Update of calibration task and replace AliTRDrunCalib.C by AddTaskTRDCalib.C (Raphaelle)
[u/mrichter/AliRoot.git] / TRD / Macros / AliTRDConfigPID.C
index 3b40b63..7462530 100644 (file)
@@ -14,6 +14,7 @@
 
 #if !defined(__CINT__) || defined(__MAKECINT__)
 #include <Riostream.h>
+#include <TPDGCode.h>
 #include <TRandom.h>
 #include <TDatime.h>
 #include <TSystem.h>
@@ -27,7 +28,7 @@
 #include "STEER/AliConfig.h"
 #include "PYTHIA6/AliDecayerPythia.h"
 #include "PYTHIA6/AliGenPythia.h"
-#include "STEER/AliMagFMaps.h"
+#include "STEER/AliMagF.h"
 #include "STRUCT/AliBODY.h"
 #include "STRUCT/AliMAG.h"
 #include "STRUCT/AliABSOv0.h"
@@ -49,7 +50,7 @@
 #include "PMD/AliPMDv1.h"
 #include "T0/AliT0v1.h"
 #include "EMCAL/AliEMCALv2.h"
-#include "ACORDE/AliACORDEv0.h"
+#include "ACORDE/AliACORDEv1.h"
 #include "VZERO/AliVZEROv6.h"
 #endif
 
@@ -84,11 +85,7 @@ enum YCut_t
 {
   kFull, kBarrel, kMuonArm
 };
-//--- Magnetic Field ---
-enum Mag_t
-{
-    k2kG, k4kG, k5kG
-};
+
 //--- Functions ---
 AliGenPythia *PythiaHVQ(PDC06Proc_t proc);
 AliGenerator *MbCocktail();
@@ -101,7 +98,7 @@ Float_t EtaToTheta(Float_t arg){
 // This part for configuration
 static DecayHvFl_t   decHvFl   = kNature; 
 static YCut_t        ycut      = kBarrel;
-static Mag_t         mag       = k5kG; 
+static AliMagF::BMap_t         mag       = k5kG; 
 //========================//
 // Set Random Number seed //
 //========================//
@@ -250,7 +247,7 @@ void Config()
   
   // electron generator
   AliGenBox *eminus= new AliGenBox(Npart);
-  eminus->SetPart(11);
+  eminus->SetPart(kElectron);
   eminus->SetMomentumRange(momen,momen);
   eminus->SetPhiRange(0,360);
   eminus->SetThetaRange(thmin,thmax);
@@ -259,7 +256,7 @@ void Config()
 
   // positron generator
   AliGenBox *eplus= new AliGenBox(Npart);
-  eplus->SetPart(-11);
+  eplus->SetPart(kPositron);
   eplus->SetMomentumRange(momen,momen);
   eplus->SetPhiRange(0,360);
   eplus->SetThetaRange(thmin,thmax);
@@ -268,7 +265,7 @@ void Config()
 
   // mu- generator
   AliGenBox *muminus= new AliGenBox(Npart);
-  muminus->SetPart(13);
+  muminus->SetPart(kMuonMinus);
   muminus->SetMomentumRange(momen,momen);
   muminus->SetPhiRange(0,360);
   muminus->SetThetaRange(thmin,thmax);
@@ -277,7 +274,7 @@ void Config()
 
   // mu+ generator
   AliGenBox *muplus= new AliGenBox(Npart);
-  muplus->SetPart(-13);
+  muplus->SetPart(kMuonPlus);
   muplus->SetMomentumRange(momen,momen);
   muplus->SetPhiRange(0,360);
   muplus->SetThetaRange(thmin,thmax);
@@ -286,7 +283,7 @@ void Config()
 
   // pi- generator
   AliGenBox *piminus= new AliGenBox(Npart);
-  piminus->SetPart(-211);
+  piminus->SetPart(kPiMinus);
   piminus->SetMomentumRange(momen,momen);
   piminus->SetPhiRange(0,360);
   piminus->SetThetaRange(thmin,thmax);
@@ -295,7 +292,7 @@ void Config()
 
   // pi+ generator
   AliGenBox *piplus= new AliGenBox(Npart);
-  piplus->SetPart(211);
+  piplus->SetPart(kPiPlus);
   piplus->SetMomentumRange(momen,momen);
   piplus->SetPhiRange(0,360);
   piplus->SetThetaRange(thmin,thmax);
@@ -304,7 +301,7 @@ void Config()
 
   // K- generator
   AliGenBox *kminus= new AliGenBox(Npart);
-  kminus->SetPart(-321);
+  kminus->SetPart(kKMinus);
   kminus->SetMomentumRange(momen,momen);
   kminus->SetPhiRange(0,360);
   kminus->SetThetaRange(thmin,thmax);
@@ -313,7 +310,7 @@ void Config()
 
   // K+ generator
   AliGenBox *kplus= new AliGenBox(Npart);
-  kplus->SetPart(321);
+  kplus->SetPart(kKPlus);
   kplus->SetMomentumRange(momen,momen);
   kplus->SetPhiRange(0,360);
   kplus->SetThetaRange(thmin,thmax);
@@ -322,7 +319,7 @@ void Config()
 
   // p generator
   AliGenBox *proton= new AliGenBox(Npart);
-  proton->SetPart(2212);
+  proton->SetPart(kProton);
   proton->SetMomentumRange(momen,momen);
   proton->SetPhiRange(0,360);
   proton->SetThetaRange(thmin,thmax);
@@ -331,7 +328,7 @@ void Config()
 
   // anti-p generator
   AliGenBox *aproton= new AliGenBox(Npart);
-  aproton->SetPart(-2212);
+  aproton->SetPart(kProtonBar);
   aproton->SetMomentumRange(momen,momen);
   aproton->SetPhiRange(0,360);
   aproton->SetThetaRange(thmin,thmax);
@@ -355,19 +352,17 @@ void Config()
 
   // FIELD
   //    
-  if (mag == k2kG) {
+  if (mag == AliMagF::k2kG) {
     comment = comment.Append(" | L3 field 0.2 T");
-  } else if (mag == k4kG) {
-    comment = comment.Append(" | L3 field 0.4 T");
-  } else if (mag == k5kG) {
+  } else if (mag == AliMagF::k5kG) {
     comment = comment.Append(" | L3 field 0.5 T");
   }
   printf("\n \n Comment: %s \n \n", comment.Data());
     
-  AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., mag);
+  AliMagF* field = new AliMagF("Maps","Maps", 1., 1., mag);
   field->SetL3ConstField(0); //Using const. field in the barrel
+  TGeoGlobalMagField::Instance()->SetField(field);
   rl->CdGAFile();
-  gAlice->SetField(field);    
 
   Int_t iABSO  = 1;
   Int_t iACORDE   = 0;
@@ -583,7 +578,7 @@ void Config()
      if (iACORDE)
     {
         //=================== ACORDE parameters ============================
-        AliACORDE *ACORDE = new AliACORDEv0("ACORDE", "normal ACORDE");
+        AliACORDE *ACORDE = new AliACORDEv1("ACORDE", "normal ACORDE");
     }
 
      if (iVZERO)