CMS 3D CMS Logo

List of all members | Public Member Functions | Private Member Functions | Private Attributes
hitfit::EtaDepResolution Class Reference

Hold on to $\eta$-dependent resolution. This class acts as a function object and returns Vector_Resolution as a function of $\eta$. It does not assume symmetry between $+\eta$ and $-\eta$. The interpretation of $\eta$ as physics $\eta$ or detector $\eta$ is left to users's implementation. More...

#include <EtaDepResolution.h>

Public Member Functions

const bool CheckEta (double eta) const
 Check is an input $\eta$ value is within the valid $\eta$-range of this instance. More...
 
bool CheckNoOverlap (const std::vector< EtaDepResElement > &v)
 Check for non-overlapping $\eta$-range between $\eta$-dependent resolution elements in a list. More...
 
 EtaDepResolution ()
 Default constructor, instantiate an EtaDepResolution object with empty list of $\eta$-dependent resolution elements. More...
 
 EtaDepResolution (const std::string &default_file)
 Instantiate an EtaDepResolution object with a filename which contains information about $\eta$-dependent resolution. More...
 
const double EtaMax () const
 Return the upper limit of the valid $\eta$-range. More...
 
const double EtaMin () const
 Return the lower limit of the valid $\eta$-range. More...
 
EtaDepResElement GetEtaDepResElement (double &eta) const
 Return the corresponding $\eta$-dependent resolution element (of type EtaDepResElement)for a value of $\eta$. More...
 
const std::vector< EtaDepResElementGetEtaDepResElement () const
 Access the internal list of $\eta$-dependent resolution elements. More...
 
Vector_Resolution GetResolution (double &eta) const
 Return the corresponding resolution for a value of $\eta$. More...
 
Vector_Resolution operator() (double &eta)
 Allow users to call this instance as a function to access the corresponding resolution for an input value of $\eta$. More...
 
std::vector< EtaDepResElement >::size_type Read (const std::string &default_file)
 Read the $\eta$-dependent resolution information from an ASCII text file. More...
 
std::vector< EtaDepResElement >::size_type Read (const Defaults_Text &defs)
 Read the $\eta$-dependent resolution information from a Defaults_Text object. More...
 
void sort ()
 
 ~EtaDepResolution ()
 

Private Member Functions

std::vector< EtaDepResElement >::const_iterator FindResolution (double &eta) const
 Internal method to return the corresponding $\eta$-dependent resolution element for a given $\eta$ value. More...
 

Private Attributes

std::vector< EtaDepResElement_EtaDepResElement
 

Detailed Description

Hold on to $\eta$-dependent resolution. This class acts as a function object and returns Vector_Resolution as a function of $\eta$. It does not assume symmetry between $+\eta$ and $-\eta$. The interpretation of $\eta$ as physics $\eta$ or detector $\eta$ is left to users's implementation.

Usage:
Users should write the $\eta$-dependent resolution in a plain ASCII text file with the following format:

etadep_etaminX = etamin

etadep_etamaxX = etamax

etadep_vecresX = vecres

where:

X is an integer starting from 0 (0, 1, 2, 3, ...). Users can write as many resolutions as long as each of them has different integer, and complete information regarding $\eta$ range an resolution.

etamin is a floating-point number, the lower limit of the valid $\eta$ range.

etamax is a floating-point number, the upper limit of the valid $\eta$ range.

vecres is a string-encoded Vector_Resolution, see the documentation for Vector_Resolution class for details.

The constructor will read the ASCII text file and read all resolutions in the file. Then it will instantiate the EtaDepResolution and sort the internal list of EtaDepResElement.

Author
Haryo Sumowidagdo Suhar.nosp@m.yo.S.nosp@m.umowi.nosp@m.dagd.nosp@m.o@cer.nosp@m.n.ch
Creation date:
June 2009.

Definition at line 86 of file EtaDepResolution.h.

Constructor & Destructor Documentation

◆ EtaDepResolution() [1/2]

hitfit::EtaDepResolution::EtaDepResolution ( )

Default constructor, instantiate an EtaDepResolution object with empty list of $\eta$-dependent resolution elements.

Definition at line 43 of file EtaDepResolution.cc.

43 {}

◆ EtaDepResolution() [2/2]

hitfit::EtaDepResolution::EtaDepResolution ( const std::string &  default_file)

Instantiate an EtaDepResolution object with a filename which contains information about $\eta$-dependent resolution.

Parameters
default_fileThe input ASCII text file.

Definition at line 45 of file EtaDepResolution.cc.

References Read().

45 { Read(default_file); }
std::vector< EtaDepResElement >::size_type Read(const std::string &default_file)
Read the -dependent resolution information from an ASCII text file.

◆ ~EtaDepResolution()

hitfit::EtaDepResolution::~EtaDepResolution ( )

Destructor.

Definition at line 47 of file EtaDepResolution.cc.

47 {}

Member Function Documentation

◆ CheckEta()

const bool hitfit::EtaDepResolution::CheckEta ( double  eta) const

