Extra header added to the list

[u/mrichter/AliRoot.git] / EMCAL / macros / Config.C

diff --git a/EMCAL/macros/Config.C b/EMCAL/macros/Config.C
old mode 100755 (executable)
new mode 100644 (file)
index 015f5f1..908b77c
--- a/EMCAL/macros/Config.C
+++ b/EMCAL/macros/Config.C
@@ -1,4 +1,3 @@
-// Config modified from production LHC12a11a
 // One can use the configuration macro in compiled mode by
 // root [0] gSystem->Load("libgeant321");
 // root [0] gSystem->SetIncludePath("-I$ROOTSYS/include -I$ALICE_ROOT/include\
@@ -135,38 +134,32 @@ void Config()
   gMC->SetCut("TOFMAX", tofmax); 
   
   
-//  int     nParticles = 1;
-//  if (gSystem->Getenv("CONFIG_NPARTICLES"))
-//  {
-//    nParticles = atoi(gSystem->Getenv("CONFIG_NPARTICLES"));
-//  }
-//  
-//  if (gSystem->Getenv("CONFIG_YEAR"))
-//  {
-//    year = atoi(gSystem->Getenv("CONFIG_YEAR"));
-//  }
-//  
-//  AliGenBox *gener = new AliGenBox(nParticles);
-//  gener->SetMomentumRange(1.,10.);
-//  
-//  if     (year == 2010)
-//    gener->SetPhiRange(80.0,120.0);
-//  else if(year == 2011)
-//    gener->SetPhiRange(80.0,180.0);
-//  else
-//    gener->SetPhiRange(80.0,190.0);
-//  
-//  gener->SetThetaRange(EtaToTheta(0.7), EtaToTheta(-0.7));
-//  
-//  gener->SetOrigin(0,0,0);        //vertex position
-//  gener->SetSigma(0,0,0);         //Sigma in (X,Y,Z) (cm) on IP position
-//  gener->SetPart(221);
-//  gener->SetDecayParticle(22);
-//  gener->SetNeutralMesonDecayInto2Photon(kTRUE);
-//  gener->SetDecayPhiRange(80,180);
-//  gener->SetDecayEtaRange(-0.7,0.7);
-//  gener->SetDecayERange(0.5,10000);
-//  gener->SetDecayPtRange(0.5,10000);
+  int     nParticles = 1;
+  if (gSystem->Getenv("CONFIG_NPARTICLES"))
+  {
+    nParticles = atoi(gSystem->Getenv("CONFIG_NPARTICLES"));
+  }
+  
+  if (gSystem->Getenv("CONFIG_YEAR"))
+  {
+    year = atoi(gSystem->Getenv("CONFIG_YEAR"));
+  }
+  
+  AliGenBox *gener = new AliGenBox(nParticles);
+  gener->SetMomentumRange(10.,10.);
+  
+  gener->SetPart(22);
+  
+  if     (year == 2010)
+    gener->SetPhiRange(80.0,120.0);
+  else if(year == 2011)
+    gener->SetPhiRange(80.0,180.0);
+  else if(year == 2012 || year == 2013)
+    gener->SetPhiRange(80.0,190.0);
+  else
+    gener->SetPhiRange(80.0,330.0); // Include DCal
+  
+  gener->SetThetaRange(EtaToTheta(0.7), EtaToTheta(-0.7));
   
 //  AliGenLib* lib   = new AliGenPHOSlib();
 //  Int_t      type  = AliGenPHOSlib::kEtaFlat;
@@ -184,28 +177,21 @@ void Config()
 //     gener->SetSigma(0,0,5.3);       //Sigma in (X,Y,Z) (cm) on IP position
 //     gener->SetForceDecay(kGammaEM);
 //
-//  gener->SetTrackingFlag(0);
-  
-  AliGenCocktail *gener = new AliGenCocktail();
-  gener->SetProjectile("A", 208, 82);
-  gener->SetTarget    ("A", 208, 82);
-  
-  // 1 Pi0 in EMCAL, 2010 configuration, 4 SM
-  AliGenParam *gEMCPi0 = GenParamCalo(1, AliGenPHOSlib::kPi0Flat, "EMCAL");
-  gener->AddGenerator(gEMCPi0,"pi0EMC", 1);
-  
-  // 1 Pi0 in PHOS
-  AliGenParam *gPHSPi0 = GenParamCalo(1, AliGenPHOSlib::kPi0Flat, "PHOS");
-  gener->AddGenerator(gPHSPi0,"pi0PHS", 1);
-  
-  // 1 Eta in EMCAL, 2010 configuration, 4 SM
-  AliGenParam *gEMCEta = GenParamCalo(1, AliGenPHOSlib::kEtaFlat, "EMCAL");
-  gener->AddGenerator(gEMCEta,"etaEMC", 1);
-  
-  // 1 Pi0 in PHOS
-  AliGenParam *gPHSEta = GenParamCalo(1, AliGenPHOSlib::kEtaFlat, "PHOS");
-  gener->AddGenerator(gPHSEta,"etaPHS", 1);
+//  //gener->SetTrackingFlag(0);
   
+//  AliGenCocktail *gener = new AliGenCocktail();
+//  gener->SetProjectile("A", 208, 82);
+//  gener->SetTarget    ("A", 208, 82);
+//  
+//  // 1 Pi0 in EMCAL, 2010 configuration, 4 SM
+//  AliGenParam *gEMCPi0 = GenParamCalo(1, AliGenPHOSlib::kPi0Flat, "EMCAL");
+//  gener->AddGenerator(gEMCPi0,"pi0EMC", 1);
+//  
+//  // 1 Eta in EMCAL, 2010 configuration, 4 SM
+//  AliGenParam *gEMCEta = GenParamCalo(1, AliGenPHOSlib::kEtaFlat, "EMCAL");
+//  gener->AddGenerator(gEMCEta,"etaEMC", 1);
+//  
+ 
   gener->Init();
   
   // 
@@ -224,7 +210,7 @@ void Config()
   Int_t   iITS   =  0;
   Int_t   iMAG   =  0;
   Int_t   iMUON  =  0;
-  Int_t   iPHOS  =  1;
+  Int_t   iPHOS  =  0;
   Int_t   iPIPE  =  0;
   Int_t   iPMD   =  0;
   Int_t   iHMPID =  0;
@@ -394,14 +380,15 @@ void Config()
   if (iEMCAL)
   {
     //=================== EMCAL parameters ============================
+    AliEMCAL *EMCAL = 0;
     if      (year == 2010)  // d phi = 40 degrees
-      AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_FIRSTYEARV1",    checkGeoAndRun);
+      EMCAL = new AliEMCALv2("EMCAL", "EMCAL_FIRSTYEARV1",    checkGeoAndRun);
     else if (year == 2011)  // d phi = 100 degrees
-      AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETEV1",     checkGeoAndRun);
-    else if (year > 2011)   // d phi = 110 degrees
-      AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETE12SMV1", checkGeoAndRun);
-    else // Old configuration with 110 degrees but not perfect geometry
-      AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETE",       checkGeoAndRun);
+      EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETEV1",     checkGeoAndRun);
+    else if (year == 2012 || year == 2013)   // d phi = 107 degrees
+      EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETE12SMV1", checkGeoAndRun);
+    else
+      EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETE12SMV1_DCAL_8SM", checkGeoAndRun); // EMCAL+DCAL dphi = 107 (EMCAL) + 33 (gap) + 67 (DCAL)
   }
   
   if (iACORDE)
@@ -427,10 +414,10 @@ Float_t EtaToTheta(Float_t 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
+  gSystem->Load("liblhapdf");      // Parton density functions                                 
+  gSystem->Load("libEGPythia6");   // TGenerator interface                                     
+  gSystem->Load("libpythia6");     // Pythia                                                   
+  gSystem->Load("libAliPythia6");  // ALICE specific
   // implementations                           
 }