EnergyLossUpdator.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/

// EnergyLossUpdator.h, (c) ATLAS Detector software
 
#ifndef TRKEXTOOLS_ENERGYLOSSUPDATOR_H
#define TRKEXTOOLS_ENERGYLOSSUPDATOR_H
 
// Gaudi
#include "AthenaBaseComps/AthAlgTool.h"
// Trk
#include "TrkEventPrimitives/PropDirection.h"
#include "TrkExInterfaces/IEnergyLossUpdator.h"
#include "TrkExUtils/MaterialInteraction.h"
#include "TrkEventPrimitives/ParticleHypothesis.h"
 
// STL
#include <utility>
 
namespace Trk {
 
class MaterialProperties;

 
class EnergyLossUpdator final
  : public AthAlgTool
  , virtual public IEnergyLossUpdator
{
 
public:
  EnergyLossUpdator(const std::string&, const std::string&, const IInterface*);

  virtual ~EnergyLossUpdator() = default;

  virtual double dEdX(const MaterialProperties& mat,
                      double p,
                      ParticleHypothesis particle = pion) const override final;

  virtual EnergyLoss energyLoss(
    const MaterialProperties& mat,
    double p,
    double pathcorrection,
    PropDirection dir = alongMomentum,
    ParticleHypothesis particle = pion,
    bool useMPV = false) const override final;

  virtual EnergyLoss updateEnergyLoss(
    EnergyLoss& eLoss,
    double caloEnergy,
    double caloEnergyError,
    double pCaloEntry,
    double momentumError,
    int& elossFlag) const override final;

  virtual void getX0ElossScales(int icalo,
                                double eta,
                                double phi,
                                double& X0Scale,
                                double& ElossScale) const override final;
 
private:
  double varianceDeltaQoverP(
    const MaterialProperties&,
    double p,
    double pathcorrection,
    PropDirection direction = alongMomentum,
    ParticleHypothesis particleHypothesis = electron) const;
 
  Trk::EnergyLoss ionizationEnergyLoss(
    const MaterialProperties& mat,
    double p,
    double pathcorrection,
    PropDirection dir = alongMomentum,
    ParticleHypothesis particle = pion) const;
 
  BooleanProperty m_detailedEloss{this, "DetailedEloss", true,
    "provide extended EnergyLoss info"};
  BooleanProperty m_optimalRadiation{this, "OptimalRadiation", true,
    "use calorimeter more optimal for radiation detection"};
};
 
} // end of namespace
 
#endif // TRKEXTOOLS_ENERGYLOSSUPDATOR_H