In vmctest/*/g4Config.C:
[u/mrichter/AliRoot.git] / test / vmctest / ppbench / g4Config.C
1 // $Id$
2 //
3 // Configuration macro for running aliroot with Geant4
4 // with primary events read from external file.
5 //
6 // By I. Hrivnacova, IPN Orsay
7         
8
9 void Config()
10 {
11   cout << "Running g4Config.C ... " << endl;
12
13   // AliRoot setup
14   //
15   gROOT->LoadMacro("$ALICE_ROOT/test/vmctest/ppbench/commonConfig.C");
16   commonConfig();
17
18   // TPC primary ionization 
19   AliTPC* TPC = (AliTPC*)gAlice->GetDetector("TPC");
20   if ( ! TPC )
21     cerr << "Cannot get TPC detector" << endl;
22   else  { 
23     cerr << "Setting TPC primary ionization" << endl;
24     TPC->SetPrimaryIonisation(); // not used with Geant3
25   }  
26
27   // Load Geant4 + Geant4 VMC libraries
28   //
29   if (gClassTable->GetID("TGeant4") == -1) {
30     // Load Geant4 libraries 
31     if (!gInterpreter->IsLoaded("$ALICE/geant4_vmc/examples/macro/g4libs.C")) {
32       gROOT->LoadMacro("$ALICE/geant4_vmc/examples/macro/g4libs.C");
33       gInterpreter->ProcessLine("g4libs()");
34     }
35   }    
36
37   // Create Geant4 VMC
38   //  
39   TGeant4 *geant4 = 0;
40   if ( ! gMC ) {
41     TG4RunConfiguration* runConfiguration 
42       = new TG4RunConfiguration("geomRoot", 
43                                 "QGSP_BERT_EMV+optical", 
44                                 "specialCuts+stackPopper+stepLimiter",
45                                  true);
46 //      = new TG4RunConfiguration("geomRootToGeant4",
47 //                                "QGSP_BERT_EMV+optical", 
48 //                                "specialCuts+stackPopper+stepLimiter",
49 //                                 true);
50
51     geant4 = new TGeant4("TGeant4", 
52                          "The Geant4 Monte Carlo : QGSP_BERT_EMV+optical", 
53                          runConfiguration);
54              // Repeat physics selection in the title; to be removed
55              // with new geant4_vmc tag (1.13)            
56
57     cout << "Geant4 has been created." << endl;
58   } 
59   else {
60     cout << "Monte Carlo has been already created." << endl;
61   }  
62
63   // Customization of Geant4 VMC
64   //
65
66     geant4->ProcessGeantCommand("/mcVerbose/all 1");  
67     geant4->ProcessGeantCommand("/mcVerbose/geometryManager 1");  
68     geant4->ProcessGeantCommand("/mcVerbose/opGeometryManager 1");  
69     geant4->ProcessGeantCommand("/mcTracking/loopVerbose 1");     
70     geant4->ProcessGeantCommand("/mcPhysics/rangeCuts 0.01 mm"); 
71     geant4->ProcessGeantCommand("/mcPhysics/selectOpProcess Scintillation");
72     geant4->ProcessGeantCommand("/mcPhysics/setOpProcessActivation false");
73     geant4->ProcessGeantCommand("/mcTracking/skipNeutrino true");
74     geant4->ProcessGeantCommand("/mcDet/setMaxStepInLowDensityMaterials 1 cm");
75
76   // Uncomment this line to get a detail info from each step 
77   // geant4->ProcessGeantCommand("/tracking/verbose 1");  
78   
79   // More info from the physics list
80   // the verbosity level is passed to all contained physics lists and their
81   // physics builders
82   geant4->ProcessGeantCommand("/mcVerbose/composedPhysicsList 2");  
83   
84   // More info from optical processes
85   //geant4->ProcessGeantCommand("/mcVerbose/opticalPhysicsList 3");  
86   
87   // More info from geometry building
88   //geant4->ProcessGeantCommand("/mcVerbose/geometryManager 1");  
89
90   // More info from setting geometry properties (in materials and surfaces)
91   // for optical physics
92   //geant4->ProcessGeantCommand("/mcVerbose/opGeometryManager 1");  
93   
94   // More info about regions construction 
95   // and conversion of VMC cuts in cuts in range per regions 
96   //geant4->ProcessGeantCommand("/mcVerbose/regionsManager 2");
97   //geant4->ProcessGeantCommand("/mcRegions/print true");
98
99   
100   // Suppress verbose info from tracks which reached maximum number of steps
101   // (default value is 30000)  
102   //geant4->ProcessGeantCommand("/mcTracking/loopVerbose 0"); 
103     
104   //
105   // Set apply cuts 
106 /*
107   geant4->ProcessGeantCommand("/run/particle/applyCuts");  
108   // geant4->ProcessGeantCommand("/mcVerbose/geometryManager 2");  
109
110   geant4->ProcessGeantCommand("/mcVerbose/composedPhysicsList 2");  
111   geant4->ProcessGeantCommand("/mcDet/volNameSeparator !");
112   geant4->ProcessGeantCommand("/mcPhysics/setStackPopperSelection e+ e- pi+ pi- kaon+ kaon- gamma");
113   //geant4->ProcessGeantCommand("/tracking/verbose 1");  
114
115   geant4->ProcessGeantCommand("/mcControl/g3Defaults");
116 !!!!n Generates warnings:
117 >>> Event 0
118 G4ProcessTable::Insert : arguments are 0 pointer 
119 G4ProcessTable::Insert : arguments are 0 pointer 
120 G4ProcessTable::Insert : arguments are 0 pointer 
121 G4ProcessTable::Insert : arguments are 0 pointer 
122 G4ProcessTable::Insert : arguments are 0 pointer 
123
124 */
125
126   // AliRoot event generator
127   // (it has to be created after MC, as it may use decayer via VMC)
128   //
129   gROOT->LoadMacro("$ALICE_ROOT/test/vmctest/ppbench/genPPbenchConfig.C");
130   gROOT->LoadMacro("$ALICE_ROOT/test/vmctest/ppbench/genExtFileConfig.C");
131             // The event generator selection (srun) is done in genPPbenchConfig.C
132             // that´s why we have to load it too
133   genExtFileConfig(srun);
134
135   cout << "Running g4Config.C finished ... " << endl;
136 }