Producer for particle flow blocks. More...

#include <PFBlockProducer.h>

Inheritance diagram for PFBlockProducer:

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

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

	~PFBlockProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
edm::InputTag	inputTagConvBremGsfRecTracks_

edm::InputTag	inputTagEGPhotons_

edm::InputTag	inputTagGsfRecTracks_

edm::InputTag	inputTagPFClustersECAL_

edm::InputTag	inputTagPFClustersHCAL_

edm::InputTag	inputTagPFClustersHFEM_

edm::InputTag	inputTagPFClustersHFHAD_

edm::InputTag	inputTagPFClustersHO_

edm::InputTag	inputTagPFClustersPS_

edm::InputTag	inputTagPFConversions_

edm::InputTag	inputTagPFNuclear_

edm::InputTag	inputTagPFV0_

edm::InputTag	inputTagRecMuons_

edm::InputTag	inputTagRecTracks_

edm::InputTag	inputTagSCBarrel_

edm::InputTag	inputTagSCEndcap_

PFBlockAlgo	pfBlockAlgo_
	Particle flow block algorithm. More...

bool	useConvBremGsfTracks_
	switch on/off Conversions Brem Recovery More...

bool	useConversions_
	switch on/off Conversions More...

bool	useEGPhotons_
	use EG photons ? More...

bool	useHO_

bool	useKDTreeTrackEcalLinker_

bool	useNuclear_
	use NuclearInteractions ? More...

bool	usePFatHLT_
	Particle Flow at HLT ? More...

bool	useSuperClusters_
	use SuperClusters ? More...

bool	useV0_
	switch on/off V0 More...

bool	verbose_
	verbose ? More...

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

typedef WorkerT< EDProducer >	WorkerType

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

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const

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

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 particle flow blocks.

This producer makes use of PFBlockAlgo, the particle flow block algorithm. Particle flow itself consists in reconstructing particles from the particle flow blocks This is done at a later stage, see PFProducer and PFAlgo.

Author: Colin Bernet

Date: April 2007

