#include <FSPFProducer.h>

Inheritance diagram for FSPFProducer:

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

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

	~FSPFProducer ()

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
	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
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
double	energy_threshold (double eta)

Private Attributes
double	barrel_th

std::vector< double >	EM_HF_ScaleFactor

double	endcap_th

double	HF_Ratio

edm::InputTag	labelPFCandidateCollection_

double	middle_th

double	par1

double	par2

edm::EDGetTokenT < reco::PFCandidateCollection >	pfCandidateToken

bool	pfPatchInHF

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

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

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

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::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

Definition at line 19 of file FSPFProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 13 of file FSPFProducer.cc.

References edm::ParameterSet::getParameter().

                                                          {
  
   labelPFCandidateCollection_ = iConfig.getParameter < edm::InputTag > ("pfCandidates");
 
   pfPatchInHF = iConfig.getParameter<bool>("pfPatchInHF");
   EM_HF_ScaleFactor = iConfig.getParameter< std::vector <double> >("EM_HF_ScaleFactor");
   HF_Ratio   = iConfig.getParameter<double>("HF_Ratio");
   par1       = iConfig.getParameter<double>("par1");
   par2       = iConfig.getParameter<double>("par2");
   barrel_th  = iConfig.getParameter<double>("barrel_th");
   endcap_th  = iConfig.getParameter<double>("endcap_th");
   middle_th  = iConfig.getParameter<double>("middle_th");
   // register products
   produces<reco::PFCandidateCollection>();
 
   // consumes
   pfCandidateToken = consumes<reco::PFCandidateCollection>(labelPFCandidateCollection_);
 }

FSPFProducer::~FSPFProducer ( )

Definition at line 32 of file FSPFProducer.cc.

32 {}

Member Function Documentation

double FSPFProducer::energy_threshold ( double eta )

private

Definition at line 109 of file FSPFProducer.cc.

References eta().

                                                 {
   if (eta<0) eta = -eta;
   if      (eta < 1.6) return barrel_th;
   else if (eta < 1.8) return middle_th;
   else                return endcap_th;
 }

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

overridevirtual

Implements edm::stream::EDProducerBase.

Definition at line 35 of file FSPFProducer.cc.

References reco::PFCandidate::egamma_HF, edm::Event::getByToken(), reco::PFCandidate::h, reco::PFCandidate::h0, reco::PFCandidate::h_HF, reco::tau::pfCandidates(), edm::Event::put(), and mathSSE::sqrt().

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

                                         {
   
   Handle < reco::PFCandidateCollection > pfCandidates;
   iEvent.getByToken (pfCandidateToken, pfCandidates);
   
   auto_ptr< reco::PFCandidateCollection >  pOutputCandidateCollection(new PFCandidateCollection);   
 
   /*  
   LogDebug("FSPFProducer")<<"START event: "
               <<iEvent.id().event()
               <<" in run "<<iEvent.id().run()<<endl;   
   */  
 
   double theNeutralFraction, px, py, pz, en;
   int n_hadron_HF = 0; 
   int n_em_HF = 0;
   int vEta=-1;
   reco::PFCandidateCollection::const_iterator  itCand =  pfCandidates->begin();
   reco::PFCandidateCollection::const_iterator  itCandEnd = pfCandidates->end();
   for( ; itCand != itCandEnd; itCand++) {
 
     // First part: create fake neutral hadrons as a function of charged hadrons and add them to the new collection
     if(itCand->particleId() == reco::PFCandidate::h){
       theNeutralFraction = par1 - par2*itCand->energy();
       if(theNeutralFraction > 0.){
     px = theNeutralFraction*itCand->px();
     py = theNeutralFraction*itCand->py();
     pz = theNeutralFraction*itCand->pz();
     en = sqrt(px*px + py*py + pz*pz); 
     if (en > energy_threshold(itCand->eta())) {
       // create a PFCandidate and add it to the particles Collection
       math::XYZTLorentzVector momentum(px,py,pz,en);
       reco::PFCandidate FakeNeutralHadron(0, momentum,  reco::PFCandidate::h0);
       pOutputCandidateCollection->push_back(FakeNeutralHadron);
     }
       }
     }
     
     // Second part: deal with HF, and put every candidate of the old collection in the new collection
     if(itCand->particleId() == reco::PFCandidate::egamma_HF){
       if (pfPatchInHF) {
     n_em_HF++;
     
     if(fabs(itCand->eta())< 4.) vEta = 0;
     else if(fabs(itCand->eta())<= 5.) vEta = 1;
     
     if (vEta==0 || vEta==1) {
       // copy these PFCandidates after the momentum rescaling
       px = EM_HF_ScaleFactor[vEta]*itCand->px();
       py = EM_HF_ScaleFactor[vEta]*itCand->py();
       pz = EM_HF_ScaleFactor[vEta]*itCand->pz();
       en = sqrt(px*px + py*py + pz*pz); 
       math::XYZTLorentzVector momentum(px,py,pz,en);
       reco::PFCandidate EMHF(itCand->charge(), momentum,  reco::PFCandidate::egamma_HF);
       if(en>0.) pOutputCandidateCollection->push_back(EMHF);
     }
       } else pOutputCandidateCollection->push_back(*itCand);
     }
     
     else if(itCand->particleId() == reco::PFCandidate::h_HF){
       if (pfPatchInHF) {
     // copy these PFCandidates to the new particles Collection only if hadron candidates are currently more than EM candidates
     n_hadron_HF++;
     if(n_em_HF < (n_hadron_HF*HF_Ratio)) pOutputCandidateCollection->push_back(*itCand);
       } else pOutputCandidateCollection->push_back(*itCand);
     }
 
     else  pOutputCandidateCollection->push_back(*itCand);
     
   }
   iEvent.put(pOutputCandidateCollection);
 }

Member Data Documentation

double FSPFProducer::barrel_th

private

Definition at line 31 of file FSPFProducer.h.

std::vector<double> FSPFProducer::EM_HF_ScaleFactor

private

Definition at line 35 of file FSPFProducer.h.

double FSPFProducer::endcap_th

private

Definition at line 31 of file FSPFProducer.h.

double FSPFProducer::HF_Ratio

private

Definition at line 34 of file FSPFProducer.h.

edm::InputTag FSPFProducer::labelPFCandidateCollection_

private

Definition at line 28 of file FSPFProducer.h.

double FSPFProducer::middle_th

private

Definition at line 31 of file FSPFProducer.h.

double FSPFProducer::par1

private

Definition at line 30 of file FSPFProducer.h.

double FSPFProducer::par2

private

Definition at line 30 of file FSPFProducer.h.

edm::EDGetTokenT<reco::PFCandidateCollection> FSPFProducer::pfCandidateToken

private

Definition at line 39 of file FSPFProducer.h.

bool FSPFProducer::pfPatchInHF

private

Definition at line 33 of file FSPFProducer.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation