STEER/STEER: Install AliRawDataHeaderSim.h for HMPID AMORE module
[u/mrichter/AliRoot.git] / macros / ConfigBeautyPPR.C
index 68d48ed..4a0a847 100644 (file)
 #include "STEER/AliMagFCM.h"
 #include "STRUCT/AliBODY.h"
 #include "STRUCT/AliMAG.h"
-#include "STRUCT/AliABSOv0.h"
-#include "STRUCT/AliDIPOv2.h"
-#include "STRUCT/AliHALL.h"
+#include "STRUCT/AliABSOv3.h"
+#include "STRUCT/AliDIPOv3.h"
+#include "STRUCT/AliHALLv3.h"
 #include "STRUCT/AliFRAMEv2.h"
-#include "STRUCT/AliSHILv2.h"
-#include "STRUCT/AliPIPEv0.h"
-#include "ITS/AliITSvPPRasymm.h"
+#include "STRUCT/AliSHILv3.h"
+#include "STRUCT/AliPIPEv3.h"
+#include "ITS/AliITSv11.h"
 #include "TPC/AliTPCv2.h"
-#include "TOF/AliTOFv5T0.h"
-#include "RICH/AliRICHv1.h"
-#include "ZDC/AliZDCv1.h"
+#include "TOF/AliTOFv6T0.h"
+#include "HMPID/AliHMPIDv3.h"
+#include "ZDC/AliZDCv3.h"
 #include "TRD/AliTRDv1.h"
 #include "FMD/AliFMDv0.h"
 #include "MUON/AliMUONv1.h"
 #include "PHOS/AliPHOSv1.h"
 #include "PMD/AliPMDv1.h"
-#include "START/AliSTARTv1.h"
-#include "CRT/AliCRTv1.h"
+#include "T0/AliT0v1.h"
+#include "ACORDE/AliACORDEv1.h"
 #endif
 
+void LoadPythia();
+
 void Config()
 {
  
@@ -46,7 +48,8 @@ void Config()
   //  gRandom->SetSeed(seed);
   gRandom->SetSeed(12345);
   cerr<<"Seed for random number generation= "<<seed<<endl; 
-
+  // Load Pythia libraries
+  LoadPythia();
   // libraries required by geant321
 #if defined(__CINT__)
   gSystem->Load("libgeant321");
@@ -57,6 +60,11 @@ void Config()
   //=======================================================================
   //  Create the output file
    
+  if(!AliCDBManager::Instance()->IsDefaultStorageSet()){
+    AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
+    AliCDBManager::Instance()->SetRun(0);
+  }
+  
   AliRunLoader* rl=0x0;
 
   cout<<"Config.C: Creating Run Loader ..."<<endl;
@@ -71,9 +79,11 @@ void Config()
   rl->SetCompressionLevel(2);
   rl->SetNumberOfEventsPerFile(3);
   gAlice->SetRunLoader(rl);
+  // gAlice->SetGeometryFromFile("geometry.root");
+  // gAlice->SetGeometryFromCDB();
 
   // Set the trigger configuration
-  gAlice->SetTriggerDescriptor("Pb-Pb");
+  AliSimulation::Instance()->SetTriggerConfig("Pb-Pb");
   cout<<"Trigger configuration is set to  Pb-Pb"<<endl;
 
   //
@@ -161,7 +171,7 @@ void Config()
   gAlice->SetField(field);    
 
   Int_t iABSO=0;
-  Int_t iCRT=0;
+  Int_t iACORDE=0;
   Int_t iDIPO=0;
   Int_t iFMD=0;
   Int_t iFRAME=0;
@@ -172,9 +182,9 @@ void Config()
   Int_t iPHOS=0;
   Int_t iPIPE=0;
   Int_t iPMD=0;
-  Int_t iRICH=0;
+  Int_t iHMPID=0;
   Int_t iSHIL=0;
-  Int_t iSTART=0;
+  Int_t iT0=0;
   Int_t iTOF=0;
   Int_t iTPC=0;
   Int_t iTRD=0;
@@ -193,19 +203,19 @@ void Config()
 
   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");
   }
 
 
@@ -213,77 +223,27 @@ void Config()
     //=================== FRAME parameters ============================
 
     AliFRAME *FRAME  = new AliFRAMEv2("FRAME","Space Frame");
-
+    FRAME->SetHoles(1);
   }
 
   if(iSHIL) {
     //=================== SHIL parameters ============================
 
-    AliSHIL *SHIL  = new AliSHILv2("SHIL","Shielding");
+    AliSHIL *SHIL  = new AliSHILv3("SHIL","Shielding");
   }
 
 
   if(iPIPE) {
     //=================== PIPE parameters ============================
 
-    AliPIPE *PIPE  = new AliPIPEv0("PIPE","Beam Pipe");
+    AliPIPE *PIPE  = new AliPIPEv3("PIPE","Beam Pipe");
   }
 
 
   if(iITS) {
+    //=================== ITS parameters ============================
 
-//=================== 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:         
-    //
-    //
-    //AliITS *ITS  = new AliITSv5symm("ITS","Updated ITS TDR detailed version with symmetric services");
-    //
-    //AliITS *ITS  = new AliITSv5asymm("ITS","Updates ITS TDR detailed version with asymmetric services");
-    //
-    AliITSvPPRasymm *ITS  = new AliITSvPPRasymm("ITS","New ITS PPR detailed version with asymmetric services");
-    ITS->SetMinorVersion(2);                                         // don't touch this parameter if you're not an ITS developer
-    ITS->SetReadDet(kFALSE);                                         // 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 [150,300]
-    ITS->SetThicknessDet2(200.);   // detector thickness on layer 2 must be in the range [150,300]
-    ITS->SetThicknessChip1(150.);  // chip thickness on layer 1 must be in the range [100,300]
-    ITS->SetThicknessChip2(150.);  // chip thickness on layer 2 must be in the range [100,300]
-    ITS->SetRails(1);             // 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(1);                // 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(1);                // 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);  
+    AliITS *ITS  = new AliITSv11("ITS","ITS v11");
   }
   
 
@@ -295,12 +255,12 @@ void Config()
 
   if(iTOF) {
     //=================== TOF parameters ============================
-    AliTOF *TOF  = new AliTOFv5T0("TOF","normal TOF");
+    AliTOF *TOF  = new AliTOFv6T0("TOF","normal TOF");
   }
 
-  if(iRICH) {
-    //=================== RICH parameters ===========================
-    AliRICH *RICH  = new AliRICHv1("RICH","normal RICH");    
+  if(iHMPID) {
+    //=================== HMPID parameters ===========================
+    AliHMPID *HMPID  = new AliHMPIDv3("HMPID","normal HMPID");    
 
   }
 
@@ -308,13 +268,13 @@ void Config()
   if(iZDC) {
     //=================== ZDC parameters ============================
 
-    AliZDC *ZDC  = new AliZDCv1("ZDC","normal ZDC");
+    AliZDC *ZDC  = new AliZDCv3("ZDC","normal ZDC");
   }
 
-  if(iCRT) {
-    //=================== CRT parameters ============================
+  if(iACORDE) {
+    //=================== ACORDE parameters ============================
 
-    AliCRT *CRT  = new AliCRTv1("CRT","normal CRT");
+    AliACORDE *ACORDE  = new AliACORDEv1("ACORDE","normal ACORDE");
   }
 
   if(iTRD) {
@@ -352,10 +312,20 @@ void Config()
 
   }
 
-  if(iSTART) {
-    //=================== START parameters ============================
-    AliSTART *START  = new AliSTARTv1("START","START Detector");
+  if(iT0) {
+    //=================== T0 parameters ============================
+    AliT0 *T0  = new AliT0v1("T0","T0 Detector");
   }
 
          
 }
+
+
+void LoadPythia()
+{
+    // Load Pythia related libraries
+    gSystem->Load("liblhapdf");      // Parton density functions
+    gSystem->Load("libEGPythia6");   // TGenerator interface
+    gSystem->Load("libpythia6");     // Pythia
+    gSystem->Load("libAliPythia6");  // ALICE specific implementations
+}