Macro checking if the current $ALICE_ROOT creates a geometry different w.r.t. to...
[u/mrichter/AliRoot.git] / macros / g4ConfigCommon.C
index e718067..69b56a4 100644 (file)
@@ -2,15 +2,25 @@
 //
 // AliRoot Configuration for running aliroot with Monte Carlo.
 // Called from g4Config.C
+//
+// By I. Hrivnacova, IPN Orsay
+
+enum PprGeo_t
+  {
+    kHoles, kNoHoles
+  };
+static PprGeo_t geo = kHoles;
 
 Float_t EtaToTheta(Float_t arg);
-static Int_t    eventsPerRun = 100;
+void    LoadPythia();
+static Int_t    eventsPerRun = 50;
 
-void ConfigCommon(Bool_t interactiveSetup)
+void ConfigCommon(Bool_t setRootGeometry = kTRUE)
 {
-  // ============================= 
-  // Root file
-  // ============================= 
+  cout << "Running ConfigCommon.C ... " << endl;
+
+  // Load Pythia libraries
+  LoadPythia();
 
   // Create the output file
   AliRunLoader* rl = 0;
@@ -24,134 +34,104 @@ void ConfigCommon(Bool_t interactiveSetup)
   rl->SetCompressionLevel(2);
   rl->SetNumberOfEventsPerFile(3);
   gAlice->SetRunLoader(rl);
+  
+  // Set Root geometry file
+  //
+  if ( setRootGeometry ) {
+    gAlice->SetRootGeometry();
+    gAlice->SetGeometryFromFile("geometry.root");
+  }
 
-  // Set External decayer
-  AliDecayer* decayer = new AliDecayerPythia();
-  decayer->SetForceDecay(kAll);
-  decayer->Init();
-  gMC->SetExternalDecayer(decayer);
-
-  // Physics process control
-  gMC ->SetProcess("DCAY",1);
-  gMC ->SetProcess("PAIR",1);
-  gMC ->SetProcess("COMP",1);
-  gMC ->SetProcess("PHOT",1);
-  gMC ->SetProcess("PFIS",0);
-  gMC ->SetProcess("DRAY",0);
-  gMC ->SetProcess("ANNI",1);
-  gMC ->SetProcess("BREM",1);
-  gMC ->SetProcess("MUNU",1);
-  //xx gMC ->SetProcess("CKOV",1);
-  gMC ->SetProcess("HADR",1); //Select pure GEANH (HADR 1) or GEANH/NUCRIN (HADR 3)
-  gMC ->SetProcess("LOSS",2);
-  gMC ->SetProcess("MULS",1);
-  //xx gMC ->SetProcess("RAYL",1);
-
-  // Energy cuts
-  // (in development)
-  Float_t cut    = 1.e-3; // 1MeV cut by default
-  Float_t tofmax = 1.e10; 
-
-  gMC ->SetCut("CUTGAM",cut);
-  gMC ->SetCut("CUTELE",cut);
-  gMC ->SetCut("CUTNEU",cut);
-  gMC ->SetCut("CUTHAD",cut);
-  gMC ->SetCut("CUTMUO",cut);
-  gMC ->SetCut("BCUTE",cut);
-  gMC ->SetCut("BCUTM",cut);
-  gMC ->SetCut("DCUTE",cut);
-  gMC ->SetCut("DCUTM",cut);
-  gMC ->SetCut("PPCUTM",cut);
-  gMC ->SetCut("TOFMAX",tofmax);
-
-  // ============================= 
-  // Event generator
-  // ============================= 
-
-  // --- Specify event type to be tracked through the ALICE setup
-  // --- All positions are in cm, angles in degrees, and P and E in GeV
-
-  Int_t nParticles = 100;
-  if (gSystem->Getenv("CONFIG_NPARTICLES")) 
-    nParticles = atoi(gSystem->Getenv("CONFIG_NPARTICLES"));
-
-  AliGenCocktail *gener = new AliGenCocktail();
-  gener->SetPhiRange(0, 360);
-  // Set pseudorapidity range from -8 to 8.
-  Float_t thmin = EtaToTheta(8);   // theta min. <---> eta max
-  Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min 
-  gener->SetThetaRange(thmin,thmax);
-  gener->SetOrigin(0, 0, 0);  //vertex position
-  gener->SetSigma(0, 0, 0);   //Sigma in (X,Y,Z) (cm) on IP position
-
-  AliGenHIJINGpara *hijingparam = new AliGenHIJINGpara(nParticles);
-  hijingparam->SetMomentumRange(0.2, 999);
-  gener->AddGenerator(hijingparam,"HIJING PARAM",1);
-
-  //    AliGenBox *genbox = new AliGenBox(nParticles);
-  //    genbox->SetPart(22);
-  //    genbox->SetPtRange(0.3, 10.00);
-  //    gener->AddGenerator(genbox,"GENBOX GAMMA for PHOS",1);
-  gener->Init();
+  // Set the trigger configuration
+  AliSimulation::Instance()->SetTriggerConfig("Pb-Pb");
+  cout<<"Trigger configuration is set to  Pb-Pb"<<endl;
 
-  // Activate this line if you want the vertex smearing to happen
-  // track by track
+    // Set Random Number seed
+    gRandom->SetSeed(123456); // Set 0 to use the currecnt time
+    AliLog::Message(AliLog::kInfo, Form("Seed for random number generation = %d",gRandom->GetSeed()), "Config.C", "Config.C", "Config()","Config.C", __LINE__);
+    int     nParticles = 100;
+    if (gSystem->Getenv("CONFIG_NPARTICLES"))
+    {
+        nParticles = atoi(gSystem->Getenv("CONFIG_NPARTICLES"));
+    }
 
-  //gener->SetVertexSmear(perTrack); 
+    AliGenCocktail *gener = new AliGenCocktail();
+    gener->SetPhiRange(0, 360);
+    // Set pseudorapidity range from -8 to 8.
+    Float_t thmin = EtaToTheta(8);   // theta min. <---> eta max
+    Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min 
+    gener->SetThetaRange(thmin,thmax);
+    gener->SetOrigin(0, 0, 0);  //vertex position
+    gener->SetSigma(0, 0, 0);   //Sigma in (X,Y,Z) (cm) on IP position
+
+    AliGenHIJINGpara *hijingparam = new AliGenHIJINGpara(nParticles);
+    hijingparam->SetMomentumRange(0.2, 999);
+    gener->AddGenerator(hijingparam,"HIJING PARAM",1);
+    gener->Init();
+
+    // 
+    // Activate this line if you want the vertex smearing to happen
+    // track by track
+    //
+    //gener->SetVertexSmear(perTrack); 
 
   // ============================= 
   // Magnetic field
   // ============================= 
 
-  // Field (L3 0.4 T)
-  AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., 1);
-  gAlice->SetField(field);    
-  
-  // Old magnetic field
-  //gAlice->SetField(-999,2);    //Specify maximum magnetic field in Tesla (neg. ==> default field)
-
-  // ============================= 
-  // Alice modules
-  // ============================= 
-
-  if (!interactiveSetup) {
-
-    // Select modules 
-  Int_t iABSO  = 1;
-  Int_t iDIPO  = 1;
-  Int_t iFMD   = 1;
-  Int_t iFRAME = 1;
-  Int_t iHALL  = 1;
-  Int_t iITS   = 1;
-  Int_t iMAG   = 1;
-  Int_t iMUON  = 1;
-  Int_t iPHOS  = 1;
-  Int_t iPIPE  = 1;
-  Int_t iPMD   = 1;
-  Int_t iRICH  = 1;
-  Int_t iSHIL  = 1;
-  Int_t iSTART = 1;
-  Int_t iTOF   = 1;
-  Int_t iTPC   = 1;
-  Int_t iTRD   = 1;
-  Int_t iZDC   = 1;
-  Int_t iEMCAL = 1;
-  Int_t iCRT   = 0;  
-  Int_t iVZERO = 1;
-
-  // ONLY FOR GEANT4
-
-  // Exclude detectors with temporary problem
-  iCRT = 0;
-  iEMCAL = 0;
-  // END OF ONLY FOR GEANT4
-
-
+    // Field (L3 0.4 T)
+    TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
+
+    Int_t   iABSO  =  1;
+    Int_t   iDIPO  =  1;
+    Int_t   iFMD   =  1;
+    Int_t   iFRAME =  1;
+    Int_t   iHALL  =  1;
+    Int_t   iITS   =  1;
+    Int_t   iMAG   =  1;
+    Int_t   iMUON  =  1;
+    Int_t   iPHOS  =  1;
+    Int_t   iPIPE  =  1;
+    Int_t   iPMD   =  1;
+    Int_t   iHMPID =  1;
+    Int_t   iSHIL  =  1;
+    Int_t   iT0    =  1;
+    Int_t   iTOF   =  1;
+    Int_t   iTPC   =  1;
+    Int_t   iTRD   =  1;
+    Int_t   iZDC   =  1;
+    Int_t   iEMCAL =  1;
+    Int_t   iACORDE = 0;
+    Int_t   iVZERO =  1;
+/*
+    Int_t   iABSO  =  0;
+    Int_t   iDIPO  =  0;
+    Int_t   iFMD   =  0;
+    Int_t   iFRAME =  0;
+    Int_t   iHALL  =  0;
+    Int_t   iITS   =  0;
+    Int_t   iMAG   =  0;
+    Int_t   iMUON  =  0;
+    Int_t   iPHOS  =  0;
+    Int_t   iPIPE  =  0;
+    Int_t   iPMD   =  0;
+    Int_t   iHMPID =  0;
+    Int_t   iSHIL  =  0;
+    Int_t   iT0    =  0;
+    Int_t   iTOF   =  0;
+    Int_t   iTPC   =  1;
+    Int_t   iTRD   =  0;
+    Int_t   iZDC   =  0;
+    Int_t   iEMCAL =  0;
+    Int_t   iACORDE = 0;
+    Int_t   iVZERO =  0;
+*/   
+// Exluded for problems
+
+    rl->CdGAFile();
     //=================== Alice BODY parameters =============================
     AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
 
-
     if (iMAG)
     {
         //=================== MAG parameters ============================
@@ -164,21 +144,21 @@ void ConfigCommon(Bool_t interactiveSetup)
     if (iABSO)
     {
         //=================== ABSO parameters ============================
-        AliABSO *ABSO = new AliABSOv0("ABSO", "Muon Absorber");
+        AliABSO *ABSO = new AliABSOv3("ABSO", "Muon Absorber");
     }
 
     if (iDIPO)
     {
         //=================== DIPO parameters ============================
 
-        AliDIPO *DIPO = new AliDIPOv2("DIPO", "Dipole version 2");
+        AliDIPO *DIPO = new AliDIPOv3("DIPO", "Dipole version 3");
     }
 
     if (iHALL)
     {
         //=================== HALL parameters ============================
 
-        AliHALL *HALL = new AliHALL("HALL", "Alice Hall");
+        AliHALL *HALL = new AliHALLv3("HALL", "Alice Hall");
     }
 
 
@@ -187,13 +167,18 @@ void ConfigCommon(Bool_t interactiveSetup)
         //=================== FRAME parameters ============================
 
         AliFRAMEv2 *FRAME = new AliFRAMEv2("FRAME", "Space Frame");
+        if (geo == kHoles) {
+         FRAME->SetHoles(1);
+       } else {
+         FRAME->SetHoles(0);
+       }
     }
 
     if (iSHIL)
     {
         //=================== SHIL parameters ============================
 
-        AliSHIL *SHIL = new AliSHILv2("SHIL", "Shielding Version 2");
+        AliSHIL *SHIL = new AliSHILv3("SHIL", "Shielding Version 3");
     }
 
 
@@ -201,63 +186,15 @@ void ConfigCommon(Bool_t interactiveSetup)
     {
         //=================== PIPE parameters ============================
 
-        AliPIPE *PIPE = new AliPIPEv0("PIPE", "Beam Pipe");
+        AliPIPE *PIPE = new AliPIPEv3("PIPE", "Beam Pipe");
     }
  
-    if(iITS) 
+    if (iITS)
     {
-      //=================== ITS parameters ============================
-      //
-      // As the innermost detector in ALICE, the Inner Tracking System "impacts" on
-      // almost all other detectors. This involves the fact that the ITS geometry
-      // still has several options to be followed in parallel in order to determine
-      // the best set-up which minimizes the induced background. All the geometries
-      // available to date are described in the following. Read carefully the comments
-      // and use the default version (the only one uncommented) unless you are making
-      // comparisons and you know what you are doing. In this case just uncomment the
-      // ITS geometry you want to use and run Aliroot.
-      //
-      // Detailed geometries:         
-      //
-      //
-      //
-/*
-       AliITSvPPRasymmFMD *ITS  = new AliITSvPPRasymmFMD("ITS","ITS PPR detailed version with asymmetric services");
-       ITS->SetMinorVersion(2);  // don't touch this parameter if you're not an ITS developer
-       ITS->SetReadDet(kTRUE);   // don't touch this parameter if you're not an ITS developer
-       //    ITS->SetWriteDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymm2.det");  // don't touch this parameter if you're not an ITS developer
-       ITS->SetThicknessDet1(200.);   // detector thickness on layer 1 must be in the range [100,300]
-       ITS->SetThicknessDet2(200.);   // detector thickness on layer 2 must be in the range [100,300]
-       ITS->SetThicknessChip1(200.);  // chip thickness on layer 1 must be in the range [150,300]
-       ITS->SetThicknessChip2(200.);  // chip thickness on layer 2 must be in the range [150,300]
-       ITS->SetRails(0);              // 1 --> rails in ; 0 --> rails out
-       ITS->SetCoolingFluid(1);       // 1 --> water ; 0 --> freon
-*/
-      //
-      // Coarse geometries (warning: no hits are produced with these coarse geometries and they unuseful 
-      // for reconstruction !):
-      //                                                       
-      //
-      AliITSvPPRcoarseasymm *ITS  = new AliITSvPPRcoarseasymm("ITS","New ITS PPR coarse version with asymmetric services");
-      ITS->SetRails(0);                // 1 --> rails in ; 0 --> rails out
-      ITS->SetSupportMaterial(0);      // 0 --> Copper ; 1 --> Aluminum ; 2 --> Carbon
-      //
-      //AliITS *ITS  = new AliITSvPPRcoarsesymm("ITS","New ITS PPR coarse version with symmetric services");
-      //ITS->SetRails(0);                // 1 --> rails in ; 0 --> rails out
-      //ITS->SetSupportMaterial(0);      // 0 --> Copper ; 1 --> Aluminum ; 2 --> Carbon
-      //                      
-      //
-      //
-      // Geant3 <-> EUCLID conversion
-      // ============================
-      //
-      // SetEUCLID is a flag to output (=1) or not to output (=0) both geometry and
-      // media to two ASCII files (called by default ITSgeometry.euc and
-      // ITSgeometry.tme) in a format understandable to the CAD system EUCLID.
-      // The default (=0) means that you dont want to use this facility.
-      //
-      ITS->SetEUCLID(0);  
-    }  
+        //=================== ITS parameters ============================
+
+       AliITS *ITS  = new AliITSv11("ITS","ITS v11");
+    }
 
     if (iTPC)
     {
@@ -265,16 +202,18 @@ void ConfigCommon(Bool_t interactiveSetup)
         AliTPC *TPC = new AliTPCv2("TPC", "Default");
     }
 
