Inheritance diagram for L1TkEmParticleProducer:

Public Types
typedef std::vector< L1TTTrackType >	L1TTTrackCollectionType

typedef TTTrack< Ref_Phase2TrackerDigi_ >	L1TTTrackType

Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	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
float	CorrectedEta (float eta, float zv) const

	L1TkEmParticleProducer (const edm::ParameterSet &)

	~L1TkEmParticleProducer () override

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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

Public Member Functions inherited from edm::global::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

std::vector< bool > const &	recordProvenanceList () 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)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~ProducerBase () noexcept(false) override

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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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 &	operator= (EDConsumerBase &&)=default

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

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
void	produce (edm::StreamID, edm::Event &, const edm::EventSetup &) const override

Private Attributes
float	chi2Max_

float	deltaZMax_

float	dRMax_

float	dRMin_

const edm::EDGetTokenT< EGammaBxCollection >	egToken_

float	etMin_

float	isoCut_

std::string	label_

bool	primaryVtxConstrain_

float	pTMinTra_

bool	relativeIsolation_

const edm::EDGetTokenT< std::vector< TTTrack< Ref_Phase2TrackerDigi_ > > >	trackToken_

const edm::EDGetTokenT< TkPrimaryVertexCollection >	vertexToken_

float	zMax_

Additional Inherited Members
Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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

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

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

Detailed Description

Definition at line 46 of file L1TkEmParticleProducer.cc.

Member Typedef Documentation

◆ L1TTTrackCollectionType

typedef std::vector<L1TTTrackType> L1TkEmParticleProducer::L1TTTrackCollectionType

Definition at line 49 of file L1TkEmParticleProducer.cc.

◆ L1TTTrackType

typedef TTTrack<Ref_Phase2TrackerDigi_> L1TkEmParticleProducer::L1TTTrackType

Definition at line 48 of file L1TkEmParticleProducer.cc.

Constructor & Destructor Documentation

◆ L1TkEmParticleProducer()

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

explicit

Definition at line 87 of file L1TkEmParticleProducer.cc.

     : egToken_(consumes<EGammaBxCollection>(iConfig.getParameter<edm::InputTag>("L1EGammaInputTag"))),
       trackToken_(consumes<std::vector<TTTrack<Ref_Phase2TrackerDigi_> > >(
           iConfig.getParameter<edm::InputTag>("L1TrackInputTag"))),
       vertexToken_(consumes<TkPrimaryVertexCollection>(iConfig.getParameter<edm::InputTag>("L1VertexInputTag"))) {
   label_ = iConfig.getParameter<std::string>("label");  // label of the collection produced
   // e.g. EG or IsoEG if all objects are kept
   // EGIsoTrk or IsoEGIsoTrk if only the EG or IsoEG
   // objects that pass a cut RelIso < isoCut_ are written
   // in the new collection.
  
   etMin_ = (float)iConfig.getParameter<double>("ETmin");
  
   // parameters for the calculation of the isolation :
   zMax_ = (float)iConfig.getParameter<double>("ZMAX");
   chi2Max_ = (float)iConfig.getParameter<double>("CHI2MAX");
   pTMinTra_ = (float)iConfig.getParameter<double>("PTMINTRA");
   dRMin_ = (float)iConfig.getParameter<double>("DRmin");
   dRMax_ = (float)iConfig.getParameter<double>("DRmax");
   primaryVtxConstrain_ = iConfig.getParameter<bool>("PrimaryVtxConstrain");
   deltaZMax_ = (float)iConfig.getParameter<double>("DeltaZMax");
   // cut applied on the isolation (if this number is <= 0, no cut is applied)
   isoCut_ = (float)iConfig.getParameter<double>("IsoCut");
   relativeIsolation_ = iConfig.getParameter<bool>("RelativeIsolation");
  
   produces<TkEmCollection>(label_);
 }

