NuclearInteractionFTFSimulator.h

Go to the documentation of this file.

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

#include <vector>

class ParticlePropagator;
class RandomEngineAndDistribution;
class G4ParticleDefinition;
class G4Track;
class G4Step;
class G4TheoFSGenerator;
class G4FTFModel;
class G4ExcitedStringDecay;
class G4LundStringFragmentation;
class G4GeneratorPrecompoundInterface;

class NuclearInteractionFTFSimulator : public MaterialEffectsSimulator
{
public:

  NuclearInteractionFTFSimulator(unsigned int distAlgo, double distCut);

  ~NuclearInteractionFTFSimulator();

private:

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

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

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

  std::vector<const G4ParticleDefinition*> theG4Hadron;
  std::vector<double> theNuclIntLength;
  std::vector<int> theId;

  G4TheoFSGenerator* theHadronicModel;
  G4FTFModel* theStringModel;
  G4ExcitedStringDecay* theStringDecay; 
  G4LundStringFragmentation* theLund;
  G4GeneratorPrecompoundInterface* theCascade; 

  G4Step* dummyStep;
  G4Track* currTrack;
  const G4ParticleDefinition* currParticle;

  G4Nucleus targetNucleus;
  G4HadProjectile theProjectile;
  G4LorentzVector curr4Mom;
  G4ThreeVector vectProj;
  G4ThreeVector theBoost;

  double theEnergyLimit;

  double theDistCut;
  double distMin;

  int numHadrons;
  int currIdx;
  unsigned int theDistAlgo;
};
#endif