-    if (iTOF)
-    {
+
+    if (iTOF) {
         //=================== TOF parameters ============================
-       AliTOF *TOF = new AliTOFv5T0("TOF", "normal TOF");
+       AliTOF *TOF = new AliTOFv6T0("TOF", "normal TOF");
     }
 
-    if (iRICH)
+
+    if (iHMPID)
     {
-        //=================== RICH parameters ===========================
-        AliRICH *RICH = new AliRICHv1("RICH", "normal RICH");
+        //=================== HMPID parameters ===========================
+        AliHMPID *HMPID = new AliHMPIDv3("HMPID", "normal HMPID");
+
     }
 
 
@@ -282,7 +221,7 @@ void ConfigCommon(Bool_t interactiveSetup)
     {
         //=================== ZDC parameters ============================
 
-        AliZDC *ZDC = new AliZDCv2("ZDC", "normal ZDC");
+        AliZDC *ZDC = new AliZDCv3("ZDC", "normal ZDC");
     }
 
     if (iTRD)
@@ -290,24 +229,18 @@ void ConfigCommon(Bool_t interactiveSetup)
         //=================== TRD parameters ============================
 
         AliTRD *TRD = new AliTRDv1("TRD", "TRD slow simulator");
-
-        // Select the gas mixture (0: 97% Xe + 3% isobutane, 1: 90% Xe + 10% CO2)
-        TRD->SetGasMix(1);
-        // Switch on TR
-        AliTRDsim *TRDsim = TRD->CreateTR();
     }
 
     if (iFMD)
     {
         //=================== FMD parameters ============================
-
-        AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
-    }
+       AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
+   }
 
     if (iMUON)
     {
         //=================== MUON parameters ===========================
-
+        // New MUONv1 version (geometry defined via builders)
         AliMUON *MUON = new AliMUONv1("MUON", "default");
     }
     //=================== PHOS parameters ===========================
