#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 noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept 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 noexcept

bool	wantsStreamRuns () const noexcept

	~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);
     }
   }
 }