Producer for PFRecTracks and PFSimParticles. More...

#include <PFSimParticleProducer.h>

Inheritance diagram for PFSimParticleProducer:

Public Types
typedef edm::Handle< reco::PFRecTrackCollection >	TrackHandle

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
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
void	getSimIDs (const TrackHandle &trackh, std::vector< unsigned > &recTrackSimID)

	PFSimParticleProducer (const edm::ParameterSet &)

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

	~PFSimParticleProducer () override

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

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

ModuleDescription const &	moduleDescription () const

	~EDProducerBase () override

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

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

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, ModuleToResolverIndicies 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 Attributes
edm::InputTag	inputTagEcalRecHitsEB_

edm::InputTag	inputTagEcalRecHitsEE_

edm::InputTag	inputTagFamosSimHits_

edm::InputTag	inputTagRecTracks_

edm::InputTag	inputTagSim_
	module label for retrieving input simtrack and simvertex More...

bool	mctruthMatchingInfo_

FSimEvent *	mySimEvent

edm::ParameterSet	particleFilter_

bool	processParticles_
	process particles on/off More...

edm::EDGetTokenT< EcalRecHitCollection >	tokenEcalRecHitsEB_

edm::EDGetTokenT< EcalRecHitCollection >	tokenEcalRecHitsEE_

edm::EDGetTokenT< edm::PCaloHitContainer >	tokenFamosSimHits_

edm::EDGetTokenT< reco::PFRecTrackCollection >	tokenRecTracks_

edm::EDGetTokenT< std::vector< SimTrack > >	tokenSim_

edm::EDGetTokenT< std::vector< SimVertex > >	tokenSimVertices_

bool	verbose_
	verbose ? More...

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

Producer for PFRecTracks and PFSimParticles.

Author: Colin Bernet

Date: April 2007

Definition at line 44 of file PFSimParticleProducer.h.

Member Typedef Documentation

typedef edm::Handle<reco::PFRecTrackCollection> PFSimParticleProducer::TrackHandle

Definition at line 53 of file PFSimParticleProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 66 of file PFSimParticleProducer.cc.

References edm::ParameterSet::getParameter(), and edm::ParameterSet::getUntrackedParameter().

 {
 
 
   processParticles_ = 
     iConfig.getUntrackedParameter<bool>("process_Particles",true);
     
   
   inputTagSim_ = iConfig.getParameter<InputTag>("sim"); 
   tokenSim_    = consumes<std::vector<SimTrack> >(inputTagSim_);
   tokenSimVertices_    = consumes<std::vector<SimVertex> >(inputTagSim_);
 
   //retrieving collections for MC Truth Matching
 
   //modif-beg
   inputTagFamosSimHits_  = iConfig.getUntrackedParameter<InputTag>("famosSimHits");
   tokenFamosSimHits_  = consumes<edm::PCaloHitContainer>(inputTagFamosSimHits_);
   mctruthMatchingInfo_ = 
     iConfig.getUntrackedParameter<bool>("MCTruthMatchingInfo",false);
   //modif-end
 
   inputTagRecTracks_ = iConfig.getParameter<InputTag>("RecTracks");
   tokenRecTracks_    = consumes<reco::PFRecTrackCollection>(inputTagRecTracks_);
 
   inputTagEcalRecHitsEB_ = iConfig.getParameter<InputTag>("ecalRecHitsEB");
   tokenEcalRecHitsEB_    = consumes<EcalRecHitCollection>(inputTagEcalRecHitsEB_);
   inputTagEcalRecHitsEE_ = iConfig.getParameter<InputTag>("ecalRecHitsEE");
   tokenEcalRecHitsEE_ 
     = consumes<EcalRecHitCollection>(inputTagEcalRecHitsEE_);
 
   verbose_ = 
     iConfig.getUntrackedParameter<bool>("verbose",false);
 
 
   // register products
   produces<reco::PFSimParticleCollection>();
   
   particleFilter_ = iConfig.getParameter<ParameterSet>
     ( "ParticleFilter" );   
   
   mySimEvent =  new FSimEvent( particleFilter_ );
 }

