Add correlation with EP
authormverweij <marta.verweij@cern.ch>
Thu, 20 Nov 2014 15:42:21 +0000 (16:42 +0100)
committermverweij <marta.verweij@cern.ch>
Thu, 20 Nov 2014 21:22:16 +0000 (22:22 +0100)
reduce output

PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskEmcalJetMass.cxx
PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskEmcalJetMass.h
PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskEmcalJetTagger.cxx
PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskEmcalJetTagger.h

index 480704e..418cf6c 100644 (file)
@@ -9,23 +9,13 @@
 #include <TH3F.h>
 #include <THnSparse.h>
 #include <TList.h>
-#include <TLorentzVector.h>
 #include <TProfile.h>
-#include <TChain.h>
-#include <TSystem.h>
-#include <TFile.h>
-#include <TKey.h>
 
-#include "AliVCluster.h"
 #include "AliVTrack.h"
 #include "AliEmcalJet.h"
 #include "AliRhoParameter.h"
 #include "AliLog.h"
 #include "AliEmcalParticle.h"
-#include "AliMCEvent.h"
-#include "AliGenPythiaEventHeader.h"
-#include "AliAODMCHeader.h"
-#include "AliMCEvent.h"
 #include "AliAnalysisManager.h"
 #include "AliJetContainer.h"
 
@@ -43,6 +33,7 @@ AliAnalysisTaskEmcalJetMass::AliAnalysisTaskEmcalJetMass() :
   fMinFractionShared(0),
   fUseUnsubJet(0),
   fJetMassType(kRaw),
+  fUseSumw2(0),
   fh3PtJet1VsMassVsLeadPtAllSel(0),
   fh3PtJet1VsMassVsLeadPtTagged(0),
   fh3PtJet1VsMassVsLeadPtTaggedMatch(0),
@@ -67,7 +58,9 @@ AliAnalysisTaskEmcalJetMass::AliAnalysisTaskEmcalJetMass() :
   fh2RatVsNConstJet1All(0),
   fh2RatVsNConstJet1Tagged(0),
   fh2RatVsNConstJet1TaggedMatch(0),
