OffsetAnalyzerDQM Class Reference

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
virtual void	analyze (edm::Event const &, edm::EventSetup const &)
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)
void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &, edm::EventSetup const &) final
void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDProducerBase () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) 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
ESProxyIndex const *	esGetTokenIndices (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::vector< ModuleDescription const * > &modules, 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
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override
void	dqmEndRun (const edm::Run &, const edm::EventSetup &) override
int	getEtaIndex (float eta)
Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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)

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::one::EDProducerBase
typedef EDProducerBase	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 16 of file OffsetAnalyzerDQM.cc.

Constructor & Destructor Documentation

OffsetAnalyzerDQM::OffsetAnalyzerDQM

(

const edm::ParameterSet &

iConfig

)

Definition at line 97 of file OffsetAnalyzerDQM.cc.

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

                                                                   {
  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];
}

Member Function Documentation

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

override

Definition at line 154 of file OffsetAnalyzerDQM.cc.

References funct::abs(), PVValHelper::eta, etabins, edm::Event::getByToken(), getEtaIndex(), mps_fire::i, 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.

                                                                            {
  //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]);
      }
    }
  }
}

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 144 of file OffsetAnalyzerDQM.cc.

References offsetPlots, and th1dPlots.

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

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

inlineoverrideprotected

Definition at line 24 of file OffsetAnalyzerDQM.cc.

{}

void OffsetAnalyzerDQM::dqmEndRun ( const edm::Run &  , const edm::EventSetup &    )

inlineoverrideprotected

Definition at line 25 of file OffsetAnalyzerDQM.cc.

References PVValHelper::eta, and getEtaIndex().

{}

int OffsetAnalyzerDQM::getEtaIndex ( float  eta )

protected

Definition at line 274 of file OffsetAnalyzerDQM.cc.

References DEFINE_FWK_MODULE, etabins, mps_fire::i, and HLT_2018_cff::nEta.

Referenced by analyze(), and dqmEndRun().

                                            {
  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;
}

Member Data Documentation

std::vector<double> OffsetAnalyzerDQM::etabins

private

Definition at line 87 of file OffsetAnalyzerDQM.cc.

Referenced by analyze(), getEtaIndex(), and OffsetAnalyzerDQM().

int OffsetAnalyzerDQM::muHigh

private

Definition at line 89 of file OffsetAnalyzerDQM.cc.

Referenced by analyze(), and OffsetAnalyzerDQM().

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

private

Definition at line 93 of file OffsetAnalyzerDQM.cc.

Referenced by analyze(), and OffsetAnalyzerDQM().

int OffsetAnalyzerDQM::npvHigh

private

Definition at line 90 of file OffsetAnalyzerDQM.cc.

Referenced by analyze(), and OffsetAnalyzerDQM().

std::string OffsetAnalyzerDQM::offsetPlotBaseName

private

Definition at line 85 of file OffsetAnalyzerDQM.cc.

Referenced by analyze(), and OffsetAnalyzerDQM().

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

private

Definition at line 82 of file OffsetAnalyzerDQM.cc.

Referenced by analyze(), bookHistograms(), and OffsetAnalyzerDQM().

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

private

Definition at line 83 of file OffsetAnalyzerDQM.cc.

Referenced by analyze(), and OffsetAnalyzerDQM().

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

private

Definition at line 94 of file OffsetAnalyzerDQM.cc.

Referenced by analyze(), and OffsetAnalyzerDQM().

std::vector<std::string> OffsetAnalyzerDQM::pftypes

private

Definition at line 86 of file OffsetAnalyzerDQM.cc.

Referenced by analyze(), and OffsetAnalyzerDQM().

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

private

Definition at line 92 of file OffsetAnalyzerDQM.cc.

Referenced by analyze(), and OffsetAnalyzerDQM().

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

private

Definition at line 81 of file OffsetAnalyzerDQM.cc.

Referenced by analyze(), bookHistograms(), and OffsetAnalyzerDQM().