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

This class implements the CaloDetDescr/ICaloRecoMaterialTool and provides the average materials needed for the LAr part of the Calorimeter/CaloTrackingGeometry volumes. More...

#include <LArRecoMaterialTool.h>

Inheritance diagram for LArRecoMaterialTool:
Collaboration diagram for LArRecoMaterialTool:

Public Member Functions

 LArRecoMaterialTool (const std::string &type, const std::string &name, const IInterface *parent)
 
 ~LArRecoMaterialTool ()
 
virtual StatusCode initialize () override
 
virtual StatusCode finalize () override
 
virtual bool get_material (CaloSubdetNames::ALIGNVOL alvol, double &mass, double &volume, double &x0, double &dEdx, double &aveA, double &aveZ) const override
 the enum ALIGNVOL covers all LAr pieces, including cryostats More...
 
virtual void get_default_material (double &density, double &x0, double &dEdx, double &aveA, double &aveZ) const override
 default material used to fill the envelope More...
 
void print () const override
 printout method 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 ()
 

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

void ScanTree ()
 Private method called at initialisation : fills all arrays by calling the "optimised" methods. More...
 
bool ScanCRYO (CaloSubdetNames::ALIGNVOL alvol, double &mass, double &x0, double &dEdx, double &aveA, double &aveZ)
 Private optimised methods for each subsystem. More...
 
bool ScanPS (CaloSubdetNames::ALIGNVOL alvol, double &mass, double &x0, double &dEdx, double &aveA, double &aveZ)
 
bool ScanEMB (CaloSubdetNames::ALIGNVOL alvol, double &mass, double &x0, double &dEdx, double &aveA, double &aveZ)
 
bool ScanEMEC (CaloSubdetNames::ALIGNVOL alvol, double &mass, double &x0, double &dEdx, double &aveA, double &aveZ)
 
bool ScanHEC (CaloSubdetNames::ALIGNVOL alvol, double &mass, double &x0, double &dEdx, double &aveA, double &aveZ)
 
bool ScanFCAL (CaloSubdetNames::ALIGNVOL alvol, double &mass, double &x0, double &dEdx, double &aveA, double &aveZ)
 
void blendGeoVolume (GeoPVConstLink &geoVolLink, std::vector< double > &volume, std::vector< double > &mass, std::vector< double > &x0, std::vector< double > &dEdX, std::vector< double > &aveA, std::vector< double > &aveZ) const
 Private method to blend GeoLogVolume. More...
 
void addMaterialFraction (const GeoLogVol &geoVol, std::vector< double > &volume, std::vector< double > &mass, std::vector< double > &x0, std::vector< double > &dEdX, std::vector< double > &aveA, std::vector< double > &aveZ) const
 Private method to add material. More...
 
void averageFraction (std::vector< double > &volumeFractions, std::vector< double > &massFractions, std::vector< double > &x0Fractions, std::vector< double > &dEdXFractions, std::vector< double > &aveAFractions, std::vector< double > &aveZFractions, double &mass, double &volume, double &x0, double &dEdx, double &aveA, double &aveZ, double &rho)
 Private method to average results stored in private arrays. More...
 
std::string map_av (CaloSubdetNames::ALIGNVOL alvol)
 
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...
 

Private Attributes

bool m_parseGeoModelForMaterial = false
 Private bool set at initialisation : decides if the GeoModel is parsed or bypassed. More...
 
const CaloCell_ID * m_calo_id
 
const CaloDM_ID * m_dm_id
 
std::vector< CaloSubdetNames::ALIGNVOL > m_alignvol_number
 
std::vector< double > m_alignvol_mass
 
std::vector< double > m_alignvol_x0
 
std::vector< double > m_alignvol_dEdX
 
std::vector< double > m_alignvol_aveA
 
std::vector< double > m_alignvol_aveZ
 
std::vector< double > m_massFractions
 
std::vector< double > m_volumeFractions
 
std::vector< double > m_x0Fractions
 
std::vector< double > m_dEdXFractions
 
std::vector< double > m_averageAFractions
 
std::vector< double > m_averageZFractions
 
std::vector< double > m_child_massFractions
 
std::vector< double > m_child_volumeFractions
 
std::vector< double > m_child_x0Fractions
 
std::vector< double > m_child_dEdXFractions
 
std::vector< double > m_child_averageAFractions
 
std::vector< double > m_child_averageZFractions
 
std::vector< double > m_child_Rho
 
std::vector< double > m_child_Volume
 
std::vector< double > m_child_Mass
 
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

This class implements the CaloDetDescr/ICaloRecoMaterialTool and provides the average materials needed for the LAr part of the Calorimeter/CaloTrackingGeometry volumes.

For this, it uses the PVLinks provided by the LArSubdetPosHelper.

Note that it means that eta/phi dependencies of the material within a given GeoModel volume is not supported : although an interface using the CaloDM identifiers is foreseen, numbers will be averaged.

Interface and implementation are preliminary !!!!

Definition at line 43 of file LArRecoMaterialTool.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

◆ LArRecoMaterialTool()

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

Definition at line 34 of file LArRecoMaterialTool.cxx.

36  :
37  AthAlgTool(type, name, parent),
38  m_calo_id(nullptr),m_dm_id(nullptr)
39 {
40  declareInterface<LArRecoMaterialTool>( this );
41 }

◆ ~LArRecoMaterialTool()

LArRecoMaterialTool::~LArRecoMaterialTool ( )

Definition at line 75 of file LArRecoMaterialTool.cxx.

76 {
77 }

Member Function Documentation

◆ addMaterialFraction()

void LArRecoMaterialTool::addMaterialFraction ( const GeoLogVol &  geoVol,
std::vector< double > &  volume,
std::vector< double > &  mass,
std::vector< double > &  x0,
std::vector< double > &  dEdX,
std::vector< double > &  aveA,
std::vector< double > &  aveZ 
) const
private

Private method to add material.

Definition at line 1069 of file LArRecoMaterialTool.cxx.

1076 {
1077  const GeoShape* childShape = geoVol.getShape();
1078  const GeoMaterial* childMaterial = geoVol.getMaterial();
1079 
1080  double childVolume = childShape->volume();
1081 
1082  double childNumberOfElements = childMaterial->getNumElements();
1083  double childX0 = childMaterial->getRadLength ();
1084  double childDeDx = childMaterial->getDeDxMin();
1085  double childRho = childMaterial->getDensity();
1086 
1087  ATH_MSG_DEBUG ("addMaterialFraction : " << geoVol.getName()
1088  << " has volume=" << childVolume << " mm**3 and "
1089  << childNumberOfElements << " elements ");
1090 
1091  for (unsigned int iEl=0; iEl<childNumberOfElements; iEl++){
1092  const GeoElement* geoEl = childMaterial->getElement(iEl);
1093  double fraction = childMaterial->getFraction(iEl);
1094 
1095  ATH_MSG_DEBUG (" direct child : " << geoEl->getName() << " fraction = " << fraction
1096  << " A= " << geoEl->getA()*(1./GeoModelKernelUnits::gram) << " Z= " << geoEl->getZ());
1097 
1098  volume.push_back(fraction*childVolume);
1099  mass.push_back(fraction*childVolume*childRho);
1100  x0.push_back(fraction*childX0);
1101  dEdX.push_back(fraction*childDeDx);
1102  aveA.push_back(fraction*(geoEl->getA()*(1./GeoModelKernelUnits::gram)));
1103  aveZ.push_back(fraction*(geoEl->getZ()));
1104  }
1105 
1106 }

◆ averageFraction()

void LArRecoMaterialTool::averageFraction ( std::vector< double > &  volumeFractions,
std::vector< double > &  massFractions,
std::vector< double > &  x0Fractions,
std::vector< double > &  dEdXFractions,
std::vector< double > &  aveAFractions,
std::vector< double > &  aveZFractions,
double &  mass,
double &  volume,
double &  x0,
double &  dEdx,
double &  aveA,
double &  aveZ,
double &  rho 
)
private

Private method to average results stored in private arrays.

Definition at line 1109 of file LArRecoMaterialTool.cxx.

