]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - macros/g4ConfigCommon.C
PAR: includes from previously enabled PARfiles
[u/mrichter/AliRoot.git] / macros / g4ConfigCommon.C
index edc27f36156ca014def1a67ef1b5b4a9a4c94c9c..2dcf4002347046145a179f08ccba07e6054ae1bd 100644 (file)
 //
 // AliRoot Configuration for running aliroot with Monte Carlo.
 // Called from g4Config.C
+//
+// By I. Hrivnacova, IPN Orsay
 
-static Int_t    eventsPerRun = 100;
-enum PprGeo_t 
-{
+enum PprGeo_t
+  {
     kHoles, kNoHoles
-};
+  };
 static PprGeo_t geo = kHoles;
 
-void ConfigCommon(Bool_t interactiveSetup)
-{
-  // ============================= 
-  // Root file
-  // ============================= 
+Float_t EtaToTheta(Float_t arg);
+void    LoadPythia();
+static Int_t    eventsPerRun = 50;
 
-  // Create the output file
-  TFile *rootfile = new TFile("galice.root","recreate");
-  rootfile->SetCompressionLevel(2);
-
-  // Set External decayer
-  AliDecayer* decayer = new AliDecayerPythia();
-  decayer->SetForceDecay(kAll);
-  decayer->Init();
-  gMC->SetExternalDecayer(decayer);
-
-  // Physics process control
-  gMC ->SetProcess("DCAY",1);
-  gMC ->SetProcess("PAIR",1);
-  gMC ->SetProcess("COMP",1);
-  gMC ->SetProcess("PHOT",1);
-  gMC ->SetProcess("PFIS",0);
-  gMC ->SetProcess("DRAY",0);
-  gMC ->SetProcess("ANNI",1);
-  gMC ->SetProcess("BREM",1);
-  gMC ->SetProcess("MUNU",1);
-  //xx gMC ->SetProcess("CKOV",1);
-  gMC ->SetProcess("HADR",1); //Select pure GEANH (HADR 1) or GEANH/NUCRIN (HADR 3)
-  gMC ->SetProcess("LOSS",2);
-  gMC ->SetProcess("MULS",1);
-  //xx gMC ->SetProcess("RAYL",1);
-
-  // Energy cuts
-  // (in development)
-  Float_t cut    = 1.e-3; // 1MeV cut by default
-  Float_t tofmax = 1.e10; 
-
-  gMC ->SetCut("CUTGAM",cut);
-  gMC ->SetCut("CUTELE",cut);
-  gMC ->SetCut("CUTNEU",cut);
-  gMC ->SetCut("CUTHAD",cut);
-  gMC ->SetCut("CUTMUO",cut);
-  gMC ->SetCut("BCUTE",cut);
-  gMC ->SetCut("BCUTM",cut);
-  gMC ->SetCut("DCUTE",cut);
-  gMC ->SetCut("DCUTM",cut);
-  gMC ->SetCut("PPCUTM",cut);
-  gMC ->SetCut("TOFMAX",tofmax);
+void ConfigCommon(Bool_t setRootGeometry = kTRUE)
+{
+  cout << "Running ConfigCommon.C ... " << endl;
 
-  // ============================= 
-  // Event generator
-  // ============================= 
+  // Load Pythia libraries
+  LoadPythia();
 
-  // --- Specify event type to be tracked through the ALICE setup
-  // --- All positions are in cm, angles in degrees, and P and E in GeV
-
-  Int_t nParticles;
-  if (gSystem->Getenv("CONFIG_NPARTICLES")) 
-    nParticles = atoi(gSystem->Getenv("CONFIG_NPARTICLES"));
-  else
-    nParticles = 100; 
-
-
-  AliGenHIJINGpara *gener = new AliGenHIJINGpara(nParticles);
-  gener->SetMomentumRange(0,999);
-  gener->SetPhiRange(0,360);
-  // Set pseudorapidity range from -8 to 8.
-  Float_t thmin = EtaToTheta(8);   // theta min. <---> eta max
-  Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min 
-  gener->SetThetaRange(thmin, thmax);
-  gener->SetOrigin(0,0,0);        //vertex position
-  gener->SetSigma(0,0,0);         //Sigma in (X,Y,Z) (cm) on IP position
-  gener->Init();
+  // Create the output file
+  AliRunLoader* rl = 0;
+  rl = AliRunLoader::Open("galice.root",
+                          AliConfig::GetDefaultEventFolderName(),
+                          "recreate");
+  if (!rl) {
+    gAlice->Fatal("Config.C","Can not instatiate the Run Loader");
+    return;
+  }
+  rl->SetCompressionLevel(2);
+  rl->SetNumberOfEventsPerFile(3);
+  gAlice->SetRunLoader(rl);
+  
+  // Set Root geometry file
   //
-  // Activate this line if you want the vertex smearing to happen
-  // track by track
+  if ( setRootGeometry ) {
+    gAlice->SetRootGeometry();
+    gAlice->SetGeometryFromFile("geometry.root");
+  }
 
-  //gener->SetVertexSmear(perTrack); 
+  // Set the trigger configuration
+  AliSimulation::Instance()->SetTriggerConfig("Pb-Pb");
+  cout<<"Trigger configuration is set to  Pb-Pb"<<endl;
 
-  // ============================= 
-  // Magnetic field
-  // ============================= 
+    // Set Random Number seed
+    gRandom->SetSeed(123456); // Set 0 to use the currecnt time
+    AliLog::Message(AliLog::kInfo, Form("Seed for random number generation = %d",gRandom->GetSeed()), "Config.C", "Config.C", "Config()","Config.C", __LINE__);
+    int     nParticles = 100;
+    if (gSystem->Getenv("CONFIG_NPARTICLES"))
+    {
+        nParticles = atoi(gSystem->Getenv("CONFIG_NPARTICLES"));
+    }
 
-  // Field (L3 0.4 T)
-  AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., 1);
-  rootfile->cd();
-  gAlice->SetField(field);    
-  
-  // Old magnetic field
-  //gAlice->SetField(-999,2);    //Specify maximum magnetic field in Tesla (neg. ==> default field)
+    AliGenCocktail *gener = new AliGenCocktail();
+    gener->SetPhiRange(0, 360);
+    // Set pseudorapidity range from -8 to 8.
+    Float_t thmin = EtaToTheta(8);   // theta min. <---> eta max
+    Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min 
+    gener->SetThetaRange(thmin,thmax);
+    gener->SetOrigin(0, 0, 0);  //vertex position
+    gener->SetSigma(0, 0, 0);   //Sigma in (X,Y,Z) (cm) on IP position
+
+    AliGenHIJINGpara *hijingparam = new AliGenHIJINGpara(nParticles);
+    hijingparam->SetMomentumRange(0.2, 999);
+    gener->AddGenerator(hijingparam,"HIJING PARAM",1);
+    gener->Init();
+
+    // 
+    // Activate this line if you want the vertex smearing to happen
+    // track by track
+    //
+    //gener->SetVertexSmear(perTrack); 
 
   // ============================= 
-  // Alice modules
+  // Magnetic field
   // ============================= 
 
-  if (!interactiveSetup) {
-
-    // Select modules 
-  Int_t iABSO=1;
-  Int_t iDIPO=1;
-  Int_t iFMD=1;
-  Int_t iFRAME=1;
-  Int_t iHALL=1;
-  Int_t iITS=1;
-  Int_t iMAG=1;
-  Int_t iMUON=1;
-  Int_t iPHOS=1;
-  Int_t iPIPE=1;
-  Int_t iPMD=1;
-  Int_t iRICH=1;
-  Int_t iSHIL=1;
-  Int_t iSTART=1;
-  Int_t iTOF=1;
-  Int_t iTPC=1;
-  Int_t iTRD=1;
-  Int_t iZDC=1;
-  Int_t iEMCAL = 1;
-  Int_t iCRT = 0;  
-  Int_t iVZERO = 1;
-
-  // ONLY FOR GEANT4
-  
-  // Exclude detectors with temporary problem
-  iCRT = 0;
-  iEMCAL = 0;
-  // END OF ONLY FOR GEANT4
-
-
+    // Field (L3 0.4 T)
+    TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
+
+    Int_t   iABSO  =  1;
+    Int_t   iDIPO  =  1;
+    Int_t   iFMD   =  1;
+    Int_t   iFRAME =  1;
+    Int_t   iHALL  =  1;
+    Int_t   iITS   =  1;
+    Int_t   iMAG   =  1;
+    Int_t   iMUON  =  1;
+    Int_t   iPHOS  =  1;
+    Int_t   iPIPE  =  1;
+    Int_t   iPMD   =  1;
+    Int_t   iHMPID =  1;
+    Int_t   iSHIL  =  1;
+    Int_t   iT0    =  1;
+    Int_t   iTOF   =  1;
+    Int_t   iTPC   =  1;
+    Int_t   iTRD   =  1;
+    Int_t   iZDC   =  1;
+    Int_t   iEMCAL =  1;
+    Int_t   iACORDE = 0;
+    Int_t   iVZERO =  1;
+/*
+    Int_t   iABSO  =  0;
+    Int_t   iDIPO  =  0;
+    Int_t   iFMD   =  0;
+    Int_t   iFRAME =  0;
+    Int_t   iHALL  =  0;
+    Int_t   iITS   =  0;
+    Int_t   iMAG   =  0;
+    Int_t   iMUON  =  0;
+    Int_t   iPHOS  =  0;
+    Int_t   iPIPE  =  0;
+    Int_t   iPMD   =  0;
+    Int_t   iHMPID =  0;
+    Int_t   iSHIL  =  0;
+    Int_t   iT0    =  0;
+    Int_t   iTOF   =  0;
+    Int_t   iTPC   =  1;
+    Int_t   iTRD   =  0;
+    Int_t   iZDC   =  0;
+    Int_t   iEMCAL =  0;
+    Int_t   iACORDE = 0;
+    Int_t   iVZERO =  0;
+*/   
+// Exluded for problems
+
+    rl->CdGAFile();
     //=================== Alice BODY parameters =============================
     AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
 
-
     if (iMAG)
     {
         //=================== MAG parameters ============================
@@ -155,21 +144,21 @@ void ConfigCommon(Bool_t interactiveSetup)
     if (iABSO)
     {
         //=================== ABSO parameters ============================
-        AliABSO *ABSO = new AliABSOv0("ABSO", "Muon Absorber");
+        AliABSO *ABSO = new AliABSOv3("ABSO", "Muon Absorber");
     }
 
     if (iDIPO)
     {
         //=================== DIPO parameters ============================
 
-        AliDIPO *DIPO = new AliDIPOv2("DIPO", "Dipole version 2");
+        AliDIPO *DIPO = new AliDIPOv3("DIPO", "Dipole version 3");
     }
 
     if (iHALL)
     {
         //=================== HALL parameters ============================
 
-        AliHALL *HALL = new AliHALL("HALL", "Alice Hall");
+        AliHALL *HALL = new AliHALLv3("HALL", "Alice Hall");
     }
 
 
@@ -178,10 +167,10 @@ void ConfigCommon(Bool_t interactiveSetup)
         //=================== FRAME parameters ============================
 
         AliFRAMEv2 *FRAME = new AliFRAMEv2("FRAME", "Space Frame");
-       if (geo == kHoles) {
-           FRAME->SetHoles(1);
+        if (geo == kHoles) {
+         FRAME->SetHoles(1);
        } else {
-           FRAME->SetHoles(0);
+         FRAME->SetHoles(0);
        }
     }
 
@@ -189,7 +178,7 @@ void ConfigCommon(Bool_t interactiveSetup)
     {
         //=================== SHIL parameters ============================
 
-        AliSHIL *SHIL = new AliSHILv2("SHIL", "Shielding");
+        AliSHIL *SHIL = new AliSHILv3("SHIL", "Shielding Version 3");
     }
 
 
@@ -197,122 +186,34 @@ void ConfigCommon(Bool_t interactiveSetup)
     {
         //=================== PIPE parameters ============================
 
-        AliPIPE *PIPE = new AliPIPEv0("PIPE", "Beam Pipe");
+        AliPIPE *PIPE = new AliPIPEv3("PIPE", "Beam Pipe");
     }
  
-  if(iITS) {
+    if (iITS)
+    {
+        //=================== ITS parameters ============================
 
-//=================== ITS parameters ============================
-    //
-    // As the innermost detector in ALICE, the Inner Tracking System "impacts" on
-    // almost all other detectors. This involves the fact that the ITS geometry
-    // still has several options to be followed in parallel in order to determine
-    // the best set-up which minimizes the induced background. All the geometries
-    // available to date are described in the following. Read carefully the comments
-    // and use the default version (the only one uncommented) unless you are making
-    // comparisons and you know what you are doing. In this case just uncomment the
-    // ITS geometry you want to use and run Aliroot.
-    //
-    // Detailed geometries:         
-    //
-    //
-    //AliITS *ITS  = new AliITSv5symm("ITS","Updated ITS TDR detailed version with symmetric services");
-    //
-    //AliITS *ITS  = new AliITSv5asymm("ITS","Updates ITS TDR detailed version with asymmetric services");
-    //
-/*
-    AliITSvPPRasymm *ITS  = new AliITSvPPRasymm("ITS","New ITS PPR detailed version with asymmetric services");
-    ITS->SetMinorVersion(2);                                    // don't touch this parameter if you're not an ITS developer
-    ITS->SetReadDet(kFALSE);                                    // don't touch this parameter if you're not an ITS developer
-    //    ITS->SetWriteDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymm2.det");  // don't touch this parameter if you're not an ITS developer
-    ITS->SetThicknessDet1(200.);   // detector thickness on layer 1 must be in the range [100,300]
-    ITS->SetThicknessDet2(200.);   // detector thickness on layer 2 must be in the range [100,300]
-    ITS->SetThicknessChip1(200.);  // chip thickness on layer 1 must be in the range [150,300]
-    ITS->SetThicknessChip2(200.);  // chip thickness on layer 2 must be in the range [150,300]
-    ITS->SetRails(0);       // 1 --> rails in ; 0 --> rails out
-    ITS->SetCoolingFluid(1);   // 1 --> water ; 0 --> freon
-    //
-    //AliITSvPPRsymm *ITS  = new AliITSvPPRsymm("ITS","New ITS PPR detailed version with symmetric services");
-    //ITS->SetMinorVersion(2);                                       // don't touch this parameter if you're not an ITS developer
-    //ITS->SetReadDet(kFALSE);                                       // don't touch this parameter if you're not an ITS developer
-    //ITS->SetWriteDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRsymm2.det"); // don't touch this parameter if you're not an ITS developer
-    //ITS->SetThicknessDet1(200.);   // detector thickness on layer 1 must be in the range [100,300]
-    //ITS->SetThicknessDet2(200.);   // detector thickness on layer 2 must be in the range [100,300]
-    //ITS->SetThicknessChip1(200.);  // chip thickness on layer 1 must be in the range [150,300]
-    //ITS->SetThicknessChip2(200.);  // chip thickness on layer 2 must be in the range [150,300]
-    //ITS->SetRails(0);              // 1 --> rails in ; 0 --> rails out
-    //ITS->SetCoolingFluid(1);       // 1 --> water ; 0 --> freon
-    //
-    //
-    // Coarse geometries (warning: no hits are produced with these coarse geometries and they unuseful 
-    // for reconstruction !):
-    //                                                     
-    //
-*/
-    AliITSvPPRcoarseasymm *ITS  = new AliITSvPPRcoarseasymm("ITS","New ITS PPR coarse version with asymmetric services");
-    ITS->SetRails(0);                // 1 --> rails in ; 0 --> rails out
-    ITS->SetSupportMaterial(0);      // 0 --> Copper ; 1 --> Aluminum ; 2 --> Carbon
-    //
-    //AliITS *ITS  = new AliITSvPPRcoarsesymm("ITS","New ITS PPR coarse version with symmetric services");
-    //ITS->SetRails(0);                // 1 --> rails in ; 0 --> rails out
-    //ITS->SetSupportMaterial(0);      // 0 --> Copper ; 1 --> Aluminum ; 2 --> Carbon
-    //                      
-    //
-    //
-    // Geant3 <-> EUCLID conversion
-    // ============================
-    //
-    // SetEUCLID is a flag to output (=1) or not to output (=0) both geometry and
-    // media to two ASCII files (called by default ITSgeometry.euc and
-    // ITSgeometry.tme) in a format understandable to the CAD system EUCLID.
-    // The default (=0) means that you dont want to use this facility.
-    //
-    ITS->SetEUCLID(0);  
-  }
-  
+       AliITS *ITS  = new AliITSv11("ITS","ITS v11");
+    }
 
     if (iTPC)
     {
-        //============================ TPC parameters ================================
-        // --- This allows the user to specify sectors for the SLOW (TPC geometry 2)
-        // --- Simulator. SecAL (SecAU) <0 means that ALL lower (upper)
-        // --- sectors are specified, any value other than that requires at least one 
-        // --- sector (lower or upper)to be specified!
-        // --- Reminder: sectors 1-24 are lower sectors (1-12 -> z>0, 13-24 -> z<0)
-        // ---           sectors 25-72 are the upper ones (25-48 -> z>0, 49-72 -> z<0)
-        // --- SecLows - number of lower sectors specified (up to 6)
-        // --- SecUps - number of upper sectors specified (up to 12)
-        // --- Sens - sensitive strips for the Slow Simulator !!!
-        // --- This does NOT work if all S or L-sectors are specified, i.e.
-        // --- if SecAL or SecAU < 0
-        //
-        //
-        //-----------------------------------------------------------------------------
-
-        //  gROOT->LoadMacro("SetTPCParam.C");
-        //  AliTPCParam *param = SetTPCParam();
+        //============================ TPC parameters ===================
         AliTPC *TPC = new AliTPCv2("TPC", "Default");
-
-        // All sectors included 
-        TPC->SetSecAL(-1);
-        TPC->SetSecAU(-1);
-
     }
 
-    if (iTOF)
-    {
-       if (geo == kHoles) {
+
+    if (iTOF) {
         //=================== TOF parameters ============================
-           AliTOF *TOF = new AliTOFv2FHoles("TOF", "TOF with Holes");
-       } else {
-           AliTOF *TOF = new AliTOFv4T0("TOF", "normal TOF");
-       }
+       AliTOF *TOF = new AliTOFv6T0("TOF", "normal TOF");
     }
 
-    if (iRICH)
+
+    if (iHMPID)
     {
-        //=================== RICH parameters ===========================
-        AliRICH *RICH = new AliRICHv3("RICH", "normal RICH");
+        //=================== HMPID parameters ===========================
+        AliHMPID *HMPID = new AliHMPIDv3("HMPID", "normal HMPID");
+
     }
 
 
@@ -320,7 +221,7 @@ void ConfigCommon(Bool_t interactiveSetup)
     {
         //=================== ZDC parameters ============================
 
-        AliZDC *ZDC = new AliZDCv2("ZDC", "normal ZDC");
+        AliZDC *ZDC = new AliZDCv3("ZDC", "normal ZDC");
     }
 
     if (iTRD)
@@ -328,41 +229,25 @@ void ConfigCommon(Bool_t interactiveSetup)
         //=================== 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 (geo == kHoles) {
-           // With hole in front of PHOS
-           TRD->SetPHOShole();
-           // With hole in front of RICH
-           TRD->SetRICHhole();
-       }
-        // Switch on TR
-        AliTRDsim *TRDsim = TRD->CreateTR();
     }
 
     if (iFMD)
     {
         //=================== FMD parameters ============================
-
-        AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
-        FMD->SetRingsSi1(256);
-        FMD->SetRingsSi2(128);
-        FMD->SetSectorsSi1(20);
-        FMD->SetSectorsSi2(24);
+       AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
    }
 
     if (iMUON)
     {
         //=================== MUON parameters ===========================
-
+        // New MUONv1 version (geometry defined via builders)
         AliMUON *MUON = new AliMUONv1("MUON", "default");
     }
     //=================== PHOS parameters ===========================
 
     if (iPHOS)
     {
-        AliPHOS *PHOS = new AliPHOSv1("PHOS", "GPS2");
+        AliPHOS *PHOS = new AliPHOSv1("PHOS", "IHEP");
     }
 
 
@@ -372,41 +257,45 @@ void ConfigCommon(Bool_t interactiveSetup)
         AliPMD *PMD = new AliPMDv1("PMD", "normal PMD");
     }
 
-    if (iSTART)
+    if (iT0)
     {
-        //=================== START parameters ============================
-        AliSTART *START = new AliSTARTv1("START", "START Detector");
+        //=================== T0 parameters ============================
+        AliT0 *T0 = new AliT0v1("T0", "T0 Detector");
     }
 
     if (iEMCAL)
     {
         //=================== EMCAL parameters ============================
-        AliEMCAL *EMCAL = new AliEMCALv1("EMCAL", "EMCALArch1a");
+        AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETE");
     }
 
-    if (iCRT)
+     if (iACORDE)
     {
-        //=================== CRT parameters ============================
-
-        AliCRT *CRT = new AliCRTv0("CRT", "normal ACORDE");
+        //=================== ACORDE parameters ============================
+        AliACORDE *ACORDE = new AliACORDEv1("ACORDE", "normal ACORDE");
     }
 
-    if (iVZERO)
+     if (iVZERO)
     {
-        //=================== CRT parameters ============================
-        AliVZERO *VZERO = new AliVZEROv2("VZERO", "normal VZERO");
+        //=================== ACORDE parameters ============================
+        AliVZERO *VZERO = new AliVZEROv7("VZERO", "normal VZERO");
     }
 
-  } // end (!isSetInteractively)
-  else {
-  
-    if (geo == kHoles) 
-      AliRunConfiguration::SwitchHoles(kTRUE);
-    else  
-      AliRunConfiguration::SwitchHoles(kFALSE);
-  }    
+     AliLog::Message(AliLog::kInfo, "End of Config", "Config.C", "Config.C", "Config()"," Config.C", __LINE__);
+
+     cout << "Running ConfigCommon.C finished ... " << endl;
+
 }
 
 Float_t EtaToTheta(Float_t arg){
   return (180./TMath::Pi())*2.*atan(exp(-arg));
 }
+
+void LoadPythia()
+{
+    // Load Pythia related libraries
+    gSystem->Load("liblhapdf");      // Parton density functions
+    gSystem->Load("libEGPythia6");   // TGenerator interface
+    gSystem->Load("libpythia6");     // Pythia
+    gSystem->Load("libAliPythia6");  // ALICE specific implementations
+}