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

Inheritance diagram for MCvsRecoVerticesAnalyzer:

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

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

TH1F *	m_hclosestvtx

TH2F *	m_hclosestvtxvsnpu

TH1F *	m_hdeltazclose

TH2F *	m_hdeltazclosevsnpu

TH1F *	m_hdeltazfirst

TH2F *	m_hdeltazfirstvsnpu

edm::EDGetTokenT< edm::HepMCProduct >	m_hepMCProductToken

const edm::ParameterSet	m_histoParameters

TH2F *	m_hrecovsmclumi2d

TProfile *	m_hrecovsmclumiprof

TProfile *	m_hrecovsmclumiweightedprof

TH2F *	m_hrecovsmcnvtx2d

TProfile *	m_hrecovsmcnvtxprof

TProfile *	m_hrecovsmcnvtxweightedprof

edm::EDGetTokenT< reco::VertexCollection >	m_recoVertexCollectionToken

const bool	m_useVisibleVertices

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

Constructor & Destructor Documentation

◆ MCvsRecoVerticesAnalyzer()

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

explicit

Definition at line 101 of file MCvsRecoVerticesAnalyzer.cc.

References edm::ParameterSet::getUntrackedParameter(), m_hclosestvtx, m_hclosestvtxvsnpu, m_hdeltazclose, m_hdeltazclosevsnpu, m_hdeltazfirst, m_hdeltazfirstvsnpu, m_histoParameters, m_hrecovsmclumi2d, m_hrecovsmclumiprof, m_hrecovsmclumiweightedprof, m_hrecovsmcnvtx2d, m_hrecovsmcnvtxprof, m_hrecovsmcnvtxweightedprof, m_useVisibleVertices, m_useweight, and TFileService::make().

     : m_useweight(iConfig.getParameter<bool>("useWeight")),
       m_useVisibleVertices(iConfig.getParameter<bool>("useVisibleVertices")),
       m_histoParameters(iConfig.getUntrackedParameter<edm::ParameterSet>("histoParameters", edm::ParameterSet())),
       m_doubleToken(consumes<double>(iConfig.getParameter<edm::InputTag>("weightProduct"))),
       m_vecPileupSummaryInfoToken(
           consumes<std::vector<PileupSummaryInfo> >(iConfig.getParameter<edm::InputTag>("pileupSummaryCollection"))),
       m_recoVertexCollectionToken(
           consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("pvCollection"))),
       m_hepMCProductToken(consumes<edm::HepMCProduct>(iConfig.getParameter<edm::InputTag>("mcTruthCollection"))) {
   //now do what ever initialization is needed
 
   if (m_useVisibleVertices)
     edm::LogInfo("UseVisibleVertices") << "Only visible vertices will be used to compute Npileup";
 
   usesResource("TFileService");
   edm::Service<TFileService> tfserv;
 
   m_hrecovsmcnvtx2d = tfserv->make<TH2F>(
       "recovsmcnvtx2d", "Number of reco vertices vs pileup interactions", 60, -0.5, 59.5, 60, -0.5, 59.5);
   m_hrecovsmcnvtx2d->GetXaxis()->SetTitle("Pileup Interactions");
   m_hrecovsmcnvtx2d->GetYaxis()->SetTitle("Reco Vertices");
   m_hrecovsmcnvtxprof =
       tfserv->make<TProfile>("recovsmcnvtxprof", "Mean number of reco vs pileup vertices", 60, -0.5, 59.5);
   m_hrecovsmcnvtxprof->GetXaxis()->SetTitle("Pileup Interactions");
   m_hrecovsmcnvtxprof->GetYaxis()->SetTitle("Reco Vertices");
 
   m_hrecovsmclumi2d = tfserv->make<TH2F>(
       "recovsmclumi2d", "Number of reco vertices vs ave pileup interactions", 200, 0., 50., 60, -0.5, 59.5);
   m_hrecovsmclumi2d->GetXaxis()->SetTitle("Average Pileup Interactions");
   m_hrecovsmclumi2d->GetYaxis()->SetTitle("Reco Vertices");
   m_hrecovsmclumiprof =
       tfserv->make<TProfile>("recovsmclumiprof", "Mean number of reco vs ave pileup vertices", 200, 0., 50.);
   m_hrecovsmclumiprof->GetXaxis()->SetTitle("Average Pileup Interactions");
   m_hrecovsmclumiprof->GetYaxis()->SetTitle("Reco Vertices");
 
   if (m_useweight) {
     m_hrecovsmcnvtxweightedprof = tfserv->make<TProfile>(
         "recovsmcnvtxweightedprof", "Mean number of reco vs pileup vertices (1-w) weight", 60, -0.5, 59.5);
     m_hrecovsmcnvtxweightedprof->GetXaxis()->SetTitle("Pileup Interactions");
     m_hrecovsmcnvtxweightedprof->GetYaxis()->SetTitle("Reco Vertices (1-w)");
 
     m_hrecovsmclumiweightedprof = tfserv->make<TProfile>(
         "recovsmclumiweightedprof", "Mean number of reco vs ave pileup vertices (1-w) weight", 200, 0., 50.);
     m_hrecovsmclumiweightedprof->GetXaxis()->SetTitle("Average Pileup Interactions");
     m_hrecovsmclumiweightedprof->GetYaxis()->SetTitle("Reco Vertices (1-w)");
   }
 
   m_hdeltazfirst = tfserv->make<TH1F>("deltazfirst",
                                       "Reco-MC vertex z position (first vertex)",
                                       m_histoParameters.getUntrackedParameter<unsigned int>("zBins", 1000),
                                       m_histoParameters.getUntrackedParameter<double>("zMin", -1.),
                                       m_histoParameters.getUntrackedParameter<double>("zMax", 1.));
   m_hdeltazfirst->GetXaxis()->SetTitle("#Delta z (cm)");
   m_hdeltazfirst->GetYaxis()->SetTitle("Events");
 
   m_hdeltazclose = tfserv->make<TH1F>("deltazclose",
                                       "Reco-MC vertex z position (closest vertex)",
                                       m_histoParameters.getUntrackedParameter<unsigned int>("zBins", 1000),
                                       m_histoParameters.getUntrackedParameter<double>("zMin", -1.),
                                       m_histoParameters.getUntrackedParameter<double>("zMax", 1.));
   m_hdeltazclose->GetXaxis()->SetTitle("#Delta z (cm)");
   m_hdeltazclose->GetYaxis()->SetTitle("Events");
 
   m_hclosestvtx = tfserv->make<TH1F>("closestvtx", "Closest reco vtx ID", 30, -0.5, 29.5);
   m_hclosestvtx->GetXaxis()->SetTitle("Vtx ID");
   m_hclosestvtx->GetYaxis()->SetTitle("Events");
 
   m_hdeltazfirstvsnpu = tfserv->make<TH2F>("deltazfirstvsnpu",
                                            "Reco-MC vertex z position (first vertex) vs Npileup",
                                            30,
                                            -0.5,
                                            29.5,
                                            m_histoParameters.getUntrackedParameter<unsigned int>("zBins", 1000),
                                            m_histoParameters.getUntrackedParameter<double>("zMin", -1.),
                                            m_histoParameters.getUntrackedParameter<double>("zMax", 1.));
   m_hdeltazfirstvsnpu->GetXaxis()->SetTitle("pileup Interactions");
   m_hdeltazfirstvsnpu->GetYaxis()->SetTitle("#Delta z (cm)");
 
   m_hdeltazclosevsnpu = tfserv->make<TH2F>("deltazclosevsnpu",
                                            "Reco-MC vertex z position (closest vertex) v Npileup",
                                            30,
                                            -0.5,
                                            29.5,
                                            m_histoParameters.getUntrackedParameter<unsigned int>("zBins", 1000),
                                            m_histoParameters.getUntrackedParameter<double>("zMin", -1.),
                                            m_histoParameters.getUntrackedParameter<double>("zMax", 1.));
   m_hdeltazclosevsnpu->GetXaxis()->SetTitle("Pileup Interactions");
   m_hdeltazclosevsnpu->GetYaxis()->SetTitle("#Delta z (cm)");
 
   m_hclosestvtxvsnpu =
       tfserv->make<TH2F>("closestvtxvsnpu", "Closest reco vtx ID vs Npileup", 30, -0.5, 29.5, 30, -0.5, 29.5);
   m_hclosestvtxvsnpu->GetXaxis()->SetTitle("Pileup Interactions");
   m_hclosestvtxvsnpu->GetYaxis()->SetTitle("Vtx ID");
 }