1118 {
1119  unsigned int parsedVolumes = massFractions.size();
1120  ATH_MSG_DEBUG (" averageFraction : will average over... "
1121  << parsedVolumes << " elements ");
1122 
1123  childMass = 0.;
1124  childVolume = 0.;
1125  childX0 = 0.;
1126  childDeDx = 0.;
1127  childAverageA = 0.;
1128  childAverageZ = 0.;
1129  childRho = 0.;
1130 
1131  for (unsigned int ivols = 0; ivols < parsedVolumes; ++ivols) {
1132  childMass += massFractions[ivols];
1133  childVolume += volumeFractions[ivols];
1134  }
1135 
1136  for (unsigned int ivolfrac = 0; ivolfrac < parsedVolumes; ++ivolfrac){
1137  double volfrac = volumeFractions[ivolfrac]/childVolume;
1138  childX0 += volfrac * x0Fractions[ivolfrac];
1139  childDeDx += volfrac * dEdXFractions[ivolfrac];
1140  double massfrac = massFractions[ivolfrac]/childMass;
1141  childAverageA += massfrac * averageAFractions[ivolfrac];
1142  childAverageZ += massfrac * averageZFractions[ivolfrac];
1143  }
1144 
1145  childRho = childMass/(GeoModelKernelUnits::gram*childVolume);
1146 
1147  ATH_MSG_DEBUG ("");
1148  ATH_MSG_DEBUG (" + averaged over " << parsedVolumes << " volumes ");
1149  ATH_MSG_DEBUG (" - volume [mm^3] : " << childVolume);
1150  ATH_MSG_DEBUG (" - mass [gram] : " << childMass*(1./GeoModelKernelUnits::gram));
1151  ATH_MSG_DEBUG (" - rho [gram/mm^3] : " << childRho);
1152  ATH_MSG_DEBUG (" - X0 : " << childX0);
1153  ATH_MSG_DEBUG (" - DeDx : " << childDeDx);
1154  ATH_MSG_DEBUG (" - AverageA : " << childAverageA);
1155  ATH_MSG_DEBUG (" - AverageZ : " << childAverageZ);
1156  ATH_MSG_DEBUG ("");
1157 
1158 }

◆ blendGeoVolume()

void LArRecoMaterialTool::blendGeoVolume ( GeoPVConstLink &  geoVolLink,
std::vector< double > &  volume,
std::vector< double > &  mass,
std::vector< double > &  x0,
std::vector< double > &  dEdX,
std::vector< double > &  aveA,
std::vector< double > &  aveZ 
) const
private

Private method to blend GeoLogVolume.

Definition at line 1048 of file LArRecoMaterialTool.cxx.

1055 {
1056 
1057  const GeoLogVol* currentVol = geoVolLink->getLogVol();
1058  if (currentVol){
1059  this->addMaterialFraction(*currentVol, volume, mass, x0, dEdX, aveA, aveZ);
1060  unsigned int currentVolChilds = geoVolLink->getNChildVols();
1061  for (unsigned int ichild = 0; ichild < currentVolChilds; ichild++){
1062  GeoPVConstLink childVolumeLink = geoVolLink->getChildVol(ichild);
1063  this->blendGeoVolume(childVolumeLink, volume, mass, x0, dEdX, aveA, aveZ);
1064  }
1065  }
1066 }

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

◆ 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

◆ finalize()

StatusCode LArRecoMaterialTool::finalize ( )
overridevirtual

Definition at line 80 of file LArRecoMaterialTool.cxx.

81 {
82  return StatusCode::SUCCESS;
83 }

◆ get_default_material()

void LArRecoMaterialTool::get_default_material ( double &  density,
double &  x0,
double &  dEdx,
double &  aveA,
double &  aveZ 
) const
overridevirtual

default material used to fill the envelope

Implements ICaloRecoMaterialTool.

Definition at line 113 of file LArRecoMaterialTool.cxx.

117 {
118  // FIXME ! should take real LAr numbers
119  density = 1.;
120  x0 = 5.;
121  dEdx = 1.1e-24;
122  aveA = 39.948;
123  aveZ = 18.;
124 
125 }

◆ get_material()

bool LArRecoMaterialTool::get_material ( CaloSubdetNames::ALIGNVOL  alvol,
double &  mass,
double &  volume,
double &  x0,
double &  dEdx,
double &  aveA,
double &  aveZ 
) const
overridevirtual

the enum ALIGNVOL covers all LAr pieces, including cryostats

Implements ICaloRecoMaterialTool.

Definition at line 87 of file LArRecoMaterialTool.cxx.

92 {
93  bool result = false;
94 
95  for ( unsigned int i = 0; i < m_alignvol_number.size(); i++ )
96  {
97  if ( m_alignvol_number[i] == alvol ) {
98  mass = m_alignvol_mass[i];
99  x0 = m_alignvol_x0[i];
100  dEdx = m_alignvol_dEdX[i];
101  aveA = m_alignvol_aveA[i];
102  aveZ = m_alignvol_aveZ[i];
103  // FIXME
104  volume = 1000.;
105  result = true;
106  }
107  }
108 
109  return result;
110 }

◆ initialize()

StatusCode LArRecoMaterialTool::initialize ( )
overridevirtual

Definition at line 44 of file LArRecoMaterialTool.cxx.

45 {
46  ATH_CHECK( detStore()->retrieve (m_calo_id, "CaloCell_ID") );
47  ATH_CHECK( detStore()->retrieve (m_dm_id, "CaloDM_ID") );
48 
49  // These vectors will contain the results :
50  m_alignvol_number.resize(0);
51  m_alignvol_mass.resize(0);
52  m_alignvol_x0.resize(0);
53  m_alignvol_dEdX.resize(0);
54  m_alignvol_aveA.resize(0);
55  m_alignvol_aveZ.resize(0);
56 
57  // ---------- Fill arrays :
58 
59  // going though the full tree is too slow to be used as default,
60  // if you want to do it set this flag to true :
61  // ( will become a property at some point ! )
62  m_parseGeoModelForMaterial = false;
63 
64  this->ScanTree();
65  //this->print();
66 
67 
68  // ---------- That's it
69 
70 
71  ATH_MSG_INFO (" LArRecoMaterialTool successfully initialized ");
72  return StatusCode::SUCCESS;
73 }

◆ 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()

const InterfaceID & LArRecoMaterialTool::interfaceID ( )
static

Definition at line 31 of file LArRecoMaterialTool.cxx.

32 { return IID_LArRecoMaterialTool; }

◆ map_av()

std::string LArRecoMaterialTool::map_av ( CaloSubdetNames::ALIGNVOL  alvol)
private

Definition at line 1161 of file LArRecoMaterialTool.cxx.

1162 {
1163  switch(alvol)
1164  {
1165  case CaloSubdetNames::LARCRYO_B:
1166  return std::string("LARCRYO_B");
1167  case CaloSubdetNames::LARCRYO_EC_POS:
1168  return std::string("LARCRYO_EC_POS");
1169  case CaloSubdetNames::LARCRYO_EC_NEG:
1170  return std::string("LARCRYO_EC_NEG");
1171  case CaloSubdetNames::SOLENOID:
1172  return std::string("SOLENOID");
1173  case CaloSubdetNames::PRESAMPLER_B_POS:
1174  return std::string("PRESAMPLER_B_POS");
1175  case CaloSubdetNames::PRESAMPLER_B_NEG:
1176  return std::string("PRESAMPLER_B_NEG");
1177  case CaloSubdetNames::EMB_POS:
1178  return std::string("EMB_POS");
1179  case CaloSubdetNames::EMB_NEG:
1180  return std::string("EMB_NEG");
1181  case CaloSubdetNames::PRESAMPLER_EC_POS:
1182  return std::string("PRESAMPLER_EC_POS");
1183  case CaloSubdetNames::PRESAMPLER_EC_NEG:
1184  return std::string("PRESAMPLER_EC_NEG");
1185  case CaloSubdetNames::EMEC_POS:
1186  return std::string("EMEC_POS");
1187  case CaloSubdetNames::EMEC_NEG:
1188  return std::string("EMEC_NEG");
1189  case CaloSubdetNames::HEC1_POS:
1190  return std::string("HEC1_POS");
1191  case CaloSubdetNames::HEC1_NEG:
1192  return std::string("HEC1_NEG");
1193  case CaloSubdetNames::HEC2_POS:
1194  return std::string("HEC2_POS");
1195  case CaloSubdetNames::HEC2_NEG:
1196  return std::string("HEC2_NEG");
1197  case CaloSubdetNames::FCAL1_POS:
1198  return std::string("FCAL1_POS");
1199  case CaloSubdetNames::FCAL1_NEG:
1200  return std::string("FCAL1_NEG");
1201  case CaloSubdetNames::FCAL2_POS:
1202  return std::string("FCAL2_POS");
1203  case CaloSubdetNames::FCAL2_NEG:
1204  return std::string("FCAL2_NEG");
1205  case CaloSubdetNames::FCAL3_POS:
1206  return std::string("FCAL3_POS");
1207  case CaloSubdetNames::FCAL3_NEG:
1208  return std::string("FCAL3_NEG");
1209  default:
1210  return std::string("");
1211  }
1212 }

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

◆ print()

void LArRecoMaterialTool::print ( ) const
overridevirtual

printout method

Implements ICaloRecoMaterialTool.

Definition at line 128 of file LArRecoMaterialTool.cxx.

129 {
130 
131  for ( unsigned int i = 0; i < m_alignvol_number.size(); i++ )
132  {
133  ATH_MSG_INFO (" + Result for ALIGNVOL " << m_alignvol_number[i]);
134  ATH_MSG_INFO (" - mass [gram] " << m_alignvol_mass[i]);
135  ATH_MSG_INFO (" - X0 " << m_alignvol_x0[i]);
136  ATH_MSG_INFO (" - DeDx " << m_alignvol_dEdX[i]);
137  ATH_MSG_INFO (" - AverageA " << m_alignvol_aveA[i]);
138  ATH_MSG_INFO (" - AverageZ " << m_alignvol_aveZ[i]);
139  }
140 }

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

◆ ScanCRYO()

bool LArRecoMaterialTool::ScanCRYO ( CaloSubdetNames::ALIGNVOL  alvol,
double &  mass,
double &  x0,
double &  dEdx,
double &  aveA,
double &  aveZ 
)
private

Private optimised methods for each subsystem.

Definition at line 191 of file LArRecoMaterialTool.cxx.

194 {
195  mass = 0.;
196  x0 = 0.;
197  dEdx = 0.;
198  aveA = 0.;
199  aveZ = 0.;
200 
201  if ( alvol != CaloSubdetNames::LARCRYO_B
202  && alvol != CaloSubdetNames::LARCRYO_EC_POS
203  && alvol != CaloSubdetNames::LARCRYO_EC_NEG
204  && alvol != CaloSubdetNames::SOLENOID ) {
205  ATH_MSG_INFO ("method ScanCryo does not support CaloSubdetNames::ALIGNVOL " << (int) alvol);
206  return false;
207  }
208 
209  // For EC, link to GeoModel does not exist (yet) : take barrel numbers
210 
211  // Shortcut :
212  if ( !m_parseGeoModelForMaterial )
213  {
214  //if ( alvol == CaloSubdetNames::LARCRYO_B || alvol == CaloSubdetNames::SOLENOID) {
215  mass = 3.61158e+07 ;
216  x0 = 89.411;
217  dEdx = 7.47033e-24 ;
218  aveA = 27.9728;
219  aveZ = 13.4253 ;
220 #if 0
221  } else {
222  mass = 3.61158e+07 ;
223  x0 = 89.411;
224  dEdx = 7.47033e-24 ;
225  aveA = 27.9728;
226  aveZ = 13.4253 ;
227  }
228 #endif
229  }
230 #if 0
231  // Dead code as long as m_parseGeoModelForMaterial is false.
232  // Comment out to avoid coverity warnings.
233  else {
234 
235  // Go for it :
236  m_child_massFractions.clear();
237  m_child_volumeFractions.clear();
238  m_child_x0Fractions.clear();
239  m_child_dEdXFractions.clear();
240  m_child_averageAFractions.clear();
241  m_child_averageZFractions.clear();
242  m_child_Rho.clear();
243  m_child_Volume.clear();
244  m_child_Mass.clear();
245 
246  double childVolume = 0.;
247  double childMass = 0.;
248  double childX0 = 0.;
249  double childDeDx = 0.;
250  double childAverageA = 0.;
251  double childAverageZ = 0.;
252  double childRho = 0.;
253 
254  const StoredPhysVol* storedPV = 0;
255  std::string key = map_av(alvol);
256  if(detStore()->contains<StoredPhysVol>(key))
257  {
258  if(detStore()->retrieve(storedPV,key)==StatusCode::FAILURE)
259  {
260  ATH_MSG_INFO ("Unable to retrieve Stored PV " << key);
261  }
262  }
263 
264  if (!storedPV) {
265  ATH_MSG_INFO ("no Volume for " << key);
266  return false;
267  }
268 
269  const GeoFullPhysVol* larTopVolLink = storedPV->getPhysVol();
270 
271  // get the number of direct childs
272  unsigned int larChilds = larTopVolLink->getNChildVols();
273 
274  for (unsigned int ichild = 0; ichild<larChilds; ++ichild){
275 
276  bool stop = false;
277 
278  // Skip some some of the children : list will get longer in future
279  if ( alvol == CaloSubdetNames::LARCRYO_B &&
280  larTopVolLink->getNameOfChildVol(ichild) == "Total LAR Volume") stop = true;
281  if ( alvol == CaloSubdetNames::SOLENOID &&
282  larTopVolLink->getNameOfChildVol(ichild) == "Total LAR Volume") stop = true;
283 
284  if (!stop) {
285 
286  ATH_MSG_DEBUG ("CaloSubdetNames::ALIGNVOL " << (int) alvol
287  <<" Direct Child " << ichild << " - with name : "
288  << larTopVolLink->getNameOfChildVol(ichild));
289 
290  PVConstLink childVolLink = larTopVolLink->getChildVol(ichild);
291 
292  m_massFractions.clear();
293  m_volumeFractions.clear();
294  m_x0Fractions.clear();
295  m_dEdXFractions.clear();
296  m_averageAFractions.clear();
297  m_averageZFractions.clear();
298 
299  this->blendGeoVolume(childVolLink,
300  m_volumeFractions,
301  m_massFractions,
302  m_x0Fractions,
303  m_dEdXFractions,
304  m_averageAFractions,
305  m_averageZFractions);
306 
307  this->averageFraction(m_volumeFractions,
308  m_massFractions,
309  m_x0Fractions,
310  m_dEdXFractions,
311  m_averageAFractions,
312  m_averageZFractions,
313  childMass, childVolume,
314  childX0, childDeDx,
315  childAverageA, childAverageZ, childRho );
316 
317  m_child_x0Fractions.push_back( childX0 );
318  m_child_dEdXFractions.push_back( childDeDx );
319  m_child_averageAFractions.push_back( childAverageA );
320  m_child_averageZFractions.push_back( childAverageZ );
321  m_child_Rho.push_back( childRho );
322  m_child_Volume.push_back( childVolume );
323  m_child_Mass.push_back( childMass );
324 
325  }
326  }
327 
328  // Now do the final average and store it :
329 
330  ATH_MSG_DEBUG ("Final average for CaloSubdetNames::ALIGNVOL " << (int) alvol
331  << " !!!!! ");
332 
333  this->averageFraction(m_child_Volume, m_child_Mass,
334  m_child_x0Fractions,
335  m_child_dEdXFractions,
336  m_child_averageAFractions,
337  m_child_averageZFractions,
338  childMass, childVolume,
339  childX0, childDeDx,
340  childAverageA, childAverageZ, childRho );
341 
342  mass = childMass/GeoModelKernelUnits::gram;
343  x0 = childX0;
344  dEdx = childDeDx;
345  aveA = childAverageA;
346  aveZ = childAverageZ;
347  }
348 #endif
349 
350  m_alignvol_number.push_back( alvol );
351  m_alignvol_mass.push_back( mass );
352  m_alignvol_x0.push_back( x0 );
353  m_alignvol_dEdX.push_back( dEdx );
354  m_alignvol_aveA.push_back( aveA );
355  m_alignvol_aveZ.push_back( aveZ );
356 
357  return true;
358 }

◆ ScanEMB()

bool LArRecoMaterialTool::ScanEMB ( CaloSubdetNames::ALIGNVOL  alvol,
double &  mass,
double &  x0,
double &  dEdx,
double &  aveA,
double &  aveZ 
)
private

Definition at line 519 of file LArRecoMaterialTool.cxx.