Definition at line 34 of file PFBlockProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 34 of file PFBlockProducer.cc.

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

                                                                {
   
 
   // use configuration file to setup input/output collection names
   inputTagRecTracks_ 
     = iConfig.getParameter<InputTag>("RecTracks");
 
   inputTagGsfRecTracks_ 
     = iConfig.getParameter<InputTag>("GsfRecTracks");
 
     inputTagConvBremGsfRecTracks_ 
     = iConfig.getParameter<InputTag>("ConvBremGsfRecTracks");
 
   inputTagRecMuons_ 
     = iConfig.getParameter<InputTag>("RecMuons");
 
   inputTagPFNuclear_ 
     = iConfig.getParameter<InputTag>("PFNuclear");
 
   inputTagPFConversions_ 
     = iConfig.getParameter<InputTag>("PFConversions");
 
   inputTagPFV0_ 
     = iConfig.getParameter<InputTag>("PFV0");
 
   inputTagPFClustersECAL_ 
     = iConfig.getParameter<InputTag>("PFClustersECAL");
 
   inputTagPFClustersHCAL_ 
     = iConfig.getParameter<InputTag>("PFClustersHCAL");
 
   inputTagPFClustersHO_ 
     = iConfig.getParameter<InputTag>("PFClustersHO");
 
   inputTagPFClustersHFEM_ 
     = iConfig.getParameter<InputTag>("PFClustersHFEM");
 
   inputTagPFClustersHFHAD_ 
     = iConfig.getParameter<InputTag>("PFClustersHFHAD");
 
   inputTagPFClustersPS_ 
     = iConfig.getParameter<InputTag>("PFClustersPS");
 
   useEGPhotons_ = iConfig.getParameter<bool>("useEGPhotons");
   
   if(useEGPhotons_) {
     inputTagEGPhotons_
       = iConfig.getParameter<InputTag>("EGPhotons");         
   }
   
   useSuperClusters_ = iConfig.existsAs<bool>("useSuperClusters") ? iConfig.getParameter<bool>("useSuperClusters") : false;
   
   if (useSuperClusters_) {
     inputTagSCBarrel_
       = iConfig.getParameter<InputTag>("SCBarrel");      
     inputTagSCEndcap_
       = iConfig.getParameter<InputTag>("SCEndcap");     
   }
 
   verbose_ = 
     iConfig.getUntrackedParameter<bool>("verbose",false);
 
   bool debug_ = 
     iConfig.getUntrackedParameter<bool>("debug",false);
 
   usePFatHLT_ = iConfig.getParameter<bool>("usePFatHLT");
 
   useNuclear_ = iConfig.getParameter<bool>("useNuclear");
 
   useConversions_ = iConfig.getParameter<bool>("useConversions");
   
   useConvBremGsfTracks_ = iConfig.getParameter<bool>("useConvBremGsfTracks");
 
   bool useConvBremPFRecTracks = iConfig.getParameter<bool>("useConvBremPFRecTracks");
 
   useV0_ = iConfig.getParameter<bool>("useV0");
 
   useHO_=  iConfig.getParameter<bool>("useHO");
 
   produces<reco::PFBlockCollection>();
   
   // Glowinski & Gouzevitch
   useKDTreeTrackEcalLinker_ = iConfig.getParameter<bool>("useKDTreeTrackEcalLinker");
   // !Glowinski & Gouzevitch
   
   // particle flow parameters  -----------------------------------
 
   std::vector<double> DPtovPtCut 
      = iConfig.getParameter<std::vector<double> >("pf_DPtoverPt_Cut");   
   if (DPtovPtCut.size()!=5)
     {
       edm::LogError("MisConfiguration")<<" vector pf_DPtoverPt_Cut has to be of size 5";
       throw;
     }
 
   std::vector<unsigned> NHitCut 
      = iConfig.getParameter<std::vector<unsigned> >("pf_NHit_Cut");   
   if (NHitCut.size()!=5)
     {
       edm::LogError("MisConfiguration")<<" vector pf_NHit_Cut has to be of size 5";
       throw;
     }
 
   bool useIterTracking
     = iConfig.getParameter<bool>("useIterTracking");
 
   int nuclearInteractionsPurity
     = iConfig.getParameter<unsigned>("nuclearInteractionsPurity");
 
   // if first parameter 0, deactivated 
   std::vector<double> EGPhotonSelectionCuts ;
 
   if (useEGPhotons_)
     EGPhotonSelectionCuts = iConfig.getParameter<std::vector<double> >("PhotonSelectionCuts");   
 
 
   if (useNuclear_){
     if (nuclearInteractionsPurity > 3 || nuclearInteractionsPurity < 1)  {
       nuclearInteractionsPurity = 1;
       edm::LogInfo("PFBlockProducer")  << "NI purity not properly implemented. Set it to the strongest level " << nuclearInteractionsPurity << endl;
     }
     vector<string> securityLevel;
     securityLevel.push_back("isNucl"); securityLevel.push_back("isNucl && isNuclLoose");  securityLevel.push_back("isNucl && isNuclLoose && isNuclKink"); 
     edm::LogInfo("PFBlockProducer")  << "NI interactions are corrected in PFlow for " << securityLevel[nuclearInteractionsPurity-1].c_str() << endl;
   }
 
 
   pfBlockAlgo_.setParameters( DPtovPtCut,
                   NHitCut,
                   useConvBremPFRecTracks,
                   useIterTracking,
                   nuclearInteractionsPurity,
                   useEGPhotons_,
                   EGPhotonSelectionCuts,
                   useSuperClusters_
                 );
   
   pfBlockAlgo_.setDebug(debug_);
 
   // Glowinski & Gouzevitch
   pfBlockAlgo_.setUseOptimization(useKDTreeTrackEcalLinker_);
   // !Glowinski & Gouzevitch
 
   // Use HO clusters for link
   pfBlockAlgo_.setHOTag(useHO_);
 
 }

