ZMuMuUserData Class Reference

Inheritance diagram for ZMuMuUserData:

Public Types
typedef math::XYZVector	Vector
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

Public Member Functions
	ZMuMuUserData (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
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
double	alpha_
EDGetTokenT< BeamSpot >	beamSpotToken_
double	beta_
int	counter
string	hltPath_
EDGetTokenT< VertexCollection >	primaryVerticesToken_
EDGetTokenT< std::vector< reco::CompositeCandidate > >	srcToken_
EDGetTokenT< GenParticleMatch >	zGenParticlesMatchToken_

Additional Inherited Members
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 29 of file ZMuMuUserData.cc.

Member Typedef Documentation

Vector

typedef math::XYZVector ZMuMuUserData::Vector

Definition at line 32 of file ZMuMuUserData.cc.

Constructor & Destructor Documentation

ZMuMuUserData()

ZMuMuUserData::ZMuMuUserData

(

const edm::ParameterSet &

cfg

)

Definition at line 46 of file ZMuMuUserData.cc.

    : srcToken_(consumes<std::vector<reco::CompositeCandidate> >(cfg.getParameter<InputTag>("src"))),
      beamSpotToken_(consumes<BeamSpot>(cfg.getParameter<InputTag>("beamSpot"))),
      primaryVerticesToken_(consumes<VertexCollection>(cfg.getParameter<InputTag>("primaryVertices"))),
      zGenParticlesMatchToken_(consumes<GenParticleMatch>(cfg.getParameter<InputTag>("zGenParticlesMatch"))),
      alpha_(cfg.getParameter<double>("alpha")),
      beta_(cfg.getParameter<double>("beta")),
      hltPath_(cfg.getParameter<std::string>("hltPath")) {
  produces<vector<pat::CompositeCandidate> >();
}

Member Function Documentation

produce()

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

overrideprivatevirtual

I fill the dau1 with positive and dau2 with negatove values for the pt, in order to flag the muons used for building zMassSa

I fill the dau1 with negatove and dau2 with positive values for the pt

Implements edm::EDProducer.

Definition at line 57 of file ZMuMuUserData.cc.

                                                          {
  Handle<std::vector<reco::CompositeCandidate> > dimuons;
  evt.getByToken(srcToken_, dimuons);
 
  Handle<BeamSpot> beamSpotHandle;
  evt.getByToken(beamSpotToken_, beamSpotHandle);
 
  Handle<VertexCollection> primaryVertices;  // Collection of primary Vertices
  evt.getByToken(primaryVerticesToken_, primaryVertices);
 
  bool isMCMatchTrue = false;
 
  Handle<GenParticleMatch> zGenParticlesMatch;
  if (evt.getByToken(zGenParticlesMatchToken_, zGenParticlesMatch)) {
    isMCMatchTrue = true;
  }
 
  //cout<<"isMCMatchTrue"<<isMCMatchTrue <<endl;
  unique_ptr<vector<pat::CompositeCandidate> > dimuonColl(new vector<pat::CompositeCandidate>());
 
  for (unsigned int i = 0; i < dimuons->size(); ++i) {
    const CompositeCandidate &z = (*dimuons)[i];
    //CandidateBaseRef zRef = dimuons ->refAt(i);
    edm::Ref<std::vector<reco::CompositeCandidate> > zRef(dimuons, i);
    pat::CompositeCandidate dimuon(z);
 
    float trueMass, truePt, trueEta, truePhi, trueY;
    if (isMCMatchTrue) {
      GenParticleRef trueZRef = (*zGenParticlesMatch)[zRef];
      //CandidateRef trueZRef = trueZIter->val;
      if (trueZRef.isNonnull()) {
        const Candidate &z = *trueZRef;
        trueMass = z.mass();
        truePt = z.pt();
        trueEta = z.eta();
        truePhi = z.phi();
        trueY = z.rapidity();
      } else {
        trueMass = -100;
        truePt = -100;
        trueEta = -100;
        truePhi = -100;
        trueY = -100;
      }
 
      dimuon.addUserFloat("TrueMass", trueMass);
      dimuon.addUserFloat("TruePt", truePt);
      dimuon.addUserFloat("TrueEta", trueEta);
      dimuon.addUserFloat("TruePhi", truePhi);
      dimuon.addUserFloat("TrueY", trueY);
    }
    const Candidate *dau1 = z.daughter(0);
    const Candidate *dau2 = z.daughter(1);
    const pat::Muon &mu1 = dynamic_cast<const pat::Muon &>(*dau1->masterClone());
    const pat::Muon &mu2 = dynamic_cast<const pat::Muon &>(*dau2->masterClone());
 
    /*cout<<"mu1 is null? "<<mu1.isMuon()<<endl;
    cout<<"mu2 is null? "<<mu2.isMuon()<<endl;
    cout<<"mu1 is global?"<<mu1.isGlobalMuon()<<endl;
    cout<<"mu2 is global?"<<mu2.isGlobalMuon()<<endl;
    */
 
    if (mu1.isGlobalMuon() == true && mu2.isGlobalMuon() == true) {
      TrackRef stAloneTrack1;
      TrackRef stAloneTrack2;
      Vector momentum;
      Candidate::PolarLorentzVector p4_1;
      double mu_mass;
      stAloneTrack1 = dau1->get<TrackRef, reco::StandAloneMuonTag>();
      stAloneTrack2 = dau2->get<TrackRef, reco::StandAloneMuonTag>();
      float zDau1SaEta = stAloneTrack1->eta();
      float zDau2SaEta = stAloneTrack2->eta();
      float zDau1SaPhi = stAloneTrack1->phi();
      float zDau2SaPhi = stAloneTrack2->phi();
      float zDau1SaPt, zDau2SaPt;
      if (counter % 2 == 0) {
        momentum = stAloneTrack1->momentum();
        p4_1 = dau2->polarP4();
        mu_mass = dau1->mass();
        zDau1SaPt = stAloneTrack1->pt();
        zDau2SaPt = -stAloneTrack2->pt();
      } else {
        momentum = stAloneTrack2->momentum();
        p4_1 = dau1->polarP4();
        mu_mass = dau2->mass();
        zDau1SaPt = -stAloneTrack1->pt();
        zDau2SaPt = stAloneTrack2->pt();
      }
 
      Candidate::PolarLorentzVector p4_2(momentum.rho(), momentum.eta(), momentum.phi(), mu_mass);
      double mass = (p4_1 + p4_2).mass();
      float zMassSa = mass;
      //cout<<"zMassSa "<<zMassSa;
      dimuon.addUserFloat("MassSa", zMassSa);
      dimuon.addUserFloat("Dau1SaPt", zDau1SaPt);
      dimuon.addUserFloat("Dau2SaPt", zDau2SaPt);
      dimuon.addUserFloat("Dau1SaPhi", zDau1SaPhi);
      dimuon.addUserFloat("Dau2SaPhi", zDau2SaPhi);
      dimuon.addUserFloat("Dau1SaEta", zDau1SaEta);
      dimuon.addUserFloat("Dau2SaEta", zDau2SaEta);
      ++counter;
    }
    dimuonColl->push_back(dimuon);
  }
 
  evt.put(std::move(dimuonColl));
}

References pat::PATObject< ObjectType >::addUserFloat(), beamSpotToken_, counter, WZMuSkim_cff::dimuons, reco::Candidate::get(), edm::Event::getByToken(), mps_fire::i, reco::Muon::isGlobalMuon(), edm::Ref< C, T, F >::isNonnull(), EgHLTOffHistBins_cfi::mass, reco::Candidate::mass(), reco::Candidate::masterClone(), eostools::move(), mu_mass, reco::Candidate::polarP4(), zMuMuMuonUserData::primaryVertices, primaryVerticesToken_, edm::Event::put(), srcToken_, z, zMuMu_SubskimPathsWithMCTruth_cff::zGenParticlesMatch, and zGenParticlesMatchToken_.

Member Data Documentation

alpha_

double ZMuMuUserData::alpha_

private

Definition at line 41 of file ZMuMuUserData.cc.

beamSpotToken_

EDGetTokenT<BeamSpot> ZMuMuUserData::beamSpotToken_

private

Definition at line 38 of file ZMuMuUserData.cc.

Referenced by produce().

beta_

double ZMuMuUserData::beta_

private

Definition at line 41 of file ZMuMuUserData.cc.

counter

int ZMuMuUserData::counter

private

Definition at line 43 of file ZMuMuUserData.cc.

Referenced by counter.Counters::__getitem__(), and produce().

hltPath_

string ZMuMuUserData::hltPath_

private

Definition at line 42 of file ZMuMuUserData.cc.

primaryVerticesToken_

EDGetTokenT<VertexCollection> ZMuMuUserData::primaryVerticesToken_

private

Definition at line 39 of file ZMuMuUserData.cc.

Referenced by produce().

srcToken_

EDGetTokenT<std::vector<reco::CompositeCandidate> > ZMuMuUserData::srcToken_

private

Definition at line 37 of file ZMuMuUserData.cc.

Referenced by produce().

zGenParticlesMatchToken_

EDGetTokenT<GenParticleMatch> ZMuMuUserData::zGenParticlesMatchToken_

private

Definition at line 40 of file ZMuMuUserData.cc.

Referenced by produce().