TPC/TPCConfig.C

   1 void Config()
   2 {
   3 //=======================================================================
   4 //  Create the output file
   5
   6 TFile *rootfile = new TFile("galice.root","recreate");
   7 rootfile->SetCompressionLevel(2);
   8
   9 //=======================================================================
  10 // ******* GEANT STEERING parameters FOR ALICE SIMULATION *******
  11 geant3->SetTRIG(1); //Number of events to be processed
  12 geant3->SetSWIT(4,1);
  13 geant3->SetDEBU(0,0,1);
  14 //geant3->SetSWIT(2,2);
  15 geant3->SetDCAY(1);
  16 geant3->SetPAIR(1);
  17 geant3->SetCOMP(1);
  18 geant3->SetPHOT(1);
  19 geant3->SetPFIS(0);
  20 geant3->SetDRAY(0);
  21 geant3->SetANNI(1);
  22 geant3->SetBREM(1);
  23 geant3->SetMUNU(1);
  24 geant3->SetCKOV(1);
  25 geant3->SetHADR(1); //Select pure GEANH (HADR 1) or GEANH/NUCRIN (HADR 3)
  26 geant3->SetLOSS(2);
  27 geant3->SetMULS(1);
  28 geant3->SetRAYL(1);
  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);
  36 //
  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
  47 gener->Init();
  48
  49 gAlice->SetField(-999,2);    //Specify maximum magnetic field in Tesla (neg. ==> default field)
  50
  51 //=================== Alice BODY parameters =============================
  52 AliBODY *BODY = new AliBODY("BODY","Alice envelop");
  53
  54
  55
  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 //-----------------------------------------------------------------------------
  69
  70 AliTPC *TPC  = new AliTPCv1("TPC","Normal TPC");
  71 TPC->SetSecAL(1);
  72 TPC->SetSecAU(1);
  73 TPC->SetSecLows(1, -1, -1, -1, -1, -1);
  74 TPC->SetSecUps(25, 26, 48, -1, -1, -1, -1, -1, -1, -1, -1, -1);
  75 TPC->SetSens(1);
  76
  77
  78
  79
  80 }