ExtEnergyLoss.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2010-2011  Belle II Collaboration                         *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Leo Piilonen                                             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#ifndef EXTENERGYLOSS_H_
#define EXTENERGYLOSS_H_
 
#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
 
#endif /* EXTENERGYLOSS_H_ */