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

Inheritance diagram for MCVerticesAnalyzer:

Public Member Functions
	MCVerticesAnalyzer (const edm::ParameterSet &)

	~MCVerticesAnalyzer () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

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

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > 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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

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::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

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

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

Definition at line 54 of file MCVerticesAnalyzer.cc.

Constructor & Destructor Documentation

◆ MCVerticesAnalyzer()

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

explicit

Definition at line 93 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
 
   usesResource("TFileService");
   edm::Service<TFileService> tfserv;
 
   m_hnvtx = tfserv->make<TH1F>("nvtx", "Number of pileup vertices", 400, -0.5, 399.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, 400, -0.5, 399.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., 400, -0.5, 399.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)", 400, -0.5, 399.5);
     m_hnvtxweight->GetXaxis()->SetTitle("Number of Interactions");
     m_hnvtxweightprof =
         tfserv->make<TProfile>("nvtxweightprof", "Mean (1-w) vs Number of pileup interactions", 400, -0.5, 399.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::~MCVerticesAnalyzer ( )

override

Definition at line 136 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

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 146 of file MCVerticesAnalyzer.cc.

References edm::HepMCProduct::GetEvent(), iEvent, 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 mps_merge::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);
     }
   }
 }

Member Data Documentation

◆ m_doubleToken

edm::EDGetTokenT<double> MCVerticesAnalyzer::m_doubleToken

private

Definition at line 66 of file MCVerticesAnalyzer.cc.

Referenced by analyze().

◆ m_hepMCProductToken

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

private

Definition at line 68 of file MCVerticesAnalyzer.cc.

Referenced by analyze().

◆ m_hlumi

TH1F* MCVerticesAnalyzer::m_hlumi

private

Definition at line 72 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

◆ m_hmainvtxx

TH1F* MCVerticesAnalyzer::m_hmainvtxx

private

Definition at line 76 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

◆ m_hmainvtxy

TH1F* MCVerticesAnalyzer::m_hmainvtxy

private

Definition at line 77 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

◆ m_hmainvtxz

TH1F* MCVerticesAnalyzer::m_hmainvtxz

private

Definition at line 78 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

◆ m_hnvtx

TH1F* MCVerticesAnalyzer::m_hnvtx

private

Definition at line 70 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

◆ m_hnvtxvsbx

TH2F* MCVerticesAnalyzer::m_hnvtxvsbx

private

Definition at line 71 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

◆ m_hnvtxvslumi

TH2F* MCVerticesAnalyzer::m_hnvtxvslumi

private

Definition at line 73 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

◆ m_hnvtxweight

TH1F* MCVerticesAnalyzer::m_hnvtxweight

private

Definition at line 74 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

◆ m_hnvtxweightprof

TProfile* MCVerticesAnalyzer::m_hnvtxweightprof

private

Definition at line 75 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

◆ m_hpileupvtxz

TH1F* MCVerticesAnalyzer::m_hpileupvtxz

private

Definition at line 79 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

◆ m_useweight

const bool MCVerticesAnalyzer::m_useweight

private

Definition at line 64 of file MCVerticesAnalyzer.cc.

Referenced by analyze(), and MCVerticesAnalyzer().

◆ m_vecPileupSummaryInfoToken

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

private

Definition at line 67 of file MCVerticesAnalyzer.cc.

Referenced by analyze().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MCVerticesAnalyzer()

◆ ~MCVerticesAnalyzer()

Member Function Documentation

◆ analyze()

Member Data Documentation

◆ 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