#include <AcceptanceHistoProducer.h>

Inheritance diagram for AcceptanceHistoProducer:

Public Member Functions
	AcceptanceHistoProducer (const edm::ParameterSet &cfg)

virtual	~AcceptanceHistoProducer ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

void	beginJob ()

Private Attributes
DQMStore *	dbe_

std::string	dqmDir_

MonitorElement *	hEtaPosEtaNeg_

MonitorElement *	hPtNegEtaNeg_

MonitorElement *	hPtPosEtaPos_

MonitorElement *	hPtPosPtNeg_

edm::InputTag	srcGenParticles_

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

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::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

Produces di-muon acceptance histograms and stores them in lumi section.

Author: Armin Burgmeier, DESY

Version

Revision:: 1.2

Id:: AcceptanceHistoProducer.h,v 1.2 2012/11/07 17:03:01 aburgmei Exp

Definition at line 25 of file AcceptanceHistoProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 9 of file AcceptanceHistoProducer.cc.

                                                                           :
   dqmDir_(cfg.getParameter<std::string>("dqmDir")),
   srcGenParticles_(cfg.getParameter<edm::InputTag>("srcGenParticles")),
   hPtPosPtNeg_(NULL), hEtaPosEtaNeg_(NULL), hPtPosEtaPos_(NULL), hPtNegEtaNeg_(NULL)
 {
 }

AcceptanceHistoProducer::~AcceptanceHistoProducer ( )

virtual

Definition at line 16 of file AcceptanceHistoProducer.cc.

17 {

18 }

Member Function Documentation

void AcceptanceHistoProducer::analyze	(	const edm::Event &	evt,
		const edm::EventSetup &	es
	)

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 40 of file AcceptanceHistoProducer.cc.

References funct::abs(), reco::LeafCandidate::charge(), reco::CompositeRefCandidateT< D >::daughter(), reco::LeafCandidate::eta(), MonitorElement::Fill(), genParticleCandidates2GenParticles_cfi::genParticles, edm::Event::getByLabel(), hEtaPosEtaNeg_, hPtNegEtaNeg_, hPtPosEtaPos_, hPtPosPtNeg_, i, NULL, reco::CompositeRefCandidateT< D >::numberOfDaughters(), reco::LeafCandidate::pdgId(), reco::LeafCandidate::pt(), srcGenParticles_, and reco::LeafCandidate::status().

 {
   assert(hPtPosPtNeg_ != NULL);
   assert(hEtaPosEtaNeg_ != NULL);
   assert(hPtPosEtaPos_ != NULL);
   assert(hPtNegEtaNeg_ != NULL);
 
   edm::Handle<reco::GenParticleCollection> genParticles;
   evt.getByLabel(srcGenParticles_, genParticles);
 
   // Look for the two muons from the matrix element (status=3)
   const reco::GenParticle* genPosMuonME = NULL;
   const reco::GenParticle* genNegMuonME = NULL;
   for(unsigned int i = 0; i < genParticles->size(); ++i)
   {
     const reco::GenParticle& part = (*genParticles)[i];
     if(abs(part.pdgId()) == 13 && part.status() == 3)
     {
       if(part.charge() < 0)
         if(!genNegMuonME)
           genNegMuonME = &part;
       if(part.charge() > 0)
         if(!genPosMuonME)
           genPosMuonME = &part;
     }
   }
 
   if(genPosMuonME && genNegMuonME)
   {
     const reco::GenParticle* genPosMuon = genPosMuonME;
     const reco::GenParticle* genNegMuon = genNegMuonME;
 
     // Follow the decay chain to find the stable decay products (status=1)
     while(genPosMuon->status() != 1 && genPosMuon->numberOfDaughters() > 0)
     {
       unsigned int i;
       for(i = 0; i < genPosMuon->numberOfDaughters(); ++i)
       {
         const reco::GenParticle* daughter = dynamic_cast<const reco::GenParticle*>(genPosMuon->daughter(i));
         if(abs(daughter->pdgId()) == 13)
           { genPosMuon = daughter; break; }
       }
 
       // No more muon daugthers? Maybe mu->e decay...
       if(i == genPosMuon->numberOfDaughters())
         break;
     }
 
     while(genNegMuon->status() != 1 && genNegMuon->numberOfDaughters() > 0)
     {
       unsigned int i;
       for(i = 0; i < genNegMuon->numberOfDaughters(); ++i)
       {
         const reco::GenParticle* daughter = dynamic_cast<const reco::GenParticle*>(genNegMuon->daughter(i));
         if(abs(daughter->pdgId()) == 13)
           { genNegMuon = daughter; break; }
       }
 
       // No more muon daugthers? Maybe mu->e decay...
       if(i == genNegMuon->numberOfDaughters())
         break;
     }
 
     // Loop might have broken earlier if there are no more daugthers, or if
     // the decay chain does not end in a muon.
     if(genPosMuon->status() == 1 && genNegMuon->status() == 1)
     {
       hPtPosPtNeg_->Fill(genPosMuon->pt(), genNegMuon->pt());
       hEtaPosEtaNeg_->Fill(genPosMuon->eta(), genNegMuon->eta());
       hPtPosEtaPos_->Fill(genPosMuon->pt(), genPosMuon->eta());
       hPtNegEtaNeg_->Fill(genNegMuon->pt(), genNegMuon->eta());
     }
   }
 }

void AcceptanceHistoProducer::beginJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 20 of file AcceptanceHistoProducer.cc.

References DQMStore::book2DD(), dbe_, dqmDir_, hEtaPosEtaNeg_, hPtNegEtaNeg_, hPtPosEtaPos_, hPtPosPtNeg_, cppFunctionSkipper::operator, DQMStore::setCurrentFolder(), MonitorElement::setLumiFlag(), and MonitorElement::setResetMe().

 {
   dbe_ = edm::Service<DQMStore>().operator->();
 
   dbe_->setCurrentFolder("MCEmbedding/ZmumuAcceptance/" + dqmDir_);
   hPtPosPtNeg_ = dbe_->book2DD("PtPosPtNeg", "Positive muon transverse momentum vs. negative muon transverse momentum;p_{T}^{+};p_{T}^{-}", 500, 0.0, 500.0, 500, 0.0, 500.0);
   hEtaPosEtaNeg_ = dbe_->book2DD("EtaPosEtaNeg", "Positive muon pseudorapdity vs. negative muon pseudorapidity;#eta^{+};#eta^{-}", 500, -2.5, 2.5, 500, -2.5, 2.5);
   hPtPosEtaPos_ = dbe_->book2DD("PtPosEtaPos", "Positive muon transverse momentum vs. positive muon pseudorapidity;p_{T}^{+};#eta^{+}", 500, 0.0, 500.0, 500, -2.5, 2.5);
   hPtNegEtaNeg_ = dbe_->book2DD("PtNegEtaNeg", "Negative muon transverse momentum vs. negative muon pseudorapidity;p_{T}^{-};#eta^{-}", 500, 0.0, 500.0, 500, -2.5, 2.5);
 
   hPtPosPtNeg_->setResetMe(true);
   hPtPosPtNeg_->setLumiFlag();
   hEtaPosEtaNeg_->setResetMe(true);
   hEtaPosEtaNeg_->setLumiFlag();
   hPtPosEtaPos_->setResetMe(true);
   hPtPosEtaPos_->setLumiFlag();
   hPtNegEtaNeg_->setResetMe(true);
   hPtNegEtaNeg_->setLumiFlag();
 }

Member Data Documentation

DQMStore* AcceptanceHistoProducer::dbe_

private

Definition at line 42 of file AcceptanceHistoProducer.h.

Referenced by beginJob().

std::string AcceptanceHistoProducer::dqmDir_

private

Definition at line 39 of file AcceptanceHistoProducer.h.

Referenced by beginJob().

MonitorElement* AcceptanceHistoProducer::hEtaPosEtaNeg_

private

Definition at line 45 of file AcceptanceHistoProducer.h.

Referenced by analyze(), and beginJob().

MonitorElement* AcceptanceHistoProducer::hPtNegEtaNeg_

private

Definition at line 47 of file AcceptanceHistoProducer.h.

Referenced by analyze(), and beginJob().

MonitorElement* AcceptanceHistoProducer::hPtPosEtaPos_

private

Definition at line 46 of file AcceptanceHistoProducer.h.

Referenced by analyze(), and beginJob().

MonitorElement* AcceptanceHistoProducer::hPtPosPtNeg_

private

Definition at line 44 of file AcceptanceHistoProducer.h.

Referenced by analyze(), and beginJob().

edm::InputTag AcceptanceHistoProducer::srcGenParticles_

private

Definition at line 40 of file AcceptanceHistoProducer.h.

Referenced by analyze().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation