#if !defined(__CINT__) || defined(__MAKECINT__)
#include <Riostream.h>
+#include <TPDGCode.h>
#include <TRandom.h>
#include <TSystem.h>
#include <TVirtualMC.h>
#include "PYTHIA6/AliDecayerPythia.h"
#include "EVGEN/AliGenCocktail.h"
#include "EVGEN/AliGenHIJINGpara.h"
-#include "STEER/AliMagFMaps.h"
+#include "STEER/AliMagF.h"
#include "STRUCT/AliBODY.h"
#include "STRUCT/AliMAG.h"
#include "STRUCT/AliABSOv3.h"
#include "STRUCT/AliFRAMEv2.h"
#include "STRUCT/AliSHILv3.h"
#include "STRUCT/AliPIPEv3.h"
-#include "ITS/AliITSvPPRasymmFMD.h"
+#include "ITS/AliITSv11Hybrid.h"
#include "TPC/AliTPCv2.h"
#include "TOF/AliTOFv6T0.h"
-#include "HMPID/AliHMPIDv2.h"
-#include "ZDC/AliZDCv2.h"
+#include "HMPID/AliHMPIDv3.h"
+#include "ZDC/AliZDCv3.h"
#include "TRD/AliTRDv1.h"
#include "FMD/AliFMDv1.h"
#include "MUON/AliMUONv1.h"
#include "PMD/AliPMDv1.h"
#include "T0/AliT0v1.h"
#include "EMCAL/AliEMCALv2.h"
-#include "ACORDE/AliACORDEv0.h"
+#include "ACORDE/AliACORDEv1.h"
#include "VZERO/AliVZEROv7.h"
#endif
Float_t EtaToTheta(Float_t arg);
+void LoadPythia();
+
void Config()
{
// Theta range given through pseudorapidity limits 22/6/2001
// Set Random Number seed
- gRandom->SetSeed(123456); // Set 0 to use the currecnt time
+ //gRandom->SetSeed(123456); // Set 0 to use the current time
+
AliLog::Message(AliLog::kInfo, Form("Seed for random number generation = %d",gRandom->GetSeed()), "Config.C", "Config.C", "Config()","Config.C", __LINE__);
-
- // libraries required by geant321
+ // Load Pythia libraries
+ LoadPythia();
+ // Libraries required by geant321
#if defined(__CINT__)
gSystem->Load("libgeant321");
#endif
new TGeant3TGeo("C++ Interface to Geant3");
- if(!AliCDBManager::Instance()->IsDefaultStorageSet()){
- AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT");
- AliCDBManager::Instance()->SetRun(0);
- }
-
AliRunLoader* rl=0x0;
AliLog::Message(AliLog::kInfo, "Creating Run Loader", "Config.C", "Config.C", "Config()"," Config.C", __LINE__);
rl->SetCompressionLevel(2);
rl->SetNumberOfEventsPerFile(3);
gAlice->SetRunLoader(rl);
+
// gAlice->SetGeometryFromFile("geometry.root");
- // gAlice->SetGeometryFromCDB();
+
+ // Uncomment if you want to load geometry from OCDB! >>>>
+/*
+ if(!AliCDBManager::Instance()->IsDefaultStorageSet()){
+ cout << "#####################################################" << endl;
+ cout << "# #" << endl;
+ cout << "# WARNING: CDB DEFAULT STORAGE NOT SET !!! #" << endl;
+ cout << "# SETTING IT TO local://$ALICE_ROOT/OCDB !!! #" << endl;
+ cout << "# #" << endl;
+ cout << "#####################################################" << endl;
+
+ AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
+ }
+
+ if(AliCDBManager::Instance()->GetRun() < 0){
+ cout << "#####################################################" << endl;
+ cout << "# #" << endl;
+ cout << "# WARNING: RUN NUMBER NOT SET !!! #" << endl;
+ cout << "# SETTING IT TO 0 !!! #" << endl;
+ cout << "# #" << endl;
+ cout << "#####################################################" << endl;
+
+ AliCDBManager::Instance()->SetRun(0);
+ }
+ gAlice->SetGeometryFromCDB();
+*/
+ // Uncomment if you want to load geometry from OCDB! <<<<
// Set the trigger configuration
- gAlice->SetTriggerDescriptor("Pb-Pb");
+ AliSimulation::Instance()->SetTriggerConfig("Pb-Pb");
cout<<"Trigger configuration is set to Pb-Pb"<<endl;
//
gener->AddGenerator(hijingparam,"HIJING PARAM",1);
// AliGenBox *genbox = new AliGenBox(nParticles);
-// genbox->SetPart(22);
+// genbox->SetPart(kGamma);
// genbox->SetPtRange(0.3, 10.00);
// gener->AddGenerator(genbox,"GENBOX GAMMA for PHOS",1);
gener->Init();
//
//gener->SetVertexSmear(perTrack);
// 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;
+ 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 = 1;
+ Int_t iVZERO = 1;
rl->CdGAFile();
//=================== Alice BODY parameters =============================
AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
//=================== FRAME parameters ============================
AliFRAMEv2 *FRAME = new AliFRAMEv2("FRAME", "Space Frame");
+ FRAME->SetHoles(1);
}
if (iSHIL)
{
//=================== ITS parameters ============================
- AliITSvPPRasymmFMD *ITS = new AliITSvPPRasymmFMD("ITS","ITS PPR detailed version with asymmetric services");
+ AliITS *ITS = new AliITSv11Hybrid("ITS","ITS v11Hybrid");
}
if (iTPC)
if (iHMPID)
{
//=================== HMPID parameters ===========================
- AliHMPID *HMPID = new AliHMPIDv2("HMPID", "normal HMPID");
+ 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)
if (iEMCAL)
{
//=================== EMCAL parameters ============================
- AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "SHISH_77_TRD1_2X2_FINAL_110DEG");
+ AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETEV1");
}
if (iACORDE)
{
//=================== ACORDE parameters ============================
- AliACORDE *ACORDE = new AliACORDEv0("ACORDE", "normal ACORDE");
+ AliACORDE *ACORDE = new AliACORDEv1("ACORDE", "normal ACORDE");
}
if (iVZERO)
{
- //=================== ACORDE parameters ============================
+ //=================== VZERO parameters ============================
AliVZERO *VZERO = new AliVZEROv7("VZERO", "normal VZERO");
}
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
+}