HectorTransport Class Reference

#include <HectorTransport.h>

Inheritance diagram for HectorTransport:

Public Member Functions
	HectorTransport (const edm::ParameterSet &ps)
void	process (const HepMC::GenEvent *ev, const edm::EventSetup &es, CLHEP::HepRandomEngine *engine) override
	~HectorTransport () override
Public Member Functions inherited from BaseProtonTransport
void	addPartToHepMC (HepMC::GenEvent *)
void	ApplyBeamCorrection (HepMC::GenParticle *p)
void	ApplyBeamCorrection (TLorentzVector &p)
	BaseProtonTransport (const edm::ParameterSet &iConfig)
double	beamEnergy ()
double	beamMomentum ()
void	clear ()
std::vector< LHCTransportLink > &	getCorrespondenceMap ()
void	setBeamEnergy (double e)
virtual	~BaseProtonTransport ()

Private Member Functions
bool	setBeamLine ()
bool	transportProton (const HepMC::GenParticle *)
	propagate the particles through a beamline to PPS More...

Private Attributes
std::unique_ptr< H_BeamLine >	m_beamline45
std::unique_ptr< H_BeamLine >	m_beamline56
double	m_fEtacut
double	m_fMomentumMin
bool	m_verbosity

Static Private Attributes
static constexpr double	fPPSBeamLineLength_ = 250.

Additional Inherited Members
Protected Types inherited from BaseProtonTransport
enum	TransportMode { TransportMode::HECTOR, TransportMode::TOTEM, TransportMode::OPTICALFUNCTIONS }
Protected Attributes inherited from BaseProtonTransport
bool	bApplyZShift
std::string	beam1Filename_
std::string	beam2Filename_
double	beamEnergy_
double	beamMomentum_
CLHEP::HepRandomEngine *	engine_
double	etaCut_
double	fBeamXatIP
double	fBeamYatIP
double	fCrossingAngle_45
double	fCrossingAngle_56
double	fPPSRegionStart_45
double	fPPSRegionStart_56
double	fVtxMeanX
double	fVtxMeanY
double	fVtxMeanZ
std::map< unsigned int, TLorentzVector >	m_beamPart
std::vector< LHCTransportLink >	m_CorrespondenceMap
double	m_sig_E
double	m_sigmaSTX
double	m_sigmaSTY
double	m_sigmaSX
double	m_sigmaSY
std::map< unsigned int, double >	m_xAtTrPoint
std::map< unsigned int, double >	m_yAtTrPoint
TransportMode	MODE
double	momentumCut_
int	NEvent
bool	verbosity_

Detailed Description

Definition at line 35 of file HectorTransport.h.

Constructor & Destructor Documentation

HectorTransport()

HectorTransport::HectorTransport

(

const edm::ParameterSet &

ps

)

Definition at line 12 of file HectorTransport.cc.

    : BaseProtonTransport(iConfig), m_beamline45(nullptr), m_beamline56(nullptr) {
  // Create LHC beam line
  MODE = TransportMode::HECTOR;  // defines the MODE for the transport
  beam1Filename_ = iConfig.getParameter<std::string>("Beam1Filename");
  beam2Filename_ = iConfig.getParameter<std::string>("Beam2Filename");
  fCrossingAngle_45 = iConfig.getParameter<double>("halfCrossingAngleSector45");
  fCrossingAngle_56 = iConfig.getParameter<double>("halfCrossingAngleSector56");
  beamEnergy_ = iConfig.getParameter<double>("BeamEnergy");
  fVtxMeanX = iConfig.getParameter<double>("VtxMeanX");
  fVtxMeanY = iConfig.getParameter<double>("VtxMeanY");
  fVtxMeanZ = iConfig.getParameter<double>("VtxMeanZ");
  m_sigmaSTX = iConfig.getParameter<double>("BeamDivergenceX");
  m_sigmaSTY = iConfig.getParameter<double>("BeamDivergenceY");
  m_sigmaSX = iConfig.getParameter<double>("BeamSigmaX");
  m_sigmaSY = iConfig.getParameter<double>("BeamSigmaY");
  m_sig_E = iConfig.getParameter<double>("BeamEnergyDispersion");
  fBeamXatIP = iConfig.getUntrackedParameter<double>("BeamXatIP", fVtxMeanX);
  fBeamYatIP = iConfig.getUntrackedParameter<double>("BeamYatIP", fVtxMeanY);
 
  //PPS
  edm::LogInfo("ProtonTransport") << "=============================================================================\n"
                                  << "             Bulding LHC Proton transporter based on HECTOR model\n"
                                  << "=============================================================================\n";
  setBeamLine();
}

