ATLAS Offline Software
Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Private Types | Private Member Functions | Static Private Member Functions | Private Attributes | List of all members
CscBipolarStripFitter Class Reference

#include <CscBipolarStripFitter.h>

Inheritance diagram for CscBipolarStripFitter:
Collaboration diagram for CscBipolarStripFitter:

Public Types

typedef std::vector< float > ChargeList
 

Public Member Functions

 CscBipolarStripFitter (const std::string &, const std::string &, const IInterface *)
 
 ~CscBipolarStripFitter ()=default
 
StatusCode initialize () override
 
Result fit (const ChargeList &charges, double samplingTime, Identifier &stripId) const
 
virtual Result fit (const ChargeList &ChargeList, double samplingTime, bool samplingPhase, Identifier &sid) const
 
virtual Result fit (const Muon::CscStripPrepData &strip) const
 
virtual Result fit (const ChargeList &ChargeList, double samplingTime, bool samplingPhase, Identifier &sid) const
 
virtual Result fit (const Muon::CscStripPrepData &strip) const
 
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & evtStore () const
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
 
virtual StatusCode sysInitialize () override
 Perform system initialization for an algorithm. More...
 
virtual StatusCode sysStart () override
 Handle START transition. More...
 
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles. More...
 
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles. More...
 
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T > &t)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc="none")
 Declare a new Gaudi property. More...
 
void updateVHKA (Gaudi::Details::PropertyBase &)
 
MsgStream & msg () const
 
MsgStream & msg (const MSG::Level lvl) const
 
bool msgLvl (const MSG::Level lvl) const
 

Static Public Member Functions

static const InterfaceID & interfaceID ()
 

Protected Member Functions

void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution More...
 
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce (T &h)
 
void extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
 Add StoreName to extra input/output deps as needed. More...
 

Private Types

typedef ServiceHandle< StoreGateSvc > StoreGateSvc_t
 

Private Member Functions

double FindPow (double z) const
 
void Derivative (double A[][3], double fp[][1], double p0[][1], int imeas, const int *meas) const
 
int TheFitter (double *x, const double ex, const double *initValues, int imeas, int *meas, int ipar, int *par, double *chi2, double *result) const
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Static Private Member Functions

static double FindInitValues (double *x, double *initValues, int *maxsample)
 
static void InvertMatrix (double matrix[][3], const int dim, const int *)
 
static void InvertSymmetric4x4 (double W[][4])
 

Private Attributes

const CscIdHelper * m_phelper
 
ToolHandle< ICscCalibTool > m_cscCalibTool
 
double m_qerr
 
double m_terr
 
double m_qerr_fail
 
double m_terr_fail
 
double m_qerrprop
 
double m_n
 
double m_n2
 
double m_zmax
 
double m_bipolarNormalization
 
double m_tsampling
 
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default) More...
 
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default) More...
 
std::vector< SG::VarHandleKeyArray * > m_vhka
 
bool m_varHandleArraysDeclared
 

Detailed Description

Definition at line 27 of file CscBipolarStripFitter.h.

Member Typedef Documentation

◆ ChargeList

typedef std::vector<float> ICscStripFitter::ChargeList
inherited

Definition at line 55 of file ICscStripFitter.h.

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ CscBipolarStripFitter()

CscBipolarStripFitter::CscBipolarStripFitter ( const std::string &  type,
const std::string &  aname,
const IInterface *  parent 
)

Definition at line 29 of file CscBipolarStripFitter.cxx.

29  :
30  AthAlgTool(type, aname, parent), m_phelper(nullptr), m_n(0), m_n2(0), m_zmax(0), m_bipolarNormalization(0), m_tsampling(0) {
31  declareInterface<ICscStripFitter>(this);
32  declareProperty("chargeError", m_qerr = 5500.0);
33  declareProperty("timeError", m_terr = 5.0);
34  declareProperty("failChargeError", m_qerr_fail = 5000.0);
35  declareProperty("failTimeError", m_terr_fail = 50.0);
36  declareProperty("chargeCalibrationError", m_qerrprop = 0.002);
37 }

◆ ~CscBipolarStripFitter()

CscBipolarStripFitter::~CscBipolarStripFitter ( )
default

Member Function Documentation

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleKeyArrayType &   
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKeyArray>

Definition at line 170 of file AthCommonDataStore.h.

172  {
173  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
174  hndl.value(),
175  hndl.documentation());
176 
177  }

◆ declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleKeyType &   
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

158  {
159  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
160  hndl.value(),
161  hndl.documentation());
162 
163  }

◆ declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleType &   
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleBase>

Definition at line 184 of file AthCommonDataStore.h.

186  {
187  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
188  hndl.value(),
189  hndl.documentation());
190  }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  t,
const SG::NotHandleType &   
)
inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

200  {
201  return PBASE::declareProperty(t);
202  }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleBase &  hndl,
const std::string &  doc,
const SG::VarHandleType &   
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

249  {
250  this->declare(hndl.vhKey());
251  hndl.vhKey().setOwner(this);
252 
253  return PBASE::declareProperty(name,hndl,doc);
254  }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleKey &  hndl,
const std::string &  doc,
const SG::VarHandleKeyType &   
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

225  {
226  this->declare(hndl);
227  hndl.setOwner(this);
228 
229  return PBASE::declareProperty(name,hndl,doc);
230  }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleKeyArray &  hndArr,
const std::string &  doc,
const SG::VarHandleKeyArrayType &   
)
inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

