CTPPSProtonReconstructionPlotter::MultiRPPlots Struct Reference

Public Member Functions
void	fill (const reco::ForwardProton &p)
	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_vtx_y_vs_xi
std::unique_ptr< TH1D >	h_chi_sq
std::unique_ptr< TH1D >	h_chi_sq_norm
std::unique_ptr< TH1D >	h_n_timing_RPs
std::unique_ptr< TH1D >	h_n_tracking_RPs
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_vtx_y
std::unique_ptr< TH1D >	h_xi
std::unique_ptr< TProfile >	p_th_x_vs_xi
std::unique_ptr< TProfile >	p_th_y_vs_xi
std::unique_ptr< TProfile >	p_vtx_y_vs_xi

Detailed Description

Definition at line 106 of file CTPPSProtonReconstructionPlotter.cc.

Constructor & Destructor Documentation

CTPPSProtonReconstructionPlotter::MultiRPPlots::MultiRPPlots ( )

inline

Definition at line 114 of file CTPPSProtonReconstructionPlotter.cc.

References lumiQTWidget::t.

                     :
        h_xi(new TH1D("", ";#xi", 100, 0., 0.3)),
        h_th_x(new TH1D("", ";#theta_{x}   (rad)", 100, -500E-6, +500E-6)),
        h_th_y(new TH1D("", ";#theta_{y}   (rad)", 100, -500E-6, +500E-6)),
        h_vtx_y(new TH1D("", ";vtx_{y}   (cm)", 100, -2., +2.)),
        h_chi_sq(new TH1D("", ";#chi^{2}", 100, 0., 0.)),
        h_chi_sq_norm(new TH1D("", ";#chi^{2}/ndf", 100, 0., 5.)),
        h_n_tracking_RPs(new TH1D("", ";n of tracking RPs", 4, -0.5, +3.5)),
        h_n_timing_RPs(new TH1D("", ";n of timing RPs", 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, -500E-3, +500E-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))
      {
        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.reset(new TH1D("", ";|t|   (GeV^2)", 100, 0., 5.));
        h_t.reset(new TH1D("", ";|t|   (GeV^2)", v_t_bin_edges.size() - 1, v_t_bin_edges.data()));
        h_t_xi_range1.reset(new TH1D("", ";|t|   (GeV^2)", v_t_bin_edges.size() - 1, v_t_bin_edges.data()));
        h_t_xi_range2.reset(new TH1D("", ";|t|   (GeV^2)", v_t_bin_edges.size() - 1, v_t_bin_edges.data()));
        h_t_xi_range3.reset(new TH1D("", ";|t|   (GeV^2)", v_t_bin_edges.size() - 1, v_t_bin_edges.data()));
        h2_t_vs_xi.reset(new TH2D("", ";#xi;|t|   (GeV^2)", 100, 0., 0.3, v_t_bin_edges.size() - 1, v_t_bin_edges.data()));
      }

Member Function Documentation

void CTPPSProtonReconstructionPlotter::MultiRPPlots::fill ( const reco::ForwardProton &  p )

inline

Definition at line 144 of file CTPPSProtonReconstructionPlotter.cc.

References reco::ForwardProton::chi2(), reco::ForwardProton::contributingLocalTracks(), TtSemiLepEvtBuilder_cfi::mt, reco::ForwardProton::ndof(), reco::ForwardProton::normalizedChi2(), CTPPSDetId::sdTimingDiamond, CTPPSDetId::sdTimingFastSilicon, CTPPSDetId::sdTrackingPixel, CTPPSDetId::sdTrackingStrip, reco::ForwardProton::t(), reco::ForwardProton::thetaX(), reco::ForwardProton::thetaY(), reco::ForwardProton::validFit(), reco::ForwardProton::vertex(), and reco::ForwardProton::xi().

      {
        if (!p.validFit())
          return;

        unsigned int n_tracking_RPs=0, n_timing_RPs=0;
        for (const auto &tr : p.contributingLocalTracks())
        {
          CTPPSDetId detId(tr->getRPId());
          if (detId.subdetId() == CTPPSDetId::sdTrackingStrip || detId.subdetId() == CTPPSDetId::sdTrackingPixel)
            n_tracking_RPs++;
          if (detId.subdetId() == CTPPSDetId::sdTimingDiamond || detId.subdetId() == CTPPSDetId::sdTimingFastSilicon)
            n_timing_RPs++;
        }

        const double th_x = p.thetaX();
        const double th_y = p.thetaY();
        const double mt = - p.t();

        h_chi_sq->Fill(p.chi2());
        if (p.ndof() > 0)
          h_chi_sq_norm->Fill(p.normalizedChi2());

        h_n_tracking_RPs->Fill(n_tracking_RPs);
        h_n_timing_RPs->Fill(n_timing_RPs);

        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);

        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());
      }

void CTPPSProtonReconstructionPlotter::MultiRPPlots::write ( ) const

inline

Definition at line 193 of file CTPPSProtonReconstructionPlotter.cc.

References CTPPSProtonReconstructionPlotter::profileToRMSGraph().

      {
        h_chi_sq->Write("h_chi_sq");
        h_chi_sq_norm->Write("h_chi_sq_norm");

        h_n_tracking_RPs->Write("h_n_tracking_RPs");
        h_n_timing_RPs->Write("h_n_timing_RPs");

        h_xi->Write("h_xi");

        h_th_x->Write("h_th_x");
        h2_th_x_vs_xi->Write("h2_th_x_vs_xi");
        p_th_x_vs_xi->Write("p_th_x_vs_xi");
        auto g_th_x_RMS_vs_xi = std::make_unique<TGraphErrors>();
        profileToRMSGraph(p_th_x_vs_xi.get(), g_th_x_RMS_vs_xi.get());
        g_th_x_RMS_vs_xi->Write("g_th_x_RMS_vs_xi");

        h_th_y->Write("h_th_y");
        h2_th_y_vs_xi->Write("h2_th_y_vs_xi");
        p_th_y_vs_xi->Write("p_th_y_vs_xi");
        auto g_th_y_RMS_vs_xi = std::make_unique<TGraphErrors>();
        profileToRMSGraph(p_th_y_vs_xi.get(), g_th_y_RMS_vs_xi.get());
        g_th_y_RMS_vs_xi->Write("g_th_y_RMS_vs_xi");

        h_vtx_y->Write("h_vtx_y");
        h2_vtx_y_vs_xi->Write("h2_vtx_y_vs_xi");
        p_vtx_y_vs_xi->Write("p_vtx_y_vs_xi");
        auto g_vtx_y_RMS_vs_xi = std::make_unique<TGraphErrors>();
        profileToRMSGraph(p_vtx_y_vs_xi.get(), g_vtx_y_RMS_vs_xi.get());
        g_vtx_y_RMS_vs_xi->Write("g_vtx_y_RMS_vs_xi");

        h_t->Scale(1., "width");

        h_t_unif->Write("h_t_unif");
        h_t->Write("h_t");
        h_t_xi_range1->Write("h_t_xi_range1");
        h_t_xi_range2->Write("h_t_xi_range2");
        h_t_xi_range3->Write("h_t_xi_range3");

        h2_t_vs_xi->Write("h2_t_vs_xi");
      }

Member Data Documentation

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

Definition at line 111 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 111 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 111 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 111 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 108 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 108 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 110 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 110 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 108 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 108 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 109 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 109 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 109 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 108 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 108 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 108 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 108 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 112 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 112 of file CTPPSProtonReconstructionPlotter.cc.

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

Definition at line 112 of file CTPPSProtonReconstructionPlotter.cc.