]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - FMD/Config.C
ACORDEv1 is now the default
[u/mrichter/AliRoot.git] / FMD / Config.C
index 39fb3d85cd81968c22e94d867f2501b007634b78..88efeace803cb351277343a0a3ed4469771ecddc 100644 (file)
@@ -8,9 +8,14 @@
 //                   -I$ALICE_ROOT -I$ALICE/geant3/TGeant3");
 // root [0] .x grun.C(1,"ConfigPPR.C++")
 //
+/** @file    Config.C
+    @author  Christian Holm Christensen <cholm@nbi.dk>
+    @date    Mon Mar 27 12:50:29 2006
+    @brief   Simulation configuration script
+*/
 #if !defined(__CINT__) || defined(__MAKECINT__)
 #include <Riostream.h>
-#include <TRandom.h>
+#include <TPDGCode.h>
 #include <TSystem.h>
 #include <TVirtualMC.h>
 #include <TGeant3.h>
@@ -40,7 +45,7 @@
 #include "ITS/AliITSvPPRasymmFMD.h"
 #include "TPC/AliTPCv2.h"
 #include "TOF/AliTOFv4T0.h"
-#include "RICH/AliRICHv1.h"
+#include "HMPID/AliHMPIDv1.h"
 #include "ZDC/AliZDCv2.h"
 #include "TRD/AliTRDv1.h"
 #include "FMD/AliFMDv1.h"
@@ -51,9 +56,9 @@
 #include "MUON/AliMUONTriggerGeometryBuilder.h"
 #include "PHOS/AliPHOSv1.h"
 #include "PMD/AliPMDv1.h"
-#include "START/AliSTARTv1.h"
-#include "EMCAL/AliEMCALv1.h"
-#include "CRT/AliCRTv0.h"
+#include "T0/AliT0v1.h"
+#include "EMCAL/AliEMCALv2.h"
+#include "ACORDE/AliACORDEv1.h"
 #include "VZERO/AliVZEROv2.h"
 #endif
 
 // 
 // Generator types 
 //
