g4Cuts.C

Go to the documentation of this file.

void SetCuts()
{
  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*/  
}