PPSAlignmentWorker::SectorData Struct Reference

Classes
struct	SlicePlots

Public Member Functions
void	init (DQMStore::IBooker &iBooker, const PPSAlignmentConfig &cfg, const SectorConfig &_scfg, const std::string &folder, bool debug)
unsigned int	process (const CTPPSLocalTrackLiteCollection &tracks, const PPSAlignmentConfig &cfg, bool debug)

Public Attributes
MonitorElement *	h2_cut_h_aft
MonitorElement *	h2_cut_h_bef
MonitorElement *	h2_cut_v_aft
MonitorElement *	h2_cut_v_bef
MonitorElement *	h_q_cut_h_aft
MonitorElement *	h_q_cut_h_bef
MonitorElement *	h_q_cut_v_aft
MonitorElement *	h_q_cut_v_bef
std::map< unsigned int, MonitorElement * >	m_h2_y_vs_x_aft_sel
std::map< unsigned int, MonitorElement * >	m_h2_y_vs_x_bef_sel
std::map< unsigned int, MonitorElement * >	m_h2_y_vs_x_mlt_sel
MonitorElement *	p_x_diffFN_vs_x_N
MonitorElement *	p_y_diffFN_vs_y_F
SectorConfig	scfg
std::map< unsigned int, SlicePlots >	x_slice_plots_F
std::map< unsigned int, SlicePlots >	x_slice_plots_N

Detailed Description

Definition at line 58 of file PPSAlignmentWorker.cc.

Member Function Documentation

init()

void PPSAlignmentWorker::SectorData::init	(	DQMStore::IBooker &	iBooker,
		const PPSAlignmentConfig &	cfg,
		const SectorConfig &	_scfg,
		const std::string &	folder,
		bool	debug
	)

