Rec::ParticleCaloCellAssociationTool Class Reference

#include <ParticleCaloCellAssociationTool.h>

Inheritance diagram for Rec::ParticleCaloCellAssociationTool:

Collaboration diagram for Rec::ParticleCaloCellAssociationTool:

Public Types
typedef std::unordered_map< size_t, std::unique_ptr< ParticleCellAssociation > >	Cache
	Method to get the ParticleCellAssociation for a given Particle. More...

Public Member Functions
	ParticleCaloCellAssociationTool (const std::string &, const std::string &, const IInterface *)
virtual	~ParticleCaloCellAssociationTool () override
virtual StatusCode	initialize () override final
virtual StatusCode	finalize () override final
virtual std::unique_ptr< ParticleCellAssociation >	particleCellAssociation (const xAOD::IParticle &particle, float dr, const CaloCellContainer *container=nullptr, const CaloExtensionCollection *extensionCache=nullptr) const override final
	Method to get the ParticleCellAssociation for a given Particle. More...
virtual ParticleCellAssociation *	particleCellAssociation (const xAOD::IParticle &particle, float dr, IParticleCaloCellAssociationTool::Cache &cache, const CaloCellContainer *container=nullptr, const CaloExtensionCollection *extensionCache=nullptr) const override final
	Method to get the ParticleCellAssociation for a given Particle. 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	getCellIntersections (const Trk::CaloExtension &caloExtension, const std::vector< const CaloCell * > &cells, ParticleCellAssociation::CellIntersections &cellIntersections) const
void	associateCells (const CaloCellContainer &container, const Trk::CaloExtension &caloExtension, float dr, std::vector< const CaloCell * > &cells) 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...

Private Attributes
ToolHandle< Trk::IParticleCaloExtensionTool >	m_caloExtensionTool
SG::ReadHandleKey< CaloCellContainer >	m_cellContainerName
Gaudi::Property< double >	m_coneSize { this, "ConeSize", 0.2 }
SG::ReadCondHandleKey< CaloDetDescrManager >	m_caloMgrKey {this,"CaloDetDescrManager", "CaloDetDescrManager"}
Trk::CaloCellSelectorLayerdR	m_defaultSelector
PathLengthUtils	m_pathLenUtil
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 32 of file ParticleCaloCellAssociationTool.h.

Member Typedef Documentation

Cache

typedef std::unordered_map<size_t,std::unique_ptr< ParticleCellAssociation > > Rec::IParticleCaloCellAssociationTool::Cache

inherited

Method to get the ParticleCellAssociation for a given Particle.

Parameters

particle	input particle
dr	cone size used for the association @ param cache cache for keeping previous results
container	cell container to be used if provided
extensionCache	A calo extension cache to be used if provided

Returns

ParticleCellAssociation* (plain ptr cache has ownership)

An alg looping over a single collection of IParticles
re-using them multiple times can use a local cache of the form std::unordered_map<size_t,std::unique_ptr<ParticleCellAssociation >>. where the key is the value of IParticle::index()

This method adds the ParticleCellAssociation to the cache look-up table which retains ownership.

Definition at line 56 of file IParticleCaloCellAssociationTool.h.

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

ParticleCaloCellAssociationTool()

Rec::ParticleCaloCellAssociationTool::ParticleCaloCellAssociationTool	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p
	)

Definition at line 22 of file ParticleCaloCellAssociationTool.cxx.

  : AthAlgTool(t, n, p)
  , m_defaultSelector(0.4)
{
  declareInterface<IParticleCaloCellAssociationTool>(this);
}

~ParticleCaloCellAssociationTool()

Rec::ParticleCaloCellAssociationTool::~ParticleCaloCellAssociationTool ( )

overridevirtualdefault

Member Function Documentation

associateCells()

void Rec::ParticleCaloCellAssociationTool::associateCells ( const CaloCellContainer &  container, const Trk::CaloExtension &  caloExtension, float  dr, std::vector< const CaloCell * > &  cells  ) const

private

Definition at line 424 of file ParticleCaloCellAssociationTool.cxx.