522 {
523  mass = 0.;
524  x0 = 0.;
525  dEdx = 0.;
526  aveA = 0.;
527  aveZ = 0.;
528 
529  if ( alvol != CaloSubdetNames::EMB_POS &&
530  alvol != CaloSubdetNames::EMB_NEG ) {
531  ATH_MSG_INFO ("method ScanEMB does not support CaloSubdetNames::ALIGNVOL " << (int) alvol);
532  return false;
533  }
534 
535  // Shortcut :
536  if (!m_parseGeoModelForMaterial)
537  {
538  mass = 9.40439e+07 ;
539  x0 = 93.9601;
540  dEdx = 6.0789e-24;
541  aveA = 98.5676;
542  aveZ = 39.7406;
543  }
544  else {
545 
546  // Here we go :
547  m_child_massFractions.clear();
548  m_child_volumeFractions.clear();
549  m_child_x0Fractions.clear();
550  m_child_dEdXFractions.clear();
551  m_child_averageAFractions.clear();
552  m_child_averageZFractions.clear();
553  m_child_Rho.clear();
554  m_child_Volume.clear();
555  m_child_Mass.clear();
556 
557  double childVolume = 0.;
558  double childMass = 0.;
559  double childX0 = 0.;
560  double childDeDx = 0.;
561  double childAverageA = 0.;
562  double childAverageZ = 0.;
563  double childRho = 0.;
564 
565  StoredPhysVol* storedPV = nullptr;
566  std::string key = map_av(alvol);
567  if(detStore()->contains<StoredPhysVol>(key))
568  {
569  if(detStore()->retrieve(storedPV,key)==StatusCode::FAILURE)
570  {
571  ATH_MSG_INFO ("Unable to retrieve Stored PV " << key);
572  }
573  }
574 
575  if (!storedPV) {
576  ATH_MSG_INFO ("no Volume for " << key);
577  mass = 0.;
578  x0 = 0.;
579  dEdx = 0.;
580  aveA = 0.;
581  aveZ = 0.;
582  return false;
583  }
584 
585  const GeoFullPhysVol* larTopVolLink = storedPV->getPhysVol();
586 
587  // get the number of direct childs
588  unsigned int larChilds = larTopVolLink->getNChildVols();
589 
590  for (unsigned int ichild = 0; ichild<larChilds; ++ichild){
591 
592  bool stop = false;
593 
594  // Skip some some of the children : list will get longer in future
595  //if ( alvol == CaloSubdetNames::LARCRYO_B &&
596  // larTopVolLink->getNameOfChildVol(ichild) == "Total LAR Volume") stop = true;
597 
598  if (!stop) {
599 
600  ATH_MSG_DEBUG ("CaloSubdetNames::ALIGNVOL " << (int) alvol
601  <<" Direct Child " << ichild << " - with name : "
602  << larTopVolLink->getNameOfChildVol(ichild));
603 
604  PVConstLink childVolLink = larTopVolLink->getChildVol(ichild);
605 
606  m_massFractions.clear();
607  m_volumeFractions.clear();
608  m_x0Fractions.clear();
609  m_dEdXFractions.clear();
610  m_averageAFractions.clear();
611  m_averageZFractions.clear();
612 
613  this->blendGeoVolume(childVolLink,
614  m_volumeFractions,
615  m_massFractions,
616  m_x0Fractions,
617  m_dEdXFractions,
618  m_averageAFractions,
619  m_averageZFractions);
620 
621  this->averageFraction(m_volumeFractions,
622  m_massFractions,
623  m_x0Fractions,
624  m_dEdXFractions,
625  m_averageAFractions,
626  m_averageZFractions,
627  childMass, childVolume,
628  childX0, childDeDx,
629  childAverageA, childAverageZ, childRho );
630 
631  m_child_x0Fractions.push_back( childX0 );
632  m_child_dEdXFractions.push_back( childDeDx );
633  m_child_averageAFractions.push_back( childAverageA );
634  m_child_averageZFractions.push_back( childAverageZ );
635  m_child_Rho.push_back( childRho );
636  m_child_Volume.push_back( childVolume );
637  m_child_Mass.push_back( childMass );
638 
639  }
640  }
641 
642  // Now do the final average and store it :
643 
644  ATH_MSG_DEBUG ("Final average for CaloSubdetNames::ALIGNVOL " << (int) alvol
645  << " !!!!! ");
646 
647  this->averageFraction(m_child_Volume, m_child_Mass,
648  m_child_x0Fractions,
649  m_child_dEdXFractions,
650  m_child_averageAFractions,
651  m_child_averageZFractions,
652  childMass, childVolume,
653  childX0, childDeDx,
654  childAverageA, childAverageZ, childRho );
655 
656  mass = childMass*(1./GeoModelKernelUnits::gram);
657  x0 = childX0;
658  dEdx = childDeDx;
659  aveA = childAverageA;
660  aveZ = childAverageZ;
661  }
662 
663  m_alignvol_number.push_back( alvol );
664  m_alignvol_mass.push_back( mass );
665  m_alignvol_x0.push_back( x0 );
666  m_alignvol_dEdX.push_back( dEdx );
667  m_alignvol_aveA.push_back( aveA );
668  m_alignvol_aveZ.push_back( aveZ );
669 
670  return true;
671 }

◆ ScanEMEC()

bool LArRecoMaterialTool::ScanEMEC ( CaloSubdetNames::ALIGNVOL  alvol,
double &  mass,
double &  x0,
double &  dEdx,
double &  aveA,
double &  aveZ 
)
private

Definition at line 674 of file LArRecoMaterialTool.cxx.

677 {
678  mass = 0.;
679  x0 = 0.;
680  dEdx = 0.;
681  aveA = 0.;
682  aveZ = 0.;
683 
684  if ( alvol != CaloSubdetNames::EMEC_POS && alvol != CaloSubdetNames::EMEC_NEG ) {
685  ATH_MSG_INFO ("method ScanEMEM does not support CaloSubdetNames::ALIGNVOL " << (int) alvol);
686  return false;
687  }
688 
689  // EMEC use a couple of custom volumes :
690  // LAr::EMEC::OuterWheel::Absorber has volume=0 mm**3 and 5 elements
691  // LAr::EMEC::OuterWheel::Electrode has volume=0 mm**3 and 5 elements
692  // ==> with volume = 0 avarages cannot work !!!
693  // Take EMB's density for now ....
694  mass = 9.40439e+07 ;
695  x0 = 93.9601;
696  dEdx = 6.0789e-24;
697  aveA = 98.5676;
698  aveZ = 39.7406;
699 
700 
701  m_alignvol_number.push_back( alvol );
702  m_alignvol_mass.push_back( mass );
703  m_alignvol_x0.push_back( x0 );
704  m_alignvol_dEdX.push_back( dEdx );
705  m_alignvol_aveA.push_back( aveA );
706  m_alignvol_aveZ.push_back( aveZ );
707 
708  return true;
709 }

◆ ScanFCAL()

bool LArRecoMaterialTool::ScanFCAL ( CaloSubdetNames::ALIGNVOL  alvol,
double &  mass,
double &  x0,
double &  dEdx,
double &  aveA,
double &  aveZ 
)
private

Definition at line 870 of file LArRecoMaterialTool.cxx.

