git://git.uio.no / u / mrichter / AliRoot.git / blame
| Commit | Line | Data |
|---|---|---|
| da32329d AM |
1 | // -*- C++ -*- |
| 2 | // $Id$ | |
| 3 | ||
| 4 | void testsl() { | |
| 5 | gSystem->AddDynamicPath("../objdir"); | |
| 6 | gSystem->Load("libStarLight"); | |
| 7 | gSystem->Load("libTStarLight.so"); | |
| 8 | ||
| 9 | TStarLight* sl = new TStarLight("starlight generator", "title", ""); | |
| 10 | ||
| 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 = 12.0 #Max value of w"); | |
| 18 | sl->SetParameter("W_MIN = 2.0 #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 )"); | |
| 29 | // is N_EVENTS valid | |
| 30 | sl->SetParameter("N_EVENTS = 20 #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"); | |
| 42 | ||
| 43 | sl->InitStarLight(); | |
| 44 | sl->PrintInputs(std::cout); | |
| 45 | ||
| 46 | sl->GetParameter("INT_PT_MAX"); | |
| 47 | ||
| 48 | TClonesArray tca("TParticle", 1000); | |
| 49 | ||
| 50 | for (Int_t counter(0); counter<20; ++counter) { | |
| 51 | Printf("--------------------------------------------------------------------------------"); | |
| 52 | sl->GenerateEvent(); | |
| 53 | sl->ImportParticles(&tca, "ALL"); | |
| 54 | for (Int_t i=0; i<tca.GetEntries(); ++i) | |
| 55 | tca.At(i)->Print(); | |
| 56 | } | |
| 57 | } |