fill average trials in a separeted histo

diff --git a/PWGJE/AliAnalysisTaskJetServices.cxx b/PWGJE/AliAnalysisTaskJetServices.cxx
index c0c0bd91d75e6616b9e67d5b3adea9578135fbd8..1d00af66124c0a107d2cd7cb6ed136e7dbcf8091 100644 (file)
--- a/PWGJE/AliAnalysisTaskJetServices.cxx
+++ b/PWGJE/AliAnalysisTaskJetServices.cxx
@@ -108,6 +108,7 @@ AliAnalysisTaskJetServices::AliAnalysisTaskJetServices():
   fTriggerName(0x0),
   fh1Xsec(0x0),
   fh1Trials(0x0),
+  fh1AvgTrials(0x0),
   fh1PtHard(0x0),
   fh1PtHardTrials(0x0),
   fh1SelectionInfoESD(0x0),
@@ -181,6 +182,7 @@ AliAnalysisTaskJetServices::AliAnalysisTaskJetServices(const char* name):
   fTriggerName(0x0),
   fh1Xsec(0x0),
   fh1Trials(0x0),
+  fh1AvgTrials(0x0),
   fh1PtHard(0x0),
   fh1PtHardTrials(0x0),
   fh1SelectionInfoESD(0x0),
@@ -250,6 +252,7 @@ Bool_t AliAnalysisTaskJetServices::Notify()
     // construct a poor man average trials 
     Float_t nEntries = (Float_t)tree->GetTree()->GetEntries();
     if(ftrials>=nEntries && nEntries>0.)fAvgTrials = ftrials/nEntries;
+    fh1Trials->Fill("#sum{avg ntrials}",ftrials); 
   }  
   return kTRUE;
 }
@@ -279,6 +282,10 @@ void AliAnalysisTaskJetServices::UserCreateOutputObjects()
   fh1Trials->GetXaxis()->SetBinLabel(1,"#sum{ntrials}");
   fHistList->Add(fh1Trials);
 
+  fh1AvgTrials = new TH1F("fh1AvgTrials","avg trials root file",1,0,1);
+  fh1AvgTrials->GetXaxis()->SetBinLabel(1,"#sum{avg ntrials}");
+  fHistList->Add(fh1AvgTrials);
+
   const Int_t nBinPt = 125;
   Double_t binLimitsPt[nBinPt+1];
   for(Int_t iPt = 0;iPt <= nBinPt;iPt++){
@@ -687,7 +694,7 @@ void AliAnalysisTaskJetServices::UserExec(Option_t */*option*/)
   Double_t ptHard = 0; 
   Double_t nTrials = 1; // Trials for MC trigger 
 
-  fh1Trials->Fill("#sum{ntrials}",fAvgTrials); 
+  fh1AvgTrials->Fill("#sum{avg ntrials}",fAvgTrials); 
   AliMCEvent* mcEvent = MCEvent();
   //    AliStack *pStack = 0; 
   if(mcEvent){

diff --git a/PWGJE/AliAnalysisTaskJetServices.h b/PWGJE/AliAnalysisTaskJetServices.h
index 4ffc606f30c5b2dc0ffa2b941d6004fd1d5e4c5c..e00b782237e1d642b3b5c552cc6040dc5a26f455 100644 (file)
--- a/PWGJE/AliAnalysisTaskJetServices.h
+++ b/PWGJE/AliAnalysisTaskJetServices.h
@@ -149,8 +149,9 @@ class AliAnalysisTaskJetServices : public AliAnalysisTaskSE
 
     TString       fNonStdFile;         // outputName for replication
     TString       *fTriggerName;  //[fNTrigger] array of trigger names
-    TProfile*     fh1Xsec;             //! pythia cross section and trials
-    TH1F*         fh1Trials;           //! trials are added
+    TProfile*     fh1Xsec;             //! pythia cross section is averaged
+    TH1F*         fh1Trials;           //! trials are added per notify...
+    TH1F*         fh1AvgTrials;        //! average trials are added per event
     TH1F*         fh1PtHard;           //! Pt har of the event...       
     TH1F*         fh1PtHardTrials;     //! Number of trials 
     TH1F*         fh1SelectionInfoESD; //! Masks that satisfy fSelectionInfo
@@ -192,7 +193,7 @@ class AliAnalysisTaskJetServices : public AliAnalysisTaskSE
     static AliAODHeader*    fgAODHeader;        //! Header for replication
     static AliAODVZERO*    fgAODVZERO;        //! vzero for replication
     static TClonesArray*  fgAODVertices;        //! primary vertex for replication
-    ClassDef(AliAnalysisTaskJetServices,15)
+    ClassDef(AliAnalysisTaskJetServices,16)
 };
  
 #endif

diff --git a/PWGJE/AliAnalysisTaskJetSpectrum2.cxx b/PWGJE/AliAnalysisTaskJetSpectrum2.cxx
index 45a79020068dfed37ecaf05560a5294b1fd07c9a..589e4fba4cf2807b6d76857658eebcbcc1e2f433 100644 (file)
--- a/PWGJE/AliAnalysisTaskJetSpectrum2.cxx
+++ b/PWGJE/AliAnalysisTaskJetSpectrum2.cxx
@@ -123,6 +123,7 @@ AliAnalysisTaskJetSpectrum2::AliAnalysisTaskJetSpectrum2():
   fTriggerName(0x0),
   fh1Xsec(0x0),
   fh1Trials(0x0),
+  fh1AvgTrials(0x0),
   fh1PtHard(0x0),
   fh1PtHardNoW(0x0),  
   fh1PtHardTrials(0x0),
@@ -227,6 +228,7 @@ AliAnalysisTaskJetSpectrum2::AliAnalysisTaskJetSpectrum2(const char* name):
   fTriggerName(0x0),
   fh1Xsec(0x0),
   fh1Trials(0x0),
+  fh1AvgTrials(0x0),
   fh1PtHard(0x0),
   fh1PtHardNoW(0x0),  
   fh1PtHardTrials(0x0),
@@ -326,6 +328,7 @@ Bool_t AliAnalysisTaskJetSpectrum2::Notify()
     // construct a poor man average trials 
     Float_t nEntries = (Float_t)tree->GetTree()->GetEntries();
     if(ftrials>=nEntries && nEntries>0.)fAvgTrials = ftrials/nEntries;
+    fh1Trials->Fill("#sum{ntrials}",ftrials); 
   }  
 
   if(fDebug)Printf("Reading File %s",fInputHandler->GetTree()->GetCurrentFile()->GetName());
