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

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)

	~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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

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
const float	chi2Max_

const float	deltaZMax_

const float	dR2Max_

const float	dR2Min_

const edm::EDGetTokenT < EGammaBxCollection >	egToken_

const float	etMin_

const float	isoCut_

const std::string	label_

const bool	primaryVtxConstrain_

const float	pTMinTra_

const bool	relativeIsolation_

const edm::EDGetTokenT < L1TTTrackCollectionType >	trackToken_

const edm::EDGetTokenT < TkPrimaryVertexCollection >	vertexToken_

const float	zMax_

Additional Inherited Members
Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

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

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

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

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

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 43 of file L1TkEmParticleProducer.cc.

Member Typedef Documentation

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

Definition at line 46 of file L1TkEmParticleProducer.cc.

typedef TTTrack<Ref_Phase2TrackerDigi_> L1TkEmParticleProducer::L1TTTrackType

Definition at line 45 of file L1TkEmParticleProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 82 of file L1TkEmParticleProducer.cc.

References label_.

     : egToken_(consumes<EGammaBxCollection>(iConfig.getParameter<edm::InputTag>("L1EGammaInputTag"))),
       trackToken_(consumes<L1TTTrackCollectionType>(iConfig.getParameter<edm::InputTag>("L1TrackInputTag"))),
       vertexToken_(consumes<TkPrimaryVertexCollection>(iConfig.getParameter<edm::InputTag>("L1VertexInputTag"))),
       label_(iConfig.getParameter<std::string>("label")),
       etMin_((float)iConfig.getParameter<double>("ETmin")),
       zMax_((float)iConfig.getParameter<double>("ZMAX")),
       chi2Max_((float)iConfig.getParameter<double>("CHI2MAX")),
       pTMinTra_((float)iConfig.getParameter<double>("PTMINTRA")),
       dR2Min_((float)iConfig.getParameter<double>("DRmin") * (float)iConfig.getParameter<double>("DRmin")),
       dR2Max_((float)iConfig.getParameter<double>("DRmax") * (float)iConfig.getParameter<double>("DRmax")),
       primaryVtxConstrain_(iConfig.getParameter<bool>("PrimaryVtxConstrain")),
       deltaZMax_((float)iConfig.getParameter<double>("DeltaZMax")),
       isoCut_((float)iConfig.getParameter<double>("IsoCut")),
       relativeIsolation_(iConfig.getParameter<bool>("RelativeIsolation")) {
   produces<TkEmCollection>(label_);
 }

L1TkEmParticleProducer::~L1TkEmParticleProducer ( )

override

Definition at line 100 of file L1TkEmParticleProducer.cc.

100 {}

Member Function Documentation

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

Definition at line 231 of file L1TkEmParticleProducer.cc.

References funct::abs(), CommonMethods::delta(), PVValHelper::eta, EtaECal, funct::exp(), log, M_PI, REcal, funct::tan(), ZEcal, and gpuVertexFinder::zv.

Referenced by produce().

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

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

static

Definition at line 271 of file L1TkEmParticleProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, and HLT_FULL_cff::InputTag.

                                                                                         {
   //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.add<string>("label", "EG");
   desc.add<edm::InputTag>("L1EGammaInputTag", edm::InputTag("simCaloStage2Digis"));
   desc.add<edm::InputTag>("L1TrackInputTag", edm::InputTag("TTTracksFromTrackletEmulation", "Level1TTTracks"));
   desc.add<edm::InputTag>("L1VertexInputTag", edm::InputTag("L1TkPrimaryVertex"));
   desc.add<double>("ETmin", -1.);
   desc.add<bool>("RelativeIsolation", true);
   desc.add<double>("IsoCut", 0.23);
   desc.add<double>("ZMAX", 25.);
   desc.add<double>("CHI2MAX", 100.);
   desc.add<double>("PTMINTRA", 2.);
   desc.add<double>("DRmin", 0.07);
   desc.add<double>("DRmax", 0.30);
   desc.add<bool>("PrimaryVtxConstrain", false);
   desc.add<double>("DeltaZMax", 0.6);
   descriptions.add("l1TkEmParticleProducer", desc);
 }

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

overrideprivatevirtual

Implements edm::global::EDProducerBase.

Definition at line 103 of file L1TkEmParticleProducer.cc.