PFBlockProducer::~PFBlockProducer ( )

Definition at line 184 of file PFBlockProducer.cc.

184 { }

Member Function Documentation

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

overridevirtual

Implements edm::EDProducer.

Definition at line 189 of file PFBlockProducer.cc.

References edm::EventID::event(), newFWLiteAna::found, edm::Event::getByLabel(), edm::EventBase::id(), LogDebug, pfConversions_cfi::pfConversions, pfV0_cfi::pfV0, edm::Event::put(), and edm::EventID::run().

                                        {
   
   LogDebug("PFBlockProducer")<<"START event: "<<iEvent.id().event()
                  <<" in run "<<iEvent.id().run()<<endl;
   
   
   // get rectracks
   
   Handle< reco::PFRecTrackCollection > recTracks;
   
   // LogDebug("PFBlockProducer")<<"get reco tracks"<<endl;
   bool found = iEvent.getByLabel(inputTagRecTracks_, recTracks);
     
   if(!found )
     LogError("PFBlockProducer")<<" cannot get rectracks: "
                    <<inputTagRecTracks_<<endl;
 
 
 
   // get GsfTracks 
   Handle< reco::GsfPFRecTrackCollection > GsfrecTracks;
 
   if(!usePFatHLT_) {
     found = iEvent.getByLabel(inputTagGsfRecTracks_,GsfrecTracks);
 
     if(!found )
       LogError("PFBlockProducer")<<" cannot get Gsfrectracks: "
                  << inputTagGsfRecTracks_ <<endl;
   }
 
   // get ConvBremGsfTracks 
   Handle< reco::GsfPFRecTrackCollection > convBremGsfrecTracks;
 
   if(useConvBremGsfTracks_) {
     found = iEvent.getByLabel(inputTagConvBremGsfRecTracks_,convBremGsfrecTracks);
 
     if(!found )
       LogError("PFBlockProducer")<<" cannot get ConvBremGsfrectracks: "
                  << inputTagConvBremGsfRecTracks_ <<endl;
   }
 
   // get recmuons
   Handle< reco::MuonCollection > recMuons;
 
   // LogDebug("PFBlockProducer")<<"get reco muons"<<endl;
   //if(!usePFatHLT_) {
   found = iEvent.getByLabel(inputTagRecMuons_, recMuons);
   
     //if(!found )
     //  LogError("PFBlockProducer")<<" cannot get recmuons: "
     //           <<inputTagRecMuons_<<endl;
 
   // get PFNuclearInteractions
   //}
   //---------- Gouzevitch
   //  Handle< reco::PFNuclearInteractionCollection > pfNuclears; 
   Handle< reco::PFDisplacedTrackerVertexCollection > pfNuclears; 
 
   if( useNuclear_ ) {
 
     found = iEvent.getByLabel(inputTagPFNuclear_, pfNuclears);
     if(!found )
       LogError("PFBlockProducer")<<" cannot get PFNuclearInteractions : "
                                <<inputTagPFNuclear_<<endl;
   }
  
   // get conversions
   Handle< reco::PFConversionCollection > pfConversions;
   if( useConversions_ ) {
     found = iEvent.getByLabel(inputTagPFConversions_, pfConversions);
     
     if(!found )
       LogError("PFBlockProducer")<<" cannot get PFConversions : "
                  <<inputTagPFConversions_<<endl;
   }
   
 
   // get V0s
   Handle< reco::PFV0Collection > pfV0;
   if( useV0_ ) {
     found = iEvent.getByLabel(inputTagPFV0_, pfV0);
     
     if(!found )
       LogError("PFBlockProducer")<<" cannot get PFV0 : "
                  <<inputTagPFV0_<<endl;
   }
 
 
   
   // get ECAL, HCAL, HO and PS clusters
   
   
   Handle< reco::PFClusterCollection > clustersECAL;
   found = iEvent.getByLabel(inputTagPFClustersECAL_, 
                 clustersECAL);      
   if(!found )
     LogError("PFBlockProducer")<<" cannot get ECAL clusters: "
                    <<inputTagPFClustersECAL_<<endl;
     
   
   Handle< reco::PFClusterCollection > clustersHCAL;
   found = iEvent.getByLabel(inputTagPFClustersHCAL_, 
                 clustersHCAL);      
   if(!found )
     LogError("PFBlockProducer")<<" cannot get HCAL clusters: "
                    <<inputTagPFClustersHCAL_<<endl;
     
   Handle< reco::PFClusterCollection > clustersHO;
   if (useHO_) {
     found = iEvent.getByLabel(inputTagPFClustersHO_, 
                   clustersHO);      
     if(!found )
       LogError("PFBlockProducer")<<" cannot get HO clusters: "
                  <<inputTagPFClustersHO_<<endl;
   }
 
   Handle< reco::PFClusterCollection > clustersHFEM;
   found = iEvent.getByLabel(inputTagPFClustersHFEM_, 
                 clustersHFEM);      
   if(!found )
     LogError("PFBlockProducer")<<" cannot get HFEM clusters: "
                    <<inputTagPFClustersHFEM_<<endl;
     
   Handle< reco::PFClusterCollection > clustersHFHAD;
   found = iEvent.getByLabel(inputTagPFClustersHFHAD_, 
                 clustersHFHAD);      
   if(!found )
     LogError("PFBlockProducer")<<" cannot get HFHAD clusters: "
                    <<inputTagPFClustersHFHAD_<<endl;
     
 
   Handle< reco::PFClusterCollection > clustersPS;
   found = iEvent.getByLabel(inputTagPFClustersPS_, 
                 clustersPS);      
   if(!found )
     LogError("PFBlockProducer")<<" cannot get PS clusters: "
                    <<inputTagPFClustersPS_<<endl;
 
   // dummy. Not used in the full framework 
   Handle< reco::PFRecTrackCollection > nuclearRecTracks;
 
   
   Handle< reco::PhotonCollection >  egPhotons;
   found = iEvent.getByLabel(inputTagEGPhotons_,
                 egPhotons);
 
   if(!found && useEGPhotons_ )
     LogError("PFBlockProducer")<<" cannot get photons" 
                    << inputTagEGPhotons_ << endl;
                    
   Handle< reco::SuperClusterCollection >  sceb;
   Handle< reco::SuperClusterCollection >  scee;
   
   if (useSuperClusters_) {
     found = iEvent.getByLabel(inputTagSCBarrel_,
                   sceb);
 
     if(!found)
       LogError("PFBlockProducer")<<" cannot get sceb" 
                 << inputTagSCBarrel_ << endl;
       
                 
     
     found = iEvent.getByLabel(inputTagSCEndcap_,
                   scee);
 
     if(!found)
       LogError("PFBlockProducer")<<" cannot get scee" 
                 << inputTagSCEndcap_ << endl;                      
                                 
   }
 
   if( usePFatHLT_  ) {
      pfBlockAlgo_.setInput( recTracks,      
                 recMuons,
                 clustersECAL,
                 clustersHCAL,
                 clustersHO,
                 clustersHFEM,
                 clustersHFHAD,
                 clustersPS);
   } else { 
     pfBlockAlgo_.setInput( recTracks, 
                GsfrecTracks,
                convBremGsfrecTracks,
                recMuons, 
                pfNuclears,
                nuclearRecTracks,
                pfConversions,
                pfV0,
                clustersECAL,
                clustersHCAL,
                clustersHO,
                clustersHFEM,
                clustersHFHAD,
                clustersPS,
                egPhotons,
                sceb,
                scee);
   }
   pfBlockAlgo_.findBlocks();
   
   if(verbose_) {
     ostringstream  str;
     str<<pfBlockAlgo_<<endl;
     LogInfo("PFBlockProducer") << str.str()<<endl;
   }    
 
   auto_ptr< reco::PFBlockCollection > 
     pOutputBlockCollection( pfBlockAlgo_.transferBlocks() ); 
   
   iEvent.put(pOutputBlockCollection);
 
   LogDebug("PFBlockProducer")<<"STOP event: "<<iEvent.id().event()
                  <<" in run "<<iEvent.id().run()<<endl;
 }

