1 enum gentype_t {hijing, gun, box, pythia, param, cocktail, fluka, halo, ntuple, scan, doublescan};
3 gentype_t gentype=param;
9 new AliGeant3("C++ Interface to Geant3");
11 //=======================================================================
12 // Create the output file
14 TFile *rootfile = new TFile("galice.root","recreate");
15 rootfile->SetCompressionLevel(2);
16 TGeant3 *geant3 = (TGeant3*)gMC;
19 //=======================================================================
20 // ******* GEANT STEERING parameters FOR ALICE SIMULATION *******
21 geant3->fGctrak->maxnst=1000000;
23 geant3->SetTRIG(1); //Number of events to be processed
24 geant3->SetSWIT(4,100);
25 geant3->SetDEBU(0,0,1);
36 geant3->SetHADR(4); //Select pure GEANH (HADR 1) or GEANH/NUCRIN (HADR 3)
40 geant3->SetAUTO(0); //Select automatic STMIN etc... calc. (AUTO 1) or manual (AUTO 0)
41 geant3->SetABAN(0); //Restore 3.16 behaviour for abandoned tracks
42 geant3->SetOPTI(2); //Select optimisation level for GEANT geometry searches (0,1,2)
43 Float_t cut = 1.e-1; // 100MeV cut by default
44 Float_t tofmax = 1.e10;
45 // GAM ELEC NHAD CHAD MUON EBREM MUHAB EDEL MUDEL MUPA TOFMAX
46 geant3->SetCUTS(1.e-4, 1.e-4, 1.e-3, 1.e-4, 1.e-4, cut, cut, cut, cut, cut, 1.e-5);
48 gAlice->TrackingLimits( 700, 2200);
51 //=======================================================================
52 // ************* STEERING parameters FOR ALICE SIMULATION **************
53 // --- Specify event type to be tracked through the ALICE setup
54 // --- All positions are in cm, angles in degrees, and P and E in GeV
59 //*********************************************
60 // Example for Fixed Particle Gun *
61 //*********************************************
62 AliGenFixed *gener = new AliGenFixed(ntracks);
63 gener->SetMomentum(3);
64 gener->SetPhiRange(0);
65 gener->SetThetaRange(0);
66 gener->SetOrigin(0,0,0); //vertex position
67 gener->SetPart(kProton); //GEANT particle type
70 //*********************************************
71 // Example for Moving Particle Gun *
72 //*********************************************
73 AliGenBox *gener = new AliGenBox(ntracks);
74 gener->SetMomentumRange(3,4);
75 gener->SetPhiRange(-360,360);
76 gener->SetThetaRange(2., 10.);
77 gener->SetOrigin(25,25,510.5);
78 gener->SetVertexSmear(perTrack);
80 gener->SetSigma(1.8, 1.8,0); //Sigma in (X,Y,Z) (cm) on IP position
81 gener->SetPart(kProton); //GEANT particle type
84 //*********************************************
85 // Scanning on a grid *
86 //*********************************************
87 AliGenScan *gener = new AliGenScan(-1);
88 gener->SetMomentumRange(20,20);
89 gener->SetPhiRange(90,90);
90 gener->SetThetaRange(0,0);
92 gener->SetSigma(3,3,0); //Sigma in (X,Y,Z) (cm) on IP position
93 gener->SetPart(kMuonMinus);
94 gener->SetRange(20, -100, 100, 20, -100, 100, 1, 500, 500);
97 //*********************************************
98 // Scanning on a grid *
99 //*********************************************
100 AliGenDoubleScan *gener = new AliGenDoubleScan(-1);
101 gener->SetMomentumRange(4,4);
102 gener->SetPhiRange(0,360);
103 gener->SetThetaRange(0,0);
105 gener->SetSigma(3,3,0); //Sigma in (X,Y,Z) (cm) on IP position
107 gener->SetRange(20, -100, 100, 20, -100, 100, 1, 500, 500);
108 gener->SetDistance(1);
113 AliGenHIJINGpara *gener = new AliGenHIJINGpara(ntracks);
114 gener->SetMomentumRange(0,999);
115 gener->SetPhiRange(0,360);
116 gener->SetThetaRange(0.104,33.52);
117 gener->SetOrigin(0,0,0); //vertex position
118 gener->SetSigma(0,0,5.6); //Sigma in (X,Y,Z) (cm) on IP position
122 //********************************************
123 // Example for Charm Production with Pythia *
124 //********************************************
126 AliGenPythia *gener = new AliGenPythia(ntracks);
127 gener->SetMomentumRange(0,999);
128 gener->SetPhiRange(0,360);
129 gener->SetThetaRange(0., 180.);
130 gener->SetYRange(-10,10);
131 gener->SetPtRange(0,100);
132 gener->SetOrigin(0,0,0); // vertex position
133 gener->SetVertexSmear(perEvent);
134 gener->SetSigma(0,0,5.6); // Sigma in (X,Y,Z) (cm) on IP position
135 // gener->SetStrucFunc(DO_Set_1);
136 gener->SetProcess(mb);
137 gener->SetEnergyCMS(5500.);
141 //*******************************************************
142 // Example for J/psi Production from Parameterisation *
143 //*******************************************************
144 AliGenParam *gener = new AliGenParam(ntracks,upsilon_p);
145 // AliGenParam *gener = new AliGenParam(ntracks, jpsi_p);
146 gener->SetMomentumRange(0,999);
147 gener->SetPhiRange(0,360);
148 gener->SetYRange(2.5,4);
149 gener->SetPtRange(0,50);
150 gener->SetOrigin(0,0,0); //vertex position
151 gener->SetSigma(0,0,0); //Sigma in (X,Y,Z) (cm) on IP position
152 gener->SetForceDecay(dimuon);
153 gener->SetCutOnChild(1);
154 gener->SetChildThetaRange(2,9);
155 gener->SetTrackingFlag(1);
160 //*******************************************************
161 // Example for a FLUKA Boundary Source *
162 //*******************************************************
163 AliGenFLUKAsource *gener = new AliGenFLUKAsource(-1);
164 gener->AddFile("$(ALICE_ROOT)/data/all32.root");
166 gener->SetPartFlag(7);
167 gener->SetMomentumRange(0,999);
168 gener->SetPhiRange(0,360);
169 gener->SetThetaRange(0., 180.);
170 gener->SetAgeMax(1.e-5);
173 gener->SetFraction(0.0315);
177 //*******************************************************
178 // Example for reading from a external file *
179 //*******************************************************
180 AliGenExtFile *gener = new AliGenExtFile(-1);
181 gener->SetFileName("$(ALICE_ROOT)/data/pbpb.root");
182 gener->SetThetaRange(0.104,33.52);
183 gener->SetOrigin(0,0,0); //vertex position
184 gener->SetSigma(0,0,5.6); //Sigma in (X,Y,Z) (cm) on IP position
185 gener->SetVertexSmear(perTrack);
186 gener->SetTrackingFlag(1);
190 //*******************************************************
191 // Example for Tunnel Halo Source *
192 //*******************************************************
193 AliGenHalo *gener = new AliGenHalo(ntracks);
194 gener->SetFileName("/h1/morsch/marsip/marsip5.mu");
198 //*******************************************************
199 // Example for a Cocktail *
200 //*******************************************************
201 AliGenCocktail *gener = new AliGenCocktail();
202 gener->SetMomentumRange(0,10);
203 gener->SetPhiRange(0,360);
204 gener->SetThetaRange(45.,135);
206 pions = new AliGenParam(100, pion_p);
207 // kaons = new AliGenParam(10 , kaon_p);
208 // protons = new AliGenParam(10 , proton_p);
209 gener->AddGenerator(pions , "Pions" , 100);
210 // gener->AddGenerator(kaons , "Kaons" , 10);
211 // gener->AddGenerator(protons, "Protons", 10);
216 // Activate this line if you want the vertex smearing to happen
221 gAlice->SetField(-2,1); //Specify maximum magnetic field in Tesla (neg. ==> default field)
229 //=================== Alice BODY parameters =============================
230 AliBODY *BODY = new AliBODY("BODY","Alice envelop");
233 //=================== MAG parameters ============================
234 // --- Start with Magnet since detector layouts may be depending ---
235 // --- on the selected Magnet dimensions ---
236 AliMAG *MAG = new AliMAG("MAG","Magnet");
240 //=================== ABSO parameters ============================
241 AliABSO *ABSO = new AliABSOv1("ABSO","Muon Absorber");
245 //=================== DIPO parameters ============================
247 AliDIPO *DIPO = new AliDIPOv2("DIPO","Dipole version 2");
252 //=================== SHIL parameters ============================
253 AliSHIL *SHIL = new AliSHILv1("SHIL","Shielding");
258 //=================== PIPE parameters ============================
259 AliPIPE *PIPE = new AliPIPEv0("PIPE","Beam Pipe");
264 //=================== MUON parameters ===========================
266 AliMUON *MUON = new AliMUONv0("MUON","normal MUON");
267 MUON->SetAcceptance(kTRUE, 0.1, 12.);
269 MUON->SetMaxStepGas(0.1);
270 MUON->SetMaxStepAlu(0.1);