References funct::abs(), HLT_FULL_cff::chi2, chi2Max_, CorrectedEta(), reco::deltaR2(), deltaZMax_, dR2Min_, egToken_, PVValHelper::eta, etMin_, edm::Event::getByToken(), isoCut_, label_, eostools::move(), phi, colinearityKinematic::Phi, primaryVtxConstrain_, pTMinTra_, edm::Event::put(), relativeIsolation_, mps_fire::result, HLT_FULL_cff::track, trackToken_, vertexToken_, z, and zMax_.

                                                                                                      {
   using namespace edm;
 
   auto result = std::make_unique<TkEmCollection>();
 
   // the L1Tracks
   edm::Handle<L1TTTrackCollectionType> L1TTTrackHandle;
   iEvent.getByToken(trackToken_, L1TTTrackHandle);
   if (!L1TTTrackHandle.isValid()) {
     LogError("L1TkEmParticleProducer") << "\nWarning: L1TTTrackCollectionType not found in the event. Exit."
                                        << std::endl;
     return;
   }
 
   // the primary vertex (used only if primaryVtxConstrain_ = true)
   float zvtxL1tk = -999;
   bool primaryVtxConstrain = primaryVtxConstrain_;
   edm::Handle<TkPrimaryVertexCollection> L1VertexHandle;
   iEvent.getByToken(vertexToken_, L1VertexHandle);
   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();
   }
 
   // Now loop over the L1EGamma objects
 
   edm::Handle<EGammaBxCollection> eGammaHandle;
   iEvent.getByToken(egToken_, eGammaHandle);
   if (!eGammaHandle.isValid()) {
     LogError("L1TkEmParticleProducer") << "\nWarning: L1EmCollection not found in the event. Exit." << std::endl;
     return;
   }
   EGammaBxCollection eGammaCollection = (*eGammaHandle.product());
   EGammaBxCollection::const_iterator egIter;
 
   int ieg = 0;
   for (egIter = eGammaCollection.begin(0); egIter != eGammaCollection.end(0); ++egIter)  // considering BX = only
   {
     edm::Ref<EGammaBxCollection> EGammaRef(eGammaHandle, ieg);
     ieg++;
 
     float et = egIter->et();
     if (et < etMin_)
       continue;
 
     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(eta, zvtxL1tk);
 
     float phi = egIter->phi();
 
     // calculate the isolation of the L1EG object with
     // respect to L1Tracks.
 
     float sumPt = 0;
     float sumPtPV = 0;
 
     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 dr2 = reco::deltaR2(Eta, Phi, eta, phi);
       if (dr2 < dR2Max_ && dr2 >= dR2Min_) {
         sumPt += Pt;
       }
 
       if (zvtxL1tk > -999 && std::abs(z - zvtxL1tk) >= deltaZMax_)
         continue;  // Now, PV constrained trackSum:
 
       dr2 = reco::deltaR2(Eta, Phi, etaPV, phi);  // recompute using the corrected eta
       if (dr2 < dR2Max_ && dr2 >= dR2Min_) {
         sumPtPV += Pt;
       }
     }  // end loop over tracks
 
     float trkisol = -999;
     float trkisolPV = -999;
     if (relativeIsolation_) {
       if (et > 0) {
         trkisol = sumPt / et;
         trkisolPV = sumPtPV / et;
       }       // relative isolation
     } else {  // absolute isolation
       trkisol = sumPt;
       trkisolPV = sumPtPV;
     }
 
     float isolation = primaryVtxConstrain ? trkisolPV : trkisol;
 
     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_);
 }

Member Data Documentation

const float L1TkEmParticleProducer::chi2Max_

private

Definition at line 68 of file L1TkEmParticleProducer.cc.

Referenced by produce().

const float L1TkEmParticleProducer::deltaZMax_

private

Definition at line 73 of file L1TkEmParticleProducer.cc.

Referenced by produce().

const float L1TkEmParticleProducer::dR2Max_

private

Definition at line 71 of file L1TkEmParticleProducer.cc.

const float L1TkEmParticleProducer::dR2Min_

private

Definition at line 70 of file L1TkEmParticleProducer.cc.

Referenced by produce().

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

private

Definition at line 60 of file L1TkEmParticleProducer.cc.

Referenced by produce().

const float L1TkEmParticleProducer::etMin_

private

Definition at line 66 of file L1TkEmParticleProducer.cc.

Referenced by produce().

const float L1TkEmParticleProducer::isoCut_

private

Definition at line 75 of file L1TkEmParticleProducer.cc.

Referenced by produce().

const std::string L1TkEmParticleProducer::label_

private

Definition at line 64 of file L1TkEmParticleProducer.cc.

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

const bool L1TkEmParticleProducer::primaryVtxConstrain_

private

Definition at line 72 of file L1TkEmParticleProducer.cc.

Referenced by produce().

const float L1TkEmParticleProducer::pTMinTra_

private

Definition at line 69 of file L1TkEmParticleProducer.cc.

Referenced by produce().

const bool L1TkEmParticleProducer::relativeIsolation_

private

Definition at line 76 of file L1TkEmParticleProducer.cc.

Referenced by produce().

const edm::EDGetTokenT<L1TTTrackCollectionType> L1TkEmParticleProducer::trackToken_

private

Definition at line 61 of file L1TkEmParticleProducer.cc.

Referenced by produce().

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

private

Definition at line 62 of file L1TkEmParticleProducer.cc.

Referenced by produce().

const float L1TkEmParticleProducer::zMax_

private

Definition at line 67 of file L1TkEmParticleProducer.cc.

Referenced by produce().

Public Types

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation