PndFtfDirect.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                 PndFtfDirect header file                  -----
// -----               Created 11/04/08  by M.Al-Turany              -----
// -------------------------------------------------------------------------

#ifndef PND_FTFDIRECT_H
#define PND_FTFDIRECT_H

#include "PndTargetGenerator.h"
#include "TF1.h"

namespace CLHEP {
class RanluxEngine;
class Hep3Vector;
} // namespace CLHEP
#include "G4Types.hh"
#include "G4String.hh"

class G4Evaporation;
class G4ExcitationHandler;
class G4PreCompoundModel;
class G4ParticleTable;
class Test30Material;
class Test30Physics;

class G4RunManager;
class G4VUserPhysicsList;
class G4DynamicParticle;
class G4VProcess;
class G4Step;
class G4Track;
class G4ParticleDefinition;
class G4Material;
class TVirtualMCStack;
class FairPrimaryGenerator;

class PndFtfDirect : public PndTargetGenerator {

 public:
  PndFtfDirect();

  PndFtfDirect(const char *configfile);
  PndFtfDirect(const PndFtfDirect &);
  PndFtfDirect(const char *particle, const char *material, int targetA, const char *generator, double mom, int seed = 0, bool noelastic = false);
  PndFtfDirect(const char *particle, const char *material, int targetA, const char *generator, double mom, bool noelastic = false);

  virtual ~PndFtfDirect();

  virtual Bool_t ReadEvent(FairPrimaryGenerator *primGen);
  Bool_t ProcessEvent(FairPrimaryGenerator *primGen);
  void SetVerbose(int v = 1) { fverbose = v; };
  void SetNoElastic(bool n = true) { fNoElastics = n; };
  void LoadG4();

 private:
  void InitZero();        // initialize numbers etc. to defaults
  void parseConfigfile(); // parse the configuration file
  void Setup();           // setup
  bool fDoSetup;
  float fMom;
  float fMode;
  double fSeed;
  double fRsigma;
  float fThtMin;

  G4RunManager *fG4RunManager;             
  G4VUserPhysicsList *fG4VUserPhysicsList; 
  std::ifstream *fin;                      
  Test30Material *fmate;                   
  Test30Physics *fphys;                    
  G4ParticleTable *fpartTable;             

  CLHEP::RanluxEngine *fdefaultEngine; 

  G4int fverbose;
  G4bool fsaverand;
  G4bool fNoElastics;
  //
  G4String fnamePart; 
  G4bool fionParticle;
  G4double fPlab; // VU
  G4double fenergy;
  G4double fsigmae;
  //
  G4int fionZ, fionA;
  G4int ftargetA;
  G4String fnameMat; 
  //
  G4bool fShen;
  G4String fnameGen; 
  //
  G4int fnevt;
  G4int fmodu;
  //
  G4double ftheStep;
  G4Material *fmaterial; 
  //
  G4bool fxsbgg;
  //
  //   // Track
  CLHEP::Hep3Vector *faPosition; 
  G4double faTime;
  CLHEP::Hep3Vector *faDirection; 
  G4double fnx, fny, fnz;

  G4ExcitationHandler *ftheDeExcitation; 
  G4Evaporation *ftheEvaporation;        
  G4PreCompoundModel *fthePreCompound;   

  G4double fsigTot;
  G4double fsigEl;
  G4double fsigIn;

  int fnpart;
  //   Root initialization
  // TFile f1("FTF.root","RECREATE","ROOT_Tree");

  Int_t factiveCnt;
  // TTree* fTree = new TTree("data","FTF Background");
  // TClonesArray* fEvt;

  // fEvt=new TClonesArray("TParticle",100);
  // fTree->Branch("Npart",&activeCnt,"Npart/I");
  // fTree->Branch("Particles",&fEvt, 32000,99);

  G4DynamicParticle *dParticle; 
  G4VProcess *proc;             
  G4Step *step;                 
  G4Track *gTrack;              
  G4ParticleDefinition *part;   

  // G4ParticleDefinition* gamma;
  G4ParticleDefinition *electron; 
  G4ParticleDefinition *proton;   
  G4ParticleDefinition *neutron;  
  G4ParticleDefinition *pin;      
  G4ParticleDefinition *pip;      
  // G4ParticleDefinition* pi0;
  G4ParticleDefinition *deu; 
  G4ParticleDefinition *tri; 
  G4ParticleDefinition *he3; 
  G4ParticleDefinition *alp; 
  // G4ParticleDefinition* ion;

  G4ParticleDefinition *anti_proton;   
  G4ParticleDefinition *anti_neutron;  
  G4ParticleDefinition *anti_deuteron; 
  G4ParticleDefinition *anti_triton;   
  G4ParticleDefinition *anti_He3;      
  G4ParticleDefinition *anti_alpha;    

  void TestPointers()
  {
    std::cout << "PndFtfDirect::fin               -> " << fin << std::endl;
    std::cout << "PndFtfDirect::fmate             -> " << fmate << std::endl;
    std::cout << "PndFtfDirect::fphys             -> " << fphys << std::endl;
    std::cout << "PndFtfDirect::fpartTable        -> " << fpartTable << std::endl;
    std::cout << "PndFtfDirect::fdefaultEngine    -> " << fdefaultEngine << std::endl;
    std::cout << "PndFtfDirect::fmaterial         -> " << fmaterial << std::endl;
    std::cout << "PndFtfDirect::faPosition        -> " << faPosition << std::endl;
    std::cout << "PndFtfDirect::faDirection       -> " << faDirection << std::endl;
    std::cout << "PndFtfDirect::ftheDeExcitation  -> " << ftheDeExcitation << std::endl;
    std::cout << "PndFtfDirect::ftheEvaporation   -> " << ftheEvaporation << std::endl;
    std::cout << "PndFtfDirect::fthePreCompound   -> " << fthePreCompound << std::endl;
    std::cout << "PndFtfDirect::dParticle         -> " << dParticle << std::endl;
    std::cout << "PndFtfDirect::proc              -> " << proc << std::endl;
    std::cout << "PndFtfDirect::step              -> " << step << std::endl;
    std::cout << "PndFtfDirect::gTrack            -> " << gTrack << std::endl;
    std::cout << "PndFtfDirect::part              -> " << part << std::endl;
    std::cout << "PndFtfDirect::electron          -> " << electron << std::endl;
    std::cout << "PndFtfDirect::proton            -> " << proton << std::endl;
    std::cout << "PndFtfDirect::neutron           -> " << neutron << std::endl;
    std::cout << "PndFtfDirect::pin               -> " << pin << std::endl;
    std::cout << "PndFtfDirect::pip               -> " << pip << std::endl;
    std::cout << "PndFtfDirect::deu               -> " << deu << std::endl;
    std::cout << "PndFtfDirect::tri               -> " << tri << std::endl;
    std::cout << "PndFtfDirect::he3               -> " << he3 << std::endl;
    std::cout << "PndFtfDirect::alp               -> " << alp << std::endl;
    std::cout << "PndFtfDirect::anti_proton       -> " << anti_proton << std::endl;
    std::cout << "PndFtfDirect::anti_neutron      -> " << anti_neutron << std::endl;
    std::cout << "PndFtfDirect::anti_deuteron     -> " << anti_deuteron << std::endl;
    std::cout << "PndFtfDirect::anti_triton       -> " << anti_triton << std::endl;
    std::cout << "PndFtfDirect::anti_He3          -> " << anti_He3 << std::endl;
    std::cout << "PndFtfDirect::anti_alpha        -> " << anti_alpha << std::endl;
  };

  ClassDef(PndFtfDirect, 1);
};

#endif