added more info and fixed calc

diff --git a/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskCLQA.cxx b/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskCLQA.cxx
index 46f99f8b2fa8109220a058dc487dc02490d6b12d..532297114dd3825f6688ee4c059998597a0be3aa 100644 (file)
--- a/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskCLQA.cxx
+++ b/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskCLQA.cxx
@@ -105,8 +105,8 @@ void AliAnalysisTaskCLQA::RunCumulants(Double_t Mmin, Double_t ptmin, Double_t p
     return;
 
   const Int_t ntracks = fTracks->GetEntries();
-  Int_t Mall=0,M=0;
-  Double_t ptmaxall=0,ptsumall=0,pt2sumall=0;
+  Int_t Mall=0,M=0,Mall2;
+  Double_t ptmaxall=0,ptsumall=0,pt2sumall=0,ptsumall2;
   Double_t tsa00=0,tsa10=0,tsa11=0;
   Double_t Q2r=0,Q2i=0;
   Double_t Q4r=0,Q4i=0;
@@ -120,6 +120,10 @@ void AliAnalysisTaskCLQA::RunCumulants(Double_t Mmin, Double_t ptmin, Double_t p
     Double_t pt = track->Pt();
     if (pt>ptmaxall)
       ptmaxall = pt;
+    if (pt>2) {
+      ptsumall2 += pt;
+      ++Mall2;
+    }
     ptsumall  +=pt;
     pt2sumall +=pt*pt;
     Double_t px = track->Px();
@@ -150,7 +154,7 @@ void AliAnalysisTaskCLQA::RunCumulants(Double_t Mmin, Double_t ptmin, Double_t p
 
   Double_t Q2abs = Q2r*Q2r+Q2i*Q2i;
   Double_t Q4abs = Q4r*Q4r+Q4i*Q4i;
-  Double_t Q42re = Q4r*Q2r*Q2r-Q4r*Q2i*Q2i+2*Q2r*Q2i;
+  Double_t Q42re = Q4r*Q2r*Q2r-Q4r*Q2i*Q2i+2*Q4i*Q2r*Q2i;
 
   Double_t tsall = -1;
   Double_t tsax = (tsa00+tsa11)*(tsa00+tsa11)-4*(tsa00*tsa11-tsa10*tsa10);
@@ -184,9 +188,11 @@ void AliAnalysisTaskCLQA::RunCumulants(Double_t Mmin, Double_t ptmin, Double_t p
                                                    5./*nSigmaDiamZ*/);
 
   fNtupCumInfo->fMall     = Mall;
+  fNtupCumInfo->fMall2    = Mall2;
   fNtupCumInfo->fPtMaxall = ptmaxall;
   fNtupCumInfo->fMPtall   = ptsumall/Mall;
   fNtupCumInfo->fMPt2all  = pt2sumall/Mall;
+  fNtupCumInfo->fMPtall2  = ptsumall2/Mall2;
   fNtupCumInfo->fTSall    = tsall;
   fNtupCumInfo->fM        = M;
   fNtupCumInfo->fQ2abs    = Q2abs;

diff --git a/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskCLQA.h b/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskCLQA.h
index 8f7f5ba33c39505d5d2afcb3972181c4cfe71162..4ffafeea39512527f30aa0fb4868949fb34c3667 100644 (file)
--- a/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskCLQA.h
+++ b/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskCLQA.h
@@ -51,10 +51,10 @@ class AliAnalysisTaskCLQA : public AliAnalysisTaskEmcalJet {
 class AliNtupCumInfo {
  public:
     AliNtupCumInfo() : fTrig(0), fRun(0), fVz(0), fIsFEC(0), fIsVSel(0), fIsP(0),
-                       fMall(0), fPtMaxall(0), fMPtall(0), fMPt2all(0), fTSall(0),
+                       fMall(0), fMall2(0), fPtMaxall(0), fMPtall(0), 
+                       fMPt2all(0), fMPtall2(0), fTSall(0),
                        fM(0), fQ2abs(0), fQ4abs(0), fQ42re(0),
-                       fPtMax(0), fMPt(0), fMPt2(0), fTS(0), 
-                       fMV0M(0), 
+                       fPtMax(0), fMPt(0), fMPt2(0), fTS(0), fMV0M(0), 
                        fCl1(0), fV0M(0), fV0MEq(0), fV0A(0), fV0AEq(0), fZNA(0) {;}
   virtual ~AliNtupCumInfo() {;}
 
@@ -66,9 +66,11 @@ class AliNtupCumInfo {
   Bool_t        fIsVSel;       // is vertex selected
   Bool_t        fIsP;          // is SPD pileup
   Int_t         fMall;         // multiplicity (tracks in eta range)
+  Int_t         fMall2;        // multiplicity (tracks above 2 GeV/c in eta range)
   Double32_t    fPtMaxall;     //[0,0,16] maximum pT
   Double32_t    fMPtall;       //[0,0,16] mean pT
   Double32_t    fMPt2all;      //[0,0,16] mean pT2
+  Double32_t    fMPtall2;      //[0,0,16] mean pT truncated above 2 GeV/c
   Double32_t    fTSall;        //[0,0,16] transverse sphericity
   Int_t         fM;            // multiplicity (tracks in pT range)
   Double_t      fQ2abs;        // Q2 absolute
@@ -86,7 +88,7 @@ class AliNtupCumInfo {
   Double32_t    fV0AEq;        //[0,0,16] class V0A Eq
   Double32_t    fZNA;          //[0,0,16] class ZNA
  
-  ClassDef(AliNtupCumInfo,1) // Cumulant storage class
+  ClassDef(AliNtupCumInfo,2) // Cumulant storage class
 };
 
 #endif