Member Data Documentation

edm::InputTag PFBlockProducer::inputTagConvBremGsfRecTracks_

private

Definition at line 49 of file PFBlockProducer.h.

edm::InputTag PFBlockProducer::inputTagEGPhotons_

private

Definition at line 60 of file PFBlockProducer.h.

edm::InputTag PFBlockProducer::inputTagGsfRecTracks_

private

Definition at line 48 of file PFBlockProducer.h.

edm::InputTag PFBlockProducer::inputTagPFClustersECAL_

private

Definition at line 52 of file PFBlockProducer.h.

edm::InputTag PFBlockProducer::inputTagPFClustersHCAL_

private

Definition at line 53 of file PFBlockProducer.h.

edm::InputTag PFBlockProducer::inputTagPFClustersHFEM_

private

Definition at line 55 of file PFBlockProducer.h.

edm::InputTag PFBlockProducer::inputTagPFClustersHFHAD_

private

Definition at line 56 of file PFBlockProducer.h.

edm::InputTag PFBlockProducer::inputTagPFClustersHO_

private

Definition at line 54 of file PFBlockProducer.h.

edm::InputTag PFBlockProducer::inputTagPFClustersPS_

private

Definition at line 57 of file PFBlockProducer.h.

edm::InputTag PFBlockProducer::inputTagPFConversions_

