Overlaps corrected, new shape of sectors
[u/mrichter/AliRoot.git] / FMD / Config.C
index 0ede4e6..54961d6 100644 (file)
     @date    Mon Mar 27 12:50:29 2006
     @brief   Simulation configuration script
 */
-#if !defined(__CINT__) || defined(__MAKECINT__)
-#include <Riostream.h>
-#include <TPDGCode.h>
-#include <TSystem.h>
-#include <TVirtualMC.h>
-#include <TGeant3.h>
-#include "STEER/AliRunLoader.h"
-#include "STEER/AliRun.h"
-#include "STEER/AliConfig.h"
-#include "STEER/AliGenerator.h"
-#include "PYTHIA6/AliDecayerPythia.h"
-#include "EVGEN/AliGenHIJINGpara.h"
-#include "THijing/AliGenHijing.h"
-#include "EVGEN/AliGenCocktail.h"
-#include "EVGEN/AliGenSlowNucleons.h"
-#include "EVGEN/AliSlowNucleonModelExp.h"
-#include "EVGEN/AliGenParam.h"
-#include "EVGEN/AliGenMUONlib.h"
-#include "EVGEN/AliGenMUONCocktail.h"
-#include "PYTHIA6/AliGenPythia.h"
-#include "STEER/AliMagFMaps.h"
-#include "STRUCT/AliBODY.h"
-#include "STRUCT/AliMAG.h"
-#include "STRUCT/AliABSOv0.h"
-#include "STRUCT/AliDIPOv2.h"
-#include "STRUCT/AliHALL.h"
-#include "STRUCT/AliFRAMEv2.h"
-#include "STRUCT/AliSHILv2.h"
-#include "STRUCT/AliPIPEv3.h"
-#include "ITS/AliITSvPPRasymmFMD.h"
-#include "TPC/AliTPCv2.h"
-#include "TOF/AliTOFv4T0.h"
-#include "HMPID/AliHMPIDv1.h"
-#include "ZDC/AliZDCv2.h"
-#include "TRD/AliTRDv1.h"
-#include "FMD/AliFMDv1.h"
-#include "MUON/AliMUONv1.h"
-#include "MUON/AliMUONSt1GeometryBuilderV2.h"
-#include "MUON/AliMUONSt2GeometryBuilder.h"
-#include "MUON/AliMUONSlatGeometryBuilder.h"
-#include "MUON/AliMUONTriggerGeometryBuilder.h"
-#include "PHOS/AliPHOSv1.h"
-#include "PMD/AliPMDv1.h"
-#include "T0/AliT0v1.h"
-#include "EMCAL/AliEMCALv2.h"
-#include "ACORDE/AliACORDEv1.h"
-#include "VZERO/AliVZEROv2.h"
-#endif
 
 //____________________________________________________________________
 // 
@@ -200,13 +152,6 @@ enum Rad_t {
 };
 
 //____________________________________________________________________
