HectorTransport Class Reference

#include <HectorTransport.h>

Inheritance diagram for HectorTransport:

Public Member Functions
	HectorTransport (const edm::ParameterSet &ps, edm::ConsumesCollector iC)
void	process (const HepMC::GenEvent *ev, const edm::EventSetup &es, CLHEP::HepRandomEngine *engine) override
	~HectorTransport () override
Public Member Functions inherited from BaseProtonTransport
void	addPartToHepMC (const HepMC::GenEvent *, 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 ()
virtual	~BaseProtonTransport ()

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

Private Attributes
const CTPPSBeamParameters *	beamParameters_ {nullptr}
edm::ESGetToken< CTPPSBeamParameters, CTPPSBeamParametersRcd >	beamParametersToken_
const BeamSpotObjects *	beamspot_ {nullptr}
edm::ESGetToken< BeamSpotObjects, BeamSpotObjectsRcd >	beamspotToken_
std::unique_ptr< H_BeamLine >	m_beamline45
std::unique_ptr< H_BeamLine >	m_beamline56

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 Member Functions inherited from BaseProtonTransport
void	setBeamFileNames (const std::string &nam1, const std::string &nam2)
void	setBeamParameters (double stx, double sty, double sx, double sy, double se)
void	setCrossingAngles (double cx45, double cx56, double cy45, double cy56)
Protected Attributes inherited from BaseProtonTransport
bool	bApplyZShift_
std::string	beam1Filename_ {""}
std::string	beam2Filename_ {""}
CLHEP::HepRandomEngine *	engine_ {nullptr}
double	etaCut_
double	fCrossingAngleX_45_ {0.0}
double	fCrossingAngleX_56_ {0.0}
double	fCrossingAngleY_45_ {0.0}
double	fCrossingAngleY_56_ {0.0}
double	fPPSRegionStart_45_
double	fPPSRegionStart_56_
std::map< unsigned int, TLorentzVector >	m_beamPart
std::vector< LHCTransportLink >	m_CorrespondenceMap
std::map< unsigned int, double >	m_xAtTrPoint
std::map< unsigned int, double >	m_yAtTrPoint
TransportMode	MODE
double	momentumCut_
bool	produceHitsRelativeToBeam_
bool	useBeamPositionFromLHCInfo_
bool	verbosity_

Detailed Description

Definition at line 39 of file HectorTransport.h.

Constructor & Destructor Documentation

HectorTransport()

HectorTransport::HectorTransport	(	const edm::ParameterSet &	ps,
		edm::ConsumesCollector	iC
	)

Definition at line 11 of file HectorTransport.cc.

References edm::ParameterSet::getParameter(), BaseProtonTransport::HECTOR, BaseProtonTransport::MODE, BaseProtonTransport::setBeamFileNames(), setBeamLine(), BaseProtonTransport::setBeamParameters(), BaseProtonTransport::setCrossingAngles(), AlCaHLTBitMon_QueryRunRegistry::string, fftjetcommon_cfi::sx, and fftjetcommon_cfi::sy.

    : BaseProtonTransport(iConfig),
      m_beamline45(nullptr),
      m_beamline56(nullptr),
      beamParametersToken_(iC.esConsumes()),
      beamspotToken_(iC.esConsumes()) {
  MODE = TransportMode::HECTOR;
  std::string s1 = iConfig.getParameter<std::string>("Beam1Filename");
  std::string s2 = iConfig.getParameter<std::string>("Beam2Filename");
  setBeamFileNames(s1, s2);
  double cax45 = iConfig.getParameter<double>("halfCrossingAngleXSector45");
  double cax56 = iConfig.getParameter<double>("halfCrossingAngleXSector56");
  double cay45 = iConfig.getParameter<double>("halfCrossingAngleYSector45");
  double cay56 = iConfig.getParameter<double>("halfCrossingAngleYSector56");
  setCrossingAngles(cax45, cax56, cay45, cay56);
  double stx = iConfig.getParameter<double>("BeamDivergenceX");
  double sty = iConfig.getParameter<double>("BeamDivergenceY");
  double sx = iConfig.getParameter<double>("BeamSigmaX");
  double sy = iConfig.getParameter<double>("BeamSigmaY");
  double se = iConfig.getParameter<double>("BeamEnergyDispersion");
  setBeamParameters(stx, sty, sx, sy, se);
  //PPS
  edm::LogVerbatim("ProtonTransport")
      << "=============================================================================\n"
      << "             Bulding LHC Proton transporter based on HECTOR model\n"
      << "=============================================================================\n";
  setBeamLine();
}

~HectorTransport()

HectorTransport::~HectorTransport ( )

override

Definition at line 39 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 44 of file HectorTransport.cc.

References beamParameters_, beamParametersToken_, beamspot_, beamspotToken_, BaseProtonTransport::clear(), BaseProtonTransport::engine_, BaseProtonTransport::etaCut_, GenParticle::GenParticle, edm::EventSetup::getData(), mps_splice::line, BaseProtonTransport::m_beamPart, BaseProtonTransport::momentumCut_, and transportProton().

                                                            {
  clear();
  engine_ = engine;  // the engine needs to be updated for each event

  beamspot_ = &iSetup.getData(beamspotToken_);
  beamParameters_ = &iSetup.getData(beamParametersToken_);

  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);
  }
}

