ZMuMuTrackUserData Class Reference

Inheritance diagram for ZMuMuTrackUserData:

Public Member Functions
	ZMuMuTrackUserData (const edm::ParameterSet &)
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
	~EDProducer () 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)

Private Member Functions
template<typename T >
vector< double >	isolation (const T &t, double ptThreshold, double etEcalThreshold, double etHcalThreshold, double dRVetoTrk, double dRTrk, double dREcal, double dRHcal, double alpha, double beta)
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
double	alpha_
EDGetTokenT< BeamSpot >	beamSpotToken_
double	beta_
double	dREcal_
double	dRHcal_
double	dRTrk_
double	dRVetoTrk_
double	etEcalThreshold_
double	etHcalThreshold_
EDGetTokenT< VertexCollection >	primaryVerticesToken_
double	ptThreshold_
EDGetTokenT< vector< pat::GenericParticle > >	srcToken_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	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
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
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 27 of file ZMuMuTrackUserData.cc.

Constructor & Destructor Documentation

ZMuMuTrackUserData()

ZMuMuTrackUserData::ZMuMuTrackUserData

(

const edm::ParameterSet &

cfg

)

Definition at line 100 of file ZMuMuTrackUserData.cc.

    : srcToken_(consumes<vector<pat::GenericParticle> >(cfg.getParameter<InputTag>("src"))),
      beamSpotToken_(consumes<BeamSpot>(cfg.getParameter<InputTag>("beamSpot"))),
      primaryVerticesToken_(consumes<VertexCollection>(cfg.getParameter<InputTag>("primaryVertices"))),
      ptThreshold_(cfg.getParameter<double>("ptThreshold")),
      etEcalThreshold_(cfg.getParameter<double>("etEcalThreshold")),
      etHcalThreshold_(cfg.getParameter<double>("etHcalThreshold")),
      dRVetoTrk_(cfg.getParameter<double>("dRVetoTrk")),
      dRTrk_(cfg.getParameter<double>("dRTrk")),
      dREcal_(cfg.getParameter<double>("dREcal")),
      dRHcal_(cfg.getParameter<double>("dRHcal")),
      alpha_(cfg.getParameter<double>("alpha")),
      beta_(cfg.getParameter<double>("beta")) {
  produces<std::vector<pat::GenericParticle> >();
}

Member Function Documentation

isolation()

template<typename T >

vector< double > ZMuMuTrackUserData::isolation ( const T &  t, double  ptThreshold, double  etEcalThreshold, double  etHcalThreshold, double  dRVetoTrk, double  dRTrk, double  dREcal, double  dRHcal, double  alpha, double  beta  )

private

Definition at line 53 of file ZMuMuTrackUserData.cc.

                                                          {
  vector<double> iso;
  const pat::IsoDeposit *trkIso = t.isoDeposit(pat::TrackIso);
  const pat::IsoDeposit *ecalIso = t.isoDeposit(pat::EcalIso);
  const pat::IsoDeposit *hcalIso = t.isoDeposit(pat::HcalIso);
 
  Direction dir = Direction(t.eta(), t.phi());
 
  pat::IsoDeposit::AbsVetos vetosTrk;
  vetosTrk.push_back(new ConeVeto(dir, dRVetoTrk));
  vetosTrk.push_back(new ThresholdVeto(ptThreshold));
 
  pat::IsoDeposit::AbsVetos vetosEcal;
  vetosEcal.push_back(new ConeVeto(dir, 0.));
  vetosEcal.push_back(new ThresholdVeto(etEcalThreshold));
 
  pat::IsoDeposit::AbsVetos vetosHcal;
  vetosHcal.push_back(new ConeVeto(dir, 0.));
  vetosHcal.push_back(new ThresholdVeto(etHcalThreshold));
 
  double isovalueTrk = (trkIso->sumWithin(dRTrk, vetosTrk));
  double isovalueEcal = (ecalIso->sumWithin(dREcal, vetosEcal));
  double isovalueHcal = (hcalIso->sumWithin(dRHcal, vetosHcal));
 
  iso.push_back(isovalueTrk);
  iso.push_back(isovalueEcal);
  iso.push_back(isovalueHcal);
 
  //double iso =  isovalueTrk + isovalueEcal + isovalueHcal;
  double combIso =
      alpha * (((1 + beta) / 2 * isovalueEcal) + ((1 - beta) / 2 * isovalueHcal)) + ((1 - alpha) * isovalueTrk);
 
  iso.push_back(combIso);
  double relIso = combIso /= t.pt();
  iso.push_back(relIso);
  return iso;
}

References zMuMuMuonUserData::alpha, zMuMuMuonUserData::beta, DeadROC_duringRun::dir, zMuMuMuonUserData::dREcal, zMuMuMuonUserData::dRHcal, zMuMuMuonUserData::dRTrk, zMuMuMuonUserData::dRVetoTrk, pat::EcalIso, zMuMuMuonUserData::etEcalThreshold, zMuMuMuonUserData::etHcalThreshold, pat::HcalIso, ewkMuLumiMonitorDQM_cfi::ptThreshold, objects.autophobj::relIso, reco::IsoDeposit::sumWithin(), submitPVValidationJobs::t, pat::TrackIso, and ZMuMuCategoriesSequences_cff::trkIso.

Referenced by produce().

produce()

void ZMuMuTrackUserData::produce ( edm::Event &  evt, const edm::EventSetup &    )

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 116 of file ZMuMuTrackUserData.cc.

                                                               {
  Handle<vector<pat::GenericParticle> > tracks;
  evt.getByToken(srcToken_, tracks);
 
  Handle<BeamSpot> beamSpotHandle;
  evt.getByToken(beamSpotToken_, beamSpotHandle);
 
  Handle<VertexCollection> primaryVertices;  // Collection of primary Vertices
  evt.getByToken(primaryVerticesToken_, primaryVertices);
 
  unique_ptr<vector<pat::GenericParticle> > tkColl(new vector<pat::GenericParticle>(*tracks));
  for (unsigned int i = 0; i < tkColl->size(); ++i) {
    pat::GenericParticle &tk = (*tkColl)[i];
    vector<double> iso = isolation(
        tk, ptThreshold_, etEcalThreshold_, etHcalThreshold_, dRVetoTrk_, dRTrk_, dREcal_, dRHcal_, alpha_, beta_);
    tk.setIsolation(pat::User1Iso, iso[0]);
    //    cout << "track User1Iso " << iso[0] << endl;
    tk.setIsolation(pat::User2Iso, iso[1]);
    //cout << "track User2Iso " << iso[1] << endl;
    tk.setIsolation(pat::User3Iso, iso[2]);
    //cout << "track User3Iso " << iso[2] << endl;
    tk.setIsolation(pat::User4Iso, iso[3]);
    //cout << "track User4Iso " << iso[3] << endl;
    tk.setIsolation(pat::User5Iso, iso[4]);
    //cout << "track User5Iso " << iso[4] << endl;
 
    float zDaudxyFromBS = -1;
    float zDaudzFromBS = -1;
    float zDaudxyFromPV = -1;
    float zDaudzFromPV = -1;
    float zDauNofMuChambers = -1;
    float zDauNofMuMatches = -1;
    float zDauChi2 = -1;
    float zDauTrkChi2 = -1;
    float zDauSaChi2 = -1;
    float zDauNofMuonHits = -1;
    float zDauNofStripHits = -1;
    float zDauNofPixelHits = -1;
    float zDauMuEnergyEm = -1;
    float zDauMuEnergyHad = -1;
 
    TrackRef muTrkRef = tk.track();
    if (muTrkRef.isNonnull()) {
      zDaudxyFromBS = muTrkRef->dxy(beamSpotHandle->position());
      zDaudzFromBS = muTrkRef->dz(beamSpotHandle->position());
      zDaudxyFromPV = muTrkRef->dxy(primaryVertices->begin()->position());
      zDaudzFromPV = muTrkRef->dz(primaryVertices->begin()->position());
      zDauChi2 = muTrkRef->normalizedChi2();
      zDauTrkChi2 = muTrkRef->normalizedChi2();
      zDauNofStripHits = muTrkRef->hitPattern().numberOfValidStripHits();
      zDauNofPixelHits = muTrkRef->hitPattern().numberOfValidPixelHits();
    }
    tk.addUserFloat("zDau_dxyFromBS", zDaudxyFromBS);
    tk.addUserFloat("zDau_dzFromBS", zDaudzFromBS);
    tk.addUserFloat("zDau_dxyFromPV", zDaudxyFromPV);
    tk.addUserFloat("zDau_dzFromPV", zDaudzFromPV);
    tk.addUserFloat("zDau_NofMuonHits", zDauNofMuonHits);
    tk.addUserFloat("zDau_TrkNofStripHits", zDauNofStripHits);
    tk.addUserFloat("zDau_TrkNofPixelHits", zDauNofPixelHits);
    tk.addUserFloat("zDau_NofMuChambers", zDauNofMuChambers);
    tk.addUserFloat("zDau_NofMuMatches", zDauNofMuMatches);
    tk.addUserFloat("zDau_Chi2", zDauChi2);
    tk.addUserFloat("zDau_TrkChi2", zDauTrkChi2);
    tk.addUserFloat("zDau_SaChi2", zDauSaChi2);
    tk.addUserFloat("zDau_MuEnergyEm", zDauMuEnergyEm);
    tk.addUserFloat("zDau_MuEnergyHad", zDauMuEnergyHad);
  }
 
  evt.put(std::move(tkColl));
}

