G4LongLivedNeutralDecay.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.                  *
 **************************************************************************/
 
// modified from GEANT4 G4Transportation class
 
#pragma once
#include "G4ios.hh"
#include "globals.hh"
#include "G4VRestDiscreteProcess.hh"
#include "G4ParticleChangeForDecay.hh"
#include "G4DecayProcessType.hh"
#include <framework/core/Module.h>
 
namespace Belle2 {
  class G4LongLivedNeutral;
  class G4LongLivedNeutralDecay : public G4VRestDiscreteProcess {
 
  public:
    G4LongLivedNeutralDecay(const G4String& processName = "LongLivedNeutralDecay");
 
    virtual ~G4LongLivedNeutralDecay();
 
  private:
    G4LongLivedNeutralDecay(const G4LongLivedNeutralDecay& right);
 
    G4LongLivedNeutralDecay& operator=(const G4LongLivedNeutralDecay& right);
 
  public:
 
    virtual G4VParticleChange* PostStepDoIt(
      const G4Track& aTrack,
      const G4Step& aStep
    ) override;
 
    virtual G4VParticleChange* AtRestDoIt(
      const G4Track& aTrack,
      const G4Step&  aStep
    ) override;
 
    virtual void BuildPhysicsTable(const G4ParticleDefinition&) override;
    virtual G4bool IsApplicable(const G4ParticleDefinition&) override;
  protected:
 
    virtual G4VParticleChange* DecayIt(
      const G4Track& aTrack,
      const G4Step&  aStep
    );
 
  public:
 
    virtual G4double AtRestGetPhysicalInteractionLength(
      const G4Track& track,
      G4ForceCondition* condition
    ) override;
 
    virtual G4double PostStepGetPhysicalInteractionLength(
      const G4Track& track,
      G4double   previousStepSize,
      G4ForceCondition* condition
    ) override;
 
  protected:
 
 
    virtual G4double GetMeanFreePath(const G4Track& aTrack,
                                     G4double   previousStepSize,
                                     G4ForceCondition* condition
                                    ) override;
    virtual G4double GetMeanLifeTime(const G4Track& aTrack,
                                     G4ForceCondition* condition
                                    ) override;
  public:
 
    virtual void StartTracking(G4Track*) override;
    virtual void EndTracking() override;
  public:
 
    G4double GetRemainderLifeTime() const;
  protected:
    G4int verboseLevel;
  protected:
 
    const G4double HighestValue; 
    G4double fRemainderLifeTime; 
    G4ParticleChangeForDecay fParticleChangeForDecay; 
  };
 
  inline
  G4VParticleChange* G4LongLivedNeutralDecay::AtRestDoIt(
    const G4Track& aTrack,
    const G4Step&   aStep
  )
  {
    return DecayIt(aTrack, aStep);
  }
 
 
  inline
  G4double G4LongLivedNeutralDecay::GetRemainderLifeTime() const
  {
    return fRemainderLifeTime;
  }
}