263  {
264 
265  // std::ostringstream ost;
266  // ost << Algorithm::name() << " VHKA declareProp: " << name
267  // << " size: " << hndArr.keys().size()
268  // << " mode: " << hndArr.mode()
269  // << " vhka size: " << m_vhka.size()
270  // << "\n";
271  // debug() << ost.str() << endmsg;
272 
273  hndArr.setOwner(this);
274  m_vhka.push_back(&hndArr);
275 
276  Gaudi::Details::PropertyBase* p = PBASE::declareProperty(name, hndArr, doc);
277  if (p != 0) {
278  p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
279  } else {
280  ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray "
281  << name);
282  }
283 
284  return p;
285 
286  }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc,
const SG::NotHandleType &   
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

337  {
338  return PBASE::declareProperty(name, property, doc);
339  }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc = "none" 
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

355  {
356  typedef typename SG::HandleClassifier<T>::type htype;
357  return declareProperty (name, property, doc, htype());
358  }

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T > &  t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145  {
146  typedef typename SG::HandleClassifier<T>::type htype;
147  return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
148  }

◆ Derivative()

void CscBipolarStripFitter::Derivative ( double  A[][3],
double  fp[][1],
double  p0[][1],
int  imeas,
const int *  meas 
) const
private

Definition at line 267 of file CscBipolarStripFitter.cxx.

267  {
268  // calculate the derivatives and the 0th order approximation
269  // around the ADC samplings
270  double norm = p0[0][0];
271  // double parm[3]={1.,p0[1][0],p0[2][0]};
272  for (int i = 0; i < imeas; i++) {
273  int ii = meas[i];
274  double z = (ii - p0[1][0]) * m_tsampling / p0[2][0];
275  double repquant = 0.;
276  double dFdzNormalized = 0.;
277  if (z > 0.) {
278  repquant = FindPow(z) * std::exp(-z) / m_bipolarNormalization;
279  dFdzNormalized = repquant * (m_n / z + z / (m_n2 + 1) - (m_n + 1.) / (m_n2 + 1.) - 1.); // repquant*(m_n/z+z/(m_n+1)-2.);
280  }
281 
282  A[ii][0] = repquant * (1. - z / (m_n2 + 1.)); // repquant*(1.-z/(m_n+1.));
283  // A[ii][0] = bipolar(&z,parm);
284  fp[ii][0] = norm * A[ii][0];
285 
286  // double normOverZmax = norm/m_bipolarNormalization;
287  double commonpart = norm * dFdzNormalized; //(z,parm);
288  A[ii][1] = commonpart * (-m_tsampling / p0[2][0]);
289  A[ii][2] = commonpart * (-z / p0[2][0]);
290  }
291  // end of derivative/zeroth order calculations
292 }

◆ detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const
inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( ) const
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase &  ExtraDeps)
protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ FindInitValues()

double CscBipolarStripFitter::FindInitValues ( double *  x,
double *  initValues,
int *  maxsample 
)
staticprivate

Definition at line 129 of file CscBipolarStripFitter.cxx.

129  {
130  // find maximum sample imax:
131  double peakingTime = -99.; // interpolated peaking time in samples
132  double amplitude = -99.; // interpolated amplitude
133  double amin, amax;
134  int imax = 0;
135  const int numSample = 4;
136  amax = x[0];
137  amin = x[0];
138  for (int j = 1; j < numSample; j++) {
139  if (amax < x[j]) {
140  amax = x[j];
141  imax = j;
142  }
143  if (amin > x[j]) amin = x[j];
144  }
145 
146  // calculate peak and amplitude:
147  (*maxsample) = imax;
148  if ((imax == 0) || (imax == numSample - 1)) // no interpolation possible
149  {
150  peakingTime = imax;
151  amplitude = amax;
152  } else // do parabolic interpolation
153  {
154  double a, b, c; // coeffients of parabola y=a*x*x+b*x+c
155  a = 0.5 * (x[imax + 1] + x[imax - 1] - 2.0 * x[imax]);
156  b = 0.5 * (x[imax + 1] - x[imax - 1]);
157  c = x[imax];
158 
159  peakingTime = -0.5 * b / a;
160  amplitude = a * peakingTime * peakingTime + b * peakingTime + c;
161  peakingTime += imax; // it was relative to imax
162  }
163 
164  initValues[0] = amplitude;
165  initValues[1] = peakingTime;
166  // - m_zmax*initValues[2]/m_tsampling;
167  return x[imax];
168 }

◆ FindPow()

double CscBipolarStripFitter::FindPow ( double  z) const
private

Definition at line 170 of file CscBipolarStripFitter.cxx.

170  {
171  double zpower = std::exp(m_n * std::log(z));
172  return zpower;
173 }

◆ fit() [1/5]

Result ICscStripFitter::fit

Definition at line 72 of file ICscStripFitter.cxx.

16  {
17  return {};
18 }

◆ fit() [2/5]

Result ICscStripFitter::fit ( const ChargeList &  ChargeList,
double  samplingTime,
bool  samplingPhase,
Identifier &  sid 
) const
virtualinherited

Reimplemented in CalibCscStripFitter.

Definition at line 16 of file ICscStripFitter.cxx.

16  {
17  return {};
18 }

◆ fit() [3/5]

Result CscBipolarStripFitter::fit ( const ChargeList &  charges,
double  samplingTime,
Identifier &  stripId 
) const

Definition at line 74 of file CscBipolarStripFitter.cxx.