PFSimParticleProducer::~PFSimParticleProducer ( )

override

Definition at line 111 of file PFSimParticleProducer.cc.

 { 
   delete mySimEvent; 
 }

Member Function Documentation

void PFSimParticleProducer::getSimIDs	(	const TrackHandle &	trackh,
		std::vector< unsigned > &	recTrackSimID
	)

Definition at line 527 of file PFSimParticleProducer.cc.

References mps_fire::i, edm::HandleBase::isValid(), FastTrackerRecHit::nSimTrackIds(), FastTrackerRecHit::simTrackId(), and reco::PFRecTrack::trackRef().

 {
 
   if( trackh.isValid() ) {
 //     cout << "Size=" << trackh->size() << endl;
 
     for(unsigned i=0;i<trackh->size(); i++) {
       
       reco::PFRecTrackRef ref( trackh,i );
       const reco::PFRecTrack& PFT   = *ref;
       const reco::TrackRef& trackref = PFT.trackRef();
 
 //       double Pt  = trackref->pt(); 
 //       double DPt = trackref->ptError();
 //       cout << " PFBlockProducer: PFrecTrack->Track Pt= " 
 //     << Pt << " DPt = " << DPt << endl;
       
       trackingRecHit_iterator rhitbeg 
     = trackref->recHitsBegin();
       trackingRecHit_iterator rhitend 
     = trackref->recHitsEnd();
       for (trackingRecHit_iterator it = rhitbeg;  
        it != rhitend; it++){
 
     if( (*it)->isValid() ){
 
       const FastTrackerRecHit * rechit 
         = (const FastTrackerRecHit*) (*it);
       
 //    cout <<  "rechit"            
 //         << " corresponding simId " 
 //         << rechit->simtrackId() 
 //         << endl;
 
       for(unsigned int st_index=0;st_index<rechit->nSimTrackIds();++st_index){
           recTrackSimID.push_back(rechit->simTrackId(st_index));
       }
       break;
     }
       }//loop track rechit
     }//loop recTracks
   }//track handle valid
 }

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

override

