Updated alignment macros and related changes (Raffaele)
[u/mrichter/AliRoot.git] / macros / g4ConfigCommon.C
index e718067..f8ccfc3 100644 (file)
@@ -4,9 +4,9 @@
 // Called from g4Config.C
 
 Float_t EtaToTheta(Float_t arg);
-static Int_t    eventsPerRun = 100;
+static Int_t    eventsPerRun = 50;
 
-void ConfigCommon(Bool_t interactiveSetup)
+void ConfigCommon()
 {
   // ============================= 
   // Root file
@@ -24,134 +24,77 @@ void ConfigCommon(Bool_t interactiveSetup)
   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"));
+    }
 
-  // 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 = 100;
-  if (gSystem->Getenv("CONFIG_NPARTICLES")) 
-    nParticles = atoi(gSystem->Getenv("CONFIG_NPARTICLES"));
-
-  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);
-
-  //    AliGenBox *genbox = new AliGenBox(nParticles);
-  //    genbox->SetPart(22);
-  //    genbox->SetPtRange(0.3, 10.00);
-  //    gener->AddGenerator(genbox,"GENBOX GAMMA for PHOS",1);
-  gener->Init();
-
-  // Activate this line if you want the vertex smearing to happen
-  // track by track
-
-  //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);
-  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 ============================
@@ -164,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");
     }
 
 
@@ -193,7 +136,7 @@ void ConfigCommon(Bool_t interactiveSetup)
     {
         //=================== SHIL parameters ============================
 
-        AliSHIL *SHIL = new AliSHILv2("SHIL", "Shielding Version 2");
+        AliSHIL *SHIL = new AliSHILv3("SHIL", "Shielding Version 3");
     }
 
 
@@ -201,63 +144,64 @@ void ConfigCommon(Bool_t interactiveSetup)
     {
         //=================== PIPE parameters ============================
 
-        AliPIPE *PIPE = new AliPIPEv0("PIPE", "Beam Pipe");
+        AliPIPE *PIPE = new AliPIPEv3("PIPE", "Beam Pipe");
     }
  
-    if(iITS) 
-    {
-      //=================== 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:         
-      //
-      //
-      //
-/*
-       AliITSvPPRasymmFMD *ITS  = new AliITSvPPRasymmFMD("ITS","ITS PPR detailed version with asymmetric services");
+    if(iITS) {
+
+    //=================== 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");
+    //
+       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(kTRUE);   // 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->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
-*/
-      //
-      // 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);  
-    }  
+       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
+    //
+    //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);  
+    }
 
     if (iTPC)
     {
@@ -265,16 +209,25 @@ void ConfigCommon(Bool_t interactiveSetup)
         AliTPC *TPC = new AliTPCv2("TPC", "Default");
     }
 
-    if (iTOF)
-    {
+
+    if (iTOF) {
         //=================== TOF parameters ============================
-       AliTOF *TOF = new AliTOFv5T0("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 AliRICHv1("RICH", "normal RICH");
+        //=================== HMPID parameters ===========================
+        AliHMPID *HMPID = new AliHMPIDv2("HMPID", "normal HMPID");
+
     }
 
 
@@ -290,24 +243,18 @@ 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);
-        // Switch on TR
-        AliTRDsim *TRDsim = TRD->CreateTR();
     }
 
     if (iFMD)
     {
         //=================== FMD parameters ============================
-
-        AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
-    }
+       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 ===========================
@@ -324,34 +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 AliEMCALv2("EMCAL", "SHISH");
+        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 AliVZEROv6("VZERO", "normal VZERO");
+        //=================== ACORDE parameters ============================
+        AliVZERO *VZERO = new AliVZEROv7("VZERO", "normal VZERO");
     }
 
-  } // end (!isSetInteractively)
+     AliLog::Message(AliLog::kInfo, "End of Config", "Config.C", "Config.C", "Config()"," Config.C", __LINE__);
 
-  cout << "End of g4ConfigCommon.C" << endl;
 }
 
 Float_t EtaToTheta(Float_t arg){