{
  const Trk::TrackParameters* pars = caloExtension.caloEntryLayerIntersection();
  if (!pars) {
    ATH_MSG_DEBUG("associateCells() - NO TrackParameters found in "
                  "caloExtension.caloEntryLayerIntersection()");
    return;
  }
 
  double eta = pars->position().eta();
  double phi = pars->position().phi();
 
  // Use Calorimeter list for CPU reasons
  SG::ReadCondHandle<CaloDetDescrManager> caloMgrHandle{m_caloMgrKey};
  const CaloDetDescrManager* caloMgr=*caloMgrHandle;
  CaloCellList myList(caloMgr,&container);
  myList.select(eta, phi, dr);
  cells.reserve(myList.ncells());
  cells.insert(cells.end(), myList.begin(), myList.end());
  ATH_MSG_DEBUG("associated cells " << cells.size() << " using cone " << dr);
}

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.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

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.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

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.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
  }

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.

  {
    return PBASE::declareProperty(t);
  }

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

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation 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.

  {
    this->declare(hndl.vhKey());
    hndl.vhKey().setOwner(this);
 
    return PBASE::declareProperty(name,hndl,doc);
  }

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

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation 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.

  {
    this->declare(hndl);
    hndl.setOwner(this);
 
    return PBASE::declareProperty(name,hndl,doc);
  }

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.

  {
 
    // std::ostringstream ost;
    // ost << Algorithm::name() << " VHKA declareProp: " << name 
    //     << " size: " << hndArr.keys().size() 
    //     << " mode: " << hndArr.mode() 
    //     << "  vhka size: " << m_vhka.size()
    //     << "\n";
    // debug() << ost.str() << endmsg;
 
    hndArr.setOwner(this);
    m_vhka.push_back(&hndArr);
 
    Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
    if (p != 0) {
      p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
    } else {
      ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                    << name);
    }
 
    return p;
 
  }

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

name	Name of the property.
property	Object holding the property value.
doc	Documentation 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.

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

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

name	Name of the property.
property	Object holding the property value.
doc	Documentation 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.

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

declareProperty() [6/6]

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

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

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

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.

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

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