74  {
75  ATH_MSG_DEBUG("Fitting with sample time " << period);
76  ATH_MSG_DEBUG(" Number of charges =" << chgs.size() << " Charges: " << chgs);
77  Result res;
78 
79  IdentifierHash stripHash;
80  if (m_phelper->get_channel_hash(stripId, stripHash)) {
81  ATH_MSG_WARNING("Unable to get CSC striphash id "
82  << " the identifier is ");
83  stripId.show();
84  }
85  ATH_MSG_DEBUG("CalibCscStripFitter:: " << stripId << " " << (unsigned int)stripHash);
86 
91  double initValues[3];
92  for (int i = 0; i < 3; ++i) { initValues[i] = 0.; }
93  // int imax = -1;
94  initValues[2] = 8.; // m_width
96 
97  double noise = m_cscCalibTool->stripNoise(stripHash);
98  res.status = 0;
99  if (m_qerr > 0)
100  res.dcharge = m_qerr;
101  else
102  res.dcharge = noise;
103 
104  res.dtime = m_terr;
105  res.charge = initValues[0];
106  res.time = initValues[1];
107 
108  initValues[1] = res.time - 2.5;
109  // int imeas =4;
110  // int meas[4] = {true,true,true,true};
111  // int ipar = 3;
112  // int par[3] = {true,true,false};
113  // double m_chi2;
114  //double result[3];
116 
117  // Add an error proportional to the charge.
118  double dqprop = m_qerrprop * res.charge;
119  res.dcharge = sqrt(res.dcharge * res.dcharge + dqprop * dqprop);
120  // Display result.
121  ATH_MSG_DEBUG(" Status: " << res.status);
122  ATH_MSG_DEBUG(" Charge: " << res.charge);
123  ATH_MSG_DEBUG(" Charge err: " << res.dcharge);
124  ATH_MSG_DEBUG(" Time: " << res.time);
125  //ATH_MSG_DEBUG(" fit : " << result[0] << " " << result[1] << " " << result[2]);
126  return res;
127 }

◆ fit() [4/5]

Result ICscStripFitter::fit

Definition at line 79 of file ICscStripFitter.cxx.

20  {
21  Identifier sid = strip.identify();
22  // IdentifierHash coll_hash = strip.collectionHash();
23 
24  Result res = fit(strip.sampleCharges(), strip.samplingTime(), strip.samplingPhase(), sid);
25  res.strip = &strip;
26  if (res.status) return res;
27  res.time += strip.timeOfFirstSample();
28  // Do we also need a phase correction here?
29  return res;
30 }

◆ fit() [5/5]

Result ICscStripFitter::fit ( const Muon::CscStripPrepData &  strip) const
virtualinherited

Definition at line 20 of file ICscStripFitter.cxx.

20  {
21  Identifier sid = strip.identify();
22  // IdentifierHash coll_hash = strip.collectionHash();
23 
24  Result res = fit(strip.sampleCharges(), strip.samplingTime(), strip.samplingPhase(), sid);
25  res.strip = &strip;
26  if (res.status) return res;
27  res.time += strip.timeOfFirstSample();
28  // Do we also need a phase correction here?
29  return res;
30 }

◆ initialize()

StatusCode CscBipolarStripFitter::initialize ( )
override

CSC calibratin tool for the Condtiions Data base access

Definition at line 41 of file CscBipolarStripFitter.cxx.

