CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData Struct Reference

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

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

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

inline

Definition at line 202 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

References CTPPSDetId::arm(), LHCOpticalFunctionsSet::exd, LHCOpticalFunctionsSet::getFcnValues(), LHCOpticalFunctionsSet::getXiValues(), 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());
    }

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

Member Data Documentation

double CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::de_x_mean

Definition at line 86 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

Referenced by write().

double CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::de_x_sigma

Definition at line 86 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

Referenced by write().

double CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::de_y_mean

Definition at line 87 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

Referenced by write().

double CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::de_y_sigma

Definition at line 87 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

Referenced by write().

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

Definition at line 137 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

Referenced by ArmData(), and write().

std::unique_ptr<TH2D> CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::h2_de_y_vs_de_x

Definition at line 84 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

Referenced by write().

std::unique_ptr<TH1D> CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::h_de_x

Definition at line 83 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

Referenced by write().

std::unique_ptr<TH1D> CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::h_de_y

Definition at line 83 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

Referenced by write().

unsigned int CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::n_events_with_tracks

Definition at line 81 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

unsigned int CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::n_exp_prot_max

Definition at line 91 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

Referenced by ArmData().

std::vector<double> CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::n_sigmas

Definition at line 89 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

Referenced by ArmData(), and write().

unsigned int CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::rpId_F

Definition at line 77 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

unsigned int CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::rpId_N

Definition at line 77 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

Referenced by UpdateOptics().

std::shared_ptr<const TSpline3> CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::s_x_to_xi_N

Definition at line 79 of file CTPPSProtonReconstructionEfficiencyEstimatorData.cc.

Referenced by UpdateOptics().