-enum Mag_t {
-  k2kG,                                //
-  k4kG,                                //
-  k5kG                         //
-};
-
-//____________________________________________________________________
 enum MC_t {
   kFLUKA, 
   kGEANT3, 
@@ -217,8 +162,8 @@ enum MC_t {
 //____________________________________________________________________
 // Functions
 Float_t       EtaToTheta(Float_t eta);
-Eg_t          LookupEG(const Char_t* name);
-AliGenerator* GeneratorFactory(EG_t eg, Rad_t rad, TString& comment);
+EG_t          LookupEG(const Char_t* name);
+AliGenerator* GeneratorFactory(EG_t eg, Rad_t rad);
 AliGenHijing* HijingStandard();
 void          ProcessEnvironmentVars(EG_t& eg, Int_t& seed);
 
@@ -228,22 +173,14 @@ Config()
 {
   //____________________________________________________________________
   // This part for configuration    
-  // EG_t  eg   = test50;
-  // EG_t  eg   = kParam_fmd;
-  // EG_t  eg   = kFMDFlat; // kParam_2000; // kPythia;
-  // EG_t  eg   = kFMDFlat;
   EG_t  eg   = kPythia6;
   Geo_t geo  = kNoHoles;
   Rad_t rad  = kGluonRadiation;
-  Mag_t mag  = k5kG;
+  AliMagF::BMap_t mag  = AliMagF::k5kG;
   Int_t seed = 12345; //Set 0 to use the current time
   MC_t  mc   = kGEANT3TGEO;
   
   //____________________________________________________________________
-  // Comment line 
-  static TString  comment;
-  
-  //____________________________________________________________________
   // Get settings from environment variables
   ProcessEnvironmentVars(eg, seed);
 
@@ -274,7 +211,7 @@ Config()
       //
       // Libraries needed by GEANT 3.21 
       //
-      gSystem->Load("$ALICE_ROOT/lib/tgt_$ALICE_TARGET/libpythia6.so");
+      gSystem->Load("$ALICE_ROOT/lib/tgt_$ALICE_TARGET/libpythia6");
       gSystem->Load("libgeant321");
       
       // 
@@ -289,10 +226,10 @@ Config()
       //
       // Libraries needed by GEANT 3.21 
       //
-      gSystem->Load("$ALICE_ROOT/lib/tgt_$ALICE_TARGET/liblhapdf.so");
-      gSystem->Load("$ALICE_ROOT/lib/tgt_$ALICE_TARGET/libpythia6.so");
-      // gSystem->Load("libEGPythia6.so"); //<- For non-debian (sigh!)
-      gSystem->Load("EGPythia6.so");
+      gSystem->Load("$ALICE_ROOT/lib/tgt_$ALICE_TARGET/liblhapdf");
+      gSystem->Load("$ALICE_ROOT/lib/tgt_$ALICE_TARGET/libpythia6");
+      gSystem->Load("libEGPythia6"); //<- For non-debian (sigh!)
+      // gSystem->Load("EGPythia6");
       gSystem->Load("libgeant321");
     
       // 
@@ -392,7 +329,7 @@ Config()
   gMC->SetProcess("MUNU",1);
   gMC->SetProcess("CKOV",1);
   gMC->SetProcess("HADR",1);
-  gMC->SetProcess("LOSS",2); // 0:none 1,3:dray 2:nodray 4:nofluct (def:2)
+  gMC->SetProcess("LOSS",1); // 0:none 1,3:dray 2:nodray 4:nofluct (def:2)
   gMC->SetProcess("MULS",1);
   gMC->SetProcess("RAYL",1);
 
@@ -414,7 +351,7 @@ Config()
   
   //__________________________________________________________________
   // Generator Configuration
-  AliGenerator* gener = GeneratorFactory(eg, rad, comment);
+  AliGenerator* gener = GeneratorFactory(eg, rad);
   gener->SetOrigin(0, 0, 0);    // vertex position
   gener->SetSigma(0, 0, 5.3);   // Sigma in (X,Y,Z) (cm) on IP position
   gener->SetCutVertexZ(1.);     // Truncate at 1 sigma
@@ -423,35 +360,12 @@ Config()
   gener->Init();
     
   //__________________________________________________________________
-  // 
-  // Comments 
-  // 
-  switch (mag) {
-  case k2kG: comment = comment.Append(" | L3 field 0.2 T"); break;
-  case k4kG: comment = comment.Append(" | L3 field 0.4 T"); break;
-  case k5kG: comment = comment.Append(" | L3 field 0.5 T"); break;
-  }
-
-  switch (rad) {
-  case kGluonRadiation: 
-    comment = comment.Append(" | Gluon Radiation On");  break;
-  default:
-    comment = comment.Append(" | Gluon Radiation Off"); break;
-  }
-
-  switch(geo) {
-  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;
-
-  //__________________________________________________________________
   // Field (L3 0.4 T)
-  AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., mag);
-  field->SetL3ConstField(0); //Using const. field in the barrel
+  AliMagF* field = new AliMagF("Maps","Maps",-1., -1., mag);
+  TGeoGlobalMagField::Instance()->SetField(field);
+
   rl->CdGAFile();
-  gAlice->SetField(field);    
+
   TFile* magF = TFile::Open("mag.root", "RECREATE");
   field->Write("mag");
   magF->Close();
@@ -474,7 +388,7 @@ Config()
   Bool_t usePMD   = kFALSE; 
   Bool_t useHMPID = kFALSE; 
   Bool_t useSHIL  = kFALSE; 
-  Bool_t useT0    = kFALSE; 
+  Bool_t useT0    = kTRUE; 
   Bool_t useTOF   = kFALSE; 
   Bool_t useTPC   = kFALSE;
   Bool_t useTRD   = kFALSE; 
@@ -500,7 +414,7 @@ Config()
   }
   if (useSHIL)   AliSHIL   *SHIL   = new AliSHILv3("SHIL", "Shielding v3");
   if (usePIPE)   AliPIPE   *PIPE   = new AliPIPEv3("PIPE", "Beam Pipe");
-  if (useITS)    AliITS    *ITS    = new AliITSv11Hybrid("ITS","ITS v11Hybrid");
+  if (useITS)    AliITS    *ITS    = new AliITSv11("ITS","ITS v11");
   if (useTPC)    AliTPC    *TPC    = new AliTPCv2("TPC", "Default");
   if (useTOF)    AliTOF    *TOF    = new AliTOFv6T0("TOF", "normal TOF");
   if (useHMPID)  AliHMPID  *HMPID  = new AliHMPIDv1("HMPID", "normal HMPID");
@@ -511,7 +425,7 @@ Config()
   if (usePHOS)   AliPHOS   *PHOS   = new AliPHOSv1("PHOS", "IHEP");
   if (usePMD)    AliPMD    *PMD    = new AliPMDv1("PMD", "normal PMD");
   if (useT0)     AliT0     *T0     = new AliT0v1("T0", "T0 Detector");
-  if (useEMCAL)  AliEMCAL  *EMCAL  = new AliEMCALv2("EMCAL", "EMCAL_COMPLETE");
+  if (useEMCAL)  AliEMCAL  *EMCAL  = new AliEMCALv2("EMCAL", "EMCAL_COMPLETEV1");
   if (useACORDE) AliACORDE *ACORDE = new AliACORDEv1("ACORDE", "normal ACORDE");
   if (useVZERO)  AliVZERO  *VZERO  = new AliVZEROv7("VZERO", "normal VZERO");
 }
@@ -535,7 +449,7 @@ LookupEG(const Char_t* name)
 
 //____________________________________________________________________  
 AliGenerator* 
-GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)  
+GeneratorFactory(EG_t eg, Rad_t rad)  
 {
   Int_t isw = 3;
   if (rad == kNoGluonRadiation) isw = 0;
@@ -555,10 +469,10 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
   case kPythia6Jets104_125:
   case kPythia6Jets125_150:
   case kPythia6Jets150_180:
-    gSystem->Load("liblhapdf.so");
-    // gSystem->Load("/usr/lib/libpythia.so");
+    gSystem->Load("liblhapdf");
+    // gSystem->Load("/usr/lib/libpythia");
     // gSystem->ListLibraries();
-    gSystem->Load("EGPythia6.so");
+    gSystem->Load("EGPythia6");
     gSystem->Load("libAliPythia6");
     break;
   }
@@ -567,7 +481,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
   switch (eg) {
   case test50:
     {
-      comment = comment.Append(":HIJINGparam test 50 particles");
       AliGenHIJINGpara *gener = new AliGenHIJINGpara(50);
       gener->SetMomentumRange(0, 999999.);
       gener->SetPhiRange(0., 360.);
@@ -580,7 +493,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kParam_8000:
     {
-      comment = comment.Append(":HIJINGparam N=8000");
       AliGenHIJINGpara *gener = new AliGenHIJINGpara(86030);
       gener->SetMomentumRange(0, 999999.);
       gener->SetPhiRange(0., 360.);
@@ -593,7 +505,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kParam_4000:
     {
-      comment = comment.Append("HIJINGparam N=4000");
       AliGenHIJINGpara *gener = new AliGenHIJINGpara(43015);
       gener->SetMomentumRange(0, 999999.);
       gener->SetPhiRange(0., 360.);
@@ -606,7 +517,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kParam_2000:
     {
-      comment = comment.Append("HIJINGparam N=2000");
       AliGenHIJINGpara *gener = new AliGenHIJINGpara(21507);
       gener->SetMomentumRange(0, 999999.);
       gener->SetPhiRange(0., 360.);
@@ -619,7 +529,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kParam_fmd:
     {
-      comment = comment.Append("HIJINGparam N=100");
       AliGenHIJINGpara *gener = new AliGenHIJINGpara(500);
       gener->SetMomentumRange(0, 999999.);
       gener->SetPhiRange(0., 360.);
@@ -635,7 +544,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     //
   case kHijing_cent1:
     {
-      comment = comment.Append("HIJING cent1");
       AliGenHijing *gener = HijingStandard();
       // impact parameter range
       gener->SetImpactParameterRange(0., 5.);
@@ -644,7 +552,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kHijing_cent2:
     {
-      comment = comment.Append("HIJING cent2");
       AliGenHijing *gener = HijingStandard();
       // impact parameter range
       gener->SetImpactParameterRange(0., 2.);
@@ -656,7 +563,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     //
   case kHijing_per1:
     {
-      comment = comment.Append("HIJING per1");
       AliGenHijing *gener = HijingStandard();
       // impact parameter range
       gener->SetImpactParameterRange(5., 8.6);
@@ -665,7 +571,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kHijing_per2:
     {
-      comment = comment.Append("HIJING per2");
       AliGenHijing *gener = HijingStandard();
       // impact parameter range
       gener->SetImpactParameterRange(8.6, 11.2);
@@ -674,7 +579,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kHijing_per3:
     {
-      comment = comment.Append("HIJING per3");
       AliGenHijing *gener = HijingStandard();
       // impact parameter range
       gener->SetImpactParameterRange(11.2, 13.2);
@@ -683,7 +587,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kHijing_per4:
     {
-      comment = comment.Append("HIJING per4");
       AliGenHijing *gener = HijingStandard();
       // impact parameter range
       gener->SetImpactParameterRange(13.2, 15.);
@@ -692,7 +595,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kHijing_per5:
     {
-      comment = comment.Append("HIJING per5");
       AliGenHijing *gener = HijingStandard();
       // impact parameter range
       gener->SetImpactParameterRange(15., 100.);
@@ -704,7 +606,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     //
   case kHijing_jj25:
     {
-      comment = comment.Append("HIJING Jet 25 GeV");
       AliGenHijing *gener = HijingStandard();
       // impact parameter range
       gener->SetImpactParameterRange(0., 5.);
@@ -721,7 +622,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
 
   case kHijing_jj50:
     {
-      comment = comment.Append("HIJING Jet 50 GeV");
       AliGenHijing *gener = HijingStandard();
       // impact parameter range
       gener->SetImpactParameterRange(0., 5.);
@@ -738,7 +638,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
 
   case kHijing_jj75:
     {
-      comment = comment.Append("HIJING Jet 75 GeV");
       AliGenHijing *gener = HijingStandard();
       // impact parameter range
       gener->SetImpactParameterRange(0., 5.);
@@ -755,7 +654,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
 
   case kHijing_jj100:
     {
-      comment = comment.Append("HIJING Jet 100 GeV");
       AliGenHijing *gener = HijingStandard();
       // impact parameter range
       gener->SetImpactParameterRange(0., 5.);
@@ -772,7 +670,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
 
   case kHijing_jj200:
     {
-      comment = comment.Append("HIJING Jet 200 GeV");
       AliGenHijing *gener = HijingStandard();
       // impact parameter range
       gener->SetImpactParameterRange(0., 5.);
@@ -791,7 +688,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     //
   case kHijing_gj25:
     {
-      comment = comment.Append("HIJING Gamma 25 GeV");
       AliGenHijing *gener = HijingStandard();
       // impact parameter range
       gener->SetImpactParameterRange(0., 5.);
@@ -808,7 +704,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
 
   case kHijing_gj50:
     {
-      comment = comment.Append("HIJING Gamma 50 GeV");
       AliGenHijing *gener = HijingStandard();
       // impact parameter range
       gener->SetImpactParameterRange(0., 5.);
@@ -825,7 +720,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
 
   case kHijing_gj75:
     {
-      comment = comment.Append("HIJING Gamma 75 GeV");
       AliGenHijing *gener = HijingStandard();
       // impact parameter range
       gener->SetImpactParameterRange(0., 5.);
@@ -842,7 +736,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
 
   case kHijing_gj100:
     {
-      comment = comment.Append("HIJING Gamma 100 GeV");
       AliGenHijing *gener = HijingStandard();
       // impact parameter range
       gener->SetImpactParameterRange(0., 5.);
@@ -859,7 +752,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
 
   case kHijing_gj200:
     {
-      comment = comment.Append("HIJING Gamma 200 GeV");
       AliGenHijing *gener = HijingStandard();
       // impact parameter range
       gener->SetImpactParameterRange(0., 5.);
@@ -875,7 +767,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kHijing_pA:
     {
-      comment = comment.Append("HIJING pA");
 
       AliGenCocktail *gener  = new AliGenCocktail();
 
@@ -912,7 +803,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kPythia6:
     {
-      comment = comment.Append(":Pythia p-p @ 14 TeV");
       AliGenPythia *gener = new AliGenPythia(-1); 
       gener->SetMomentumRange(0,999999);
       gener->SetThetaRange(0., 180.);
@@ -925,7 +815,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kPythia6Jets20_24:
     {
-      comment = comment.Append(":Pythia jets 20-24 GeV @ 5.5 TeV");
       AliGenPythia * gener = new AliGenPythia(-1);
       gener->SetEnergyCMS(5500.);//        Centre of mass energy
       gener->SetProcess(kPyJets);//        Process type
@@ -944,7 +833,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kPythia6Jets24_29:
     {
-      comment = comment.Append(":Pythia jets 24-29 GeV @ 5.5 TeV");
       AliGenPythia * gener = new AliGenPythia(-1);
       gener->SetEnergyCMS(5500.);//        Centre of mass energy
       gener->SetProcess(kPyJets);//        Process type
@@ -963,7 +851,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kPythia6Jets29_35:
     {
-      comment = comment.Append(":Pythia jets 29-35 GeV @ 5.5 TeV");
       AliGenPythia * gener = new AliGenPythia(-1);
       gener->SetEnergyCMS(5500.);//        Centre of mass energy
       gener->SetProcess(kPyJets);//        Process type
@@ -982,7 +869,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kPythia6Jets35_42:
     {
-      comment = comment.Append(":Pythia jets 35-42 GeV @ 5.5 TeV");
       AliGenPythia * gener = new AliGenPythia(-1);
       gener->SetEnergyCMS(5500.);//        Centre of mass energy
       gener->SetProcess(kPyJets);//        Process type
@@ -1001,7 +887,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kPythia6Jets42_50:
     {
-      comment = comment.Append(":Pythia jets 42-50 GeV @ 5.5 TeV");
       AliGenPythia * gener = new AliGenPythia(-1);
       gener->SetEnergyCMS(5500.);//        Centre of mass energy
       gener->SetProcess(kPyJets);//        Process type
@@ -1020,7 +905,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kPythia6Jets50_60:
     {
-      comment = comment.Append(":Pythia jets 50-60 GeV @ 5.5 TeV");
       AliGenPythia * gener = new AliGenPythia(-1);
       gener->SetEnergyCMS(5500.);//        Centre of mass energy
       gener->SetProcess(kPyJets);//        Process type
@@ -1039,7 +923,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kPythia6Jets60_72:
     {
-      comment = comment.Append(":Pythia jets 60-72 GeV @ 5.5 TeV");
       AliGenPythia * gener = new AliGenPythia(-1);
       gener->SetEnergyCMS(5500.);//        Centre of mass energy
       gener->SetProcess(kPyJets);//        Process type
@@ -1058,7 +941,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kPythia6Jets72_86:
     {
-      comment = comment.Append(":Pythia jets 72-86 GeV @ 5.5 TeV");
       AliGenPythia * gener = new AliGenPythia(-1);
       gener->SetEnergyCMS(5500.);//        Centre of mass energy
       gener->SetProcess(kPyJets);//        Process type
@@ -1077,7 +959,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kPythia6Jets86_104:
     {
-      comment = comment.Append(":Pythia jets 86-104 GeV @ 5.5 TeV");
       AliGenPythia * gener = new AliGenPythia(-1);
       gener->SetEnergyCMS(5500.);//        Centre of mass energy
       gener->SetProcess(kPyJets);//        Process type
@@ -1096,7 +977,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kPythia6Jets104_125:
     {
-      comment = comment.Append(":Pythia jets 105-125 GeV @ 5.5 TeV");
       AliGenPythia * gener = new AliGenPythia(-1);
       gener->SetEnergyCMS(5500.);//        Centre of mass energy
       gener->SetProcess(kPyJets);//        Process type
@@ -1115,7 +995,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kPythia6Jets125_150:
     {
-      comment = comment.Append(":Pythia jets 125-150 GeV @ 5.5 TeV");
       AliGenPythia * gener = new AliGenPythia(-1);
       gener->SetEnergyCMS(5500.);//        Centre of mass energy
       gener->SetProcess(kPyJets);//        Process type
@@ -1134,7 +1013,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kPythia6Jets150_180:
     {
-      comment = comment.Append(":Pythia jets 150-180 GeV @ 5.5 TeV");
       AliGenPythia * gener = new AliGenPythia(-1);
       gener->SetEnergyCMS(5500.);//        Centre of mass energy
       gener->SetProcess(kPyJets);//        Process type
@@ -1153,7 +1031,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kD0PbPb5500:
     {
-      comment = comment.Append(" D0 in Pb-Pb at 5.5 TeV");
       AliGenPythia * gener = new AliGenPythia(10);
       gener->SetProcess(kPyD0PbPbMNR);
       gener->SetStrucFunc(kCTEQ4L);
@@ -1170,7 +1047,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kCharmSemiElPbPb5500:
     {
-      comment = comment.Append(" Charm in Pb-Pb at 5.5 TeV");
       AliGenPythia * gener = new AliGenPythia(10);
       gener->SetProcess(kPyCharmPbPbMNR);
       gener->SetStrucFunc(kCTEQ4L);
@@ -1186,7 +1062,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kBeautySemiElPbPb5500:
     {
-      comment = comment.Append(" Beauty in Pb-Pb at 5.5 TeV");
       AliGenPythia *gener = new AliGenPythia(10);
       gener->SetProcess(kPyBeautyPbPbMNR);
       gener->SetStrucFunc(kCTEQ4L);
@@ -1202,7 +1077,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kCocktailTRD:
     {
-      comment = comment.Append(" Cocktail for TRD at 5.5 TeV");
       AliGenCocktail *gener  = new AliGenCocktail();
 
       AliGenParam *jpsi = new AliGenParam(10,
@@ -1248,7 +1122,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kPyJJ:
     {
-      comment = comment.Append(" Jet-jet at 5.5 TeV");
       AliGenPythia *gener = new AliGenPythia(-1);
       gener->SetEnergyCMS(5500.);
       gener->SetProcess(kPyJets);
@@ -1262,7 +1135,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kPyGJ:
     {
-      comment = comment.Append(" Gamma-jet at 5.5 TeV");
       AliGenPythia *gener = new AliGenPythia(-1);
       gener->SetEnergyCMS(5500.);
       gener->SetProcess(kPyDirectGamma);
@@ -1277,7 +1149,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kMuonCocktailCent1:
     {
-      comment = comment.Append(" Muon Cocktail Cent1");
       AliGenMUONCocktail * gener = new AliGenMUONCocktail();
       gener->SetPtRange(1.0,100.);       // Transverse momentum range   
       gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
@@ -1292,7 +1163,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kMuonCocktailPer1:
     {
-      comment = comment.Append(" Muon Cocktail Per1");
       AliGenMUONCocktail * gener = new AliGenMUONCocktail();
       gener->SetPtRange(1.0,100.);       // Transverse momentum range   
       gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
@@ -1307,7 +1177,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kMuonCocktailPer4:
     {
-      comment = comment.Append(" Muon Cocktail Per4");
       AliGenMUONCocktail * gener = new AliGenMUONCocktail();
       gener->SetPtRange(1.0,100.);       // Transverse momentum range   
       gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
@@ -1322,7 +1191,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kMuonCocktailCent1HighPt:
     {
-      comment = comment.Append(" Muon Cocktail HighPt Cent1");
       AliGenMUONCocktail * gener = new AliGenMUONCocktail();
       gener->SetPtRange(1.0,100.);       // Transverse momentum range   
       gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
@@ -1337,7 +1205,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kMuonCocktailPer1HighPt :
     {
-      comment = comment.Append(" Muon Cocktail HighPt Per1");
       AliGenMUONCocktail * gener = new AliGenMUONCocktail();
       gener->SetPtRange(1.0,100.);       // Transverse momentum range   
       gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
@@ -1352,7 +1219,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kMuonCocktailPer4HighPt:
     {
-      comment = comment.Append(" Muon Cocktail HighPt Per4");
       AliGenMUONCocktail * gener = new AliGenMUONCocktail();
       gener->SetPtRange(1.0,100.);       // Transverse momentum range   
       gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
@@ -1367,7 +1233,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kMuonCocktailCent1Single:
     {
-      comment = comment.Append(" Muon Cocktail Single Cent1");
       AliGenMUONCocktail * gener = new AliGenMUONCocktail();
       gener->SetPtRange(1.0,100.);       // Transverse momentum range   
       gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
@@ -1382,7 +1247,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kMuonCocktailPer1Single :
     {
-      comment = comment.Append(" Muon Cocktail Single Per1");
       AliGenMUONCocktail * gener = new AliGenMUONCocktail();
       gener->SetPtRange(1.0,100.);       // Transverse momentum range   
       gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
@@ -1397,7 +1261,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kMuonCocktailPer4Single:
     {
-      comment = comment.Append(" Muon Cocktail Single Per4");
       AliGenMUONCocktail * gener = new AliGenMUONCocktail();
       gener->SetPtRange(1.0,100.);       // Transverse momentum range   
       gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
@@ -1412,7 +1275,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kFMD1Flat: 
     {
-      comment = comment.Append(" Flat in FMD1 range");
       AliGenBox* gener = new AliGenBox(2000);
       gener->SetPart(kPiPlus);
       gener->SetMomentumRange(3,4);
@@ -1423,8 +1285,7 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kFMD2Flat: 
     {
-      comment = comment.Append(" Flat in FMD2 range");
-      AliGenBox* gener = new AliGenBox(2000);
+      AliGenBox* gener = new AliGenBox(100);
       gener->SetPart(kPiPlus);
       gener->SetMomentumRange(3,4);
       gener->SetPhiRange(0, 360);
@@ -1434,7 +1295,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kFMD3Flat: 
     {
-      comment = comment.Append(" Flat in FMD3 range");
       AliGenBox* gener = new AliGenBox(2000);
       gener->SetPart(kPiPlus);
       gener->SetMomentumRange(3,4);
@@ -1445,7 +1305,6 @@ GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)
     break;
   case kFMDFlat:
     {
-      comment = comment.Append(" Flat in FMD range");
       AliGenCocktail* gener = new AliGenCocktail();
       gener->SetMomentumRange(3,4);
       gener->SetPhiRange(0, 360);