db/d68/NuclearInteractionFTFSimulator_8h_source.html

#ifndef FastSimulation_MaterialEffects_NuclearInteractionFTFSimulator_H

#define FastSimulation_MaterialEffects_NuclearInteractionFTFSimulator_H


#include "FastSimulation/MaterialEffects/interface/MaterialEffectsSimulator.h"


#include "G4Nucleus.hh"

#include "G4HadProjectile.hh"

#include "G4LorentzVector.hh"

#include "G4ThreeVector.hh"


class ParticlePropagator;

class RandomEngineAndDistribution;

class G4ParticleDefinition;

class G4Track;

class G4Step;

class G4TheoFSGenerator;

class G4FTFModel;

class G4ExcitedStringDecay;

class G4LundStringFragmentation;

class G4GeneratorPrecompoundInterface;

class G4CascadeInterface;

class G4DiffuseElastic;

class G4PhysicsLogVector;


static const int numHadrons = 30;

static const int npoints = 21;


class NuclearInteractionFTFSimulator : public MaterialEffectsSimulator {

public:

  NuclearInteractionFTFSimulator(unsigned int distAlgo, double distCut, double elimit, double eth);


  ~NuclearInteractionFTFSimulator() override;


private:

  void compute(ParticlePropagator& Particle, RandomEngineAndDistribution const*) override;


  void saveDaughter(ParticlePropagator& Particle, const G4LorentzVector& lv, int pdgid);


  double distanceToPrimary(const RawParticle& Particle, const RawParticle& aDaughter) const;


  static const G4ParticleDefinition* theG4Hadron[numHadrons];

  static int theId[numHadrons];

  G4PhysicsLogVector* vect;


  G4TheoFSGenerator* theHadronicModel;

  G4FTFModel* theStringModel;

  G4ExcitedStringDecay* theStringDecay;

  G4LundStringFragmentation* theLund;

  G4GeneratorPrecompoundInterface* theCascade;


  G4CascadeInterface* theBertiniCascade;

  G4DiffuseElastic* theDiffuseElastic;


  G4Step* dummyStep;

  G4Track* currTrack;

  const G4ParticleDefinition* currParticle;


  G4Nucleus targetNucleus;

  G4HadProjectile theProjectile;

  G4LorentzVector curr4Mom;

  G4ThreeVector vectProj;

  G4ThreeVector theBoost;


  double theBertiniLimit;

  double theEnergyLimit;


  double theDistCut;

  double distMin;


  double intLengthElastic;

  double intLengthInelastic;


  int currIdx;

  size_t index;

  unsigned int theDistAlgo;

};

#endif