Inheritance diagram for OffsetAnalyzerDQM:

Classes
struct	Plot1D

struct	PlotProfile

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

	OffsetAnalyzerDQM (const edm::ParameterSet &)

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

void	beginStream (edm::StreamID id) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Protected Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

int	getEtaIndex (float eta)

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Attributes
std::vector< double >	etabins

int	muHigh

edm::EDGetTokenT< edm::View< PileupSummaryInfo > >	muToken

int	npvHigh

std::string	offsetPlotBaseName

std::map< std::string, PlotProfile >	offsetPlots

std::map< int, std::string >	pdgMap

edm::EDGetTokenT< edm::View< pat::PackedCandidate > >	pfToken

std::vector< std::string >	pftypes

edm::EDGetTokenT< edm::View< reco::Vertex > >	pvToken

std::map< std::string, Plot1D >	th1dPlots

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 16 of file OffsetAnalyzerDQM.cc.

Constructor & Destructor Documentation

◆ OffsetAnalyzerDQM()

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

Definition at line 95 of file OffsetAnalyzerDQM.cc.

                                                                    {
   offsetPlotBaseName = iConfig.getParameter<std::string>("offsetPlotBaseName");
  
   pvToken = consumes<edm::View<reco::Vertex>>(iConfig.getParameter<edm::InputTag>("pvTag"));
   muToken = consumes<edm::View<PileupSummaryInfo>>(iConfig.getParameter<edm::InputTag>("muTag"));
   pfToken = consumes<edm::View<pat::PackedCandidate>>(iConfig.getParameter<edm::InputTag>("pfTag"));
  
   etabins = iConfig.getParameter<std::vector<double>>("etabins");
   pftypes = iConfig.getParameter<std::vector<std::string>>("pftypes");
  
   muHigh = iConfig.getUntrackedParameter<int>("muHigh");
   npvHigh = iConfig.getUntrackedParameter<int>("npvHigh");
  
   //initialize offset plots
   const auto& offset_psets = iConfig.getParameter<std::vector<edm::ParameterSet>>("offsetPlots");
   for (auto& pset : offset_psets) {
     std::string name = pset.getParameter<std::string>("name");
     std::string title = pset.getParameter<std::string>("title");
     std::string dir = pset.getParameter<std::string>("dir");
     std::vector<double> vx = pset.getParameter<std::vector<double>>("vx");
     int ny = pset.getParameter<uint32_t>("ny");
     double y0 = pset.getParameter<double>("y0");
     double y1 = pset.getParameter<double>("y1");
  
     offsetPlots[name] = PlotProfile(name, title, dir, 0, 0, 0, vx, ny, y0, y1);
   }
  
   //initialize th1d
   const auto& th1d_psets = iConfig.getParameter<std::vector<edm::ParameterSet>>("th1dPlots");
   for (auto& pset : th1d_psets) {
     std::string name = pset.getParameter<std::string>("name");
     std::string title = pset.getParameter<std::string>("title");
     std::string dir = pset.getParameter<std::string>("dir");
     int nx = pset.getParameter<uint32_t>("nx");
     double x0 = pset.getParameter<double>("x0");
     double x1 = pset.getParameter<double>("x1");
  
     th1dPlots[name] = Plot1D(name, title, dir, nx, x0, x1);
   }
  
   //create pdg map
   std::vector<uint32_t> pdgKeys = iConfig.getParameter<std::vector<uint32_t>>("pdgKeys");
   std::vector<std::string> pdgStrs = iConfig.getParameter<std::vector<std::string>>("pdgStrs");
   for (int i = 0, n = pdgKeys.size(); i < n; i++)
     pdgMap[pdgKeys[i]] = pdgStrs[i];
 }

References DeadROC_duringRun::dir, etabins, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, muHigh, muToken, dqmiodumpmetadata::n, Skims_PA_cff::name, npvHigh, offsetPlotBaseName, offsetPlots, offsetAnalyzerDQM_cff::pdgKeys, pdgMap, offsetAnalyzerDQM_cff::pdgStrs, pfToken, pftypes, Plot1D, muonDTDigis_cfi::pset, pvToken, AlCaHLTBitMon_QueryRunRegistry::string, th1dPlots, overlapproblemtsosanalyzer_cfi::title, testProducerWithPsetDescEmpty_cfi::x1, and testProducerWithPsetDescEmpty_cfi::y1.

Member Function Documentation

