Setting a trigger configuration (C.Cheshkov)

diff --git a/macros/Config_PythiaHeavyFlavours.C b/macros/Config_PythiaHeavyFlavours.C
index d28b551863934aa62c1544b44edca642ac808d0f..177e075c0fd04b30c7a95714b8786d176550a5e6 100644 (file)
--- a/macros/Config_PythiaHeavyFlavours.C
+++ b/macros/Config_PythiaHeavyFlavours.C
@@ -17,7 +17,7 @@
 #include <TDatime.h>
 #include <TSystem.h>
 #include <TVirtualMC.h>
-#include <TGeant3.h>
+#include <TGeant3TGeo.h>
 #include "STEER/AliRunLoader.h"
 #include "STEER/AliRun.h"
 #include "STEER/AliConfig.h"
@@ -70,6 +70,14 @@ enum Mag_t
 {
     k2kG, k4kG, k5kG
 };
+//--- Trigger config ---
+enum TrigConf_t
+{
+    kDefaultPPTrig, kDefaultPbPbTrig
+};
+const char * TrigConfName[] = {
+    "p-p","Pb-Pb"
+};
 //--- Functions ---
 AliGenPythia *PythiaHVQ(ProcessHvFl_t proc);
 
@@ -79,6 +87,7 @@ static ProcessHvFl_t procHvFl = kCharmPbPb5500;
 static DecayHvFl_t   decHvFl  = kNature; 
 static YCut_t        ycut     = kFull;
 static Mag_t         mag      = k5kG; 
+static TrigConf_t    trig     = kDefaultPbPbTrig; // default PbPb trigger configuration
 // nEvts = -1  : you get 1 QQbar pair and all the fragmentation and 
 //               decay chain
 // nEvts = N>0 : you get N charm / beauty Hadrons 
@@ -152,6 +161,10 @@ void Config()
   rl->SetNumberOfEventsPerFile(3);
   gAlice->SetRunLoader(rl);
 
+  // Set the trigger configuration
+  gAlice->SetTriggerDescriptor(TrigConfName[trig]);
+  cout<<"Trigger configuration is set to  "<<TrigConfName[trig]<<endl;
+
   //
   //=======================================================================
   // ************* STEERING parameters FOR ALICE SIMULATION **************