ATLAS Offline Software
Classes | 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
Rec::MuidCaloEnergyParam Class Reference

#include <MuidCaloEnergyParam.h>

Inheritance diagram for Rec::MuidCaloEnergyParam:
Collaboration diagram for Rec::MuidCaloEnergyParam:

Classes

struct  BinsWeights
 

Public Member Functions

 MuidCaloEnergyParam (const std::string &type, const std::string &name, const IInterface *parent)
 
 ~MuidCaloEnergyParam (void)
 
StatusCode initialize ()
 
CaloEnergy * meanParametrizedEnergy (double trackMomentum, double eta, double phi) const
 IMuidCaloEnergyParam interface: to get the mean muon energy loss in the calorimeter (with asymmetric errors describing the whole distribution) More...
 
CaloEnergy * mopParametrizedEnergy (double trackMomentum, double eta, double phi) const
 IMuidCaloEnergyParam interface: to get the most probable (peak value) energy loss in the calorimeter (with asymmetric errors describing the whole distribution)
More...
 
CaloEnergy * mopPeakEnergy (double trackMomentum, double eta, double phi) const
 IMuidCaloEnergyParam interface: to get the most probable (peak value) energy loss in the calorimeter (with asymmetric errors describing 2 sigma around the peak region)
More...
 
CaloEnergy * mopDepositedEnergy (double trackMomentum, double eta, double phi) const
 IMuidCaloEnergyParam interface: to get the most probable (peak value) energy loss and error in the calorimeter. More...
 
double x0mapInertMaterial (double) const
 IMuidCaloEnergyParam interface: returns the fraction of inert material between InDet and MS. More...
 
double x0mapEmMaterial (double) const
 IMuidCaloEnergyParam interface: returns the fraction of em calorimeter inert material between InDet and MS. More...
 
double x0mapHecMaterial (double) const
 IMuidCaloEnergyParam interface: returns the fraction of hec calorimeter inert material between InDet and MS. More...
 
double caloCompartmentDepth (int, int) const
 IMuidCaloEnergyParam interface: depth of calo compartment. More...
 
double emMopFraction (double eta) const
 IMuidCaloEnergyParam interface: fraction of the full mop deposition corresponding to the em. More...
 
double em1MopFraction (double eta) const
 IMuidCaloEnergyParam interface: fraction of the full mop deposition corresponding to the em presampler + compartment #1. More...
 
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 ()
 AlgTool and IAlgTool interface methods. More...
 

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

BinsWeights etaFixedBin (double) const
 
BinsWeights etaVariableBin (double) const
 
double meanEnergyLoss (double, double) const
 
std::pair< double, double > meanEnergyLossError (double, double) const
 
double mopEnergyLoss (double, double) const
 
double mopEnergyLossError (double, double) const
 
double symmetricMopEnergyLoss (double, double) 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 caloCompartmentDepthTile (int)
 
static double caloCompartmentDepthLArHEC (int)
 
static double caloCompartmentDepthLArEM (int)
 

Private Attributes

bool m_cosmics
 
double m_smoothingFraction
 
double m_binWidth
 
double m_etaOffset
 
double m_inverseWidth
 
double m_meanEnergyLossP0 [26] {}
 
double m_meanEnergyLossP1 [26] {}
 
double m_meanEnergyLossP2 [26] {}
 
double m_meanEnergyLossErrorLeftP0 [26] {}
 
double m_meanEnergyLossErrorLeftP1 [26] {}
 
double m_meanEnergyLossErrorRightP0 [26] {}
 
double m_meanEnergyLossErrorRightP1 [26] {}
 
double m_mopEnergyLossP0 [26] {}
 
double m_mopEnergyLossP1 [26] {}
 
double m_mopEnergyLossP2 [26] {}
 
double m_mopEnergyLossErrorP0 [26] {}
 
double m_mopEnergyLossErrorP1 [26] {}
 
double m_mopSymmetricEnergyLossP0 [26] {}
 
double m_mopSymmetricEnergyLossP1 [26] {}
 
double m_mopSymmetricEnergyLossP2 [26] {}
 
double m_emMopFraction [26] {}
 
double m_etaGranularity
 
double m_dead_eta [277] {}
 
double m_dead_x01 [277] {}
 
double m_dead_x02 [277] {}
 
double m_dead_x0tot [277] {}
 
double m_dead_cryo0 [277] {}
 
double m_dead_cryo1 [277] {}
 
double m_dead_cryo2 [277] {}
 
double m_dead_girder [277] {}
 
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 28 of file MuidCaloEnergyParam.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ MuidCaloEnergyParam()

MuidCaloEnergyParam::MuidCaloEnergyParam ( const std::string &  type,
const std::string &  name,
const IInterface *  parent 
)

Definition at line 31 of file MuidCaloEnergyParam.cxx.

31  :
32  AthAlgTool(type, name, parent),
33  m_cosmics(false),
34  m_smoothingFraction(0.5),
35  m_binWidth(0.),
36  m_etaOffset(0.),
37  m_inverseWidth(0.),
38  m_etaGranularity(0.) {
39  declareInterface<IMuidCaloEnergyParam>(this);
40  declareProperty("Cosmics", m_cosmics);
41  declareProperty("SmoothingFraction", m_smoothingFraction);
42  }

◆ ~MuidCaloEnergyParam()

MuidCaloEnergyParam::~MuidCaloEnergyParam ( void  )
default

Member Function Documentation

◆ caloCompartmentDepth()

double MuidCaloEnergyParam::caloCompartmentDepth ( int  ,
int   
) const
virtual

IMuidCaloEnergyParam interface: depth of calo compartment.

Implements Rec::IMuidCaloEnergyParam.

Definition at line 658 of file MuidCaloEnergyParam.cxx.

658  {
659  /* ------------------
660  icalo calorimeter
661  0 Tile
662  1 LArHEC
663  2 LArEM
664  ------------------ */
665  if (icalo == 0)
666  return caloCompartmentDepthTile(icomp);
667  else if (icalo == 1)
668  return caloCompartmentDepthLArHEC(icomp);
669  else if (icalo == 2)
670  return caloCompartmentDepthLArEM(icomp);
671  else
672  ATH_MSG_WARNING(" Incorrect selection of calorimeter... returning 0.");
673  return 0.;
674  }

◆ caloCompartmentDepthLArEM()

double MuidCaloEnergyParam::caloCompartmentDepthLArEM ( int  icomp)
staticprivate

Definition at line 695 of file MuidCaloEnergyParam.cxx.

695  {
696  // depth in radiation lengths (from Liquid Argon TDR p 5) eta = 0
697  constexpr double comp[4] = {0.2, 4.3, 16., 2.};
698  double total = 0.;
699  for (int i = 0; i < 4; i++) total += comp[i];
700  if (icomp >= 0 && icomp < 4) { return comp[icomp] / total; }
701  return 0.;
702  }

◆ caloCompartmentDepthLArHEC()

double MuidCaloEnergyParam::caloCompartmentDepthLArHEC ( int  icomp)
staticprivate

Definition at line 685 of file MuidCaloEnergyParam.cxx.

685  {
686  // depth in absorption lengths (from Liquid Argon TDR Table 1-3 p 11)
687  constexpr double comp[4] = {1.4, 2.9, 5.7 / 2., 5.7 / 2.}; // divide 5.7 rad lengths equally to 3rd and 4th sampling
688  double total = 0.;
689  for (int i = 0; i < 4; i++) total += comp[i];
690  const double inv_total = 1. / total;
691  if (icomp >= 0 && icomp < 4) { return comp[icomp] * inv_total; }
692  return 0.;
693  }

◆ caloCompartmentDepthTile()

double MuidCaloEnergyParam::caloCompartmentDepthTile ( int  icomp)
staticprivate

Definition at line 676 of file MuidCaloEnergyParam.cxx.

676  {
677  // depth in absorption lengths (from Tile TDR p 22) eta = 0
678  constexpr double comp[4] = {1.5, 4.2, 1.9, 0.}; // assign 0. to ITC for the moment
679  double total = 0.;
680  for (int i = 0; i < 4; i++) total += comp[i];
681  if (icomp >= 0 && icomp < 4) { return comp[icomp] / total; }
682  return 0.;
683  }

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

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

◆ em1MopFraction()

double MuidCaloEnergyParam::em1MopFraction ( double  eta) const
virtual

IMuidCaloEnergyParam interface: fraction of the full mop deposition corresponding to the em presampler + compartment #1.

Implements Rec::IMuidCaloEnergyParam.

Definition at line 711 of file MuidCaloEnergyParam.cxx.

711  {
712  // FIXME!!
713  double em1Fraction = 0.1 * emMopFraction(eta);
714  return em1Fraction;
715  }

◆ emMopFraction()

double MuidCaloEnergyParam::emMopFraction ( double  eta) const
virtual

IMuidCaloEnergyParam interface: fraction of the full mop deposition corresponding to the em.

Implements Rec::IMuidCaloEnergyParam.

Definition at line 704 of file MuidCaloEnergyParam.cxx.

704  {
705  int ieta = static_cast<int>(std::abs(eta) / 0.10);
706  if (ieta > 25) ieta = 25;
707 
708  return m_emMopFraction[ieta];
709  }

◆ etaFixedBin()

MuidCaloEnergyParam::BinsWeights MuidCaloEnergyParam::etaFixedBin ( double  eta) const
private

Definition at line 460 of file MuidCaloEnergyParam.cxx.

460  {
461  MuidCaloEnergyParam::BinsWeights BW{};
462  eta = std::abs(eta);
463  double offsetEta1 = eta - m_etaOffset;
464  double offsetEta2 = eta + m_etaOffset;
465  BW.etaWeight2 = 0.5;
466  if (offsetEta2 < m_binWidth) // first (i.e. central) bin #0 - no interpolation
467  {
468  BW.etaBin1 = 0;
469  BW.etaBin2 = 0;
470  } else if (offsetEta1 >= 2.50) // last bin #25 (eta > 2.5) - no interpolation
471  {
472  BW.etaBin1 = 25;
473  BW.etaBin2 = 25;
474  } else if (m_smoothingFraction > 0.) {
475  double eta1 = floor(m_inverseWidth * offsetEta1);
476  BW.etaBin1 = static_cast<int>(eta1);
477  double eta2 = floor(m_inverseWidth * offsetEta2);
478  BW.etaBin2 = static_cast<int>(eta2);
479  if (BW.etaBin2 > BW.etaBin1) BW.etaWeight2 = (m_inverseWidth * offsetEta2 - eta2) / m_smoothingFraction;
480  } else {
481  double eta1 = floor(m_inverseWidth * eta);
482  BW.etaBin1 = static_cast<int>(eta1);
483  BW.etaBin2 = BW.etaBin1;
484  }
485 
486  BW.etaWeight1 = 1. - BW.etaWeight2;
487  return BW;
488  }

◆ etaVariableBin()

MuidCaloEnergyParam::BinsWeights MuidCaloEnergyParam::etaVariableBin ( double  eta) const
private

Definition at line 490 of file MuidCaloEnergyParam.cxx.

