Bug fix in multiplicity limits.
[u/mrichter/AliRoot.git] / macros / g4Config.C
1 // $Id$
2 //
3 // Configuration macro for running aliroot with Geant4.
4 // All AliRoot specifics are defined in g4ConfigCommon.C
5 // 
6 // In order to run aliroot with Geant4, you have to Initialize
7 // AliRun with this g4Config.C
8 //
9 // aliroot
10 // root [0] gAlice->Init("g4Config.C");        
11 // root [1] gAlice->Run();      
12 //
13 // You can also start from a mini GUI - g4menu.C.
14 // See description in this macro.
15 //
16 // By I. Hrivnacova, IPN Orsay
17         
18
19 void Config()
20 {
21   cout << "Running g4Config.C ... " << endl;
22
23   // AliRoot setup
24   //
25   gROOT->LoadMacro("$ALICE_ROOT/macros/g4ConfigCommon.C");
26   ConfigCommon();
27
28   // Load Geant4 + Geant4 VMC libraries
29   //
30   if (gClassTable->GetID("TGeant4") == -1) {
31     // Load Geant4 libraries 
32     if (!gInterpreter->IsLoaded("$ALICE/geant4_vmc/examples/macro/g4libs.C")) {
33       gROOT->LoadMacro("$ALICE/geant4_vmc/examples/macro/g4libs.C");
34       gInterpreter->ProcessLine("g4libs()");
35     }
36   }    
37
38   // Create Geant4 VMC
39   //  
40   TGeant4 *geant4 = 0;
41   if (!gMC) {
42     TG4RunConfiguration* runConfiguration 
43       = new TG4RunConfiguration("geomRoot", 
44                                 "emStandard+optical", 
45                                 "specialCuts+specialControls+stackPopper",
46                                  true);
47 //      = new TG4RunConfiguration("geomRootToGeant4",
48 //                                "emStandard+optical", 
49 //                                "specialCuts+specialControls+stackPopper",
50 //                                 true);
51       
52
53     geant4 = new TGeant4("TGeant4", "The Geant4 Monte Carlo", runConfiguration);
54     cout << "Geant4 has been created." << endl;
55   } else {
56     cout << "Monte Carlo has been already created." << endl;
57   }  
58     
59   //
60   // Set External decayer
61   TVirtualMCDecayer *decayer = new AliDecayerPythia();
62
63   decayer->SetForceDecay(kAll);
64   decayer->Init();
65   gMC->SetExternalDecayer(decayer);
66   
67   //=======================================================================
68   // ************* STEERING parameters FOR ALICE SIMULATION **************
69   // --- Specify event type to be tracked through the ALICE setup
70   // --- All positions are in cm, angles in degrees, and P and E in GeV
71
72   gMC->SetProcess("DCAY",1);
73   gMC->SetProcess("PAIR",1);
74   gMC->SetProcess("COMP",1);
75   gMC->SetProcess("PHOT",1);
76   gMC->SetProcess("PFIS",0);
77   gMC->SetProcess("DRAY",0);
78   gMC->SetProcess("ANNI",1);
79   gMC->SetProcess("BREM",1);
80   gMC->SetProcess("MUNU",1);
81   gMC->SetProcess("CKOV",1);
82   gMC->SetProcess("HADR",1);
83   gMC->SetProcess("LOSS",2);
84   gMC->SetProcess("MULS",1);
85   //gMC->SetProcess("RAYL",1);
86
87   Float_t cut = 1.e-3;        // 1MeV cut by default
88   Float_t tofmax = 1.e10;
89
90   gMC->SetCut("CUTGAM", cut);
91   gMC->SetCut("CUTELE", cut);
92   gMC->SetCut("CUTNEU", cut);
93   gMC->SetCut("CUTHAD", cut);
94   gMC->SetCut("CUTMUO", cut);
95   gMC->SetCut("BCUTE",  cut); 
96   gMC->SetCut("BCUTM",  cut); 
97   gMC->SetCut("DCUTE",  cut); 
98   gMC->SetCut("DCUTM",  cut); 
99   gMC->SetCut("PPCUTM", cut);
100   gMC->SetCut("TOFMAX", tofmax); 
101
102   // Set apply cuts 
103   geant4->ProcessGeantCommand("/run/particle/applyCuts");  
104   // geant4->ProcessGeantCommand("/mcVerbose/geometryManager 2");  
105   geant4->ProcessGeantCommand("/mcVerbose/composedPhysicsList 2");  
106   geant4->ProcessGeantCommand("/mcDet/volNameSeparator !");
107   geant4->ProcessGeantCommand("/mcControl/g3Defaults");
108   geant4->ProcessGeantCommand("/mcPhysics/setStackPopperSelection e+ e- pi+ pi- kaon+ kaon- gamma");
109   // geant4->ProcessGeantCommand("/tracking/verbose 1");  
110
111   cout << "Running g4Config.C finished ... " << endl;
112 }