SetCuts.C

Go to the documentation of this file.

/** Configuration macro for setting common cuts and processes for G3, G4 and Fluka
    specific cuts and processes to g3 or g4 should be set in the g3Config.C, g4Config.C or flConfig.C
**/
#include "TVirtualMC.h"
 
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); /** pair production*/
    gMC->SetProcess("COMP", 1); /**Compton scattering*/
    gMC->SetProcess("PHOT", 1); /** photo electric effect */
    gMC->SetProcess("PFIS", 0); /**photofission*/
    gMC->SetProcess("DRAY", 1); /**delta-ray*/
    gMC->SetProcess("ANNI", 1); /**annihilation*/
    gMC->SetProcess("BREM", 1); /**bremsstrahlung*/
    gMC->SetProcess("HADR", 3); /**hadronic process*/
    gMC->SetProcess("MUNU", 1); /**muon nuclear interaction*/
    gMC->SetProcess("DCAY", 1); /**decay*/
    gMC->SetProcess("LOSS", 1); /**energy loss*/
    gMC->SetProcess("MULS", 1); /**multiple scattering*/
 
    Double_t cut1 = 1.0E-3;    // GeV --> 1 MeV
    Double_t tofmax = 1.E10;   // seconds
    cout << "SetCuts Macro: Setting cuts.." << endl;
 
    gMC->SetCut("CUTGAM", cut1);   /** gammas (GeV)*/
    gMC->SetCut("CUTELE", cut1);   /** electrons (GeV)*/
    gMC->SetCut("CUTNEU", cut1);   /** neutral hadrons (GeV)*/
    gMC->SetCut("CUTHAD", cut1);   /** charged hadrons (GeV)*/
    gMC->SetCut("CUTMUO", cut1);   /** muons (GeV)*/
    gMC->SetCut("BCUTE", cut1);    /** electron bremsstrahlung (GeV)*/
    gMC->SetCut("BCUTM", cut1);    /** muon and hadron bremsstrahlung(GeV)*/
    gMC->SetCut("DCUTE", cut1);    /** delta-rays by electrons (GeV)*/
    gMC->SetCut("DCUTM", cut1);    /** delta-rays by muons (GeV)*/
    gMC->SetCut("PPCUTM", cut1);   /** direct pair production by muons (GeV)*/
    gMC->SetCut("TOFMAX", tofmax); /**time of flight cut in seconds*/
}