g4Cuts.C File Reference

Go to the source code of this file.

Functions
void	SetCuts ()

Function Documentation

void SetCuts

(

)

Configuration macro for setting common cuts and processes for G3, G4 and Fluka (M. Al-Turany 27.03.2008) specific cuts and processes to g3 or g4 should be set in the g3Config.C, g4Config.C or flConfig.C

pair production

Compton scattering

photo electric effect

photofission

delta-ray

annihilation

bremsstrahlung

hadronic process

muon nuclear interaction

decay

energy loss

multiple scattering

gammas (GeV)

electrons (GeV)

neutral hadrons (GeV)

charged hadrons (GeV)

muons (GeV)

electron bremsstrahlung (GeV)

muon and hadron bremsstrahlung(GeV)

delta-rays by electrons (GeV)

delta-rays by muons (GeV)

direct pair production by muons (GeV)

time of flight cut in seconds

Definition at line 7 of file g4Cuts.C.

References Double_t.

{
  cout << "SetCuts Macro: Setting Processes.." <<endl;
   
  // ------>>>> IMPORTANT!!!!
  // For a correct comparison between GEANE and MC (pull distributions) 
  // or for a simulation without the generation of secondary particles:
  // 1. set LOSS = 2, DRAY = 0, BREM = 1
  // 2. set the following cut values: CUTGAM, CUTELE, CUTNEU, CUTHAD, CUTMUO = 1 MeV or less
  //                                  BCUTE, BCUTM, DCUTE, DCUTM, PPCUTM = 10 TeV
  // (For an explanation of the chosen values, please refer to the GEANT User's Guide
  // or to message #5362 in the PandaRoot Forum >> Monte Carlo Engines >> g3Config.C thread)
  // 
  // The default settings refer to a complete simulation which generates and follows also the secondary particles.


  gMC->SetProcess("PAIR",1); 
  gMC->SetProcess("COMP",1); 
  gMC->SetProcess("PHOT",1); 
  gMC->SetProcess("PFIS",0); 
  gMC->SetProcess("DRAY",1); 
  gMC->SetProcess("ANNI",1); 
  gMC->SetProcess("BREM",1); 
  gMC->SetProcess("HADR",1); 
  gMC->SetProcess("MUNU",1); 
  gMC->SetProcess("DCAY",1);          //ww 1
  gMC->SetProcess("LOSS",1);    //ww 1
  gMC->SetProcess("MULS",1); 

  Double_t cut1 = 1.0E-3;         // GeV --> 1 MeV
  Double_t cutb = 1.0E4;          // GeV --> 10 TeV
  Double_t tofmax = 1.E10;        // seconds
  cout << "SetCuts Macro: Setting cuts.." <<endl;
  
  gMC->SetCut("CUTGAM",cut1);   
  gMC->SetCut("CUTELE",cut1);   
  gMC->SetCut("CUTNEU",cut1);   
  gMC->SetCut("CUTHAD",cut1);   
  gMC->SetCut("CUTMUO",cut1);   
  gMC->SetCut("BCUTE",cut1);    
  gMC->SetCut("BCUTM",cut1);     
  gMC->SetCut("DCUTE",cut1);    
  gMC->SetCut("DCUTM",cut1);    
  gMC->SetCut("PPCUTM",cut1);   
  gMC->SetCut("TOFMAX",tofmax); 

  // //for pdf generation
  // gMC->SetProcess("PAIR",0); /** pair production*/
  // gMC->SetProcess("COMP",0); /**Compton scattering*/
  // gMC->SetProcess("PHOT",0); /** photo electric effect */
  // gMC->SetProcess("PFIS",0); /**photofission*/
  // gMC->SetProcess("DRAY",0); /**delta-ray*/
  // gMC->SetProcess("ANNI",0); /**annihilation*/
  // gMC->SetProcess("BREM",0); /**bremsstrahlung*/ //ww 1
  // gMC->SetProcess("HADR",0); /**hadronic process*/ //ww 1
  // gMC->SetProcess("MUNU",0); /**muon nuclear interaction*/
  // gMC->SetProcess("DCAY",0); /**decay*/         //ww 1
  // gMC->SetProcess("LOSS",0); /**energy loss*/   //ww 1
  // gMC->SetProcess("MULS",0); /**multiple scattering*/  
}