// gAlice->SetConfigFunction("Config('$HOME','box');");
// april 3rd: added L3 magnet
+void LoadPythia();
+
void Config(char directory[100]="", char option[6]="trg1mu")
{
static Int_t sseed = 0; // Set 0 to use the current time
//=====================================================================
- // Libraries required by geant321
+ LoadPythia();
gSystem->Load("libgeant321.so");
new TGeant3TGeo("C++ Interface to Geant3");
+
//=======================================================================
if(!AliCDBManager::Instance()->IsDefaultStorageSet()){
- AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT");
+ AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
AliCDBManager::Instance()->SetRun(0);
}
// External Radius (in cm)
// 183., 245., 395., 560., 563., 850., 900.
//=======================================================================
+ AliGenPythia* PythiaForMUONCocktail(Decay_t dt)
+ {
+ AliGenPythia *pythia = new AliGenPythia(1);
+ pythia->SetProcess(kPyMbMSEL1);
+ pythia->SetStrucFunc(kCTEQ5L);
+ pythia->SetEnergyCMS(14000.);
+ pythia->SetForceDecay(dt);
+ pythia->SetPtRange(0.,100.);
+ pythia->SetYRange(-8.,8.);
+ pythia->SetPhiRange(0.,360.);
+ pythia->SetPtHard(2.76,-1.0);
+ pythia->SwitchHFOff();
+ return pythia;
+ }
if (!strcmp(option,"trg2mu")) {
AliGenMUONCocktailpp *gener = new AliGenMUONCocktailpp();
gener->SetOrigin(0.,0.,0.);
gener->SetSigma(0.,0.,5.);
gener->SetVertexSmear(kPerEvent);
+ Decay_t dt = gener->GetDecayModePythia(dt);
+ AliGenPythia* pythia = PythiaForMUONCocktail(dt);
+ pythia->Init();
+ gener->AddGenerator(pythia,"Pythia",1);
gener->Init();
}
gener->SetOrigin(0.,0.,0.);
gener->SetSigma(0.,0.,5.);
gener->SetVertexSmear(kPerEvent);
+ Decay_t dt = gener->GetDecayModePythia();
+ AliGenPythia* pythia = PythiaForMUONCocktail(dt);
+ pythia->Init();
+ gener->AddGenerator(pythia,"Pythia",1);
gener->Init();
}
//=============================================================
// Field (L3 0.5 T) outside dimuon spectrometer
- AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., AliMagFMaps::k5kG);
+ AliMagF* field = new AliMagF("Maps","Maps",2,1.,1., 10.,AliMagF::k5kG);
field->SetL3ConstField(0); // Using const. field in the barrel
- gAlice->SetField(field);
+ TGeoGlobalMagField::Instance()->SetField(field);
Int_t iITS = 1;
Int_t iZDC = 1;
AliSHIL *SHIL = new AliSHILv3("SHIL", "Shielding Version 3");
//=================== ITS parameters =============================
if(iITS) {
- AliITSvPPRasymmFMD *ITS = new AliITSvPPRasymmFMD("ITS","ITS PPR detailed version with asymmetric services");
+ AliITS *ITS = new AliITSv11Hybrid("ITS","ITS v11Hybrid");
}
//=================== ZDC parameters =============================
if(iZDC) {
- AliZDC *ZDC = new AliZDCv2("ZDC", "normal ZDC");
+ AliZDC *ZDC = new AliZDCv3("ZDC", "normal ZDC");
}
//=================== FMD parameters =============================
if(iFMD) {
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
+}
+
+