873 {
874  mass = 0.;
875  x0 = 0.;
876  dEdx = 0.;
877  aveA = 0.;
878  aveZ = 0.;
879 
880  if ( alvol != CaloSubdetNames::FCAL1_POS &&
881  alvol != CaloSubdetNames::FCAL1_NEG &&
882  alvol != CaloSubdetNames::FCAL2_POS &&
883  alvol != CaloSubdetNames::FCAL2_NEG &&
884  alvol != CaloSubdetNames::FCAL3_POS &&
885  alvol != CaloSubdetNames::FCAL3_NEG ) {
886  ATH_MSG_INFO ("method ScanEMB does not support CaloSubdetNames::ALIGNVOL "
887  << (int) alvol);
888  return false;
889  }
890 
891  // Shortcut :
892  if (!m_parseGeoModelForMaterial)
893  {
894  if (alvol == CaloSubdetNames::FCAL1_POS || alvol == CaloSubdetNames::FCAL1_NEG) {
895  mass = 47825.7;
896  x0 = 127.627;
897  dEdx = 4.36158e-24;
898  aveA = 42.5374;
899  aveZ = 19.2601;
900  }
901  else if (alvol == CaloSubdetNames::FCAL2_POS || alvol == CaloSubdetNames::FCAL2_NEG) {
902  mass = 1.77759e+06;
903  x0 = 81.6747;
904  dEdx = 1.9157e-23;
905  aveA = 159.937;
906  aveZ = 64.5724;
907  }
908  else if (alvol == CaloSubdetNames::FCAL3_POS || alvol == CaloSubdetNames::FCAL3_NEG) {
909  mass = 1.80016e+06;
910  x0 = 82.9447;
911  dEdx = 1.86512e-23;
912  aveA = 158.884;
913  aveZ = 64.1669;
914  }
915 
916  }
917 
918  else {
919 
920  // Here we go :
921  m_child_massFractions.clear();
922  m_child_volumeFractions.clear();
923  m_child_x0Fractions.clear();
924  m_child_dEdXFractions.clear();
925  m_child_averageAFractions.clear();
926  m_child_averageZFractions.clear();
927  m_child_Rho.clear();
928  m_child_Volume.clear();
929  m_child_Mass.clear();
930 
931  double childVolume = 0.;
932  double childMass = 0.;
933  double childX0 = 0.;
934  double childDeDx = 0.;
935  double childAverageA = 0.;
936  double childAverageZ = 0.;
937  double childRho = 0.;
938 
939  StoredPhysVol* storedPV = nullptr;
940  std::string key = map_av(alvol);
941  if(detStore()->contains<StoredPhysVol>(key))
942  {
943  if(detStore()->retrieve(storedPV,key)==StatusCode::FAILURE)
944  {
945  ATH_MSG_INFO ("Unable to retrieve Stored PV " << key);
946  }
947  }
948 
949  if (!storedPV) {
950  ATH_MSG_INFO ("no Volume for " << key);
951  mass = 0.;
952  x0 = 0.;
953  dEdx = 0.;
954  aveA = 0.;
955  aveZ = 0.;
956  return false;
957  }
958 
959  const GeoFullPhysVol* larTopVolLink = storedPV->getPhysVol();
960 
961  // get the number of direct childs
962  unsigned int larChilds = larTopVolLink->getNChildVols();
963 
964  for (unsigned int ichild = 0; ichild<larChilds; ++ichild){
965 
966  bool stop = false;
967 
968  // Skip some some of the children : list will get longer in future
969  //if ( alvol == CaloSubdetNames::LARCRYO_B &&
970  // larTopVolLink->getNameOfChildVol(ichild) == "Total LAR Volume") stop = true;
971 
972  if (!stop) {
973 
974  ATH_MSG_DEBUG ("CaloSubdetNames::ALIGNVOL " << (int) alvol
975  <<" Direct Child " << ichild << " - with name : "
976  << larTopVolLink->getNameOfChildVol(ichild));
977 
978  PVConstLink childVolLink = larTopVolLink->getChildVol(ichild);
979 
980  m_massFractions.clear();
981  m_volumeFractions.clear();
982  m_x0Fractions.clear();
983  m_dEdXFractions.clear();
984  m_averageAFractions.clear();
985  m_averageZFractions.clear();
986 
987  this->blendGeoVolume(childVolLink,
988  m_volumeFractions,
989  m_massFractions,
990  m_x0Fractions,
991  m_dEdXFractions,
992  m_averageAFractions,
993  m_averageZFractions);
994 
995  this->averageFraction(m_volumeFractions,
996  m_massFractions,
997  m_x0Fractions,
998  m_dEdXFractions,
999  m_averageAFractions,
1000  m_averageZFractions,
1001  childMass, childVolume,
1002  childX0, childDeDx,
1003  childAverageA, childAverageZ, childRho );
1004 
1005  m_child_x0Fractions.push_back( childX0 );
1006  m_child_dEdXFractions.push_back( childDeDx );
1007  m_child_averageAFractions.push_back( childAverageA );
1008  m_child_averageZFractions.push_back( childAverageZ );
1009  m_child_Rho.push_back( childRho );
1010  m_child_Volume.push_back( childVolume );
1011  m_child_Mass.push_back( childMass );
1012 
1013  }
1014  }
1015 
1016  // Now do the final average and store it :
1017 
1018  ATH_MSG_DEBUG ("Final average for CaloSubdetNames::ALIGNVOL " << (int) alvol
1019  << " !!!!! ");
1020 
1021  this->averageFraction(m_child_Volume, m_child_Mass,
1022  m_child_x0Fractions,
1023  m_child_dEdXFractions,
1024  m_child_averageAFractions,
1025  m_child_averageZFractions,
1026  childMass, childVolume,
1027  childX0, childDeDx,
1028  childAverageA, childAverageZ, childRho );
1029 
1030  mass = childMass*(1./GeoModelKernelUnits::gram);
1031  x0 = childX0;
1032  dEdx = childDeDx;
1033  aveA = childAverageA;
1034  aveZ = childAverageZ;
1035  }
1036 
1037  m_alignvol_number.push_back( alvol );
1038  m_alignvol_mass.push_back( mass );
1039  m_alignvol_x0.push_back( x0 );
1040  m_alignvol_dEdX.push_back( dEdx );
1041  m_alignvol_aveA.push_back( aveA );
1042  m_alignvol_aveZ.push_back( aveZ );
1043 
1044  return true;
1045 }

◆ ScanHEC()

bool LArRecoMaterialTool::ScanHEC ( CaloSubdetNames::ALIGNVOL  alvol,
double &  mass,
double &  x0,
double &  dEdx,
double &  aveA,
double &  aveZ 
)
private

Definition at line 712 of file LArRecoMaterialTool.cxx.

715 {
716  mass = 0.;
717  x0 = 0.;
718  dEdx = 0.;
719  aveA = 0.;
720  aveZ = 0.;
721 
722  if ( alvol != CaloSubdetNames::HEC1_POS &&
723  alvol != CaloSubdetNames::HEC1_NEG &&
724  alvol != CaloSubdetNames::HEC2_POS &&
725  alvol != CaloSubdetNames::HEC2_NEG ) {
726  ATH_MSG_INFO ("method ScanEMB does not support CaloSubdetNames::ALIGNVOL "
727  << (int) alvol);
728  return false;
729  }
730 
731  // Shortcut :
732  if (!m_parseGeoModelForMaterial)
733  {
734  mass = 6.06935e+06;
735  x0 = 91.1845;
736  dEdx = 1.14141e-23;
737  aveA = 58.9808;
738  aveZ = 26.8777;
739  }
740  else {
741 
742  // Here we go :
743  m_child_massFractions.clear();
744  m_child_volumeFractions.clear();
745  m_child_x0Fractions.clear();
746  m_child_dEdXFractions.clear();
747  m_child_averageAFractions.clear();
748  m_child_averageZFractions.clear();
749  m_child_Rho.clear();
750  m_child_Volume.clear();
751  m_child_Mass.clear();
752 
753  double childVolume = 0.;
754  double childMass = 0.;
755  double childX0 = 0.;
756  double childDeDx = 0.;
757  double childAverageA = 0.;
758  double childAverageZ = 0.;
759  double childRho = 0.;
760 
761  StoredPhysVol* storedPV = nullptr;
762  std::string key = map_av(alvol);
763  if(detStore()->contains<StoredPhysVol>(key))
764  {
765  if(detStore()->retrieve(storedPV,key)==StatusCode::FAILURE)
766  {
767  ATH_MSG_INFO ("Unable to retrieve Stored PV " << key);
768  }
769  }
770 
771  if (!storedPV) {
772  ATH_MSG_INFO ("no Volume for " << key);
773  mass = 0.;
774  x0 = 0.;
775  dEdx = 0.;
776  aveA = 0.;
777  aveZ = 0.;
778  return false;
779  }
780 
781  const GeoFullPhysVol* larTopVolLink = storedPV->getPhysVol();
782 
783  // get the number of direct childs
784  unsigned int larChilds = larTopVolLink->getNChildVols();
785 
786  for (unsigned int ichild = 0; ichild<larChilds; ++ichild){
787 
788  bool stop = false;
789 
790  // Skip some some of the children : list will get longer in future
791  //if ( alvol == CaloSubdetNames::LARCRYO_B &&
792  // larTopVolLink->getNameOfChildVol(ichild) == "Total LAR Volume") stop = true;
793 
794  if (!stop) {
795 
796  ATH_MSG_DEBUG ("CaloSubdetNames::ALIGNVOL " << (int) alvol
797  <<" Direct Child " << ichild << " - with name : "
798  << larTopVolLink->getNameOfChildVol(ichild));
799 
800  PVConstLink childVolLink = larTopVolLink->getChildVol(ichild);
801 
802  m_massFractions.clear();
803  m_volumeFractions.clear();
804  m_x0Fractions.clear();
805  m_dEdXFractions.clear();
806  m_averageAFractions.clear();
807  m_averageZFractions.clear();
808 
809  this->blendGeoVolume(childVolLink,
810  m_volumeFractions,
811  m_massFractions,
812  m_x0Fractions,
813  m_dEdXFractions,
814  m_averageAFractions,
815  m_averageZFractions);
816 
817  this->averageFraction(m_volumeFractions,
818  m_massFractions,
819  m_x0Fractions,
820  m_dEdXFractions,
821  m_averageAFractions,
822  m_averageZFractions,
823  childMass, childVolume,
824  childX0, childDeDx,
825  childAverageA, childAverageZ, childRho );
826 
827  m_child_x0Fractions.push_back( childX0 );
828  m_child_dEdXFractions.push_back( childDeDx );
829  m_child_averageAFractions.push_back( childAverageA );
830  m_child_averageZFractions.push_back( childAverageZ );
831  m_child_Rho.push_back( childRho );
832  m_child_Volume.push_back( childVolume );
833  m_child_Mass.push_back( childMass );
834 
835  }
836  }
837 
838  // Now do the final average and store it :
839 
840  ATH_MSG_DEBUG ("Final average for CaloSubdetNames::ALIGNVOL " << (int) alvol
841  << " !!!!! ");
842 
843  this->averageFraction(m_child_Volume, m_child_Mass,
844  m_child_x0Fractions,
845  m_child_dEdXFractions,
846  m_child_averageAFractions,
847  m_child_averageZFractions,
848  childMass, childVolume,
849  childX0, childDeDx,
850  childAverageA, childAverageZ, childRho );
851 
852  mass = childMass*(1./GeoModelKernelUnits::gram);
853  x0 = childX0;
854  dEdx = childDeDx;
855  aveA = childAverageA;
856  aveZ = childAverageZ;
857  }
858 
859  m_alignvol_number.push_back( alvol );
860  m_alignvol_mass.push_back( mass );
861  m_alignvol_x0.push_back( x0 );
862  m_alignvol_dEdX.push_back( dEdx );
863  m_alignvol_aveA.push_back( aveA );
864  m_alignvol_aveZ.push_back( aveZ );
865 
866  return true;
867 }

