5 gSystem->AddDynamicPath("../../../objdir/lib/tgt_linux/");
6 gSystem->Load("libStarLight");
7 gSystem->Load("libAliStarLight.so");
9 TStarLight* sl = new TStarLight("starlight generator", "title", "");
11 sl->SetParameter("BEAM_1_Z = 82 #Z of projectile");
12 sl->SetParameter("BEAM_1_A = 208 #A of projectile");
13 sl->SetParameter("BEAM_2_Z = 82 #Z of target");
14 sl->SetParameter("BEAM_2_A = 208 #A of target");
15 sl->SetParameter("BEAM_1_GAMMA = 1470 #Gamma of the colliding ions");
16 sl->SetParameter("BEAM_2_GAMMA = 1470 #Gamma of the colliding ions");
17 sl->SetParameter("W_MAX = 15 #Max value of w");
18 sl->SetParameter("W_MIN = -1 #Min value of w");
19 sl->SetParameter("W_N_BINS = 40 #Bins i w");
20 sl->SetParameter("RAP_MAX = 8. #max y");
21 sl->SetParameter("RAP_N_BINS = 80 #Bins i y");
22 sl->SetParameter("CUT_PT = 0 #Cut in pT? 0 = (no, 1 = yes)");
23 sl->SetParameter("PT_MIN = 1.0 #Minimum pT in GeV");
24 sl->SetParameter("PT_MAX = 3.0 #Maximum pT in GeV");
25 sl->SetParameter("CUT_ETA = 0 #Cut in pseudorapidity? (0 = no, 1 = yes)");
26 sl->SetParameter("ETA_MIN = -10 #Minimum pseudorapidity");
27 sl->SetParameter("ETA_MAX = 10 #Maximum pseudorapidity");
28 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 )");
30 sl->SetParameter("N_EVENTS = 1000 #Number of events");
31 sl->SetParameter("PROD_PID = 113 #Channel of interest (not relevant for photonuclear processes)");
32 sl->SetParameter("RND_SEED = 34533 #Random number seed");
33 sl->SetParameter("OUTPUT_FORMAT = 2 #Form of the output");
34 sl->SetParameter("BREAKUP_MODE = 5 #Controls the nuclear breakup");
35 sl->SetParameter("INTERFERENCE = 0 #Interference (0 = off, 1 = on)");
36 sl->SetParameter("IF_STRENGTH = 1. #% of intefernce (0.0 - 0.1)");
37 sl->SetParameter("COHERENT = 1 #Coherent=1,Incoherent=0");
38 sl->SetParameter("INCO_FACTOR = 1. #percentage of incoherence");
39 sl->SetParameter("BFORD = 9.5 #");
40 sl->SetParameter("INT_PT_MAX = 0.24 #Maximum pt considered, when interference is turned on");
41 sl->SetParameter("INT_PT_N_BINS = 120 #Number of pt bins when interference is turned on");
44 sl->PrintInputs(std::cout);
46 sl->GetParameter("INT_PT_MAX");
48 TClonesArray tca("TParticle", 1000);
50 for (Int_t counter(0); counter<1000; ++counter) {
51 Printf("--------------------------------------------------------------------------------");
53 sl->ImportParticles(&tca, "ALL");
54 for (Int_t i=0; i<tca.GetEntries(); ++i)