CTPPSProtonReconstructionPlotter::MultiRPPlots Struct Reference

Public Member Functions
void	fill (const reco::ForwardProton &p, unsigned int nTracks, bool n1f1)
	MultiRPPlots ()
void	write () const

Public Attributes
std::unique_ptr< TH2D >	h2_t_vs_xi
std::unique_ptr< TH2D >	h2_th_x_vs_xi
std::unique_ptr< TH2D >	h2_th_y_vs_xi
std::unique_ptr< TH2D >	h2_timing_tracks_vs_prot_mult
std::unique_ptr< TH2D >	h2_vtx_y_vs_xi
std::unique_ptr< TH2D >	h2_x_timing_vs_x_tracking_ClCo
std::unique_ptr< TH2D >	h2_y_vs_x_tt0_ClCo
std::unique_ptr< TH2D >	h2_y_vs_x_tt1_ClCo
std::unique_ptr< TH2D >	h2_y_vs_x_ttm_ClCo
std::unique_ptr< TH1D >	h_chi_sq
std::unique_ptr< TH1D >	h_chi_sq_norm
std::unique_ptr< TH1D >	h_de_x_match_timing_vs_tracking
std::unique_ptr< TH1D >	h_de_x_match_timing_vs_tracking_ClCo
std::unique_ptr< TH1D >	h_de_x_rel_timing_vs_tracking
std::unique_ptr< TH1D >	h_de_x_rel_timing_vs_tracking_ClCo
std::unique_ptr< TH1D >	h_de_x_timing_vs_tracking
std::unique_ptr< TH1D >	h_de_x_timing_vs_tracking_ClCo
std::unique_ptr< TH1D >	h_log_chi_sq
std::unique_ptr< TH1D >	h_multiplicity
std::unique_ptr< TH1D >	h_n_contrib_timing_tracks
std::unique_ptr< TH1D >	h_n_contrib_tracking_tracks
std::unique_ptr< TH1D >	h_t
std::unique_ptr< TH1D >	h_t_unif
std::unique_ptr< TH1D >	h_t_xi_range1
std::unique_ptr< TH1D >	h_t_xi_range2
std::unique_ptr< TH1D >	h_t_xi_range3
std::unique_ptr< TH1D >	h_th_x
std::unique_ptr< TH1D >	h_th_y
std::unique_ptr< TH1D >	h_time
std::unique_ptr< TH1D >	h_time_unc
std::unique_ptr< TH1D >	h_vtx_y
std::unique_ptr< TH1D >	h_xi
std::unique_ptr< TH1D >	h_xi_n1f1
std::map< unsigned int, TH1D * >	m_h_xi_nTracks
std::unique_ptr< TProfile >	p_th_x_vs_xi
std::unique_ptr< TProfile >	p_th_y_vs_xi
std::unique_ptr< TProfile >	p_time_unc_vs_x_ClCo
std::unique_ptr< TProfile >	p_time_unc_vs_xi
std::unique_ptr< TProfile >	p_vtx_y_vs_xi

Detailed Description

Definition at line 156 of file CTPPSProtonReconstructionPlotter.cc.

Constructor & Destructor Documentation

MultiRPPlots()

CTPPSProtonReconstructionPlotter::MultiRPPlots::MultiRPPlots ( )

inline

