AliGenerator* gener = 0x0;
if (proc == kStarlight) {
- gSystem->Load("libStarLight");
- gSystem->Load("libAliStarLight.so");
+ TString ap="../../build/lib/tgt_linux/";
+ gSystem->Load(ap+"libStarLight");
+ gSystem->Load(ap+"libAliStarLight");
AliGenStarLight* sl = new AliGenStarLight(1000*1000);
sl->SetParameter("CUT_ETA = 0 #Cut in pseudorapidity? (0 = no, 1 = yes)");
sl->SetParameter("ETA_MIN = -10 #Minimum pseudorapidity");
sl->SetParameter("ETA_MAX = 10 #Maximum pseudorapidity");
- sl->SetParameter("PROD_MODE = 2 #gg or gP switch (1 = 2-photon, 2 = coherent vector meson (narrow), 3 = coherent vector meson (wide), # 4 = incoherent vector meson, 5 = A+A DPMJet single, 6 = A+A DPMJet double, 7 = p+A DPMJet single, 8 = p+A Pythia single )");
- sl->SetParameter("PROD_PID = 113 #Channel of interest (not relevant for photonuclear processes)");
+ sl->SetParameter("PROD_MODE = 4 #gg or gP switch (1 = 2-photon, 2 = coherent vector meson (narrow), 3 = coherent vector meson (wide), # 4 = incoherent vector meson, 5 = A+A DPMJet single, 6 = A+A DPMJet double, 7 = p+A DPMJet single, 8 = p+A Pythia single )");
+ sl->SetParameter("PROD_PID = 913 #Channel of interest (not relevant for photonuclear processes)");
sl->SetParameter("RND_SEED = 34533 #Random number seed");
sl->SetParameter("BREAKUP_MODE = 5 #Controls the nuclear breakup");
sl->SetParameter("INTERFERENCE = 0 #Interference (0 = off, 1 = on)");
sl->SetParameter("INT_PT_MAX = 0.24 #Maximum pt considered, when interference is turned on");
sl->SetParameter("INT_PT_N_BINS = 120 #Number of pt bins when interference is turned on");
+ sl->SetRapidityMotherRange(-1., 1.);
+
sl->Init();
sl->GetTStarLight()->PrintInputs(std::cout);
gener->SetCutVertexZ(1.); // Truncate at 1 sigma
gener->SetVertexSmear(kPerEvent);
gener->SetTrackingFlag(1);
- gener->SetYRange(-0.9, 0.9);
+// gener->SetYRange(-0.9, 0.9);
gener->Init();
}