References BaseProtonTransport::beam1Filename_, BaseProtonTransport::beam2Filename_, BaseProtonTransport::beamEnergy_, BaseProtonTransport::fBeamXatIP, BaseProtonTransport::fBeamYatIP, BaseProtonTransport::fCrossingAngle_45, BaseProtonTransport::fCrossingAngle_56, BaseProtonTransport::fVtxMeanX, BaseProtonTransport::fVtxMeanY, BaseProtonTransport::fVtxMeanZ, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), BaseProtonTransport::HECTOR, BaseProtonTransport::m_sig_E, BaseProtonTransport::m_sigmaSTX, BaseProtonTransport::m_sigmaSTY, BaseProtonTransport::m_sigmaSX, BaseProtonTransport::m_sigmaSY, BaseProtonTransport::MODE, setBeamLine(), and AlCaHLTBitMon_QueryRunRegistry::string.

~HectorTransport()

HectorTransport::~HectorTransport ( )

override

Definition at line 10 of file HectorTransport.cc.

{};

Member Function Documentation

process()

void HectorTransport::process ( const HepMC::GenEvent *  ev, const edm::EventSetup &  es, CLHEP::HepRandomEngine *  engine  )

overridevirtual

Implements BaseProtonTransport.

Definition at line 42 of file HectorTransport.cc.

                                                             {
  this->clear();
 
  engine_ = _engine;  // the engine needs to be updated for each event
 
  for (HepMC::GenEvent::particle_const_iterator eventParticle = evt->particles_begin();
       eventParticle != evt->particles_end();
       ++eventParticle) {
    if (!((*eventParticle)->status() == 1 && (*eventParticle)->pdg_id() == 2212))
      continue;
 
    if (!(fabs((*eventParticle)->momentum().eta()) > etaCut_ && fabs((*eventParticle)->momentum().pz()) > momentumCut_))
      continue;  // discard protons outside kinematic acceptance
 
    unsigned int line = (*eventParticle)->barcode();
    HepMC::GenParticle* gpart = (*eventParticle);
 
    if (gpart->pdg_id() != 2212)
      continue;  // only transport stable protons
    if (gpart->status() != 1)
      continue;
    if (m_beamPart.find(line) != m_beamPart.end())
      continue;  // assures this protons has not been already propagated
 
    transportProton(gpart);
  }
}

References BaseProtonTransport::clear(), BaseProtonTransport::engine_, BaseProtonTransport::etaCut_, GenParticle::GenParticle, mps_splice::line, BaseProtonTransport::m_beamPart, BaseProtonTransport::momentumCut_, and transportProton().

setBeamLine()

bool HectorTransport::setBeamLine ( )

private

