ATLAS Offline Software
RtResolutionFromPoints.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
3 */
4 
6 
8 #include "GaudiKernel/MsgStream.h"
12 
13 using namespace MuonCalib;
14 
15 //*****************************************************************************
16 
17 //::::::::::::::::::::::::
18 //:: METHOD get_min_max ::
19 //::::::::::::::::::::::::
20 void RtResolutionFromPoints::get_min_max(const std::vector<SamplePoint> &sample_points, double &x_min, double &x_max) {
21  for (unsigned int k = 0; k < sample_points.size(); k++) {
22  if (k == 0) {
23  x_min = sample_points[k].x1();
24  x_max = x_min;
25  }
26  if (x_min > sample_points[k].x1()) { x_min = sample_points[k].x1(); }
27  if (x_max < sample_points[k].x1()) { x_max = sample_points[k].x1(); }
28  }
29 
30  return;
31 }
32 
33 //*****************************************************************************
34 
35 //:::::::::::::::::::::::::::::::::::::
36 //:: METHOD getRtResolutionChebyshev ::
37 //:::::::::::::::::::::::::::::::::::::
38 RtResolutionChebyshev RtResolutionFromPoints::getRtResolutionChebyshev(const std::vector<SamplePoint> &sample_points,
39  const unsigned int &order) {
41  // VARIABLES //
43  MsgStream log(Athena::getMessageSvc(), "RtResolutionChebyshev");
44 
45  std::vector<double> res_param(order + 3); // input parameters of RtChebyshev
46  BaseFunctionFitter fitter(order + 1); // Chebyshev fitter
47  ChebyshevPolynomial chebyshev; // Chebyshev polynomial
48  std::vector<SamplePoint> my_points(sample_points); // copy of the sample to add reduced times
49 
51  // GET THE MINIMUM AND MAXIMUM TIMES AND CALCULATE REDUCED TIMES //
53  get_min_max(sample_points, res_param[0], res_param[1]);
54 
55  double mean(0.5 * (res_param[1] + res_param[0]));
56  double length(res_param[1] - res_param[0]);
57  for (unsigned int k = 0; k < my_points.size(); k++) { my_points[k].set_x1(2 * (sample_points[k].x1() - mean) / length); }
58 
60  // PERFORM A CHEBYSHEV FIT TO THE SAMPLE POINTS //
62  fitter.fit_parameters(my_points, 1, sample_points.size(), chebyshev);
63  for (unsigned int k = 0; k < order + 1; k++) { res_param[k + 2] = fitter.coefficients()[k]; }
64 
66  // CREATE AN RtChebyshev OBJECT WITH THE CORRECT PARAMETERS //
68  RtResolutionChebyshev rt_res_chebyshev(res_param);
69 
70  return rt_res_chebyshev;
71 }
72 
73 //*****************************************************************************
74 
75 //::::::::::::::::::::::::::::::::::
76 //:: METHOD getRtResolutionLookUp ::
77 //::::::::::::::::::::::::::::::::::
78 RtResolutionLookUp RtResolutionFromPoints::getRtResolutionLookUp(const std::vector<SamplePoint> &sample_points) {
80  // VARIABLES //
82  RtResolutionChebyshev res(getRtResolutionChebyshev(sample_points, 8)); // auxiliary resolution-t
83  unsigned int nb_points(100); // number of (r, t) points
84  double bin_width((res.tUpper() - res.tLower()) / static_cast<double>(nb_points)); // step size
85  std::vector<double> res_param(nb_points + 2); // r-t parameters
86 
88  // CREATE AN RtRelationLookUp OBJECT WITH THE CORRECT PARAMETERS //
90  res_param[0] = res.tLower();
91  res_param[1] = bin_width;
92  for (unsigned int k = 0; k < nb_points; k++) { res_param[k + 2] = res.resolution(res.tLower() + k * bin_width); }
93  RtResolutionLookUp rt_res_relation_look_up(res_param);
94 
95  return rt_res_relation_look_up;
96 }
LArSamples::FitterData::fitter
const ShapeFitter * fitter
Definition: ShapeFitter.cxx:23
plotBeamSpotCompare.x1
x1
Definition: plotBeamSpotCompare.py:216
mean
void mean(std::vector< double > &bins, std::vector< double > &values, const std::vector< std::string > &files, const std::string &histname, const std::string &tplotname, const std::string &label="")
Definition: dependence.cxx:254
getMessageSvc.h
singleton-like access to IMessageSvc via open function and helper
MuonCalib::ChebyshevPolynomial
Definition: ChebyshevPolynomial.h:17
MuonCalib::RtResolutionFromPoints::getRtResolutionChebyshev
static RtResolutionChebyshev getRtResolutionChebyshev(const std::vector< SamplePoint > &sample_points, const unsigned int &order)
< get an RtResolutionChebyshev resembling the sigma(t) function as described by the sample points in ...
Definition: RtResolutionFromPoints.cxx:38
CheckAppliedSFs.bin_width
bin_width
Definition: CheckAppliedSFs.py:242
MuonCalib::RtResolutionFromPoints::get_min_max
static void get_min_max(const std::vector< SamplePoint > &sample_points, double &x_min, double &x_max)
Definition: RtResolutionFromPoints.cxx:20
Athena::getMessageSvc
IMessageSvc * getMessageSvc(bool quiet=false)
Definition: getMessageSvc.cxx:20
MuonCalib::BaseFunctionFitter
Definition: BaseFunctionFitter.h:39
mc.order
order
Configure Herwig7.
Definition: mc.Herwig7_Dijet.py:12
res
std::pair< std::vector< unsigned int >, bool > res
Definition: JetGroupProductTest.cxx:14
MuonCalib
CscCalcPed - algorithm that finds the Cathode Strip Chamber pedestals from an RDO.
Definition: CscCalcPed.cxx:22
BaseFunctionFitter.h
MuonCalib::RtResolutionFromPoints::getRtResolutionLookUp
static RtResolutionLookUp getRtResolutionLookUp(const std::vector< SamplePoint > &sample_points)
Definition: RtResolutionFromPoints.cxx:78
RtResolutionFromPoints.h
MuonCalib::RtResolutionChebyshev
Definition: RtResolutionChebyshev.h:34
ChebyshevPolynomial.h
python.CaloCondTools.log
log
Definition: CaloCondTools.py:20
PolygonBase.h
length
double length(const pvec &v)
Definition: FPGATrackSimLLPDoubletHoughTransformTool.cxx:26
fitman.k
k
Definition: fitman.py:528
MuonCalib::RtResolutionLookUp
Equidistant look up table for resolution tables with the time as key.
Definition: RtResolutionLookUp.h:23