41  {
42  ATH_MSG_DEBUG("Initializing " << name());
43 
44  ATH_MSG_DEBUG("Properties for " << name() << ":");
45  ATH_MSG_DEBUG(" Charge error: " << m_qerr);
46  ATH_MSG_DEBUG(" Time error: " << m_terr);
47  ATH_MSG_DEBUG(" Fail charge error: " << m_qerr_fail);
48  ATH_MSG_DEBUG(" Fail time error: " << m_terr_fail);
49  ATH_MSG_DEBUG(" Charge calibration error: " << m_qerrprop);
50 
52  ATH_CHECK_RECOVERABLE(m_cscCalibTool.retrieve());
53 
54  const MuonGM::MuonDetectorManager *muDetMgr = nullptr;
55  ATH_CHECK_RECOVERABLE(detStore()->retrieve(muDetMgr));
56  ATH_MSG_DEBUG("Retrieved geometry.");
57 
58  m_phelper = muDetMgr->cscIdHelper();
59 
60  m_n = m_cscCalibTool->getNumberOfIntegration(); // 12.;
61  m_n2 = m_cscCalibTool->getNumberOfIntegration2(); // 11.66;
62  m_zmax = m_n + m_n2 + 2. - sqrt((m_n + 2. + m_n2) * (m_n + 2. + m_n2) - 4. * m_n * (m_n2 + 1));
63  m_zmax *= 0.5;
64  // m_n+1 - sqrt(m_n+1.0);
65  m_bipolarNormalization = FindPow(m_zmax) * (1 - m_zmax / (m_n2 + 1)) * std::exp(-m_zmax);
66  // FindPow(m_zmax)*(1-m_zmax/(m_n+1))*std::exp(-m_zmax);
67  m_tsampling = m_cscCalibTool->getSamplingTime() / 2.; // 50/2 = 25.;
68 
69  return StatusCode::SUCCESS;
70 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const
overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ interfaceID()

static const InterfaceID& ICscStripFitter::interfaceID ( )
inlinestaticinherited

Must declare this, with name of interface

Definition at line 59 of file ICscStripFitter.h.

59  {
60  static const InterfaceID IID_ICscStripFitter("ICscStripFitter", 1, 0);
61  return IID_ICscStripFitter;
62  }

◆ InvertMatrix()

void CscBipolarStripFitter::InvertMatrix ( double  matrix[][3],
const int  dim,
const int *  correspdim 
)
staticprivate

Definition at line 175 of file CscBipolarStripFitter.cxx.

175  {
176  // invert 2x2 or 3x3 symmetric matrix
177  if (dim == 1) {
178  int ii = correspdim[0];
179  matrix[ii][ii] = 1. / matrix[ii][ii];
180 
181  } else if (dim == 2) {
182  int ii = correspdim[0];
183  int jj = correspdim[1];
184  double determinant = -matrix[jj][ii] * matrix[ii][jj] + matrix[ii][ii] * matrix[jj][jj];
185  // if(std::abs(determinant)<1.e-13)
186  // std::cout<<" zero determinant "<<std::endl;
187  double i00 = matrix[ii][ii];
188  matrix[ii][ii] = matrix[jj][jj] / determinant;
189  matrix[jj][jj] = i00 / determinant;
190  matrix[ii][jj] = -matrix[ii][jj] / determinant;
191  matrix[jj][ii] = matrix[ii][jj];
192 
193  } else if (dim == 3) {
194  double sm12 = matrix[0][1] * matrix[0][1];
195  double sm13 = matrix[0][2] * matrix[0][2];
196  double sm23 = matrix[1][2] * matrix[1][2];
197  double determinant = matrix[0][0] * matrix[1][1] * matrix[2][2] - matrix[0][0] * sm23 - sm12 * matrix[2][2] +
198  2. * matrix[0][1] * matrix[0][2] * matrix[1][2] - matrix[1][1] * sm13;
199 
200  // if(std::abs(determinant)<1.e-13)
201  // std::cout << "zero determinant"<<std::endl;
202  double i00 = matrix[1][1] * matrix[2][2] - sm23;
203  double i11 = matrix[0][0] * matrix[2][2] - sm13;
204  double i22 = matrix[0][0] * matrix[1][1] - sm12;
205 
206  matrix[1][0] = (matrix[0][2] * matrix[1][2] - matrix[2][2] * matrix[0][1]) / determinant;
207  matrix[2][0] = (matrix[0][1] * matrix[1][2] - matrix[0][2] * matrix[1][1]) / determinant;
208  matrix[2][1] = (matrix[0][1] * matrix[0][2] - matrix[0][0] * matrix[1][2]) / determinant;
209  matrix[0][1] = matrix[1][0];
210  matrix[0][2] = matrix[2][0];
211  matrix[1][2] = matrix[2][1];
212  matrix[0][0] = i00 / determinant;
213  matrix[1][1] = i11 / determinant;
214  matrix[2][2] = i22 / determinant;
215 
216  } else {
217  // int ierr;
218  // matrix.invert(ierr);
219  printf("this is not a 1x1 or 2x2 or 3x3 matrix\n");
220  }
221 }

◆ InvertSymmetric4x4()

void CscBipolarStripFitter::InvertSymmetric4x4 ( double  W[][4])
staticprivate

Definition at line 223 of file CscBipolarStripFitter.cxx.

223  {
224  double Determinant =
225  W[0][3] * W[0][3] * W[1][2] * W[1][2] - 2. * W[0][2] * W[0][3] * W[1][2] * W[1][3] + W[0][2] * W[0][2] * W[1][3] * W[1][3] -
226  W[0][3] * W[0][3] * W[1][1] * W[2][2] + 2. * W[0][1] * W[0][3] * W[1][3] * W[2][2] - W[0][0] * W[1][3] * W[1][3] * W[2][2] +
227  2. * W[0][2] * W[0][3] * W[1][1] * W[2][3] - 2. * W[0][1] * W[0][3] * W[1][2] * W[2][3] -
228  2. * W[0][1] * W[0][2] * W[1][3] * W[2][3] + 2. * W[0][0] * W[1][2] * W[1][3] * W[2][3] + W[0][1] * W[0][1] * W[2][3] * W[2][3] -
229  W[0][0] * W[1][1] * W[2][3] * W[2][3] - W[0][2] * W[0][2] * W[1][1] * W[3][3] + 2. * W[0][1] * W[0][2] * W[1][2] * W[3][3] -
230  W[0][0] * W[1][2] * W[1][2] * W[3][3] - W[0][1] * W[0][1] * W[2][2] * W[3][3] + W[0][0] * W[1][1] * W[2][2] * W[3][3];
231 
232  W[0][1] = W[3][0] * W[3][1] * W[2][2] - W[3][0] * W[2][1] * W[3][2] - W[2][0] * W[3][1] * W[3][2] + W[1][0] * W[3][2] * W[3][2] +
233  W[2][0] * W[2][1] * W[3][3] - W[1][0] * W[2][2] * W[3][3];
234  W[0][2] = -W[3][0] * W[2][1] * W[3][1] + W[2][0] * W[3][1] * W[3][1] + W[3][0] * W[1][1] * W[3][2] - W[1][0] * W[3][1] * W[3][2] -
235  W[2][0] * W[1][1] * W[3][3] + W[1][0] * W[2][1] * W[3][3];
236  W[0][3] = W[3][0] * W[2][1] * W[2][1] - W[2][0] * W[2][1] * W[3][1] - W[3][0] * W[1][1] * W[2][2] + W[1][0] * W[3][1] * W[2][2] +
237  W[2][0] * W[1][1] * W[3][2] - W[1][0] * W[2][1] * W[3][2];
238  W[1][2] = W[3][0] * W[3][0] * W[2][1] - W[2][0] * W[3][0] * W[3][1] - W[1][0] * W[3][0] * W[3][2] + W[0][0] * W[3][1] * W[3][2] +
239  W[1][0] * W[2][0] * W[3][3] - W[0][0] * W[2][1] * W[3][3];
240 
241  W[1][3] = -W[2][0] * W[3][0] * W[2][1] + W[2][0] * W[2][0] * W[3][1] + W[1][0] * W[3][0] * W[2][2] - W[0][0] * W[3][1] * W[2][2] -
242  W[1][0] * W[2][0] * W[3][2] + W[0][0] * W[2][1] * W[3][2];
243  W[2][3] = W[2][0] * W[3][0] * W[1][1] - W[1][0] * W[3][0] * W[2][1] - W[1][0] * W[2][0] * W[3][1] + W[0][0] * W[2][1] * W[3][1] +
244  W[1][0] * W[1][0] * W[3][2] - W[0][0] * W[1][1] * W[3][2];
245 
246  double W00 = -W[3][1] * W[3][1] * W[2][2] + 2. * W[2][1] * W[3][1] * W[3][2] - W[1][1] * W[3][2] * W[3][2] -
247  W[2][1] * W[2][1] * W[3][3] + W[1][1] * W[2][2] * W[3][3];
248  double W11 = -W[3][0] * W[3][0] * W[2][2] + 2. * W[2][0] * W[3][0] * W[3][2] - W[0][0] * W[3][2] * W[3][2] -
249  W[2][0] * W[2][0] * W[3][3] + W[0][0] * W[2][2] * W[3][3];
250  double W22 = -W[3][0] * W[3][0] * W[1][1] + 2. * W[1][0] * W[3][0] * W[3][1] - W[0][0] * W[3][1] * W[3][1] -
251  W[1][0] * W[1][0] * W[3][3] + W[0][0] * W[1][1] * W[3][3];
252  double W33 = -W[2][0] * W[2][0] * W[1][1] + 2. * W[1][0] * W[2][0] * W[2][1] - W[0][0] * W[2][1] * W[2][1] -
253  W[1][0] * W[1][0] * W[2][2] + W[0][0] * W[1][1] * W[2][2];
254 
255  for (int i = 0; i < 3; i++) {
256  for (int j = 1; j < 4; j++) {
257  if (i >= j) continue;
258  W[i][j] = W[i][j] / Determinant;
259  W[j][i] = W[i][j];
260  }
261  }
262  W[0][0] = W00 / Determinant;
263  W[1][1] = W11 / Determinant;
264  W[2][2] = W22 / Determinant;
265  W[3][3] = W33 / Determinant;
266 }

◆ msg() [1/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

24  {
25  return this->msgStream();
26  }

◆ msg() [2/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( const MSG::Level  lvl) const
inlineinherited

Definition at line 27 of file AthCommonMsg.h.

27  {
28  return this->msgStream(lvl);
29  }

◆ msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level  lvl) const
inlineinherited

Definition at line 30 of file AthCommonMsg.h.

30  {
31  return this->msgLevel(lvl);
32  }

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const
overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T &  h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381  {
382  h.renounce();
383  PBASE::renounce (h);
384  }

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray &  handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364  {
365  handlesArray.renounce();
366  }

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )
overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )
overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ TheFitter()

