Inheritance diagram for ZMuMuMuonUserData:

Public Member Functions
	ZMuMuMuonUserData (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

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

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

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_

string	hltPath_

EDGetTokenT< VertexCollection >	primaryVerticesToken_

double	ptThreshold_

EDGetTokenT< vector< pat::Muon > >	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	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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)

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<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 28 of file ZMuMuMuonUserData.cc.

Constructor & Destructor Documentation

ZMuMuMuonUserData::ZMuMuMuonUserData ( const edm::ParameterSet & cfg )

Definition at line 110 of file ZMuMuMuonUserData.cc.

     : srcToken_(consumes<vector<pat::Muon> >(cfg.getParameter<InputTag>("src"))),
       beamSpotToken_(consumes<BeamSpot>(cfg.getParameter<InputTag>("beamSpot"))),
       primaryVerticesToken_(consumes<VertexCollection>(cfg.getParameter<InputTag>("primaryVertices"))),
       alpha_(cfg.getParameter<double>("alpha")),
       beta_(cfg.getParameter<double>("beta")),
       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")),
       hltPath_(cfg.getParameter<std::string>("hltPath")) {
   produces<std::vector<pat::Muon> >();
 }

Member Function Documentation

template<typename T >

vector< double > ZMuMuMuonUserData::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 55 of file ZMuMuMuonUserData.cc.

References zMuMuMuonUserData::beta, DeadROC_duringRun::dir, pat::EcalIso, pat::HcalIso, objects.autophobj::relIso, reco::IsoDeposit::sumWithin(), pat::TrackIso, and ZMuMuCategoriesSequences_cff::trkIso.

Referenced by produce().

                                                          {
   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);
   //cout<<"isoTrk"<<isovalueTrk<<" "<<t.trackIso()<<endl;
   iso.push_back(isovalueEcal);
   //cout<<"isoEcal"<<isovalueEcal<<" "<<t.ecalIso()<<endl;
   iso.push_back(isovalueHcal);
   //cout<<"isoHcal"<<isovalueHcal<<" "<<t.hcalIso()<<endl;
   //double isovalueTrk  = t.trackIso();
   //double isovalueEcal = t.ecalIso();
   //double isovalueHcal = t.hcalIso();
 
   //double iso =  isovalueTrk + isovalueEcal + isovalueHcal;
   double combIso =
       alpha * (((1 + beta) / 2 * isovalueEcal) + ((1 - beta) / 2 * isovalueHcal)) + ((1 - alpha) * isovalueTrk);
   iso.push_back(combIso);
   //cout<<"combIso"<<iso[3]<<endl;
 
   double relIso = combIso /= t.pt();
   iso.push_back(relIso);
   //cout<<"relIso"<<iso[4]<<endl;
   return iso;
 }

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

overrideprivate

Definition at line 127 of file ZMuMuMuonUserData.cc.

References pat::PATObject< ObjectType >::addUserFloat(), alpha_, beamSpotToken_, beta_, reco::Muon::calEnergy(), DEFINE_FWK_MODULE, dREcal_, dRHcal_, dRTrk_, dRVetoTrk_, reco::MuonEnergy::em, etEcalThreshold_, etHcalThreshold_, edm::Event::getByToken(), pat::Muon::globalTrack(), reco::MuonEnergy::had, hltPath_, mps_fire::i, pat::Muon::innerTrack(), reco::Muon::isGlobalMuon(), edm::Ref< C, T, F >::isNonnull(), isolation(), reco::Muon::isStandAloneMuon(), reco::Muon::isTrackerMuon(), visualization-live-secondInstance_cfg::m, eostools::move(), PDWG_BPHSkim_cff::muons, reco::Muon::numberOfChambers(), reco::Muon::numberOfMatches(), pat::Muon::outerTrack(), reco::BeamSpot::position(), zMuMuMuonUserData::primaryVertices, primaryVerticesToken_, ptThreshold_, edm::Event::put(), pat::Lepton< LeptonType >::setIsolation(), srcToken_, pat::PATObject< ObjectType >::triggerObjectMatchesByPath(), pat::User1Iso, pat::User2Iso, pat::User3Iso, pat::User4Iso, and pat::User5Iso.

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

                                                               {
   Handle<vector<pat::Muon> > muons;
   evt.getByToken(srcToken_, muons);
 
   Handle<BeamSpot> beamSpotHandle;
   evt.getByToken(beamSpotToken_, beamSpotHandle);
 
   Handle<VertexCollection> primaryVertices;  // Collection of primary Vertices
   evt.getByToken(primaryVerticesToken_, primaryVertices);
 
   unique_ptr<vector<pat::Muon> > muonColl(new vector<pat::Muon>(*muons));
   for (unsigned int i = 0; i < muonColl->size(); ++i) {
     pat::Muon &m = (*muonColl)[i];
     //pat::Muon *mu = new pat::Muon(m);
     vector<double> iso = isolation(
         m, ptThreshold_, etEcalThreshold_, etHcalThreshold_, dRVetoTrk_, dRTrk_, dREcal_, dRHcal_, alpha_, beta_);
     m.setIsolation(pat::User1Iso, iso[0]);
     //cout<<"muon User1Iso "<<iso[0]<<endl;
     m.setIsolation(pat::User2Iso, iso[1]);
     //cout<<"iso2 "<<iso[1]<<endl;
     m.setIsolation(pat::User3Iso, iso[2]);
     //cout<<"iso3 "<<iso[2]<<endl;
     m.setIsolation(pat::User4Iso, iso[3]);
     //cout<<"iso4 "<<iso[3]<<endl;
     m.setIsolation(pat::User5Iso, iso[4]);
     //cout<<"iso5 "<<iso[4]<<endl;
     float zDauMuEnergyEm = m.calEnergy().em;
     float zDauMuEnergyHad = m.calEnergy().had;
     TrackRef muGlbRef = m.globalTrack();
     TrackRef muTrkRef = m.innerTrack();
     TrackRef muSaRef = m.outerTrack();
     float zDaudxyFromBS = -1;
     float zDaudzFromBS = -1;
     float zDaudxyFromPV = -1;
     float zDaudzFromPV = -1;
     int zDauNofMuChambers = m.numberOfChambers();
     int zDauNofMuMatches = m.numberOfMatches();
     //  for the following variables looking at global/trk and sta at the same time
     float zDauChi2 = -1;
     float zDauTrkChi2 = -1;
     float zDauSaChi2 = -1;
     float zDauNofMuonHits = -1;
     float zDauSaNofMuonHits = -1;
     float zDauNofStripHits = -1;
     float zDauTrkNofStripHits = -1;
     float zDauNofPixelHits = -1;
     float zDauTrkNofPixelHits = -1;
 
     if (muGlbRef.isNonnull() && m.isGlobalMuon() == true) {
       zDaudxyFromBS = muGlbRef->dxy(beamSpotHandle->position());
       zDaudzFromBS = muGlbRef->dz(beamSpotHandle->position());
       zDaudxyFromPV = muGlbRef->dxy(primaryVertices->begin()->position());
       zDaudzFromPV = muGlbRef->dz(primaryVertices->begin()->position());
       zDauChi2 = muGlbRef->normalizedChi2();
       zDauTrkChi2 = muTrkRef->normalizedChi2();
       zDauSaChi2 = muSaRef->normalizedChi2();
       zDauNofMuonHits = muGlbRef->hitPattern().numberOfValidMuonHits();
       zDauSaNofMuonHits = muSaRef->hitPattern().numberOfValidMuonHits();
       zDauNofStripHits = muGlbRef->hitPattern().numberOfValidStripHits();
       zDauTrkNofStripHits = muTrkRef->hitPattern().numberOfValidStripHits();
       zDauNofPixelHits = muGlbRef->hitPattern().numberOfValidPixelHits();
       zDauTrkNofPixelHits = muTrkRef->hitPattern().numberOfValidPixelHits();
     } else if (muSaRef.isNonnull() && m.isStandAloneMuon() == true) {
       zDaudxyFromBS = muSaRef->dxy(beamSpotHandle->position());
       zDaudzFromBS = muSaRef->dz(beamSpotHandle->position());
       zDaudxyFromPV = muSaRef->dxy(primaryVertices->begin()->position());
       zDaudzFromPV = muSaRef->dz(primaryVertices->begin()->position());
       zDauSaChi2 = muSaRef->normalizedChi2();
       zDauSaNofMuonHits = muSaRef->hitPattern().numberOfValidMuonHits();
 
     } else if (muTrkRef.isNonnull() && m.isTrackerMuon() == true) {
       zDaudxyFromBS = muTrkRef->dxy(beamSpotHandle->position());
       zDaudzFromBS = muTrkRef->dz(beamSpotHandle->position());
       zDaudxyFromPV = muTrkRef->dxy(primaryVertices->begin()->position());
       zDaudzFromPV = muTrkRef->dz(primaryVertices->begin()->position());
       zDauTrkChi2 = muTrkRef->normalizedChi2();
       zDauTrkNofStripHits = muTrkRef->hitPattern().numberOfValidStripHits();
       zDauTrkNofPixelHits = muTrkRef->hitPattern().numberOfValidPixelHits();
     }
 
     const pat::TriggerObjectStandAloneCollection muHLTMatches = m.triggerObjectMatchesByPath(hltPath_);
     float muHLTBit;
     int dimTrig = muHLTMatches.size();
     if (dimTrig != 0) {
       muHLTBit = 1;
     } else {
       muHLTBit = 0;
     }
     m.addUserFloat("zDau_dxyFromBS", zDaudxyFromBS);
     m.addUserFloat("zDau_dzFromBS", zDaudzFromBS);
     m.addUserFloat("zDau_dxyFromPV", zDaudxyFromPV);
     m.addUserFloat("zDau_dzFromPV", zDaudzFromPV);
     m.addUserFloat("zDau_HLTBit", muHLTBit);
     m.addUserFloat("zDau_dzFromPV", zDaudzFromPV);
     m.addUserFloat("zDau_Chi2", zDauChi2);
     m.addUserFloat("zDau_TrkChi2", zDauTrkChi2);
     m.addUserFloat("zDau_SaChi2", zDauSaChi2);
     m.addUserFloat("zDau_NofMuonHits", zDauNofMuonHits);
     m.addUserFloat("zDau_SaNofMuonHits", zDauSaNofMuonHits);
     m.addUserFloat("zDau_NofStripHits", zDauNofStripHits);
     m.addUserFloat("zDau_TrkNofStripHits", zDauTrkNofStripHits);
     m.addUserFloat("zDau_NofPixelHits", zDauNofPixelHits);
     m.addUserFloat("zDau_TrkNofPixelHits", zDauTrkNofPixelHits);
     m.addUserFloat("zDau_NofMuChambers", zDauNofMuChambers);
     m.addUserFloat("zDau_NofMuMatches", zDauNofMuMatches);
     m.addUserFloat("zDau_MuEnergyEm", zDauMuEnergyEm);
     m.addUserFloat("zDau_MuEnergyHad", zDauMuEnergyHad);
   }
 
   evt.put(std::move(muonColl));
 }

Member Data Documentation

double ZMuMuMuonUserData::alpha_

private

Definition at line 38 of file ZMuMuMuonUserData.cc.

Referenced by produce().

EDGetTokenT<BeamSpot> ZMuMuMuonUserData::beamSpotToken_

private

Definition at line 36 of file ZMuMuMuonUserData.cc.

Referenced by produce().

double ZMuMuMuonUserData::beta_

private

Definition at line 38 of file ZMuMuMuonUserData.cc.

Referenced by produce().

double ZMuMuMuonUserData::dREcal_

private

Definition at line 39 of file ZMuMuMuonUserData.cc.

Referenced by produce().

double ZMuMuMuonUserData::dRHcal_

private

Definition at line 39 of file ZMuMuMuonUserData.cc.

Referenced by produce().

double ZMuMuMuonUserData::dRTrk_

private

Definition at line 39 of file ZMuMuMuonUserData.cc.

Referenced by produce().

double ZMuMuMuonUserData::dRVetoTrk_

private

Definition at line 39 of file ZMuMuMuonUserData.cc.

Referenced by produce().

double ZMuMuMuonUserData::etEcalThreshold_

private

Definition at line 39 of file ZMuMuMuonUserData.cc.

Referenced by produce().

double ZMuMuMuonUserData::etHcalThreshold_

private

Definition at line 39 of file ZMuMuMuonUserData.cc.

Referenced by produce().

string ZMuMuMuonUserData::hltPath_

private

Definition at line 40 of file ZMuMuMuonUserData.cc.

Referenced by produce().

EDGetTokenT<VertexCollection> ZMuMuMuonUserData::primaryVerticesToken_

private

Definition at line 37 of file ZMuMuMuonUserData.cc.

Referenced by produce().

double ZMuMuMuonUserData::ptThreshold_

private

Definition at line 39 of file ZMuMuMuonUserData.cc.

Referenced by produce().

EDGetTokenT<vector<pat::Muon> > ZMuMuMuonUserData::srcToken_

private

Definition at line 35 of file ZMuMuMuonUserData.cc.

Referenced by produce().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation