Classes
struct	EffPlots

Public Member Functions
	ArmData ()

void	UpdateOptics (const LHCInterpolatedOpticalFunctionsSetCollection &ofc)

void	write () const

Public Attributes
double	de_x_mean

double	de_x_sigma

double	de_y_mean

double	de_y_sigma

std::map< unsigned int, std::map< unsigned int, EffPlots > >	effPlots

std::unique_ptr< TH2D >	h2_de_y_vs_de_x

std::unique_ptr< TH1D >	h_de_x

std::unique_ptr< TH1D >	h_de_y

unsigned int	n_events_with_tracks

unsigned int	n_exp_prot_max

std::vector< double >	n_sigmas

unsigned int	rpId_F

unsigned int	rpId_N

std::shared_ptr< const TSpline3 >	s_x_to_xi_N

Detailed Description

Definition at line 76 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

Constructor & Destructor Documentation

◆ ArmData()

CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::ArmData ( )

inline

Definition at line 139 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

References effPlots, n_exp_prot_max, and n_sigmas.

         : n_events_with_tracks(0),
 
           h_de_x(new TH1D("", ";x_{F} - x_{N}   (mm)", 200, -1., +1.)),
           h_de_y(new TH1D("", ";y_{F} - y_{N}   (mm)", 200, -1., +1.)),
           h2_de_y_vs_de_x(new TH2D("", ";x_{F} - x_{N}   (mm);y_{F} - y_{N}   (mm)", 100, -1., +1., 100, -1., +1.)),
 
           de_x_mean(0.),
           de_x_sigma(0.),
           de_y_mean(0.),
           de_y_sigma(0.) {
       for (unsigned int nep = 0; nep <= n_exp_prot_max; ++nep) {
         for (unsigned int nsi = 0; nsi < n_sigmas.size(); ++nsi) {
           effPlots[nep][nsi] = EffPlots();
         }
       }
     }

Member Function Documentation

◆ UpdateOptics()

void CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::UpdateOptics ( const LHCInterpolatedOpticalFunctionsSetCollection & ofc )

inline

Definition at line 202 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

References CTPPSDetId::arm(), protons_cff::decRPId, LHCOpticalFunctionsSet::exd, haddnano::of, AlCaHLTBitMon_ParallelJobs::p, CTPPSDetId::rp(), rpId_N, s_x_to_xi_N, and CTPPSDetId::station().

                                                                                {
       const LHCInterpolatedOpticalFunctionsSet *of = nullptr;
 
       for (const auto &p : ofc) {
         CTPPSDetId rpId(p.first);
         unsigned int decRPId = rpId.arm() * 100 + rpId.station() * 10 + rpId.rp();
 
         if (decRPId == rpId_N) {
           of = &p.second;
           break;
         }
       }
 
       if (!of) {
         edm::LogError("ArmData::UpdateOptics") << "Cannot find optical functions for RP " << rpId_N;
         return;
       }
 
       std::vector<double> xiValues =
           of->getXiValues();  // local copy made since the TSpline constructor needs non-const parameters
       std::vector<double> xDValues = of->getFcnValues()[LHCOpticalFunctionsSet::exd];
       s_x_to_xi_N = std::make_shared<TSpline3>("", xDValues.data(), xiValues.data(), xiValues.size());
     }

◆ write()

void CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::write ( ) const

inline

Definition at line 157 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

References visDQMUpload::buf, de_x_mean, de_x_sigma, de_y_mean, de_y_sigma, effPlots, g, h2_de_y_vs_de_x, h_de_x, h_de_y, and n_sigmas.

                        {
       h_de_x->Write("h_de_x");
       h_de_y->Write("h_de_y");
       h2_de_y_vs_de_x->Write("h2_de_y_vs_de_x");
 
       char buf[100];
 
       for (const auto &n_si : n_sigmas) {
         sprintf(buf, "g_de_y_vs_de_x_n_si=%.1f", n_si);
         TGraph *g = new TGraph();
         g->SetName(buf);
 
         g->SetPoint(0, de_x_mean - n_si * de_x_sigma, de_y_mean - n_si * de_y_sigma);
         g->SetPoint(1, de_x_mean + n_si * de_x_sigma, de_y_mean - n_si * de_y_sigma);
         g->SetPoint(2, de_x_mean + n_si * de_x_sigma, de_y_mean + n_si * de_y_sigma);
         g->SetPoint(3, de_x_mean - n_si * de_x_sigma, de_y_mean + n_si * de_y_sigma);
         g->SetPoint(4, de_x_mean - n_si * de_x_sigma, de_y_mean - n_si * de_y_sigma);
 
         g->Write();
       }
 
       TDirectory *d_top = gDirectory;
 
       for (const auto &nep_p : effPlots) {
         if (nep_p.first == 0)
           sprintf(buf, "exp prot all");
         else
           sprintf(buf, "exp prot %u", nep_p.first);
 
         TDirectory *d_nep = d_top->mkdir(buf);
 
         for (const auto &nsi_p : nep_p.second) {
           sprintf(buf, "nsi = %.1f", n_sigmas[nsi_p.first]);
 
           TDirectory *d_nsi = d_nep->mkdir(buf);
 
           gDirectory = d_nsi;
 
           nsi_p.second.write();
         }
       }
 
       gDirectory = d_top;
     }