Changed call to TClonesArray::GetEntries() to TClonesArray::GetEntriesFast() for...

[u/mrichter/AliRoot.git] / FMD / Config.C

diff --git a/FMD/Config.C b/FMD/Config.C
index 32f1c10311573c0fbacda46eda8e1735730c44d7..71a5925faa003046dfc14d3ad3fa4c714f8f5bde 100644 (file)
--- a/FMD/Config.C
+++ b/FMD/Config.C
@@ -15,7 +15,7 @@
 */
 #if !defined(__CINT__) || defined(__MAKECINT__)
 #include <Riostream.h>
-#include <TRandom.h>
+#include <TPDGCode.h>
 #include <TSystem.h>
 #include <TVirtualMC.h>
 #include <TGeant3.h>
@@ -45,7 +45,7 @@
 #include "ITS/AliITSvPPRasymmFMD.h"
 #include "TPC/AliTPCv2.h"
 #include "TOF/AliTOFv4T0.h"
-#include "RICH/AliRICHv1.h"
+#include "HMPID/AliHMPIDv1.h"
 #include "ZDC/AliZDCv2.h"
 #include "TRD/AliTRDv1.h"
 #include "FMD/AliFMDv1.h"
@@ -56,9 +56,9 @@
 #include "MUON/AliMUONTriggerGeometryBuilder.h"
 #include "PHOS/AliPHOSv1.h"
 #include "PMD/AliPMDv1.h"
-#include "START/AliSTARTv1.h"
+#include "T0/AliT0v1.h"
 #include "EMCAL/AliEMCALv1.h"
-#include "CRT/AliCRTv0.h"
+#include "ACORDE/AliACORDEv0.h"
 #include "VZERO/AliVZEROv2.h"
 #endif
 
@@ -228,9 +228,10 @@ Config()
 {
   //____________________________________________________________________
   // This part for configuration    
-  //static EG_t eg = test50;
-  //EG_t  eg   = kParam_fmd;
+  // EG_t  eg   = test50;
+  // EG_t  eg   = kParam_fmd;
   EG_t  eg   = kParam_2000; // kPythia;
+  // EG_t  eg   = kFMDFlat;
   Geo_t geo  = kNoHoles;
   Rad_t rad  = kGluonRadiation;
   Mag_t mag  = k5kG;
@@ -286,6 +287,7 @@ Config()
       //
       // Libraries needed by GEANT 3.21 
       //
+      gSystem->Load("EGPythia6.so");
       gSystem->Load("libgeant321");
     
       // 
@@ -385,7 +387,7 @@ Config()
   gMC->SetProcess("MUNU",1);
   gMC->SetProcess("CKOV",1);
   gMC->SetProcess("HADR",1);
-  gMC->SetProcess("LOSS",2);
+  gMC->SetProcess("LOSS",2); // 0:none 1,3:dray 2:nodray 4:nofluct (def:2)
   gMC->SetProcess("MULS",1);
   gMC->SetProcess("RAYL",1);
 
@@ -433,8 +435,8 @@ Config()
   }
 
   switch(geo) {
-  case kHoles: comment = comment.Append(" | Holes for PHOS/RICH"); break;
-  default:     comment = comment.Append(" | No holes for PHOS/RICH"); break;
+  case kHoles: comment = comment.Append(" | Holes for PHOS/HMPID"); break;
+  default:     comment = comment.Append(" | No holes for PHOS/HMPID"); break;
   }
 
   std::cout << "\n\n Comment: " << comment << "\n" << std::endl;
@@ -451,7 +453,7 @@ Config()
   // Used detectors 
   // 
   Bool_t useABSO  = kFALSE; 
-  Bool_t useCRT   = kFALSE; 
+  Bool_t useACORDE   = kFALSE; 
   Bool_t useDIPO  = kFALSE; 
   Bool_t useFMD   = kTRUE; 
   Bool_t useFRAME = kFALSE; 
@@ -462,9 +464,9 @@ Config()
   Bool_t usePHOS  = kFALSE; 
   Bool_t usePIPE  = kFALSE; 
   Bool_t usePMD   = kFALSE; 
-  Bool_t useRICH  = kFALSE; 
+  Bool_t useHMPID  = kFALSE; 
   Bool_t useSHIL  = kFALSE; 
-  Bool_t useSTART = kFALSE; 
+  Bool_t useT0 = kFALSE; 
   Bool_t useTOF   = kFALSE; 
   Bool_t useTPC   = kFALSE;
   Bool_t useTRD   = kFALSE; 
@@ -629,9 +631,9 @@ Config()
     AliTOF *TOF = new AliTOFv4T0("TOF", "normal TOF");
   }
 
-  if (useRICH) {
-    // ================== RICH parameters ============================
-    AliRICH *RICH = new AliRICHv1("RICH", "normal RICH");
+  if (useHMPID) {
+    // ================== HMPID parameters ============================
+    AliHMPID *HMPID = new AliHMPIDv1("HMPID", "normal HMPID");
 
   }
 
@@ -649,8 +651,8 @@ Config()
     if (geo == kHoles) {
       // With hole in front of PHOS
       TRD->SetPHOShole();
-      // With hole in front of RICH
-      TRD->SetRICHhole();
+      // With hole in front of HMPID
+      TRD->SetHMPIDhole();
     }
     // Switch on TR
     AliTRDsim *TRDsim = TRD->CreateTR();
@@ -659,7 +661,6 @@ Config()
   if (useFMD) {
     // =================== FMD parameters ============================
     AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
-    AliLog::SetModuleDebugLevel("FMD", 1);
     // FMD->UseDetailed(kFALSE);
     // FMD->UseAssembly();
     // FMD->UseOld();
@@ -684,9 +685,9 @@ Config()
     AliPMD *PMD = new AliPMDv1("PMD", "normal PMD");
   }
 
-  if (useSTART) {
-    // =================== START parameters ==========================
-    AliSTART *START = new AliSTARTv1("START", "START Detector");
+  if (useT0) {
+    // =================== T0 parameters ==========================
+    AliT0 *T0 = new AliT0v1("T0", "T0 Detector");
   }
 
   if (useEMCAL) {
@@ -694,9 +695,9 @@ Config()
     AliEMCAL *EMCAL = new AliEMCALv1("EMCAL", "EMCAL_55_25");
   }
 
-  if (useCRT) {
-    // =================== CRT parameters ============================
-    AliCRT *CRT = new AliCRTv0("CRT", "normal ACORDE");
+  if (useACORDE) {
+    // =================== ACORDE parameters ============================
+    AliACORDE *ACORDE = new AliACORDEv0("ACORDE", "normal ACORDE");
   }
 
   if (useVZERO) {
@@ -1581,7 +1582,7 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     {
       comment = comment.Append(" Flat in FMD1 range");
       AliGenBox* gener = new AliGenBox(2000);
-      gener->SetPart(211);
+      gener->SetPart(kPiPlus);
       gener->SetMomentumRange(3,4);
       gener->SetPhiRange(0, 360);
       gener->SetThetaRange(0.77, 3.08);
@@ -1592,7 +1593,7 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     {
       comment = comment.Append(" Flat in FMD2 range");
       AliGenBox* gener = new AliGenBox(2000);
-      gener->SetPart(211);
+      gener->SetPart(kPiPlus);
       gener->SetMomentumRange(3,4);
       gener->SetPhiRange(0, 360);
       gener->SetThetaRange(2.95, 20.42);
@@ -1603,7 +1604,7 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     {
       comment = comment.Append(" Flat in FMD3 range");
       AliGenBox* gener = new AliGenBox(2000);
-      gener->SetPart(211);
+      gener->SetPart(kPiPlus);
       gener->SetMomentumRange(3,4);
       gener->SetPhiRange(0, 360);
       gener->SetThetaRange(155.97, 176.73);
@@ -1613,18 +1614,20 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
   case kFMDFlat:
     {
       comment = comment.Append(" Flat in FMD range");
-      AliGenCocktail* gener = AliGenCocktail("FMD cocktail");
-      gener->SetPart(211);
+      AliGenCocktail* gener = new AliGenCocktail();
       gener->SetMomentumRange(3,4);
       gener->SetPhiRange(0, 360);
       AliGenBox* gener3 = new AliGenBox(2000);
       gener3->SetThetaRange(155.97, 176.73);
+      gener3->SetPart(kPiPlus);
       gener->AddGenerator(gener3, "FMD3", .33);
       AliGenBox* gener2 = new AliGenBox(2000);
       gener2->SetThetaRange(2.95, 20.42);
+      gener2->SetPart(kPiPlus);
       gener->AddGenerator(gener2, "FMD2", .33);
       AliGenBox* gener1 = new AliGenBox(2000);
       gener1->SetThetaRange(0.77, 3.08);
+      gener1->SetPart(kPiPlus);
       gener->AddGenerator(gener1, "FMD1", .34);
       gGener = gener;
     }