STARLIGHT/starlight/objdir/testsl.C

   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 }