setBeamLine()

bool HectorTransport::setBeamLine ( )

private

Definition at line 187 of file HectorTransport.cc.

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().

                                  {
  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::make_unique<H_BeamLine>(-1, fPPSBeamLineLength_ + 0.1);
    m_beamline45->fill(b2.fullPath(), 1, "IP5");
    m_beamline56 = std::make_unique<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";
    edm::LogVerbatim("HectorTransportSetup") << "===================================================================\n"
                                             << "                  Forward beam line elements \n";
    m_beamline45->showElements();
    edm::LogVerbatim("HectorTransportSetup") << "===================================================================\n"
                                             << "                 Backward beam line elements \n";
    m_beamline56->showElements();
  }
  return true;
}

transportProton()

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

private

propagate the particles through a beamline to PPS

Definition at line 75 of file HectorTransport.cc.

References BaseProtonTransport::beamEnergy(), BaseProtonTransport::beamMomentum(), beamParameters_, beamspot_, ALCARECOTkAlJpsiMuMu_cff::charge, cm_to_mm, funct::cos(), MillePedeFileConverter_cfg::e, BaseProtonTransport::fCrossingAngleX_45_, BaseProtonTransport::fCrossingAngleX_56_, BaseProtonTransport::fCrossingAngleY_45_, BaseProtonTransport::fCrossingAngleY_56_, BaseProtonTransport::fPPSRegionStart_45_, BaseProtonTransport::fPPSRegionStart_56_, CTPPSBeamParameters::getVtxOffsetX45(), CTPPSBeamParameters::getVtxOffsetY45(), mps_splice::line, LogDebug, m_beamline45, m_beamline56, BaseProtonTransport::m_beamPart, BaseProtonTransport::m_xAtTrPoint, BaseProtonTransport::m_yAtTrPoint, EgHLTOffHistBins_cfi::mass, mm_to_um, funct::pow(), BaseProtonTransport::produceHitsRelativeToBeam_, ProtonMassSQ, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, mathSSE::sqrt(), funct::tan(), theta(), um_to_mm, urad, BaseProtonTransport::useBeamPositionFromLHCInfo_, BaseProtonTransport::verbosity_, BeamSpotObjects::x(), and BeamSpotObjects::y().

Referenced by process().

                                                                   {
  edm::LogVerbatim("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.;

  // remove the CMS vertex shift
  double vtxXoffset;
  double vtxYoffset;
  if (useBeamPositionFromLHCInfo_) {
    vtxXoffset = beamParameters_->getVtxOffsetX45() * cm_to_mm;
    vtxYoffset = beamParameters_->getVtxOffsetY45() * cm_to_mm;
  } else {
    vtxXoffset = beamspot_->x() * cm_to_mm;
    vtxYoffset = beamspot_->y() * cm_to_mm;
  }

  // Momentum in LHC ref. frame
  px = -gpart->momentum().px();
  py = gpart->momentum().py();
  pz = -gpart->momentum().pz();
  e = gpart->momentum().e();

  int direction = (pz > 0) ? 1 : -1;  // in relation to LHC ref frame

  double XforPosition = -gpart->production_vertex()->position().x();  //mm in the ref. frame of LHC
  double YforPosition = gpart->production_vertex()->position().y();   //mm
  double ZforPosition = -gpart->production_vertex()->position().z();  //mm
  double fCrossingAngleX = (pz < 0) ? fCrossingAngleX_45_ : fCrossingAngleX_56_;
  double fCrossingAngleY = (pz < 0) ? fCrossingAngleY_45_ : fCrossingAngleY_56_;
  //
  H_BeamParticle h_p(mass, charge);
  h_p.set4Momentum(px, py, pz, e);
  //
  // shift the starting position of the track to Z=0 if configured so (all the variables below are still in cm)
  XforPosition = XforPosition - ZforPosition * (px / pz + fCrossingAngleX * urad);  // theta ~=tan(theta)
  YforPosition = YforPosition - ZforPosition * (py / pz + fCrossingAngleY * urad);
  XforPosition -= (-vtxXoffset);  // X was already in the LHC ref. frame
  YforPosition -= vtxYoffset;

  // set position, but do not invert the coordinate for the angles (TX,TY) because it is done by set4Momentum
  h_p.setPosition(XforPosition * mm_to_um, YforPosition * mm_to_um, h_p.getTX(), h_p.getTY(), 0.);

  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);
  x1_ctpps = h_p.getX();
  y1_ctpps = h_p.getY();

  double thx = h_p.getTX();
  double thy = h_p.getTY();

  if (produceHitsRelativeToBeam_) {
    H_BeamParticle p_beam(mass, charge);
    p_beam.set4Momentum(0., 0., beamMomentum(), beamEnergy());
    p_beam.setPosition(0., 0., fCrossingAngleX * urad, fCrossingAngleY * urad, 0.);
    p_beam.computePath(&*beamline);
    thx -= p_beam.getTX();
    thy -= p_beam.getTY();
    x1_ctpps -= p_beam.getX();
    y1_ctpps -= p_beam.getY();
    edm::LogVerbatim("HectorTransportEventProcessing")
        << "Shifting the hit relative to beam :  X = " << p_beam.getX() << "(mm)      Y = " << p_beam.getY() << "(mm)";
  }

  double partP = sqrt(pow(h_p.getE(), 2) - ProtonMassSQ);
  double theta = sqrt(thx * thx + thy * thy) * urad;

  // copy the kinematic changing to CMS ref. frame, only the negative Pz needs to be changed
  TLorentzVector p_out(-tan(thx * urad) * partP * cos(theta),
                       tan(thy * urad) * partP * cos(theta),
                       -direction * partP * cos(theta),
                       h_p.getE());

  m_beamPart[line] = p_out;
  m_xAtTrPoint[line] = -x1_ctpps * um_to_mm;  // move to CMS ref. frame
  m_yAtTrPoint[line] = y1_ctpps * um_to_mm;
  if (verbosity_)
    LogDebug("HectorTransportEventProcessing")
        << "HectorTransport:filterPPS: barcode = " << line << " x=  " << x1_ctpps << " y= " << y1_ctpps;

  return true;
}

Member Data Documentation

beamParameters_

const CTPPSBeamParameters* HectorTransport::beamParameters_ {nullptr}

private

Definition at line 62 of file HectorTransport.h.

Referenced by process(), and transportProton().

beamParametersToken_

edm::ESGetToken<CTPPSBeamParameters, CTPPSBeamParametersRcd> HectorTransport::beamParametersToken_

private

Definition at line 59 of file HectorTransport.h.

Referenced by process().

beamspot_

const BeamSpotObjects* HectorTransport::beamspot_ {nullptr}

private

Definition at line 63 of file HectorTransport.h.

Referenced by process(), and transportProton().

beamspotToken_

edm::ESGetToken<BeamSpotObjects, BeamSpotObjectsRcd> HectorTransport::beamspotToken_

private

Definition at line 60 of file HectorTransport.h.

Referenced by process().

fPPSBeamLineLength_

constexpr double HectorTransport::fPPSBeamLineLength_ = 250.

staticprivate

Definition at line 47 of file HectorTransport.h.

Referenced by setBeamLine().

m_beamline45

std::unique_ptr<H_BeamLine> HectorTransport::m_beamline45

private

Definition at line 56 of file HectorTransport.h.

Referenced by setBeamLine(), and transportProton().

m_beamline56

std::unique_ptr<H_BeamLine> HectorTransport::m_beamline56

private

Definition at line 57 of file HectorTransport.h.

Referenced by setBeamLine(), and transportProton().