Definition at line 174 of file HectorTransport.cc.

                                  {
  edm::FileInPath b1(beam1Filename_.c_str());
  edm::FileInPath b2(beam2Filename_.c_str());
 
  // construct beam line for PPS (forward 1 backward 2):
  if (fPPSBeamLineLength_ > 0.) {
    m_beamline45 = std::unique_ptr<H_BeamLine>(new H_BeamLine(-1, fPPSBeamLineLength_ + 0.1));  // (direction, length)
    m_beamline45->fill(b2.fullPath(), 1, "IP5");
    m_beamline56 = std::unique_ptr<H_BeamLine>(new H_BeamLine(1, fPPSBeamLineLength_ + 0.1));  //
    m_beamline56->fill(b1.fullPath(), 1, "IP5");
  } else {
    if (verbosity_)
      LogDebug("HectorTransportSetup") << "HectorTransport: WARNING: lengthctpps=  " << fPPSBeamLineLength_;
    return false;
  }
  if (verbosity_) {
    edm::LogInfo("HectorTransportSetup") << "===================================================================\n"
                                         << " * * * * * * * * * * * * * * * * * * * * * * * * * * * *           \n"
                                         << " *                                                         *       \n"
                                         << " *                   --<--<--  A fast simulator --<--<--     *     \n"
                                         << " *                 | --<--<--     of particle   --<--<--     *     \n"
                                         << " *  ----HECTOR----<                                          *     \n"
                                         << " *                 | -->-->-- transport through-->-->--      *     \n"
                                         << " *                   -->-->-- generic beamlines -->-->--     *     \n"
                                         << " *                                                           *     \n"
                                         << " * JINST 2:P09005 (2007)                                     *     \n"
                                         << " *      X Rouby, J de Favereau, K Piotrzkowski (CP3)         *     \n"
                                         << " *       http://www.fynu.ucl.ac.be/hector.html               *     \n"
                                         << " *                                                           *     \n"
                                         << " * Center for Cosmology, Particle Physics and Phenomenology  *     \n"
                                         << " *              Universite catholique de Louvain             *     \n"
                                         << " *                 Louvain-la-Neuve, Belgium                 *     \n"
                                         << " *                                                         *       \n"
                                         << " * * * * * * * * * * * * * * * * * * * * * * * * * * * *           \n"
                                         << " HectorTransport configuration: \n"
                                         << " Beam line length      = " << fPPSBeamLineLength_ << "\n"
                                         << " PPS Region Start 44   =  " << fPPSRegionStart_45 << "\n"
                                         << " PPS Region Start 56   =  " << fPPSRegionStart_56 << "\n"
                                         << "===================================================================\n";
  }
  if (verbosity_) {
    edm::LogInfo("HectorTransportSetup") << "====================================================================\n"
                                         << "                  Forward beam line elements \n";
    m_beamline45->showElements();
    edm::LogInfo("HectorTransportSetup") << "====================================================================\n"
                                         << "                 Backward beam line elements \n";
    m_beamline56->showElements();
  }
  return true;
}

References b1, b2, BaseProtonTransport::beam1Filename_, BaseProtonTransport::beam2Filename_, fPPSBeamLineLength_, BaseProtonTransport::fPPSRegionStart_45, BaseProtonTransport::fPPSRegionStart_56, LogDebug, m_beamline45, m_beamline56, and BaseProtonTransport::verbosity_.

Referenced by HectorTransport().

transportProton()

bool HectorTransport::transportProton ( const HepMC::GenParticle *  gpart )

private

propagate the particles through a beamline to PPS

Definition at line 71 of file HectorTransport.cc.

                                                                   {
  edm::LogInfo("ProtonTransport") << "Starting proton transport using HECTOR method\n";
 
  double px, py, pz, e;
  unsigned int line = (gpart)->barcode();
 
  double mass = gpart->generatedMass();
  double charge = 1.;
 
  px = gpart->momentum().px();
  py = gpart->momentum().py();
  pz = gpart->momentum().pz();
  e = gpart->momentum().e();
 
  e = sqrt(pow(mass, 2) + pow(px, 2) + pow(py, 2) + pow(pz, 2));
 
  int direction = (pz > 0) ? 1 : -1;
 
  // Apply Beam and Crossing Angle Corrections
  TLorentzVector p_out(px, py, pz, e);
  PPSTools::LorentzBoost(p_out,
                         "LAB",
                         {fCrossingAngle_56,  //Beam1
                          fCrossingAngle_45,  //Beam2
                          beamMomentum_,
                          beamEnergy_});
 
  ApplyBeamCorrection(p_out);
 
  // from mm to cm
  double XforPosition = gpart->production_vertex()->position().x() / cm;  //cm
  double YforPosition = gpart->production_vertex()->position().y() / cm;  //cm
  double ZforPosition = gpart->production_vertex()->position().z() / cm;  //cm
 
  H_BeamParticle h_p(mass, charge);
  h_p.set4Momentum(-direction * p_out.Px(), p_out.Py(), fabs(p_out.Pz()), p_out.E());
  // shift the beam position to the given beam position at IP (in cm)
  XforPosition = (XforPosition - fVtxMeanX) + fBeamXatIP * mm_to_cm;
  YforPosition = (YforPosition - fVtxMeanY) + fBeamYatIP * mm_to_cm;
  //ZforPosition stays the same, move the closed orbit only in the X,Y plane
  //
  // shift the starting position of the track to Z=0 if configured so (all the variables below are still in cm)
  if (bApplyZShift) {
    double fCrossingAngle = (p_out.Pz() > 0) ? fCrossingAngle_45 : -fCrossingAngle_56;
    XforPosition = XforPosition +
                   (tan((long double)fCrossingAngle * urad) - ((long double)p_out.Px()) / ((long double)p_out.Pz())) *
                       ZforPosition;
    YforPosition = YforPosition - ((long double)p_out.Py()) / ((long double)p_out.Pz()) * ZforPosition;
    ZforPosition = 0.;
  }
 
  //
  // set position, but do not invert the coordinate for the angles (TX,TY) because it is done by set4Momentum
  h_p.setPosition(-direction * XforPosition * cm_to_um,
                  YforPosition * cm_to_um,
                  h_p.getTX(),
                  h_p.getTY(),
                  -direction * ZforPosition * cm_to_m);
  bool is_stop;
  float x1_ctpps;
  float y1_ctpps;
 
  H_BeamLine* _beamline = nullptr;
  double _targetZ = 0;
  switch (direction) {
    case -1:
      _beamline = &*m_beamline56;  // negative side propagation
      _targetZ = fPPSRegionStart_56;
      break;
    case 1:
      _beamline = &*m_beamline45;
      _targetZ = fPPSRegionStart_45;
      break;
  }
  // insert protection for NULL beamlines here
  h_p.computePath(&*_beamline);
  is_stop = h_p.stopped(&*_beamline);
  if (verbosity_)
    LogDebug("HectorTransportEventProcessing")
        << "HectorTransport:filterPPS: barcode = " << line << " is_stop=  " << is_stop;
  if (is_stop) {
    return false;
  }
  //
  //propagating
  //
  h_p.propagate(_targetZ);
 
  p_out = PPSTools::HectorParticle2LorentzVector(h_p, direction);
 
  p_out.SetPx(direction * p_out.Px());
  x1_ctpps = direction * h_p.getX() * um_to_mm;
  y1_ctpps = h_p.getY() * um_to_mm;
 
  if (verbosity_)
    LogDebug("HectorTransportEventProcessing")
        << "HectorTransport:filterPPS: barcode = " << line << " x=  " << x1_ctpps << " y= " << y1_ctpps;
 
  m_beamPart[line] = p_out;
  m_xAtTrPoint[line] = x1_ctpps;
  m_yAtTrPoint[line] = y1_ctpps;
  return true;
}