Definition at line 116 of file PPSAlignmentWorker.cc.

                                                      {
  scfg = _scfg;
 
  // binning
  const double bin_size_x = cfg.binning().bin_size_x_;
  const unsigned int n_bins_x = cfg.binning().n_bins_x_;
 
  const double pixel_x_offset = cfg.binning().pixel_x_offset_;
 
  const double x_min_pix = pixel_x_offset, x_max_pix = pixel_x_offset + n_bins_x * bin_size_x;
  const double x_min_str = 0., x_max_str = n_bins_x * bin_size_x;
 
  const unsigned int n_bins_y = cfg.binning().n_bins_y_;
  const double y_min = cfg.binning().y_min_, y_max = cfg.binning().y_max_;
 
  // hit distributions
  iBooker.setCurrentFolder(folder + "/" + scfg.name_ + "/before selection/" + scfg.rp_N_.name_);
  m_h2_y_vs_x_bef_sel[scfg.rp_N_.id_] =
      iBooker.book2DD("h2_y_vs_x", ";x;y", n_bins_x, x_min_str, x_max_str, n_bins_y, y_min, y_max);
  iBooker.setCurrentFolder(folder + "/" + scfg.name_ + "/before selection/" + scfg.rp_F_.name_);
  m_h2_y_vs_x_bef_sel[scfg.rp_F_.id_] =
      iBooker.book2DD("h2_y_vs_x", ";x;y", n_bins_x, x_min_pix, x_max_pix, n_bins_y, y_min, y_max);
 
  iBooker.setCurrentFolder(folder + "/" + scfg.name_ + "/multiplicity selection/" + scfg.rp_N_.name_);
  m_h2_y_vs_x_mlt_sel[scfg.rp_N_.id_] =
      iBooker.book2DD("h2_y_vs_x", ";x;y", n_bins_x, x_min_str, x_max_str, n_bins_y, y_min, y_max);
  iBooker.setCurrentFolder(folder + "/" + scfg.name_ + "/multiplicity selection/" + scfg.rp_F_.name_);
  m_h2_y_vs_x_mlt_sel[scfg.rp_F_.id_] =
      iBooker.book2DD("h2_y_vs_x", ";x;y", n_bins_x, x_min_pix, x_max_pix, n_bins_y, y_min, y_max);
 
  iBooker.setCurrentFolder(folder + "/" + scfg.name_ + "/after selection/" + scfg.rp_N_.name_);
  m_h2_y_vs_x_aft_sel[scfg.rp_N_.id_] =
      iBooker.book2DD("h2_y_vs_x", ";x;y", n_bins_x, x_min_str, x_max_str, n_bins_y, y_min, y_max);
  iBooker.setCurrentFolder(folder + "/" + scfg.name_ + "/after selection/" + scfg.rp_F_.name_);
  m_h2_y_vs_x_aft_sel[scfg.rp_F_.id_] =
      iBooker.book2DD("h2_y_vs_x", ";x;y", n_bins_x, x_min_pix, x_max_pix, n_bins_y, y_min, y_max);
 
  // cut plots
  iBooker.setCurrentFolder(folder + "/" + scfg.name_ + "/cuts/cut_h");
  h_q_cut_h_bef = iBooker.book1DD("h_q_cut_h_bef", ";cq_h", 400, -2., 2.);
  h_q_cut_h_aft = iBooker.book1DD("h_q_cut_h_aft", ";cq_h", 400, -2., 2.);
  h2_cut_h_bef =
      iBooker.book2DD("h2_cut_h_bef", ";x_up;x_dw", n_bins_x, x_min_str, x_max_str, n_bins_x, x_min_pix, x_max_pix);
  h2_cut_h_aft =
      iBooker.book2DD("h2_cut_h_aft", ";x_up;x_dw", n_bins_x, x_min_str, x_max_str, n_bins_x, x_min_pix, x_max_pix);
 
  iBooker.setCurrentFolder(folder + "/" + scfg.name_ + "/cuts/cut_v");
  h_q_cut_v_bef = iBooker.book1DD("h_q_cut_v_bef", ";cq_v", 400, -2., 2.);
  h_q_cut_v_aft = iBooker.book1DD("h_q_cut_v_aft", ";cq_v", 400, -2., 2.);
  h2_cut_v_bef = iBooker.book2DD("h2_cut_v_bef", ";y_up;y_dw", n_bins_y, y_min, y_max, n_bins_y, y_min, y_max);
  h2_cut_v_aft = iBooker.book2DD("h2_cut_v_aft", ";y_up;y_dw", n_bins_y, y_min, y_max, n_bins_y, y_min, y_max);
 
  // near-far plots
  iBooker.setCurrentFolder(folder + "/" + scfg.name_ + "/near_far");
  auto *tmp = new TProfile("", ";x_{N};x_{F} - x_{N}", 100, 0., 20.);
  p_x_diffFN_vs_x_N = iBooker.bookProfile("p_x_diffFN_vs_x_N", tmp);
 
  for (int i = 0; i < scfg.rp_N_.x_slice_n_; i++) {
    const double xMin = scfg.rp_N_.x_slice_min_ + i * scfg.rp_N_.x_slice_w_;
    const double xMax = scfg.rp_N_.x_slice_min_ + (i + 1) * scfg.rp_N_.x_slice_w_;
 
    char buf[100];
    sprintf(buf, "%.1f-%.1f", xMin, xMax);
 
    iBooker.setCurrentFolder(folder + "/" + scfg.name_ + "/near_far/x slices, N/" + buf);
    x_slice_plots_N.insert({i, SlicePlots(iBooker, debug)});
  }
 
  for (int i = 0; i < scfg.rp_F_.x_slice_n_; i++) {
    const double xMin = scfg.rp_F_.x_slice_min_ + i * scfg.rp_F_.x_slice_w_;
    const double xMax = scfg.rp_F_.x_slice_min_ + (i + 1) * scfg.rp_F_.x_slice_w_;
 
    char buf[100];
    sprintf(buf, "%.1f-%.1f", xMin, xMax);
 
    iBooker.setCurrentFolder(folder + "/" + scfg.name_ + "/near_far/x slices, F/" + buf);
    x_slice_plots_F.insert({i, SlicePlots(iBooker, debug)});
  }
}

