From: mvl <marco.van.leeuwen@cern.ch>
Date: Wed, 23 Apr 2014 12:59:11 +0000 (+0200)
Subject: From Filip: Fill more histos for normalisation
X-Git-Url: http://git.uio.no/git/?a=commitdiff_plain;h=29ac2c9bebcba0cce52681f2842e0a93d2da5285;p=u%2Fmrichter%2FAliRoot.git

From Filip: Fill more histos for normalisation
---

diff --git a/PWGJE/AliAnalysisTaskJetCorePP.cxx b/PWGJE/AliAnalysisTaskJetCorePP.cxx
index c57c389c4fe..4171d2215fc 100644
--- a/PWGJE/AliAnalysisTaskJetCorePP.cxx
+++ b/PWGJE/AliAnalysisTaskJetCorePP.cxx
@@ -62,6 +62,7 @@
 #include "AliAODJet.h"
 #include "AliVVertex.h"
 #include "AliAnalysisTaskJetCorePP.h"
+#include "AliHeader.h" //KF//
 
 using std::cout;
 using std::endl;
@@ -907,20 +908,6 @@ void AliAnalysisTaskJetCorePP::UserExec(Option_t *)
 {
    //User Exec
    
-   if(fIsKine){ //Fill Xsection and number of trials
-      Float_t xsection = 0;
-      Float_t trials  = 0;
- 
-      AliMCEvent *mcEvent = MCEvent();
-      AliGenPythiaEventHeader *genPH = dynamic_cast<AliGenPythiaEventHeader*>(mcEvent->GenEventHeader());
-      if(genPH){
-         xsection = genPH->GetXsection();
-         trials  = genPH->Trials();
-      }
-      fh1Xsec->Fill("<#sigma>",xsection);
-      fh1Trials->Fill("#sum{ntrials}",trials);
-   }//KF//
-
 
    //Event loop
    Double_t eventW  = 1.0;
@@ -966,8 +953,25 @@ void AliAnalysisTaskJetCorePP::UserExec(Option_t *)
          PostData(1, fOutputList);
          return;
       }
+ 
+      Float_t xsection = 0;
+      Float_t trials  = 0;
+
+      AliGenPythiaEventHeader *genPH =
+         dynamic_cast<AliGenPythiaEventHeader*> (fMcEvent->GenEventHeader()); 
+      if(genPH){
+         xsection = genPH->GetXsection();
+         trials   = genPH->Trials();
+         ptHard   = genPH->GetPtHard();
+      }
+      fh1Xsec->Fill("<#sigma>",xsection);
+      fh1Trials->Fill("#sum{ntrials}",trials);
+      fh1PtHard->Fill(ptHard,eventW);
+      fh1PtHardNoW->Fill(ptHard,1);
+      fh1PtHardTrials->Fill(ptHard,trials);
    }
 
+
    fESD = dynamic_cast<AliESDEvent*>(InputEvent());
    if(!fESD){
       if(fDebug>1) AliError("ESD not available");