FairRoot
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
SetCuts.C
Go to the documentation of this file.
1 /********************************************************************************
2  * Copyright (C) 2014 GSI Helmholtzzentrum fuer Schwerionenforschung GmbH *
3  * *
4  * This software is distributed under the terms of the *
5  * GNU Lesser General Public Licence (LGPL) version 3, *
6  * copied verbatim in the file "LICENSE" *
7  ********************************************************************************/
8 
14 void SetCuts()
15 {
16  // cout << "SetCuts Macro: Setting Processes.." <<endl;
17 
18  // ------>>>> IMPORTANT!!!!
19  // For a correct comparison between GEANE and MC (pull distributions)
20  // or for a simulation without the generation of secondary particles:
21  // 1. set LOSS = 2, DRAY = 0, BREM = 1
22  // 2. set the following cut values: CUTGAM, CUTELE, CUTNEU, CUTHAD, CUTMUO = 1 MeV or less
23  // BCUTE, BCUTM, DCUTE, DCUTM, PPCUTM = 10 TeV
24  // (For an explanation of the chosen values, please refer to the GEANT User's Guide
25  // or to message #5362 in the PandaRoot Forum >> Monte Carlo Engines >> g3Config.C thread)
26  //
27  // The default settings refer to a complete simulation which generates and follows also the secondary particles.
28  TVirtualMC *MC = TVirtualMC::GetMC();
29 
30  MC->SetProcess("PAIR", 1);
31  MC->SetProcess("COMP", 1);
32  MC->SetProcess("PHOT", 1);
33  MC->SetProcess("PFIS", 0);
34  MC->SetProcess("DRAY", 1);
35  MC->SetProcess("ANNI", 1);
36  MC->SetProcess("BREM", 1);
37  MC->SetProcess("HADR", 1);
38  MC->SetProcess("MUNU", 1);
39  MC->SetProcess("DCAY", 1);
40  MC->SetProcess("LOSS", 1);
41  MC->SetProcess("MULS", 1);
43  Double_t cut1 = 1.0E-3; // GeV --> 1 MeV
44  Double_t cutb = 1.0E4; // GeV --> 10 TeV
45  Double_t tofmax = 1.E10; // seconds
46  // cout << "SetCuts Macro: Setting cuts.." <<endl;
47 
48  MC->SetCut("CUTGAM", cut1);
49  MC->SetCut("CUTELE", cut1);
50  MC->SetCut("CUTNEU", cut1);
51  MC->SetCut("CUTHAD", cut1);
52  MC->SetCut("CUTMUO", cut1);
53  MC->SetCut("BCUTE", cut1);
54  MC->SetCut("BCUTM", cut1);
55  MC->SetCut("DCUTE", cut1);
56  MC->SetCut("DCUTM", cut1);
57  MC->SetCut("PPCUTM", cut1);
58  MC->SetCut("TOFMAX", tofmax);
59 }
SetCuts
void SetCuts()
Definition: SetCuts.C:14