ExtEnergyLoss.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#pragma once
 
#include <globals.hh>
#include <G4VContinuousProcess.hh>
 
class G4ParticleDefinition;
class G4Track;
class G4Step;
 
namespace Belle2 {
 
  namespace Simulation {
 
    class EnergyLossForExtrapolator;

    class ExtEnergyLoss: public G4VContinuousProcess {
 
    public:

      ExtEnergyLoss(const G4String& processName = "ExtEnergyLoss", G4ProcessType aType = fElectromagnetic);

      virtual ~ExtEnergyLoss();

      inline G4bool IsApplicable(const G4ParticleDefinition&);

      G4double GetContinuousStepLimit(const G4Track&, G4double, G4double, G4double&);

      G4VParticleChange* AlongStepDoIt(const G4Track&, const G4Step&);

      inline G4double GetStepLimit();

      inline void SetStepLimit(G4double);
 
    private:

      EnergyLossForExtrapolator* m_energyLossForExtrapolator;

      G4double m_StepLimit;
 
    };
 
    inline G4bool ExtEnergyLoss::IsApplicable(const G4ParticleDefinition& aParticle)
    {
      // This process applies to all geant4e-specific charged particles.
      // (Verified that it is only called for these - but let's check anyway.)
      return ((aParticle.GetPDGCharge() != 0) && (aParticle.GetParticleName().compare(0, 4, "g4e_") == 0));
    }
 
    inline G4double ExtEnergyLoss::GetStepLimit()
    {
      return m_StepLimit;
    }
 
    inline void ExtEnergyLoss::SetStepLimit(G4double val)
    {
      m_StepLimit = val;
    }
 
  } //end namespace Simulation

} //end Belle2 namespace