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

Inheritance diagram for CTPPSSimHitProducer:

Public Types
typedef CLHEP::HepLorentzVector	LorentzVector

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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

	~CTPPSSimHitProducer ()

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

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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

virtual void	endStream () override

virtual 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
Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 58 of file CTPPSSimHitProducer.cc.

Member Typedef Documentation

typedef CLHEP::HepLorentzVector CTPPSSimHitProducer::LorentzVector

Definition at line 62 of file CTPPSSimHitProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 89 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 ( )

Definition at line 102 of file CTPPSSimHitProducer.cc.

 {
 
 }

Member Function Documentation

void CTPPSSimHitProducer::beginStream ( edm::StreamID )

overrideprivatevirtual

Reimplemented from edm::stream::EDProducerBase.

Definition at line 216 of file CTPPSSimHitProducer.cc.

217 {

218 }

void CTPPSSimHitProducer::endStream ( )

overrideprivatevirtual

Reimplemented from edm::stream::EDProducerBase.

Definition at line 222 of file CTPPSSimHitProducer.cc.

References DEFINE_FWK_MODULE.

222 {

223 }

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

overrideprivatevirtual

Definition at line 114 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(), gen::n, AlCaHLTBitMon_ParallelJobs::p, sysUtil::pid, edm::Event::put(), MetAnalyzer::pv(), and s_to_ns.

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;
 
         // 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();
             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
                 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
                 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");
 }