References dqm::implementation::IBooker::book1DD(), dqm::implementation::IBooker::book2DD(), dqm::implementation::IBooker::bookProfile(), visDQMUpload::buf, looper::cfg, debug, printsummarytable::folder, h2_cut_h_aft, h2_cut_h_bef, h2_cut_v_aft, h2_cut_v_bef, h_q_cut_h_aft, h_q_cut_h_bef, h_q_cut_v_aft, h_q_cut_v_bef, mps_fire::i, RPConfig::id_, m_h2_y_vs_x_aft_sel, m_h2_y_vs_x_bef_sel, m_h2_y_vs_x_mlt_sel, RPConfig::name_, SectorConfig::name_, p_x_diffFN_vs_x_N, SectorConfig::rp_F_, SectorConfig::rp_N_, scfg, dqm::implementation::NavigatorBase::setCurrentFolder(), createJobs::tmp, RPConfig::x_slice_min_, RPConfig::x_slice_n_, x_slice_plots_F, x_slice_plots_N, RPConfig::x_slice_w_, multiplicitycorr_cfi::xMax, and photonAnalyzer_cfi::xMin.

Referenced by PPSAlignmentWorker::bookHistograms().

process()

unsigned int PPSAlignmentWorker::SectorData::process	(	const CTPPSLocalTrackLiteCollection &	tracks,
		const PPSAlignmentConfig &	cfg,
		bool	debug
	)