References chi2Max_, deltaZMax_, dRMax_, dRMin_, etMin_, dqmMemoryStats::float, edm::ParameterSet::getParameter(), isoCut_, label_, primaryVtxConstrain_, pTMinTra_, relativeIsolation_, AlCaHLTBitMon_QueryRunRegistry::string, and zMax_.

◆ ~L1TkEmParticleProducer()

L1TkEmParticleProducer::~L1TkEmParticleProducer ( )

override

Definition at line 115 of file L1TkEmParticleProducer.cc.

115 {}

Member Function Documentation

◆ CorrectedEta()

float L1TkEmParticleProducer::CorrectedEta	(	float	eta,
		float	zv
	)		const

Definition at line 253 of file L1TkEmParticleProducer.cc.

                                                                     {
   // Correct the eta of the L1EG object once we know the zvertex
  
   bool IsBarrel = (std::abs(eta) < EtaECal);
  
   float theta = 2. * atan(exp(-eta));
   if (theta < 0)
     theta = theta + M_PI;
   float tantheta = tan(theta);
  
   float delta;
   if (IsBarrel) {
     delta = REcal / tantheta;
   } else {
     if (theta > 0)
       delta = ZEcal;
     if (theta < 0)
       delta = -ZEcal;
   }
  
   float tanthetaprime = delta * tantheta / (delta - zv);
  
   float thetaprime = atan(tanthetaprime);
   if (thetaprime < 0)
     thetaprime = thetaprime + M_PI;
  
   float etaprime = -log(tan(thetaprime / 2.));
   return etaprime;
 }

References funct::abs(), dumpMFGeometry_cfg::delta, PVValHelper::eta, EtaECal, JetChargeProducer_cfi::exp, dqm-mbProfile::log, M_PI, REcal, funct::tan(), theta(), and ZEcal.

Referenced by produce().

◆ fillDescriptions()