Check is an input $\eta$ value is within the valid $\eta$-range of this instance.

Parameters
etaThe $\eta$ value to be check.
Return:
true if eta is within the valid $\eta$-range.
false if is not within the valid $\eta$-range.

Definition at line 112 of file EtaDepResolution.cc.

References _EtaDepResElement, PVValHelper::eta, and FindResolution().

112 { return FindResolution(eta) != _EtaDepResElement.end(); }
std::vector< EtaDepResElement > _EtaDepResElement
std::vector< EtaDepResElement >::const_iterator FindResolution(double &eta) const
Internal method to return the corresponding -dependent resolution element for a given value...

◆ CheckNoOverlap()

bool hitfit::EtaDepResolution::CheckNoOverlap ( const std::vector< EtaDepResElement > &  v)

Check for non-overlapping $\eta$-range between $\eta$-dependent resolution elements in a list.

Parameters
vThe list of $\eta$-dependent resolution elements to check.
Return:
true if there is no overlap.
false if there is overlap.

Definition at line 87 of file EtaDepResolution.cc.

References mps_fire::i, dqmiolumiharvest::j, and findQualityFiles::v.

Referenced by Read().

87  {
88  for (std::vector<EtaDepResElement>::size_type i = 0; i != v.size(); i++) {
89  for (std::vector<EtaDepResElement>::size_type j = i + 1; j != v.size(); j++) {
90  if (v[i].IsOverlap(v[j])) {
91  return false;
92  }
93  }
94  }
95  return true;
96  }
uint16_t size_type

◆ EtaMax()

const double hitfit::EtaDepResolution::EtaMax ( ) const

Return the upper limit of the valid $\eta$-range.

Definition at line 105 of file EtaDepResolution.cc.

References _EtaDepResElement.

105  {
106  if (!(_EtaDepResElement.empty())) {
107  return _EtaDepResElement.back().EtaMax();
108  }
109  return -999.; // a ridiculously negative large number
110  }
std::vector< EtaDepResElement > _EtaDepResElement

◆ EtaMin()

const double hitfit::EtaDepResolution::EtaMin ( ) const

Return the lower limit of the valid $\eta$-range.

Definition at line 98 of file EtaDepResolution.cc.

References _EtaDepResElement.

98  {
99  if (!(_EtaDepResElement.empty())) {
100  return _EtaDepResElement.front().EtaMin();
101  }
102  return 999.; // a ridiculously positive large number
103  }
std::vector< EtaDepResElement > _EtaDepResElement

◆ FindResolution()

std::vector< EtaDepResElement >::const_iterator hitfit::EtaDepResolution::FindResolution ( double &  eta) const
private

Internal method to return the corresponding $\eta$-dependent resolution element for a given $\eta$ value.

Parameters
etaThe value of $\eta$ whose corresponding resolution is to be found.
Return:
An iterator (pointer) to the corresponding $\eta$-dependent resolution element.

Definition at line 31 of file EtaDepResolution.cc.

References _EtaDepResElement, and PVValHelper::eta.

Referenced by CheckEta(), GetEtaDepResElement(), and GetResolution().

31  {
32  for (std::vector<EtaDepResElement>::const_iterator res = _EtaDepResElement.begin(); res != _EtaDepResElement.end();
33  ++res) {
34  if (res->IsInInterval(eta) || res->IsOnEdge(eta)) {
35  return res;
36  }
37  }
38  return _EtaDepResElement.end();
39  }
Definition: Electron.h:6
std::vector< EtaDepResElement > _EtaDepResElement

◆ GetEtaDepResElement() [1/2]

EtaDepResElement hitfit::EtaDepResolution::GetEtaDepResElement ( double &  eta) const

Return the corresponding $\eta$-dependent resolution element (of type EtaDepResElement)for a value of $\eta$.

Parameters
etaThe $\eta$ value for which the corresponding element
Return:
The corresponding element if the input $\eta$ value is valid. Throw a runtime error if the $\eta$ value is invalid.

Definition at line 128 of file EtaDepResolution.cc.

References PVValHelper::eta, and FindResolution().

128 { return *(FindResolution(eta)); }
std::vector< EtaDepResElement >::const_iterator FindResolution(double &eta) const
Internal method to return the corresponding -dependent resolution element for a given value...

◆ GetEtaDepResElement() [2/2]

const std::vector< EtaDepResElement > hitfit::EtaDepResolution::GetEtaDepResElement ( ) const

Access the internal list of $\eta$-dependent resolution elements.

Return:
The list of $\eta$-dependent resolution elements.

Definition at line 130 of file EtaDepResolution.cc.

References _EtaDepResElement.

130 { return _EtaDepResElement; }
std::vector< EtaDepResElement > _EtaDepResElement

◆ GetResolution()

Vector_Resolution hitfit::EtaDepResolution::GetResolution ( double &  eta) const

Return the corresponding resolution for a value of $\eta$.

Parameters
etaThe $\eta$ value for which the corresponding resolution is desired.
Return:
The corresponding resolution if the input $\eta$ value is valid. Throw a runtime error if the $\eta$ value is invalid.

Definition at line 116 of file EtaDepResolution.cc.

References _EtaDepResElement, PVValHelper::eta, and FindResolution().

Referenced by operator()().

116  {
117  std::vector<EtaDepResElement>::const_iterator etaDepResEleVecIter = FindResolution(eta);
118  if (etaDepResEleVecIter != _EtaDepResElement.end()) {
119  return etaDepResEleVecIter->GetResolution();
120  }
121 
122  std::stringstream message;
123  message << "Error, the given eta value : " << eta << " is not inside the valid eta range!";
124 
125  throw std::runtime_error(message.str());
126  }
std::vector< EtaDepResElement > _EtaDepResElement
std::vector< EtaDepResElement >::const_iterator FindResolution(double &eta) const
Internal method to return the corresponding -dependent resolution element for a given value...

◆ operator()()

Vector_Resolution hitfit::EtaDepResolution::operator() ( double &  eta)

Allow users to call this instance as a function to access the corresponding resolution for an input value of $\eta$.

Parameters
etaThe $\eta$ value for which the corresponding resolution is desired.
Return:
The corresponding resolution if the input $\eta$ value is valid. Throw a runtime error if the $\eta$ value is invalid.

Definition at line 114 of file EtaDepResolution.cc.

References PVValHelper::eta, and GetResolution().

114 { return GetResolution(eta); }
Vector_Resolution GetResolution(double &eta) const
Return the corresponding resolution for a value of .

◆ Read() [1/2]

std::vector< EtaDepResElement >::size_type hitfit::EtaDepResolution::Read ( const std::string &  default_file)

Read the $\eta$-dependent resolution information from an ASCII text file.

Parameters
default_fileThe ASCII text file to read.
Return:
The number of $\eta$-dependent resolution element read from the file.

Definition at line 49 of file EtaDepResolution.cc.

References _EtaDepResElement.

Referenced by EtaDepResolution().

49  {
50  const Defaults_Text defs(default_file);
51  Read(defs);
52  return _EtaDepResElement.size();
53  }
std::vector< EtaDepResElement >::size_type Read(const std::string &default_file)
Read the -dependent resolution information from an ASCII text file.
std::vector< EtaDepResElement > _EtaDepResElement

◆ Read() [2/2]

std::vector< EtaDepResElement >::size_type hitfit::EtaDepResolution::Read ( const Defaults_Text defs)

Read the $\eta$-dependent resolution information from a Defaults_Text object.

Parameters
defsThe Defaults_Text object to read.
Return:
The number of $\eta$-dependent resolution element read from the file.

Definition at line 55 of file EtaDepResolution.cc.

References _EtaDepResElement, CheckNoOverlap(), TrackingDataMCValidation_Standalone_cff::etamax, TrackingDataMCValidation_Standalone_cff::etamin, hitfit::Defaults_Text::exists(), hitfit::Defaults_Text::get_float(), hitfit::Defaults_Text::get_string(), mps_fire::i, and sort().

55  {
56  _EtaDepResElement.clear();
57 
59  std::ostringstream os_etamin;
60  std::ostringstream os_etamax;
61  std::ostringstream os_res;
62 
63  os_etamin << "etadep_etamin" << i;
64  os_etamax << "etadep_etamax" << i;
65  os_res << "etadep_vecres" << i;
66 
67  if (defs.exists(os_etamin.str()) && defs.exists(os_etamax.str()) && defs.exists(os_res.str())) {
68  double etamin = defs.get_float(os_etamin.str());
69  double etamax = defs.get_float(os_etamax.str());
70  Vector_Resolution res(defs.get_string(os_res.str()));
71  _EtaDepResElement.push_back(EtaDepResElement(etamin, etamax, res));
72 
73  } else {
74  break;
75  }
76  }
77 
79  sort();
80  } else {
81  _EtaDepResElement.clear();
82  }
83 
84  return _EtaDepResElement.size();
85  }
bool CheckNoOverlap(const std::vector< EtaDepResElement > &v)
Check for non-overlapping -range between -dependent resolution elements in a list.
uint16_t size_type
Definition: Electron.h:6
std::vector< EtaDepResElement > _EtaDepResElement

◆ sort()

void hitfit::EtaDepResolution::sort ( )

Sort the internal list of $\eta$-dependent resolution elements.

Definition at line 41 of file EtaDepResolution.cc.

References _EtaDepResElement.

Referenced by Read().

41 { std::stable_sort(_EtaDepResElement.begin(), _EtaDepResElement.end()); }
std::vector< EtaDepResElement > _EtaDepResElement

Member Data Documentation

◆ _EtaDepResElement

std::vector<EtaDepResElement> hitfit::EtaDepResolution::_EtaDepResElement
private

List of $\eta$-dependent resolution elements.

Definition at line 91 of file EtaDepResolution.h.

Referenced by CheckEta(), EtaMax(), EtaMin(), FindResolution(), GetEtaDepResElement(), GetResolution(), Read(), and sort().