@@ -412,6 +415,9 @@ void AliAnalysisTaskJetSpectrum2::UserCreateOutputObjects()
   fh1Trials = new TH1F("fh1Trials","trials root file",1,0,1);
   fh1Trials->GetXaxis()->SetBinLabel(1,"#sum{ntrials}");
   fHistList->Add(fh1Trials);
+  fh1AvgTrials = new TH1F("fh1AvgTrials","trials root file",1,0,1);
+  fh1AvgTrials->GetXaxis()->SetBinLabel(1,"#sum{avg ntrials}");
+  fHistList->Add(fh1AvgTrials);
   fh1PtHard = new TH1F("fh1PtHard","PYTHIA Pt hard;p_{T,hard}",nBinPt,binLimitsPt);
   fHistList->Add(fh1PtHard);
   fh1PtHardNoW = new TH1F("fh1PtHardNoW","PYTHIA Pt hard no weight;p_{T,hard}",nBinPt,binLimitsPt);
@@ -806,7 +812,7 @@ void AliAnalysisTaskJetSpectrum2::UserExec(Option_t */*option*/){
   Double_t eventW = 1;
   Double_t ptHard = 0; 
   Double_t nTrials = 1; // Trials for MC trigger 
-  fh1Trials->Fill("#sum{ntrials}",fAvgTrials); 
+  fh1AvgTrials->Fill("#sum{avg ntrials}",fAvgTrials); 
 
   // Getting some global properties
   fCentrality = GetCentrality();

diff --git a/PWGJE/AliAnalysisTaskJetSpectrum2.h b/PWGJE/AliAnalysisTaskJetSpectrum2.h
index 483f7b2f7b8acd4baefd1ca15b69209f93917421..362d24e1504a333ae9b6e96274b9a01e913a9810 100644 (file)
--- a/PWGJE/AliAnalysisTaskJetSpectrum2.h
+++ b/PWGJE/AliAnalysisTaskJetSpectrum2.h
@@ -203,6 +203,7 @@ class AliAnalysisTaskJetSpectrum2 : public AliAnalysisTaskSE
 
     TProfile*     fh1Xsec;   //! pythia cross section and trials
     TH1F*         fh1Trials; //! trials are added
+    TH1F*         fh1AvgTrials; //! trials are added
     TH1F*         fh1PtHard;  //! Pt har of the event...       
     TH1F*         fh1PtHardNoW;  //! Pt har of the event without weigt       
     TH1F*         fh1PtHardTrials;  //! Number of trials 
@@ -254,7 +255,7 @@ class AliAnalysisTaskJetSpectrum2 : public AliAnalysisTaskSE
     TList *fHistList;                  //! Output list
    
 
-    ClassDef(AliAnalysisTaskJetSpectrum2, 22); // Analysis task for standard jet analysis
+    ClassDef(AliAnalysisTaskJetSpectrum2, 23); // Analysis task for standard jet analysis
 };
  
 #endif