Updated alignment macros and related changes (Raffaele)
[u/mrichter/AliRoot.git] / macros / g4ConfigCommon.C
index d06d335..f8ccfc3 100644 (file)
 // AliRoot Configuration for running aliroot with Monte Carlo.
 // Called from g4Config.C
 
-static Int_t    eventsPerRun = 100;
-enum PprGeo_t 
-{
-    kHoles, kNoHoles
-};
-static PprGeo_t geo = kHoles;
+Float_t EtaToTheta(Float_t arg);
+static Int_t    eventsPerRun = 50;
 
-void ConfigCommon(Bool_t interactiveSetup)
+void ConfigCommon()
 {
   // ============================= 
   // Root file
   // ============================= 
 
   // 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);
-
-  // ============================= 
-  // Event generator
-  // ============================= 
-
-  // --- 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();
-  //
-  // Activate this line if you want the vertex smearing to happen
-  // track by track
+  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);
+//  gAlice->SetRootGeometry();
+//  gAlice->SetGeometryFileName("geometry.root");
+
+  // Set the trigger configuration
+  gAlice->SetTriggerDescriptor("Pb-Pb");
+  cout<<"Trigger configuration is set to  Pb-Pb"<<endl;
+
+    // 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 = 500;
+    if (gSystem->Getenv("CONFIG_NPARTICLES"))
+    {
+        nParticles = atoi(gSystem->Getenv("CONFIG_NPARTICLES"));
+    }
 
-  //gener->SetVertexSmear(perTrack); 
+///*
+    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); 
 
   // ============================= 
   // Magnetic field
   // ============================= 
 
-  // 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)
-
-  // ============================= 
-  // Alice modules
-  // ============================= 
-
-  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)
+    AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., 1);
+    gAlice->SetField(field);    
+
+    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;
+
+    rl->CdGAFile();
     //=================== Alice BODY parameters =============================
     AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
 
-
     if (iMAG)
     {
         //=================== MAG parameters ============================
@@ -155,21 +107,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,18 +130,13 @@ void ConfigCommon(Bool_t interactiveSetup)
         //=================== FRAME parameters ============================
 
         AliFRAMEv2 *FRAME = new AliFRAMEv2("FRAME", "Space Frame");
-       if (geo == kHoles) {
-           FRAME->SetHoles(1);
-       } else {
-           FRAME->SetHoles(0);
-       }
     }
 
     if (iSHIL)
     {
         //=================== SHIL parameters ============================
 
-        AliSHIL *SHIL = new AliSHILv2("SHIL", "Shielding");
+        AliSHIL *SHIL = new AliSHILv3("SHIL", "Shielding Version 3");
     }
 
 
@@ -197,12 +144,12 @@ 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
@@ -220,38 +167,24 @@ void ConfigCommon(Bool_t interactiveSetup)
     //
     //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
+       AliITSvPPRasymmFMD *ITS  = new AliITSvPPRasymmFMD("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
-    //
-    //
+       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(150.);  // chip thickness on layer 1 must be in the range [150,300]
+       ITS->SetThicknessChip2(150.);  // 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
+    //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
@@ -267,55 +200,33 @@ void ConfigCommon(Bool_t interactiveSetup)
     // 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);  
-  }
-  
+       ITS->SetEUCLID(0);  
+    }
 
     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");
+       // Partial geometry: modules at 2,3,4,6,7,11,12,14,15,16
+       // starting at 6h in positive direction
+       //      Int_t TOFSectors[18]={-1,-1,0,0,0,-1,0,0,-1,-1,-1,0,0,-1,0,0,0,0};
+       // Partial geometry: modules at 1,2,6,7,9,10,11,12,15,16,17
+       // (ALICE numbering convention)
+               Int_t TOFSectors[18]={-1,0,0,-1,-1,-1,0,0,-1,0,0,0,0,-1,-1,0,0,0};
+       TOF->SetTOFSectors(TOFSectors);
     }
 
-    if (iRICH)
+
+    if (iHMPID)
     {
-        //=================== RICH parameters ===========================
-/*
-        AliRICH *RICH = new AliRICHv3("RICH", "normal RICH");
-*/
-        AliRICH *RICH = new AliRICHv1("RICH", "normal RICH");
+        //=================== HMPID parameters ===========================
+        AliHMPID *HMPID = new AliHMPIDv2("HMPID", "normal HMPID");
 
     }
 
@@ -332,41 +243,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");
     }
 
 
@@ -376,39 +271,32 @@ 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", "SHISH_77_TRD1_2X2_FINAL_110DEG");
     }
 
-    if (iCRT)
+     if (iACORDE)
     {
-        //=================== CRT parameters ============================
-
-        AliCRT *CRT = new AliCRTv0("CRT", "normal ACORDE");
+        //=================== ACORDE parameters ============================
+        AliACORDE *ACORDE = new AliACORDEv0("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__);
+
 }
 
 Float_t EtaToTheta(Float_t arg){