updated STARLIGHT to r191 of http://starlight.hepforge.org/svn/trunk

[u/mrichter/AliRoot.git] / STARLIGHT / starlight / TStarLight / testsl.C

diff --git a/STARLIGHT/starlight/TStarLight/testsl.C b/STARLIGHT/starlight/TStarLight/testsl.C
index 1e70384950dcd6e28c484600f56337553cd8421e..5f4c6f71761adec750a5931548c519776be94725 100644 (file)
--- a/STARLIGHT/starlight/TStarLight/testsl.C
+++ b/STARLIGHT/starlight/TStarLight/testsl.C
@@ -37,7 +37,7 @@ void testsl() {
   sl->SetParameter("ETA_MIN      =  -10    #Minimum pseudorapidity");
   sl->SetParameter("ETA_MAX      =   10    #Maximum pseudorapidity");
   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     =  113    #Channel of interest (not relevant for photonuclear processes)");
+  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)");
@@ -52,25 +52,27 @@ void testsl() {
   TClonesArray tca("TParticle", 100);
 
   TLorentzVector v[2], vSum;
-  TH1* hM  = new TH1D("hM",  "STARLIGHT;M#(){#pi^{+}#pi^{-}}",     100, 0., 2.);
-  TH1* hPt = new TH1D("hPt", "STARLIGHT;P_{T}#(){#pi^{+}#pi^{-}}", 100, 0., 1.);
-  TH1* hY  = new TH1D("hY",  "STARLIGHT;Y#(){#pi^{+}#pi^{-}}",     160,-8., 8.);
+  TH1* hM  = new TH1D("hM",  "STARLIGHT;M#(){#pi^{+}#pi^{-}}",     80, 0., 2.);
+  TH1* hPt = new TH1D("hPt", "STARLIGHT;P_{T}#(){#pi^{+}#pi^{-}}", 80, 0., 2.);
+  TH1* hY  = new TH1D("hY",  "STARLIGHT;Y#(){#pi^{+}#pi^{-}}",    100,-10., 10.);
 
   std::ofstream ofs("rho0.txt");
   TParticle *p;
-  for (Int_t counter(0); counter<1000;) {
+  for (Int_t counter(0); counter<20000;) {
     sl->GenerateEvent();    
     sl->BoostEvent();    
     sl->ImportParticles(&tca, "ALL");
-    Bool_t isOK = kTRUE;
-    for (Int_t i=0; i<tca.GetEntries() && isOK; ++i) {
+    Bool_t genOK = kTRUE;
+    TLorentzVector vSum;
+    for (Int_t i=0; i<tca.GetEntries() && genOK; ++i) {
       p = (TParticle*)tca.At(i);
       p->Momentum(v[i]);
-//       isOK = TMath::Abs(v[i].Rapidity()) <= 0.9;
+      vSum += v[i];
+//       genOK = TMath::Abs(v[i].Rapidity()) <= 1.5;
     }
     tca.Clear();
-    if (!isOK) continue;
-    Printf("counter, %d", counter, isOK);
+    if (!genOK) continue;
+    Printf("counter, %d", counter, genOK);
     ++counter;
     vSum = v[0] + v[1];
     ofs << std::fixed << std::setprecision(4)