LHCOpticalFunctionsSet Class Reference

Set of optical functions corresponding to one scoring plane along LHC. More...

#include <LHCOpticalFunctionsSet.h>

Inheritance diagram for LHCOpticalFunctionsSet:

Public Types
enum	{   evx, eLx, e14, exd,   evpx, eLpx, e24, expd,   e32, evy, eLy, eyd,   e42, evpy, eLpy, eypd,   nFunctions }
	indices for m_fcn_values and m_splines data members More...

Public Member Functions
const std::vector< std::vector< double > > &	getFcnValues () const
double	getScoringPlaneZ () const
	returns the position of the scoring plane (LHC/TOTEM convention) More...
const std::vector< double > &	getXiValues () const
	LHCOpticalFunctionsSet ()=default
	LHCOpticalFunctionsSet (const std::string &fileName, const std::string &directoryName, double z)
	fills m_*_values fields from a ROOT file More...
	~LHCOpticalFunctionsSet ()=default

Protected Attributes
std::vector< std::vector< double > >	m_fcn_values
	length unit cm More...
std::vector< double >	m_xi_values
double	m_z
	position of the scoring plane, in LHC/TOTEM convention, cm More...

Private Member Functions
template<class Archive >
void	serialize (Archive &ar, const unsigned int version)

Friends
class	boost::serialization::access
template<typename CondSerializationT , typename Enabled >
struct	cond::serialization::access

Detailed Description

Set of optical functions corresponding to one scoring plane along LHC.

Definition at line 13 of file LHCOpticalFunctionsSet.h.

Member Enumeration Documentation

anonymous enum

indices for m_fcn_values and m_splines data members

Enumerator
evx
eLx
e14
exd
evpx
eLpx
e24
expd
e32
evy
eLy
eyd
e42
evpy
eLpy
eypd
nFunctions

Definition at line 17 of file LHCOpticalFunctionsSet.h.

{ evx, eLx, e14, exd, evpx, eLpx, e24, expd, e32, evy, eLy, eyd, e42, evpy, eLpy, eypd, nFunctions };

Constructor & Destructor Documentation

LHCOpticalFunctionsSet::LHCOpticalFunctionsSet ( )

default

LHCOpticalFunctionsSet::LHCOpticalFunctionsSet	(	const std::string &	fileName,
		const std::string &	directoryName,
		double	z
	)

fills m_*_values fields from a ROOT file

Definition at line 12 of file LHCOpticalFunctionsSet.cc.

References e14, e24, e32, e42, eLpx, eLpy, eLx, eLy, evpx, evpy, evx, evy, Exception, exd, expd, eyd, eypd, cuy::graphs, m_fcn_values, m_xi_values, nFunctions, pi, AlCaHLTBitMon_QueryRunRegistry::string, GlobalPosition_Frontier_DevDB_cff::tag, and hybridSuperClusters_cfi::xi.

                                                                                                                  :
  m_z(z)
{
  TFile *f_in = TFile::Open(fileName.c_str());
  if (f_in == nullptr)
    throw cms::Exception("LHCOpticalFunctionsSet") << "Cannot open file " << fileName << ".";

  std::vector<TGraph*> graphs(nFunctions);
  for (unsigned int fi = 0; fi < nFunctions; ++fi)
  {
    std::string tag;
    if (fi == evx) tag = "v_x";
    else if (fi == eLx) tag = "L_x";
    else if (fi == e14) tag = "E_14";
    else if (fi == exd) tag = "x_D";
    else if (fi == evpx) tag = "vp_x";
    else if (fi == eLpx) tag = "Lp_x";
    else if (fi == e24) tag = "E_24";
    else if (fi == expd) tag = "xp_D";
    else if (fi == e32) tag = "E_32";
    else if (fi == evy) tag = "v_y";
    else if (fi == eLy) tag = "L_y";
    else if (fi == eyd) tag = "y_D";
    else if (fi == e42) tag = "E_42";
    else if (fi == evpy) tag = "vp_y";
    else if (fi == eLpy) tag = "Lp_y";
    else if (fi == eypd) tag = "yp_D";
    else
      throw cms::Exception("LHCOpticalFunctionsSet") << "Invalid tag for optical functions: \"" << fi << "\"";

    std::string objPath = directoryName + "/g_" + tag + "_vs_xi";
    auto gr_obj = dynamic_cast<TGraph*>( f_in->Get(objPath.c_str()) );
    if (!gr_obj)
      throw cms::Exception("LHCOpticalFunctionsSet") << "Cannot load object " << objPath << " from file " << fileName << ".";

    graphs[fi] = gr_obj;
  }

  const unsigned int num_xi_vals = graphs[0]->GetN();
  m_xi_values.resize(num_xi_vals);

  m_fcn_values.resize(nFunctions);

  for (unsigned int fi = 0; fi < nFunctions; ++fi)
    m_fcn_values[fi].resize(num_xi_vals);

  for (unsigned int pi = 0; pi < num_xi_vals; ++pi)
  {
    const double xi = graphs[0]->GetX()[pi];
    m_xi_values[pi] = xi;

    for (unsigned int fi = 0; fi < m_fcn_values.size(); ++fi)
      m_fcn_values[fi][pi] = graphs[fi]->Eval(xi);
  }

  delete f_in;
}

LHCOpticalFunctionsSet::~LHCOpticalFunctionsSet ( )

default

Member Function Documentation

const std::vector<std::vector<double> >& LHCOpticalFunctionsSet::getFcnValues ( ) const

inline

Definition at line 30 of file LHCOpticalFunctionsSet.h.

References m_fcn_values.

Referenced by ProtonReconstructionAlgorithm::init(), and CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::UpdateOptics().

{ return m_fcn_values; }

double LHCOpticalFunctionsSet::getScoringPlaneZ ( ) const

inline

returns the position of the scoring plane (LHC/TOTEM convention)

Definition at line 27 of file LHCOpticalFunctionsSet.h.

References m_z.

{ return m_z; }

const std::vector<double>& LHCOpticalFunctionsSet::getXiValues ( ) const

inline

Definition at line 29 of file LHCOpticalFunctionsSet.h.

References m_xi_values.

Referenced by ProtonReconstructionAlgorithm::init(), and CTPPSProtonReconstructionEfficiencyEstimatorData::ArmData::UpdateOptics().

{ return m_xi_values; }

template<class Archive >

void LHCOpticalFunctionsSet::serialize ( Archive &  ar, const unsigned int  version  )

private

Friends And Related Function Documentation

friend class boost::serialization::access

friend

Definition at line 39 of file LHCOpticalFunctionsSet.h.

template<typename CondSerializationT , typename Enabled >

friend struct cond::serialization::access

friend

Definition at line 39 of file LHCOpticalFunctionsSet.h.

Member Data Documentation

std::vector<std::vector<double> > LHCOpticalFunctionsSet::m_fcn_values

protected

length unit cm

Definition at line 37 of file LHCOpticalFunctionsSet.h.

Referenced by getFcnValues(), LHCInterpolatedOpticalFunctionsSet::initializeSplines(), LHCOpticalFunctionsSet(), and CTPPSInterpolatedOpticalFunctionsESSource::produce().

std::vector<double> LHCOpticalFunctionsSet::m_xi_values

protected

Definition at line 36 of file LHCOpticalFunctionsSet.h.

Referenced by getXiValues(), LHCInterpolatedOpticalFunctionsSet::initializeSplines(), LHCOpticalFunctionsSet(), and CTPPSInterpolatedOpticalFunctionsESSource::produce().

double LHCOpticalFunctionsSet::m_z

protected

position of the scoring plane, in LHC/TOTEM convention, cm

Definition at line 34 of file LHCOpticalFunctionsSet.h.

Referenced by getScoringPlaneZ(), and CTPPSInterpolatedOpticalFunctionsESSource::produce().