◆ ScanPS()

bool LArRecoMaterialTool::ScanPS ( CaloSubdetNames::ALIGNVOL  alvol,
double &  mass,
double &  x0,
double &  dEdx,
double &  aveA,
double &  aveZ 
)
private

Definition at line 361 of file LArRecoMaterialTool.cxx.

364 {
365  mass = 0.;
366  x0 = 0.;
367  dEdx = 0.;
368  aveA = 0.;
369  aveZ = 0.;
370 
371  if ( alvol != CaloSubdetNames::PRESAMPLER_B_POS
372  && alvol != CaloSubdetNames::PRESAMPLER_B_NEG
373  && alvol != CaloSubdetNames::PRESAMPLER_EC_POS
374  && alvol != CaloSubdetNames::PRESAMPLER_EC_NEG
375  ) {
376  ATH_MSG_INFO ("method ScanCryo does not support CaloSubdetNames::ALIGNVOL " << (int) alvol);
377  return false;
378  }
379 
380  // for EMEC take Barrel numbers
381 
382  // Shortcut :
383  if (!m_parseGeoModelForMaterial ||
384  alvol == CaloSubdetNames::PRESAMPLER_EC_POS ||
385  alvol == CaloSubdetNames::PRESAMPLER_EC_NEG)
386  {
387  mass = 1.80785e+06 ;
388  x0 = 141.837;
389  dEdx = 3.55847e-24;
390  aveA = 39.5576;
391  aveZ = 17.8274;
392  }
393 
394  else {
395 
396  // Go for it :
397  m_child_massFractions.clear();
398  m_child_volumeFractions.clear();
399  m_child_x0Fractions.clear();
400  m_child_dEdXFractions.clear();
401  m_child_averageAFractions.clear();
402  m_child_averageZFractions.clear();
403  m_child_Rho.clear();
404  m_child_Volume.clear();
405  m_child_Mass.clear();
406 
407  double childVolume = 0.;
408  double childMass = 0.;
409  double childX0 = 0.;
410  double childDeDx = 0.;
411  double childAverageA = 0.;
412  double childAverageZ = 0.;
413  double childRho = 0.;
414 
415  StoredPhysVol* storedPV = nullptr;
416  std::string key = map_av(alvol);
417  if(detStore()->contains<StoredPhysVol>(key))
418  {
419  if(detStore()->retrieve(storedPV,key)==StatusCode::FAILURE)
420  {
421  ATH_MSG_INFO ("Unable to retrieve Stored PV " << key);
422  }
423  }
424 
425  if (!storedPV) {
426  ATH_MSG_INFO ("no Volume for " << key);
427  return false;
428  }
429 
430  const GeoFullPhysVol* larTopVolLink = storedPV->getPhysVol();
431 
432  // get the number of direct childs
433  unsigned int larChilds = larTopVolLink->getNChildVols();
434 
435  for (unsigned int ichild = 0; ichild<larChilds; ++ichild){
436 
437  bool stop = false;
438 
439  // Skip some some of the children : list will get longer in future
440  //if ( alvol == CaloSubdetNames::LARCRYO_B &&
441  // larTopVolLink->getNameOfChildVol(ichild) == "Total LAR Volume") stop = true;
442 
443  if (!stop) {
444 
445  ATH_MSG_DEBUG ("CaloSubdetNames::ALIGNVOL " << (int) alvol
446  <<" Direct Child " << ichild << " - with name : "
447  << larTopVolLink->getNameOfChildVol(ichild));
448 
449  PVConstLink childVolLink = larTopVolLink->getChildVol(ichild);
450 
451  m_massFractions.clear();
452  m_volumeFractions.clear();
453  m_x0Fractions.clear();
454  m_dEdXFractions.clear();
455  m_averageAFractions.clear();
456  m_averageZFractions.clear();
457 
458  this->blendGeoVolume(childVolLink,
459  m_volumeFractions,
460  m_massFractions,
461  m_x0Fractions,
462  m_dEdXFractions,
463  m_averageAFractions,
464  m_averageZFractions);
465 
466  this->averageFraction(m_volumeFractions,
467  m_massFractions,
468  m_x0Fractions,
469  m_dEdXFractions,
470  m_averageAFractions,
471  m_averageZFractions,
472  childMass, childVolume,
473  childX0, childDeDx,
474  childAverageA, childAverageZ, childRho );
475 
476  m_child_x0Fractions.push_back( childX0 );
477  m_child_dEdXFractions.push_back( childDeDx );
478  m_child_averageAFractions.push_back( childAverageA );
479  m_child_averageZFractions.push_back( childAverageZ );
480  m_child_Rho.push_back( childRho );
481  m_child_Volume.push_back( childVolume );
482  m_child_Mass.push_back( childMass );
483 
484  }
485  }
486 
487  // Now do the final average and store it :
488 
489  ATH_MSG_DEBUG ("Final average for CaloSubdetNames::ALIGNVOL " << (int) alvol
490  << " !!!!! ");
491 
492  this->averageFraction(m_child_Volume, m_child_Mass,
493  m_child_x0Fractions,
494  m_child_dEdXFractions,
495  m_child_averageAFractions,
496  m_child_averageZFractions,
497  childMass, childVolume,
498  childX0, childDeDx,
499  childAverageA, childAverageZ, childRho );
500 
501  mass = childMass*(1./GeoModelKernelUnits::gram);
502  x0 = childX0;
503  dEdx = childDeDx;
504  aveA = childAverageA;
505  aveZ = childAverageZ;
506  }
507 
508  m_alignvol_number.push_back( alvol );
509  m_alignvol_mass.push_back( mass );
510  m_alignvol_x0.push_back( x0 );
511  m_alignvol_dEdX.push_back( dEdx );
512  m_alignvol_aveA.push_back( aveA );
513  m_alignvol_aveZ.push_back( aveZ );
514 
515  return true;
516 }

◆ ScanTree()

void LArRecoMaterialTool::ScanTree ( )
private

Private method called at initialisation : fills all arrays by calling the "optimised" methods.

Definition at line 143 of file LArRecoMaterialTool.cxx.