Definition at line 179 of file CTPPSProtonReconstructionPlotter.cc.

        : h_multiplicity(new TH1D("", ";reconstructed protons per event", 11, -0.5, 10.5)),
          h_xi(new TH1D("", ";#xi", 100, 0., 0.3)),
          h_th_x(new TH1D("", ";#theta_{x}   (rad)", 250, -500E-6, +500E-6)),
          h_th_y(new TH1D("", ";#theta_{y}   (rad)", 500, -1000E-6, +1000E-6)),
          h_vtx_y(new TH1D("", ";vtx_{y}   (cm)", 100, -100E-3, +100E-3)),
          h_chi_sq(new TH1D("", ";#chi^{2}", 100, 0., 10.)),
          h_log_chi_sq(new TH1D("", ";log_{10} #chi^{2}", 100, -20., 5.)),
          h_chi_sq_norm(new TH1D("", ";#chi^{2}/ndf", 100, 0., 5.)),
          h_time(new TH1D("", ";time   (ns)", 100, -2., +2.)),
          h_time_unc(new TH1D("", ";time unc   (ns)", 100, -1., +1.)),
          p_time_unc_vs_x_ClCo(new TProfile("", ";x_tracking   (mm);time unc   (ns)", 100, 0., 30.)),
          p_time_unc_vs_xi(new TProfile("", ";xi;time unc   (ns)", 100, 0., 0.3)),
          h_n_contrib_tracking_tracks(new TH1D("", ";n of contrib. tracking tracks per reco proton", 4, -0.5, +3.5)),
          h_n_contrib_timing_tracks(new TH1D("", ";n of contrib. timing tracks per reco proton", 4, -0.5, +3.5)),
          h2_th_x_vs_xi(new TH2D("", ";#xi;#theta_{x}   (rad)", 100, 0., 0.3, 100, -500E-6, +500E-6)),
          h2_th_y_vs_xi(new TH2D("", ";#xi;#theta_{y}   (rad)", 100, 0., 0.3, 100, -500E-6, +500E-6)),
          h2_vtx_y_vs_xi(new TH2D("", ";#xi;vtx_{y}   (cm)", 100, 0., 0.3, 100, -100E-3, +100E-3)),
          p_th_x_vs_xi(new TProfile("", ";#xi;#theta_{x}   (rad)", 100, 0., 0.3)),
          p_th_y_vs_xi(new TProfile("", ";#xi;#theta_{y}   (rad)", 100, 0., 0.3)),
          p_vtx_y_vs_xi(new TProfile("", ";#xi;vtx_{y}   (cm)", 100, 0., 0.3)),
          h2_timing_tracks_vs_prot_mult(
              new TH2D("", ";reco protons per event;timing tracks per event", 11, -0.5, 10.5, 11, -0.5, 10.5)),
          h_xi_n1f1(new TH1D("", ";#xi", 100, 0., 0.3)),
 
          h2_x_timing_vs_x_tracking_ClCo(
              new TH2D("", ";x_tracking   (mm);x_timing   (mm)", 100, 0., 20., 100, 0., 20.)),
          h_de_x_timing_vs_tracking(new TH1D("", ";#Delta x   (mm)", 200, -1., +1.)),
          h_de_x_rel_timing_vs_tracking(new TH1D("", ";#Delta x / #sigma(x)", 200, -20., +20.)),
          h_de_x_match_timing_vs_tracking(new TH1D("", ";match between tracking and timing tracks", 2, -0.5, +1.5)),
          h_de_x_timing_vs_tracking_ClCo(new TH1D("", ";#Delta x   (mm)", 200, -1., +1.)),
          h_de_x_rel_timing_vs_tracking_ClCo(new TH1D("", ";#Delta x / #sigma(x)", 200, -20., +20.)),
          h_de_x_match_timing_vs_tracking_ClCo(
              new TH1D("", ";match between tracking and timing tracks", 2, -0.5, +1.5)),
 
          h2_y_vs_x_tt0_ClCo(new TH2D("", ";x   (mm);y   (mm)", 100, -5., 25., 100, -15., +15.)),
          h2_y_vs_x_tt1_ClCo(new TH2D("", ";x   (mm);y   (mm)", 100, -5., 25., 100, -15., +15.)),
          h2_y_vs_x_ttm_ClCo(new TH2D("", ";x   (mm);y   (mm)", 100, -5., 25., 100, -15., +15.)) {
      std::vector<double> v_t_bin_edges;
      for (double t = 0; t <= 5.;) {
        v_t_bin_edges.push_back(t);
        const double de_t = 0.05 + 0.09 * t + 0.02 * t * t;
        t += de_t;
      }
      h_t_unif = std::make_unique<TH1D>("", ";|t|   (GeV^2)", 100, 0., 5.);
      h_t = std::make_unique<TH1D>("", ";|t|   (GeV^2)", v_t_bin_edges.size() - 1, v_t_bin_edges.data());
      h_t_xi_range1 = std::make_unique<TH1D>("", ";|t|   (GeV^2)", v_t_bin_edges.size() - 1, v_t_bin_edges.data());
      h_t_xi_range2 = std::make_unique<TH1D>("", ";|t|   (GeV^2)", v_t_bin_edges.size() - 1, v_t_bin_edges.data());
      h_t_xi_range3 = std::make_unique<TH1D>("", ";|t|   (GeV^2)", v_t_bin_edges.size() - 1, v_t_bin_edges.data());
      h2_t_vs_xi = std::make_unique<TH2D>(
          "", ";#xi;|t|   (GeV^2)", 100, 0., 0.3, v_t_bin_edges.size() - 1, v_t_bin_edges.data());
 
      for (unsigned int n = 2; n <= 10; ++n)
        m_h_xi_nTracks[n] = new TH1D(*h_xi);
    }
 