Definition at line 200 of file PPSAlignmentWorker.cc.

                                                                 {
  CTPPSLocalTrackLiteCollection tracksUp, tracksDw;
 
  for (const auto &tr : tracks) {
    CTPPSDetId rpId(tr.rpId());
    unsigned int rpDecId = rpId.arm() * 100 + rpId.station() * 10 + rpId.rp();
 
    if (rpDecId != scfg.rp_N_.id_ && rpDecId != scfg.rp_F_.id_)
      continue;
 
    double x = tr.x();
    double y = tr.y();
 
    // re-build track object
    CTPPSLocalTrackLite trCorr(tr.rpId(),
                               x,
                               0.,
                               y,
                               0.,
                               tr.tx(),
                               tr.txUnc(),
                               tr.ty(),
                               tr.tyUnc(),
                               tr.chiSquaredOverNDF(),
                               tr.pixelTrackRecoInfo(),
                               tr.numberOfPointsUsedForFit(),
                               tr.time(),
                               tr.timeUnc());
 
    // store corrected track into the right collection
    if (rpDecId == scfg.rp_N_.id_)
      tracksUp.push_back(trCorr);
    if (rpDecId == scfg.rp_F_.id_)
      tracksDw.push_back(trCorr);
  }
 
  // update plots before selection
  for (const auto &tr : tracksUp)
    m_h2_y_vs_x_bef_sel[scfg.rp_N_.id_]->Fill(tr.x(), tr.y());
 
  for (const auto &tr : tracksDw)
    m_h2_y_vs_x_bef_sel[scfg.rp_F_.id_]->Fill(tr.x(), tr.y());
 
  // skip crowded events
  if (tracksUp.size() > cfg.maxRPTracksSize())
    return 0;
 
  if (tracksDw.size() > cfg.maxRPTracksSize())
    return 0;
 
  // update plots with multiplicity selection
  for (const auto &tr : tracksUp)
    m_h2_y_vs_x_mlt_sel[scfg.rp_N_.id_]->Fill(tr.x(), tr.y());
 
  for (const auto &tr : tracksDw)
    m_h2_y_vs_x_mlt_sel[scfg.rp_F_.id_]->Fill(tr.x(), tr.y());
 
  // do the selection
  unsigned int pairsSelected = 0;
 
  for (const auto &trUp : tracksUp) {
    for (const auto &trDw : tracksDw) {
      h2_cut_h_bef->Fill(trUp.x(), trDw.x());
      h2_cut_v_bef->Fill(trUp.y(), trDw.y());
 
      const double cq_h = trDw.x() + scfg.cut_h_a_ * trUp.x() + scfg.cut_h_c_;
      h_q_cut_h_bef->Fill(cq_h);
      const bool cv_h = (std::fabs(cq_h) < cfg.n_si() * scfg.cut_h_si_);
 
      const double cq_v = trDw.y() + scfg.cut_v_a_ * trUp.y() + scfg.cut_v_c_;
      h_q_cut_v_bef->Fill(cq_v);
      const bool cv_v = (std::fabs(cq_v) < cfg.n_si() * scfg.cut_v_si_);
 
      bool cutsPassed = true;
      if (scfg.cut_h_apply_)
        cutsPassed &= cv_h;
      if (scfg.cut_v_apply_)
        cutsPassed &= cv_v;
 
      if (cutsPassed) {
        pairsSelected++;
 
        h_q_cut_h_aft->Fill(cq_h);
        h_q_cut_v_aft->Fill(cq_v);
 
        h2_cut_h_aft->Fill(trUp.x(), trDw.x());
        h2_cut_v_aft->Fill(trUp.y(), trDw.y());
 
        m_h2_y_vs_x_aft_sel[scfg.rp_N_.id_]->Fill(trUp.x(), trUp.y());
        m_h2_y_vs_x_aft_sel[scfg.rp_F_.id_]->Fill(trDw.x(), trDw.y());
 
        p_x_diffFN_vs_x_N->Fill(trUp.x(), trDw.x() - trUp.x());
 
        int idx = (trUp.x() - scfg.rp_N_.x_slice_min_) / scfg.rp_N_.x_slice_w_;
        if (idx >= 0 && idx < scfg.rp_N_.x_slice_n_) {
          x_slice_plots_N[idx].h_y->Fill(trUp.y());
          x_slice_plots_N[idx].p_y_diffFN_vs_y->Fill(trUp.y(), trDw.y() - trUp.y());
          if (debug)
            x_slice_plots_N[idx].h2_y_diffFN_vs_y->Fill(trUp.y(), trDw.y() - trUp.y());
        }
 
        idx = (trDw.x() - scfg.rp_F_.x_slice_min_) / scfg.rp_F_.x_slice_w_;
        if (idx >= 0 && idx < scfg.rp_F_.x_slice_n_) {
          x_slice_plots_F[idx].h_y->Fill(trDw.y());
          x_slice_plots_F[idx].p_y_diffFN_vs_y->Fill(trDw.y(), trDw.y() - trUp.y());
          if (debug)
            x_slice_plots_F[idx].h2_y_diffFN_vs_y->Fill(trDw.y(), trDw.y() - trUp.y());
        }
      }
    }
  }
 
  return pairsSelected;
}

References looper::cfg, SectorConfig::cut_h_a_, SectorConfig::cut_h_apply_, SectorConfig::cut_h_c_, SectorConfig::cut_h_si_, SectorConfig::cut_v_a_, SectorConfig::cut_v_apply_, SectorConfig::cut_v_c_, SectorConfig::cut_v_si_, debug, dqm::impl::MonitorElement::Fill(), h2_cut_h_aft, h2_cut_h_bef, h2_cut_v_aft, h2_cut_v_bef, h_q_cut_h_aft, h_q_cut_h_bef, h_q_cut_v_aft, h_q_cut_v_bef, RPConfig::id_, heavyIonCSV_trainingSettings::idx, if(), m_h2_y_vs_x_aft_sel, m_h2_y_vs_x_bef_sel, m_h2_y_vs_x_mlt_sel, p_x_diffFN_vs_x_N, SectorConfig::rp_F_, SectorConfig::rp_N_, profile_2016_postTS2_cff::rpId, scfg, tracks, x, RPConfig::x_slice_min_, RPConfig::x_slice_n_, x_slice_plots_F, x_slice_plots_N, RPConfig::x_slice_w_, and y.

Referenced by PPSAlignmentWorker::analyze().

Member Data Documentation

h2_cut_h_aft

MonitorElement * PPSAlignmentWorker::SectorData::h2_cut_h_aft

Definition at line 70 of file PPSAlignmentWorker.cc.

Referenced by init(), and process().

h2_cut_h_bef

MonitorElement* PPSAlignmentWorker::SectorData::h2_cut_h_bef

Definition at line 70 of file PPSAlignmentWorker.cc.

Referenced by init(), and process().

h2_cut_v_aft

MonitorElement * PPSAlignmentWorker::SectorData::h2_cut_v_aft

Definition at line 73 of file PPSAlignmentWorker.cc.

Referenced by init(), and process().

h2_cut_v_bef

MonitorElement* PPSAlignmentWorker::SectorData::h2_cut_v_bef

Definition at line 73 of file PPSAlignmentWorker.cc.

Referenced by init(), and process().

h_q_cut_h_aft

MonitorElement * PPSAlignmentWorker::SectorData::h_q_cut_h_aft

Definition at line 69 of file PPSAlignmentWorker.cc.

Referenced by init(), and process().

h_q_cut_h_bef

MonitorElement* PPSAlignmentWorker::SectorData::h_q_cut_h_bef

Definition at line 69 of file PPSAlignmentWorker.cc.

Referenced by init(), and process().

h_q_cut_v_aft

MonitorElement * PPSAlignmentWorker::SectorData::h_q_cut_v_aft

Definition at line 72 of file PPSAlignmentWorker.cc.

Referenced by init(), and process().

h_q_cut_v_bef

MonitorElement* PPSAlignmentWorker::SectorData::h_q_cut_v_bef

Definition at line 72 of file PPSAlignmentWorker.cc.

Referenced by init(), and process().

m_h2_y_vs_x_aft_sel

std::map<unsigned int, MonitorElement *> PPSAlignmentWorker::SectorData::m_h2_y_vs_x_aft_sel

Definition at line 66 of file PPSAlignmentWorker.cc.

Referenced by init(), and process().

m_h2_y_vs_x_bef_sel

std::map<unsigned int, MonitorElement *> PPSAlignmentWorker::SectorData::m_h2_y_vs_x_bef_sel

Definition at line 62 of file PPSAlignmentWorker.cc.

Referenced by init(), and process().

m_h2_y_vs_x_mlt_sel

std::map<unsigned int, MonitorElement *> PPSAlignmentWorker::SectorData::m_h2_y_vs_x_mlt_sel

Definition at line 64 of file PPSAlignmentWorker.cc.

Referenced by init(), and process().

p_x_diffFN_vs_x_N

MonitorElement* PPSAlignmentWorker::SectorData::p_x_diffFN_vs_x_N

Definition at line 76 of file PPSAlignmentWorker.cc.

Referenced by init(), and process().

p_y_diffFN_vs_y_F

MonitorElement* PPSAlignmentWorker::SectorData::p_y_diffFN_vs_y_F

Definition at line 77 of file PPSAlignmentWorker.cc.

scfg

SectorConfig PPSAlignmentWorker::SectorData::scfg

Definition at line 59 of file PPSAlignmentWorker.cc.

Referenced by init(), and process().

x_slice_plots_F

std::map<unsigned int, SlicePlots> PPSAlignmentWorker::SectorData::x_slice_plots_F

Definition at line 88 of file PPSAlignmentWorker.cc.

Referenced by init(), and process().

x_slice_plots_N

std::map<unsigned int, SlicePlots> PPSAlignmentWorker::SectorData::x_slice_plots_N

Definition at line 88 of file PPSAlignmentWorker.cc.

Referenced by init(), and process().