]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - macros/Config.C
dding blank lines
[u/mrichter/AliRoot.git] / macros / Config.C
index 28e80da27b0d3ab3f48575e88167a06df0891bc3..dcee82aebecbfc3e4fa3614c0b69b182b5e23ce1 100644 (file)
@@ -6,6 +6,7 @@
 
 #if !defined(__CINT__) || defined(__MAKECINT__)
 #include <Riostream.h>
+#include <TPDGCode.h>
 #include <TRandom.h>
 #include <TSystem.h>
 #include <TVirtualMC.h>
 #include "STRUCT/AliFRAMEv2.h"
 #include "STRUCT/AliSHILv3.h"
 #include "STRUCT/AliPIPEv3.h"
-#include "ITS/AliITSvPPRasymmFMD.h"
+#include "ITS/AliITSv11Hybrid.h"
 #include "TPC/AliTPCv2.h"
 #include "TOF/AliTOFv6T0.h"
-#include "HMPID/AliHMPIDv2.h"
-#include "ZDC/AliZDCv2.h"
+#include "HMPID/AliHMPIDv3.h"
+#include "ZDC/AliZDCv3.h"
 #include "TRD/AliTRDv1.h"
 #include "FMD/AliFMDv1.h"
 #include "MUON/AliMUONv1.h"
 #include "VZERO/AliVZEROv7.h"
 #endif
 
-enum PprGeo_t
-  {
-    kHoles, kNoHoles
-  };
-static PprGeo_t geo = kHoles;
-
 Float_t EtaToTheta(Float_t arg);
+void    LoadPythia();
+
 
 void Config()
 {
@@ -55,22 +52,19 @@ void Config()
     // Theta range given through pseudorapidity limits 22/6/2001
 
     // Set Random Number seed
-  gRandom->SetSeed(123456); // Set 0 to use the currecnt time
+  //gRandom->SetSeed(123456); // Set 0 to use the current time
+  
   AliLog::Message(AliLog::kInfo, Form("Seed for random number generation = %d",gRandom->GetSeed()), "Config.C", "Config.C", "Config()","Config.C", __LINE__);
 
-
-   // libraries required by geant321
+  // Load Pythia libraries
+  LoadPythia();
+  // Libraries required by geant321
 #if defined(__CINT__)
     gSystem->Load("libgeant321");
 #endif
 
     new     TGeant3TGeo("C++ Interface to Geant3");
 
-    if(!AliCDBManager::Instance()->IsDefaultStorageSet()){
-      AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT");
-      AliCDBManager::Instance()->SetRun(0);
-    }
-
     AliRunLoader* rl=0x0;
 
     AliLog::Message(AliLog::kInfo, "Creating Run Loader", "Config.C", "Config.C", "Config()"," Config.C", __LINE__);
@@ -86,8 +80,35 @@ void Config()
     rl->SetCompressionLevel(2);
     rl->SetNumberOfEventsPerFile(3);
     gAlice->SetRunLoader(rl);
+    
     // gAlice->SetGeometryFromFile("geometry.root");
-    // gAlice->SetGeometryFromCDB();
+
+    // Uncomment if you want to load geometry from OCDB!   >>>>
+/*    
+    if(!AliCDBManager::Instance()->IsDefaultStorageSet()){
+        cout << "#####################################################" << endl;
+        cout << "#                                                   #" << endl;
+        cout << "#     WARNING: CDB DEFAULT STORAGE NOT SET !!!      #" << endl;
+        cout << "#     SETTING IT TO local://$ALICE_ROOT !!!         #" << endl;
+        cout << "#                                                   #" << endl;
+        cout << "#####################################################" << endl;
+          
+         AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT");
+    }
+    
+    if(AliCDBManager::Instance()->GetRun() < 0){
+        cout << "#####################################################" << endl;
+        cout << "#                                                   #" << endl;
+        cout << "#     WARNING: RUN NUMBER NOT SET !!!               #" << endl;
+        cout << "#     SETTING IT TO 0 !!!                           #" << endl;
+        cout << "#                                                   #" << endl;
+        cout << "#####################################################" << endl;
+          
+         AliCDBManager::Instance()->SetRun(0);
+    }
+    gAlice->SetGeometryFromCDB();
+*/
+    // Uncomment if you want to load geometry from OCDB!   <<<<
 
     // Set the trigger configuration
     gAlice->SetTriggerDescriptor("Pb-Pb");
@@ -157,7 +178,7 @@ void Config()
     gener->AddGenerator(hijingparam,"HIJING PARAM",1);
 
 //    AliGenBox *genbox = new AliGenBox(nParticles);
-//    genbox->SetPart(22);
+//    genbox->SetPart(kGamma);
 //    genbox->SetPtRange(0.3, 10.00);
 //    gener->AddGenerator(genbox,"GENBOX GAMMA for PHOS",1);
     gener->Init();
@@ -233,11 +254,7 @@ void Config()
         //=================== FRAME parameters ============================
 
         AliFRAMEv2 *FRAME = new AliFRAMEv2("FRAME", "Space Frame");
-        if (geo == kHoles) {
-         FRAME->SetHoles(1);
-       } else {
-         FRAME->SetHoles(0);
-       }
+       FRAME->SetHoles(1);
     }
 
     if (iSHIL)
@@ -259,7 +276,7 @@ void Config()
     {
         //=================== ITS parameters ============================
 
-       AliITSvPPRasymmFMD *ITS  = new AliITSvPPRasymmFMD("ITS","ITS PPR detailed version with asymmetric services");
+       AliITS *ITS  = new AliITSv11Hybrid("ITS","ITS v11Hybrid");
     }
 
     if (iTPC)
@@ -278,7 +295,7 @@ void Config()
     if (iHMPID)
     {
         //=================== HMPID parameters ===========================
-        AliHMPID *HMPID = new AliHMPIDv2("HMPID", "normal HMPID");
+        AliHMPID *HMPID = new AliHMPIDv3("HMPID", "normal HMPID");
 
     }
 
@@ -287,7 +304,7 @@ void Config()
     {
         //=================== ZDC parameters ============================
 
-        AliZDC *ZDC = new AliZDCv2("ZDC", "normal ZDC");
+        AliZDC *ZDC = new AliZDCv3("ZDC", "normal ZDC");
     }
 
     if (iTRD)
@@ -332,7 +349,7 @@ void Config()
     if (iEMCAL)
     {
         //=================== EMCAL parameters ============================
-        AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "SHISH_77_TRD1_2X2_FINAL_110DEG");
+        AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETE");
     }
 
      if (iACORDE)
@@ -354,3 +371,13 @@ void Config()
 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
+}