int CscBipolarStripFitter::TheFitter ( double *  x,
const double  ex,
const double *  initValues,
int  imeas,
int *  meas,
int  ipar,
int *  par,
double *  chi2,
double *  result 
) const
private

Definition at line 294 of file CscBipolarStripFitter.cxx.

295  {
296  // maximum iterations
297  const int maxIter = 7;
298  // tolerances
299  double fitTolerance0 = 0.1;
300  double fitTolerance1 = 0.01;
301  // CLHEP::HepMatrix p0(3,1,0); // the matrix of the initial fit parameters
302  double p0[3][1];
303  for (int j = 0; j < 3; j++) p0[j][0] = initValues[j];
304 
305  // CLHEP::HepMatrix m(4,1,0); // the matrix of ADC measurements (samples: 0,1,2,3)
306  double m[4][1];
307  // CLHEP::HepMatrix W(4,4,0); // the error matrix of the ADC measurements
308  double W[4][4];
309  for (int i = 0; i < 4; i++) {
310  m[i][0] = x[i];
311  W[i][i] = ex * ex;
312  }
313  // covariances
314  W[0][1] = 0.03 * ex * ex;
315  W[0][2] = -0.411 * ex * ex;
316  W[0][3] = -0.188 * ex * ex;
317  W[1][2] = 0.0275 * ex * ex;
318  W[1][3] = -0.4303 * ex * ex;
319  W[2][3] = 0. * ex * ex;
320  W[1][0] = W[0][1];
321  W[2][0] = W[0][2];
322  W[3][0] = W[0][3];
323  W[2][1] = W[1][2];
324  W[3][1] = W[1][3];
325  W[3][2] = W[2][3];
326 
327  // WW.invert(ierr);
328  InvertSymmetric4x4(W);
329 
330  // Taylor std::expansion of the bipolar pulse model around the
331  // samplings : F(x) = F(p0) + A *(p-p0) + higher.order
332  // CLHEP::HepMatrix fp(4,1,0); // the matrix of 0th order approximation
333  double fp[4][1];
334  // CLHEP::HepMatrix A(4,3,0); // the matrix of derivatives
335  double A[4][3];
336  for (int i = 0; i < 4; ++i) {
337  fp[i][0] = 0.;
338  for (int j = 0; j < 3; ++j) { A[i][j] = 0.; }
339  }
340  // remarks :
341  // if the pulse peaks in the last sampling fit with a constant shaping time
342  // if the pulse peaks in the first sampling fit without using the last sampling
343  // (too large contribution to the chi2
344  int counter = 0;
345  bool converged = false;
346  double amplitudeChangeOld = 0.;
347  bool diverganceCandidate = false;
348  // CLHEP::HepMatrix weight(3,3,1); // weight matrix allocated once
349  // the non-fitted parts are taken care appropriately
350  // at least if the fitting parameters or measurements
351  // don't change during the fitting procedure
352  double weight[3][3];
353 
354  // CLHEP::HepMatrix residFactor(3,4,0); // residFactor allocated once
355  double residFactor[3][4];
356 
357  for (int i = 0; i < 3; ++i) {
358  for (int j = 0; j < 4; ++j) {
359  if (j < 3) { weight[i][j] = 0.; }
360  residFactor[i][j] = 0.;
361  }
362  }
363  weight[0][0] = 1.;
364  weight[1][1] = 1.;
365  weight[2][2] = 1.;
366 
367  while (!converged && counter < maxIter) // fit loop
368  {
369  Derivative(A, fp, p0, imeas, meas); // calculate the matrix of derivatives and 0th order approximation
370  // matrix multiplication
371  // the weight matrix is symmetric
372  // weight= A.T()*W*A;//.assign(A.T()*W*A);
373 
374  double helpmatrix[4][3];
375  for (int i = 0; i < 4; ++i) {
376  for (int j = 0; j < 3; ++j) { helpmatrix[i][j] = 0.; }
377  }
378 
379  for (int i = 0; i < imeas; i++) {
380  int ii = meas[i];
381  for (int j = 0; j < ipar; j++) {
382  int jj = par[j];
383  for (int k = 0; k < imeas; k++) {
384  int kk = meas[k];
385  helpmatrix[ii][jj] += W[ii][kk] * A[kk][jj];
386  }
387  }
388  }
389  for (int i = 0; i < ipar; i++) {
390  int ii = par[i];
391  for (int j = i; j < ipar; j++) {
392  int jj = par[j];
393  weight[ii][jj] = 0.;
394  for (int k = 0; k < imeas; k++) {
395  int kk = meas[k];
396  weight[ii][jj] += A[kk][ii] * helpmatrix[kk][jj]; // A[kk][ii]*A[kk][jj];
397  }
398  // weight[ii][jj]*=W[0][0];
399  weight[jj][ii] = weight[ii][jj];
400  }
401  }
402  // weight.invert(ierr); // inversion of weight matrix
403  // hand-made inversion of 2x2 or 3x3 symmetric matrix
404  InvertMatrix(weight, ipar, par);
405 
406  // calculate W*(A.T()*W)
407  // residFactor = weight*(A.T()*W);
408  double helpmatrix2[3][4];
409  for (int i = 0; i < 3; ++i) {
410  for (int j = 0; j < 4; ++j) { helpmatrix2[i][j] = 0.; }
411  }
412 
413  for (int i = 0; i < ipar; i++) {
414  int ii = par[i];
415  for (int j = 0; j < imeas; j++) {
416  int jj = meas[j];
417  for (int k = 0; k < imeas; k++) {
418  int kk = meas[k];
419  helpmatrix2[ii][jj] += A[kk][ii] * W[kk][jj];
420  }
421  }
422  }
423 
424  for (int i = 0; i < ipar; i++) {
425  int ii = par[i];
426  for (int j = 0; j < imeas; j++) {
427  int jj = meas[j];
428  residFactor[ii][jj] = 0.;
429  for (int k = 0; k < ipar; k++) {
430  int kk = par[k];
431  residFactor[ii][jj] += weight[ii][kk] * helpmatrix2[kk][jj];
432  }
433  // residFactor[ii][jj]*=W[0][0];
434  }
435  }
436 
437  double paramDiff[3];
438  for (int i = 0; i < 3; ++i) { paramDiff[i] = 0.; }
439 
440  for (int i = 0; i < ipar; i++) {
441  int ii = par[i];
442  // estimation of new parameters
443  // paramDiff[i][0] += (weight*(A.T()*W)*(m-fp))[i][0];
444  for (int j = 0; j < imeas; j++) {
445  int jj = meas[j];
446  paramDiff[ii] += residFactor[ii][jj] * (m[jj][0] - fp[jj][0]);
447  }
448  p0[ii][0] += paramDiff[ii];
449  }
450  std::cout << "##### " << p0[0][0] << " " << p0[1][0] << " " << p0[2][0] << std::endl;
451  // if the parameters are not physical, keep them sensible
452  // if peaking time less than -0.5
453  double peakingTime = p0[1][0] + m_zmax * p0[2][0] / m_tsampling;
454  if (peakingTime < -0.5 || peakingTime > 3.) p0[1][0] = initValues[1];
455 
456  if (p0[0][0] < 0.) p0[0][0] = initValues[0];
457  double amplitudeChangeNew = std::abs(paramDiff[0]);
458  if (std::abs(paramDiff[0]) < fitTolerance0 && std::abs(paramDiff[1]) < fitTolerance1) {
459  converged = true;
460  // calculate chi2
461  // (m-fp).T()*W*(m-fp)
462  double residual[4];
463  for (int i = 0; i < 4; ++i) { residual[i] = 0.; }
464 
465  for (int i = 0; i < imeas; i++) {
466  int ii = meas[i];
467  residual[i] = m[ii][0] - fp[ii][0];
468  }
469 
470  double tmpChi2 = 0.;
471  double helpmatrixchi2[4][1];
472  for (int i = 0; i < 4; ++i) { helpmatrixchi2[i][0] = 0.; }
473  for (int i = 0; i < imeas; i++) {
474  int ii = meas[i];
475  for (int k = 0; k < imeas; k++) {
476  int kk = meas[k];
477  helpmatrixchi2[ii][0] += W[ii][kk] * residual[kk];
478  }
479  }
480  for (int k = 0; k < imeas; k++) {
481  int kk = meas[k];
482  tmpChi2 += residual[kk] * helpmatrixchi2[kk][0];
483  // std::cout << residual[kk] << " ";//*W[kk][kk]*residual[kk]<< " ";
484  }
485  // std::cout<<std::endl;
486  (*chi2) = tmpChi2;
487  } else if (counter > 4 && (amplitudeChangeNew > 4. * amplitudeChangeOld)) {
488  if (diverganceCandidate) {
489  // diverging fit
490  // return parabola interpolation
491  printf("%3.2f %3.2f %3.2f %3.2f\n ", x[0], x[1], x[2], x[3]);
492  printf("Diverging fit\n");
493  return 4;
494  } else
495  diverganceCandidate = true;
496  }
497  if (p0[0][0] < 0.) {
498  // negative amplitude
499  // fit diverged
500  // return parabola
501  return 4;
502  }
503  // if after a couple of iterations the amplitude is low
504  // reduce the tolerances (or the maximum iterations)
505  // low amplitude pulses tend to oscillate and exhaust all iterations
506  if (p0[0][0] < 20.) {
507  fitTolerance0 = 0.1;
508  fitTolerance1 = 0.05;
509  }
510  amplitudeChangeOld = amplitudeChangeNew;
511  counter++;
512  }
513  /*
514  std::cout << " Error matrix "<<std::endl;
515  for(int i=0;i<3;i++)
516  {
517  for(int j=0;j<3;j++)
518  std::cout << weight[i][j]<<"\t";
519  std::cout<<std::endl;
520  }
521  */
522 
523  result[0] = p0[0][0];
524  result[1] = m_zmax * p0[2][0] / m_tsampling + p0[1][0];
525  result[2] = p0[2][0];
526 
527  if (counter == maxIter) return 3;
528  return 0;
529 }

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase &  )
inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

