#include <PFNuclearProducer.h>

Inheritance diagram for PFNuclearProducer:

Public Member Functions
	PFNuclearProducer (const edm::ParameterSet &)
	Constructor. More...

	~PFNuclearProducer () override
	Destructor. More...

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Private Member Functions
void	beginRun (const edm::Run &, const edm::EventSetup &) override

void	endRun (const edm::Run &, const edm::EventSetup &) override

void	produce (edm::Event &, const edm::EventSetup &) override
	Produce the PFRecTrack collection. More...

Private Attributes
double	likelihoodCut_

std::vector< edm::EDGetTokenT< reco::NuclearInteractionCollection > >	nuclearContainers_

PFTrackTransformer *	pfTransformer_
	PFTrackTransformer. More...

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

Definition at line 12 of file PFNuclearProducer.h.

Constructor & Destructor Documentation

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

explicit

Constructor.

Definition at line 10 of file PFNuclearProducer.cc.

References edm::ParameterSet::getParameter(), mps_fire::i, likelihoodCut_, and nuclearContainers_.

                                                                :
   pfTransformer_(nullptr)
 {
   produces<reco::PFRecTrackCollection>();
   produces<reco::PFNuclearInteractionCollection>();
 
   std::vector<edm::InputTag> tags= 
     iConfig.getParameter< vector < InputTag > >("nuclearColList");
 
   for (unsigned int i=0;i<tags.size();++i) 
     nuclearContainers_.push_back(consumes<reco::NuclearInteractionCollection>(tags[i]));
 
   likelihoodCut_
      = iConfig.getParameter<double>("likelihoodCut");
 }

PFNuclearProducer::~PFNuclearProducer ( )

override

Destructor.

Definition at line 26 of file PFNuclearProducer.cc.

References pfTransformer_.

 {
   delete pfTransformer_;
 }

Member Function Documentation

void PFNuclearProducer::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 79 of file PFNuclearProducer.cc.

References edm::EventSetup::get(), MagneticField::inTesla(), seedCreatorFromRegionConsecutiveHitsEDProducer_cff::magneticField, PFTrackTransformer::OnlyProp(), and pfTransformer_.

 {
   ESHandle<MagneticField> magneticField;
   iSetup.get<IdealMagneticFieldRecord>().get(magneticField);
   pfTransformer_= new PFTrackTransformer(math::XYZVector(magneticField->inTesla(GlobalPoint(0,0,0))));
   pfTransformer_->OnlyProp();
 }

void PFNuclearProducer::endRun	(	const edm::Run &	run,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 90 of file PFNuclearProducer.cc.

References pfTransformer_.

                                         {
   delete pfTransformer_;
   pfTransformer_=nullptr;
 }

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

overrideprivate

Produce the PFRecTrack collection.

Definition at line 32 of file PFNuclearProducer.cc.

References PFTrackTransformer::addPoints(), edm::Event::getByToken(), edm::Event::getRefBeforePut(), icoll, reco::PFRecTrack::KF, likelihood(), likelihoodCut_, eostools::move(), nuclearContainers_, pfTransformer_, edm::Handle< T >::product(), edm::RefVector< C, T, F >::push_back(), and edm::Event::put().

 {
   typedef reco::NuclearInteraction::trackRef_iterator trackRef_iterator;
 
   //create the empty collections 
   auto pfNuclearColl = std::make_unique<reco::PFNuclearInteractionCollection>();
   auto pfNuclearRecTrackColl = std::make_unique<reco::PFRecTrackCollection>();
   
   reco::PFRecTrackRefProd pfTrackRefProd = iEvent.getRefBeforePut<reco::PFRecTrackCollection>();
   int hid=0;
 
   // loop on the nuclear interaction collections
   for (unsigned int istr=0; istr<nuclearContainers_.size();istr++){
     
     Handle<reco::NuclearInteractionCollection> nuclCollH;
     iEvent.getByToken(nuclearContainers_[istr], nuclCollH);
     const reco::NuclearInteractionCollection& nuclColl = *(nuclCollH.product());
 
     // loop on all NuclearInteraction 
     for( unsigned int icoll=0; icoll < nuclColl.size(); icoll++) {
 
       if( nuclColl[icoll].likelihood() < likelihoodCut_) continue;
 
       reco::PFRecTrackRefVector pfRecTkcoll;
 
       // convert the secondary tracks
       for(trackRef_iterator it = nuclColl[icoll].secondaryTracks_begin(); it!=nuclColl[icoll].secondaryTracks_end(); it++){
         reco::PFRecTrack pftrack( (*it)->charge(), 
                     reco::PFRecTrack::KF, 
                     it->key(), (reco::TrackRef)((*it).castTo<reco::TrackRef>()) );
         Trajectory FakeTraj;
         bool valid = pfTransformer_->addPoints( pftrack, **it, FakeTraj);
         if(valid) {
       pfRecTkcoll.push_back(reco::PFRecTrackRef( pfTrackRefProd, hid++ ));  
           pfNuclearRecTrackColl->push_back(pftrack);
         }
       }
       reco::NuclearInteractionRef niRef(nuclCollH, icoll);
       pfNuclearColl->push_back( reco::PFNuclearInteraction( niRef, pfRecTkcoll ));
     }
   }
   iEvent.put(std::move(pfNuclearRecTrackColl));
   iEvent.put(std::move(pfNuclearColl));
 }