macros/g4ConfigCommon.C

   1 // $Id$
   2 //
   3 // AliRoot Configuration for running aliroot with Monte Carlo.
   4 // Called from g4Config.C
   5 //
   6 // By I. Hrivnacova, IPN Orsay
   7
   8 enum PprGeo_t
   9   {
  10     kHoles, kNoHoles
  11   };
  12 static PprGeo_t geo = kHoles;
  13
  14 Float_t EtaToTheta(Float_t arg);
  15 void    LoadPythia();
  16 static Int_t    eventsPerRun = 50;
  17
  18 void ConfigCommon(Bool_t setRootGeometry = kTRUE)
  19 {
  20   cout << "Running ConfigCommon.C ... " << endl;
  21
  22   // Load Pythia libraries
  23   LoadPythia();
  24
  25   // Create the output file
  26   AliRunLoader* rl = 0;
  27   rl = AliRunLoader::Open("galice.root",
  28                            AliConfig::GetDefaultEventFolderName(),
  29                            "recreate");
  30   if (!rl) {
  31     gAlice->Fatal("Config.C","Can not instatiate the Run Loader");
  32     return;
  33   }
  34   rl->SetCompressionLevel(2);
  35   rl->SetNumberOfEventsPerFile(3);
  36   gAlice->SetRunLoader(rl);
  37
  38   // Set Root geometry file
  39   //
  40   if ( setRootGeometry ) {
  41     gAlice->SetRootGeometry();
  42     gAlice->SetGeometryFromFile("geometry.root");
  43   }
  44
  45   // Set the trigger configuration
  46   AliSimulation::Instance()->SetTriggerConfig("Pb-Pb");
  47   cout<<"Trigger configuration is set to  Pb-Pb"<<endl;
  48
  49     // Set Random Number seed
  50     gRandom->SetSeed(123456); // Set 0 to use the currecnt time
  51     AliLog::Message(AliLog::kInfo, Form("Seed for random number generation = %d",gRandom->GetSeed()), "Config.C", "Config.C", "Config()","Config.C", __LINE__);
  52     int     nParticles = 100;
  53     if (gSystem->Getenv("CONFIG_NPARTICLES"))
  54     {
  55         nParticles = atoi(gSystem->Getenv("CONFIG_NPARTICLES"));
  56     }
  57
  58     AliGenCocktail *gener = new AliGenCocktail();
  59     gener->SetPhiRange(0, 360);
  60     // Set pseudorapidity range from -8 to 8.
  61     Float_t thmin = EtaToTheta(8);   // theta min. <---> eta max
  62     Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min
  63     gener->SetThetaRange(thmin,thmax);
  64     gener->SetOrigin(0, 0, 0);  //vertex position
  65     gener->SetSigma(0, 0, 0);   //Sigma in (X,Y,Z) (cm) on IP position
  66
  67     AliGenHIJINGpara *hijingparam = new AliGenHIJINGpara(nParticles);
  68     hijingparam->SetMomentumRange(0.2, 999);
  69     gener->AddGenerator(hijingparam,"HIJING PARAM",1);
  70     gener->Init();
  71
  72     //
  73     // Activate this line if you want the vertex smearing to happen
  74     // track by track
  75     //
  76     //gener->SetVertexSmear(perTrack);
  77
  78   // =============================
  79   // Magnetic field
  80   // =============================
  81
  82     // Field (L3 0.4 T)
  83     TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
  84
  85     Int_t   iABSO  =  1;
  86     Int_t   iDIPO  =  1;
  87     Int_t   iFMD   =  1;
  88     Int_t   iFRAME =  1;
  89     Int_t   iHALL  =  1;
  90     Int_t   iITS   =  1;
  91     Int_t   iMAG   =  1;
  92     Int_t   iMUON  =  1;
  93     Int_t   iPHOS  =  1;
  94     Int_t   iPIPE  =  1;
  95     Int_t   iPMD   =  1;
  96     Int_t   iHMPID =  1;
  97     Int_t   iSHIL  =  1;
  98     Int_t   iT0    =  1;
  99     Int_t   iTOF   =  1;
 100     Int_t   iTPC   =  1;
 101     Int_t   iTRD   =  1;
 102     Int_t   iZDC   =  1;
 103     Int_t   iEMCAL =  1;
 104     Int_t   iACORDE = 0;
 105     Int_t   iVZERO =  1;
 106 /*
 107     Int_t   iABSO  =  0;
 108     Int_t   iDIPO  =  0;
 109     Int_t   iFMD   =  0;
 110     Int_t   iFRAME =  0;
 111     Int_t   iHALL  =  0;
 112     Int_t   iITS   =  0;
 113     Int_t   iMAG   =  0;
 114     Int_t   iMUON  =  0;
 115     Int_t   iPHOS  =  0;
 116     Int_t   iPIPE  =  0;
 117     Int_t   iPMD   =  0;
 118     Int_t   iHMPID =  0;
 119     Int_t   iSHIL  =  0;
 120     Int_t   iT0    =  0;
 121     Int_t   iTOF   =  0;
 122     Int_t   iTPC   =  1;
 123     Int_t   iTRD   =  0;
 124     Int_t   iZDC   =  0;
 125     Int_t   iEMCAL =  0;
 126     Int_t   iACORDE = 0;
 127     Int_t   iVZERO =  0;
 128 */
 129 // Exluded for problems
 130
 131     rl->CdGAFile();
 132     //=================== Alice BODY parameters =============================
 133     AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
 134
 135     if (iMAG)
 136     {
 137         //=================== MAG parameters ============================
 138         // --- Start with Magnet since detector layouts may be depending ---
 139         // --- on the selected Magnet dimensions ---
 140         AliMAG *MAG = new AliMAG("MAG", "Magnet");
 141     }
 142
 143
 144     if (iABSO)
 145     {
 146         //=================== ABSO parameters ============================
 147         AliABSO *ABSO = new AliABSOv3("ABSO", "Muon Absorber");
 148     }
 149
 150     if (iDIPO)
 151     {
 152         //=================== DIPO parameters ============================
 153
 154         AliDIPO *DIPO = new AliDIPOv3("DIPO", "Dipole version 3");
 155     }
 156
 157     if (iHALL)
 158     {
 159         //=================== HALL parameters ============================
 160
 161         AliHALL *HALL = new AliHALLv3("HALL", "Alice Hall");
 162     }
 163
 164
 165     if (iFRAME)
 166     {
 167         //=================== FRAME parameters ============================
 168
 169         AliFRAMEv2 *FRAME = new AliFRAMEv2("FRAME", "Space Frame");
 170         if (geo == kHoles) {
 171           FRAME->SetHoles(1);
 172         } else {
 173           FRAME->SetHoles(0);
 174         }
 175     }
 176
 177     if (iSHIL)
 178     {
 179         //=================== SHIL parameters ============================
 180
 181         AliSHIL *SHIL = new AliSHILv3("SHIL", "Shielding Version 3");
 182     }
 183
 184
 185     if (iPIPE)
 186     {
 187         //=================== PIPE parameters ============================
 188
 189         AliPIPE *PIPE = new AliPIPEv3("PIPE", "Beam Pipe");
 190     }
 191
 192     if (iITS)
 193     {
 194         //=================== ITS parameters ============================
 195
 196         AliITS *ITS  = new AliITSv11("ITS","ITS v11");
 197     }
 198
 199     if (iTPC)
 200     {
 201         //============================ TPC parameters ===================
 202         AliTPC *TPC = new AliTPCv2("TPC", "Default");
 203     }
 204
 205
 206     if (iTOF) {
 207         //=================== TOF parameters ============================
 208         AliTOF *TOF = new AliTOFv6T0("TOF", "normal TOF");
 209     }
 210
 211
 212     if (iHMPID)
 213     {
 214         //=================== HMPID parameters ===========================
 215         AliHMPID *HMPID = new AliHMPIDv3("HMPID", "normal HMPID");
 216
 217     }
 218
 219
 220     if (iZDC)
 221     {
 222         //=================== ZDC parameters ============================
 223
 224         AliZDC *ZDC = new AliZDCv3("ZDC", "normal ZDC");
 225     }
 226
 227     if (iTRD)
 228     {
 229         //=================== TRD parameters ============================
 230
 231         AliTRD *TRD = new AliTRDv1("TRD", "TRD slow simulator");
 232     }
 233
 234     if (iFMD)
 235     {
 236         //=================== FMD parameters ============================
 237         AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
 238    }
 239
 240     if (iMUON)
 241     {
 242         //=================== MUON parameters ===========================
 243         // New MUONv1 version (geometry defined via builders)
 244         AliMUON *MUON = new AliMUONv1("MUON", "default");
 245     }
 246     //=================== PHOS parameters ===========================
 247
 248     if (iPHOS)
 249     {
 250         AliPHOS *PHOS = new AliPHOSv1("PHOS", "IHEP");
 251     }
 252
 253
 254     if (iPMD)
 255     {
 256         //=================== PMD parameters ============================
 257         AliPMD *PMD = new AliPMDv1("PMD", "normal PMD");
 258     }
 259
 260     if (iT0)
 261     {
 262         //=================== T0 parameters ============================
 263         AliT0 *T0 = new AliT0v1("T0", "T0 Detector");
 264     }
 265
 266     if (iEMCAL)
 267     {
 268         //=================== EMCAL parameters ============================
 269         AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETE");
 270     }
 271
 272      if (iACORDE)
 273     {
 274         //=================== ACORDE parameters ============================
 275         AliACORDE *ACORDE = new AliACORDEv1("ACORDE", "normal ACORDE");
 276     }
 277
 278      if (iVZERO)
 279     {
 280         //=================== ACORDE parameters ============================
 281         AliVZERO *VZERO = new AliVZEROv7("VZERO", "normal VZERO");
 282     }
 283
 284      AliLog::Message(AliLog::kInfo, "End of Config", "Config.C", "Config.C", "Config()"," Config.C", __LINE__);
 285
 286      cout << "Running ConfigCommon.C finished ... " << endl;
 287
 288 }
 289
 290 Float_t EtaToTheta(Float_t arg){
 291   return (180./TMath::Pi())*2.*atan(exp(-arg));
 292 }
 293
 294 void LoadPythia()
 295 {
 296     // Load Pythia related libraries
 297     gSystem->Load("liblhapdf.so");      // Parton density functions
 298     gSystem->Load("libEGPythia6.so");   // TGenerator interface
 299     gSystem->Load("libpythia6.so");     // Pythia
 300     gSystem->Load("libAliPythia6.so");  // ALICE specific implementations
 301 }