308  {
309  // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310  // << " size: " << m_vhka.size() << endmsg;
311  for (auto &a : m_vhka) {
312  std::vector<SG::VarHandleKey*> keys = a->keys();
313  for (auto k : keys) {
314  k->setOwner(this);
315  }
316  }
317  }

Member Data Documentation

◆ m_bipolarNormalization

double CscBipolarStripFitter::m_bipolarNormalization
private

Definition at line 74 of file CscBipolarStripFitter.h.

◆ m_cscCalibTool

ToolHandle<ICscCalibTool> CscBipolarStripFitter::m_cscCalibTool
private
Initial value:
{
this,
"cscCalibTool",
"",
}

Definition at line 51 of file CscBipolarStripFitter.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_n

double CscBipolarStripFitter::m_n
private

Definition at line 71 of file CscBipolarStripFitter.h.

◆ m_n2

double CscBipolarStripFitter::m_n2
private

Definition at line 72 of file CscBipolarStripFitter.h.

◆ m_phelper

const CscIdHelper* CscBipolarStripFitter::m_phelper
private

Definition at line 48 of file CscBipolarStripFitter.h.

◆ m_qerr

double CscBipolarStripFitter::m_qerr
private

Definition at line 66 of file CscBipolarStripFitter.h.

◆ m_qerr_fail