References h2_t_vs_xi, h_t, h_t_unif, h_t_xi_range1, h_t_xi_range2, h_t_xi_range3, h_xi, m_h_xi_nTracks, dqmiodumpmetadata::n, and submitPVValidationJobs::t.

Member Function Documentation

fill()

void CTPPSProtonReconstructionPlotter::MultiRPPlots::fill ( const reco::ForwardProton &  p, unsigned int  nTracks, bool  n1f1  )

inline

Definition at line 235 of file CTPPSProtonReconstructionPlotter.cc.

                          : p.contributingLocalTracks()) {
        CTPPSDetId detId(tr->rpId());
        if (detId.subdetId() == CTPPSDetId::sdTrackingStrip || detId.subdetId() == CTPPSDetId::sdTrackingPixel)
          n_contrib_tracking_tracks++;
        if (detId.subdetId() == CTPPSDetId::sdTimingDiamond || detId.subdetId() == CTPPSDetId::sdTimingFastSilicon)
          n_contrib_timing_tracks++;
      }
 
      const double th_x = p.thetaX();
      const double th_y = p.thetaY();
      const double mt = -p.t();
 
      h_chi_sq->Fill(p.chi2());
      h_log_chi_sq->Fill(log10(p.chi2()));
      if (p.ndof() > 0)
        h_chi_sq_norm->Fill(p.normalizedChi2());
 
      h_n_contrib_tracking_tracks->Fill(n_contrib_tracking_tracks);
      h_n_contrib_timing_tracks->Fill(n_contrib_timing_tracks);
 
      h_xi->Fill(p.xi());
 
      h_th_x->Fill(th_x);
      h_th_y->Fill(th_y);
 
      h_vtx_y->Fill(p.vertex().y());
 
      h_t_unif->Fill(mt);
      h_t->Fill(mt);
      if (p.xi() > 0.04 && p.xi() < 0.07)
        h_t_xi_range1->Fill(mt);
      if (p.xi() > 0.07 && p.xi() < 0.10)
        h_t_xi_range2->Fill(mt);
      if (p.xi() > 0.10 && p.xi() < 0.13)
        h_t_xi_range3->Fill(mt);
 
      if (p.timeError() > 0.) {
        h_time->Fill(p.time());
        h_time_unc->Fill(p.timeError());
        //p_time_unc_vs_x_ClCo filled in ClCo code below
        p_time_unc_vs_xi->Fill(p.xi(), p.timeError());
      }
 
      h2_th_x_vs_xi->Fill(p.xi(), th_x);
      h2_th_y_vs_xi->Fill(p.xi(), th_y);
      h2_vtx_y_vs_xi->Fill(p.xi(), p.vertex().y());
      h2_t_vs_xi->Fill(p.xi(), mt);
 
      p_th_x_vs_xi->Fill(p.xi(), th_x);
      p_th_y_vs_xi->Fill(p.xi(), th_y);
      p_vtx_y_vs_xi->Fill(p.xi(), p.vertex().y());
 
      auto it = m_h_xi_nTracks.find(nTracks);
      if (it != m_h_xi_nTracks.end())
        it->second->Fill(p.xi());
 
      if (n1f1)
        h_xi_n1f1->Fill(p.xi());
    }
 

References CTPPSDetId::sdTimingDiamond, CTPPSDetId::sdTimingFastSilicon, CTPPSDetId::sdTrackingPixel, and CTPPSDetId::sdTrackingStrip.

write()

void CTPPSProtonReconstructionPlotter::MultiRPPlots::write ( ) const

inline

Definition at line 300 of file CTPPSProtonReconstructionPlotter.cc.

                        : m_h_xi_nTracks) {
        char buf[100];
        sprintf(buf, "h_xi_nTracks_%u", p.first);
        p.second->Write(buf);
      }
 
      gDirectory = d_top;
 
      h_xi_n1f1->Write("h_xi_n1f1");
 
      h2_x_timing_vs_x_tracking_ClCo->Write("h2_x_timing_vs_x_tracking_ClCo");
 
      h_de_x_timing_vs_tracking->Write("h_de_x_timing_vs_tracking");
      h_de_x_rel_timing_vs_tracking->Write("h_de_x_rel_timing_vs_tracking");
      h_de_x_match_timing_vs_tracking->Write("h_de_x_match_timing_vs_tracking");
 
      h_de_x_timing_vs_tracking_ClCo->Write("h_de_x_timing_vs_tracking_ClCo");
      h_de_x_rel_timing_vs_tracking_ClCo->Write("h_de_x_rel_timing_vs_tracking_ClCo");
      h_de_x_match_timing_vs_tracking_ClCo->Write("h_de_x_match_timing_vs_tracking_ClCo");
 
      h2_y_vs_x_tt0_ClCo->Write("h2_y_vs_x_tt0_ClCo");
      h2_y_vs_x_tt1_ClCo->Write("h2_y_vs_x_tt1_ClCo");
      h2_y_vs_x_ttm_ClCo->Write("h2_y_vs_x_ttm_ClCo");
    }
  };

Member Data Documentation

h2_t_vs_xi

std::unique_ptr<TH2D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h2_t_vs_xi

Definition at line 163 of file CTPPSProtonReconstructionPlotter.cc.

Referenced by MultiRPPlots().

h2_th_x_vs_xi

std::unique_ptr<TH2D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h2_th_x_vs_xi

Definition at line 163 of file CTPPSProtonReconstructionPlotter.cc.

h2_th_y_vs_xi

std::unique_ptr<TH2D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h2_th_y_vs_xi

Definition at line 163 of file CTPPSProtonReconstructionPlotter.cc.

h2_timing_tracks_vs_prot_mult

std::unique_ptr<TH2D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h2_timing_tracks_vs_prot_mult

Definition at line 166 of file CTPPSProtonReconstructionPlotter.cc.

h2_vtx_y_vs_xi

std::unique_ptr<TH2D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h2_vtx_y_vs_xi

Definition at line 163 of file CTPPSProtonReconstructionPlotter.cc.

h2_x_timing_vs_x_tracking_ClCo

std::unique_ptr<TH2D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h2_x_timing_vs_x_tracking_ClCo

Definition at line 171 of file CTPPSProtonReconstructionPlotter.cc.

h2_y_vs_x_tt0_ClCo

std::unique_ptr<TH2D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h2_y_vs_x_tt0_ClCo

Definition at line 177 of file CTPPSProtonReconstructionPlotter.cc.

h2_y_vs_x_tt1_ClCo

std::unique_ptr<TH2D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h2_y_vs_x_tt1_ClCo

Definition at line 177 of file CTPPSProtonReconstructionPlotter.cc.

h2_y_vs_x_ttm_ClCo

std::unique_ptr<TH2D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h2_y_vs_x_ttm_ClCo

Definition at line 177 of file CTPPSProtonReconstructionPlotter.cc.

h_chi_sq

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_chi_sq

Definition at line 158 of file CTPPSProtonReconstructionPlotter.cc.

h_chi_sq_norm

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_chi_sq_norm

Definition at line 158 of file CTPPSProtonReconstructionPlotter.cc.

h_de_x_match_timing_vs_tracking

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_de_x_match_timing_vs_tracking

Definition at line 173 of file CTPPSProtonReconstructionPlotter.cc.

h_de_x_match_timing_vs_tracking_ClCo

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_de_x_match_timing_vs_tracking_ClCo

Definition at line 174 of file CTPPSProtonReconstructionPlotter.cc.

h_de_x_rel_timing_vs_tracking

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_de_x_rel_timing_vs_tracking

Definition at line 173 of file CTPPSProtonReconstructionPlotter.cc.

h_de_x_rel_timing_vs_tracking_ClCo

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_de_x_rel_timing_vs_tracking_ClCo

