@date Mon Mar 27 12:50:29 2006
@brief Simulation configuration script
*/
-#if !defined(__CINT__) || defined(__MAKECINT__)
-#include <Riostream.h>
-#include <TPDGCode.h>
-#include <TSystem.h>
-#include <TVirtualMC.h>
-#include <TGeant3.h>
-#include "STEER/AliRunLoader.h"
-#include "STEER/AliRun.h"
-#include "STEER/AliConfig.h"
-#include "STEER/AliGenerator.h"
-#include "PYTHIA6/AliDecayerPythia.h"
-#include "EVGEN/AliGenHIJINGpara.h"
-#include "THijing/AliGenHijing.h"
-#include "EVGEN/AliGenCocktail.h"
-#include "EVGEN/AliGenSlowNucleons.h"
-#include "EVGEN/AliSlowNucleonModelExp.h"
-#include "EVGEN/AliGenParam.h"
-#include "EVGEN/AliGenMUONlib.h"
-#include "EVGEN/AliGenMUONCocktail.h"
-#include "PYTHIA6/AliGenPythia.h"
-#include "STEER/AliMagFMaps.h"
-#include "STRUCT/AliBODY.h"
-#include "STRUCT/AliMAG.h"
-#include "STRUCT/AliABSOv0.h"
-#include "STRUCT/AliDIPOv2.h"
-#include "STRUCT/AliHALL.h"
-#include "STRUCT/AliFRAMEv2.h"
-#include "STRUCT/AliSHILv2.h"
-#include "STRUCT/AliPIPEv0.h"
-#include "ITS/AliITSvPPRasymmFMD.h"
-#include "TPC/AliTPCv2.h"
-#include "TOF/AliTOFv4T0.h"
-#include "HMPID/AliHMPIDv1.h"
-#include "ZDC/AliZDCv2.h"
-#include "TRD/AliTRDv1.h"
-#include "FMD/AliFMDv1.h"
-#include "MUON/AliMUONv1.h"
-#include "MUON/AliMUONSt1GeometryBuilderV2.h"
-#include "MUON/AliMUONSt2GeometryBuilder.h"
-#include "MUON/AliMUONSlatGeometryBuilder.h"
-#include "MUON/AliMUONTriggerGeometryBuilder.h"
-#include "PHOS/AliPHOSv1.h"
-#include "PMD/AliPMDv1.h"
-#include "T0/AliT0v1.h"
-#include "EMCAL/AliEMCALv1.h"
-#include "ACORDE/AliACORDEv0.h"
-#include "VZERO/AliVZEROv2.h"
-#endif
//____________________________________________________________________
//
//____________________________________________________________________
// Functions
Float_t EtaToTheta(Float_t eta);
-Eg_t LookupEG(const Char_t* name);
+EG_t LookupEG(const Char_t* name);
AliGenerator* GeneratorFactory(EG_t eg, Rad_t rad, TString& comment);
AliGenHijing* HijingStandard();
void ProcessEnvironmentVars(EG_t& eg, Int_t& seed);
// This part for configuration
// EG_t eg = test50;
// EG_t eg = kParam_fmd;
- EG_t eg = kParam_2000; // kPythia;
+ // EG_t eg = kFMDFlat; // kParam_2000; // kPythia;
// EG_t eg = kFMDFlat;
+ EG_t eg = kPythia6;
+ // EG_t eg = kFMD2Flat;
Geo_t geo = kNoHoles;
Rad_t rad = kGluonRadiation;
Mag_t mag = k5kG;
//
// Libraries needed by GEANT 3.21
//
+ gSystem->Load("$ALICE_ROOT/lib/tgt_$ALICE_TARGET/libpythia6.so");
gSystem->Load("libgeant321");
//
//
// Libraries needed by GEANT 3.21
//
+ gSystem->Load("$ALICE_ROOT/lib/tgt_$ALICE_TARGET/liblhapdf.so");
+ gSystem->Load("$ALICE_ROOT/lib/tgt_$ALICE_TARGET/libpythia6.so");
+ // gSystem->Load("libEGPythia6.so"); //<- For non-debian (sigh!)
gSystem->Load("EGPythia6.so");
gSystem->Load("libgeant321");
gMC->SetProcess("MUNU",1);
gMC->SetProcess("CKOV",1);
gMC->SetProcess("HADR",1);
- gMC->SetProcess("LOSS",2); // 0:none 1,3:dray 2:nodray 4:nofluct (def:2)
+ gMC->SetProcess("LOSS",1); // 0:none 1,3:dray 2:nodray 4:nofluct (def:2)
gMC->SetProcess("MULS",1);
gMC->SetProcess("RAYL",1);
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 useACORDE = kFALSE;
+ Bool_t useACORDE= kFALSE;
Bool_t useDIPO = kFALSE;
Bool_t useFMD = kTRUE;
Bool_t useFRAME = kFALSE;
Bool_t useMAG = kFALSE;
Bool_t useMUON = kFALSE;
Bool_t usePHOS = kFALSE;
- Bool_t usePIPE = kFALSE;
+ Bool_t usePIPE = kTRUE;
Bool_t usePMD = kFALSE;
- Bool_t useHMPID = kFALSE;
+ Bool_t useHMPID = kFALSE;
Bool_t useSHIL = kFALSE;
- Bool_t useT0 = kFALSE;
+ Bool_t useT0 = kFALSE;
Bool_t useTOF = kFALSE;
Bool_t useTPC = kFALSE;
Bool_t useTRD = kFALSE;
AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
- if (useMAG) {
- // =================== 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 ===========================
- AliABSO *ABSO = new AliABSOv0("ABSO", "Muon Absorber");
- }
-
- if (useDIPO) {
- // =================== DIPO parameters ===========================
- AliDIPO *DIPO = new AliDIPOv2("DIPO", "Dipole version 2");
- }
-
- if (useHALL) {
- // =================== HALL parameters ===========================
- AliHALL *HALL = new AliHALL("HALL", "Alice Hall");
- }
-
-
+ if (useMAG) AliMAG *MAG = new AliMAG("MAG", "Magnet");
+ if (useABSO) AliABSO *ABSO = new AliABSOv3("ABSO", "Muon Absorber");
+ if (useDIPO) AliDIPO *DIPO = new AliDIPOv3("DIPO", "Dipole version 3");
+ if (useHALL) AliHALL *HALL = new AliHALLv3("HALL", "Alice Hall");
if (useFRAME) {
- // ================== FRAME parameters ===========================
AliFRAMEv2 *FRAME = new AliFRAMEv2("FRAME", "Space Frame");
switch (geo) {
case kHoles: FRAME->SetHoles(1); break;
default: FRAME->SetHoles(0); break;
}
}
-
- if (useSHIL) {
- // ================== SHIL parameters ============================
- AliSHIL *SHIL = new AliSHILv2("SHIL", "Shielding Version 2");
- }
-
-
- if (usePIPE) {
- // ================== PIPE parameters ============================
- AliPIPE *PIPE = new AliPIPEv0("PIPE", "Beam Pipe");
- }
-
- if (useITS) {
- // =================== 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");
- // don't touch this parameter if you're not an ITS developer
- 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");
- // detector thickness on layer 1 must be in the range [100,300]
- ITS->SetThicknessDet1(200.);
- // detector thickness on layer 2 must be in the range [100,300]
- ITS->SetThicknessDet2(200.);
- // chip thickness on layer 1 must be in the range [150,300]
- ITS->SetThicknessChip1(200.);
- // chip thickness on layer 2 must be in the range [150,300]
- ITS->SetThicknessChip2(200.);
- // 1 --> rails in ; 0 --> rails out
- ITS->SetRails(0);
- // 1 --> water ; 0 --> freon
- ITS->SetCoolingFluid(1);
-
- // 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");
- // 1 --> rails in ; 0 --> rails out
- // ITS->SetRails(0);
- // 0 --> Copper ; 1 --> Aluminum ; 2 --> Carbon
- // ITS->SetSupportMaterial(0);
- //
- // AliITS *ITS =
- // new AliITSvPPRcoarsesymm("ITS","New ITS PPR coarse version "
- // "with symmetric services");
- // 1 --> rails in ; 0 --> rails out
- // ITS->SetRails(0);
- // 0 --> Copper ; 1 --> Aluminum ; 2 --> Carbon
- // ITS->SetSupportMaterial(0);
- //
- // 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 (useTPC) {
- // =================== 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();
- AliTPC *TPC = new AliTPCv2("TPC", "Default");
- }
-
- if (useTOF) {
- // ================== TOF parameters =============================
- AliTOF *TOF = new AliTOFv4T0("TOF", "normal TOF");
- }
-
- if (useHMPID) {
- // ================== HMPID parameters ============================
- AliHMPID *HMPID = new AliHMPIDv1("HMPID", "normal HMPID");
-
- }
-
- if (useZDC) {
- // ================== ZDC parameters =============================
- AliZDC *ZDC = new AliZDCv2("ZDC", "normal ZDC");
- }
-
- if (useTRD) {
- // ================== 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 HMPID
- TRD->SetHMPIDhole();
- }
- // Switch on TR
- AliTRDsim *TRDsim = TRD->CreateTR();
- }
-
- if (useFMD) {
- // =================== FMD parameters ============================
- AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
- // FMD->UseDetailed(kFALSE);
- // FMD->UseAssembly();
- // FMD->UseOld();
- }
-
- if (useMUON) {
- // =================== 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));
- }
-
- if (usePHOS) {
- // =================== PHOS parameters ===========================
- AliPHOS *PHOS = new AliPHOSv1("PHOS", "IHEP");
- }
-
- if (usePMD) {
- // =================== PMD parameters ============================
- AliPMD *PMD = new AliPMDv1("PMD", "normal PMD");
- }
-
- if (useT0) {
- // =================== T0 parameters ==========================
- AliT0 *T0 = new AliT0v1("T0", "T0 Detector");
- }
-
- if (useEMCAL) {
- // =================== EMCAL parameters ==========================
- AliEMCAL *EMCAL = new AliEMCALv1("EMCAL", "EMCAL_55_25");
- }
-
- if (useACORDE) {
- // =================== ACORDE parameters ============================
- AliACORDE *ACORDE = new AliACORDEv0("ACORDE", "normal ACORDE");
- }
-
- if (useVZERO) {
- // =================== V0 parameters =============================
- AliVZERO *VZERO = new AliVZEROv3("VZERO", "normal VZERO");
- }
+ if (useSHIL) AliSHIL *SHIL = new AliSHILv3("SHIL", "Shielding v3");
+ if (usePIPE) AliPIPE *PIPE = new AliPIPEv3("PIPE", "Beam Pipe");
+ if (useITS) AliITS *ITS = new AliITSv11Hybrid("ITS","ITS v11Hybrid");
+ if (useTPC) AliTPC *TPC = new AliTPCv2("TPC", "Default");
+ if (useTOF) AliTOF *TOF = new AliTOFv6T0("TOF", "normal TOF");
+ if (useHMPID) AliHMPID *HMPID = new AliHMPIDv1("HMPID", "normal HMPID");
+ if (useZDC) AliZDC *ZDC = new AliZDCv3("ZDC", "normal ZDC");
+ if (useTRD) AliTRD *TRD = new AliTRDv1("TRD", "TRD slow simulator");
+ if (useFMD) AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
+ if (useMUON) AliMUON *MUON = new AliMUONv1("MUON", "default");
+ if (usePHOS) AliPHOS *PHOS = new AliPHOSv1("PHOS", "IHEP");
+ if (usePMD) AliPMD *PMD = new AliPMDv1("PMD", "normal PMD");
+ if (useT0) AliT0 *T0 = new AliT0v1("T0", "T0 Detector");
+ if (useEMCAL) AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETE");
+ if (useACORDE) AliACORDE *ACORDE = new AliACORDEv1("ACORDE", "normal ACORDE");
+ if (useVZERO) AliVZERO *VZERO = new AliVZEROv7("VZERO", "normal VZERO");
}
//____________________________________________________________________
Int_t isw = 3;
if (rad == kNoGluonRadiation) isw = 0;
+ // Possibly load libAliPythia6
+ switch (eg) {
+ case kPythia6:
+ case kPythia6Jets20_24:
+ case kPythia6Jets24_29:
+ case kPythia6Jets29_35:
+ case kPythia6Jets35_42:
+ case kPythia6Jets42_50:
+ case kPythia6Jets50_60:
+ case kPythia6Jets60_72:
+ case kPythia6Jets72_86:
+ case kPythia6Jets86_104:
+ case kPythia6Jets104_125:
+ case kPythia6Jets125_150:
+ case kPythia6Jets150_180:
+ gSystem->Load("liblhapdf.so");
+ // gSystem->Load("/usr/lib/libpythia.so");
+ // gSystem->ListLibraries();
+ gSystem->Load("EGPythia6.so");
+ gSystem->Load("libAliPythia6");
+ break;
+ }
AliGenerator * gGener = 0;
switch (eg) {
case kFMD2Flat:
{
comment = comment.Append(" Flat in FMD2 range");
- AliGenBox* gener = new AliGenBox(2000);
+ AliGenBox* gener = new AliGenBox(10);
gener->SetPart(kPiPlus);
gener->SetMomentumRange(3,4);
gener->SetPhiRange(0, 360);