144 {
145  double mass = 0.;
146  double x0 = 0.;
147  double dEdx = 0.;
148  double aveA = 0.;
149  double aveZ = 0.;
150 
151  bool result = false;
152 
153  ATH_MSG_INFO (" Entering method ScanTree() ");
154 
155  result = this->ScanCRYO(CaloSubdetNames::LARCRYO_B,mass,x0,dEdx,aveA,aveZ);
156  result = this->ScanCRYO(CaloSubdetNames::LARCRYO_EC_POS,mass,x0,dEdx,aveA,aveZ);
157  result = this->ScanCRYO(CaloSubdetNames::LARCRYO_EC_NEG,mass,x0,dEdx,aveA,aveZ);
158  result = this->ScanCRYO(CaloSubdetNames::SOLENOID,mass,x0,dEdx,aveA,aveZ);
159 
160  result = this->ScanPS(CaloSubdetNames::PRESAMPLER_B_POS,mass,x0,dEdx,aveA,aveZ);
161  result = this->ScanPS(CaloSubdetNames::PRESAMPLER_B_NEG,mass,x0,dEdx,aveA,aveZ);
162 
163  result = this->ScanEMB(CaloSubdetNames::EMB_POS,mass,x0,dEdx,aveA,aveZ);
164  result = this->ScanEMB(CaloSubdetNames::EMB_NEG,mass,x0,dEdx,aveA,aveZ);
165 
166  result = this->ScanPS(CaloSubdetNames::PRESAMPLER_EC_POS,mass,x0,dEdx,aveA,aveZ);
167  result = this->ScanPS(CaloSubdetNames::PRESAMPLER_EC_NEG,mass,x0,dEdx,aveA,aveZ);
168 
169  result = this->ScanEMEC(CaloSubdetNames::EMEC_POS,mass,x0,dEdx,aveA,aveZ);
170  result = this->ScanEMEC(CaloSubdetNames::EMEC_NEG,mass,x0,dEdx,aveA,aveZ);
171 
172  result = this->ScanHEC(CaloSubdetNames::HEC1_POS,mass,x0,dEdx,aveA,aveZ);
173  result = this->ScanHEC(CaloSubdetNames::HEC1_NEG,mass,x0,dEdx,aveA,aveZ);
174  result = this->ScanHEC(CaloSubdetNames::HEC2_POS,mass,x0,dEdx,aveA,aveZ);
175  result = this->ScanHEC(CaloSubdetNames::HEC2_NEG,mass,x0,dEdx,aveA,aveZ);
176 
177  result = this->ScanFCAL(CaloSubdetNames::FCAL1_POS,mass,x0,dEdx,aveA,aveZ);
178  result = this->ScanFCAL(CaloSubdetNames::FCAL1_NEG,mass,x0,dEdx,aveA,aveZ);
179  result = this->ScanFCAL(CaloSubdetNames::FCAL2_POS,mass,x0,dEdx,aveA,aveZ);
180  result = this->ScanFCAL(CaloSubdetNames::FCAL2_NEG,mass,x0,dEdx,aveA,aveZ);
181  result = this->ScanFCAL(CaloSubdetNames::FCAL3_POS,mass,x0,dEdx,aveA,aveZ);
182  result = this->ScanFCAL(CaloSubdetNames::FCAL3_NEG,mass,x0,dEdx,aveA,aveZ);
183 
184  if (!result)
185  ATH_MSG_INFO (" If seems to have failed ?");
186 
187  return;
188 }

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

Member Data Documentation

◆ m_alignvol_aveA

std::vector<double> LArRecoMaterialTool::m_alignvol_aveA
private

Definition at line 152 of file LArRecoMaterialTool.h.

◆ m_alignvol_aveZ

std::vector<double> LArRecoMaterialTool::m_alignvol_aveZ
private

Definition at line 153 of file LArRecoMaterialTool.h.

◆ m_alignvol_dEdX

std::vector<double> LArRecoMaterialTool::m_alignvol_dEdX
private

Definition at line 151 of file LArRecoMaterialTool.h.

◆ m_alignvol_mass

std::vector<double> LArRecoMaterialTool::m_alignvol_mass
private

Definition at line 149 of file LArRecoMaterialTool.h.

◆ m_alignvol_number

std::vector<CaloSubdetNames::ALIGNVOL> LArRecoMaterialTool::m_alignvol_number
private

Definition at line 148 of file LArRecoMaterialTool.h.

◆ m_alignvol_x0

std::vector<double> LArRecoMaterialTool::m_alignvol_x0
private

Definition at line 150 of file LArRecoMaterialTool.h.

◆ m_averageAFractions

std::vector<double> LArRecoMaterialTool::m_averageAFractions
private

Definition at line 160 of file LArRecoMaterialTool.h.

◆ m_averageZFractions

std::vector<double> LArRecoMaterialTool::m_averageZFractions
private

Definition at line 161 of file LArRecoMaterialTool.h.

◆ m_calo_id

const CaloCell_ID* LArRecoMaterialTool::m_calo_id
private

Definition at line 144 of file LArRecoMaterialTool.h.

◆ m_child_averageAFractions

std::vector<double> LArRecoMaterialTool::m_child_averageAFractions
private

Definition at line 167 of file LArRecoMaterialTool.h.

◆ m_child_averageZFractions

std::vector<double> LArRecoMaterialTool::m_child_averageZFractions
private

Definition at line 168 of file LArRecoMaterialTool.h.

◆ m_child_dEdXFractions

std::vector<double> LArRecoMaterialTool::m_child_dEdXFractions
private

Definition at line 166 of file LArRecoMaterialTool.h.

◆ m_child_Mass

std::vector<double> LArRecoMaterialTool::m_child_Mass
private

Definition at line 171 of file LArRecoMaterialTool.h.

◆ m_child_massFractions

std::vector<double> LArRecoMaterialTool::m_child_massFractions
private

Definition at line 163 of file LArRecoMaterialTool.h.

◆ m_child_Rho

std::vector<double> LArRecoMaterialTool::m_child_Rho
private

Definition at line 169 of file LArRecoMaterialTool.h.

◆ m_child_Volume

std::vector<double> LArRecoMaterialTool::m_child_Volume
private

Definition at line 170 of file LArRecoMaterialTool.h.

◆ m_child_volumeFractions

std::vector<double> LArRecoMaterialTool::m_child_volumeFractions
private

Definition at line 164 of file LArRecoMaterialTool.h.

◆ m_child_x0Fractions

std::vector<double> LArRecoMaterialTool::m_child_x0Fractions
private

Definition at line 165 of file LArRecoMaterialTool.h.

◆ m_dEdXFractions

std::vector<double> LArRecoMaterialTool::m_dEdXFractions
private

Definition at line 159 of file LArRecoMaterialTool.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_dm_id

const CaloDM_ID* LArRecoMaterialTool::m_dm_id
private

Definition at line 145 of file LArRecoMaterialTool.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_massFractions

std::vector<double> LArRecoMaterialTool::m_massFractions
private

Definition at line 156 of file LArRecoMaterialTool.h.

◆ m_parseGeoModelForMaterial

bool LArRecoMaterialTool::m_parseGeoModelForMaterial = false
private

Private bool set at initialisation : decides if the GeoModel is parsed or bypassed.

Definition at line 79 of file LArRecoMaterialTool.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_volumeFractions

std::vector<double> LArRecoMaterialTool::m_volumeFractions
private

Definition at line 157 of file LArRecoMaterialTool.h.

◆ m_x0Fractions

std::vector<double> LArRecoMaterialTool::m_x0Fractions
private

Definition at line 158 of file LArRecoMaterialTool.h.

