In vmctest/*/g4Config.C:
[u/mrichter/AliRoot.git] / test / vmctest / ppbench / g4Config.C
CommitLineData
e6dc94be 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
9void Config()
10{
11 cout << "Running g4Config.C ... " << endl;
12
13 // AliRoot setup
14 //
15 gROOT->LoadMacro("$ALICE_ROOT/test/vmctest/ppbench/commonConfig.C");
c05b5bd1 16 commonConfig();
e6dc94be 17
d111ffce 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
e6dc94be 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",
c05b5bd1 48// "specialCuts+stackPopper+stepLimiter",
e6dc94be 49// true);
e6dc94be 50
e2750a44 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
e6dc94be 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
e6dc94be 66 geant4->ProcessGeantCommand("/mcVerbose/all 1");
d111ffce 67 geant4->ProcessGeantCommand("/mcVerbose/geometryManager 1");
e6dc94be 68 geant4->ProcessGeantCommand("/mcVerbose/opGeometryManager 1");
d111ffce 69 geant4->ProcessGeantCommand("/mcTracking/loopVerbose 1");
e6dc94be 70 geant4->ProcessGeantCommand("/mcPhysics/rangeCuts 0.01 mm");
e6dc94be 71 geant4->ProcessGeantCommand("/mcPhysics/selectOpProcess Scintillation");
72 geant4->ProcessGeantCommand("/mcPhysics/setOpProcessActivation false");
73 geant4->ProcessGeantCommand("/mcTracking/skipNeutrino true");
d111ffce 74 geant4->ProcessGeantCommand("/mcDet/setMaxStepInLowDensityMaterials 1 cm");
e6dc94be 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
118G4ProcessTable::Insert : arguments are 0 pointer
119G4ProcessTable::Insert : arguments are 0 pointer
120G4ProcessTable::Insert : arguments are 0 pointer
121G4ProcessTable::Insert : arguments are 0 pointer
122G4ProcessTable::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