490  {
491  eta = std::abs(eta);
492  if (eta < 1.05) // barrel: contains the variable bin region
493  {
494  if (eta < 0.75) // identical to fixed binning
495  {
496  return etaFixedBin(eta);
497  } else if (eta < 0.80) // wide bin#7
498  {
499  return etaFixedBin(0.75);
500  } else if (eta < 1.00) // similar to fixed binning (with half bin offset)
501  {
502  return etaFixedBin(eta - 0.5 * m_binWidth);
503  } else // wide bin#9
504  {
505  return etaFixedBin(0.95);
506  }
507  } else // endcap: fixed binning (with full bin offset)
508  {
509  return etaFixedBin(eta - m_binWidth);
510  }
511  }

◆ 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

◆ initialize()

StatusCode MuidCaloEnergyParam::initialize ( )

Definition at line 48 of file MuidCaloEnergyParam.cxx.

48  {
49  ATH_MSG_DEBUG(" Initializing ");
50 
51  // set up eta bin width and precomputed smoothing constants
52  m_binWidth = 0.1;
53  if (m_smoothingFraction < 0.) m_smoothingFraction = 0.;
54  if (m_smoothingFraction > 1.) m_smoothingFraction = 1.;
55  m_etaOffset = 0.5 * m_binWidth * m_smoothingFraction;
56  m_inverseWidth = 1. / m_binWidth;
57 
58  // fix to constant term for ATLAS-GEO-10-00-00 --> derived using a pt=5GeV
59  // single muon sample comparing the actual eloss and the expected according
60  // to the parametrisation after comparison with data, this is only to be
61  // applied to barrel region
62  constexpr double fix[26] = {0.21892, 0.13412, 0.064942, 0.154388, 0.16903, 0.212424, 0.246268, 0.234471, 0.23724,
63  0.209337, 0.186304, 0.100489, 0.0740633, 0.0319171, 0.0352365, 0.0218442, 0.19691, 0.124486,
64  0.0882254, 0.132367, 0.181742, 0.213001, 0.192657, 0.219348, 0.261353, 0.142727};
65 
66  // Initialization of the mean energy deposition arrays
67  //(function [0]+[1]*Log(0.0067*p*p)+[2]*p with GeV units)
68  constexpr double meanEnergyLoss_p0[26] = {2.7909, 2.82994, 2.91726, 3.00969, 3.09838, 3.25242, 3.32701, 2.99668, 3.05648,
69  3.50282, 4.41603, 4.64789, 3.90948, 3.29712, 2.86013, 2.8295, 3.13659, 3.6972,
70  3.69033, 3.76311, 3.62784, 3.60054, 3.67615, 3.57151, 3.59287, 3.40605};
71  constexpr double meanEnergyLoss_p1[26] = {0.0497463, 0.0440959, 0.0331098, 0.0565649, 0.0721456, 0.0948259, 0.0825271,
72  0.0543386, 0.0790654, 0.0494549, 0.118818, 0.139855, 0.00587402, 0.0589792,
73  0.0461476, 0.0285649, 0.0460483, 0.130774, 0.0999097, 0.126468, 0.0302406,
74  0.0645617, 0.0416555, 0.0503394, -0.00448376, 0.135223};
75  constexpr double meanEnergyLoss_p2[26] = {0.0102245, 0.0113924, 0.0090614, 0.0105087, 0.0104329, 0.010759, 0.0103308,
76  0.011115, 0.00944668, 0.0126312, 0.0127012, 0.0162508, 0.0131209, 0.00842482,
77  0.00889747, 0.00755592, 0.00944529, 0.0106429, 0.0111037, 0.0113795, 0.0129818,
78  0.0128469, 0.0138801, 0.015034, 0.0135201, 0.0109176};
79  // tabulated coefficients for left assymetric error: pol1 with GeV units
80  constexpr double meanEnergyLoss_left_p0[26] = {0.290254, 0.253085, 0.252006, 0.279803, 0.284166, 0.269267, 0.370476,
81  0.309723, 0.339191, 0.392122, 0.504243, 0.452033, 0.389579, 0.325196,
82  0.241338, 0.230355, 0.319177, 0.341208, 0.262456, 0.300411, 0.255884,
83  0.257605, 0.272939, 0.247617, 0.267175, 0.377071};
84  constexpr double meanEnergyLoss_left_p1[26] = {0.00570519, 0.00594466, 0.00519917, 0.0054666, 0.00552764, 0.00631414, 0.00539921,
85  0.00573473, 0.00508457, 0.00624818, 0.00691296, 0.00864038, 0.00640888, 0.00488251,
86  0.00507001, 0.00462523, 0.00529166, 0.00594752, 0.00640289, 0.00655527, 0.00631339,
87  0.00643479, 0.00749501, 0.00766192, 0.00637578, 0.00567618};
88  // tabulated coefficients for right assymetric error: pol1 with GeV units
89  constexpr double meanEnergyLoss_right_p0[26] = {0.478493, 0.587466, 0.678961, 0.723157, 0.536312, 0.60896, 0.570827,
90  0.641256, 0.713216, 0.648768, 0.684193, 0.472149, 0.697594, 0.648713,
91  0.577021, 0.52603, 0.536587, 0.573851, 0.603372, 0.817339, 0.631103,
92  0.491342, 0.546488, 0.46241, 0.46241, 0.951949};
93  constexpr double meanEnergyLoss_right_p1[26] = {0.0250521, 0.0284547, 0.0202018, 0.0202702, 0.0257285, 0.0230366, 0.0240065,
94  0.0209904, 0.0220515, 0.0296286, 0.0314476, 0.0452839, 0.0343686, 0.0197284,
95  0.0217681, 0.018329, 0.0193914, 0.0305364, 0.0258112, 0.0272972, 0.0279697,
96  0.0228007, 0.0268285, 0.0312118, 0.0312118, 0.0205223};
97 
98  for (int i = 0; i < 26; ++i) {
99  if (i < 13) {
100  m_meanEnergyLossP0[i] = meanEnergyLoss_p0[i] + fix[i];
101  } else {
102  m_meanEnergyLossP0[i] = meanEnergyLoss_p0[i];
103  }
104  m_meanEnergyLossP1[i] = meanEnergyLoss_p1[i];
105  m_meanEnergyLossP2[i] = meanEnergyLoss_p2[i];
106  m_meanEnergyLossErrorLeftP0[i] = meanEnergyLoss_left_p0[i];
107  m_meanEnergyLossErrorLeftP1[i] = meanEnergyLoss_left_p1[i];
108  m_meanEnergyLossErrorRightP0[i] = meanEnergyLoss_right_p0[i];
109  m_meanEnergyLossErrorRightP1[i] = meanEnergyLoss_right_p1[i];
110  }
111 
112  // Initialization of the mop energy deposition arrays
113  // tabulated coefficients for MOP energy loss and for its error
114  // (function [0]+[1]*Log(0.0067*p*p)+[2]*p with GeV units)
115  constexpr double mopEnergyLoss_p0[26] = {2.57291, 2.57581, 2.65592, 2.73852, 2.81242, 2.94167, 3.01323, 2.72711, 2.75402,
116  3.17494, 4.07206, 4.28525, 3.58247, 2.99343, 2.59807, 2.56842, 2.83854, 3.4031,
117  3.40046, 3.3927, 3.38412, 3.38123, 3.4004, 3.37481, 3.3171, 3.28507};
118  constexpr double mopEnergyLoss_p1[26] = {0.0464874, 0.0213596, 0.0113974, 0.0314924, 0.0283828, 0.0363767, 0.0405601,
119  0.0211039, 0.0273863, 0.0351794, 0.0381777, 0.079885, 0.0176414, 0.0280359,
120  0.0173828, 0.0280728, 0.0366678, 0.0416954, 0.029448, 0.0300134, 0.0262876,
121  0.0229897, 0.0386583, 0.0305293, 0.0214848, 0.0316924};
122  constexpr double mopEnergyLoss_p2[26] = {0.00172448, 0.00225414, 0.00225085, 0.00196091, 0.00229175, 0.00232876, 0.00248703,
123  0.00239371, 0.00204022, 0.00238259, 0.00387238, 0.00398731, 0.0030019, 0.00220271,
124  0.00173978, 0.00153614, 0.00224971, 0.00291137, 0.00297603, 0.00313257, 0.00330651,
125  0.00319543, 0.00287297, 0.00333396, 0.00301831, 0.0027259};
126 
127  constexpr double mopEnergyLoss_sigma_p0[26] = {0.129825, 0.120123, 0.110096, 0.129288, 0.13066, 0.128621, 0.139855,
128  0.136144, 0.143521, 0.192355, 0.22055, 0.24224, 0.195584, 0.147344,
129  0.121476, 0.112492, 0.175842, 0.127176, 0.12163, 0.123569, 0.116129,
130  0.118597, 0.126694, 0.129753, 0.128578, 0.109232};
131  constexpr double mopEnergyLoss_sigma_p1[26] = {0.00117338, 0.00114104, 0.00116575, 0.00106617, 0.00121612, 0.0013242, 0.00133937,
132  0.00125934, 0.00114077, 0.00122792, 0.00190011, 0.00203161, 0.00127366, 0.00110695,
133  0.000916009, 0.000920279, 0.00131875, 0.00147309, 0.0015842, 0.0016292, 0.00166247,
134  0.00151376, 0.00149124, 0.00163482, 0.00156624, 0.00148398};
135 
136  for (int i = 0; i < 26; ++i) {
137  if (i < 13) {
138  m_mopEnergyLossP0[i] = mopEnergyLoss_p0[i] + fix[i];
139  } else {
140  m_mopEnergyLossP0[i] = mopEnergyLoss_p0[i];
141  }
142  m_mopEnergyLossP1[i] = mopEnergyLoss_p1[i];
143  m_mopEnergyLossP2[i] = mopEnergyLoss_p2[i];
144  m_mopEnergyLossErrorP0[i] = mopEnergyLoss_sigma_p0[i];
145  m_mopEnergyLossErrorP1[i] = mopEnergyLoss_sigma_p1[i];
146  }
147 
148  // not sure the mop symmetric is still needed...
149  constexpr double mopSymmetricEnergyLoss_p0[26] = {2.70345, 2.68656, 2.75213, 2.85555, 2.95905, 3.09844, 3.19676, 2.87172, 2.92657,
150  3.38863, 4.28484, 4.54497, 3.77225, 3.18316, 2.72541, 2.68793, 3.08594, 3.52879,
151  3.5142, 3.48087, 3.50184, 3.4893, 3.52484, 3.48467, 3.44116, 3.43226};
152  constexpr double mopSymmetricEnergyLoss_p1[26] = {
153  0.0388166, 0.00777747, 0.0261084, 0.0406787, 0.0322767, 0.038223, 0.0409763, 0.00096891, 0.0439047,
154  0.0430484, 0.0419927, 0.0766375, -0.000214255, 0.0420977, -0.0143302, -0.0022226, 0.0574891, 0.0260815,
155  -0.00223455, 0.0016692, -0.00404322, -0.00620859, 0.00603695, -0.0271089, 0.0386722, 0.0273018};
156  constexpr double mopSymmetricEnergyLoss_p2[26] = {
157  0.00263823, 0.00360428, 0.00293494, 0.00235128, 0.00245869, 0.00304957, 0.00317534, 0.00428433, 0.00231015,
158  0.0029778, 0.00518089, 0.00522931, 0.00446165, 0.0023943, 0.00365962, 0.00310547, 0.00236659, 0.00405503,
159  0.00488483, 0.00499606, 0.0050601, 0.0045985, 0.00515479, 0.00607963, 0.00405829, 0.00440702};
160  for (int i = 0; i < 26; ++i) {
161  m_mopSymmetricEnergyLossP0[i] = mopSymmetricEnergyLoss_p0[i];
162  m_mopSymmetricEnergyLossP1[i] = mopSymmetricEnergyLoss_p1[i];
163  m_mopSymmetricEnergyLossP2[i] = mopSymmetricEnergyLoss_p2[i];
164  }
165 
166  // Initialize the em energy deposit fraction
167  constexpr double emMopFraction[26] = {0.100105, 0.101756, 0.1033, 0.103466, 0.10543, 0.103732, 0.106252, 0.120923, 0.116363,
168  0.116984, 0.118591, 0.104436, 0.100136, 0.0934386, 0.104249, 0.122974, 0.122026, 0.113485,
169  0.0978895, 0.0967849, 0.100603, 0.101762, 0.104887, 0.106119, 0.104617, 0.100805};
170  for (int i = 0; i < 26; ++i) { m_emMopFraction[i] = emMopFraction[i]; }
171 
172  // Initialization of the dead material map arrays
173  constexpr double dead_eta[277] = {
174  0.005, 0.015, 0.025, 0.035, 0.045, 0.055, 0.065, 0.075, 0.085, 0.095, 0.105, 0.115, 0.125, 0.135, 0.145, 0.155, 0.165, 0.175,
175  0.185, 0.195, 0.205, 0.215, 0.225, 0.235, 0.245, 0.255, 0.265, 0.275, 0.285, 0.295, 0.305, 0.315, 0.325, 0.335, 0.345, 0.355,
176  0.365, 0.375, 0.385, 0.395, 0.405, 0.415, 0.425, 0.435, 0.445, 0.455, 0.465, 0.475, 0.485, 0.495, 0.505, 0.515, 0.525, 0.535,
177  0.545, 0.555, 0.565, 0.575, 0.585, 0.595, 0.605, 0.615, 0.625, 0.635, 0.645, 0.655, 0.665, 0.675, 0.685, 0.695, 0.705, 0.715,
178  0.725, 0.735, 0.745, 0.755, 0.765, 0.775, 0.785, 0.795, 0.805, 0.815, 0.825, 0.835, 0.845, 0.855, 0.865, 0.875, 0.885, 0.895,
179  0.905, 0.915, 0.925, 0.935, 0.945, 0.955, 0.965, 0.975, 0.985, 0.995, 1.005, 1.015, 1.025, 1.035, 1.045, 1.055, 1.065, 1.075,
180  1.085, 1.095, 1.105, 1.115, 1.125, 1.135, 1.145, 1.155, 1.165, 1.175, 1.185, 1.195, 1.205, 1.215, 1.225, 1.235, 1.245, 1.255,
181  1.265, 1.275, 1.285, 1.295, 1.305, 1.315, 1.325, 1.335, 1.345, 1.355, 1.365, 1.375, 1.385, 1.395, 1.405, 1.415, 1.425, 1.435,
182  1.445, 1.455, 1.465, 1.475, 1.485, 1.495, 1.505, 1.515, 1.525, 1.535, 1.545, 1.555, 1.565, 1.575, 1.585, 1.595, 1.605, 1.615,
183  1.625, 1.635, 1.645, 1.655, 1.665, 1.675, 1.685, 1.695, 1.705, 1.715, 1.725, 1.735, 1.745, 1.755, 1.765, 1.775, 1.785, 1.795,
184  1.805, 1.815, 1.825, 1.835, 1.845, 1.855, 1.865, 1.875, 1.885, 1.895, 1.905, 1.915, 1.925, 1.935, 1.945, 1.955, 1.965, 1.975,
185  1.985, 1.995, 2.005, 2.015, 2.025, 2.035, 2.045, 2.055, 2.065, 2.075, 2.085, 2.095, 2.105, 2.115, 2.125, 2.135, 2.145, 2.155,
186  2.165, 2.175, 2.185, 2.195, 2.205, 2.215, 2.225, 2.235, 2.245, 2.255, 2.265, 2.275, 2.285, 2.295, 2.305, 2.315, 2.325, 2.335,
187  2.345, 2.355, 2.365, 2.375, 2.385, 2.395, 2.405, 2.415, 2.425, 2.435, 2.445, 2.455, 2.465, 2.475, 2.485, 2.495, 2.505, 2.515,
188  2.525, 2.535, 2.545, 2.555, 2.565, 2.575, 2.585, 2.595, 2.605, 2.615, 2.625, 2.635, 2.645, 2.655, 2.665, 2.675, 2.685, 2.695,
189  2.705, 2.715, 2.725, 2.735, 2.745, 2.755, 2.765};
190  constexpr double dead_x01[277] = {
191  1.8469, 1.7213, 1.7217, 1.7222, 1.7366, 1.7402, 1.7667, 1.7847, 1.7927, 1.7968, 1.8054, 1.7959, 1.7737, 1.7626, 1.7475, 1.7259,
192  1.7211, 1.7032, 1.6776, 1.6731, 1.6688, 1.6516, 1.6474, 1.6497, 1.6394, 1.6384, 1.6322, 1.6288, 1.6187, 1.6251, 1.6474, 1.6547,
193  1.6712, 1.6758, 1.6901, 1.6935, 1.6883, 1.6966, 1.7005, 1.6978, 1.7044, 1.7184, 1.7257, 1.7322, 1.7327, 1.7402, 1.7483, 1.7332,
194  1.7477, 1.7599, 1.783, 1.7907, 1.8299, 1.8343, 1.8427, 1.8513, 1.8271, 1.8276, 1.8345, 1.8702, 1.8821, 1.8889, 1.9048, 1.9538,
195  2.0558, 2.1443, 2.2201, 2.2642, 2.2656, 2.2472, 2.2776, 2.2798, 2.3388, 2.3809, 2.3763, 2.3869, 2.4254, 2.4758, 2.5282, 2.5364,
196  2.5637, 2.5917, 2.5999, 2.6376, 2.6386, 2.6423, 2.6479, 2.674, 2.676, 2.6806, 2.7002, 2.6985, 2.704, 2.6898, 2.7206, 2.7414,
197  2.7647, 2.7856, 2.7867, 2.8078, 2.8187, 2.8461, 2.8572, 2.8544, 2.8487, 2.8796, 2.9303, 2.9401, 2.928, 2.9376, 2.9437, 3.0083,
198  2.9823, 2.9945, 2.9809, 3.0485, 3.0494, 3.0595, 3.0642, 3.1098, 3.1232, 3.1306, 3.0995, 3.1497, 3.1399, 3.1285, 3.1489, 3.1748,
199  3.2728, 3.2219, 3.2252, 3.2826, 3.3505, 3.4002, 3.4492, 3.4564, 3.4622, 3.4631, 3.4628, 3.5027, 3.5506, 3.6616, 3.8859, 4.0115,
200  4.1715, 4.3591, 4.4486, 4.4909, 4.5445, 4.5871, 4.6396, 4.398, 3.6563, 4.672, 3.9715, 3.0968, 2.7105, 2.7465, 2.7719, 2.7361,
201  2.8049, 2.7947, 2.7961, 2.8328, 2.8847, 2.8935, 3.1969, 3.5616, 3.9473, 3.3614, 2.5151, 2.3433, 2.3557, 2.0129, 2.0082, 2.0186,
202  1.9549, 1.9559, 1.873, 1.9516, 1.9924, 2.0921, 2.1637, 2.0858, 2.0419, 1.9849, 2.0455, 1.9273, 1.8368, 1.7655, 1.779, 1.726,
203  1.724, 1.7487, 1.7251, 1.6957, 1.7127, 1.6828, 1.6412, 1.6795, 1.6421, 1.631, 1.6366, 1.6835, 1.6884, 1.6811, 1.6724, 1.6534,
204  1.6404, 1.6498, 1.6268, 1.6205, 1.6453, 1.657, 1.6372, 1.6172, 1.5865, 1.561, 1.5822, 1.6654, 1.6161, 1.57, 1.5173, 1.4629,
205  1.5068, 1.4886, 1.496, 1.5121, 1.4521, 1.4178, 1.4192, 1.4095, 1.42, 1.4151, 1.4129, 1.3897, 1.4006, 1.3767, 1.3668, 1.352,
206  1.3715, 1.3606, 1.3526, 1.3903, 1.4008, 1.3913, 1.363, 1.3655, 1.3634, 1.366, 1.3687, 1.3626, 1.3925, 1.3764, 1.3442, 1.3366,
207  1.3301, 1.3277, 1.3036, 1.304, 1.3068, 1.3114, 1.303, 1.345, 1.3567, 1.3424, 1.3092, 1.3122, 1.3137, 1.3231, 1.3236, 1.2951,
208  1.2928, 1.3456, 1.8221, 2.0029, 2.0075};
209  constexpr double dead_x02[277] = {
210  27.3438, 27.2208, 27.2263, 27.243, 31.8853, 28.1255, 27.5756, 27.0132, 27.0413, 27.0637, 27.1041, 27.1223, 27.1304, 27.1563,
211  27.1702, 27.1865, 27.2222, 27.8551, 28.488, 30.206, 28.5837, 28.0049, 27.4459, 27.5156, 27.5572, 27.6216, 27.6856, 27.7451,
212  27.8147, 27.8949, 27.9914, 28.078, 28.1787, 28.2678, 28.3671, 28.4603, 29.1869, 29.9369, 31.3911, 29.6731, 29.6085, 29.0849,
213  29.2069, 29.3229, 29.4405, 29.5618, 29.6891, 29.7915, 29.9375, 30.0806, 30.2322, 30.3764, 30.5561, 30.6942, 30.8476, 31.6978,
214  31.8183, 33.6375, 33.149, 32.3431, 32.5023, 31.9797, 32.1619, 32.3803, 32.6547, 32.9188, 33.1731, 33.3989, 33.5851, 33.7547,
215  33.9763, 34.1728, 34.4293, 34.676, 35.3402, 35.8687, 36.8943, 37.3228, 37.0672, 35.8788, 30.0848, 30.2997, 30.4957, 30.7282,
216  30.917, 31.1227, 31.3252, 31.5502, 31.7627, 31.9752, 32.2057, 32.4178, 32.6401, 32.8458, 33.9831, 34.4918, 34.9247, 35.4172,
217  35.0033, 35.2645, 34.6033, 34.8723, 35.1359, 35.8106, 37.0775, 37.4231, 37.5118, 37.6158, 38.0132, 38.5355, 38.9488, 40.1903,
218  40.5554, 39.3216, 40.47, 40.5203, 40.3146, 40.0119, 40.0371, 40.0744, 40.0102, 40.0563, 39.9966, 40.0213, 40.0384, 39.9564,
219  38.3485, 37.0049, 36.011, 34.3675, 32.746, 31.6373, 30.2222, 28.8545, 31.2748, 30.5887, 29.1005, 28.1896, 28.9692, 29.0944,
220  29.5132, 29.8182, 30.2405, 30.4836, 30.6649, 31.4598, 31.8968, 30.9419, 32.2847, 32.2785, 29.4552, 28.8184, 27.8418, 26.9469,
221  26.5213, 26.0314, 26.0441, 26.4897, 26.9313, 27.3156, 27.7043, 27.9047, 28.1472, 28.4043, 28.6613, 28.8634, 29.3494, 29.8874,
222  30.4377, 30.0112, 29.3163, 29.2913, 29.4448, 29.2385, 29.3655, 29.5045, 29.5655, 29.6873, 29.7209, 29.9143, 30.0671, 30.2626,
223  30.4579, 30.4513, 30.5094, 30.5531, 30.7121, 30.6904, 30.6949, 30.7176, 30.8232, 30.9326, 31.0953, 31.376, 31.6177, 31.8488,
224  32.1461, 32.4569, 32.7735, 33.1196, 33.4516, 33.7533, 34.0885, 34.4624, 34.7467, 34.9918, 35.2164, 35.4157, 35.6083, 35.812,
225  35.9719, 36.1537, 36.3483, 36.5131, 36.65, 36.7811, 36.8951, 37.0119, 37.1696, 37.3852, 37.4656, 37.5443, 37.6143, 37.6783,
226  37.839, 37.9339, 38.0519, 38.1968, 38.2431, 38.3125, 38.4156, 38.5057, 38.6145, 38.7049, 38.797, 38.8656, 38.9668, 39.031,
227  39.108, 39.177, 39.2796, 39.2613, 39.3312, 39.4452, 39.6109, 39.6561, 39.7265, 39.7582, 39.7801, 39.8325, 40.0373, 30.445,
228  30.4259, 30.3236, 30.335, 30.4099, 30.4834, 30.4998, 30.549, 30.6477, 30.6725, 30.7221, 30.7792, 30.8878, 30.9677, 31.025,
229  31.0644, 31.1425, 31.2212, 31.3097, 31.392, 31.4477, 31.5309, 31.6727, 32.2409, 32.5151, 32.6164};
230  constexpr double dead_x0tot[277] = {
231  108.064, 114.015, 102.328, 103.979, 108.575, 99.904, 107.609, 114.499, 110.58, 105.503, 111.6, 109.618, 107.369, 108.663,
232  108.656, 108.724, 107.947, 111.123, 113.801, 114.927, 110.352, 113.048, 109.105, 112.754, 110.352, 112.808, 113.401, 111.948,
233  113.811, 111.447, 114.069, 115.163, 115.728, 113.215, 116.055, 115.433, 116.233, 118.157, 119.787, 118.035, 117.5, 117.333,
234  116.777, 119.61, 119.342, 121.258, 119.331, 121.887, 122.346, 122.314, 122.692, 122.296, 123.452, 124.217, 125.029, 126.766,
235  127.559, 129.726, 129.524, 130.828, 128.679, 129.477, 129.923, 131.236, 132.793, 132.212, 132.597, 133.999, 131.29, 129.22,
236  126.212, 126.427, 127.45, 124.279, 119.29, 115.547, 114.276, 116.323, 117.272, 117.058, 112.01, 112.727, 114.061, 115.384,
237  115.753, 114.081, 112.184, 111.034, 109.526, 109.854, 111.763, 113.335, 114.989, 116.799, 119.382, 121.044, 122.683, 123.843,
238  123.798, 124.567, 122.571, 122.951, 124.179, 124.977, 128.058, 131.607, 135.404, 139.132, 142.167, 146.221, 150.063, 154.225,
239  158.159, 160.614, 164.357, 168.218, 171.728, 174.947, 177.505, 181.286, 183.258, 186.901, 192.562, 189.001, 189.092, 185.551,
240  181.202, 177.128, 168.268, 159.719, 156.125, 151.878, 149.137, 144.408, 141.419, 138.169, 133.679, 130.259, 128.061, 126.093,
241  123.544, 121.135, 119.087, 116.403, 113.911, 112.066, 109.821, 105.606, 104.274, 101.77, 100.493, 99.9533, 98.9401, 98.3901,
242  97.6717, 97.3916, 97.613, 98.0975, 98.6056, 99.3877, 99.9002, 100.231, 100.728, 101.391, 102.93, 105.539, 106.374, 107.412,
243  108.678, 108.733, 111.575, 110.964, 113.989, 117.078, 120.737, 124.857, 129.56, 134.323, 138.298, 142.47, 146.55, 149.97,
244  153.422, 153.182, 152.819, 152.886, 152.96, 152.822, 152.713, 152.624, 152.621, 152.623, 152.681, 152.86, 153.001, 153.134,
245  153.335, 153.551, 153.775, 154.031, 154.275, 154.49, 154.74, 155.03, 155.233, 155.398, 155.545, 155.667, 155.785, 155.915,
246  156.003, 156.114, 156.24, 156.337, 156.407, 156.474, 156.524, 162.991, 163.283, 163.632, 163.845, 164.055, 164.256, 164.851,
247  165.542, 165.713, 165.796, 165.9, 166.052, 166.246, 166.473, 166.686, 166.917, 167.129, 167.341, 167.53, 167.75, 167.932,
248  168.201, 168.312, 168.131, 168.228, 168.412, 168.64, 168.918, 169.075, 169.257, 169.399, 169.53, 169.691, 170.004, 160.519,
249  160.606, 160.61, 160.726, 160.905, 161.082, 161.201, 161.352, 161.552, 161.677, 161.826, 161.982, 162.189, 162.366, 162.52,
250  162.655, 162.829, 163.002, 163.184, 163.36, 163.508, 163.683, 163.916, 164.574, 164.938, 165.128};
251  constexpr double dead_cryo0[277] = {
252  0.6732, 0.6734, 0.6735, 0.6737, 0.6739, 0.6743, 0.6746, 0.6751, 0.6757, 0.6763, 0.677, 0.6777, 0.6785, 0.6794, 0.6804, 0.6813,
253  0.6825, 0.6836, 0.6848, 0.6861, 0.6874, 0.6889, 0.6904, 0.6919, 0.6936, 0.6952, 0.697, 0.6988, 0.7007, 0.7028, 0.7048, 0.707,
254  0.7091, 0.7114, 0.7137, 0.7162, 0.7186, 0.7212, 0.7237, 0.7265, 0.7293, 0.7321, 0.735, 0.7379, 0.741, 0.7442, 0.7474, 0.7507,
255  0.754, 0.7575, 0.7609, 0.7646, 0.7682, 0.7719, 0.7758, 0.7796, 0.7836, 0.7876, 0.7918, 0.796, 0.8002, 0.8046, 0.8091, 0.8136,
256  0.8182, 0.8229, 0.8277, 0.8326, 0.8375, 0.8425, 0.8476, 0.8528, 0.8581, 0.8634, 0.869, 0.8744, 0.88, 0.8857, 0.8916, 0.8975,
257  0.9034, 0.9095, 0.9157, 0.9219, 0.9283, 0.9347, 0.9412, 0.9479, 0.9546, 0.9614, 0.9683, 0.9753, 0.9824, 0.9896, 0.997, 1.0043,
258  1.0119, 1.0194, 1.0271, 1.0349, 1.0429, 1.0509, 1.059, 1.0672, 1.0756, 1.084, 1.0925, 1.1012, 1.1099, 1.1189, 1.1278, 1.137,
259  1.1461, 1.1556, 1.1649, 1.1746, 1.1842, 1.1941, 1.2039, 1.214, 1.2242, 1.2344, 1.2449, 1.2553, 1.266, 1.2768, 1.2877, 1.2988,
260  1.31, 1.3212, 1.3327, 1.3443, 1.3559, 1.3678, 1.3798, 1.3919, 1.4041, 1.4165, 1.429, 1.4418, 1.4545, 1.4675, 1.4807, 1.4938,
261  1.5073, 1.5208, 1.5346, 1.5485, 1.5625, 1.5766, 1.5601, 1.2354, 0.4486, 0.5128, 0.3507, 0.1844, 0.1841, 0.1837, 0.1835, 0.1832,
262  0.1828, 0.1826, 0.1823, 0.182, 0.1817, 0.1814, 0.1812, 0.181, 0.1807, 0.1804, 0.1802, 0.1799, 0.1797, 0.1795, 0.1793, 0.1791,
263  0.1789, 0.1787, 0.1784, 0.1782, 0.178, 0.1778, 0.1777, 0.1775, 0.1773, 0.1771, 0.177, 0.1768, 0.1766, 0.1764, 0.1763, 0.1762,
264  0.176, 0.1759, 0.1758, 0.1756, 0.1754, 0.1753, 0.1752, 0.1751, 0.1749, 0.1748, 0.1746, 0.1745, 0.1745, 0.1743, 0.1742, 0.174,
265  0.174, 0.1739, 0.1737, 0.1736, 0.1735, 0.1735, 0.1734, 0.1732, 0.1731, 0.1731, 0.173, 0.1729, 0.1728, 0.1728, 0.1727, 0.1726,
266  0.1725, 0.1725, 0.1724, 0.1722, 0.1722, 0.1721, 0.172, 0.172, 0.1719, 0.1719, 0.1718, 0.1717, 0.1716, 0.1716, 0.1716, 0.1715,
267  0.1714, 0.1714, 0.1713, 0.1713, 0.1712, 0.1711, 0.1711, 0.1711, 0.1711, 0.171, 0.1709, 0.1709, 0.1709, 0.1708, 0.1708, 0.1708,
268  0.1707, 0.1706, 0.1706, 0.1706, 0.1706, 0.1705, 0.1704, 0.1704, 0.1704, 0.1703, 0.1703, 0.1703, 0.1702, 0.1702, 0.1702, 0.1702,
269  0.1701, 0.1701, 0.1701, 0.17, 0.17};
270  constexpr double dead_cryo1[277] = {
271  1.6192, 1.6194, 1.6197, 1.6203, 6.3509, 2.5841, 1.9874, 1.3905, 1.3916, 1.3929, 1.3943, 1.3957, 1.3975, 1.3992, 1.4013, 1.4033,
272  1.4055, 2.0138, 2.6243, 4.3032, 2.6344, 2.0294, 1.4221, 1.4251, 1.4284, 1.432, 1.4357, 1.4394, 1.4433, 1.4475, 1.4516, 1.4561,
273  1.4605, 1.4652, 1.4702, 1.475, 2.117, 2.7636, 4.1287, 2.3159, 2.1483, 1.5079, 1.5139, 1.5199, 1.5261, 1.5326, 1.5393, 1.546,
274  1.553, 1.5599, 1.5673, 1.5747, 1.5821, 1.5899, 1.5977, 2.297, 2.3084, 3.9766, 3.3404, 2.3449, 2.345, 1.6573, 1.6664, 1.6757,
275  1.6853, 1.6949, 1.7047, 1.7148, 1.7249, 1.7352, 1.7457, 1.7564, 1.7673, 1.7784, 2.2539, 2.576, 3.3648, 3.5384, 3.0253, 2.6439,
276  1.8606, 1.8733, 1.8859, 1.8989, 1.9119, 1.9251, 1.9387, 1.9521, 1.966, 1.98, 1.9943, 2.0088, 2.0235, 2.0383, 2.937, 3.214,
277  3.4082, 3.6655, 3.026, 3.0489, 2.148, 2.1643, 2.1812, 2.198, 3.6572, 3.7318, 3.5266, 3.3742, 3.534, 3.7906, 3.7431, 4.3293,
278  4.4603, 2.9506, 3.8466, 3.5582, 3.078, 2.4882, 2.2279, 1.7277, 1.2731, 1.0272, 0.7105, 0.3933, 0.2912, 0.2902, 0.2893, 0.2883,
279  0.2875, 0.2866, 0.2857, 0.2849, 0.284, 0.2834, 2.6585, 2.6785, 2.6989, 2.7197, 2.7405, 2.7614, 2.7828, 2.8045, 2.8265, 2.8488,
280  2.8712, 2.8939, 2.9114, 2.7992, 2.9324, 3.4047, 0.414, 0.252, 0.1168, 0.9172, 0.9156, 0.914, 0.9123, 0.9108, 0.9091, 0.9076,
281  0.9062, 0.9048, 0.9033, 0.9018, 0.9006, 0.8993, 0.8979, 0.8966, 0.8955, 0.8942, 0.893, 0.8918, 0.8907, 0.8896, 0.8886, 0.8874,
282  0.8865, 0.8853, 0.8843, 0.8833, 0.8823, 0.9288, 0.9279, 0.927, 0.926, 0.9251, 0.9242, 0.9233, 0.9224, 0.9217, 0.9208, 0.9199,
283  0.9192, 0.9184, 0.9176, 0.9169, 0.9162, 0.9154, 0.9148, 0.9141, 0.9135, 0.9128, 0.9121, 0.9115, 0.9109, 0.9103, 0.9097, 0.9091,
284  0.9086, 0.9079, 0.9074, 0.9069, 0.9064, 0.9059, 0.9053, 0.9048, 0.9043, 0.9038, 0.9034, 0.903, 0.9025, 0.9022, 0.9017, 0.9013,
285  0.9009, 0.9004, 0.9002, 0.8996, 0.8994, 0.899, 0.8986, 0.8982, 0.8979, 0.8976, 0.8971, 0.897, 0.8966, 0.8963, 0.896, 0.8957,
286  0.8954, 0.8951, 0.8948, 0.8945, 0.8942, 0.8941, 0.8937, 0.8934, 0.8932, 0.8929, 0.8927, 0.8924, 0.8922, 0.892, 0.8917, 0.8915,
287  0.8914, 0.8912, 0.891, 0.8907, 0.8905, 0.8903, 0.8902, 0.89, 0.8898, 0.8897, 0.8894, 0.8893, 0.889, 0.8889, 0.8888, 0.8886,
288  0.8886, 0.8883, 0.8882, 0.888, 0.8879};
289  constexpr double dead_cryo2[277] = {
290  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
291  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
292  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
293  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
294  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
295  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2895, 0.9101, 1.5272, 1.8249, 1.8279, 1.8312,
296  1.8348, 1.8385, 1.8424, 1.8466, 1.851, 1.8555, 1.8604, 1.5309, 0, 0, 0, 0, 0, 0, 0, 0,
297  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0763, 0.6379, 0.8382, 1.0003, 1.2294,
298  1.5831, 2.0439, 2.3796, 2.6456, 2.9827, 3.3447, 1.141, 1.138, 1.1349, 1.132, 1.1291, 1.1264, 1.1236, 1.1208, 1.1183, 1.1156,
299  1.1132, 1.1108, 1.1084, 1.106, 1.1037, 1.1014, 4.514, 4.4934, 4.5226, 3.5091, 3.5412, 3.5737, 3.6066, 3.6398, 3.6734, 3.7074,
300  3.7418, 3.7764, 3.8115, 4.423, 5.2744, 6.1316, 6.6266, 6.6847, 6.7435, 6.803, 9.6759, 9.0086, 8.1971, 7.7761, 7.5583, 7.7543,
301  8.5776, 9.3644, 9.4149, 9.466, 9.4281, 8.6493, 7.8651, 7.5418, 7.2376, 7.3746, 7.4009, 7.3938, 7.3869, 7.3801, 7.3735, 7.3671,
302  7.3608, 7.3546, 7.3485, 7.3426, 7.3367, 7.331, 7.3254, 7.3199, 7.3146, 7.3093, 7.3042, 7.2991, 7.2943, 7.2894, 7.2846, 7.28,
303  7.2754, 7.271, 7.2666, 7.2624, 7.2582, 7.2541, 7.2501, 7.2461, 7.2423, 7.2386, 7.2348, 7.2312, 7.2276, 7.2242, 7.2207, 7.2174,
304  7.2141, 7.211, 7.2078, 7.2047, 7.2017, 7.1987, 7.1958, 7.1931, 7.1902, 7.1875, 7.1848, 7.1822, 7.1796, 7.1772, 7.1747, 7.1723,
305  7.1699, 7.1676, 7.1653, 7.1632, 7.161, 7.1588, 7.1567, 7.1547, 7.1527, 7.1507, 7.1487, 7.1469, 7.145, 7.1433, 7.1415, 7.1397,
306  7.138, 7.1364, 7.1347, 7.1331, 7.1315, 7.13, 7.1284, 7.1269, 7.1255, 7.1241, 7.1227, 7.1213, 7.12, 7.1187, 7.1175, 7.1162,
307  7.1149, 7.1137, 7.1126, 7.1114, 7.1102};
308  constexpr double dead_girder[277] = {
309  11.9514, 11.9525, 11.9549, 11.9584, 11.9632, 11.9693, 11.9764, 11.9848, 11.9943, 12.0051, 12.0171, 12.0303, 12.0446, 12.0602,
310  12.077, 12.095, 12.1142, 12.1346, 12.1563, 12.1791, 12.2032, 12.2285, 12.255, 12.2827, 12.3116, 12.3418, 12.3732, 12.4059,
311  12.4398, 12.475, 12.5113, 12.549, 12.5879, 12.6281, 12.6695, 12.7122, 12.7561, 12.8014, 12.8479, 12.8957, 12.9448, 12.9952,
312  13.0469, 13.0998, 13.1541, 13.2098, 13.2667, 13.325, 13.3845, 13.4455, 13.5078, 13.5714, 13.6364, 13.7028, 13.7705, 13.8396,
313  13.91, 13.9818, 14.0552, 14.1298, 14.206, 14.2835, 12.3788, 7.5701, 7.4957, 6.6037, 3.8717, 0, 0, 0,
314  0, 3.5585, 7.2483, 7.505, 8.4513, 10.9561, 15.539, 15.7235, 15.8265, 15.9311, 16.0372, 16.1449, 16.2543, 16.3653,
315  16.4779, 16.5921, 16.7081, 16.8256, 16.9449, 17.0659, 17.1886, 17.313, 17.4392, 17.567, 17.6966, 17.8281, 17.9612, 18.0962,
316  18.233, 18.3716, 18.5121, 18.6543, 18.7985, 18.9446, 19.0925, 19.2424, 19.3942, 19.5479, 19.7036, 19.8612, 20.0208, 20.1824,
317  20.3461, 20.5118, 20.6795, 20.8493, 21.0212, 16.6712, 11.697, 11.2949, 10.114, 9.906, 8.5687, 3.7101, 0, 0,
318  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
319  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
320  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
321  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
322  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
323  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
324  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
325  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
326  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
327  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
328  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
329  for (int i = 0; i < 277; i++) {
330  m_dead_eta[i] = dead_eta[i];
331  m_dead_x01[i] = dead_x01[i];
332  m_dead_x02[i] = dead_x02[i];
333  m_dead_x0tot[i] = dead_x0tot[i];
334  m_dead_cryo0[i] = dead_cryo0[i];
335  m_dead_cryo1[i] = dead_cryo1[i];
336  m_dead_cryo2[i] = dead_cryo2[i];
337  m_dead_girder[i] = dead_girder[i];
338  }
339  m_etaGranularity = 0.01;
340  //---- End of material map arrays initialization -----//
341 
342  return StatusCode::SUCCESS;
343  }

◆ 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& Rec::IMuidCaloEnergyParam::interfaceID ( )
inlinestaticinherited

AlgTool and IAlgTool interface methods.

Definition at line 37 of file IMuidCaloEnergyParam.h.

37 { return IID_IMuidCaloEnergyParam; }

◆ meanEnergyLoss()

double MuidCaloEnergyParam::meanEnergyLoss ( double  eta,
double  momentum 
) const
private

Definition at line 513 of file MuidCaloEnergyParam.cxx.

513  {
514  // Mean Energy Loss parametrization with asymmetric gaussian fit
515  // tabulated coefficients for function [0]+[1]*Log(0.0067*p*p)+[2]*p with GeV units
516  MuidCaloEnergyParam::BinsWeights BW = etaVariableBin(eta);
517  double p0 = BW.etaWeight1 * m_meanEnergyLossP0[BW.etaBin1] + BW.etaWeight2 * m_meanEnergyLossP0[BW.etaBin2];
518  double p1 = BW.etaWeight1 * m_meanEnergyLossP1[BW.etaBin1] + BW.etaWeight2 * m_meanEnergyLossP1[BW.etaBin2];
519  double p2 = BW.etaWeight1 * m_meanEnergyLossP2[BW.etaBin1] + BW.etaWeight2 * m_meanEnergyLossP2[BW.etaBin2];
520  double parametrisedDeposit = p0 * Units::GeV + p1 * Units::GeV * std::log(0.0067 * std::pow(momentum * PerGeV, 2)) + p2 * momentum;
521 
522  // additional offset from high-statistics Z->mumu MC (measured by Peter K 30/11/2011, 03/2012)
523  // NOTE: constant bin width used here
524  BW = etaFixedBin(eta);
525  // mc11 tune
526  // double fixFromPeter[26] = { -0.203781 , -0.146775 , -0.0952301 , -0.186141 , -0.180687 ,
527  // -0.225095 , -0.195912 , -0.221158 , -0.198556 , -0.271551 ,
528  // -0.200737 , -0.151854 , -0.0799826 , -0.237688 , -0.255126 ,
529  // -0.313466 , -0.35815 , -0.334661 , -0.458948 , -0.469967 ,
530  // -0.516059 , -0.522978 , -0.53896 , -0.550683 , -0.51722 ,
531  // -0.543858 };
532  // mc12 tune
533  constexpr double fixFromPeter[26] = {-0.159195, -0.039180, 0.071644, -0.118753, -0.144227, -0.182778, -0.158818,
534  -0.188853, -0.180997, -0.269571, -0.227767, -0.191740, -0.121575, -0.304938,
535  -0.296636, -0.312171, -0.352982, -0.303020, -0.357250, -0.396249, -0.433535,
536  -0.417214, -0.424093, -0.415557, -0.325372, -0.553009};
537  double fix = BW.etaWeight1 * fixFromPeter[BW.etaBin1] + BW.etaWeight2 * fixFromPeter[BW.etaBin2];
538  parametrisedDeposit += fix * Units::GeV;
539 
540  return parametrisedDeposit;
541  }

◆ meanEnergyLossError()

std::pair< double, double > MuidCaloEnergyParam::meanEnergyLossError ( double  eta,
double  momentum 
) const
private

Definition at line 543 of file MuidCaloEnergyParam.cxx.

543  {
544  MuidCaloEnergyParam::BinsWeights BW = etaVariableBin(eta);
545  double left_p0 = BW.etaWeight1 * m_meanEnergyLossErrorLeftP0[BW.etaBin1] + BW.etaWeight2 * m_meanEnergyLossErrorLeftP0[BW.etaBin2];
546  double left_p1 = BW.etaWeight1 * m_meanEnergyLossErrorLeftP1[BW.etaBin1] + BW.etaWeight2 * m_meanEnergyLossErrorLeftP1[BW.etaBin2];
547  double right_p0 =
548  BW.etaWeight1 * m_meanEnergyLossErrorRightP0[BW.etaBin1] + BW.etaWeight2 * m_meanEnergyLossErrorRightP0[BW.etaBin2];
549  double right_p1 =
550  BW.etaWeight1 * m_meanEnergyLossErrorRightP1[BW.etaBin1] + BW.etaWeight2 * m_meanEnergyLossErrorRightP1[BW.etaBin2];
551  double sigma_left = left_p0 * Units::GeV + left_p1 * momentum;
552  double sigma_right = right_p0 * Units::GeV + right_p1 * momentum;
553 
554  return std::make_pair(sigma_left, sigma_right);
555  }

◆ meanParametrizedEnergy()

CaloEnergy * MuidCaloEnergyParam::meanParametrizedEnergy ( double  trackMomentum,
double  eta,
double  phi 
) const
virtual

IMuidCaloEnergyParam interface: to get the mean muon energy loss in the calorimeter (with asymmetric errors describing the whole distribution)

Implements Rec::IMuidCaloEnergyParam.

Definition at line 344 of file MuidCaloEnergyParam.cxx.

344  {
345  CaloEnergy* caloEnergy = nullptr;
346  double meanEloss = meanEnergyLoss(eta, trackMomentum);
347  std::pair<double, double> parametrizedErrors = meanEnergyLossError(eta, trackMomentum);
348 
349  ATH_MSG_DEBUG(" Mean_deposited_energy = " << meanEloss / Units::GeV);
350 
351  // cosmics: energy sign flip in upper hemisphere
352  if (m_cosmics && phi > 0.) {
353  caloEnergy = new CaloEnergy(-meanEloss, 0.5 * (parametrizedErrors.first + parametrizedErrors.second), parametrizedErrors.second,
354  parametrizedErrors.first, CaloEnergy::Parametrized);
355  } else {
356  caloEnergy = new CaloEnergy(meanEloss, 0.5 * (parametrizedErrors.first + parametrizedErrors.second), parametrizedErrors.first,
357  parametrizedErrors.second, CaloEnergy::Parametrized);
358  }
359  // set parametrized energy loss in em
360  double emLoss = caloEnergy->deltaE() * emMopFraction(eta);
361  caloEnergy->set_fsrCandidateEnergy(emLoss);
362  return caloEnergy;
363  }

◆ mopDepositedEnergy()

CaloEnergy * MuidCaloEnergyParam::mopDepositedEnergy ( double  trackMomentum,
double  eta,
double  phi 
) const
virtual

IMuidCaloEnergyParam interface: to get the most probable (peak value) energy loss and error in the calorimeter.

Implements Rec::IMuidCaloEnergyParam.

Definition at line 434 of file MuidCaloEnergyParam.cxx.

434  {
435  CaloEnergy* caloEnergy = nullptr;
436 
437  std::pair<double, double> eta_phi;
438  eta_phi.first = eta;
439  eta_phi.second = phi;
440 
441  double mop_eloss = mopEnergyLoss(eta_phi.first, trackMomentum);
442 
443  // double mop_eloss = symmetricMopEnergyLoss(eta_phi.first,trackMomentum);
444  double mop_eloss_error = mopEnergyLossError(eta_phi.first, trackMomentum);
445 
446  ATH_MSG_DEBUG(" Most probable deposited energy... " << mop_eloss / Units::GeV);
447 
448  // cosmics: energy sign flip in upper hemisphere
449  if (m_cosmics && phi > 0.) {
450  caloEnergy = new CaloEnergy(-mop_eloss, mop_eloss_error, mop_eloss_error, mop_eloss_error, CaloEnergy::Parametrized);
451  } else {
452  caloEnergy = new CaloEnergy(mop_eloss, mop_eloss_error, mop_eloss_error, mop_eloss_error, CaloEnergy::Parametrized);
453  }
454  // set parametrized energy loss in em
455  double emLoss = caloEnergy->deltaE() * emMopFraction(eta);
456  caloEnergy->set_fsrCandidateEnergy(emLoss);
457  return caloEnergy;
458  }

◆ mopEnergyLoss()

double MuidCaloEnergyParam::mopEnergyLoss ( double  eta,
double  momentum 
) const
private

Definition at line 557 of file MuidCaloEnergyParam.cxx.

557  {
558  // Mop Energy Loss parametrization with asymmetric gaussian fit in the mop region
559  // tabulated coefficients for function [0]+[1]*Log(0.0067*p*p)+[2]*p with GeV units
560  MuidCaloEnergyParam::BinsWeights BW = etaVariableBin(eta);
561  double p0 = BW.etaWeight1 * m_mopEnergyLossP0[BW.etaBin1] + BW.etaWeight2 * m_mopEnergyLossP0[BW.etaBin2];
562  double p1 = BW.etaWeight1 * m_mopEnergyLossP1[BW.etaBin1] + BW.etaWeight2 * m_mopEnergyLossP1[BW.etaBin2];
563  double p2 = BW.etaWeight1 * m_mopEnergyLossP2[BW.etaBin1] + BW.etaWeight2 * m_mopEnergyLossP2[BW.etaBin2];
564  double parametrisedDeposit = p0 * Units::GeV + p1 * Units::GeV * std::log(0.0067 * std::pow(momentum * PerGeV, 2)) + p2 * momentum;
565 
566  // additional offset from high-statistics Z->mumu MC (measured by Peter K 30/11/2011)
567  // NOTE: constant bin width used here
568  BW = etaFixedBin(eta);
569  // mc11 tune
570  // double fixFromPeter[26] = { -0.203781 , -0.146775 , -0.0952301 , -0.186141 , -0.180687 ,
571  // -0.225095 , -0.195912 , -0.221158 , -0.198556 , -0.271551 ,
572  // -0.200737 , -0.151854 , -0.0799826 , -0.237688 , -0.255126 ,
573  // -0.313466 , -0.35815 , -0.334661 , -0.458948 , -0.469967 ,
574  // -0.516059 , -0.522978 , -0.53896 , -0.550683 , -0.51722 ,
575  // -0.543858 };
576  // mc12 tune
577  constexpr double fixFromPeter[26] = {-0.159195, -0.039180, 0.071644, -0.118753, -0.144227, -0.182778, -0.158818,
578  -0.188853, -0.180997, -0.269571, -0.227767, -0.191740, -0.121575, -0.304938,
579  -0.296636, -0.312171, -0.352982, -0.303020, -0.357250, -0.396249, -0.433535,
580  -0.417214, -0.424093, -0.415557, -0.325372, -0.553009};
581 
582  double fix = BW.etaWeight1 * fixFromPeter[BW.etaBin1] + BW.etaWeight2 * fixFromPeter[BW.etaBin2];
583  parametrisedDeposit += fix * Units::GeV;
584 
585  return parametrisedDeposit;
586  }

◆ mopEnergyLossError()

double MuidCaloEnergyParam::mopEnergyLossError ( double  eta,
double  momentum 
) const
private

Definition at line 588 of file MuidCaloEnergyParam.cxx.

588  {
589  eta = std::abs(eta);
590  MuidCaloEnergyParam::BinsWeights BW = etaVariableBin(eta);
591  double sigma_p0 = BW.etaWeight1 * m_mopEnergyLossErrorP0[BW.etaBin1] + BW.etaWeight2 * m_mopEnergyLossErrorP0[BW.etaBin2];
592  double sigma_p1 = BW.etaWeight1 * m_mopEnergyLossErrorP1[BW.etaBin1] + BW.etaWeight2 * m_mopEnergyLossErrorP1[BW.etaBin2];
593  double sigma = sigma_p0 * Units::GeV + sigma_p1 * momentum;
594 
595  // fix for mc11 as underestimate in barrel, overestimate in endcap
596  if (eta < 1.7) {
597  sigma *= 1.1;
598  } else {
599  sigma *= 0.9;
600  }
601 
602  return sigma;
603  }

◆ mopParametrizedEnergy()

CaloEnergy * MuidCaloEnergyParam::mopParametrizedEnergy ( double  trackMomentum,
double  eta,
double  phi 
) const
virtual

IMuidCaloEnergyParam interface: to get the most probable (peak value) energy loss in the calorimeter (with asymmetric errors describing the whole distribution)

Implements Rec::IMuidCaloEnergyParam.

Definition at line 365 of file MuidCaloEnergyParam.cxx.

365  {
366  CaloEnergy* caloEnergy = nullptr;
367 
368  // mop and mean energy deposition
369  double mopEloss = mopEnergyLoss(eta, trackMomentum);
370  double meanEloss = meanEnergyLoss(eta, trackMomentum);
371 
372  ATH_MSG_DEBUG(" Most probable (mean) deposited energy " << mopEloss / Units::GeV << " (" << meanEloss / Units::GeV << ")");
373 
374  // mop energy deposition uncertainty (Gaussian fit around mop value +- 2 sigma)
375  double mopError = mopEnergyLossError(eta, trackMomentum);
376 
377  // scale low transverse momentum to allow for fluctuations
378  double pt = trackMomentum / std::cosh(eta);
379  double scale = 1.;
380  if (pt < FifteenGeV) scale = (2.5 - pt * (1.5 / (FifteenGeV)));
381 
382  // asymmetric error to allow for Landau
383 
384  constexpr double alpha = 0.2;
385  constexpr double beta = 300. * Units::MeV;
386  double mopErrorPlus = mopError + alpha * (meanEloss - mopEloss) + beta;
387 
388  // cosmics: energy sign flip in upper hemisphere
389  if (m_cosmics && phi > 0.) {
390  caloEnergy = new CaloEnergy(-mopEloss, scale * mopError, scale * mopErrorPlus, scale * mopError, CaloEnergy::Parametrized);
391  } else {
392  caloEnergy = new CaloEnergy(mopEloss, scale * mopError, scale * mopError, scale * mopErrorPlus, CaloEnergy::Parametrized);
393  }
394  // set parametrized energy loss in em
395  double emLoss = caloEnergy->deltaE() * emMopFraction(eta);
396  caloEnergy->set_fsrCandidateEnergy(emLoss);
397  return caloEnergy;
398  }

◆ mopPeakEnergy()

CaloEnergy * MuidCaloEnergyParam::mopPeakEnergy ( double  trackMomentum,
double  eta,
double  phi 
) const
virtual

IMuidCaloEnergyParam interface: to get the most probable (peak value) energy loss in the calorimeter (with asymmetric errors describing 2 sigma around the peak region)

Implements Rec::IMuidCaloEnergyParam.

Definition at line 400 of file MuidCaloEnergyParam.cxx.

400  {
401  CaloEnergy* caloEnergy = nullptr;
402 
403  // mop and mean energy deposition
404  double mopEloss = mopEnergyLoss(eta, trackMomentum);
405  double meanEloss = meanEnergyLoss(eta, trackMomentum);
406 
407  ATH_MSG_DEBUG(" Most probable (mean) deposited energy " << mopEloss / Units::GeV << " (" << meanEloss / Units::GeV << ")");
408 
409  // mop energy deposition uncertainty (Gaussian fit around mop value +- 2 sigma)
410  double mopError = mopEnergyLossError(eta, trackMomentum);
411 
412  // scale low transverse momentum to allow for fluctuations
413  double pt = trackMomentum / std::cosh(eta);
414  double scale = 1.;
415  if (pt < FifteenGeV) scale = (2.5 - pt * (1.5 / (FifteenGeV)));
416 
417  // asymmetric error to allow for Landau
418  constexpr double alpha = 0.1;
419  constexpr double beta = 150. * Units::MeV;
420  double mopErrorPlus = mopError + alpha * (meanEloss - mopEloss) + beta;
421 
422  // cosmics: energy sign flip in upper hemisphere
423  if (m_cosmics && phi > 0.) {
424  caloEnergy = new CaloEnergy(-mopEloss, scale * mopError, scale * mopErrorPlus, scale * mopError, CaloEnergy::MOP);
425  } else {
426  caloEnergy = new CaloEnergy(mopEloss, scale * mopError, scale * mopError, scale * mopErrorPlus, CaloEnergy::MOP);
427  }
428  // set parametrized energy loss in em
429  double emLoss = caloEnergy->deltaE() * emMopFraction(eta);
430  caloEnergy->set_fsrCandidateEnergy(emLoss);
431  return caloEnergy;
432  }

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

◆ symmetricMopEnergyLoss()

double MuidCaloEnergyParam::symmetricMopEnergyLoss ( double  eta,
double  momentum 
) const
private

Definition at line 605 of file MuidCaloEnergyParam.cxx.

605  {
606  // Mop Energy Loss parametrization with symmetric gaussian fit in the mop region
607  // new coefficients with symmetric fits
608  MuidCaloEnergyParam::BinsWeights BW = etaVariableBin(eta);
609  double p0 = BW.etaWeight1 * m_mopSymmetricEnergyLossP0[BW.etaBin1] + BW.etaWeight2 * m_mopSymmetricEnergyLossP0[BW.etaBin2];
610  double p1 = BW.etaWeight1 * m_mopSymmetricEnergyLossP1[BW.etaBin1] + BW.etaWeight2 * m_mopSymmetricEnergyLossP1[BW.etaBin2];
611  double p2 = BW.etaWeight1 * m_mopSymmetricEnergyLossP2[BW.etaBin1] + BW.etaWeight2 * m_mopSymmetricEnergyLossP2[BW.etaBin2];
612  double parametrisedDeposit = p0 * Units::GeV + p1 * Units::GeV * log(0.0067 * std::pow(momentum * PerGeV, 1)) + p2 * momentum;
613  return parametrisedDeposit;
614  }

◆ 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.

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

◆ x0mapEmMaterial()

double MuidCaloEnergyParam::x0mapEmMaterial ( double  eta) const
virtual

IMuidCaloEnergyParam interface: returns the fraction of em calorimeter inert material between InDet and MS.

Implements Rec::IMuidCaloEnergyParam.

Definition at line 630 of file MuidCaloEnergyParam.cxx.

630  {
631  int index = static_cast<int>(floor(std::abs(eta) / m_etaGranularity));
632  if (index >= 277) index = 276;
633 
634  double x01 = m_dead_x01[index];
635  double x02 = m_dead_x02[index];
636  double x0tot = m_dead_x0tot[index];
637  double x0cryo1 = m_dead_cryo1[index];
638  // double x0cryo2 = m_dead_cryo2[index];
639  // double x0girder = m_dead_girder[index];
640  double sum = x02 - x01 - x0cryo1;
641  return sum / x0tot;
642  }

◆ x0mapHecMaterial()

double MuidCaloEnergyParam::x0mapHecMaterial ( double  eta) const
virtual

IMuidCaloEnergyParam interface: returns the fraction of hec calorimeter inert material between InDet and MS.

Implements Rec::IMuidCaloEnergyParam.

Definition at line 644 of file MuidCaloEnergyParam.cxx.

644  {
645  int index = static_cast<int>(floor(std::abs(eta) / m_etaGranularity));
646  if (index >= 277) index = 276;
647 
648  // double x01 = m_dead_x01[index];
649  double x02 = m_dead_x02[index];
650  double x0tot = m_dead_x0tot[index];
651  // double x0cryo1 = m_dead_cryo1[index];
652  double x0cryo2 = m_dead_cryo2[index];
653  double x0girder = m_dead_girder[index];
654  double sum = x0tot - x02 - x0cryo2 - x0girder;
655  return sum / x0tot;
656  }

◆ x0mapInertMaterial()

double MuidCaloEnergyParam::x0mapInertMaterial ( double  eta) const
virtual

IMuidCaloEnergyParam interface: returns the fraction of inert material between InDet and MS.

Implements Rec::IMuidCaloEnergyParam.

Definition at line 616 of file MuidCaloEnergyParam.cxx.

616  {
617  int index = static_cast<int>(floor(std::abs(eta) / m_etaGranularity));
618  if (index >= 277) index = 276;
619 
620  double x01 = m_dead_x01[index];
621  // double x02 = m_dead_x02[index];
622  double x0tot = m_dead_x0tot[index];
623  double x0cryo1 = m_dead_cryo1[index];
624  double x0cryo2 = m_dead_cryo2[index];
625  double x0girder = m_dead_girder[index];
626  double sum = x01 + x0cryo1 + x0cryo2 + x0girder;
627  return sum / x0tot;
628  }

Member Data Documentation

◆ m_binWidth

double Rec::MuidCaloEnergyParam::m_binWidth
private

Definition at line 89 of file MuidCaloEnergyParam.h.

◆ m_cosmics

bool Rec::MuidCaloEnergyParam::m_cosmics
private

Definition at line 85 of file MuidCaloEnergyParam.h.

◆ m_dead_cryo0

double Rec::MuidCaloEnergyParam::m_dead_cryo0[277] {}
private

Definition at line 120 of file MuidCaloEnergyParam.h.

◆ m_dead_cryo1

double Rec::MuidCaloEnergyParam::m_dead_cryo1[277] {}
private

Definition at line 121 of file MuidCaloEnergyParam.h.

◆ m_dead_cryo2

double Rec::MuidCaloEnergyParam::m_dead_cryo2[277] {}
private

Definition at line 122 of file MuidCaloEnergyParam.h.

◆ m_dead_eta

double Rec::MuidCaloEnergyParam::m_dead_eta[277] {}
private

Definition at line 116 of file MuidCaloEnergyParam.h.

◆ m_dead_girder

double Rec::MuidCaloEnergyParam::m_dead_girder[277] {}
private

Definition at line 123 of file MuidCaloEnergyParam.h.

◆ m_dead_x01

double Rec::MuidCaloEnergyParam::m_dead_x01[277] {}
private

Definition at line 117 of file MuidCaloEnergyParam.h.

◆ m_dead_x02

double Rec::MuidCaloEnergyParam::m_dead_x02[277] {}
private

Definition at line 118 of file MuidCaloEnergyParam.h.

◆ m_dead_x0tot

double Rec::MuidCaloEnergyParam::m_dead_x0tot[277] {}
private

Definition at line 119 of file MuidCaloEnergyParam.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_emMopFraction

double Rec::MuidCaloEnergyParam::m_emMopFraction[26] {}
private

Definition at line 112 of file MuidCaloEnergyParam.h.

◆ m_etaGranularity

double Rec::MuidCaloEnergyParam::m_etaGranularity
private

Definition at line 115 of file MuidCaloEnergyParam.h.

◆ m_etaOffset

double Rec::MuidCaloEnergyParam::m_etaOffset
private

Definition at line 90 of file MuidCaloEnergyParam.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_inverseWidth

double Rec::MuidCaloEnergyParam::m_inverseWidth
private

Definition at line 91 of file MuidCaloEnergyParam.h.

◆ m_meanEnergyLossErrorLeftP0

double Rec::MuidCaloEnergyParam::m_meanEnergyLossErrorLeftP0[26] {}
private

Definition at line 97 of file MuidCaloEnergyParam.h.

◆ m_meanEnergyLossErrorLeftP1

double Rec::MuidCaloEnergyParam::m_meanEnergyLossErrorLeftP1[26] {}
private

Definition at line 98 of file MuidCaloEnergyParam.h.

◆ m_meanEnergyLossErrorRightP0

double Rec::MuidCaloEnergyParam::m_meanEnergyLossErrorRightP0[26] {}
private

Definition at line 99 of file MuidCaloEnergyParam.h.

◆ m_meanEnergyLossErrorRightP1

double Rec::MuidCaloEnergyParam::m_meanEnergyLossErrorRightP1[26] {}
private

Definition at line 100 of file MuidCaloEnergyParam.h.

◆ m_meanEnergyLossP0

double Rec::MuidCaloEnergyParam::m_meanEnergyLossP0[26] {}
private

Definition at line 94 of file MuidCaloEnergyParam.h.

◆ m_meanEnergyLossP1

double Rec::MuidCaloEnergyParam::m_meanEnergyLossP1[26] {}
private

Definition at line 95 of file MuidCaloEnergyParam.h.

◆ m_meanEnergyLossP2

double Rec::MuidCaloEnergyParam::m_meanEnergyLossP2[26] {}
private

Definition at line 96 of file MuidCaloEnergyParam.h.

◆ m_mopEnergyLossErrorP0

double Rec::MuidCaloEnergyParam::m_mopEnergyLossErrorP0[26] {}
private

Definition at line 105 of file MuidCaloEnergyParam.h.

◆ m_mopEnergyLossErrorP1

double Rec::MuidCaloEnergyParam::m_mopEnergyLossErrorP1[26] {}
private

Definition at line 106 of file MuidCaloEnergyParam.h.

◆ m_mopEnergyLossP0

double Rec::MuidCaloEnergyParam::m_mopEnergyLossP0[26] {}
private

Definition at line 102 of file MuidCaloEnergyParam.h.

◆ m_mopEnergyLossP1

double Rec::MuidCaloEnergyParam::m_mopEnergyLossP1[26] {}
private

Definition at line 103 of file MuidCaloEnergyParam.h.

◆ m_mopEnergyLossP2

double Rec::MuidCaloEnergyParam::m_mopEnergyLossP2[26] {}
private

Definition at line 104 of file MuidCaloEnergyParam.h.

◆ m_mopSymmetricEnergyLossP0

double Rec::MuidCaloEnergyParam::m_mopSymmetricEnergyLossP0[26] {}
private

Definition at line 108 of file MuidCaloEnergyParam.h.

◆ m_mopSymmetricEnergyLossP1

double Rec::MuidCaloEnergyParam::m_mopSymmetricEnergyLossP1[26] {}
private

Definition at line 109 of file MuidCaloEnergyParam.h.

◆ m_mopSymmetricEnergyLossP2

double Rec::MuidCaloEnergyParam::m_mopSymmetricEnergyLossP2[26] {}
private

Definition at line 110 of file MuidCaloEnergyParam.h.

◆ m_smoothingFraction

double Rec::MuidCaloEnergyParam::m_smoothingFraction
private

Definition at line 86 of file MuidCaloEnergyParam.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.

The documentation for this class was generated from the following files:
str::trackMomentum
const std::string trackMomentum
Definition: BTagTrackIpAccessor.cxx:15
Rec::MuidCaloEnergyParam::m_dead_girder
double m_dead_girder[277]
Definition: MuidCaloEnergyParam.h:123
plotting.yearwise_luminosity_vs_mu.comp
comp
Definition: yearwise_luminosity_vs_mu.py:24
Rec::MuidCaloEnergyParam::m_cosmics
bool m_cosmics
Definition: MuidCaloEnergyParam.h:85
Rec::MuidCaloEnergyParam::mopEnergyLossError
double mopEnergyLossError(double, double) const
Definition: MuidCaloEnergyParam.cxx:588
pdg_comparison.sigma
sigma
Definition: pdg_comparison.py:324
Rec::MuidCaloEnergyParam::m_inverseWidth
double m_inverseWidth
Definition: MuidCaloEnergyParam.h:91
python.PerfMonSerializer.p
def p
Definition: PerfMonSerializer.py:743
ParticleGun_SamplingFraction.eta2
eta2
Definition: ParticleGun_SamplingFraction.py:96
phi
Scalar phi() const
phi method
Definition: AmgMatrixBasePlugin.h:64
Rec::MuidCaloEnergyParam::m_mopEnergyLossP0
double m_mopEnergyLossP0[26]
Definition: MuidCaloEnergyParam.h:102
CaloEnergy
class extending the basic Trk::EnergyLoss to describe the measured or parameterised muon energy loss ...
Definition: CaloEnergy.h:28
Rec::MuidCaloEnergyParam::caloCompartmentDepthTile
static double caloCompartmentDepthTile(int)
Definition: MuidCaloEnergyParam.cxx:676
eta
Scalar eta() const
pseudorapidity method
Definition: AmgMatrixBasePlugin.h:79
index
Definition: index.py:1
AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
Rec::MuidCaloEnergyParam::m_meanEnergyLossP1
double m_meanEnergyLossP1[26]
Definition: MuidCaloEnergyParam.h:95
Rec::MuidCaloEnergyParam::emMopFraction
double emMopFraction(double eta) const
IMuidCaloEnergyParam interface: fraction of the full mop deposition corresponding to the em.
Definition: MuidCaloEnergyParam.cxx:704
conifer::pow
constexpr int pow(int x)
Definition: conifer.h:20
Rec::MuidCaloEnergyParam::m_mopEnergyLossErrorP1
double m_mopEnergyLossErrorP1[26]
Definition: MuidCaloEnergyParam.h:106
python.SystemOfUnits.MeV
int MeV
Definition: SystemOfUnits.py:154
test_pyathena.pt
pt
Definition: test_pyathena.py:11
Rec::MuidCaloEnergyParam::etaFixedBin
BinsWeights etaFixedBin(double) const
Definition: MuidCaloEnergyParam.cxx:460
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
xAOD::eta1
setEt setPhi setE277 setWeta2 eta1
Definition: TrigEMCluster_v1.cxx:41
Rec::MuidCaloEnergyParam::m_dead_cryo2
double m_dead_cryo2[277]
Definition: MuidCaloEnergyParam.h:122
Rec::MuidCaloEnergyParam::m_meanEnergyLossErrorRightP0
double m_meanEnergyLossErrorRightP0[26]
Definition: MuidCaloEnergyParam.h:99
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
yodamerge_tmp.scale
scale
Definition: yodamerge_tmp.py:138
Rec::MuidCaloEnergyParam::m_meanEnergyLossErrorLeftP0
double m_meanEnergyLossErrorLeftP0[26]
Definition: MuidCaloEnergyParam.h:97
Rec::MuidCaloEnergyParam::m_etaGranularity
double m_etaGranularity
Definition: MuidCaloEnergyParam.h:115
Rec::MuidCaloEnergyParam::m_dead_x01
double m_dead_x01[277]
Definition: MuidCaloEnergyParam.h:117
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
Rec::MuidCaloEnergyParam::m_binWidth
double m_binWidth
Definition: MuidCaloEnergyParam.h:89
AthCommonDataStore
Definition: AthCommonDataStore.h:52
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
convertTimingResiduals.sum
sum
Definition: convertTimingResiduals.py:55
ParticleGun_EoverP_Config.momentum
momentum
Definition: ParticleGun_EoverP_Config.py:63
lumiFormat.i
int i
Definition: lumiFormat.py:92
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
Rec::MuidCaloEnergyParam::meanEnergyLossError
std::pair< double, double > meanEnergyLossError(double, double) const
Definition: MuidCaloEnergyParam.cxx:543
Rec::MuidCaloEnergyParam::m_mopEnergyLossErrorP0
double m_mopEnergyLossErrorP0[26]
Definition: MuidCaloEnergyParam.h:105
Trk::EnergyLoss::deltaE
double deltaE() const
returns the
test_pyathena.parent
parent
Definition: test_pyathena.py:15
Rec::MuidCaloEnergyParam::m_meanEnergyLossErrorLeftP1
double m_meanEnergyLossErrorLeftP1[26]
Definition: MuidCaloEnergyParam.h:98
Rec::MuidCaloEnergyParam::m_dead_eta
double m_dead_eta[277]
Definition: MuidCaloEnergyParam.h:116
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
Rec::MuidCaloEnergyParam::m_smoothingFraction
double m_smoothingFraction
Definition: MuidCaloEnergyParam.h:86
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
Rec::MuidCaloEnergyParam::caloCompartmentDepthLArEM
static double caloCompartmentDepthLArEM(int)
Definition: MuidCaloEnergyParam.cxx:695
Rec::MuidCaloEnergyParam::etaVariableBin
BinsWeights etaVariableBin(double) const
Definition: MuidCaloEnergyParam.cxx:490
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
Rec::MuidCaloEnergyParam::m_dead_cryo0
double m_dead_cryo0[277]
Definition: MuidCaloEnergyParam.h:120
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:192
Rec::MuidCaloEnergyParam::caloCompartmentDepthLArHEC
static double caloCompartmentDepthLArHEC(int)
Definition: MuidCaloEnergyParam.cxx:685
Rec::MuidCaloEnergyParam::m_dead_cryo1
double m_dead_cryo1[277]
Definition: MuidCaloEnergyParam.h:121
Rec::MuidCaloEnergyParam::m_mopSymmetricEnergyLossP0
double m_mopSymmetricEnergyLossP0[26]
Definition: MuidCaloEnergyParam.h:108
Rec::MuidCaloEnergyParam::m_emMopFraction
double m_emMopFraction[26]
Definition: MuidCaloEnergyParam.h:112
DeMoScan.index
string index
Definition: DeMoScan.py:362
CaloEnergy::set_fsrCandidateEnergy
void set_fsrCandidateEnergy(const float fs)
FSR Candidate Energy.
Definition: CaloEnergy.h:118
a
TList * a
Definition: liststreamerinfos.cxx:10
h
CaloEnergy::Parametrized
@ Parametrized
Definition: CaloEnergy.h:43
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
python.CaloScaleNoiseConfig.type
type
Definition: CaloScaleNoiseConfig.py:78
Rec::MuidCaloEnergyParam::m_meanEnergyLossErrorRightP1
double m_meanEnergyLossErrorRightP1[26]
Definition: MuidCaloEnergyParam.h:100
Rec::MuidCaloEnergyParam::m_etaOffset
double m_etaOffset
Definition: MuidCaloEnergyParam.h:90
Rec::MuidCaloEnergyParam::m_mopSymmetricEnergyLossP2
double m_mopSymmetricEnergyLossP2[26]
Definition: MuidCaloEnergyParam.h:110
python.CaloCondTools.log
log
Definition: CaloCondTools.py:20
Rec::MuidCaloEnergyParam::m_mopEnergyLossP1
double m_mopEnergyLossP1[26]
Definition: MuidCaloEnergyParam.h:103
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:616
Rec::MuidCaloEnergyParam::m_mopEnergyLossP2
double m_mopEnergyLossP2[26]
Definition: MuidCaloEnergyParam.h:104
declareProperty
#define declareProperty(n, p, h)
Definition: BaseFakeBkgTool.cxx:15
Rec::MuidCaloEnergyParam::m_mopSymmetricEnergyLossP1
double m_mopSymmetricEnergyLossP1[26]
Definition: MuidCaloEnergyParam.h:109
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:790
CaloEnergy::MOP
@ MOP
Definition: CaloEnergy.h:43
MuonParameters::beta
@ beta
Definition: MuonParamDefs.h:144
Rec::MuidCaloEnergyParam::mopEnergyLoss
double mopEnergyLoss(double, double) const
Definition: MuidCaloEnergyParam.cxx:557
GeV
#define GeV
Definition: CaloTransverseBalanceVecMon.cxx:30
Rec::MuidCaloEnergyParam::m_meanEnergyLossP2
double m_meanEnergyLossP2[26]
Definition: MuidCaloEnergyParam.h:96
Rec::MuidCaloEnergyParam::m_dead_x02
double m_dead_x02[277]
Definition: MuidCaloEnergyParam.h:118
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
fitman.k
k
Definition: fitman.py:528
Rec::MuidCaloEnergyParam::m_dead_x0tot
double m_dead_x0tot[277]
Definition: MuidCaloEnergyParam.h:119
Rec::MuidCaloEnergyParam::m_meanEnergyLossP0
double m_meanEnergyLossP0[26]
Definition: MuidCaloEnergyParam.h:94
Rec::MuidCaloEnergyParam::meanEnergyLoss
double meanEnergyLoss(double, double) const
Definition: MuidCaloEnergyParam.cxx:513
Generated on Fri Jun 28 2024 22:27:01 for ATLAS Offline Software by   doxygen 1.8.18