TotemRPDQMHarvester Class Reference

Inheritance diagram for TotemRPDQMHarvester:

Public Member Functions
	TotemRPDQMHarvester (const edm::ParameterSet &ps)
	~TotemRPDQMHarvester () override

Protected Member Functions
void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override
void	dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, const edm::LuminosityBlock &, const edm::EventSetup &) override

Private Member Functions
void	MakeHitNumberRatios (unsigned int id, DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter)
void	MakePlaneEfficiencyHistograms (unsigned int id, DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, bool &rpPlotInitialized)

Detailed Description

Definition at line 20 of file TotemRPDQMHarvester.cc.

Constructor & Destructor Documentation

TotemRPDQMHarvester::TotemRPDQMHarvester

(

const edm::ParameterSet &

ps

)

Definition at line 50 of file TotemRPDQMHarvester.cc.

{}

TotemRPDQMHarvester::~TotemRPDQMHarvester ( )

override

Definition at line 54 of file TotemRPDQMHarvester.cc.

{}

Member Function Documentation

void TotemRPDQMHarvester::dqmEndJob ( DQMStore::IBooker &  , DQMStore::IGetter &    )

inlineoverrideprotected

Definition at line 26 of file TotemRPDQMHarvester.cc.

{}

void TotemRPDQMHarvester::dqmEndLuminosityBlock ( DQMStore::IBooker &  ibooker, DQMStore::IGetter &  igetter, const edm::LuminosityBlock &  , const edm::EventSetup &    )

overrideprotected

Definition at line 181 of file TotemRPDQMHarvester.cc.

                                                                       {
  // loop over arms
  for (unsigned int arm = 0; arm < 2; arm++) {
    // loop over stations
    for (unsigned int st = 0; st < 3; st += 2) {
      // loop over RPs
      for (unsigned int rp = 0; rp < 6; ++rp) {
        if (st == 2) {
          // unit 220-nr is not equipped
          if (rp <= 2)
            continue;

          // RP 220-fr-hr contains pixels
          if (rp == 3)
            continue;
        }

        TotemRPDetId rpId(arm, st, rp);

        MakeHitNumberRatios(rpId, ibooker, igetter);

        bool rpPlotInitialized = false;

        // loop over planes
        for (unsigned int pl = 0; pl < 10; ++pl) {
          TotemRPDetId plId(arm, st, rp, pl);

          MakePlaneEfficiencyHistograms(plId, ibooker, igetter, rpPlotInitialized);
        }
      }
    }
  }
}

void TotemRPDQMHarvester::MakeHitNumberRatios ( unsigned int  id, DQMStore::IBooker &  ibooker, DQMStore::IGetter &  igetter  )

private

Definition at line 58 of file TotemRPDQMHarvester.cc.

References dqm::implementation::IBooker::book1D(), c, dqm::implementation::IGetter::get(), dqm::impl::MonitorElement::getBinContent(), dqm::impl::MonitorElement::getNbinsX(), dqm::impl::MonitorElement::getNbinsY(), dqm::legacy::MonitorElement::getTH1F(), dqm::legacy::MonitorElement::getTH2F(), CTPPSDetId::nFull, CTPPSDetId::nPath, fed_dqm_sourceclient-live_cfg::path, alignCSCRings::r, CTPPSDetId::rpName(), alignCSCRings::s, dqm::impl::MonitorElement::setBinContent(), dqm::implementation::NavigatorBase::setCurrentFolder(), and runGCPTkAlMap::title.

                                                                                                                 {
  // get source histogram
  string path;
  TotemRPDetId(id).rpName(path, TotemRPDetId::nPath);

  MonitorElement *activity = igetter.get(path + "/activity in planes (2D)");

  if (!activity)
    return;

  // book new histogram, if not yet done
  const string hit_ratio_name = "hit ratio in hot spot";
  MonitorElement *hit_ratio = igetter.get(path + "/" + hit_ratio_name);

  if (hit_ratio == nullptr) {
    ibooker.setCurrentFolder(path);
    string title;
    TotemRPDetId(id).rpName(title, TotemRPDetId::nFull);
    hit_ratio = ibooker.book1D(hit_ratio_name, title + ";plane;N_hits(320<strip<440) / N_hits(all)", 10, -0.5, 9.5);
  } else {
    hit_ratio->getTH1F()->Reset();
  }

  // calculate ratios
  TAxis *y_axis = activity->getTH2F()->GetYaxis();
  for (int bix = 1; bix <= activity->getNbinsX(); ++bix) {
    double S_full = 0., S_sel = 0.;
    for (int biy = 1; biy <= activity->getNbinsY(); ++biy) {
      double c = activity->getBinContent(bix, biy);
      double s = y_axis->GetBinCenter(biy);

      S_full += c;

      if (s > 320. && s < 440.)
        S_sel += c;
    }

    double r = (S_full > 0.) ? S_sel / S_full : 0.;

    hit_ratio->setBinContent(bix, r);
  }
}

