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