ZMassHistogrammer Class Reference

Inheritance diagram for ZMassHistogrammer:

Public Member Functions
	ZMassHistogrammer (const edm::ParameterSet &pset)
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

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

Private Attributes
edm::InputTag	gen_
TH1F *	h_mZ_
TH1F *	h_mZMC_
edm::InputTag	z_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
CurrentProcessingContext const *	currentContext () const
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 5 of file ZMassHistogrammer.cc.

Constructor & Destructor Documentation

ZMassHistogrammer::ZMassHistogrammer

(

const edm::ParameterSet &

pset

)

Definition at line 29 of file ZMassHistogrammer.cc.

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

                                                             :
  z_(pset.getParameter<InputTag>("z")),
  gen_(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  )

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 38 of file ZMassHistogrammer.cc.

References reco::Candidate::charge(), gather_cfg::cout, reco::Candidate::daughter(), relval_steps::gen(), gen_, h_mZ_, h_mZMC_, i, reco::Candidate::mass(), reco::Candidate::numberOfDaughters(), reco::Candidate::pdgId(), reco::Candidate::status(), detailsBasic3DVector::z, and z_.

                                                                                 { 
  cout << ">>> Z Mass analyze" << endl;
  Handle<CandidateCollection> z;
  Handle<CandidateCollection> gen;
  event.getByLabel(z_, z);
  event.getByLabel(gen_, 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;
      }
    }
  }
}

Member Data Documentation

edm::InputTag ZMassHistogrammer::gen_

private

Definition at line 10 of file ZMassHistogrammer.cc.

Referenced by analyze().

TH1F* ZMassHistogrammer::h_mZ_

private

Definition at line 11 of file ZMassHistogrammer.cc.

Referenced by analyze(), and ZMassHistogrammer().

TH1F * ZMassHistogrammer::h_mZMC_

private

Definition at line 11 of file ZMassHistogrammer.cc.

Referenced by analyze(), and ZMassHistogrammer().

edm::InputTag ZMassHistogrammer::z_

private

Definition at line 10 of file ZMassHistogrammer.cc.

Referenced by analyze().