Inheritance diagram for ZMassHistogrammer:

Public Member Functions
	ZMassHistogrammer (const edm::ParameterSet &pset)

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () 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

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)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &setup) override

Private Attributes
edm::EDGetTokenT< reco::CandidateView >	genToken_

TH1F *	h_mZ_

TH1F *	h_mZMC_

edm::EDGetTokenT< reco::CandidateView >	zToken_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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

Constructor & Destructor Documentation

ZMassHistogrammer::ZMassHistogrammer ( const edm::ParameterSet & pset )

Definition at line 30 of file ZMassHistogrammer.cc.

References gather_cfg::cout, h_mZ_, h_mZMC_, and TFileService::make().

                                                              :
   zToken_(consumes<reco::CandidateView>(pset.getParameter<InputTag>("z"))),
   genToken_(consumes<reco::CandidateView>(pset.getParameter<InputTag>("gen"))) {
   cout << ">>> Z Mass constructor" << endl;
   Service<TFileService> fs;
   h_mZ_ = fs->make<TH1F>("ZMass", "Z mass (GeV/c^{2})", 100,  0, 200);
   h_mZMC_ = fs->make<TH1F>("ZMCMass", "Z MC mass (GeV/c^{2})", 100,  0, 200);
 }

Member Function Documentation

void ZMassHistogrammer::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

overrideprivate

Definition at line 39 of file ZMassHistogrammer.cc.

References reco::Candidate::charge(), gather_cfg::cout, reco::Candidate::daughter(), DEFINE_FWK_MODULE, relval_steps::gen(), genToken_, h_mZ_, h_mZMC_, mps_fire::i, reco::Candidate::mass(), reco::Candidate::numberOfDaughters(), reco::Candidate::pdgId(), edm::OwnVector< T, P >::size(), reco::Candidate::status(), z, and zToken_.

                                                                                  {
   cout << ">>> Z Mass analyze" << endl;
   Handle<CandidateCollection> z;
   Handle<CandidateCollection> gen;
   event.getByToken(zToken_, z);
   event.getByToken(genToken_, gen);
   for(unsigned int i = 0; i < z->size(); ++i) {
     const Candidate &zCand = (*z)[i];
     h_mZ_->Fill(zCand.mass());
   }
   for(unsigned int i = 0; i < gen->size(); ++i) {
     const Candidate &genCand = (*gen)[i];
     if((genCand.pdgId() == 23) && (genCand.status() == 2)) //this is an intermediate Z0
       cout << ">>> intermediate Z0 found, with " << genCand.numberOfDaughters()
        << " daughters" << endl;
     if((genCand.pdgId() == 23)&&(genCand.status() == 3)) { //this is a Z0
       cout << ">>> Z0 found, with " << genCand.numberOfDaughters()
        << " daughters" << endl;
       h_mZMC_->Fill(genCand.mass());
       if(genCand.numberOfDaughters() == 3) {//Z0 decays in mu+ mu-, the 3rd daughter is the same Z0
     const Candidate * dauGen0 = genCand.daughter(0);
     const Candidate * dauGen1 = genCand.daughter(1);
     const Candidate * dauGen2 = genCand.daughter(2);
     cout << ">>> daughter MC 0 PDG Id " << dauGen0->pdgId()
          << ", status " << dauGen0->status()
          << ", charge " << dauGen0->charge()
          << endl;
     cout << ">>> daughter MC 1 PDG Id " << dauGen1->pdgId()
          << ", status " << dauGen1->status()
          << ", charge " << dauGen1->charge()
          << endl;
     cout << ">>> daughter MC 2 PDG Id " << dauGen2->pdgId()
          << ", status " << dauGen2->status()
          << ", charge " << dauGen2->charge() << endl;
       }
     }
   }
 }