Definition at line 116 of file PFSimParticleProducer.cc.

 {  
   // init Particle data table (from Pythia)
   edm::ESHandle < HepPDT::ParticleDataTable > pdt;
   //  edm::ESHandle < DefaultConfig::ParticleDataTable > pdt;
   iSetup.getData(pdt);
   mySimEvent->initializePdt(&(*pdt));
 
   ParticleTable::Sentry ptable(mySimEvent->theTable());
   LogDebug("PFSimParticleProducer")<<"START event: "<<iEvent.id().event()
                    <<" in run "<<iEvent.id().run()<<endl;
  
   //MC Truth Matching only with Famos and UnFoldedMode option to true!!
 
   //vector to store the trackIDs of rectracks corresponding
   //to the simulated particle.
   std::vector<unsigned> recTrackSimID;
   
   //In order to know which simparticule contribute to 
   //a given Ecal RecHit energy, we need to access 
   //the PCAloHit from FastSim. 
   
   typedef std::pair<double, unsigned> hitSimID;
   typedef std::list< std::pair<double, unsigned> >::iterator ITM;
   std::vector< std::list <hitSimID> > caloHitsEBID(62000);
   std::vector<double> caloHitsEBTotE(62000,0.0);
   
   if(mctruthMatchingInfo_){
      
     //getting the PCAloHit
     edm::Handle<edm::PCaloHitContainer> pcalohits;
     //   bool found_phit 
     //     = iEvent.getByLabel("famosSimHits","EcalHitsEB",
     //          pcalohits);  
     //modif-beg
     bool found_phit 
       = iEvent.getByToken(tokenFamosSimHits_,pcalohits);
     //modif-end
     
     if(!found_phit) {
       ostringstream err;
       err<<"could not find pcaloHit "<<"famosSimHits:EcalHitsEB";
       LogError("PFSimParticleProducer")<<err.str()<<endl;
       
       throw cms::Exception( "MissingProduct", err.str());
     }
     else {
       assert( pcalohits.isValid() );
       
       //     cout << "PFSimParticleProducer: number of pcalohits="
       //     << pcalohits->size() << endl;
       
       edm::PCaloHitContainer::const_iterator it    
     = pcalohits.product()->begin();
       edm::PCaloHitContainer::const_iterator itend 
     = pcalohits.product()->end();
       
       //loop on the PCaloHit from FastSim Calorimetry
       for(;it!=itend;++it)
     {
       EBDetId detid(it->id());
       
       //    cout << detid << " " << detid.rawId()
       //         << " " <<  detid.hashedIndex() 
       //         << " " << it->energy() 
       //         << " " << it->id() 
       //         << " trackId=" 
       //         << it->geantTrackId() 
       //         << endl;
       
       if(it->energy() > 0.0) {
         std::pair<double, unsigned> phitsimid
           = make_pair(it->energy(),it->geantTrackId());
         caloHitsEBID[detid.hashedIndex()].push_back(phitsimid);
         caloHitsEBTotE[detid.hashedIndex()] 
           += it->energy(); //summing pcalhit energy   
       }//energy > 0
       
     }//loop PcaloHits    
     }//pcalohit handle access
     
     //Retrieving the PFRecTrack collection for
     //Monte Carlo Truth Matching tool
     Handle< reco::PFRecTrackCollection > recTracks;
     try{      
       LogDebug("PFSimParticleProducer")<<"getting PFRecTracks"<<endl;
       iEvent.getByToken(tokenRecTracks_, recTracks);
       
     } catch (cms::Exception& err) { 
       LogError("PFSimParticleProducer")<<err
                        <<" cannot get collection "
                        <<"particleFlowBlock"<<":"
                        <<""
                        <<endl;
     }//pfrectrack handle access
     
     //getting the simID corresponding to 
     //each of the PFRecTracks
     getSimIDs( recTracks, recTrackSimID );
     
   }//mctruthMatchingInfo_ //modif
   
   // deal with true particles 
   if( processParticles_) {
 
     auto pOutputPFSimParticleCollection = std::make_unique<reco::PFSimParticleCollection>(); 
 
     Handle<vector<SimTrack> > simTracks;
     bool found = iEvent.getByToken(tokenSim_,simTracks);
     if(!found) {
 
       ostringstream err;
       err<<"cannot find sim tracks: "<<inputTagSim_;
       LogError("PFSimParticleProducer")<<err.str()<<endl;
       
       throw cms::Exception( "MissingProduct", err.str());
     }
       
     
     
     Handle<vector<SimVertex> > simVertices;
     found = iEvent.getByToken(tokenSimVertices_,simVertices);
     if(!found) {
       LogError("PFSimParticleProducer")
     <<"cannot find sim vertices: "<<inputTagSim_<<endl;
       return;
     }
 
       
 
     //     for(unsigned it = 0; it<simTracks->size(); it++ ) {
     //       cout<<"\t track "<< (*simTracks)[it]<<" "
     //    <<(*simTracks)[it].momentum().vect().perp()<<" "
     //    <<(*simTracks)[it].momentum().e()<<endl;
     //     }
 
     mySimEvent->fill( *simTracks, *simVertices );
       
     if(verbose_) 
       mySimEvent->print();
 
     // const std::vector<FSimTrack>& fsimTracks = *(mySimEvent->tracks() );
     for(unsigned i=0; i<mySimEvent->nTracks(); i++) {
     
       const FSimTrack& fst = mySimEvent->track(i);
 
       int motherId = -1;
       if( ! fst.noMother() ) // a mother exist
     motherId = fst.mother().id();
 
       //This is finding out the simID corresponding 
       //to the recTrack
 //       cout << "Particle " << i 
 //     << " " << fst.genpartIndex() 
 //     << " -------------------------------------" << endl;
 
       //GETTING THE TRACK ID
       unsigned         recTrackID = 99999;
       vector<unsigned> recHitContrib; //modif
       vector<double>   recHitContribFrac; //modif
 
       if(mctruthMatchingInfo_){ //modif
 
     for(unsigned lo=0; lo<recTrackSimID.size(); 
         lo++) {
       if( i == recTrackSimID[lo] ) {
 //      cout << "Corresponding Rec Track " 
 //       << lo << endl;
         recTrackID = lo;
       }//match track
     }//loop rectrack
 //  if( recTrackID == 99999 ) 
 //    cout << "Sim Track not reconstructed pT=" <<  
 //      fst.momentum().pt() << endl;
     
     // get the ecalBarrel rechits for MC truth matching tool
     edm::Handle<EcalRecHitCollection> rhcHandle;
     bool found = iEvent.getByToken(tokenEcalRecHitsEB_, 
                        rhcHandle);
     if(!found) {
       ostringstream err;
       err<<"could not find rechits "<< inputTagEcalRecHitsEB_;
       LogError("PFSimParticleProducer")<<err.str()<<endl;
       
       throw cms::Exception( "MissingProduct", err.str());
     }
     else {
       assert( rhcHandle.isValid() );
 //    cout << "PFSimParticleProducer: number of rechits="
 //         << rhcHandle->size() << endl;
       
       EBRecHitCollection::const_iterator it_rh    
         = rhcHandle.product()->begin();
       EBRecHitCollection::const_iterator itend_rh 
         = rhcHandle.product()->end();
       
       for(;it_rh!=itend_rh;++it_rh)
         {
           unsigned rhit_hi 
         = EBDetId(it_rh->id()).hashedIndex();
           EBDetId detid(it_rh->id());
 //      cout << detid << " " << detid.rawId()
 //       << " " <<  detid.hashedIndex() 
 //       << " " << it_rh->energy() << endl;
         
           ITM it_phit    = caloHitsEBID[rhit_hi].begin();
           ITM itend_phit = caloHitsEBID[rhit_hi].end();    
           for(;it_phit!=itend_phit;++it_phit)
         {
           if(i == it_phit->second)
             {
               //Alex (08/10/08) TO BE REMOVED, eliminating
               //duplicated rechits
               bool alreadyin = false;
               for( unsigned ihit = 0; ihit < recHitContrib.size(); 
                ++ihit )
             if(detid.rawId() == recHitContrib[ihit]) 
               alreadyin = true;
               
               if(!alreadyin){       
             double pcalofraction = 0.0;
             if(caloHitsEBTotE[rhit_hi] != 0.0)
               pcalofraction 
                 = (it_phit->first/caloHitsEBTotE[rhit_hi])*100.0;
             
             //store info
             recHitContrib.push_back(it_rh->id());
             recHitContribFrac.push_back(pcalofraction);
               }//selected rechits       
             }//matching
         }//loop pcalohit
           
         }//loop rechits
       
     }//getting the rechits
 
       }//mctruthMatchingInfo_ //modif
 
 //       cout << "This particule has " << recHitContrib.size() 
 //     << " rechit contribution" << endl;
 //       for( unsigned ih = 0; ih < recHitContrib.size(); ++ih )
 //  cout << recHitContrib[ih] 
 //       << " f=" << recHitContribFrac[ih] << " ";
 //       cout << endl;
 
       reco::PFSimParticle particle(  fst.charge(), 
                      fst.type(), 
                      fst.id(), 
                      motherId,
                      fst.daughters(),
                      recTrackID,
                      recHitContrib,
                      recHitContribFrac);
 
 
       const FSimVertex& originVtx = fst.vertex();
 
       math::XYZPoint          posOrig( originVtx.position().x(), 
                        originVtx.position().y(), 
                        originVtx.position().z() );
 
       math::XYZTLorentzVector momOrig( fst.momentum().px(), 
                        fst.momentum().py(), 
                        fst.momentum().pz(), 
                        fst.momentum().e() );
       reco::PFTrajectoryPoint 
     pointOrig(-1, 
           reco::PFTrajectoryPoint::ClosestApproach,
           posOrig, momOrig);
 
       // point 0 is origin vertex
       particle.addPoint(pointOrig);
     
 
       if( ! fst.noEndVertex() ) {
     const FSimVertex& endVtx = fst.endVertex();
     
     math::XYZPoint          posEnd( endVtx.position().x(), 
                     endVtx.position().y(), 
                     endVtx.position().z() );
     //       cout<<"end vertex : "
     //    <<endVtx.position().x()<<" "
     //    <<endVtx.position().y()<<endl;
     
     math::XYZTLorentzVector momEnd;
     
     reco::PFTrajectoryPoint 
       pointEnd( -1, 
             reco::PFTrajectoryPoint::BeamPipeOrEndVertex,
             posEnd, momEnd);
     
     particle.addPoint(pointEnd);
       }
       else { // add a dummy point
     reco::PFTrajectoryPoint dummy;
     particle.addPoint(dummy); 
       }
 
 
       if( fst.onLayer1() ) { // PS layer1
     const RawParticle& rp = fst.layer1Entrance();
       
     math::XYZPoint posLayer1( rp.x(), rp.y(), rp.z() );
     math::XYZTLorentzVector momLayer1( rp.px(), rp.py(), rp.pz(), rp.e() );
     reco::PFTrajectoryPoint layer1Pt(-1, reco::PFTrajectoryPoint::PS1, 
                      posLayer1, momLayer1);
     
     particle.addPoint( layer1Pt ); 
 
     // extrapolate to cluster depth
       }
       else { // add a dummy point
     reco::PFTrajectoryPoint dummy;
     particle.addPoint(dummy); 
       }
 
       if( fst.onLayer2() ) { // PS layer2
     const RawParticle& rp = fst.layer2Entrance();
       
     math::XYZPoint posLayer2( rp.x(), rp.y(), rp.z() );
     math::XYZTLorentzVector momLayer2( rp.px(), rp.py(), rp.pz(), rp.e() );
     reco::PFTrajectoryPoint layer2Pt(-1, reco::PFTrajectoryPoint::PS2, 
                      posLayer2, momLayer2);
     
     particle.addPoint( layer2Pt ); 
 
     // extrapolate to cluster depth
       }
       else { // add a dummy point
     reco::PFTrajectoryPoint dummy;
     particle.addPoint(dummy); 
       }
 
       if( fst.onEcal() ) {
     const RawParticle& rp = fst.ecalEntrance();
     
     math::XYZPoint posECAL( rp.x(), rp.y(), rp.z() );
     math::XYZTLorentzVector momECAL( rp.px(), rp.py(), rp.pz(), rp.e() );
     reco::PFTrajectoryPoint ecalPt(-1, 
                        reco::PFTrajectoryPoint::ECALEntrance, 
                        posECAL, momECAL);
     
     particle.addPoint( ecalPt ); 
 
     // extrapolate to cluster depth
       }
       else { // add a dummy point
     reco::PFTrajectoryPoint dummy;
     particle.addPoint(dummy); 
       }
     
       // add a dummy point for ECAL Shower max
       reco::PFTrajectoryPoint dummy;
       particle.addPoint(dummy); 
 
       if( fst.onHcal() ) {
 
     const RawParticle& rpin = fst.hcalEntrance();
     
     math::XYZPoint posHCALin( rpin.x(), rpin.y(), rpin.z() );
     math::XYZTLorentzVector momHCALin( rpin.px(), rpin.py(), rpin.pz(), 
                        rpin.e() );
     reco::PFTrajectoryPoint hcalPtin(-1, 
                      reco::PFTrajectoryPoint::HCALEntrance,
                      posHCALin, momHCALin);
     
     particle.addPoint( hcalPtin ); 
 
     const RawParticle& rpout = fst.hcalExit();
     
     math::XYZPoint posHCALout( rpout.x(), rpout.y(), rpout.z() );
     math::XYZTLorentzVector momHCALout( rpout.px(), rpout.py(), rpout.pz(),
                         rpout.e() );
     reco::PFTrajectoryPoint 
       hcalPtout(-1, reco::PFTrajectoryPoint::HCALExit, 
              posHCALout, momHCALout);
     
     particle.addPoint( hcalPtout );     
 
     const RawParticle& rpho = fst.hoEntrance();
     
     math::XYZPoint posHOEntrance( rpho.x(), rpho.y(), rpho.z() );
     math::XYZTLorentzVector momHOEntrance( rpho.px(), rpho.py(), rpho.pz(),
                         rpho.e() );
     reco::PFTrajectoryPoint 
       hoPtin(-1, reco::PFTrajectoryPoint::HOLayer, 
          posHOEntrance, momHOEntrance);
     
     particle.addPoint( hoPtin );    
 
 
 
 
       }
       else { // add a dummy point
     reco::PFTrajectoryPoint dummy;
     particle.addPoint(dummy); 
       }
           
       pOutputPFSimParticleCollection->push_back( particle );
     }
     
     iEvent.put(std::move(pOutputPFSimParticleCollection));
   }
 
   
   LogDebug("PFSimParticleProducer")<<"STOP event: "<<iEvent.id().event()
                    <<" in run "<<iEvent.id().run()<<endl;
 }

