L1TkEmParticleProducer Class Reference

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 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
template<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...
template<Transition B>
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)
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 ()
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()
template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()
template<Transition Tr = Transition::Event>
auto	produces () noexcept
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 ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
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

L1TTTrackCollectionType

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

Definition at line 46 of file L1TkEmParticleProducer.cc.

L1TTTrackType

typedef TTTrack<Ref_Phase2TrackerDigi_> L1TkEmParticleProducer::L1TTTrackType

Definition at line 45 of file L1TkEmParticleProducer.cc.

Constructor & Destructor Documentation

L1TkEmParticleProducer()

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::~L1TkEmParticleProducer ( )

override

Definition at line 100 of file L1TkEmParticleProducer.cc.

{}

Member Function Documentation

CorrectedEta()

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

Definition at line 237 of file L1TkEmParticleProducer.cc.

References funct::abs(), dumpMFGeometry_cfg::delta, PVValHelper::eta, EtaECal, JetChargeProducer_cfi::exp, dqm-mbProfile::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;
}

fillDescriptions()

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

static

Definition at line 277 of file L1TkEmParticleProducer.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and HLT_2022v11_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);
}

produce()

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(), BXVector< T >::begin(), hltPixelTracks_cff::chi2, chi2Max_, CorrectedEta(), reco::deltaR2(), deltaZMax_, dR2Min_, egToken_, BXVector< T >::end(), EgHLTOffHistBins_cfi::et, PVValHelper::eta, etMin_, iEvent, isoCut_, edm::HandleBase::isValid(), label_, eostools::move(), phi, VtxSmearedParameters_cfi::Phi, L1TkHTMissProducer_cfi::primaryVtxConstrain, primaryVtxConstrain_, edm::Handle< T >::product(), pTMinTra_, relativeIsolation_, mps_fire::result, TtFullHadEvtBuilder_cfi::sumPt, HLT_2022v11_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();
  }

  if (!L1TTTrackHandle.isValid()) {
    LogError("L1TkEmParticleProducer") << "\nWarning: L1TTTrackCollectionType not found in the event. Exit."
                                       << std::endl;
    return;
  }

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

chi2Max_

const float L1TkEmParticleProducer::chi2Max_

private

Definition at line 68 of file L1TkEmParticleProducer.cc.

Referenced by produce().

deltaZMax_

const float L1TkEmParticleProducer::deltaZMax_

private

Definition at line 73 of file L1TkEmParticleProducer.cc.

Referenced by produce().

dR2Max_

const float L1TkEmParticleProducer::dR2Max_

private

Definition at line 71 of file L1TkEmParticleProducer.cc.

dR2Min_

const float L1TkEmParticleProducer::dR2Min_

private

Definition at line 70 of file L1TkEmParticleProducer.cc.

Referenced by produce().

egToken_

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

private

Definition at line 60 of file L1TkEmParticleProducer.cc.

Referenced by produce().

etMin_

const float L1TkEmParticleProducer::etMin_

private

Definition at line 66 of file L1TkEmParticleProducer.cc.

Referenced by produce().

isoCut_

const float L1TkEmParticleProducer::isoCut_

private

Definition at line 75 of file L1TkEmParticleProducer.cc.

Referenced by produce().

label_

const std::string L1TkEmParticleProducer::label_

private

Definition at line 64 of file L1TkEmParticleProducer.cc.

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

primaryVtxConstrain_

const bool L1TkEmParticleProducer::primaryVtxConstrain_

private

Definition at line 72 of file L1TkEmParticleProducer.cc.

Referenced by produce().

pTMinTra_

const float L1TkEmParticleProducer::pTMinTra_

private

Definition at line 69 of file L1TkEmParticleProducer.cc.

Referenced by produce().

relativeIsolation_

const bool L1TkEmParticleProducer::relativeIsolation_

private

Definition at line 76 of file L1TkEmParticleProducer.cc.

Referenced by produce().

trackToken_

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

private

Definition at line 61 of file L1TkEmParticleProducer.cc.

Referenced by produce().

vertexToken_

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

private

Definition at line 62 of file L1TkEmParticleProducer.cc.

Referenced by produce().

zMax_

const float L1TkEmParticleProducer::zMax_

private

Definition at line 67 of file L1TkEmParticleProducer.cc.

Referenced by produce().