◆ ~MCvsRecoVerticesAnalyzer()

MCvsRecoVerticesAnalyzer::~MCvsRecoVerticesAnalyzer ( )

override

Definition at line 197 of file MCvsRecoVerticesAnalyzer.cc.

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

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 207 of file MCvsRecoVerticesAnalyzer.cc.

References funct::abs(), edm::HepMCProduct::GetEvent(), iEvent, m_doubleToken, m_hclosestvtx, m_hclosestvtxvsnpu, m_hdeltazclose, m_hdeltazclosevsnpu, m_hdeltazfirst, m_hdeltazfirstvsnpu, m_hepMCProductToken, m_hrecovsmclumi2d, m_hrecovsmclumiprof, m_hrecovsmclumiweightedprof, m_hrecovsmcnvtx2d, m_hrecovsmcnvtxprof, m_hrecovsmcnvtxweightedprof, m_recoVertexCollectionToken, m_useVisibleVertices, m_useweight, m_vecPileupSummaryInfoToken, mps_merge::weight, and z.

                                                                                           {
   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);
 
   // look for the intime PileupSummaryInfo
 
   std::vector<PileupSummaryInfo>::const_iterator pileupinfo;
 
   for (pileupinfo = pileupinfos->begin(); pileupinfo != pileupinfos->end(); ++pileupinfo) {
     if (pileupinfo->getBunchCrossing() == 0)
       break;
   }
 
   //
 
   edm::Handle<reco::VertexCollection> pvcoll;
   iEvent.getByToken(m_recoVertexCollectionToken, pvcoll);
 
   //
 
   if (pileupinfo->getBunchCrossing() != 0) {
     edm::LogError("NoInTimePileUpInfo") << "Cannot find the in-time pileup info " << pileupinfo->getBunchCrossing();
 
   } else {
     int npileup = pileupinfo->getPU_NumInteractions();
 
     if (m_useVisibleVertices)
       npileup = pileupinfo->getPU_zpositions().size();
 
     m_hrecovsmcnvtx2d->Fill(npileup, pvcoll->size(), weight);
     m_hrecovsmcnvtxprof->Fill(npileup, pvcoll->size(), weight);
 
     m_hrecovsmclumi2d->Fill(pileupinfo->getTrueNumInteractions(), pvcoll->size(), weight);
     m_hrecovsmclumiprof->Fill(pileupinfo->getTrueNumInteractions(), pvcoll->size(), weight);
 
     if (m_useweight) {
       m_hrecovsmcnvtxweightedprof->Fill(npileup, pvcoll->size(), 1. - weight);
       m_hrecovsmclumiweightedprof->Fill(pileupinfo->getTrueNumInteractions(), pvcoll->size(), 1. - weight);
     }
     //
 
     edm::Handle<edm::HepMCProduct> EvtHandle;
     iEvent.getByToken(m_hepMCProductToken, EvtHandle);
 
     const HepMC::GenEvent* Evt = EvtHandle->GetEvent();
 
     // compute the difference between the main interaction vertex z position and the first vertex of the collection
 
     if (!pvcoll->empty()) {
       if (!(*pvcoll)[0].isFake()) {
         // get the first vertex
         if (Evt->vertices_begin() != Evt->vertices_end()) {
           m_hdeltazfirst->Fill((*pvcoll)[0].z() - (*Evt->vertices_begin())->point3d().z() / 10., weight);
           m_hdeltazfirstvsnpu->Fill(npileup, (*pvcoll)[0].z() - (*Evt->vertices_begin())->point3d().z() / 10., weight);
         }
       }
     }
 
     // compute the difference between the main interaction vertex z position and the closest reco vertex
 
     double minabsdist = -1.;
     double mindist = -999.;
     int closestvtx = -1;
 
     for (unsigned int ivtx = 0; ivtx < pvcoll->size(); ++ivtx) {
       if (closestvtx < 0 ||
           minabsdist > std::abs((*pvcoll)[ivtx].z() - (*Evt->vertices_begin())->point3d().z() / 10.)) {
         mindist = (*pvcoll)[ivtx].z() - (*Evt->vertices_begin())->point3d().z() / 10.;
         closestvtx = ivtx;
         minabsdist = std::abs(mindist);
       }
     }
     if (closestvtx >= 0) {
       m_hdeltazclose->Fill(mindist, weight);
       m_hdeltazclosevsnpu->Fill(npileup, mindist, weight);
       m_hclosestvtx->Fill(closestvtx, weight);
       m_hclosestvtxvsnpu->Fill(npileup, closestvtx, weight);
     }
   }
 }