Partial TOF and TRD geometry as described by Karel in his message to the physics...

[u/mrichter/AliRoot.git] / macros / Config_PDC06.C

diff --git a/macros/Config_PDC06.C b/macros/Config_PDC06.C
index 0841e529b0bffeaa47519cf4765a85de5cbdaed0..dd551519de64dfb10ee14c8e61b9bba937361878 100644 (file)
--- a/macros/Config_PDC06.C
+++ b/macros/Config_PDC06.C
@@ -46,7 +46,7 @@
 #include "START/AliSTARTv1.h"
 #include "EMCAL/AliEMCALv2.h"
 #include "CRT/AliCRTv0.h"
 #include "START/AliSTARTv1.h"
 #include "EMCAL/AliEMCALv2.h"
 #include "CRT/AliCRTv0.h"

-#include "VZERO/AliVZEROv6.h"
+#include "VZERO/AliVZEROv7.h"

 #endif
 
 
 #endif
 
 


@@ -85,6 +85,17 @@ enum Mag_t
 {
     k2kG, k4kG, k5kG
 };
 {
     k2kG, k4kG, k5kG
 };

+
+//--- Trigger config ---
+enum TrigConf_t
+{
+    kDefaultPPTrig, kDefaultPbPbTrig
+};
+
+const char * TrigConfName[] = {
+    "p-p","Pb-Pb"
+};
+

 //--- Functions ---
 AliGenPythia *PythiaHVQ(PDC06Proc_t proc);
 AliGenerator *MbCocktail();
 //--- Functions ---
 AliGenPythia *PythiaHVQ(PDC06Proc_t proc);
 AliGenerator *MbCocktail();


@@ -96,6 +107,7 @@ static PDC06Proc_t   proc     = kPyOmegaPlus;
 static DecayHvFl_t   decHvFl  = kNature; 
 static YCut_t        ycut     = kFull;
 static Mag_t         mag      = k5kG; 
 static DecayHvFl_t   decHvFl  = kNature; 
 static YCut_t        ycut     = kFull;
 static Mag_t         mag      = k5kG; 

+static TrigConf_t    trig     = kDefaultPPTrig; // default pp trigger configuration

 //========================//
 // Set Random Number seed //
 //========================//
 //========================//
 // Set Random Number seed //
 //========================//


@@ -168,6 +180,10 @@ void Config()
   rl->SetCompressionLevel(2);
   rl->SetNumberOfEventsPerFile(1000);
   gAlice->SetRunLoader(rl);
   rl->SetCompressionLevel(2);
   rl->SetNumberOfEventsPerFile(1000);
   gAlice->SetRunLoader(rl);

+  
+  // Set the trigger configuration
+  gAlice->SetTriggerDescriptor(TrigConfName[trig]);
+  cout<<"Trigger configuration is set to  "<<TrigConfName[trig]<<endl;

 
   //
   //=======================================================================
 
   //
   //=======================================================================


@@ -306,7 +322,7 @@ void Config()
   printf("\n \n Diamond size x-y: %10.3e z: %10.3e\n \n", sigmaxy, sigmaz);
     
   gener->SetSigma(sigmaxy, sigmaxy, sigmaz);      // Sigma in (X,Y,Z) (cm) on IP position
   printf("\n \n Diamond size x-y: %10.3e z: %10.3e\n \n", sigmaxy, sigmaz);
     
   gener->SetSigma(sigmaxy, sigmaxy, sigmaz);      // Sigma in (X,Y,Z) (cm) on IP position

-  gener->SetCutVertexZ(1.);        // Truncate at 1 sigma
+  gener->SetCutVertexZ(3.);        // Truncate at 3 sigma

   gener->SetVertexSmear(kPerEvent);
 
   gener->Init();
   gener->SetVertexSmear(kPerEvent);
 
   gener->Init();


@@ -330,7 +346,7 @@ void Config()
 
 
   Int_t iABSO  = 1;
 
 
   Int_t iABSO  = 1;

-  Int_t iCRT   = 1;
+  Int_t iCRT   = 0;

   Int_t iDIPO  = 1;
   Int_t iEMCAL = 1;
   Int_t iFMD   = 1;
   Int_t iDIPO  = 1;
   Int_t iEMCAL = 1;
   Int_t iFMD   = 1;


@@ -474,6 +490,10 @@ void Config()
     if (iTOF) {
         //=================== TOF parameters ============================
        AliTOF *TOF = new AliTOFv5T0("TOF", "normal TOF");
     if (iTOF) {
         //=================== TOF parameters ============================
        AliTOF *TOF = new AliTOFv5T0("TOF", "normal TOF");

+       // Partial geometry: modules at 2,3,4,6,7,11,12,14,15,16
+       // starting at 6h in positive direction
+       Int_t TOFSectors[18]={-1,-1,0,0,0,-1,0,0,-1,-1,-1,0,0,-1,0,0,0,-1};
+       TOF->SetTOFSectors(TOFSectors);

     }
 
 
     }
 
 


@@ -497,11 +517,19 @@ void Config()
         //=================== TRD parameters ============================
 
         AliTRD *TRD = new AliTRDv1("TRD", "TRD slow simulator");
         //=================== 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();
+        AliTRDgeometry *geoTRD = TRD->GetGeometry();
+       // Partial geometry: modules at 2,3,4,6,12,14,15,16
+       // starting at 6h in positive direction
+       geoTRD->SetSMstatus(0,0);
+        geoTRD->SetSMstatus(1,0);
+        geoTRD->SetSMstatus(5,0);
+        geoTRD->SetSMstatus(7,0);
+        geoTRD->SetSMstatus(8,0);
+        geoTRD->SetSMstatus(9,0);
+        geoTRD->SetSMstatus(10,0);
+        geoTRD->SetSMstatus(11,0);
+        geoTRD->SetSMstatus(13,0);
+        geoTRD->SetSMstatus(17,0);

     }
 
     if (iFMD)
     }
 
     if (iFMD)


@@ -539,7 +567,7 @@ void Config()
     if (iEMCAL)
     {
         //=================== EMCAL parameters ============================
     if (iEMCAL)
     {
         //=================== EMCAL parameters ============================

-        AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "SHISH");
+        AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "SHISH_77_TRD1_2X2_FINAL_110DEG");

     }
 
      if (iCRT)
     }
 
      if (iCRT)


@@ -551,7 +579,7 @@ void Config()
      if (iVZERO)
     {
         //=================== CRT parameters ============================
      if (iVZERO)
     {
         //=================== CRT parameters ============================

-        AliVZERO *VZERO = new AliVZEROv6("VZERO", "normal VZERO");
+        AliVZERO *VZERO = new AliVZEROv7("VZERO", "normal VZERO");

     }
 }
 //
     }
 }
 //