◆ analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 152 of file OffsetAnalyzerDQM.cc.

                                                                             {
   //npv//
   edm::Handle<edm::View<reco::Vertex>> vertexHandle;
   iEvent.getByToken(pvToken, vertexHandle);
  
   unsigned int nPVall = vertexHandle->size();
   bool isGoodPV[nPVall];
   for (size_t i = 0; i < nPVall; ++i)
     isGoodPV[i] = false;
  
   int npv = 0;
   for (unsigned int i = 0; i < nPVall; i++) {
     const auto& pv = vertexHandle->at(i);
  
     if (!pv.isFake() && pv.ndof() >= 4 && fabs(pv.z()) <= 24.0 && fabs(pv.position().rho()) <= 2.0) {
       npv++;
       isGoodPV[i] = true;
     }
   }
   th1dPlots["npv"].fill(npv);
   int npv_in_range = npv;
   if (npv_in_range < 0)
     npv_in_range = 0;
   else if (npv_in_range >= npvHigh)
     npv_in_range = npvHigh - 1;  // make sure int_mu won't lead to non-existing ME
  
   //mu//
   int int_mu = -1;
   edm::Handle<edm::View<PileupSummaryInfo>> muHandle;
   if (iEvent.getByToken(muToken, muHandle)) {
     const auto& summary = *muHandle;
     auto it = std::find_if(summary.begin(), summary.end(), [](const auto& s) { return s.getBunchCrossing() == 0; });
  
     if (it->getBunchCrossing() != 0) {
       edm::LogError("OffsetAnalyzerDQM") << "Cannot find the in-time pileup info " << it->getBunchCrossing();
     } else {
       float mu = it->getTrueNumInteractions();
       th1dPlots["mu"].fill(mu);
       int_mu = mu + 0.5;
     }
   }
   if (int_mu >= muHigh)
     int_mu = muHigh - 1;  // make sure int_mu won't lead to non-existing ME
  
   //create map of pftypes vs total energy / eta
   std::map<std::string, std::vector<double>> m_pftype_etaE;
   int nEta = etabins.size() - 1;
   for (const auto& pftype : pftypes)
     m_pftype_etaE[pftype].assign(nEta, 0.0);
  
   //pf particles//
   edm::Handle<edm::View<pat::PackedCandidate>> pfHandle;
   iEvent.getByToken(pfToken, pfHandle);
  
   for (unsigned int i = 0, n = pfHandle->size(); i < n; i++) {
     const auto& cand = pfHandle->at(i);
  
     int etaIndex = getEtaIndex(cand.eta());
     std::string pftype = pdgMap[abs(cand.pdgId())];
     if (etaIndex == -1 || pftype.empty())
       continue;
  
     if (pftype == "chm") {  //check charged hadrons ONLY
       bool attached = false;
  
       for (unsigned int ipv = 0; ipv < nPVall && !attached; ipv++) {
         if (isGoodPV[ipv] && cand.fromPV(ipv) == 3)
           attached = true;  //pv used in fit
       }
       if (!attached)
         pftype = "chu";  //unmatched charged hadron
     }
     /*
         reco::TrackRef candTrkRef( cand.trackRef() );
         if ( pftype == "chm" && !candTrkRef.isNull() ) { //check charged hadrons ONLY
             bool attached = false;
  
             for (auto ipv=vertexHandle->begin(), endpv=vertexHandle->end(); ipv != endpv && !attached; ++ipv) {
                 if ( !ipv->isFake() && ipv->ndof() >= 4 && fabs(ipv->z()) < 24 ) { //must be attached to a good pv
  
                     for(auto ivtrk=ipv->tracks_begin(), endvtrk=ipv->tracks_end(); ivtrk != endvtrk && !attached; ++ivtrk) {
                         reco::TrackRef pvTrkRef(ivtrk->castTo<reco::TrackRef>());
                         if (pvTrkRef == candTrkRef) attached = true;
                     }
                 }
             }
             if (!attached) pftype = "chu"; //unmatched charged hadron
         }
 */
     m_pftype_etaE[pftype][etaIndex] += cand.et();
   }
  
   for (const auto& pair : m_pftype_etaE) {
     std::string pftype = pair.first;
     std::vector<double> etaE = pair.second;
  
     std::string offset_name_npv = offsetPlotBaseName + "_npv" + std::to_string(npv_in_range) + "_" + pftype;
     if (offsetPlots.find(offset_name_npv) == offsetPlots.end())
       return;  //npv is out of range ()
  
     for (int i = 0; i < nEta; i++) {
       double eta = 0.5 * (etabins[i] + etabins[i + 1]);
       offsetPlots[offset_name_npv].fill2D(eta, etaE[i]);
     }
  
     if (int_mu != -1) {
       std::string offset_name_mu = offsetPlotBaseName + "_mu" + std::to_string(int_mu) + "_" + pftype;
       if (offsetPlots.find(offset_name_mu) == offsetPlots.end())
         return;  //mu is out of range
  
       for (int i = 0; i < nEta; i++) {
         double eta = 0.5 * (etabins[i] + etabins[i + 1]);
         offsetPlots[offset_name_mu].fill2D(eta, etaE[i]);
       }
     }
   }
 }

References funct::abs(), PVValHelper::eta, etabins, getEtaIndex(), mps_fire::i, iEvent, amptDefaultParameters_cff::mu, muHigh, muToken, dqmiodumpmetadata::n, HLT_2018_cff::nEta, npvHigh, offsetPlotBaseName, offsetPlots, pdgMap, pfToken, pftypes, MetAnalyzer::pv(), pvToken, alignCSCRings::s, AlCaHLTBitMon_QueryRunRegistry::string, edmLumisInFiles::summary, and th1dPlots.

◆ bookHistograms()

void OffsetAnalyzerDQM::bookHistograms	(	DQMStore::IBooker &	booker,
		edm::Run const &	,
		edm::EventSetup const &
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 142 of file OffsetAnalyzerDQM.cc.

                                                                                                    {
   //std::cout << "OffsetAnalyzerDQM booking offset histograms" << std::endl;
   for (auto& pair : offsetPlots) {
     pair.second.book(booker);
   }
   for (auto& pair : th1dPlots) {
     pair.second.book(booker);
   }
 }

References offsetPlots, and th1dPlots.

◆ getEtaIndex()

int OffsetAnalyzerDQM::getEtaIndex ( float eta )

protected

Definition at line 272 of file OffsetAnalyzerDQM.cc.

                                             {
   int nEta = etabins.size() - 1;
  
   for (int i = 0; i < nEta; i++) {
     if (etabins[i] <= eta && eta < etabins[i + 1])
       return i;
   }
   if (eta == etabins[nEta])
     return nEta - 1;
   else
     return -1;
 }