MCVerticesAnalyzer Class Reference

#include <TrackingPFG/PileUp/src/MCVerticesAnalyzer.cc>

Inheritance diagram for MCVerticesAnalyzer:

Public Member Functions
	MCVerticesAnalyzer (const edm::ParameterSet &)
	~MCVerticesAnalyzer ()
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
std::vector< ConsumesInfo >	consumesInfo () const
	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
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
virtual void	beginRun (const edm::Run &, const edm::EventSetup &)
virtual void	endJob ()
virtual void	endRun (const edm::Run &, const edm::EventSetup &)

Private Attributes
edm::EDGetTokenT< double >	m_doubleToken
edm::EDGetTokenT < edm::HepMCProduct >	m_hepMCProductToken
TH1F *	m_hlumi
TH1F *	m_hmainvtxx
TH1F *	m_hmainvtxy
TH1F *	m_hmainvtxz
TH1F *	m_hnvtx
TH2F *	m_hnvtxvsbx
TH2F *	m_hnvtxvslumi
TH1F *	m_hnvtxweight
TProfile *	m_hnvtxweightprof
TH1F *	m_hpileupvtxz
const bool	m_useweight
edm::EDGetTokenT< std::vector < PileupSummaryInfo > >	m_vecPileupSummaryInfoToken

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

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 57 of file MCVerticesAnalyzer.cc.

Constructor & Destructor Documentation

MCVerticesAnalyzer::MCVerticesAnalyzer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 103 of file MCVerticesAnalyzer.cc.

References m_hlumi, m_hmainvtxx, m_hmainvtxy, m_hmainvtxz, m_hnvtx, m_hnvtxvsbx, m_hnvtxvslumi, m_hnvtxweight, m_hnvtxweightprof, m_hpileupvtxz, m_useweight, and TFileService::make().

  : m_useweight( iConfig.getParameter< bool >( "useWeight" ) )
  , m_doubleToken( consumes< double >( iConfig.getParameter< edm::InputTag >( "weightProduct" ) ) )
  , m_vecPileupSummaryInfoToken( consumes< std::vector<PileupSummaryInfo> >( iConfig.getParameter< edm::InputTag >( "pileupSummaryCollection" ) ) )
  , m_hepMCProductToken( consumes< edm::HepMCProduct >( iConfig.getParameter< edm::InputTag >( "mcTruthCollection" ) ) )
{
   //now do what ever initialization is needed



  edm::Service<TFileService> tfserv;

  m_hnvtx = tfserv->make<TH1F>("nvtx","Number of pileup vertices",60,-0.5,59.5);
  m_hnvtx->GetXaxis()->SetTitle("Number of Interactions");

  m_hnvtxvsbx = tfserv->make<TH2F>("nvtxvsbx","Number of pileup vertices vs BX",9,-4.5,4.5,60,-0.5,59.5);
  m_hnvtxvsbx->GetXaxis()->SetTitle("BX number");
  m_hnvtxvsbx->GetYaxis()->SetTitle("Number of Interactions");

  m_hlumi = tfserv->make<TH1F>("lumi","BX luminosity*xsect",200,0.,50.);
  m_hlumi->GetXaxis()->SetTitle("Average Number of Interactions");

  m_hnvtxvslumi = tfserv->make<TH2F>("nvtxvslumi","Npileup vs BX luminosity*xsect",200,0.,50.,60,-0.5,59.5);
  m_hnvtxvslumi->GetXaxis()->SetTitle("Average Number of Interactions");  m_hnvtxvslumi->GetYaxis()->SetTitle("Number of Interactions");

  if(m_useweight) {
    m_hnvtxweight = tfserv->make<TH1F>("nvtxweight","Number of pileup vertices (1-w)",60,-0.5,59.5);
    m_hnvtxweight->GetXaxis()->SetTitle("Number of Interactions");
    m_hnvtxweightprof = tfserv->make<TProfile>("nvtxweightprof","Mean (1-w) vs Number of pileup interactions",60,-0.5,59.5);
    m_hnvtxweightprof->GetXaxis()->SetTitle("Number of Interactions");
  }

  m_hmainvtxx = tfserv->make<TH1F>("mainvtxx","Main vertex x position",200,-.5,.5);
  m_hmainvtxx->GetXaxis()->SetTitle("X (cm)");
  m_hmainvtxy = tfserv->make<TH1F>("mainvtxy","Main vertex y position",200,-.5,.5);
  m_hmainvtxy->GetXaxis()->SetTitle("Y (cm)");
  m_hmainvtxz = tfserv->make<TH1F>("mainvtxz","Main vertex z position",600,-30.,30.);
  m_hmainvtxz->GetXaxis()->SetTitle("Z (cm)");
  m_hpileupvtxz = tfserv->make<TH1F>("pileupvtxz","PileUp vertices z position",600,-30.,30.);
  m_hpileupvtxz->GetXaxis()->SetTitle("Z (cm)");

}

MCVerticesAnalyzer::~MCVerticesAnalyzer

(

)

Definition at line 147 of file MCVerticesAnalyzer.cc.

{

   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)

}

Member Function Documentation

void MCVerticesAnalyzer::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 162 of file MCVerticesAnalyzer.cc.

References edm::Event::getByToken(), edm::HandleBase::isValid(), m_doubleToken, m_hepMCProductToken, m_hlumi, m_hmainvtxx, m_hmainvtxy, m_hmainvtxz, m_hnvtx, m_hnvtxvsbx, m_hnvtxvslumi, m_hnvtxweight, m_hnvtxweightprof, m_hpileupvtxz, m_useweight, m_vecPileupSummaryInfoToken, and histoStyle::weight.

{

   double weight = 1.;

   if(m_useweight) {
     edm::Handle<double> weightprod;
     iEvent.getByToken( m_doubleToken, weightprod );

     weight = *weightprod;

   }


   edm::Handle<std::vector<PileupSummaryInfo> >  pileupinfos;
   iEvent.getByToken( m_vecPileupSummaryInfoToken, pileupinfos );

   //

   if(pileupinfos.isValid()) {

  // look for the intime PileupSummaryInfo

     std::vector<PileupSummaryInfo>::const_iterator pileupinfo;

     for(pileupinfo = pileupinfos->begin(); pileupinfo != pileupinfos->end() ; ++pileupinfo) {
       m_hnvtxvsbx->Fill(pileupinfo->getBunchCrossing(),pileupinfo->getPU_NumInteractions(),weight);
     } 


     for(pileupinfo = pileupinfos->begin(); pileupinfo != pileupinfos->end() ; ++pileupinfo) {
       if(pileupinfo->getBunchCrossing()==0) break;
     } 
 
     //

     if(pileupinfo->getBunchCrossing()!=0) {
       edm::LogError("NoInTimePileUpInfo") << "Cannot find the in-time pileup info " << pileupinfo->getBunchCrossing();
     }
     else {

       m_hlumi->Fill(pileupinfo->getTrueNumInteractions(),weight);
       m_hnvtx->Fill(pileupinfo->getPU_NumInteractions(),weight);
       m_hnvtxvslumi->Fill(pileupinfo->getTrueNumInteractions(),pileupinfo->getPU_NumInteractions(),weight);

       if(m_useweight) {
     m_hnvtxweight->Fill(pileupinfo->getPU_NumInteractions(),1.-weight);
     m_hnvtxweightprof->Fill(pileupinfo->getPU_NumInteractions(),1.-weight);
       }
       
       const std::vector<float>& zpositions = pileupinfo->getPU_zpositions();
       
       for(std::vector<float>::const_iterator zpos = zpositions.begin() ; zpos != zpositions.end() ; ++zpos) {
     
     m_hpileupvtxz->Fill(*zpos,weight);
     
       }
     }
   }
   // main interaction part

   edm::Handle< edm::HepMCProduct > EvtHandle ;
   iEvent.getByToken( m_hepMCProductToken, EvtHandle );

   if(EvtHandle.isValid()) {

     const HepMC::GenEvent* Evt = EvtHandle->GetEvent();

   // get the first vertex
     
     if(Evt->vertices_begin() != Evt->vertices_end()) {

       m_hmainvtxx->Fill((*Evt->vertices_begin())->point3d().x()/10.,weight);
       m_hmainvtxy->Fill((*Evt->vertices_begin())->point3d().y()/10.,weight);
       m_hmainvtxz->Fill((*Evt->vertices_begin())->point3d().z()/10.,weight);

     }
   }
}

void MCVerticesAnalyzer::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 256 of file MCVerticesAnalyzer.cc.

{
}

void MCVerticesAnalyzer::beginRun ( const edm::Run &  iRun, const edm::EventSetup &    )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 243 of file MCVerticesAnalyzer.cc.

{
}

void MCVerticesAnalyzer::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 262 of file MCVerticesAnalyzer.cc.

{
}

void MCVerticesAnalyzer::endRun ( const edm::Run &  iRun, const edm::EventSetup &    )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 248 of file MCVerticesAnalyzer.cc.

{
}

Member Data Documentation

edm::EDGetTokenT< double > MCVerticesAnalyzer::m_doubleToken

private

Definition at line 75 of file MCVerticesAnalyzer.cc.

Referenced by analyze().

edm::EDGetTokenT< edm::HepMCProduct > MCVerticesAnalyzer::m_hepMCProductToken

private

Definition at line 77 of file MCVerticesAnalyzer.cc.

Referenced by analyze().

TH1F* MCVerticesAnalyzer::m_hlumi

private

Definition at line 81 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

TH1F* MCVerticesAnalyzer::m_hmainvtxx

private

Definition at line 85 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

TH1F* MCVerticesAnalyzer::m_hmainvtxy

private

Definition at line 86 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

TH1F* MCVerticesAnalyzer::m_hmainvtxz

private

Definition at line 87 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

TH1F* MCVerticesAnalyzer::m_hnvtx

private

Definition at line 79 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

TH2F* MCVerticesAnalyzer::m_hnvtxvsbx

private

Definition at line 80 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

TH2F* MCVerticesAnalyzer::m_hnvtxvslumi

private

Definition at line 82 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

TH1F* MCVerticesAnalyzer::m_hnvtxweight

private

Definition at line 83 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

TProfile* MCVerticesAnalyzer::m_hnvtxweightprof

private

Definition at line 84 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

TH1F* MCVerticesAnalyzer::m_hpileupvtxz

private

Definition at line 88 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

const bool MCVerticesAnalyzer::m_useweight

private

Definition at line 73 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

edm::EDGetTokenT< std::vector<PileupSummaryInfo> > MCVerticesAnalyzer::m_vecPileupSummaryInfoToken

private

Definition at line 76 of file MCVerticesAnalyzer.cc.

Referenced by analyze().