TotemTransport Class Reference

#include <TotemTransport.h>

Inheritance diagram for TotemTransport:

Public Member Functions
void	process (const HepMC::GenEvent *ev, const edm::EventSetup &es, CLHEP::HepRandomEngine *engine) override
	TotemTransport (const edm::ParameterSet &ps)
	~TotemTransport () 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
LHCOpticsApproximator *	ReadParameterization (const std::string &, const std::string &)
bool	transportProton (HepMC::GenParticle *)

Private Attributes
LHCOpticsApproximator *	m_aprox_ip_150_l_
LHCOpticsApproximator *	m_aprox_ip_150_r_
double	m_beampipe_aperture_radius_
std::string	m_model_ip_150_l_name_
std::string	m_model_ip_150_r_name_

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 19 of file TotemTransport.h.

Constructor & Destructor Documentation

TotemTransport()

TotemTransport::TotemTransport

(

const edm::ParameterSet &

ps

)

Definition at line 9 of file TotemTransport.cc.

References BaseProtonTransport::beam1Filename_, BaseProtonTransport::beam2Filename_, edm::errors::Configuration, Exception, BaseProtonTransport::fPPSRegionStart_56_, edm::ParameterSet::getParameter(), m_aprox_ip_150_l_, m_aprox_ip_150_r_, m_model_ip_150_l_name_, m_model_ip_150_r_name_, BaseProtonTransport::MODE, ReadParameterization(), BaseProtonTransport::setBeamFileNames(), BaseProtonTransport::setBeamParameters(), BaseProtonTransport::setCrossingAngles(), AlCaHLTBitMon_QueryRunRegistry::string, fftjetcommon_cfi::sx, fftjetcommon_cfi::sy, and BaseProtonTransport::TOTEM.

    : BaseProtonTransport(iConfig),
      m_model_ip_150_r_name_(iConfig.getParameter<std::string>("Model_IP_150_R_Name")),
      m_model_ip_150_l_name_(iConfig.getParameter<std::string>("Model_IP_150_L_Name")),
      m_beampipe_aperture_radius_(iConfig.getParameter<double>("BeampipeApertureRadius")) {
  MODE = TransportMode::TOTEM;
  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");
  setCrossingAngles(cax45, cax56, 0.0, 0.0);
  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);

  if (fPPSRegionStart_56_ > 0)
    fPPSRegionStart_56_ *= -1;  // make sure sector 56 has negative position, as TOTEM convention

  edm::LogVerbatim("TotemTransport")
      << "=============================================================================\n"
      << "             Bulding LHC Proton transporter based on TOTEM model\n"
      << "=============================================================================\n";

  m_aprox_ip_150_r_ = ReadParameterization(m_model_ip_150_r_name_, beam1Filename_);
  m_aprox_ip_150_l_ = ReadParameterization(m_model_ip_150_l_name_, beam2Filename_);

  if (m_aprox_ip_150_r_ == nullptr || m_aprox_ip_150_l_ == nullptr) {
    edm::LogError("TotemTransport") << "Parameterisation " << m_model_ip_150_r_name_ << " or " << m_model_ip_150_l_name_
                                    << " missing in file. Cannot proceed. ";
    throw edm::Exception(edm::errors::Configuration) << "TotemTransport is not properly initialized";
  }
  edm::LogVerbatim("TotemTransport") << "Parameterizations read from file, pointers:" << m_aprox_ip_150_r_ << " "
                                     << m_aprox_ip_150_l_ << " ";
}

~TotemTransport()

TotemTransport::~TotemTransport ( )

override

Definition at line 47 of file TotemTransport.cc.

{}

Member Function Documentation

process()

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

overridevirtual

Implements BaseProtonTransport.

Definition at line 52 of file TotemTransport.cc.

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

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

  HepMC::GenEvent evt(*ievt);

  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())  // assures this protons has not been already propagated
      continue;

    transportProton(gpart);
  }
}

ReadParameterization()

LHCOpticsApproximator * TotemTransport::ReadParameterization ( const std::string &  m_model_name, const std::string &  rootfile  )

private

Definition at line 181 of file TotemTransport.cc.

References f, MillePedeFileConverter_cfg::fileName, and doHarvest::rootfile.

Referenced by TotemTransport().

                                                                                       {
  edm::FileInPath fileName(rootfile.c_str());
  TFile* f = TFile::Open(fileName.fullPath().c_str(), "read");
  if (!f) {
    edm::LogError("TotemTransport") << "File " << fileName << " not found. Exiting.";
    return nullptr;
  }
  edm::LogVerbatim("TotemTransport") << "Root file opened, pointer:" << f;

  // read parametrization
  LHCOpticsApproximator* aprox = (LHCOpticsApproximator*)f->Get(m_model_name.c_str());
  f->Close();
  return aprox;
}

transportProton()

bool TotemTransport::transportProton ( HepMC::GenParticle *  in_trk )

private

Definition at line 86 of file TotemTransport.cc.

References BaseProtonTransport::ApplyBeamCorrection(), MillePedeFileConverter_cfg::e, BaseProtonTransport::fCrossingAngleX_45_, BaseProtonTransport::fCrossingAngleX_56_, BaseProtonTransport::fPPSRegionStart_45_, BaseProtonTransport::fPPSRegionStart_56_, mps_splice::line, LogDebug, m_aprox_ip_150_l_, m_aprox_ip_150_r_, BaseProtonTransport::m_beamPart, m_beampipe_aperture_radius_, BaseProtonTransport::m_xAtTrPoint, BaseProtonTransport::m_yAtTrPoint, ProtonMassSQ, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, mathSSE::sqrt(), LHCOpticsApproximator::Transport_m_GeV(), urad, and BaseProtonTransport::verbosity_.

Referenced by process().

                                                             {
  //
  edm::LogVerbatim("TotemTransport") << "Starting proton transport using TOTEM method\n";
  //
  ApplyBeamCorrection(in_trk);

  const HepMC::GenVertex* in_pos = in_trk->production_vertex();
  const HepMC::FourVector in_mom = in_trk->momentum();
  //
  // ATTENTION: HepMC uses mm, vertex config of CMS uses cm and SimTransport uses mm
  //
  double in_position[3] = {in_pos->position().x(), in_pos->position().y(), in_pos->position().z()};  //in LHC ref. frame

  double crossingAngleX = (in_mom.z() > 0) ? fCrossingAngleX_45_ : fCrossingAngleX_56_;

  // Move the position to z=0. Do it in the CMS ref frame. Totem parameterization does the rotation internatlly
  in_position[0] =
      in_position[0] - in_position[2] * (in_mom.x() / in_mom.z() - crossingAngleX * urad);  // in CMS ref. frame
  in_position[1] = in_position[1] - in_position[2] * (in_mom.y() / (in_mom.z()));
  in_position[2] = 0.;
  double in_momentum[3] = {in_mom.x(), in_mom.y(), in_mom.z()};
  double out_position[3];
  double out_momentum[3];
  edm::LogVerbatim("TotemTransport") << "before transport ->"
                                     << " position: " << in_position[0] << ", " << in_position[1] << ", "
                                     << in_position[2] << " momentum: " << in_momentum[0] << ", " << in_momentum[1]
                                     << ", " << in_momentum[2];

  LHCOpticsApproximator* approximator = nullptr;
  double zin;
  double zout;
  if (in_mom.z() > 0) {
    approximator = m_aprox_ip_150_l_;
    zin = 0.0;  // Totem propagations assumes the starting point at 0 (zero)
    zout = fPPSRegionStart_45_;
  } else {
    approximator = m_aprox_ip_150_r_;
    zin = 0.0;  // Totem propagations assumes the starting point at 0 (zero)
    zout = fPPSRegionStart_56_;
  }

  bool invert_beam_coord_system =
      true;  // it doesn't matter the option here, it is hard coded as TRUE inside LHCOpticsApproximator!

  bool tracked = approximator->Transport_m_GeV(
      in_position, in_momentum, out_position, out_momentum, invert_beam_coord_system, zout - zin);

  if (!tracked)
    return false;

  edm::LogVerbatim("TotemTransport") << "after transport -> "
                                     << "position: " << out_position[0] << ", " << out_position[1] << ", "
                                     << out_position[2] << "momentum: " << out_momentum[0] << ", " << out_momentum[1]
                                     << ", " << out_momentum[2];

  if (out_position[0] * out_position[0] + out_position[1] * out_position[1] >
      m_beampipe_aperture_radius_ * m_beampipe_aperture_radius_) {
    edm::LogVerbatim("TotemTransport") << "Proton ouside beampipe\n"
                                       << "===== END Transport "
                                       << "====================";
    return false;
  }

  using CLHEP::meter;
  TVector3 out_pos(out_position[0] * meter, out_position[1] * meter, out_position[2] * meter);
  TVector3 out_mom(out_momentum[0], out_momentum[1], out_momentum[2]);

  if (verbosity_) {
    LogDebug("TotemTransport") << "output -> position: ";
    out_pos.Print();
    LogDebug("TotemTransport") << " momentum: ";
    out_mom.Print();
  }

  double px = -out_momentum[0];  // calculates px by means of TH_X, which is in the LHC ref. frame.
  double py = out_momentum[1];   // this need to be checked again, since it seems an invertion is occuring in the prop.
  double pz =
      out_momentum[2];  // totem calculates output pz already in the CMS ref. frame, it doesn't need to be converted
  double e = sqrt(px * px + py * py + pz * pz + ProtonMassSQ);
  TLorentzVector p_out(px, py, pz, e);
  double x1_ctpps = -out_position[0] * meter;  // Totem parameterization uses meter, one need it in millimeter
  double y1_ctpps = out_position[1] * meter;

  unsigned int line = in_trk->barcode();

  if (verbosity_)
    LogDebug("TotemTransport") << "TotemTransport: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;
}

Member Data Documentation

m_aprox_ip_150_l_

LHCOpticsApproximator* TotemTransport::m_aprox_ip_150_l_

private

Definition at line 31 of file TotemTransport.h.

Referenced by TotemTransport(), and transportProton().

m_aprox_ip_150_r_

LHCOpticsApproximator* TotemTransport::m_aprox_ip_150_r_

private

Definition at line 30 of file TotemTransport.h.

Referenced by TotemTransport(), and transportProton().

m_beampipe_aperture_radius_

double TotemTransport::m_beampipe_aperture_radius_

private

Definition at line 36 of file TotemTransport.h.

Referenced by transportProton().

m_model_ip_150_l_name_

std::string TotemTransport::m_model_ip_150_l_name_

private

Definition at line 34 of file TotemTransport.h.

Referenced by TotemTransport().

m_model_ip_150_r_name_

std::string TotemTransport::m_model_ip_150_r_name_

private

Definition at line 33 of file TotemTransport.h.

Referenced by TotemTransport().