double CscBipolarStripFitter::m_qerr_fail
private

Definition at line 68 of file CscBipolarStripFitter.h.

◆ m_qerrprop

double CscBipolarStripFitter::m_qerrprop
private

Definition at line 70 of file CscBipolarStripFitter.h.

◆ m_terr

double CscBipolarStripFitter::m_terr
private

Definition at line 67 of file CscBipolarStripFitter.h.

◆ m_terr_fail

double CscBipolarStripFitter::m_terr_fail
private

Definition at line 69 of file CscBipolarStripFitter.h.

◆ m_tsampling

double CscBipolarStripFitter::m_tsampling
private

Definition at line 75 of file CscBipolarStripFitter.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_zmax

double CscBipolarStripFitter::m_zmax
private

Definition at line 73 of file CscBipolarStripFitter.h.

The documentation for this class was generated from the following files:
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
CscBipolarStripFitter::m_zmax
double m_zmax
Definition: CscBipolarStripFitter.h:73
python.AtlRunQueryAMI.period
period
Definition: AtlRunQueryAMI.py:225
CscBipolarStripFitter::fit
Result fit(const ChargeList &charges, double samplingTime, Identifier &stripId) const
Definition: CscBipolarStripFitter.cxx:74
CscBipolarStripFitter::m_cscCalibTool
ToolHandle< ICscCalibTool > m_cscCalibTool
Definition: CscBipolarStripFitter.h:51
PlotCalibFromCool.norm
norm
Definition: PlotCalibFromCool.py:100
CscBipolarStripFitter::m_terr_fail
double m_terr_fail
Definition: CscBipolarStripFitter.h:69
yodamerge_tmp.dim
dim
Definition: yodamerge_tmp.py:239
get_generator_info.result
result
Definition: get_generator_info.py:21
SiliconTech::strip
@ strip
python.SystemOfUnits.m
int m
Definition: SystemOfUnits.py:91
CscBipolarStripFitter::m_n2
double m_n2
Definition: CscBipolarStripFitter.h:72
ClusterSeg::residual
@ residual
Definition: ClusterNtuple.h:20
CscBipolarStripFitter::m_tsampling
double m_tsampling
Definition: CscBipolarStripFitter.h:75
JetTiledMap::W
@ W
Definition: TiledEtaPhiMap.h:44
AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
ICscStripFitter::fit
virtual Result fit(const ChargeList &ChargeList, double samplingTime, bool samplingPhase, Identifier &sid) const
Definition: ICscStripFitter.cxx:16
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
ATH_CHECK_RECOVERABLE
#define ATH_CHECK_RECOVERABLE
Evaluate an expression and check for errors.
Definition: AthCheckMacros.h:48
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
drawFromPickle.exp
exp
Definition: drawFromPickle.py:36
x
#define x
AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore
const ServiceHandle< StoreGateSvc > & detStore() const
The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:95
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
python.Dumpers.aname
string aname
Definition: Dumpers.py:5546
dqt_zlumi_pandas.weight
int weight
Definition: dqt_zlumi_pandas.py:189
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:27
python.changerun.kk
list kk
Definition: changerun.py:41
A
CscBipolarStripFitter::m_n
double m_n
Definition: CscBipolarStripFitter.h:71
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
AthCommonDataStore
Definition: AthCommonDataStore.h:52
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
lumiFormat.i
int i
Definition: lumiFormat.py:85
z
#define z
trigmenu_modify_prescale_json.fp
fp
Definition: trigmenu_modify_prescale_json.py:53
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
CscBipolarStripFitter::m_qerrprop
double m_qerrprop
Definition: CscBipolarStripFitter.h:70
res
std::pair< std::vector< unsigned int >, bool > res
Definition: JetGroupProductTest.cxx:14
CscBipolarStripFitter::m_bipolarNormalization
double m_bipolarNormalization
Definition: CscBipolarStripFitter.h:74
test_pyathena.parent
parent
Definition: test_pyathena.py:15
CscBipolarStripFitter::Derivative
void Derivative(double A[][3], double fp[][1], double p0[][1], int imeas, const int *meas) const
Definition: CscBipolarStripFitter.cxx:267
CscBipolarStripFitter::m_qerr
double m_qerr
Definition: CscBipolarStripFitter.h:66
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
CscBipolarStripFitter::InvertSymmetric4x4
static void InvertSymmetric4x4(double W[][4])
Definition: CscBipolarStripFitter.cxx:223
imax
int imax(int i, int j)
Definition: TileLaserTimingTool.cxx:33
Identifier::show
void show() const
Print out in hex form.
Definition: Identifier.cxx:30
AthAlgTool::AthAlgTool
AthAlgTool()
Default constructor:
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
MuonIdHelper::get_channel_hash
virtual int get_channel_hash(const Identifier &id, IdentifierHash &hash_id) const
Definition: MuonIdHelper.cxx:139
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
Result
Definition: fbtTestBasics.cxx:49
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:221
createCoolChannelIdFile.par
par
Definition: createCoolChannelIdFile.py:29
CscBipolarStripFitter::m_qerr_fail
double m_qerr_fail
Definition: CscBipolarStripFitter.h:68
plotBeamSpotMon.b
b
Definition: plotBeamSpotMon.py:77
CscBipolarStripFitter::InvertMatrix
static void InvertMatrix(double matrix[][3], const int dim, const int *)
Definition: CscBipolarStripFitter.cxx:175
CscBipolarStripFitter::m_phelper
const CscIdHelper * m_phelper
Definition: CscBipolarStripFitter.h:48
python.testIfMatch.matrix
matrix
Definition: testIfMatch.py:66
a
TList * a
Definition: liststreamerinfos.cxx:10
h
MuonGM::MuonDetectorManager
The MuonDetectorManager stores the transient representation of the Muon Spectrometer geometry and pro...
Definition: MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonDetectorManager.h:50
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
python.CaloScaleNoiseConfig.type
type
Definition: CaloScaleNoiseConfig.py:78
python.CaloCondTools.log
log
Definition: CaloCondTools.py:20
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:623
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:798
MuonGM::MuonDetectorManager::cscIdHelper
const CscIdHelper * cscIdHelper() const
Definition: MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonDetectorManager.h:232
CscBipolarStripFitter::m_terr
double m_terr
Definition: CscBipolarStripFitter.h:67
IdentifierHash
This is a "hash" representation of an Identifier. This encodes a 32 bit index which can be used to lo...
Definition: IdentifierHash.h:25
test_pyathena.counter
counter
Definition: test_pyathena.py:15
CscBipolarStripFitter::FindPow
double FindPow(double z) const
Definition: CscBipolarStripFitter.cxx:170
python.compressB64.c
def c
Definition: compressB64.py:93
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
WriteCellNoiseToCool.noise
noise
Definition: WriteCellNoiseToCool.py:380
TRTCalib_cfilter.p0
p0
Definition: TRTCalib_cfilter.py:129
fitman.k
k
Definition: fitman.py:528
Identifier
Definition: IdentifierFieldParser.cxx:14
Generated on Thu Nov 7 2024 21:40:54 for ATLAS Offline Software by   doxygen 1.8.18