3 //=======================================================================
4 // Create the output file
6 TFile *rootfile = new TFile("galice.root","recreate");
7 rootfile->SetCompressionLevel(2);
9 //=======================================================================
10 // ******* GEANT STEERING parameters FOR ALICE SIMULATION *******
11 geant3->SetTRIG(1); //Number of events to be processed
13 geant3->SetDEBU(0,0,1);
14 //geant3->SetSWIT(2,2);
25 geant3->SetHADR(1); //Select pure GEANH (HADR 1) or GEANH/NUCRIN (HADR 3)
29 geant3->SetAUTO(1); //Select automatic STMIN etc... calc. (AUTO 1) or manual (AUTO 0)
30 geant3->SetABAN(0); //Restore 3.16 behaviour for abandoned tracks
31 geant3->SetOPTI(2); //Select optimisation level for GEANT geometry searches (0,1,2)
32 Float_t cut = 1.e-3; // 1MeV cut by default
33 Float_t tofmax = 1.e10;
34 // GAM ELEC NHAD CHAD MUON EBREM MUHAB EDEL MUDEL MUPA TOFMAX
35 geant3->SetCUTS(cut,cut, cut, cut, cut, cut, cut, cut, cut, cut, tofmax);
37 //=======================================================================
38 // ************* STEERING parameters FOR ALICE SIMULATION **************
39 // --- Specify event type to be tracked through the ALICE setup
40 // --- All positions are in cm, angles in degrees, and P and E in GeV
41 AliGenHIJINGpara *gener = new AliGenHIJINGpara(1000);
42 gener->SetMomentumRange(0,100);
43 gener->SetPhiRange(0,360);
44 gener->SetThetaRange(20,160);
45 gener->SetOrigin(0,0,0); //vertex position
46 gener->SetSigma(0,0,0); //Sigma in (X,Y,Z) (cm) on IP position
49 gAlice->SetField(-999,2); //Specify maximum magnetic field in Tesla (neg. ==> default field)
51 //=================== Alice BODY parameters =============================
52 AliBODY *BODY = new AliBODY("BODY","Alice envelop");
56 //============================ TPC parameters ================================
57 // --- This allows the user to specify sectors for the SLOW (TPC geometry 2)
58 // --- Simulator. TPCSECL (TPSECU) <0 means that ALL lower (upper)
59 // --- sectors are specified, any value other than that requires at least one
60 // --- sector (lower or upper)to be specified!
61 // --- Reminder: sectors 1-24 are lower sectors (1-12 -> z>0, 13-24 -> z<0)
62 // --- sectors 25-72 are the upper ones (25-48 -> z>0, 49-72 -> z<0)
63 // --- TPCLOWS - number of lower sectors specified (up to 6)
64 // --- TPCUPS - number of upper sectors specified (up to 12)
65 // --- TPCSENS - sensitive strips for the Slow Simulator !!!
66 // --- This does NOT work if all S or L-sectors are specified, i.e.
67 // --- if TPCSECL or TPCSECU < 0
68 //-----------------------------------------------------------------------------
70 AliTPC *TPC = new AliTPCv1("TPC","Normal TPC");
73 TPC->SetSecLows(1, -1, -1, -1, -1, -1);
74 TPC->SetSecUps(25, 26, 48, -1, -1, -1, -1, -1, -1, -1, -1, -1);