-  fh2EtMassOverEtRSq(0)
+  fh3JetPtVsMassVsEPRelAllSel(0),
+  fh3JetPtVsMassVsEPRelTagged(0),
+  fh3JetPtVsMassVsEPRelTaggedMatch(0)
 {
   // Default constructor.
 
@@ -84,20 +77,22 @@ AliAnalysisTaskEmcalJetMass::AliAnalysisTaskEmcalJetMass() :
   fh2MassVsNConstJet1Tagged            = new TH2F*[fNcentBins];
   fh2MassVsNConstJet1TaggedMatch       = new TH2F*[fNcentBins];
 
-  fh3PtJet1VsRatVsLeadPtAllSel        = new TH3F*[fNcentBins];
-  fh3PtJet1VsRatVsLeadPtTagged        = new TH3F*[fNcentBins];
-  fh3PtJet1VsRatVsLeadPtTaggedMatch   = new TH3F*[fNcentBins];
-  fpPtVsRatJet1All                    = new TProfile*[fNcentBins];
-  fpPtVsRatJet1Tagged                 = new TProfile*[fNcentBins];
-  fpPtVsRatJet1TaggedMatch            = new TProfile*[fNcentBins];
-  fh2RatVsAreaJet1All                 = new TH2F*[fNcentBins];
-  fh2RatVsAreaJet1Tagged              = new TH2F*[fNcentBins];
-  fh2RatVsAreaJet1TaggedMatch         = new TH2F*[fNcentBins];
-  fh2RatVsNConstJet1All               = new TH2F*[fNcentBins];
-  fh2RatVsNConstJet1Tagged            = new TH2F*[fNcentBins];
-  fh2RatVsNConstJet1TaggedMatch       = new TH2F*[fNcentBins];
-
-  fh2EtMassOverEtRSq                   = new TH2F*[fNcentBins];
+  fh3PtJet1VsRatVsLeadPtAllSel         = new TH3F*[fNcentBins];
+  fh3PtJet1VsRatVsLeadPtTagged         = new TH3F*[fNcentBins];
+  fh3PtJet1VsRatVsLeadPtTaggedMatch    = new TH3F*[fNcentBins];
+  fpPtVsRatJet1All                     = new TProfile*[fNcentBins];
+  fpPtVsRatJet1Tagged                  = new TProfile*[fNcentBins];
+  fpPtVsRatJet1TaggedMatch             = new TProfile*[fNcentBins];
+  fh2RatVsAreaJet1All                  = new TH2F*[fNcentBins];
+  fh2RatVsAreaJet1Tagged               = new TH2F*[fNcentBins];
+  fh2RatVsAreaJet1TaggedMatch          = new TH2F*[fNcentBins];
+  fh2RatVsNConstJet1All                = new TH2F*[fNcentBins];
+  fh2RatVsNConstJet1Tagged             = new TH2F*[fNcentBins];
+  fh2RatVsNConstJet1TaggedMatch        = new TH2F*[fNcentBins];
+
+  fh3JetPtVsMassVsEPRelAllSel          = new TH3F*[fNcentBins];
+  fh3JetPtVsMassVsEPRelTagged          = new TH3F*[fNcentBins];
+  fh3JetPtVsMassVsEPRelTaggedMatch     = new TH3F*[fNcentBins];
 
   for (Int_t i = 0; i < fNcentBins; i++) {
     fh3PtJet1VsMassVsLeadPtAllSel[i]        = 0;
@@ -113,20 +108,22 @@ AliAnalysisTaskEmcalJetMass::AliAnalysisTaskEmcalJetMass() :
     fh2MassVsNConstJet1Tagged[i]            = 0;
     fh2MassVsNConstJet1TaggedMatch[i]       = 0;
 
-    fh3PtJet1VsRatVsLeadPtAllSel[i]        = 0;
-    fh3PtJet1VsRatVsLeadPtTagged[i]        = 0;
-    fh3PtJet1VsRatVsLeadPtTaggedMatch[i]   = 0;
-    fpPtVsRatJet1All[i]                    = 0;
-    fpPtVsRatJet1Tagged[i]                 = 0;
-    fpPtVsRatJet1TaggedMatch[i]            = 0;
-    fh2RatVsAreaJet1All[i]                 = 0;
-    fh2RatVsAreaJet1Tagged[i]              = 0;
-    fh2RatVsAreaJet1TaggedMatch[i]         = 0;
-    fh2RatVsNConstJet1All[i]               = 0;
-    fh2RatVsNConstJet1Tagged[i]            = 0;
-    fh2RatVsNConstJet1TaggedMatch[i]       = 0;
-
-    fh2EtMassOverEtRSq[i]                  = 0;
+    fh3PtJet1VsRatVsLeadPtAllSel[i]         = 0;
+    fh3PtJet1VsRatVsLeadPtTagged[i]         = 0;
+    fh3PtJet1VsRatVsLeadPtTaggedMatch[i]    = 0;
+    fpPtVsRatJet1All[i]                     = 0;
+    fpPtVsRatJet1Tagged[i]                  = 0;
+    fpPtVsRatJet1TaggedMatch[i]             = 0;
+    fh2RatVsAreaJet1All[i]                  = 0;
+    fh2RatVsAreaJet1Tagged[i]               = 0;
+    fh2RatVsAreaJet1TaggedMatch[i]          = 0;
+    fh2RatVsNConstJet1All[i]                = 0;
+    fh2RatVsNConstJet1Tagged[i]             = 0;
+    fh2RatVsNConstJet1TaggedMatch[i]        = 0;
+
+    fh3JetPtVsMassVsEPRelAllSel[i]          = 0;
+    fh3JetPtVsMassVsEPRelTagged[i]          = 0;
+    fh3JetPtVsMassVsEPRelTaggedMatch[i]     = 0;
   }
 
   SetMakeGeneralHistograms(kTRUE);
@@ -140,6 +137,7 @@ AliAnalysisTaskEmcalJetMass::AliAnalysisTaskEmcalJetMass(const char *name) :
   fMinFractionShared(0),
   fUseUnsubJet(0),
   fJetMassType(kRaw),
+  fUseSumw2(0),
   fh3PtJet1VsMassVsLeadPtAllSel(0),
   fh3PtJet1VsMassVsLeadPtTagged(0),
   fh3PtJet1VsMassVsLeadPtTaggedMatch(0),
@@ -164,7 +162,9 @@ AliAnalysisTaskEmcalJetMass::AliAnalysisTaskEmcalJetMass(const char *name) :
   fh2RatVsNConstJet1All(0),
   fh2RatVsNConstJet1Tagged(0),
   fh2RatVsNConstJet1TaggedMatch(0),
-  fh2EtMassOverEtRSq(0)
+  fh3JetPtVsMassVsEPRelAllSel(0),
+  fh3JetPtVsMassVsEPRelTagged(0),
+  fh3JetPtVsMassVsEPRelTaggedMatch(0)
 {
   // Standard constructor.
 
@@ -194,7 +194,9 @@ AliAnalysisTaskEmcalJetMass::AliAnalysisTaskEmcalJetMass(const char *name) :
   fh2RatVsNConstJet1Tagged             = new TH2F*[fNcentBins];
   fh2RatVsNConstJet1TaggedMatch        = new TH2F*[fNcentBins];
 
-  fh2EtMassOverEtRSq                   = new TH2F*[fNcentBins];
+  fh3JetPtVsMassVsEPRelAllSel          = new TH3F*[fNcentBins];
+  fh3JetPtVsMassVsEPRelTagged          = new TH3F*[fNcentBins];
+  fh3JetPtVsMassVsEPRelTaggedMatch     = new TH3F*[fNcentBins];
 
   for (Int_t i = 0; i < fNcentBins; i++) {
     fh3PtJet1VsMassVsLeadPtAllSel[i]        = 0;
@@ -223,7 +225,9 @@ AliAnalysisTaskEmcalJetMass::AliAnalysisTaskEmcalJetMass(const char *name) :
     fh2RatVsNConstJet1Tagged[i]             = 0;
     fh2RatVsNConstJet1TaggedMatch[i]        = 0;
 
-    fh2EtMassOverEtRSq[i]                   = 0;
+    fh3JetPtVsMassVsEPRelAllSel[i]          = 0;
+    fh3JetPtVsMassVsEPRelTagged[i]          = 0;
+    fh3JetPtVsMassVsEPRelTaggedMatch[i]     = 0;
   }
 
   SetMakeGeneralHistograms(kTRUE);
@@ -389,22 +393,33 @@ void AliAnalysisTaskEmcalJetMass::UserCreateOutputObjects()
     fh2RatVsNConstJet1TaggedMatch[i] = new TH2F(histName.Data(),histTitle.Data(),nBinsR,minR,maxR,nBinsNConst,minNConst,maxNConst);
     fOutput->Add(fh2RatVsNConstJet1TaggedMatch[i]);
 
+    histName = TString::Format("fh3JetPtVsMassVsEPRelAllSel_%d",i);
+    histTitle = TString::Format("%s;#it{p}_{T,jet1};#it{M}_{jet1};#varphi_{jet}-#Psi",histName.Data());
+    fh3JetPtVsMassVsEPRelAllSel[i] = new TH3F(histName.Data(),histTitle.Data(),nBinsPt,minPt,maxPt,nBinsM,minM,maxM,100,0.,TMath::Pi());
+    fOutput->Add(fh3JetPtVsMassVsEPRelAllSel[i]);
 
-    histName = TString::Format("fh2EtMassOverEtRSq_%d",i);
-    histTitle = TString::Format("%s;#it{E}_{T};(#it{M}/(#it{E}_{T}#it{R}))^{2}",histName.Data());
-    fh2EtMassOverEtRSq[i] = new TH2F(histName.Data(),histTitle.Data(),nBinsPt,minPt,maxPt,100,0.,1.);
-    fOutput->Add(fh2EtMassOverEtRSq[i]);
+    histName = TString::Format("fh3JetPtVsMassVsEPRelTagged_%d",i);
+    histTitle = TString::Format("%s;#it{p}_{T,jet1};#it{M}_{jet1};#varphi_{jet}-#Psi",histName.Data());
+    fh3JetPtVsMassVsEPRelTagged[i] =  new TH3F(histName.Data(),histTitle.Data(),nBinsPt,minPt,maxPt,nBinsM,minM,maxM,100,0.,TMath::Pi());
+    fOutput->Add(fh3JetPtVsMassVsEPRelTagged[i]);
+
+    histName = TString::Format("fh3JetPtVsMassVsEPRelTaggedMatch_%d",i);
+    histTitle = TString::Format("%s;#it{p}_{T,jet1};#it{M}_{jet1};#varphi_{jet}-#Psi",histName.Data());
+    fh3JetPtVsMassVsEPRelTaggedMatch[i] =  new TH3F(histName.Data(),histTitle.Data(),nBinsPt,minPt,maxPt,nBinsM,minM,maxM,100,0.,TMath::Pi());
+    fOutput->Add(fh3JetPtVsMassVsEPRelTaggedMatch[i]);
   }
 
-  // =========== Switch on Sumw2 for all histos ===========
-  for (Int_t i=0; i<fOutput->GetEntries(); ++i) {
-    TH1 *h1 = dynamic_cast<TH1*>(fOutput->At(i));
-    if (h1){
-      h1->Sumw2();
-      continue;
+  if(fUseSumw2) {
+    // =========== Switch on Sumw2 for all histos ===========
+    for (Int_t i=0; i<fOutput->GetEntries(); ++i) {
+      TH1 *h1 = dynamic_cast<TH1*>(fOutput->At(i));
+      if (h1){
+       h1->Sumw2();
+       continue;
+      }
+      THnSparse *hn = dynamic_cast<THnSparse*>(fOutput->At(i));
+      if(hn)hn->Sumw2();
     }
-    THnSparse *hn = dynamic_cast<THnSparse*>(fOutput->At(i));
-    if(hn)hn->Sumw2();
   }
 
   TH1::AddDirectory(oldStatus);
@@ -436,6 +451,9 @@ Bool_t AliAnalysisTaskEmcalJetMass::FillHistograms()
       Double_t mJet1 = GetJetMass(jet1);
       Double_t rat = -1.;
       if(ptJet1<0. || ptJet1>0.) rat = mJet1/ptJet1;
+      Double_t ep = jet1->Phi() - fEPV0;
+      while (ep < 0) ep += TMath::Pi();
+      while (ep >= TMath::Pi()) ep -= TMath::Pi();
 
       Double_t fraction = -1.;
       if(fUseUnsubJet) {
@@ -462,6 +480,7 @@ Bool_t AliAnalysisTaskEmcalJetMass::FillHistograms()
       fpPtVsMassJet1All[fCentBin]->Fill(ptJet1,mJet1);
       fh2MassVsAreaJet1All[fCentBin]->Fill(mJet1,jet1->Area());
       fh2MassVsNConstJet1All[fCentBin]->Fill(mJet1,jet1->GetNumberOfConstituents());
+      fh3JetPtVsMassVsEPRelAllSel[fCentBin]->Fill(ptJet1,mJet1,ep);
 
       fh3PtJet1VsRatVsLeadPtAllSel[fCentBin]->Fill(ptJet1,rat,maxTrackPt);
       fpPtVsRatJet1All[fCentBin]->Fill(ptJet1,rat);
@@ -475,6 +494,7 @@ Bool_t AliAnalysisTaskEmcalJetMass::FillHistograms()
       fpPtVsMassJet1Tagged[fCentBin]->Fill(ptJet1,mJet1);
       fh2MassVsAreaJet1Tagged[fCentBin]->Fill(mJet1,jet1->Area());
       fh2MassVsNConstJet1Tagged[fCentBin]->Fill(mJet1,jet1->GetNumberOfConstituents());
+      fh3JetPtVsMassVsEPRelTagged[fCentBin]->Fill(ptJet1,mJet1,ep);
 
       fh3PtJet1VsRatVsLeadPtTagged[fCentBin]->Fill(ptJet1,rat,maxTrackPt);
       fpPtVsRatJet1Tagged[fCentBin]->Fill(ptJet1,rat);
@@ -483,26 +503,21 @@ Bool_t AliAnalysisTaskEmcalJetMass::FillHistograms()
 
       //matching
       if(fMinFractionShared>0. && fraction>fMinFractionShared) {
+       AliEmcalJet *jetM = jet1->ClosestJet();
+       maxTrackPt = jetM->MaxTrackPt();
        fh3PtJet1VsMassVsLeadPtTaggedMatch[fCentBin]->Fill(ptJet1,mJet1,maxTrackPt);
        fpPtVsMassJet1TaggedMatch[fCentBin]->Fill(ptJet1,mJet1);
        fh2MassVsAreaJet1TaggedMatch[fCentBin]->Fill(mJet1,jet1->Area());
        fh2MassVsNConstJet1TaggedMatch[fCentBin]->Fill(mJet1,jet1->GetNumberOfConstituents());
+       fh3JetPtVsMassVsEPRelTaggedMatch[fCentBin]->Fill(ptJet1,mJet1,ep);
 
        fh3PtJet1VsRatVsLeadPtTaggedMatch[fCentBin]->Fill(ptJet1,rat,maxTrackPt);
        fpPtVsRatJet1TaggedMatch[fCentBin]->Fill(ptJet1,rat);
        fh2RatVsAreaJet1TaggedMatch[fCentBin]->Fill(rat,jet1->Area());
        fh2RatVsNConstJet1TaggedMatch[fCentBin]->Fill(rat,jet1->GetNumberOfConstituents());
       }
-      
-      Double_t Et2 = mJet1*mJet1 + ptJet1*ptJet1;
-      Double_t Et = 0.;    Double_t massOverEtR = 0.;
-      if(Et2>0.) Et = TMath::Sqrt(Et2);
-      if((Et*jetCont->GetJetRadius())>0.) 
-       massOverEtR = mJet1/(Et*jetCont->GetJetRadius());
-      fh2EtMassOverEtRSq[fCentBin]->Fill(Et,massOverEtR*massOverEtR);
     }
   }
-  
   return kTRUE;
 }
 
index 6403368..782345d 100644 (file)
@@ -32,6 +32,7 @@ class AliAnalysisTaskEmcalJetMass : public AliAnalysisTaskEmcalJet {
   void SetJetContainerUnsub(Int_t c)                                 { fContainerUnsub    = c   ; }
   void SetMinFractionShared(Double_t f, Bool_t useUnsubJet = kFALSE) { fMinFractionShared = f   ; fUseUnsubJet = useUnsubJet; }
   void SetJetMassType(JetMassType t)                                 { fJetMassType       = t   ; }
+  void SetUseSumw2(Bool_t b)                                         { fUseSumw2          = b   ; }
 
  protected:
   Bool_t                              RetrieveEventObjects();
@@ -40,11 +41,13 @@ class AliAnalysisTaskEmcalJetMass : public AliAnalysisTaskEmcalJet {
 
   Double_t                            GetJetMass(AliEmcalJet *jet);
  
+ private:
   Int_t                               fContainerBase;              // jets to be analyzed
   Int_t                               fContainerUnsub;             // unsubtracted jets
   Double_t                            fMinFractionShared;          // only fill histos for jets if shared fraction larger than X
   Bool_t                              fUseUnsubJet;                // calc fraction of unsubtracted jet (relevant for constituent subtraction)
   JetMassType                         fJetMassType;                // jet mass type to be used
+  Bool_t                              fUseSumw2;                   // activate sumw2 for output histograms
 
   TH3F            **fh3PtJet1VsMassVsLeadPtAllSel;         //!all jets after std selection pt vs mass vs leading track pt
   TH3F            **fh3PtJet1VsMassVsLeadPtTagged;         //!tagged jets pt vs mass vs leading track pt
@@ -72,13 +75,14 @@ class AliAnalysisTaskEmcalJetMass : public AliAnalysisTaskEmcalJet {
   TH2F            **fh2RatVsNConstJet1Tagged;              //!mass/pt vs number of constituents of tagged jets
   TH2F            **fh2RatVsNConstJet1TaggedMatch;         //!mass/pt vs number of constituents of tagged jets matched to MC
 
-  TH2F            **fh2EtMassOverEtRSq;                    //!Et vs (M/Et*R)^2
+  TH3F            **fh3JetPtVsMassVsEPRelAllSel;           //!jet pt vs mass vs (phi-ep) for all jets
+  TH3F            **fh3JetPtVsMassVsEPRelTagged;           //!jet pt vs mass vs (phi-ep) for tagged jets
+  TH3F            **fh3JetPtVsMassVsEPRelTaggedMatch;      //!jet pt vs mass vs (phi-ep) for tagged matched jets
 
- private:
   AliAnalysisTaskEmcalJetMass(const AliAnalysisTaskEmcalJetMass&);            // not implemented
   AliAnalysisTaskEmcalJetMass &operator=(const AliAnalysisTaskEmcalJetMass&); // not implemented
 
-  ClassDef(AliAnalysisTaskEmcalJetMass, 8)
+  ClassDef(AliAnalysisTaskEmcalJetMass, 9)
 };
 #endif
 
index 47b4efd..0a7c4c2 100644 (file)
@@ -45,6 +45,7 @@ AliAnalysisTaskEmcalJetTagger::AliAnalysisTaskEmcalJetTagger() :
   fContainerBase(0),
   fContainerTag(1),
   fMinFractionShared(0),
+  fUseSumw2(0),
   fMatchingDone(0),
   fh3PtJet1VsDeltaEtaDeltaPhi(0),
   fh2PtJet1VsDeltaR(0),
@@ -88,6 +89,7 @@ AliAnalysisTaskEmcalJetTagger::AliAnalysisTaskEmcalJetTagger(const char *name) :
   fContainerBase(0),
   fContainerTag(1),
   fMinFractionShared(0),
+  fUseSumw2(0),
   fMatchingDone(0),
   fh3PtJet1VsDeltaEtaDeltaPhi(0),
   fh2PtJet1VsDeltaR(0),
@@ -140,7 +142,7 @@ void AliAnalysisTaskEmcalJetTagger::UserCreateOutputObjects()
   TH1::AddDirectory(kFALSE);
 
   const Int_t nBinsPt          = 250;
-  const Int_t nBinsDPhi        = 100;
+  const Int_t nBinsDPhi        = 72;
   const Int_t nBinsDEta        = 100;
   const Int_t nBinsDR          = 50;
   const Int_t nBinsFraction    = 101;
@@ -205,15 +207,17 @@ void AliAnalysisTaskEmcalJetTagger::UserCreateOutputObjects()
   fh3PtJetAreaDRConst = new TH3F("fh3PtJetAreaDRConst","fh3PtJetAreaDRConst;pT;A;#Delta R",nBinsPt,minPt,maxPt,100,0.,1.,100,0.,1.);
   fOutput->Add(fh3PtJetAreaDRConst);
 
-  // =========== Switch on Sumw2 for all histos ===========
-  for (Int_t i=0; i<fOutput->GetEntries(); ++i) {
-    TH1 *h1 = dynamic_cast<TH1*>(fOutput->At(i));
-    if (h1){
-      h1->Sumw2();
-      continue;
+  if(fUseSumw2) {
+    // =========== Switch on Sumw2 for all histos ===========
+    for (Int_t i=0; i<fOutput->GetEntries(); ++i) {
+      TH1 *h1 = dynamic_cast<TH1*>(fOutput->At(i));
+      if (h1){
+       h1->Sumw2();
+       continue;
+      }
+      THnSparse *hn = dynamic_cast<THnSparse*>(fOutput->At(i));
+      if(hn)hn->Sumw2();
     }
-    THnSparse *hn = dynamic_cast<THnSparse*>(fOutput->At(i));
-    if(hn)hn->Sumw2();
   }
 
   TH1::AddDirectory(oldStatus);
index e3fac35..bfefea0 100644 (file)
@@ -41,6 +41,8 @@ class AliAnalysisTaskEmcalJetTagger : public AliAnalysisTaskEmcalJet {
 
   void SetMinFractionShared(Double_t f)                         { fMinFractionShared = f; }
 
+  void SetUseSumw2(Bool_t b)                                    { fUseSumw2 = b;}
+
   //Getters
   Double_t GetFractionSharedPt(const AliEmcalJet *jet1, const AliEmcalJet *jet2) const;
 
@@ -57,16 +59,15 @@ class AliAnalysisTaskEmcalJetTagger : public AliAnalysisTaskEmcalJet {
   void     MatchJetsGeo(Int_t c1 = -1, Int_t c2 = -1, Int_t iDebug = 0, Float_t maxDist = 0.3, Int_t type = 2, Bool_t bReset = kTRUE);
   void     ResetTagging(const Int_t c);
   
+ private:
   JetTaggingType                      fJetTaggingType;             // jet matching type
   JetTaggingMethod                    fJetTaggingMethod;           // jet matching method
-
   Int_t                               fContainerBase;              // jets to be tagged
   Int_t                               fContainerTag;               // jets used for tagging
-
   Double_t                            fMinFractionShared;          // only fill histos for jets if shared fraction larger than X
-  
- private:
-  Bool_t            fMatchingDone;                // flag to indicate if matching is done or not
+  Bool_t                              fUseSumw2;                   // activate sumw2 for output histograms
+  Bool_t                              fMatchingDone;               // flag to indicate if matching is done or not
+
   TH3F            **fh3PtJet1VsDeltaEtaDeltaPhi;  //!pt jet 1 vs deta vs dphi
   TH2F            **fh2PtJet1VsDeltaR;            //!pt jet 1 vs dR
   TH2F            **fh2PtJet2VsFraction;          //!pt jet 1 vs shared fraction
@@ -83,7 +84,7 @@ class AliAnalysisTaskEmcalJetTagger : public AliAnalysisTaskEmcalJet {
   AliAnalysisTaskEmcalJetTagger(const AliAnalysisTaskEmcalJetTagger&);            // not implemented
   AliAnalysisTaskEmcalJetTagger &operator=(const AliAnalysisTaskEmcalJetTagger&); // not implemented
 
-  ClassDef(AliAnalysisTaskEmcalJetTagger, 5)
+  ClassDef(AliAnalysisTaskEmcalJetTagger, 6)
 };
 #endif