MuNuclearSplittingProcessXSecBias.h

Go to the documentation of this file.

#ifndef LARG4_MNXS_H
#define LARG4_MNXS_H

#include "Geant4/G4ForceCondition.hh"
#include "Geant4/G4GPILSelection.hh"
#include "Geant4/G4Types.hh"
#include "Geant4/G4VParticleChange.hh"
#include "Geant4/G4VProcess.hh"
#include "Geant4/G4WrapperProcess.hh"
#include "Geant4/Randomize.hh"
#include <cmath>
class G4Step;
class G4Track;

namespace larg4 {

  class MuNuclearSplittingProcessXSecBias : public G4WrapperProcess {
    // Override PostStepDoIt method
  public:
    MuNuclearSplittingProcessXSecBias(){};
    ~MuNuclearSplittingProcessXSecBias(){};

    void SetNSplit(G4int nTrx, G4int xB = 0, G4double xFac = 1)
    {
      fNSplit = nTrx, eFactor = (G4double)xFac, xBiasMode = xB;
    };
    void SetIsActive(G4bool doIt) { fActive = doIt; };

    G4VParticleChange* PostStepDoIt(const G4Track& track, const G4Step& step);
    G4VParticleChange* AlongStepDoIt(const G4Track& track, const G4Step& step);
    // GPIL    //////////////
    virtual G4double AlongStepGetPhysicalInteractionLength(const G4Track& track,
                                                           G4double previousStepSize,
                                                           G4double currentMinimumStep,
                                                           G4double& proposedSafety,
                                                           G4GPILSelection* selection);
    virtual G4double PostStepGetPhysicalInteractionLength(const G4Track& track,
                                                          G4double previousStepSize,
                                                          G4ForceCondition* condition);

  protected:
    virtual void ResetNumberOfInteractionLengthLeft()
    {
      G4VProcess::theNumberOfInteractionLengthLeft = -std::log(G4UniformRand());
      theInitialNumberOfInteractionLength = G4VProcess::theNumberOfInteractionLengthLeft;
    }

  private:
    // Data members
    G4int fNSplit;
    G4bool fActive;
    G4int xBiasMode;
    G4double eFactor; // enhancement factor to the cross-setion

    G4VParticleChange fParticleChange;
    // weight change applied at AlongStepDoIt()
    G4double wc;
    G4double theInitialNumberOfInteractionLength;

    G4double XBiasSurvivalProbability();
    G4double XBiasSecondaryWeight();
    G4double GetTotalNumberOfInteractionLengthTraversed();
  };

} // end namespace

#endif // MNSP