Definition at line 174 of file CTPPSProtonReconstructionPlotter.cc.

h_de_x_timing_vs_tracking

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_de_x_timing_vs_tracking

Definition at line 173 of file CTPPSProtonReconstructionPlotter.cc.

h_de_x_timing_vs_tracking_ClCo

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_de_x_timing_vs_tracking_ClCo

Definition at line 174 of file CTPPSProtonReconstructionPlotter.cc.

h_log_chi_sq

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_log_chi_sq

Definition at line 158 of file CTPPSProtonReconstructionPlotter.cc.

h_multiplicity

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_multiplicity

Definition at line 157 of file CTPPSProtonReconstructionPlotter.cc.

h_n_contrib_timing_tracks

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_n_contrib_timing_tracks

Definition at line 162 of file CTPPSProtonReconstructionPlotter.cc.

h_n_contrib_tracking_tracks

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_n_contrib_tracking_tracks

Definition at line 162 of file CTPPSProtonReconstructionPlotter.cc.

h_t

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_t

Definition at line 158 of file CTPPSProtonReconstructionPlotter.cc.

Referenced by MultiRPPlots().

h_t_unif

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_t_unif

Definition at line 158 of file CTPPSProtonReconstructionPlotter.cc.

Referenced by MultiRPPlots().

h_t_xi_range1

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_t_xi_range1

Definition at line 159 of file CTPPSProtonReconstructionPlotter.cc.

Referenced by MultiRPPlots().

h_t_xi_range2

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_t_xi_range2

Definition at line 159 of file CTPPSProtonReconstructionPlotter.cc.

Referenced by MultiRPPlots().

h_t_xi_range3

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_t_xi_range3

Definition at line 159 of file CTPPSProtonReconstructionPlotter.cc.

Referenced by MultiRPPlots().

h_th_x

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_th_x

Definition at line 158 of file CTPPSProtonReconstructionPlotter.cc.

h_th_y

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_th_y

Definition at line 158 of file CTPPSProtonReconstructionPlotter.cc.

h_time

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_time

Definition at line 160 of file CTPPSProtonReconstructionPlotter.cc.

h_time_unc

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_time_unc

Definition at line 160 of file CTPPSProtonReconstructionPlotter.cc.

h_vtx_y

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_vtx_y

Definition at line 158 of file CTPPSProtonReconstructionPlotter.cc.

h_xi

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_xi

Definition at line 158 of file CTPPSProtonReconstructionPlotter.cc.

Referenced by MultiRPPlots().

h_xi_n1f1

std::unique_ptr<TH1D> CTPPSProtonReconstructionPlotter::MultiRPPlots::h_xi_n1f1

Definition at line 169 of file CTPPSProtonReconstructionPlotter.cc.

m_h_xi_nTracks

std::map<unsigned int, TH1D *> CTPPSProtonReconstructionPlotter::MultiRPPlots::m_h_xi_nTracks

Definition at line 168 of file CTPPSProtonReconstructionPlotter.cc.

Referenced by MultiRPPlots().

p_th_x_vs_xi

std::unique_ptr<TProfile> CTPPSProtonReconstructionPlotter::MultiRPPlots::p_th_x_vs_xi

Definition at line 164 of file CTPPSProtonReconstructionPlotter.cc.

p_th_y_vs_xi

std::unique_ptr<TProfile> CTPPSProtonReconstructionPlotter::MultiRPPlots::p_th_y_vs_xi

Definition at line 164 of file CTPPSProtonReconstructionPlotter.cc.

p_time_unc_vs_x_ClCo

std::unique_ptr<TProfile> CTPPSProtonReconstructionPlotter::MultiRPPlots::p_time_unc_vs_x_ClCo

Definition at line 161 of file CTPPSProtonReconstructionPlotter.cc.

p_time_unc_vs_xi

std::unique_ptr<TProfile> CTPPSProtonReconstructionPlotter::MultiRPPlots::p_time_unc_vs_xi

Definition at line 161 of file CTPPSProtonReconstructionPlotter.cc.

p_vtx_y_vs_xi

std::unique_ptr<TProfile> CTPPSProtonReconstructionPlotter::MultiRPPlots::p_vtx_y_vs_xi

Definition at line 164 of file CTPPSProtonReconstructionPlotter.cc.