void TotemRPDQMHarvester::MakePlaneEfficiencyHistograms ( unsigned int  id, DQMStore::IBooker &  ibooker, DQMStore::IGetter &  igetter, bool &  rpPlotInitialized  )

private

Definition at line 103 of file TotemRPDQMHarvester.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), postValidation_cfi::efficiency, dqm::implementation::IGetter::get(), dqm::impl::MonitorElement::getBinContent(), dqm::impl::MonitorElement::getNbinsX(), dqm::legacy::MonitorElement::getTH1(), dqm::legacy::MonitorElement::getTH1F(), dqm::legacy::MonitorElement::getTH2F(), CTPPSDetId::nFull, CTPPSDetId::nPath, pileupDistInMC::num, AlCaHLTBitMon_ParallelJobs::p, fed_dqm_sourceclient-live_cfg::path, TotemRPDetId::plane(), TotemRPDetId::planeName(), CTPPSDetId::rpId(), CTPPSDetId::rpName(), dqm::impl::MonitorElement::setBinContent(), dqm::impl::MonitorElement::setBinError(), dqm::implementation::NavigatorBase::setCurrentFolder(), mathSSE::sqrt(), and runGCPTkAlMap::title.

                                                                                 {
  TotemRPDetId detId(id);

  // get source histograms
  string path;
  detId.planeName(path, TotemRPDetId::nPath);

  MonitorElement *efficiency_num = igetter.get(path + "/efficiency num");
  MonitorElement *efficiency_den = igetter.get(path + "/efficiency den");

  if (!efficiency_num || !efficiency_den)
    return;

  // book new plane histogram, if not yet done
  const string efficiency_name = "efficiency";
  MonitorElement *efficiency = igetter.get(path + "/" + efficiency_name);

  if (efficiency == nullptr) {
    string title;
    detId.planeName(title, TotemRPDetId::nFull);
    TAxis *axis = efficiency_den->getTH1()->GetXaxis();
    ibooker.setCurrentFolder(path);
    efficiency = ibooker.book1D(
        efficiency_name, title + ";track position   (mm)", axis->GetNbins(), axis->GetXmin(), axis->GetXmax());
  } else {
    efficiency->getTH1F()->Reset();
  }

  // book new RP histogram, if not yet done
  CTPPSDetId rpId = detId.rpId();
  rpId.rpName(path, TotemRPDetId::nPath);
  const string rp_efficiency_name = "plane efficiency";
  MonitorElement *rp_efficiency = igetter.get(path + "/" + rp_efficiency_name);

  if (rp_efficiency == nullptr) {
    string title;
    rpId.rpName(title, TotemRPDetId::nFull);
    TAxis *axis = efficiency_den->getTH1()->GetXaxis();
    ibooker.setCurrentFolder(path);
    rp_efficiency = ibooker.book2D(rp_efficiency_name,
                                   title + ";plane;track position   (mm)",
                                   10,
                                   -0.5,
                                   9.5,
                                   axis->GetNbins(),
                                   axis->GetXmin(),
                                   axis->GetXmax());
    rpPlotInitialized = true;
  } else {
    if (!rpPlotInitialized)
      rp_efficiency->getTH2F()->Reset();
    rpPlotInitialized = true;
  }

  // calculate and fill efficiencies
  for (signed int bi = 1; bi <= efficiency->getNbinsX(); bi++) {
    double num = efficiency_num->getBinContent(bi);
    double den = efficiency_den->getBinContent(bi);

    if (den > 0) {
      double p = num / den;
      double p_unc = sqrt(p * (1. - p) / den);
      efficiency->setBinContent(bi, p);
      efficiency->setBinError(bi, p_unc);

      int pl_bi = detId.plane() + 1;
      rp_efficiency->setBinContent(pl_bi, bi, p);
    } else {
      efficiency->setBinContent(bi, 0.);
    }
  }
}