References BaseProtonTransport::ApplyBeamCorrection(), BaseProtonTransport::bApplyZShift, BaseProtonTransport::beamEnergy_, BaseProtonTransport::beamMomentum_, ALCARECOTkAlJpsiMuMu_cff::charge, cm_to_m, cm_to_um, MillePedeFileConverter_cfg::e, BaseProtonTransport::fBeamXatIP, BaseProtonTransport::fBeamYatIP, BaseProtonTransport::fCrossingAngle_45, BaseProtonTransport::fCrossingAngle_56, BaseProtonTransport::fPPSRegionStart_45, BaseProtonTransport::fPPSRegionStart_56, BaseProtonTransport::fVtxMeanX, BaseProtonTransport::fVtxMeanY, PPSTools::HectorParticle2LorentzVector(), mps_splice::line, LogDebug, PPSTools::LorentzBoost(), m_beamline45, m_beamline56, BaseProtonTransport::m_beamPart, BaseProtonTransport::m_xAtTrPoint, BaseProtonTransport::m_yAtTrPoint, EgHLTOffHistBins_cfi::mass, mm_to_cm, funct::pow(), multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, mathSSE::sqrt(), funct::tan(), um_to_mm, urad, and BaseProtonTransport::verbosity_.

Referenced by process().

Member Data Documentation

fPPSBeamLineLength_

constexpr double HectorTransport::fPPSBeamLineLength_ = 250.

staticconstexprprivate

Definition at line 44 of file HectorTransport.h.

Referenced by setBeamLine().

m_beamline45

std::unique_ptr<H_BeamLine> HectorTransport::m_beamline45

private

Definition at line 58 of file HectorTransport.h.

Referenced by setBeamLine(), and transportProton().

m_beamline56

std::unique_ptr<H_BeamLine> HectorTransport::m_beamline56

private

Definition at line 59 of file HectorTransport.h.

Referenced by setBeamLine(), and transportProton().

m_fEtacut

double HectorTransport::m_fEtacut

private

Definition at line 54 of file HectorTransport.h.

m_fMomentumMin

double HectorTransport::m_fMomentumMin

private

Definition at line 55 of file HectorTransport.h.

m_verbosity

bool HectorTransport::m_verbosity

private

Definition at line 43 of file HectorTransport.h.