-enum PprRun_t {
+enum EG_t {
   test50,
   kParam_8000,                 //
   kParam_4000,                 //
   kParam_2000,                 //
+  kParam_fmd,                  //
   kHijing_cent1,               //
   kHijing_cent2,               //
   kHijing_per1,                        //
@@ -112,18 +118,23 @@ enum PprRun_t {
   kMuonCocktailCent1Single,    //
   kMuonCocktailPer1Single,     //
   kMuonCocktailPer4Single,
-  kRunMax
+  kFMD1Flat, 
+  kFMD2Flat, 
+  kFMD3Flat,
+  kFMDFlat,
+  kEgMax
 };
 
 //____________________________________________________________________
 // 
 // Generator types names
 //
-const char* pprRunName[kRunMax] = {
+const char* egName[kEgMax] = {
   "test50",
   "kParam_8000",               //
   "kParam_4000",               //
   "kParam_2000",               //
+  "kParam_fmd",                        //
   "kHijing_cent1",             //
   "kHijing_cent2",             //
   "kHijing_per1",              //
@@ -169,23 +180,27 @@ const char* pprRunName[kRunMax] = {
   "kMuonCocktailPer4HighPt",   //
   "kMuonCocktailCent1Single",  //
   "kMuonCocktailPer1Single",   //
-  "kMuonCocktailPer4Single"
+  "kMuonCocktailPer4Single",
+  "kFMD1Flat",
+  "kFMD2Flat",
+  "kFMD3Flat",
+  "kFMDFlat"
 };
 
 //____________________________________________________________________
-enum PprGeo_t {
+enum Geo_t {
   kHoles,                      //
   kNoHoles                     //
 };
 
 //____________________________________________________________________
-enum PprRad_t {
+enum Rad_t {
   kGluonRadiation,             //
   kNoGluonRadiation            //
 };
 
 //____________________________________________________________________
-enum PprMag_t {
+enum Mag_t {
   k2kG,                                //
   k4kG,                                //
   k5kG                         //
@@ -195,48 +210,50 @@ enum PprMag_t {
 enum MC_t {
   kFLUKA, 
   kGEANT3, 
-  kGEANT4
+  kGEANT4, 
+  kGEANT3TGEO,
 };
 
-//____________________________________________________________________
-// This part for configuration    
-//static PprRun_t srun = test50;
-static PprRun_t srun  = kParam_8000;
-static PprGeo_t sgeo  = kNoHoles;
-static PprRad_t srad  = kGluonRadiation;
-static PprMag_t smag  = k5kG;
-static Int_t    sseed = 12345; //Set 0 to use the current time
-static MC_t     smc   = kGEANT3;
-
-
-//____________________________________________________________________
-// Comment line 
-static TString  comment;
-
 //____________________________________________________________________
 // Functions
-Float_t       EtaToTheta(Float_t arg);
-AliGenerator* GeneratorFactory(PprRun_t srun);
+Float_t       EtaToTheta(Float_t eta);
+Eg_t          LookupEG(const Char_t* name);
+AliGenerator* GeneratorFactory(EG_t eg, Rad_t rad, TString& comment);
 AliGenHijing* HijingStandard();
-void          ProcessEnvironmentVars();
+void          ProcessEnvironmentVars(EG_t& eg, Int_t& seed);
 
 //____________________________________________________________________
 void 
 Config()
 {
-  // ThetaRange is (0., 180.). It was (0.28,179.72) 7/12/00 09:00
-  // Theta range given through pseudorapidity limits 22/6/2001
-
+  //____________________________________________________________________
+  // This part for configuration    
+  // EG_t  eg   = test50;
+  // EG_t  eg   = kParam_fmd;
+  EG_t  eg   = kFMDFlat; // kParam_2000; // kPythia;
+  // EG_t  eg   = kFMDFlat;
+  Geo_t geo  = kNoHoles;
+  Rad_t rad  = kGluonRadiation;
+  Mag_t mag  = 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();
+  ProcessEnvironmentVars(eg, seed);
 
+  //____________________________________________________________________
   // Set Random Number seed
-  gRandom->SetSeed(sseed);
+  gRandom->SetSeed(seed);
   cout<<"Seed for random number generation= "<<gRandom->GetSeed()<<endl; 
 
 
   //__________________________________________________________________
-  switch (smc) {
+  switch (mc) {
   case kFLUKA: 
     // 
     // libraries required by fluka21
@@ -244,7 +261,7 @@ Config()
     gSystem->Load("libGeom");
     cout << "\t* Loading TFluka..." << endl;  
     gSystem->Load("libTFluka");    
-    
+    gSystem->MakeDirectory("peg");
     // 
     // FLUKA MC
     //
@@ -252,15 +269,34 @@ Config()
     new TFluka("C++ Interface to Fluka", 0/*verbosity*/);
     break;
   case kGEANT3: 
-    //
-    // Libraries needed by GEANT 3.21 
-    //
-    gSystem->Load("libgeant321");
+    {
+      //
+      // Libraries needed by GEANT 3.21 
+      //
+      gSystem->Load("libgeant321");
+      
+      // 
+      // GEANT 3.21 MC 
+      // 
+      TGeant3* gmc = new TGeant3("C++ Interface to Geant3");
+      gmc->SetSWIT(4, 1000);
+    }
+    break;
+  case kGEANT3TGEO:
+    {
+      //
+      // Libraries needed by GEANT 3.21 
+      //
+      gSystem->Load("EGPythia6.so");
+      gSystem->Load("libgeant321");
     
-    // 
-    // GEANT 3.21 MC 
-    // 
-    new TGeant3("C++ Interface to Geant3");
+      // 
+      // GEANT 3.21 MC 
+      // 
+      TGeant3TGeo* gmc  = new TGeant3TGeo("C++ Interface to Geant3");
+      gmc->SetSWIT(4, 1000);
+      Printf("Making a TGeant3TGeo objet");
+    }
     break;
   default:
     gAlice->Fatal("Config.C", "No MC type chosen");
@@ -284,7 +320,7 @@ Config()
 
   //__________________________________________________________________
   // For FLUKA 
-  switch (smc) {
+  switch (mc) {
   case kFLUKA: 
     {
       //
@@ -295,41 +331,51 @@ Config()
        ((TFluka*)gMC)->SetGeneratePemf(kFALSE);
       else
        ((TFluka*)gMC)->SetGeneratePemf(kTRUE);
+      TString flupro(gSystem->Getenv("FLUPRO"));
+      if (flupro.IsNull()) 
+       Fatal("Config.C", "Environment variable FLUPRO not set");
+#if 0
+      char* files[] = { "brems_fin.bin", "cohff.bin", "elasct.bin", 
+                       "gxsect.bin", "nuclear.bin", "sigmapi.bin", 
+                       0 };
+      char* file = files[0];
+      while (file) {
+       TString which(gSystem->Which(".", file));
+       if (which.IsNull()) {
+         if (gSystem->Symlink(Form("%s/%s", flupro.Data(), file), file)!=0) 
+           Fatal("Config.C", "Couldn't link $(FLUPRO)/%s -> .", file);
+       }
+       file++;
+      }
+#endif
+      TString neuxsc(gSystem->Which(".", "neuxsc.bin"));
+      if (neuxsc.IsNull()) 
+       gSystem->Symlink(Form("%s/neuxsc_72.bin", flupro.Data()), 
+                        "neuxsc.bin"); 
+      gSystem->CopyFile("$(FLUPRO)/random.dat", "old.seed", kTRUE);
     }
     break;
   }
 
+  //__________________________________________________________________
   //
   // Set External decayer
+#if 0
   AliDecayer *decayer = new AliDecayerPythia();
-
-
-  switch (srun) {
-  case kD0PbPb5500:
-    decayer->SetForceDecay(kHadronicD);
-    break;
-  case kCharmSemiElPbPb5500:
-    decayer->SetForceDecay(kSemiElectronic);
-    break;
-  case kBeautySemiElPbPb5500:
-    decayer->SetForceDecay(kSemiElectronic);
-    break;
-  default:
-    decayer->SetForceDecay(kAll);
-    break;
+  switch (eg) {
+  case kD0PbPb5500:           decayer->SetForceDecay(kHadronicD);      break;
+  case kCharmSemiElPbPb5500:  decayer->SetForceDecay(kSemiElectronic); break;
+  case kBeautySemiElPbPb5500: decayer->SetForceDecay(kSemiElectronic); break;
+  default:                    decayer->SetForceDecay(kAll);            break;
   }
   decayer->Init();
   gMC->SetExternalDecayer(decayer);
+#endif
 
-  //
-  //
-  //=======================================================================
-  //
-  //=======================================================================
-  // ************* STEERING parameters FOR ALICE SIMULATION **************
-  // --- Specify event type to be tracked through the ALICE setup
-  // --- All positions are in cm, angles in degrees, and P and E in GeV
-
+  //__________________________________________________________________
+  // *********** STEERING parameters FOR ALICE SIMULATION ************
+  // - Specify event type to be tracked through the ALICE setup
+  // - All positions are in cm, angles in degrees, and P and E in GeV 
   gMC->SetProcess("DCAY",1);
   gMC->SetProcess("PAIR",1);
   gMC->SetProcess("COMP",1);
@@ -341,7 +387,7 @@ Config()
   gMC->SetProcess("MUNU",1);
   gMC->SetProcess("CKOV",1);
   gMC->SetProcess("HADR",1);
-  gMC->SetProcess("LOSS",2);
+  gMC->SetProcess("LOSS",2); // 0:none 1,3:dray 2:nodray 4:nofluct (def:2)
   gMC->SetProcess("MULS",1);
   gMC->SetProcess("RAYL",1);
 
@@ -360,10 +406,10 @@ Config()
   gMC->SetCut("PPCUTM", cut);
   gMC->SetCut("TOFMAX", tofmax); 
 
-
+  
+  //__________________________________________________________________
   // Generator Configuration
-  gAlice->SetDebug(0);
-  AliGenerator* gener = GeneratorFactory(srun);
+  AliGenerator* gener = GeneratorFactory(eg, rad, comment);
   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
@@ -371,54 +417,59 @@ Config()
   gener->SetTrackingFlag(1);
   gener->Init();
     
+  //__________________________________________________________________
   // 
   // Comments 
   // 
-  switch (smag) {
+  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 (srad) {
+  switch (rad) {
   case kGluonRadiation: 
     comment = comment.Append(" | Gluon Radiation On");  break;
   default:
     comment = comment.Append(" | Gluon Radiation Off"); break;
   }
 
-  switch(sgeo) {
-  case kHoles: comment = comment.Append(" | Holes for PHOS/RICH"); break;
-  default:     comment = comment.Append(" | No holes for PHOS/RICH"); 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., smag);
+  AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., mag);
   field->SetL3ConstField(0); //Using const. field in the barrel
   rl->CdGAFile();
   gAlice->SetField(field);    
+  TFile* magF = TFile::Open("mag.root", "RECREATE");
+  field->Write("mag");
+  magF->Close();
 
+  //__________________________________________________________________
   // 
   // Used detectors 
   // 
   Bool_t useABSO  = kFALSE; 
-  Bool_t useCRT   = kFALSE; 
+  Bool_t useACORDE= kFALSE; 
   Bool_t useDIPO  = kFALSE; 
   Bool_t useFMD   = kTRUE; 
   Bool_t useFRAME = kFALSE; 
   Bool_t useHALL  = kFALSE; 
-  Bool_t useITS   = kFALSE; 
+  Bool_t useITS   = kFALSE;
   Bool_t useMAG   = kFALSE; 
   Bool_t useMUON  = kFALSE; 
   Bool_t usePHOS  = kFALSE; 
   Bool_t usePIPE  = kFALSE; 
   Bool_t usePMD   = kFALSE; 
-  Bool_t useRICH  = kFALSE; 
+  Bool_t useHMPID = kFALSE; 
   Bool_t useSHIL  = kFALSE; 
-  Bool_t useSTART = kFALSE; 
+  Bool_t useT0    = kFALSE; 
   Bool_t useTOF   = kFALSE; 
   Bool_t useTPC   = kFALSE;
   Bool_t useTRD   = kFALSE; 
@@ -427,57 +478,57 @@ Config()
   Bool_t useVZERO = kFALSE;
 
   cout << "\t* Creating the detectors ..." << endl;
-  //=================== Alice BODY parameters =============================
-  //=================== Alice BODY parameters =============================
+  // ================= Alice BODY parameters =========================
   AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
   
   
   if (useMAG) {
-    //=================== MAG parameters ============================
+    // =================== MAG parameters ============================
     // Start with Magnet since detector layouts may be depending on
     // the selected Magnet dimensions 
     AliMAG *MAG = new AliMAG("MAG", "Magnet");
   }
 
   if (useABSO) {
-    //=================== ABSO parameters ============================
+    // =================== ABSO parameters ===========================
     AliABSO *ABSO = new AliABSOv0("ABSO", "Muon Absorber");
   }
 
   if (useDIPO) {
-    //=================== DIPO parameters ============================
-    
+    // =================== DIPO parameters ===========================
     AliDIPO *DIPO = new AliDIPOv2("DIPO", "Dipole version 2");
   }
 
   if (useHALL) {
-    //=================== HALL parameters ============================
+    // =================== HALL parameters ===========================
     AliHALL *HALL = new AliHALL("HALL", "Alice Hall");
   }
 
 
   if (useFRAME) {
-    //=================== FRAME parameters ============================
+    // ================== FRAME parameters ===========================
     AliFRAMEv2 *FRAME = new AliFRAMEv2("FRAME", "Space Frame");
-    switch (sgeo) {
+    switch (geo) {
     case kHoles: FRAME->SetHoles(1); break;
     default:     FRAME->SetHoles(0); break;
     }
   }
 
   if (useSHIL) {
-    //=================== SHIL parameters ============================
+    // ================== SHIL parameters ============================
     AliSHIL *SHIL = new AliSHILv2("SHIL", "Shielding Version 2");
   }
 
 
   if (usePIPE) {
-    //=================== PIPE parameters ============================
+    // ================== PIPE parameters ============================
     AliPIPE *PIPE = new AliPIPEv0("PIPE", "Beam Pipe");
   }
   
   if (useITS) {
-    //=================== ITS parameters ============================
+     AliITS *ITS  = new AliITSv11Hybrid("ITS","ITS v11Hybrid");
+#if 0
+    // =================== ITS parameters ============================
     //
     // As the innermost detector in ALICE, the Inner Tracking System
     // "impacts" on almost all other detectors. This involves the fact
@@ -552,10 +603,11 @@ Config()
     // want to use this facility.
     //
     ITS->SetEUCLID(0);
+#endif
   }
 
   if (useTPC) {
-    //============================ TPC parameters ====================
+    // =================== TPC parameters ============================
     //
     // This allows the user to specify sectors for the SLOW (TPC
     // geometry 2) Simulator. SecAL (SecAU) <0 means that ALL lower
@@ -578,85 +630,84 @@ Config()
     //  gROOT->LoadMacro("SetTPCParam.C");
     //  AliTPCParam *param = SetTPCParam();
     AliTPC *TPC = new AliTPCv2("TPC", "Default");
-    
-    // All sectors included
-    TPC->SetSecAL(-1);
-    TPC->SetSecAU(-1);
   }
 
   if (useTOF) {
-    //=================== TOF parameters ============================
+    // ================== TOF parameters =============================
     AliTOF *TOF = new AliTOFv4T0("TOF", "normal TOF");
   }
 
-  if (useRICH) {
-    //=================== RICH parameters ===========================
-    AliRICH *RICH = new AliRICHv1("RICH", "normal RICH");
+  if (useHMPID) {
+    // ================== HMPID parameters ============================
+    AliHMPID *HMPID = new AliHMPIDv1("HMPID", "normal HMPID");
 
   }
 
   if (useZDC) {
-    //=================== ZDC parameters ============================
+    // ================== ZDC parameters =============================
     AliZDC *ZDC = new AliZDCv2("ZDC", "normal ZDC");
   }
 
   if (useTRD) {
-    //=================== TRD parameters ============================
+    // ================== 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);
-    if (sgeo == kHoles) {
+    if (geo == kHoles) {
       // With hole in front of PHOS
       TRD->SetPHOShole();
-      // With hole in front of RICH
-      TRD->SetRICHhole();
+      // With hole in front of HMPID
+      TRD->SetHMPIDhole();
     }
     // Switch on TR
     AliTRDsim *TRDsim = TRD->CreateTR();
   }
 
   if (useFMD) {
-    //=================== FMD parameters ============================
+    // =================== FMD parameters ============================
     AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
+    // FMD->UseDetailed(kFALSE);
+    // FMD->UseAssembly();
+    // FMD->UseOld();
   }
 
   if (useMUON) {
-    //=================== MUON parameters ===========================
+    // =================== MUON parameters ===========================
     AliMUON *MUON = new AliMUONv1("MUON", "default");
-    MUON->AddGeometryBuilder(new AliMUONSt1GeometryBuilder(MUON));
-    MUON->AddGeometryBuilder(new AliMUONSt2GeometryBuilder(MUON));
-    MUON->AddGeometryBuilder(new AliMUONSlatGeometryBuilder(MUON));
-    MUON->AddGeometryBuilder(new AliMUONTriggerGeometryBuilder(MUON));
+    // MUON->AddGeometryBuilder(new AliMUONSt1GeometryBuilder(MUON));
+    // MUON->AddGeometryBuilder(new AliMUONSt2GeometryBuilder(MUON));
+    // MUON->AddGeometryBuilder(new AliMUONSlatGeometryBuilder(MUON));
+    // MUON->AddGeometryBuilder(new AliMUONTriggerGeometryBuilder(MUON));
   }
 
   if (usePHOS) {
-    //=================== PHOS parameters ===========================
+    // =================== PHOS parameters ===========================
     AliPHOS *PHOS = new AliPHOSv1("PHOS", "IHEP");
   }
 
   if (usePMD) {
-    //=================== PMD parameters ============================
+    // =================== PMD parameters ============================
     AliPMD *PMD = new AliPMDv1("PMD", "normal PMD");
   }
 
-  if (useSTART) {
-    //=================== START parameters ============================
-    AliSTART *START = new AliSTARTv1("START", "START Detector");
+  if (useT0) {
+    // =================== T0 parameters ==========================
+    AliT0 *T0 = new AliT0v1("T0", "T0 Detector");
   }
 
   if (useEMCAL) {
-    //=================== EMCAL parameters ============================
-    AliEMCAL *EMCAL = new AliEMCALv1("EMCAL", "EMCAL_55_25");
+    // =================== EMCAL parameters ==========================
+    AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETE");
   }
 
-  if (useCRT) {
-    //=================== CRT parameters ============================
-    AliCRT *CRT = new AliCRTv0("CRT", "normal ACORDE");
+  if (useACORDE) {
+    // =================== ACORDE parameters ============================
+    AliACORDE *ACORDE = new AliACORDEv1("ACORDE", "normal ACORDE");
   }
 
   if (useVZERO) {
-    //=================== CRT parameters ============================
+    // =================== V0 parameters =============================
     AliVZERO *VZERO = new AliVZEROv3("VZERO", "normal VZERO");
   }
 }
@@ -664,20 +715,30 @@ Config()
 //____________________________________________________________________
 Float_t EtaToTheta(Float_t arg)
 {
-  return (180./TMath::Pi())*2.*atan(exp(-arg));
+  return (180./TMath::Pi())*2.*TMath::ATan(TMath::Exp(-arg));
 }
 
-
 //____________________________________________________________________
+Int_t 
+LookupEG(const Char_t* name) 
+{
+  TString n(name);
+  for (Int_t i = 0; i < kEgMax; i++) {
+    if (n == egName[i]) return i;
+  }
+  return -1;
+}
+
+//____________________________________________________________________  
 AliGenerator* 
-GeneratorFactory(PprRun_t srun) 
+GeneratorFactory(EG_t eg, Rad_t rad, TString& comment)  
 {
   Int_t isw = 3;
-  if (srad == kNoGluonRadiation) isw = 0;
+  if (rad == kNoGluonRadiation) isw = 0;
   
   
   AliGenerator * gGener = 0;
-  switch (srun) {
+  switch (eg) {
   case test50:
     {
       comment = comment.Append(":HIJINGparam test 50 particles");
@@ -730,6 +791,19 @@ GeneratorFactory(PprRun_t srun)
       gGener=gener;
     }
     break;
+  case kParam_fmd:
+    {
+      comment = comment.Append("HIJINGparam N=100");
+      AliGenHIJINGpara *gener = new AliGenHIJINGpara(500);
+      gener->SetMomentumRange(0, 999999.);
+      gener->SetPhiRange(0., 360.);
+      // Set pseudorapidity range from -8 to 8.
+      Float_t thmin = EtaToTheta(6);   // theta min. <---> eta max
+      Float_t thmax = EtaToTheta(2);  // theta max. <---> eta min 
+      gener->SetThetaRange(thmin,thmax);
+      gGener=gener;
+    }
+    break;
     //
     //  Hijing Central
     //
@@ -1510,6 +1584,61 @@ GeneratorFactory(PprRun_t srun)
       gGener=gener;
     }
     break;
+  case kFMD1Flat: 
+    {
+      comment = comment.Append(" Flat in FMD1 range");
+      AliGenBox* gener = new AliGenBox(2000);
+      gener->SetPart(kPiPlus);
+      gener->SetMomentumRange(3,4);
+      gener->SetPhiRange(0, 360);
+      gener->SetThetaRange(0.77, 3.08);
+      gGener = gener;
+    }
+    break;
+  case kFMD2Flat: 
+    {
+      comment = comment.Append(" Flat in FMD2 range");
+      AliGenBox* gener = new AliGenBox(2000);
+      gener->SetPart(kPiPlus);
+      gener->SetMomentumRange(3,4);
+      gener->SetPhiRange(0, 360);
+      gener->SetThetaRange(2.95, 20.42);
+      gGener = gener;
+    }
+    break;
+  case kFMD3Flat: 
+    {
+      comment = comment.Append(" Flat in FMD3 range");
+      AliGenBox* gener = new AliGenBox(2000);
+      gener->SetPart(kPiPlus);
+      gener->SetMomentumRange(3,4);
+      gener->SetPhiRange(0, 360);
+      gener->SetThetaRange(155.97, 176.73);
+      gGener = gener;
+    }
+    break;
+  case kFMDFlat:
+    {
+      comment = comment.Append(" Flat in FMD range");
+      AliGenCocktail* gener = new AliGenCocktail();
+      gener->SetMomentumRange(3,4);
+      gener->SetPhiRange(0, 360);
+      AliGenBox* gener3 = new AliGenBox(2000);
+      gener3->SetThetaRange(155.97, 176.73);
+      gener3->SetPart(kPiPlus);
+      gener->AddGenerator(gener3, "FMD3", .33);
+      AliGenBox* gener2 = new AliGenBox(2000);
+      gener2->SetThetaRange(2.95, 20.42);
+      gener2->SetPart(kPiPlus);
+      gener->AddGenerator(gener2, "FMD2", .33);
+      AliGenBox* gener1 = new AliGenBox(2000);
+      gener1->SetThetaRange(0.77, 3.08);
+      gener1->SetPart(kPiPlus);
+      gener->AddGenerator(gener1, "FMD1", .34);
+      gGener = gener;
+    }
+    break;
+    
   default: break;
   }
   return gGener;
@@ -1545,21 +1674,16 @@ HijingStandard()
 
 //____________________________________________________________________
 void 
-ProcessEnvironmentVars()
+ProcessEnvironmentVars(EG_t& eg, Int_t& seed)
 {
   // Run type
   if (gSystem->Getenv("CONFIG_RUN_TYPE")) {
-    for (Int_t iRun = 0; iRun < kRunMax; iRun++) {
-      if (strcmp(gSystem->Getenv("CONFIG_RUN_TYPE"), pprRunName[iRun])==0) {
-       srun = (PprRun_t)iRun;
-       cout<<"Run type set to "<<pprRunName[iRun]<<endl;
-      }
-    }
+    Int_t eg1 = LookupEG(gSystem->Getenv("CONFIG_RUN_TYPE"));
+    if  (eg1 >= 0) eg = EG_t(eg1);
   }
-
   // Random Number seed
   if (gSystem->Getenv("CONFIG_SEED")) {
-    sseed = atoi(gSystem->Getenv("CONFIG_SEED"));
+    seed = atoi(gSystem->Getenv("CONFIG_SEED"));
   }
 }