FullModelHadronicProcess.h

Go to the documentation of this file.

#ifndef SimG4Core_CustomPhysics_FullModelHadronicProcess_h
#define SimG4Core_CustomPhysics_FullModelHadronicProcess_h 1

#include "globals.hh"
#include "G4VDiscreteProcess.hh"
#include "G4EnergyRangeManager.hh"
#include "G4Nucleus.hh"
#include "G4ReactionProduct.hh"
#include <vector>
#include "G4HadronicException.hh"

#include "SimG4Core/CustomPhysics/interface/FullModelReactionDynamics.h"

class G4ProcessHelper;

class FullModelHadronicProcess : public G4VDiscreteProcess
{
public:
  
  FullModelHadronicProcess(G4ProcessHelper * aHelper, 
               const G4String& processName = "FullModelHadronicProcess");
    
  ~FullModelHadronicProcess() override;

  
  G4bool IsApplicable(const G4ParticleDefinition& aP) override;
  
  G4VParticleChange *PostStepDoIt(const G4Track &aTrack, const G4Step &aStep) override;

protected:
  
  const G4ParticleDefinition *theParticle;
  G4ParticleDefinition *newParticle;
  
  G4ParticleChange theParticleChange;
  
private:    
  
  virtual G4double GetMicroscopicCrossSection( const G4DynamicParticle *aParticle, 
                           const G4Element *anElement, 
                           G4double aTemp );

  G4double GetMeanFreePath(const G4Track& aTrack, G4double, G4ForceCondition*) override;
  
  void CalculateMomenta( G4FastVector<G4ReactionProduct,MYGHADLISTSIZE> &vec,
             G4int &vecLen,
             const G4HadProjectile *originalIncident,
             const G4DynamicParticle *originalTarget,
             G4ReactionProduct &modifiedOriginal,
             G4Nucleus &targetNucleus,
             G4ReactionProduct &currentParticle,
             G4ReactionProduct &targetParticle,
             G4bool &incidentHasChanged,
             G4bool &targetHasChanged,
             G4bool quasiElastic );

  G4bool MarkLeadingStrangeParticle(
                    const G4ReactionProduct &currentParticle,
                    const G4ReactionProduct &targetParticle,
                    G4ReactionProduct &leadParticle );
    
  void Rotate(G4FastVector<G4ReactionProduct,MYGHADLISTSIZE> &vec, G4int &vecLen);

  const G4DynamicParticle* FindRhadron(G4ParticleChange*);

  G4ProcessHelper* theHelper;
  G4bool toyModel;
  G4double cache;
  G4ThreeVector what;
};

#endif