MuNuclearSplittingProcessXSecBias.h

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
/// \file  MuNuclearSplittingProcessXSecBias.h
/// \brief Check of Geant4 to run the LArSoft detector simulation
///
/// \author  brebel@fnal.gov
////////////////////////////////////////////////////////////////////////

#ifndef LARG4_MNXS_H
#define LARG4_MNXS_H

#include <cmath>
#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"
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