The documentation for this class was generated from the following files:
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
LArRecoMaterialTool::m_child_averageZFractions
std::vector< double > m_child_averageZFractions
Definition: LArRecoMaterialTool.h:168
LArRecoMaterialTool::ScanEMEC
bool ScanEMEC(CaloSubdetNames::ALIGNVOL alvol, double &mass, double &x0, double &dEdx, double &aveA, double &aveZ)
Definition: LArRecoMaterialTool.cxx:674
CaloSubdetNames::HEC2_POS
@ HEC2_POS
Definition: CaloSubdetNames.h:42
get_generator_info.result
result
Definition: get_generator_info.py:21
python.PerfMonSerializer.p
def p
Definition: PerfMonSerializer.py:743
LArRecoMaterialTool::m_alignvol_mass
std::vector< double > m_alignvol_mass
Definition: LArRecoMaterialTool.h:149
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
LArRecoMaterialTool::m_child_Mass
std::vector< double > m_child_Mass
Definition: LArRecoMaterialTool.h:171
AthCommonDataStore::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
CaloSubdetNames::EMB_POS
@ EMB_POS
Definition: CaloSubdetNames.h:38
LArRecoMaterialTool::m_dm_id
const CaloDM_ID * m_dm_id
Definition: LArRecoMaterialTool.h:145
CaloSubdetNames::PRESAMPLER_B_NEG
@ PRESAMPLER_B_NEG
Definition: CaloSubdetNames.h:37
CaloSubdetNames::FCAL2_POS
@ FCAL2_POS
Definition: CaloSubdetNames.h:44
LArRecoMaterialTool::m_child_Volume
std::vector< double > m_child_Volume
Definition: LArRecoMaterialTool.h:170
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
PixelModuleFeMask_create_db.stop
int stop
Definition: PixelModuleFeMask_create_db.py:76
python.SystemOfUnits.gram
int gram
Definition: SystemOfUnits.py:165
LArRecoMaterialTool::m_averageZFractions
std::vector< double > m_averageZFractions
Definition: LArRecoMaterialTool.h:161
LArRecoMaterialTool::m_child_volumeFractions
std::vector< double > m_child_volumeFractions
Definition: LArRecoMaterialTool.h:164
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
LArRecoMaterialTool::m_calo_id
const CaloCell_ID * m_calo_id
Definition: LArRecoMaterialTool.h:144
LArRecoMaterialTool::ScanPS
bool ScanPS(CaloSubdetNames::ALIGNVOL alvol, double &mass, double &x0, double &dEdx, double &aveA, double &aveZ)
Definition: LArRecoMaterialTool.cxx:361
CaloSubdetNames::HEC1_NEG
@ HEC1_NEG
Definition: CaloSubdetNames.h:41
dqt_zlumi_pandas.mass
mass
Definition: dqt_zlumi_pandas.py:170
CaloSubdetNames::FCAL1_POS
@ FCAL1_POS
Definition: CaloSubdetNames.h:43
LArRecoMaterialTool::ScanTree
void ScanTree()
Private method called at initialisation : fills all arrays by calling the "optimised" methods.
Definition: LArRecoMaterialTool.cxx:143
CaloSubdetNames::EMEC_POS
@ EMEC_POS
Definition: CaloSubdetNames.h:40
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
LArRecoMaterialTool::m_alignvol_aveZ
std::vector< double > m_alignvol_aveZ
Definition: LArRecoMaterialTool.h:153
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
CaloSubdetNames::HEC2_NEG
@ HEC2_NEG
Definition: CaloSubdetNames.h:42
CaloSubdetNames::FCAL3_POS
@ FCAL3_POS
Definition: CaloSubdetNames.h:45
StoredPhysVol
Definition: StoredPhysVol.h:27
CaloSubdetNames::PRESAMPLER_EC_POS
@ PRESAMPLER_EC_POS
Definition: CaloSubdetNames.h:39
CaloSubdetNames::PRESAMPLER_B_POS
@ PRESAMPLER_B_POS
Definition: CaloSubdetNames.h:37
CaloSubdetNames::EMEC_NEG
@ EMEC_NEG
Definition: CaloSubdetNames.h:40
StoredPhysVol::getPhysVol
GeoFullPhysVol * getPhysVol()
Destructor.
Definition: StoredPhysVol.cxx:20
LArRecoMaterialTool::ScanEMB
bool ScanEMB(CaloSubdetNames::ALIGNVOL alvol, double &mass, double &x0, double &dEdx, double &aveA, double &aveZ)
Definition: LArRecoMaterialTool.cxx:519
AthCommonDataStore
Definition: AthCommonDataStore.h:52
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
CaloSubdetNames::LARCRYO_B
@ LARCRYO_B
Definition: CaloSubdetNames.h:36
lumiFormat.i
int i
Definition: lumiFormat.py:92
LArRecoMaterialTool::m_child_averageAFractions
std::vector< double > m_child_averageAFractions
Definition: LArRecoMaterialTool.h:167
LArRecoMaterialTool::m_child_x0Fractions
std::vector< double > m_child_x0Fractions
Definition: LArRecoMaterialTool.h:165
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
LArRecoMaterialTool::m_alignvol_aveA
std::vector< double > m_alignvol_aveA
Definition: LArRecoMaterialTool.h:152
CaloSubdetNames::LARCRYO_EC_POS
@ LARCRYO_EC_POS
Definition: CaloSubdetNames.h:36
CaloSubdetNames::FCAL3_NEG
@ FCAL3_NEG
Definition: CaloSubdetNames.h:45
test_pyathena.parent
parent
Definition: test_pyathena.py:15
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
LArRecoMaterialTool::ScanFCAL
bool ScanFCAL(CaloSubdetNames::ALIGNVOL alvol, double &mass, double &x0, double &dEdx, double &aveA, double &aveZ)
Definition: LArRecoMaterialTool.cxx:870
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
AthAlgTool::AthAlgTool
AthAlgTool()
Default constructor:
LArRecoMaterialTool::m_child_Rho
std::vector< double > m_child_Rho
Definition: LArRecoMaterialTool.h:169
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
CaloSubdetNames::HEC1_POS
@ HEC1_POS
Definition: CaloSubdetNames.h:41
LArRecoMaterialTool::m_dEdXFractions
std::vector< double > m_dEdXFractions
Definition: LArRecoMaterialTool.h:159
LArRecoMaterialTool::m_massFractions
std::vector< double > m_massFractions
Definition: LArRecoMaterialTool.h:156
CaloSubdetNames::PRESAMPLER_EC_NEG
@ PRESAMPLER_EC_NEG
Definition: CaloSubdetNames.h:39
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:192
CaloSubdetNames::EMB_NEG
@ EMB_NEG
Definition: CaloSubdetNames.h:38
CaloSubdetNames::FCAL2_NEG
@ FCAL2_NEG
Definition: CaloSubdetNames.h:44
CaloSubdetNames::SOLENOID
@ SOLENOID
Definition: CaloSubdetNames.h:46
LArRecoMaterialTool::ScanCRYO
bool ScanCRYO(CaloSubdetNames::ALIGNVOL alvol, double &mass, double &x0, double &dEdx, double &aveA, double &aveZ)
Private optimised methods for each subsystem.
Definition: LArRecoMaterialTool.cxx:191
LArRecoMaterialTool::m_child_massFractions
std::vector< double > m_child_massFractions
Definition: LArRecoMaterialTool.h:163
LArRecoMaterialTool::m_averageAFractions
std::vector< double > m_averageAFractions
Definition: LArRecoMaterialTool.h:160
LArRecoMaterialTool::m_alignvol_dEdX
std::vector< double > m_alignvol_dEdX
Definition: LArRecoMaterialTool.h:151
LArRecoMaterialTool::blendGeoVolume
void blendGeoVolume(GeoPVConstLink &geoVolLink, std::vector< double > &volume, std::vector< double > &mass, std::vector< double > &x0, std::vector< double > &dEdX, std::vector< double > &aveA, std::vector< double > &aveZ) const
Private method to blend GeoLogVolume.
Definition: LArRecoMaterialTool.cxx:1048
LArRecoMaterialTool::map_av
std::string map_av(CaloSubdetNames::ALIGNVOL alvol)
Definition: LArRecoMaterialTool.cxx:1161
a
TList * a
Definition: liststreamerinfos.cxx:10
LArRecoMaterialTool::m_parseGeoModelForMaterial
bool m_parseGeoModelForMaterial
Private bool set at initialisation : decides if the GeoModel is parsed or bypassed.
Definition: LArRecoMaterialTool.h:79
LArRecoMaterialTool::ScanHEC
bool ScanHEC(CaloSubdetNames::ALIGNVOL alvol, double &mass, double &x0, double &dEdx, double &aveA, double &aveZ)
Definition: LArRecoMaterialTool.cxx:712
h
CaloSubdetNames::FCAL1_NEG
@ FCAL1_NEG
Definition: CaloSubdetNames.h:43
python.CaloScaleNoiseConfig.type
type
Definition: CaloScaleNoiseConfig.py:78
LArRecoMaterialTool::m_child_dEdXFractions
std::vector< double > m_child_dEdXFractions
Definition: LArRecoMaterialTool.h:166
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:616
declareProperty
#define declareProperty(n, p, h)
Definition: BaseFakeBkgTool.cxx:15
LArRecoMaterialTool::m_alignvol_x0
std::vector< double > m_alignvol_x0
Definition: LArRecoMaterialTool.h:150
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:790
LArRecoMaterialTool::m_x0Fractions
std::vector< double > m_x0Fractions
Definition: LArRecoMaterialTool.h:158
LArRecoMaterialTool::addMaterialFraction
void addMaterialFraction(const GeoLogVol &geoVol, std::vector< double > &volume, std::vector< double > &mass, std::vector< double > &x0, std::vector< double > &dEdX, std::vector< double > &aveA, std::vector< double > &aveZ) const
Private method to add material.
Definition: LArRecoMaterialTool.cxx:1069
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
LArRecoMaterialTool::averageFraction
void averageFraction(std::vector< double > &volumeFractions, std::vector< double > &massFractions, std::vector< double > &x0Fractions, std::vector< double > &dEdXFractions, std::vector< double > &aveAFractions, std::vector< double > &aveZFractions, double &mass, double &volume, double &x0, double &dEdx, double &aveA, double &aveZ, double &rho)
Private method to average results stored in private arrays.
Definition: LArRecoMaterialTool.cxx:1109
CaloSubdetNames::LARCRYO_EC_NEG
@ LARCRYO_EC_NEG
Definition: CaloSubdetNames.h:36
fitman.k
k
Definition: fitman.py:528
LArRecoMaterialTool::m_volumeFractions
std::vector< double > m_volumeFractions
Definition: LArRecoMaterialTool.h:157
LArRecoMaterialTool::m_alignvol_number
std::vector< CaloSubdetNames::ALIGNVOL > m_alignvol_number
Definition: LArRecoMaterialTool.h:148
mapkey::key
key
Definition: TElectronEfficiencyCorrectionTool.cxx:37
Generated on Fri Jun 28 2024 21:52:22 for ATLAS Offline Software by   doxygen 1.8.18