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