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 70 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

Constructor & Destructor Documentation

◆ ArmData()

CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::ArmData ( )

inline

Definition at line 114 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

                                   {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();
         }
       }
     }
  
     void write() const {

References effPlots, n_exp_prot_max, and n_sigmas.

Member Function Documentation

◆ UpdateOptics()

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

inline

Definition at line 177 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

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

References protons_cff::decRPId, AlCaHLTBitMon_ParallelJobs::p, profile_2016_postTS2_cff::rpId, and rpId_N.

◆ write()

void CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::write ( ) const

inline

Definition at line 132 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

                             : 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;
     }
  
     void UpdateOptics(const LHCInterpolatedOpticalFunctionsSetCollection &ofc) {