{ 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 Rec::ParticleCaloCellAssociationTool::finalize ( )

finaloverridevirtual

Definition at line 51 of file ParticleCaloCellAssociationTool.cxx.

{
  return StatusCode::SUCCESS;
}

getCellIntersections()

void Rec::ParticleCaloCellAssociationTool::getCellIntersections ( const Trk::CaloExtension &  caloExtension, const std::vector< const CaloCell * > &  cells, ParticleCellAssociation::CellIntersections &  cellIntersections  ) const

private

Definition at line 143 of file ParticleCaloCellAssociationTool.cxx.

{
  // use 3D pathlength in cells
 
  bool use3D = true;
 
  cellIntersections.reserve(extension.caloLayerIntersections().size() * 1.3);
 
  CaloExtensionHelpers::EntryExitLayerMap entryExitLayerMap;
  CaloExtensionHelpers::entryExitLayerMap(extension, entryExitLayerMap);
  ATH_MSG_DEBUG("EntryExitLayerMap " << entryExitLayerMap.size());
 
  CaloExtensionHelpers::ScalarLayerMap eLossLayerMap, pathLenLayerMap;
  CaloExtensionHelpers::eLossLayerMap(extension, eLossLayerMap);
  CaloExtensionHelpers::pathLenLayerMap(extension, pathLenLayerMap);
 
  ATH_MSG_DEBUG("Getting cells intersections using cells " << cells.size());
  for (const auto* cell : cells) {
    // get sampling and look up entry/exit points
    CaloSampling::CaloSample sample = cell->caloDDE()->getSampling();
 
    auto pos = entryExitLayerMap.find(sample);
    if (pos == entryExitLayerMap.end())
      continue;
    //
    // pos2 and weight2 are introduced because the PreSamplerB has sometimes a
    // very small size
    //         PresamplerB and EMB1 are merged
    //
    auto pos2 = pos;
    if (sample == CaloSampling::PreSamplerB) {
      pos2 = entryExitLayerMap.find(CaloSampling::EMB1);
      if (pos2 == entryExitLayerMap.end()) {
        pos2 = pos;
      }
    }
    double path = 0.;
 
    double drFix = cell->caloDDE()->dr();
    double dzFix = cell->caloDDE()->dz();
    //      double dphi = cell->caloDDE()->dphi();
 
    int isample = cell->caloDDE()->getSampling();
    bool barrel = false;
    if (cell->caloDDE()->getSubCalo() == CaloCell_ID::TILE)
      barrel = true;
    if (sample == CaloSampling::PreSamplerB || sample == CaloSampling::EMB1 ||
        sample == CaloSampling::EMB2 || sample == CaloSampling::EMB3)
      barrel = true;
 
    double drTG = fabs((pos->second.first - pos2->second.second).perp());
    double dzTG = fabs((pos->second.first - pos2->second.second).z());
 
    if (barrel)
      ATH_MSG_VERBOSE(" barrel cell sampling "
                      << cell->caloDDE()->getSampling() << " dr "
                      << cell->caloDDE()->dr() << " drTG " << drTG);
    if (!barrel)
      ATH_MSG_VERBOSE(" endcap cell sampling "
                      << cell->caloDDE()->getSampling() << " dz "
                      << cell->caloDDE()->dz() << " dzTG " << dzTG);
 
    if (drFix == 0.) {
      // recalculate the r values from the other cells
      // BUG/FEATURE: extract dr from cell container for sampling 4 5 6 7 needed
      // EME BUG/FEATURE: extract dr from cell container for sampling 8 9 10 11
      // needed HEC
      if (cell->caloDDE()->deta() > 0) {
        double theta = atan2(cell->caloDDE()->r(), cell->z());
        double dtheta =
          2 * cell->caloDDE()->deta() * sin(theta / 2.) * cos(theta / 2);
        if (theta + dtheta < M_PI) {
          double dr =
            fabs(cell->z() * tan(theta + dtheta) - cell->z() * tan(theta));
          drFix = fabs(dr);
          double detaCheck =
            -log(tan((theta + dtheta) / 2.)) + log(tan((theta) / 2.));
          ATH_MSG_VERBOSE(" FIX cell sampling "
                          << cell->caloDDE()->getSampling() << " deta "
                          << cell->caloDDE()->deta() << " detaCheck "
                          << detaCheck << " drFix " << drFix);
        } else {
          ATH_MSG_WARNING(" FIXR cell sampling failed: theta "
                          << theta << " dtheta " << dtheta << " sum/pi "
                          << (theta + dtheta) * M_1_PI << " deta "
                          << cell->caloDDE()->deta());
        }
        //          ATH_MSG_VERBOSE(" FIX cell sampling deta " << deta << "
        //          dtheta " << dtheta  << "  scale " <<  scale << " theta " <<
        //          theta );
      } else {
        double drMin = 100000.;
        int dscut = 1;
        if (!barrel)
          dscut = 0;
        const CaloCell* cellFound = nullptr;
        for (const auto* celln : cells) {
          if (cell == celln)
            continue;
          if (cell->caloDDE()->getSubCalo() == celln->caloDDE()->getSubCalo()) {
            int dsample = isample - celln->caloDDE()->getSampling();
            if (abs(dsample) == dscut) {
              double drNew = fabs(cell->caloDDE()->r() - celln->caloDDE()->r());
              if (drNew < 1)
                continue;
              if (drNew < drMin) {
                drMin = drNew;
                cellFound = celln;
              }
            }
          }
        }
        drFix = drMin;
        ATH_MSG_VERBOSE(" Problem cell sampling "
                        << cell->caloDDE()->getSampling() << " x "
                        << cell->caloDDE()->x() << " y " << cell->caloDDE()->y()
                        << " z " << cell->caloDDE()->z() << " dr "
                        << cell->caloDDE()->dr() << " drFix " << drFix
                        << " drTG " << drTG);
        if (cellFound)
          ATH_MSG_VERBOSE(" cellFound sampling "
                          << cellFound->caloDDE()->getSampling() << " x "
                          << cellFound->caloDDE()->x() << " y "
                          << cellFound->caloDDE()->y() << " z "
                          << cellFound->caloDDE()->z() << " dr "
                          << cellFound->caloDDE()->dr() << " dscut " << dscut
                          << " drFix " << drFix);
      }
    }
 
    if (dzFix == 0.) {
      // recalculate z values from the other cells
      // BUG/FEATURE: extract dz from cell container for sampling 0 1 2 3 needed
      // EMB
      if (cell->caloDDE()->deta() > 0) {
        double theta = atan2(cell->caloDDE()->r(), cell->z());
        double dtheta =
          2 * cell->caloDDE()->deta() * sin(theta / 2.) * cos(theta / 2);
        if (theta + dtheta < M_PI) {
          double dz = fabs(cell->caloDDE()->r() / tan(theta + dtheta) -
                           cell->caloDDE()->r() / tan(theta));
          dzFix = dz;
        } else {
          ATH_MSG_WARNING(" FIXZ cell sampling failed: theta "
                          << theta << " dtheta " << dtheta << " sum/pi "
                          << (theta + dtheta) * M_1_PI << " deta "
                          << cell->caloDDE()->deta());
        }
        double detaCheck =
          -log(tan((theta + dtheta) / 2.)) + log(tan((theta) / 2.));
        ATH_MSG_VERBOSE(" Fix cell sampling "
                        << cell->caloDDE()->getSampling() << " deta "
                        << cell->caloDDE()->deta() << " detaCheck  "
                        << detaCheck << " dtheta " << dtheta << " dzFix "
                        << dzFix);
      } else {
        double dzMin = 100000.;
        int dscut = 1;
        if (barrel)
          dscut = 0;
        const CaloCell* cellFound = nullptr;
        for (const auto* celln : cells) {
          if (cell == celln)
            continue;
          if (cell->caloDDE()->getSubCalo() == celln->caloDDE()->getSubCalo()) {
            int isample2 = celln->caloDDE()->getSampling();
            if (abs(isample - isample2) == dscut) {
              double dzNew = fabs(cell->caloDDE()->z() - celln->caloDDE()->z());
              if (dzNew < 1)
                continue;
              if (dzNew < dzMin) {
                dzMin = dzNew;
                cellFound = celln;
              }
            }
          }
        }
        dzFix = dzMin;
        ATH_MSG_VERBOSE(" Problem cell sampling "
                        << cell->caloDDE()->getSampling() << " x "
                        << cell->caloDDE()->x() << " y " << cell->caloDDE()->y()
                        << " z " << cell->caloDDE()->z() << " dz "
                        << cell->caloDDE()->dz() << " dzFix " << dzFix
                        << " dzTG " << dzTG);
        if (cellFound)
          ATH_MSG_VERBOSE(" cellFound sampling "
                          << cellFound->caloDDE()->getSampling() << " x "
                          << cellFound->caloDDE()->x() << " y "
                          << cellFound->caloDDE()->y() << " z "
                          << cellFound->caloDDE()->z() << " dz "
                          << cellFound->caloDDE()->dz() << " dscut " << dscut
                          << " dzFix " << dzFix);
      }
    }
    //
    // always use fixed values that correspond to the Calorimeter Tracking
    // Geometry these are different from the CaloCell values
    //
 
    if (cell->energy() > 50.)
      ATH_MSG_DEBUG(" cell sampling and size "
                    << cell->caloDDE()->getSampling() << " cell energy "
                    << cell->energy() << " dzFix " << dzFix << " dzTG " << dzTG
                    << " drFix " << drFix << " drTG " << drTG << " barrel "
                    << barrel);
 
    if (!barrel)
      dzFix = dzTG;
    if (barrel)
      drFix = drTG;
 
    if (use3D) {
      // m_pathLenUtil.pathInsideCell( *cell, entryExitLayerMap);
      double pathInMM = PathLengthUtils::get3DPathLength(
        *cell, pos->second.first, pos2->second.second, drFix, dzFix);
      double totpath = (pos->second.first - pos2->second.second).mag();
      path = totpath != 0 ? pathInMM / totpath : 0.;
      if (path > 0 || cell->energy() > 50.) {
        ATH_MSG_DEBUG(" cell sampling and size "
                      << cell->caloDDE()->getSampling() << " cell energy "
                      << cell->energy() << " drFix " << drFix << " dzFix "
                      << dzFix << " path " << path << " length TG " << totpath);
        ATH_MSG_DEBUG(" cell dr " << cell->caloDDE()->dr() << " cell dz "
                                  << cell->caloDDE()->dz() << " deta "
                                  << cell->caloDDE()->deta());
      }
    }
 
    double path2 = 0.;
 
    if (!use3D)
      path2 = pathInsideCell(*cell, pos->second.first, pos2->second.second);
 
    if (path2 <= 0. && path <= 0.)
      continue;
 
    // auto entrancePair = entryExitLayerMap.find(entranceID);
    auto eLossPair = eLossLayerMap.find(sample);
    double eLoss = 0.;
    //
    // Just store total expected eloss
    //
    if (eLossPair != eLossLayerMap.end()) {
      eLoss = eLossPair->second;
      if (sample == CaloSampling::PreSamplerB) {
        auto eLossPair2 = eLossLayerMap.find(CaloSampling::EMB1);
        if (eLossPair2 != eLossLayerMap.end()) {
          eLoss = 0.5 * (eLossPair->second) + 0.5 * (eLossPair2->second);
        }
      } else if (sample == CaloSampling::EMB1) {
        auto eLossPair2 = eLossLayerMap.find(CaloSampling::PreSamplerB);
        if (eLossPair2 != eLossLayerMap.end()) {
          eLoss = 0.5 * (eLossPair->second) + 0.5 * (eLossPair2->second);
        }
      }
    } // IF
 
    ATH_MSG_DEBUG(" PATH3D = "
                  << path << " PATH2D = " << path2 << " eLoss " << eLoss
                  << " cell energy " << (cell)->energy() << " radius "
                  << cell->caloDDE()->r() << " phi " << cell->caloDDE()->phi()
                  << " dr " << cell->caloDDE()->dr() << " dphi "
                  << cell->caloDDE()->dphi() << " x " << cell->caloDDE()->x()
                  << " y " << cell->caloDDE()->y() << " z "
                  << cell->caloDDE()->z() << " dx " << cell->caloDDE()->dx()
                  << " dy " << cell->caloDDE()->dy() << " dz "
                  << cell->caloDDE()->dz() << " volume "
                  << cell->caloDDE()->volume());
 
    cellIntersections.emplace_back(
      cell, new ParticleCellIntersection(*cell, eLoss, use3D ? path : path2));
  }
  ATH_MSG_DEBUG(" added cell intersections  " << cellIntersections.size());
}

initialize()

StatusCode Rec::ParticleCaloCellAssociationTool::initialize ( )

finaloverridevirtual

Definition at line 35 of file ParticleCaloCellAssociationTool.cxx.

{
  /* Retrieve track extrapolator from ToolService */
  ATH_CHECK(m_caloExtensionTool.retrieve());
 
  m_defaultSelector.setConeSize(m_coneSize);
 
  if (!m_cellContainerName.key().empty()) {
    ATH_CHECK(m_cellContainerName.initialize());
  }
  
  ATH_CHECK(m_caloMgrKey.initialize());
  return StatusCode::SUCCESS;
}

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 & Rec::IParticleCaloCellAssociationTool::interfaceID ( )

inlinestaticinherited

Definition at line 65 of file IParticleCaloCellAssociationTool.h.

                                                                        {
  return IID_IParticleCaloCellAssociationTool; 
}

msg() [1/2]

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msg() [2/2]

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

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                  {
    return this->msgStream(lvl);
  }

msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

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.

particleCellAssociation() [1/2]

std::unique_ptr< ParticleCellAssociation > Rec::ParticleCaloCellAssociationTool::particleCellAssociation ( const xAOD::IParticle &  particle, float  dr, const CaloCellContainer *  container = nullptr, const CaloExtensionCollection *  extensionCache = nullptr  ) const

finaloverridevirtual

Method to get the ParticleCellAssociation for a given Particle.

Parameters

particle	input particle
dr	cone size used for the association
container	cell container to be used if provided

Returns

std::unique_ptr<ParticleCellAssociation>

Implements Rec::IParticleCaloCellAssociationTool.

Definition at line 57 of file ParticleCaloCellAssociationTool.cxx.

{
  ATH_MSG_DEBUG(" particleCellAssociation: ptr " << &particle << " dr " << dr);
  // get the extrapolation into the calo
  std::unique_ptr<const Trk::CaloExtension> caloExtensionUPtr;
  const Trk::CaloExtension* caloExtension = nullptr;
  if (extensionCache)
    caloExtension =
      m_caloExtensionTool->caloExtension(particle, *extensionCache);
  else {
    caloExtensionUPtr = m_caloExtensionTool->caloExtension(
      Gaudi::Hive::currentContext(), particle);
    caloExtension = caloExtensionUPtr.get();
  }
  if (!caloExtension) {
    ATH_MSG_DEBUG("Failed to get calo extension");
    return nullptr;
  }
  if (caloExtension->caloLayerIntersections().empty()) {
    ATH_MSG_DEBUG(
      "Received a caloExtension object without track extrapolation");
    return nullptr;
  }
  // retrieve the cell container if not provided, return false it retrieval
  // failed
  if (!container) {
    if (m_cellContainerName.key().empty()) {
      ATH_MSG_DEBUG("Failed to get calo cell container");
      return nullptr;
    }
    SG::ReadHandle<CaloCellContainer> cccHdl(m_cellContainerName);
    container = &(*cccHdl);
  }
  std::vector<const CaloCell*> cells;
  // update cone size in case it is smaller than the default
  if (dr < m_coneSize) {
    dr = m_coneSize;
  }
  associateCells(*container, *caloExtension, dr, cells);
 
  // get cell intersections
  ParticleCellAssociation::CellIntersections cellIntersections;
  getCellIntersections(*caloExtension, cells, cellIntersections);
  if (!caloExtensionUPtr)
    // Have to manually copy the calo extension object. Clearly the class
    // wants to be shared through a shared_ptr but this clearly is not an
    // option
    caloExtensionUPtr = std::make_unique<Trk::CaloExtension>(
      caloExtension->caloEntryLayerIntersection()
        ? caloExtension->caloEntryLayerIntersection()->clone()
        : nullptr,
      caloExtension->muonEntryLayerIntersection()
        ? caloExtension->muonEntryLayerIntersection()->clone()
        : nullptr,
      std::vector<Trk::CurvilinearParameters>(
        caloExtension->caloLayerIntersections()));
  return std::make_unique<ParticleCellAssociation>(caloExtensionUPtr.release(),
                                                   std::move(cells),
                                                   dr,
                                                   std::move(cellIntersections),
                                                   container);
}

particleCellAssociation() [2/2]

ParticleCellAssociation * Rec::ParticleCaloCellAssociationTool::particleCellAssociation ( const xAOD::IParticle &  particle, float  dr, IParticleCaloCellAssociationTool::Cache &  cache, const CaloCellContainer *  container = nullptr, const CaloExtensionCollection *  extensionCache = nullptr  ) const

finaloverridevirtual

Method to get the ParticleCellAssociation for a given Particle.

Parameters

particle	input particle
dr	cone size used for the association @ param cache cache for keeping previous results
container	cell container to be used if provided

Returns

ParticleCellAssociation* (plain ptr cache has ownership)

An alg looping over a single collection of IParticles re-using them multiple times can use a local cache of the form std::unordered_map<size_t,std::unique_ptr<ParticleCellAssociation >>. where the key is the value of IParticle::index()

This method adds the ParticleCellAssociation to the cache look-up table which retains ownership.

Implements Rec::IParticleCaloCellAssociationTool.

Definition at line 125 of file ParticleCaloCellAssociationTool.cxx.

{
  /*if not there , default ctor for unique_ptr (nullptr)*/
  std::unique_ptr<ParticleCellAssociation>& association =
    cache[particle.index()];
  if (association == nullptr) {
    association =
      particleCellAssociation(particle, dr, container, extensionCache);
  }
  return association.get();
}

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.

  {
    h.renounce();
    PBASE::renounce (h);
  }

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.

                                                          {
    handlesArray.renounce();
  }

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.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

m_caloExtensionTool

ToolHandle<Trk::IParticleCaloExtensionTool> Rec::ParticleCaloCellAssociationTool::m_caloExtensionTool

private

Initial value:

{
    this,
    "ParticleCaloExtensionTool",
    ""
  }

Definition at line 94 of file ParticleCaloCellAssociationTool.h.

m_caloMgrKey

SG::ReadCondHandleKey<CaloDetDescrManager> Rec::ParticleCaloCellAssociationTool::m_caloMgrKey {this,"CaloDetDescrManager", "CaloDetDescrManager"}

private

Definition at line 106 of file ParticleCaloCellAssociationTool.h.

m_cellContainerName

SG::ReadHandleKey<CaloCellContainer> Rec::ParticleCaloCellAssociationTool::m_cellContainerName

private

Initial value:

{ this,
                                                            "CaloCellContainer",
                                                            "AllCalo" }

Definition at line 100 of file ParticleCaloCellAssociationTool.h.

m_coneSize

Gaudi::Property<double> Rec::ParticleCaloCellAssociationTool::m_coneSize { this, "ConeSize", 0.2 }

private

Definition at line 104 of file ParticleCaloCellAssociationTool.h.

m_defaultSelector

Trk::CaloCellSelectorLayerdR Rec::ParticleCaloCellAssociationTool::m_defaultSelector

private

Definition at line 108 of file ParticleCaloCellAssociationTool.h.

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_pathLenUtil

PathLengthUtils Rec::ParticleCaloCellAssociationTool::m_pathLenUtil

private

Definition at line 110 of file ParticleCaloCellAssociationTool.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:

• ParticleCaloCellAssociationTool.h
• ParticleCaloCellAssociationTool.cxx