@@ -324,36 +257,45 @@ void ConfigCommon(Bool_t interactiveSetup)
         AliPMD *PMD = new AliPMDv1("PMD", "normal PMD");
     }
 
-    if (iSTART)
+    if (iT0)
     {
-        //=================== START parameters ============================
-        AliSTART *START = new AliSTARTv1("START", "START Detector");
+        //=================== T0 parameters ============================
+        AliT0 *T0 = new AliT0v1("T0", "T0 Detector");
     }
 
     if (iEMCAL)
     {
         //=================== EMCAL parameters ============================
-        AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "SHISH");
+        AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETE");
     }
 
-    if (iCRT)
+     if (iACORDE)
     {
-        //=================== CRT parameters ============================
-
-        AliCRT *CRT = new AliCRTv0("CRT", "normal ACORDE");
+        //=================== ACORDE parameters ============================
+        AliACORDE *ACORDE = new AliACORDEv1("ACORDE", "normal ACORDE");
     }
 
-    if (iVZERO)
+     if (iVZERO)
     {
-        //=================== CRT parameters ============================
-        AliVZERO *VZERO = new AliVZEROv6("VZERO", "normal VZERO");
+        //=================== ACORDE parameters ============================
+        AliVZERO *VZERO = new AliVZEROv7("VZERO", "normal VZERO");
     }
 
-  } // end (!isSetInteractively)
+     AliLog::Message(AliLog::kInfo, "End of Config", "Config.C", "Config.C", "Config()"," Config.C", __LINE__);
+
+     cout << "Running ConfigCommon.C finished ... " << endl;
 
-  cout << "End of g4ConfigCommon.C" << endl;
 }
 
 Float_t EtaToTheta(Float_t arg){
   return (180./TMath::Pi())*2.*atan(exp(-arg));
 }
+
+void LoadPythia()
+{
+    // Load Pythia related libraries
+    gSystem->Load("liblhapdf.so");      // Parton density functions
+    gSystem->Load("libEGPythia6.so");   // TGenerator interface
+    gSystem->Load("libpythia6.so");     // Pythia
+    gSystem->Load("libAliPythia6.so");  // ALICE specific implementations
+}