Member Data Documentation

edm::InputTag PFSimParticleProducer::inputTagEcalRecHitsEB_

private

Definition at line 74 of file PFSimParticleProducer.h.

edm::InputTag PFSimParticleProducer::inputTagEcalRecHitsEE_

private

Definition at line 76 of file PFSimParticleProducer.h.

edm::InputTag PFSimParticleProducer::inputTagFamosSimHits_

private

Definition at line 68 of file PFSimParticleProducer.h.

edm::InputTag PFSimParticleProducer::inputTagRecTracks_

private

Definition at line 72 of file PFSimParticleProducer.h.

edm::InputTag PFSimParticleProducer::inputTagSim_

private

module label for retrieving input simtrack and simvertex

Definition at line 61 of file PFSimParticleProducer.h.

bool PFSimParticleProducer::mctruthMatchingInfo_

private

Definition at line 67 of file PFSimParticleProducer.h.

FSimEvent* PFSimParticleProducer::mySimEvent

private

Definition at line 82 of file PFSimParticleProducer.h.

edm::ParameterSet PFSimParticleProducer::particleFilter_

private

Definition at line 81 of file PFSimParticleProducer.h.

bool PFSimParticleProducer::processParticles_

private

process particles on/off

Definition at line 87 of file PFSimParticleProducer.h.

edm::EDGetTokenT<EcalRecHitCollection> PFSimParticleProducer::tokenEcalRecHitsEB_

private

Definition at line 75 of file PFSimParticleProducer.h.

edm::EDGetTokenT<EcalRecHitCollection> PFSimParticleProducer::tokenEcalRecHitsEE_

private

Definition at line 77 of file PFSimParticleProducer.h.

edm::EDGetTokenT<edm::PCaloHitContainer> PFSimParticleProducer::tokenFamosSimHits_

private

Definition at line 69 of file PFSimParticleProducer.h.

edm::EDGetTokenT<reco::PFRecTrackCollection> PFSimParticleProducer::tokenRecTracks_

private

Definition at line 73 of file PFSimParticleProducer.h.

edm::EDGetTokenT<std::vector<SimTrack> > PFSimParticleProducer::tokenSim_

private

Definition at line 62 of file PFSimParticleProducer.h.

edm::EDGetTokenT<std::vector<SimVertex> > PFSimParticleProducer::tokenSimVertices_

private

Definition at line 63 of file PFSimParticleProducer.h.

bool PFSimParticleProducer::verbose_

private

verbose ?

Definition at line 90 of file PFSimParticleProducer.h.

Public Types

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation