Updated ctiteria, renamed files, modified printout (Jacek, Marian)

[u/mrichter/AliRoot.git] / EVGEN / TestPrimaries.C

diff --git a/EVGEN/TestPrimaries.C b/EVGEN/TestPrimaries.C
index 51ac7003fbf40fdfee72b4595c568df84bed8bdc..8c46fa73335b9d91e10fb1a16b8b3c2d979b6aad 100644 (file)
--- a/EVGEN/TestPrimaries.C
+++ b/EVGEN/TestPrimaries.C
@@ -13,12 +13,21 @@ void TestPrimaries(Int_t evNumber1=0, Int_t evNumber2=0)
 //  Create some histograms
 
     TH1F *thetaH   =  new TH1F("thetaH","Theta distribution",180,0,180);
+    thetaH->Sumw2();
     TH1F *phiH     =  new TH1F("phiH","Phi distribution"  ,180,-180,180);
+    phiH->Sumw2();
+    TH1F *phiH1    =  new TH1F("phiH1","Phi distribution"  ,180,-180,180);
+    phiH1->Sumw2();
     TH1F *etaH     =  new TH1F("etaH","Pseudorapidity",120,-12,12);
+    etaH->Sumw2();
     TH1F *yH       =  new TH1F("yH","Rapidity distribution",120,-12,12);
+    yH->Sumw2();
     TH1F *eH       =  new TH1F("eH","Energy distribution",100,0,100);
+    eH->Sumw2();
     TH1F *eetaH    =  new TH1F("eetaH","Pseudorapidity",120,0,12);
+    eetaH->Sumw2();
     TH1F *ptH      =  new TH1F("ptH","Pt distribution",150,0,15);
+    ptH->Sumw2();
 //
 //   Loop over events 
 //
@@ -55,7 +64,10 @@ void TestPrimaries(Int_t evNumber1=0, Int_t evNumber2=0)
            if (DataBase->GetParticle(mpart)->Charge() == 0) continue;
            Float_t wgt = 1./(Float_t ((evNumber2-evNumber1)+1.));
            thetaH->Fill(theta*180./TMath::Pi(),wgt);
-           phiH->Fill(phi*180./TMath::Pi(),wgt);
+           if (pT<2)
+             phiH->Fill(phi*180./TMath::Pi(),wgt);
+           else
+             phiH1->Fill(phi*180./TMath::Pi(),wgt);
            etaH->Fill(eta,5.*wgt);    
            eetaH->Fill(eta,E);    
            yH->Fill(y,5.*wgt);      
@@ -75,14 +87,16 @@ void TestPrimaries(Int_t evNumber1=0, Int_t evNumber2=0)
     TCanvas *c2 = new TCanvas("c2","Canvas 1",400,10,600,700);
     c2->Divide(2,2);
     c2->cd(1); phiH->Draw();
-    c2->cd(2); thetaH->Draw();
-    c2->cd(3); eetaH->Draw();
+    c2->cd(2); phiH1->Draw();
+    c2->cd(3); thetaH->Draw();
+    c2->cd(4); eetaH->Draw();
 
     TFile* f = new TFile("kineH.root", "update");
     ptH->Write();
     etaH->Write();
     yH->Write();
     phiH->Write();
+    phiH1->Write();
     f->Close();