References pat::PATObject< ObjectType >::addUserFloat(), alpha_, beamSpotToken_, beta_, dREcal_, dRHcal_, dRTrk_, dRVetoTrk_, etEcalThreshold_, etHcalThreshold_, edm::Event::getByToken(), mps_fire::i, edm::Ref< C, T, F >::isNonnull(), isolation(), eostools::move(), reco::BeamSpot::position(), zMuMuMuonUserData::primaryVertices, primaryVerticesToken_, ptThreshold_, edm::Event::put(), pat::GenericParticle::setIsolation(), srcToken_, pat::GenericParticle::track(), PDWG_EXOHSCP_cff::tracks, pat::User1Iso, pat::User2Iso, pat::User3Iso, pat::User4Iso, and pat::User5Iso.

Member Data Documentation

alpha_

double ZMuMuTrackUserData::alpha_

private

Definition at line 38 of file ZMuMuTrackUserData.cc.

Referenced by produce().

beamSpotToken_

EDGetTokenT<BeamSpot> ZMuMuTrackUserData::beamSpotToken_

private

Definition at line 35 of file ZMuMuTrackUserData.cc.

Referenced by produce().

beta_

double ZMuMuTrackUserData::beta_

private

Definition at line 38 of file ZMuMuTrackUserData.cc.

Referenced by produce().

dREcal_

double ZMuMuTrackUserData::dREcal_

private

Definition at line 37 of file ZMuMuTrackUserData.cc.

Referenced by produce().

dRHcal_

double ZMuMuTrackUserData::dRHcal_

private

Definition at line 37 of file ZMuMuTrackUserData.cc.

Referenced by produce().

dRTrk_

double ZMuMuTrackUserData::dRTrk_

private

Definition at line 37 of file ZMuMuTrackUserData.cc.

Referenced by produce().

dRVetoTrk_

double ZMuMuTrackUserData::dRVetoTrk_

private

Definition at line 37 of file ZMuMuTrackUserData.cc.

Referenced by produce().

etEcalThreshold_

double ZMuMuTrackUserData::etEcalThreshold_

private

Definition at line 37 of file ZMuMuTrackUserData.cc.

Referenced by produce().

etHcalThreshold_

double ZMuMuTrackUserData::etHcalThreshold_

private

Definition at line 37 of file ZMuMuTrackUserData.cc.

Referenced by produce().

primaryVerticesToken_

EDGetTokenT<VertexCollection> ZMuMuTrackUserData::primaryVerticesToken_

private

Definition at line 36 of file ZMuMuTrackUserData.cc.

Referenced by produce().

ptThreshold_

double ZMuMuTrackUserData::ptThreshold_

private

Definition at line 37 of file ZMuMuTrackUserData.cc.

Referenced by produce().

srcToken_

EDGetTokenT<vector<pat::GenericParticle> > ZMuMuTrackUserData::srcToken_

private

Definition at line 34 of file ZMuMuTrackUserData.cc.

Referenced by produce().