#include <FastSimulation/CTPPSSimHitProducer/plugins/CTPPSSimHitProducer.cc>

Inheritance diagram for CTPPSSimHitProducer:

Public Member Functions
	CTPPSSimHitProducer (const edm::ParameterSet &)

	~CTPPSSimHitProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Member Functions
void	beginStream (edm::StreamID) override

void	endStream () override

void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
edm::Handle< edm::HepMCProduct >	EvtHandle

double	fz_timing

double	fz_tracker1

double	fz_tracker2

edm::EDGetTokenT< edm::HepMCProduct >	mcEventToken

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 57 of file CTPPSSimHitProducer.cc.

Constructor & Destructor Documentation

CTPPSSimHitProducer::CTPPSSimHitProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 84 of file CTPPSSimHitProducer.cc.

References fz_timing, fz_tracker1, fz_tracker2, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mcEventToken, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                        {
   produces<edm::PSimHitContainer>("CTPPSHits");
   // consumes
   mcEventToken =
       mayConsume<edm::HepMCProduct>(iConfig.getUntrackedParameter<edm::InputTag>("MCEvent", std::string("")));
 
   // Read the position of the trackers and timing detectors
   fz_tracker1 = iConfig.getParameter<double>("Z_Tracker1");
   fz_tracker2 = iConfig.getParameter<double>("Z_Tracker2");
   fz_timing = iConfig.getParameter<double>("Z_Timing");
 }

CTPPSSimHitProducer::~CTPPSSimHitProducer ( )

override

Definition at line 96 of file CTPPSSimHitProducer.cc.

96 {}

Member Function Documentation

void CTPPSSimHitProducer::beginStream ( edm::StreamID )

overrideprivate

Definition at line 210 of file CTPPSSimHitProducer.cc.

210 {}

void CTPPSSimHitProducer::endStream ( )

overrideprivate

Definition at line 213 of file CTPPSSimHitProducer.cc.

References DEFINE_FWK_MODULE.

213 {}

void CTPPSSimHitProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 103 of file CTPPSSimHitProducer.cc.

References funct::abs(), c_light_s, EvtHandle, fz_timing, fz_tracker2, edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), mps_fire::i, m_to_mm, mcEventToken, mm_to_m, eostools::move(), dqmiodumpmetadata::n, AlCaHLTBitMon_ParallelJobs::p, edm::Event::put(), MetAnalyzer::pv(), s_to_ns, and bphysicsOniaDQM_cfi::vertex.

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

                                                                                {
   using namespace edm;
 
   std::vector<PSimHit> theCTPPSHits;
   iEvent.getByToken(mcEventToken, EvtHandle);
 
   const HepMC::GenEvent* Evt = EvtHandle->GetEvent();
   std::vector<math::XYZTLorentzVector> protonCTPPS;
   protonCTPPS.clear();
   for (HepMC::GenEvent::vertex_const_iterator ivtx = Evt->vertices_begin(); ivtx != Evt->vertices_end(); ivtx++) {
     if ((*ivtx)->id() != 0)
       continue;
     double prim_vtxZ = (*ivtx)->position().z() * mm_to_m;  //in meters
     // Get the vertices at the entrance of CTPPS and get the protons coming out of them (propagated by Hector)
     for (HepMC::GenVertex::particles_out_const_iterator i = (*ivtx)->particles_out_const_begin();
          i != (*ivtx)->particles_out_const_end();
          ++i) {
       int pid = (*i)->pdg_id();
       if (pid != 2212)
         continue;
 
       HepMC::GenVertex* pv = (*i)->production_vertex();
       const HepMC::FourVector& vertex = pv->position();
       const HepMC::FourVector p((*i)->momentum());
       protonCTPPS.push_back(math::XYZTLorentzVector(p.x(), p.y(), p.z(), p.t()));
 
       LocalPoint initialPosition_tr1, initialPosition_tr2, initialPosition_tof;
 
       double t0_tr1 = 0., t0_tr2 = 0., t0_tof = 0.;
       //  t0_* has dimensions of mm
       //  Convert to ns for internal calculations.
       const double c_light_s = 2.99792458e+11;  // mm/s;
       const double s_to_ns = 1.e9;
       const double m_to_mm = 1.e3;
       double x_tr1 = 0., x_tr2 = 0., x_tof = 0., y_tr1 = 0., y_tr2 = 0., y_tof = 0., z_tr1 = 0.;
       double z_tr2 = fz_tracker2;  //m
       double z_tof = fz_timing;    //m
       int Direction = 0;
 
       // Read the vertex made by SimTransport at the entrance of det1
       if (std::abs(vertex.eta()) > 8. && (*i)->status() == 1) {
         if (vertex.z() > 0)
           Direction = 1;
         else if (vertex.z() < 0)
           Direction = -1;
 
         //Get the global coordinates at Tracker1, equal to those of the vertex at CTPPS
         x_tr1 = vertex.x();
         y_tr1 = vertex.y();
         z_tr1 = vertex.z() * mm_to_m;  //move z from mm to meters
         Local3DPoint xyzzy_tr1(x_tr1, y_tr1, z_tr1);
         initialPosition_tr1 = xyzzy_tr1;
         t0_tr1 = vertex.t() / c_light_s * s_to_ns;
         //Get the global coordinates at Tracker2 by propagating as a straight line from Tracker1
         t0_tr2 = z_tr2 * m_to_mm / c_light_s * s_to_ns;  //discuss latter if needs to be corrected with vertex position
         t0_tr2 = t0_tr1 + (z_tr2 - z_tr1) * m_to_mm / c_light_s * s_to_ns;  //corrected with vertex position
         z_tr2 *= Direction;
         x_tr2 = x_tr1 + (p.x() / p.z()) * (z_tr2 - z_tr1) * m_to_mm;
         y_tr2 = y_tr1 + (p.y() / p.z()) * (z_tr2 - z_tr1) * m_to_mm;
         Local3DPoint xyzzy_tr2(x_tr2, y_tr2, z_tr2);
         initialPosition_tr2 = xyzzy_tr2;
         //Propagate as a straight line from Tracker1
         t0_tof = z_tof * m_to_mm / c_light_s * s_to_ns;  //discuss latter if needs to be corrected with vertex position
         t0_tof = (z_tof - prim_vtxZ) * m_to_mm / c_light_s * s_to_ns;  //corrected with vertex position
         z_tof *= Direction;
         x_tof = x_tr1 + (p.x() / p.z()) * (z_tof - z_tr1) * m_to_mm;
         y_tof = y_tr1 + (p.y() / p.z()) * (z_tof - z_tr1) * m_to_mm;
         Local3DPoint xyzzy_tof(x_tof, y_tof, z_tof);
         initialPosition_tof = xyzzy_tof;
 
         // DetId codification for PSimHit from CTPPSPixel- It will be replaced by CTPPSDetId
         // 2014314496 -> Tracker1 zPositive
         // 2014838784 -> Tracker2 zPositive
         // 2046820352 -> Timing   zPositive
         // 2031091712 -> Tracker1 zNegative
         // 2031616000 -> Tracker2 zNegative
         // 2063597568 -> Timing   zNegative
 
         if (Direction > 0.) {
           PSimHit hit_tr1(xyzzy_tr1, xyzzy_tr1, 0., t0_tr1, 0., pid, 2014314496, 0, 0., 0., 2);
           PSimHit hit_tr2(xyzzy_tr2, xyzzy_tr2, 0., t0_tr2, 0., pid, 2014838784, 0, 0., 0., 2);
           PSimHit hit_tof(xyzzy_tof, xyzzy_tof, 0., t0_tof, 0., pid, 2046820352, 0, 0., 0., 2);
           theCTPPSHits.push_back(hit_tr1);
           theCTPPSHits.push_back(hit_tr2);
           theCTPPSHits.push_back(hit_tof);
         }
         if (Direction < 0.) {
           PSimHit hit_tr1(xyzzy_tr1, xyzzy_tr1, 0., t0_tr1, 0., pid, 2031091712, 0, 0., 0., 2);
           PSimHit hit_tr2(xyzzy_tr2, xyzzy_tr2, 0., t0_tr2, 0., pid, 2031616000, 0, 0., 0., 2);
           PSimHit hit_tof(xyzzy_tof, xyzzy_tof, 0., t0_tof, 0., pid, 2063597568, 0, 0., 0., 2);
           theCTPPSHits.push_back(hit_tr1);
           theCTPPSHits.push_back(hit_tr2);
           theCTPPSHits.push_back(hit_tof);
         }
       }
     }
   }
   std::unique_ptr<edm::PSimHitContainer> pctpps(new edm::PSimHitContainer);
   int n = 0;
   for (std::vector<PSimHit>::const_iterator i = theCTPPSHits.begin(); i != theCTPPSHits.end(); i++) {
     pctpps->push_back(*i);
     n += 1;
   }
   iEvent.put(std::move(pctpps), "CTPPSHits");
 }