void L1TkEmParticleProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 284 of file L1TkEmParticleProducer.cc.

                                                                                         {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

◆ produce()

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

overrideprivatevirtual

Implements edm::global::EDProducerBase.

Definition at line 118 of file L1TkEmParticleProducer.cc.

                                                                                                      {
   using namespace edm;
  
   auto result = std::make_unique<TkEmCollection>();
  
   // the L1EGamma objects
   edm::Handle<EGammaBxCollection> eGammaHandle;
   iEvent.getByToken(egToken_, eGammaHandle);
   EGammaBxCollection eGammaCollection = (*eGammaHandle.product());
   EGammaBxCollection::const_iterator egIter;
  
   // the L1Tracks
   edm::Handle<std::vector<TTTrack<Ref_Phase2TrackerDigi_> > > L1TTTrackHandle;
   iEvent.getByToken(trackToken_, L1TTTrackHandle);
  
   // the primary vertex (used only if primaryVtxConstrain_ = true)
   float zvtxL1tk = -999;
   edm::Handle<TkPrimaryVertexCollection> L1VertexHandle;
   iEvent.getByToken(vertexToken_, L1VertexHandle);
   bool primaryVtxConstrain = primaryVtxConstrain_;
   if (!L1VertexHandle.isValid()) {
     LogWarning("L1TkEmParticleProducer")
         << "Warning: TkPrimaryVertexCollection not found in the event. Won't use any PrimaryVertex constraint."
         << std::endl;
     primaryVtxConstrain = false;
   } else {
     std::vector<TkPrimaryVertex>::const_iterator vtxIter = L1VertexHandle->begin();
     // by convention, the first vertex in the collection is the one that should
     // be used by default
     zvtxL1tk = vtxIter->zvertex();
   }
  
   if (!L1TTTrackHandle.isValid()) {
     LogError("L1TkEmParticleProducer") << "\nWarning: L1TTTrackCollectionType not found in the event. Exit."
                                        << std::endl;
     return;
   }
  
   // Now loop over the L1EGamma objects
  
   if (!eGammaHandle.isValid()) {
     LogError("L1TkEmParticleProducer") << "\nWarning: L1EmCollection not found in the event. Exit." << std::endl;
     return;
   }
  
   int ieg = 0;
   for (egIter = eGammaCollection.begin(0); egIter != eGammaCollection.end(0); ++egIter)  // considering BX = only
   {
     edm::Ref<EGammaBxCollection> EGammaRef(eGammaHandle, ieg);
     ieg++;
  
     float eta = egIter->eta();
     // The eta of the L1EG object is seen from (0,0,0).
     // if primaryVtxConstrain_ = true, and for the PV constrained iso, use the zvtxL1tk to correct the eta(L1EG)
     // that is used in the calculation of DeltaR.
     float etaPV = CorrectedEta((float)eta, zvtxL1tk);
  
     float phi = egIter->phi();
     float et = egIter->et();
  
     if (et < etMin_)
       continue;
  
     // calculate the isolation of the L1EG object with
     // respect to L1Tracks.
  
     float trkisol = -999;
     float sumPt = 0;
     float sumPtPV = 0;
     float trkisolPV = -999;
  
     for (const auto& track : *L1TTTrackHandle) {
       float Pt = track.momentum().perp();
       float Eta = track.momentum().eta();
       float Phi = track.momentum().phi();
       float z = track.POCA().z();
       if (fabs(z) > zMax_)
         continue;
       if (Pt < pTMinTra_)
         continue;
       float chi2 = track.chi2();
       if (chi2 > chi2Max_)
         continue;
  
       float dr = reco::deltaR(Eta, Phi, eta, phi);
       if (dr < dRMax_ && dr >= dRMin_) {
         sumPt += Pt;
       }
  
       if (zvtxL1tk > -999 && std::abs(z - zvtxL1tk) >= deltaZMax_)
         continue;  // Now, PV constrained trackSum:
  
       dr = reco::deltaR(Eta, Phi, etaPV, phi);  // recompute using the corrected eta
  
       if (dr < dRMax_ && dr >= dRMin_) {
         sumPtPV += Pt;
       }
     }  // end loop over tracks
  
     if (relativeIsolation_) {
       if (et > 0) {
         trkisol = sumPt / et;
         trkisolPV = sumPtPV / et;
       }       // relative isolation
     } else {  // absolute isolation
       trkisol = sumPt;
       trkisolPV = sumPtPV;
     }
  
     float isolation = trkisol;
     if (primaryVtxConstrain) {
       isolation = trkisolPV;
     }
  
     const math::XYZTLorentzVector P4 = egIter->p4();
     TkEm trkEm(P4, EGammaRef, trkisol, trkisolPV);
  
     if (isoCut_ <= 0) {
       // write the L1TkEm particle to the collection,
       // irrespective of its relative isolation
       result->push_back(trkEm);
     } else {
       // the object is written to the collection only
       // if it passes the isolation cut
       if (isolation <= isoCut_)
         result->push_back(trkEm);
     }
  
   }  // end loop over EGamma objects
  
   iEvent.put(std::move(result), label_);
 }

References funct::abs(), BXVector< T >::begin(), hltPixelTracks_cff::chi2, chi2Max_, CorrectedEta(), reco::deltaR(), deltaZMax_, flavorHistoryFilter_cfi::dr, dRMin_, egToken_, BXVector< T >::end(), EgHLTOffHistBins_cfi::et, PVValHelper::eta, etMin_, iEvent, isoCut_, edm::HandleBase::isValid(), label_, eostools::move(), phi, VtxSmearedParameters_cfi::Phi, primaryVtxConstrain_, edm::Handle< T >::product(), pTMinTra_, relativeIsolation_, mps_fire::result, TtFullHadEvtBuilder_cfi::sumPt, HLT_FULL_cff::track, trackToken_, vertexToken_, z, and zMax_.

Member Data Documentation

◆ chi2Max_

float L1TkEmParticleProducer::chi2Max_

private

Definition at line 68 of file L1TkEmParticleProducer.cc.

Referenced by L1TkEmParticleProducer(), and produce().

◆ deltaZMax_

float L1TkEmParticleProducer::deltaZMax_

private

Definition at line 74 of file L1TkEmParticleProducer.cc.

Referenced by L1TkEmParticleProducer(), and produce().

◆ dRMax_

float L1TkEmParticleProducer::dRMax_

private

Definition at line 70 of file L1TkEmParticleProducer.cc.

Referenced by L1TkEmParticleProducer().

◆ dRMin_

float L1TkEmParticleProducer::dRMin_

private

Definition at line 69 of file L1TkEmParticleProducer.cc.

Referenced by L1TkEmParticleProducer(), and produce().

◆ egToken_

const edm::EDGetTokenT<EGammaBxCollection> L1TkEmParticleProducer::egToken_

private

Definition at line 79 of file L1TkEmParticleProducer.cc.

Referenced by produce().

◆ etMin_

float L1TkEmParticleProducer::etMin_

private

Definition at line 65 of file L1TkEmParticleProducer.cc.

Referenced by L1TkEmParticleProducer(), and produce().

◆ isoCut_

float L1TkEmParticleProducer::isoCut_

private

Definition at line 76 of file L1TkEmParticleProducer.cc.

Referenced by L1TkEmParticleProducer(), and produce().

◆ label_

std::string L1TkEmParticleProducer::label_

private

Definition at line 63 of file L1TkEmParticleProducer.cc.

Referenced by SequenceTypes.DummyModule::__repr__(), Mixins._Labelable::_findDependencies(), L1TkEmParticleProducer(), produce(), and Mixins._Labelable::setLabel().

◆ primaryVtxConstrain_

bool L1TkEmParticleProducer::primaryVtxConstrain_

private

Definition at line 72 of file L1TkEmParticleProducer.cc.

Referenced by L1TkEmParticleProducer(), and produce().

◆ pTMinTra_

float L1TkEmParticleProducer::pTMinTra_

private

Definition at line 71 of file L1TkEmParticleProducer.cc.

Referenced by L1TkEmParticleProducer(), and produce().

◆ relativeIsolation_

bool L1TkEmParticleProducer::relativeIsolation_

private

Definition at line 77 of file L1TkEmParticleProducer.cc.

Referenced by L1TkEmParticleProducer(), and produce().

◆ trackToken_

const edm::EDGetTokenT<std::vector<TTTrack<Ref_Phase2TrackerDigi_> > > L1TkEmParticleProducer::trackToken_

private

Definition at line 80 of file L1TkEmParticleProducer.cc.

Referenced by produce().

◆ vertexToken_

const edm::EDGetTokenT<TkPrimaryVertexCollection> L1TkEmParticleProducer::vertexToken_

private

Definition at line 81 of file L1TkEmParticleProducer.cc.

Referenced by produce().

◆ zMax_

float L1TkEmParticleProducer::zMax_

private

Definition at line 67 of file L1TkEmParticleProducer.cc.

Referenced by L1TkEmParticleProducer(), and produce().

Public Types

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ L1TTTrackCollectionType

◆ L1TTTrackType

Constructor & Destructor Documentation

◆ L1TkEmParticleProducer()

◆ ~L1TkEmParticleProducer()

Member Function Documentation

◆ CorrectedEta()

◆ fillDescriptions()

◆ produce()

Member Data Documentation

◆ chi2Max_

◆ deltaZMax_

◆ dRMax_

◆ dRMin_

◆ egToken_

◆ etMin_

◆ isoCut_

◆ label_

◆ primaryVtxConstrain_

◆ pTMinTra_

◆ relativeIsolation_

◆ trackToken_

◆ vertexToken_

◆ zMax_