//
// AliRoot Configuration for running aliroot with Monte Carlo.
// Called from g4Config.C
+//
+// By I. Hrivnacova, IPN Orsay
-Float_t EtaToTheta(Float_t arg);
-static Int_t eventsPerRun = 100;
-enum PprGeo_t
-{
+enum PprGeo_t
+ {
kHoles, kNoHoles
-};
+ };
static PprGeo_t geo = kHoles;
-void ConfigCommon(Bool_t interactiveSetup)
+Float_t EtaToTheta(Float_t arg);
+void LoadPythia();
+static Int_t eventsPerRun = 50;
+
+void ConfigCommon(Bool_t setRootGeometry = kTRUE)
{
- // =============================
- // Root file
- // =============================
+ cout << "Running ConfigCommon.C ... " << endl;
+
+ // Load Pythia libraries
+ LoadPythia();
// Create the output file
AliRunLoader* rl = 0;
rl->SetCompressionLevel(2);
rl->SetNumberOfEventsPerFile(3);
gAlice->SetRunLoader(rl);
+
+ // Set Root geometry file
+ //
+ if ( setRootGeometry ) {
+ gAlice->SetRootGeometry();
+ gAlice->SetGeometryFromFile("geometry.root");
+ }
- // 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();
+ // Set the trigger configuration
+ gAlice->SetTriggerDescriptor("Pb-Pb");
+ cout<<"Trigger configuration is set to Pb-Pb"<<endl;
- // Activate this line if you want the vertex smearing to happen
- // track by track
+ // 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"));
+ }
- //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)
+ 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 ============================
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");
}
//=================== 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);
}
}
{
//=================== SHIL parameters ============================
- AliSHIL *SHIL = new AliSHILv2("SHIL", "Shielding Version 2");
+ AliSHIL *SHIL = new AliSHILv3("SHIL", "Shielding Version 3");
}
{
//=================== PIPE parameters ============================
- AliPIPE *PIPE = new AliPIPEv0("PIPE", "Beam Pipe");
+ AliPIPE *PIPE = new AliPIPEv3("PIPE", "Beam Pipe");
}
- if(iITS)
+ 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");
- 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->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 parameters ============================
+
+ AliITS *ITS = new AliITSv11Hybrid("ITS","ITS v11Hybrid");
+ }
if (iTPC)
{
AliTPC *TPC = new AliTPCv2("TPC", "Default");
}
- if (iTOF)
- {
+
+ if (iTOF) {
//=================== TOF parameters ============================
- AliTOF *TOF = new AliTOFv4T0("TOF", "normal TOF");
+ AliTOF *TOF = new AliTOFv6T0("TOF", "normal TOF");
}
- if (iRICH)
+
+ if (iHMPID)
{
- //=================== RICH parameters ===========================
- AliRICH *RICH = new AliRICHv1("RICH", "normal RICH");
+ //=================== HMPID parameters ===========================
+ AliHMPID *HMPID = new AliHMPIDv3("HMPID", "normal HMPID");
+
}
{
//=================== ZDC parameters ============================
- AliZDC *ZDC = new AliZDCv2("ZDC", "normal ZDC");
+ AliZDC *ZDC = new AliZDCv3("ZDC", "normal ZDC");
}
if (iTRD)
//=================== 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");
- }
+ AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
+ }
if (iMUON)
{
//=================== MUON parameters ===========================
-
+ // New MUONv1 version (geometry defined via builders)
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));
}
//=================== PHOS parameters ===========================
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", "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 AliVZEROv3("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;
- cout << "End of g4ConfigCommon.C" << endl;
}
Float_t EtaToTheta(Float_t arg){
return (180./TMath::Pi())*2.*atan(exp(-arg));
}
+
+void LoadPythia()
+{
+ // Load Pythia related libraries
+ gSystem->Load("liblhapdf.so"); // Parton density functions
+ gSystem->Load("libEGPythia6.so"); // TGenerator interface
+ gSystem->Load("libpythia6.so"); // Pythia
+ gSystem->Load("libAliPythia6.so"); // ALICE specific implementations
+}