private

Definition at line 58 of file PFBlockProducer.h.

edm::InputTag PFBlockProducer::inputTagPFNuclear_

private

Definition at line 51 of file PFBlockProducer.h.

edm::InputTag PFBlockProducer::inputTagPFV0_

private

Definition at line 59 of file PFBlockProducer.h.

edm::InputTag PFBlockProducer::inputTagRecMuons_

private

Definition at line 50 of file PFBlockProducer.h.

edm::InputTag PFBlockProducer::inputTagRecTracks_

private

Definition at line 47 of file PFBlockProducer.h.

edm::InputTag PFBlockProducer::inputTagSCBarrel_

private

Definition at line 61 of file PFBlockProducer.h.

edm::InputTag PFBlockProducer::inputTagSCEndcap_

private

Definition at line 62 of file PFBlockProducer.h.

PFBlockAlgo PFBlockProducer::pfBlockAlgo_

private

Particle flow block algorithm.

Definition at line 97 of file PFBlockProducer.h.

bool PFBlockProducer::useConvBremGsfTracks_

private

switch on/off Conversions Brem Recovery

Definition at line 83 of file PFBlockProducer.h.

bool PFBlockProducer::useConversions_

private

switch on/off Conversions

Definition at line 80 of file PFBlockProducer.h.

bool PFBlockProducer::useEGPhotons_

private

use EG photons ?

Definition at line 74 of file PFBlockProducer.h.

bool PFBlockProducer::useHO_

private

Definition at line 65 of file PFBlockProducer.h.

bool PFBlockProducer::useKDTreeTrackEcalLinker_

private

Definition at line 93 of file PFBlockProducer.h.

bool PFBlockProducer::useNuclear_

private

use NuclearInteractions ?

Definition at line 71 of file PFBlockProducer.h.

bool PFBlockProducer::usePFatHLT_

private

Particle Flow at HLT ?

Definition at line 89 of file PFBlockProducer.h.

bool PFBlockProducer::useSuperClusters_

private

use SuperClusters ?

Definition at line 77 of file PFBlockProducer.h.

bool PFBlockProducer::useV0_

private

switch on/off V0

Definition at line 86 of file PFBlockProducer.h.

bool PFBlockProducer::verbose_

private

verbose ?

Definition at line 68 of file PFBlockProducer.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation