Changes by Salvatore
authorloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 3 Jan 2013 09:37:04 +0000 (09:37 +0000)
committerloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 3 Jan 2013 09:37:04 +0000 (09:37 +0000)
PWGJE/EMCALJetTasks/AliAnalysisTaskDeltaPt.cxx
PWGJE/EMCALJetTasks/AliAnalysisTaskDeltaPt.h
PWGJE/EMCALJetTasks/AliAnalysisTaskRhoBase.cxx
PWGJE/EMCALJetTasks/AliAnalysisTaskRhoBase.h
PWGJE/EMCALJetTasks/AliEmcalJetTask.cxx
PWGJE/EMCALJetTasks/AliEmcalJetTask.h
PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskSAQA.cxx
PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C
PWGJE/EMCALJetTasks/macros/AddTaskJetRespPtHard.C
PWGJE/EMCALJetTasks/macros/AddTaskSAJF.C
PWGJE/EMCALJetTasks/macros/AddTaskSAQA.C

index 8f9a6d3..3f3e700 100644 (file)
@@ -56,13 +56,13 @@ AliAnalysisTaskDeltaPt::AliAnalysisTaskDeltaPt() :
     fHistDeltaPtRCExLJ[i] = 0;
     fHistDeltaPtRCRand[i] = 0;
     fHistEmbNotFoundPhiEta[i] = 0;
-    fHistEmbJetsPt[i] = 0;
-    fHistEmbJetsCorrPt[i] = 0;
-    fHistEmbJetsArea[i] = 0;
+    fHistEmbJetsPtArea[i] = 0;
+    fHistEmbJetsCorrPtArea[i] = 0;
     fHistEmbPartPt[i] = 0;
     fHistDistEmbPartJetAxis[i] = 0;
+    fHistEmbBkgArea[i] = 0;
     fHistRhoVSEmbBkg[i] = 0;
-    fHistDeltaPtEmb[i] = 0;
+    fHistDeltaPtEmbArea[i] = 0;
   }
 
   SetMakeGeneralHistograms(kTRUE);
@@ -102,13 +102,13 @@ AliAnalysisTaskDeltaPt::AliAnalysisTaskDeltaPt(const char *name) :
     fHistDeltaPtRCExLJ[i] = 0;
     fHistDeltaPtRCRand[i] = 0;
     fHistEmbNotFoundPhiEta[i] = 0;
-    fHistEmbJetsPt[i] = 0;
-    fHistEmbJetsCorrPt[i] = 0;
-    fHistEmbJetsArea[i] = 0;
+    fHistEmbJetsPtArea[i] = 0;
+    fHistEmbJetsCorrPtArea[i] = 0;
     fHistEmbPartPt[i] = 0;
     fHistDistEmbPartJetAxis[i] = 0;
+    fHistEmbBkgArea[i] = 0;
     fHistRhoVSEmbBkg[i] = 0;
-    fHistDeltaPtEmb[i] = 0;
+    fHistDeltaPtEmbArea[i] = 0;
   }
 
   SetMakeGeneralHistograms(kTRUE);
@@ -133,7 +133,7 @@ void AliAnalysisTaskDeltaPt::UserCreateOutputObjects()
   fOutput->Add(fHistRCPhiEta);
 
   fHistRCPtExLJVSDPhiLJ = new TH2F("fHistRCPtExLJVSDPhiLJ","fHistRCPtExLJVSDPhiLJ", fNbins, fMinBinPt, fMaxBinPt, 128, -1.6, 4.8);
-  fHistRCPtExLJVSDPhiLJ->GetXaxis()->SetTitle("rigid cone p_{T} [GeV/c]");
+  fHistRCPtExLJVSDPhiLJ->GetXaxis()->SetTitle("rigid cone #it{p}_{T} (GeV/#it{c})");
   fHistRCPtExLJVSDPhiLJ->GetYaxis()->SetTitle("#Delta#phi");
   fOutput->Add(fHistRCPtExLJVSDPhiLJ);
 
@@ -155,78 +155,71 @@ void AliAnalysisTaskDeltaPt::UserCreateOutputObjects()
     histname = "fHistRCPt_";
     histname += i;
     fHistRCPt[i] = new TH1F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt * 2);
-    fHistRCPt[i]->GetXaxis()->SetTitle("rigid cone p_{T} [GeV/c]");
+    fHistRCPt[i]->GetXaxis()->SetTitle("rigid cone #it{p}_{T} (GeV/#it{c})");
     fHistRCPt[i]->GetYaxis()->SetTitle("counts");
     fOutput->Add(fHistRCPt[i]);
 
     histname = "fHistRCPtExLJ_";
     histname += i;
     fHistRCPtExLJ[i] = new TH1F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt * 2);
-    fHistRCPtExLJ[i]->GetXaxis()->SetTitle("rigid cone p_{T}^{RC} [GeV/c]");
+    fHistRCPtExLJ[i]->GetXaxis()->SetTitle("rigid cone #it{p}_{T}^{RC} (GeV/#it{c})");
     fHistRCPtExLJ[i]->GetYaxis()->SetTitle("counts");
     fOutput->Add(fHistRCPtExLJ[i]);
 
     histname = "fHistRCPtRand_";
     histname += i;
     fHistRCPtRand[i] = new TH1F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt * 2);
-    fHistRCPtRand[i]->GetXaxis()->SetTitle("rigid cone p_{T}^{RC} [GeV/c]");
+    fHistRCPtRand[i]->GetXaxis()->SetTitle("rigid cone #it{p}_{T}^{RC} (GeV/#it{c})");
     fHistRCPtRand[i]->GetYaxis()->SetTitle("counts");
     fOutput->Add(fHistRCPtRand[i]);
 
     histname = "fHistRhoVSRCPt_";
     histname += i;
-    fHistRhoVSRCPt[i] = new TH2F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt*2, fNbins, fMinBinPt, fMaxBinPt);
-    fHistRhoVSRCPt[i]->GetXaxis()->SetTitle("A#rho [GeV/c]");
-    fHistRhoVSRCPt[i]->GetYaxis()->SetTitle("rigid cone p_{T} [GeV/c]");
+    fHistRhoVSRCPt[i] = new TH2F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt, fNbins, fMinBinPt, fMaxBinPt);
+    fHistRhoVSRCPt[i]->GetXaxis()->SetTitle("A#rho (GeV/#it{c})");
+    fHistRhoVSRCPt[i]->GetYaxis()->SetTitle("rigid cone #it{p}_{T} (GeV/#it{c})");
     fOutput->Add(fHistRhoVSRCPt[i]);
 
     histname = "fHistDeltaPtRC_";
     histname += i;
     fHistDeltaPtRC[i] = new TH1F(histname.Data(), histname.Data(), fNbins * 2, -fMaxBinPt, fMaxBinPt);
-    fHistDeltaPtRC[i]->GetXaxis()->SetTitle("#deltap_{T}^{RC} [GeV/c]");
+    fHistDeltaPtRC[i]->GetXaxis()->SetTitle("#delta#it{p}_{T}^{RC} (GeV/#it{c})");
     fHistDeltaPtRC[i]->GetYaxis()->SetTitle("counts");
     fOutput->Add(fHistDeltaPtRC[i]);
 
     histname = "fHistDeltaPtRCExLJ_";
     histname += i;
     fHistDeltaPtRCExLJ[i] = new TH1F(histname.Data(), histname.Data(), fNbins * 2, -fMaxBinPt, fMaxBinPt);
-    fHistDeltaPtRCExLJ[i]->GetXaxis()->SetTitle("#deltap_{T}^{RC} [GeV/c]");
+    fHistDeltaPtRCExLJ[i]->GetXaxis()->SetTitle("#delta#it{p}_{T}^{RC} (GeV/#it{c})");
     fHistDeltaPtRCExLJ[i]->GetYaxis()->SetTitle("counts");
     fOutput->Add(fHistDeltaPtRCExLJ[i]);
 
     histname = "fHistDeltaPtRCRand_";
     histname += i;
     fHistDeltaPtRCRand[i] = new TH1F(histname.Data(), histname.Data(), fNbins * 2, -fMaxBinPt, fMaxBinPt);
-    fHistDeltaPtRCRand[i]->GetXaxis()->SetTitle("#deltap_{T}^{RC} [GeV/c]");
+    fHistDeltaPtRCRand[i]->GetXaxis()->SetTitle("#delta#it{p}_{T}^{RC} (GeV/#it{c})");
     fHistDeltaPtRCRand[i]->GetYaxis()->SetTitle("counts");
     fOutput->Add(fHistDeltaPtRCRand[i]);
 
     if (!fEmbJetsName.IsNull()) {
-      histname = "fHistEmbJetsPt_";
+      histname = "fHistEmbJetsPtArea_";
       histname += i;
-      fHistEmbJetsPt[i] = new TH1F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt);
-      fHistEmbJetsPt[i]->GetXaxis()->SetTitle("embedded jet p_{T}^{raw} [GeV/c]");
-      fHistEmbJetsPt[i]->GetYaxis()->SetTitle("counts");
-      fOutput->Add(fHistEmbJetsPt[i]);
+      fHistEmbJetsPtArea[i] = new TH2F(histname.Data(), histname.Data(), 40, 0, fJetRadius * fJetRadius * TMath::Pi() * 3, fNbins, fMinBinPt, fMaxBinPt);
+      fHistEmbJetsPtArea[i]->GetXaxis()->SetTitle("area");
+      fHistEmbJetsPtArea[i]->GetYaxis()->SetTitle("embedded jet #it{p}_{T}^{raw} (GeV/#it{c})");
+      fOutput->Add(fHistEmbJetsPtArea[i]);
 
-      histname = "fHistEmbJetsCorrPt_";
+      histname = "fHistEmbJetsCorrPtArea_";
       histname += i;
-      fHistEmbJetsCorrPt[i] = new TH1F(histname.Data(), histname.Data(), fNbins * 2, -fMaxBinPt, fMaxBinPt);
-      fHistEmbJetsCorrPt[i]->GetXaxis()->SetTitle("embedded jet p_{T}^{corr} [GeV/c]");
-      fHistEmbJetsCorrPt[i]->GetYaxis()->SetTitle("counts");
-      fOutput->Add(fHistEmbJetsCorrPt[i]);
-
-      histname = "fHistEmbJetsArea_";
-      histname += i;
-      fHistEmbJetsArea[i] = new TH1F(histname.Data(), histname.Data(), 40, 0, fJetRadius * fJetRadius * TMath::Pi() * 3);
-      fHistEmbJetsArea[i]->GetXaxis()->SetTitle("area");
-      fHistEmbJetsArea[i]->GetYaxis()->SetTitle("counts");
-      fOutput->Add(fHistEmbJetsArea[i]);
+      fHistEmbJetsCorrPtArea[i] = new TH2F(histname.Data(), histname.Data(), 40, 0, fJetRadius * fJetRadius * TMath::Pi() * 3, fNbins * 2, -fMaxBinPt, fMaxBinPt);
+      fHistEmbJetsCorrPtArea[i]->GetXaxis()->SetTitle("area");
+      fHistEmbJetsCorrPtArea[i]->GetYaxis()->SetTitle("embedded jet #it{p}_{T}^{corr} (GeV/#it{c})");
+      fOutput->Add(fHistEmbJetsCorrPtArea[i]);
 
       histname = "fHistEmbPartPt_";
       histname += i;
       fHistEmbPartPt[i] = new TH1F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt);
-      fHistEmbPartPt[i]->GetXaxis()->SetTitle("embedded particle p_{T}^{emb} [GeV/c]");
+      fHistEmbPartPt[i]->GetXaxis()->SetTitle("embedded particle it{p}_{T}^{emb} (GeV/#it{c})");
       fHistEmbPartPt[i]->GetYaxis()->SetTitle("counts");
       fOutput->Add(fHistEmbPartPt[i]);
 
@@ -244,19 +237,26 @@ void AliAnalysisTaskDeltaPt::UserCreateOutputObjects()
       fHistEmbNotFoundPhiEta[i]->GetYaxis()->SetTitle("#phi");
       fOutput->Add(fHistEmbNotFoundPhiEta[i]);
 
+      histname = "fHistEmbBkgArea_";
+      histname += i;
+      fHistEmbBkgArea[i] = new TH2F(histname.Data(), histname.Data(), 40, 0, fJetRadius * fJetRadius * TMath::Pi() * 3, fNbins, fMinBinPt, fMaxBinPt);
+      fHistEmbBkgArea[i]->GetXaxis()->SetTitle("area");
+      fHistEmbBkgArea[i]->GetYaxis()->SetTitle("background of embedded track (GeV/#it{c})");
+      fOutput->Add(fHistEmbBkgArea[i]);
+
       histname = "fHistRhoVSEmbBkg_";
       histname += i;
-      fHistRhoVSEmbBkg[i] = new TH2F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt*2, fNbins, fMinBinPt, fMaxBinPt);
-      fHistRhoVSEmbBkg[i]->GetXaxis()->SetTitle("A#rho [GeV/c]");
-      fHistRhoVSEmbBkg[i]->GetYaxis()->SetTitle("background of embedded track [GeV/c]");
+      fHistRhoVSEmbBkg[i] = new TH2F(histname.Data(), histname.Data(), fNbins, fMinBinPt, fMaxBinPt, fNbins, fMinBinPt, fMaxBinPt);
+      fHistRhoVSEmbBkg[i]->GetXaxis()->SetTitle("A#rho (GeV/#it{c})");
+      fHistRhoVSEmbBkg[i]->GetYaxis()->SetTitle("background of embedded track (GeV/#it{c})");
       fOutput->Add(fHistRhoVSEmbBkg[i]);
       
-      histname = "fHistDeltaPtEmb_";
+      histname = "fHistDeltaPtEmbArea_";
       histname += i;
-      fHistDeltaPtEmb[i] = new TH1F(histname.Data(), histname.Data(), fNbins * 2, -fMaxBinPt, fMaxBinPt);
-      fHistDeltaPtEmb[i]->GetXaxis()->SetTitle("#deltap_{T}^{emb} [GeV/c]");
-      fHistDeltaPtEmb[i]->GetYaxis()->SetTitle("counts");
-      fOutput->Add(fHistDeltaPtEmb[i]);
+      fHistDeltaPtEmbArea[i] = new TH2F(histname.Data(), histname.Data(), 40, 0, fJetRadius * fJetRadius * TMath::Pi() * 3, fNbins * 2, -fMaxBinPt, fMaxBinPt);
+      fHistDeltaPtEmbArea[i]->GetXaxis()->SetTitle("area");
+      fHistDeltaPtEmbArea[i]->GetYaxis()->SetTitle("#delta#it{p}_{T}^{emb} (GeV/#it{c})");
+      fOutput->Add(fHistDeltaPtEmbArea[i]);
     }
   }
 
@@ -281,17 +281,17 @@ Bool_t AliAnalysisTaskDeltaPt::FillHistograms()
   
   const Float_t rcArea = fJetRadius * fJetRadius * TMath::Pi();
 
-  // Simple random cones
   for (Int_t i = 0; i < fRCperEvent; i++) {
+    // Simple random cones
     Float_t RCpt = 0;
     Float_t RCeta = 0;
     Float_t RCphi = 0;
     GetRandomCone(RCpt, RCeta, RCphi, 0);
     if (RCpt > 0) {
       fHistRCPhiEta->Fill(RCeta, RCphi);
-      fHistRhoVSRCPt[fCentBin]->Fill(fRhoVal, RCpt);
+      fHistRhoVSRCPt[fCentBin]->Fill(fRhoVal * rcArea, RCpt);
 
-      fHistRCPt[fCentBin]->Fill(RCpt / rcArea);
+      fHistRCPt[fCentBin]->Fill(RCpt);
       fHistDeltaPtRC[fCentBin]->Fill(RCpt - rcArea * fRhoVal);
     }
   
@@ -307,7 +307,7 @@ Bool_t AliAnalysisTaskDeltaPt::FillHistograms()
        if (dphi < -1.6) dphi += TMath::Pi() * 2; 
        fHistRCPtExLJVSDPhiLJ->Fill(RCpt, dphi);
       }
-      fHistRCPtExLJ[fCentBin]->Fill(RCpt / rcArea);
+      fHistRCPtExLJ[fCentBin]->Fill(RCpt);
       fHistDeltaPtRCExLJ[fCentBin]->Fill(RCpt - rcArea * fRhoVal);
     }
     
@@ -317,7 +317,7 @@ Bool_t AliAnalysisTaskDeltaPt::FillHistograms()
     RCphi = 0;
     GetRandomCone(RCpt, RCeta, RCphi, 0, fRandTracks, fRandCaloClusters);
     if (RCpt > 0) {
-      fHistRCPtRand[fCentBin]->Fill(RCpt / rcArea);
+      fHistRCPtRand[fCentBin]->Fill(RCpt);
       fHistDeltaPtRCRand[fCentBin]->Fill(RCpt - rcArea * fRhoVal);
     }  
   }
@@ -367,13 +367,14 @@ Bool_t AliAnalysisTaskDeltaPt::FillHistograms()
     
     fHistDistEmbPartJetAxis[fCentBin]->Fill(distProbeJet);
 
-    fHistEmbJetsPt[fCentBin]->Fill(embJet->Pt());
-    fHistEmbJetsCorrPt[fCentBin]->Fill(embJet->Pt() - fRhoVal * embJet->Area());
-    fHistEmbJetsArea[fCentBin]->Fill(embJet->Area());
+    fHistEmbJetsPtArea[fCentBin]->Fill(embJet->Area(), embJet->Pt());
+    fHistEmbJetsCorrPtArea[fCentBin]->Fill(embJet->Area(), embJet->Pt() - fRhoVal * embJet->Area());
     fHistEmbJetPhiEta->Fill(embJet->Eta(), embJet->Phi());
 
-    fHistDeltaPtEmb[fCentBin]->Fill(embJet->Pt() - embJet->Area() * fRhoVal - probePt);
-    fHistRhoVSEmbBkg[fCentBin]->Fill(fRhoVal, embJet->Pt() - probePt);
+    fHistEmbBkgArea[fCentBin]->Fill(embJet->Area(), embJet->Pt() - probePt);
+    fHistRhoVSEmbBkg[fCentBin]->Fill(fRhoVal * embJet->Area(), embJet->Pt() - probePt);
+    fHistDeltaPtEmbArea[fCentBin]->Fill(embJet->Area(), embJet->Pt() - embJet->Area() * fRhoVal - probePt);
+
   }
   else {
     if (fEmbTracks)
index 63a59bc..f4d36fd 100644 (file)
@@ -70,20 +70,20 @@ class AliAnalysisTaskDeltaPt : public AliAnalysisTaskEmcalJet {
 
   // Jet embedding
   TH2F                       *fHistEmbNotFoundPhiEta[4];   //!Phi-Eta of "not found" embedded particles
-  TH1F                       *fHistEmbJetsPt[4];           //!Pt distribution of embedded jets
-  TH1F                       *fHistEmbJetsCorrPt[4];       //!Pt-rho*A distribution of embedded jets
-  TH1F                       *fHistEmbJetsArea[4];         //!Area distribution of embedded jets
+  TH2F                       *fHistEmbJetsPtArea[4];       //!Pt vs. area of embedded jets
+  TH2F                       *fHistEmbJetsCorrPtArea[4];   //!Pt-rho*A vs. area of embedded jets
   TH1F                       *fHistEmbPartPt[4];           //!Pt distribution of embedded particle
   TH2F                       *fHistEmbJetPhiEta;           //!Phi-Eta distribution of embedded jets
   TH2F                       *fHistEmbPartPhiEta;          //!Phi-Eta distribution of embedded particles
   TH1F                       *fHistDistEmbPartJetAxis[4];  //!Distance between embedded particle and jet axis
+  TH2F                       *fHistEmbBkgArea[4];          //!Pt(embjet) - Pt(embtrack) vs. area of embedded jets
   TH2F                       *fHistRhoVSEmbBkg[4];         //!Area(embjet) * rho vs. Pt(embjet) - Pt(embtrack)
-  TH1F                       *fHistDeltaPtEmb[4];          //!deltaPt = Pt(embjet) - Area(embjet) * rho - Pt(embtrack)
+  TH2F                       *fHistDeltaPtEmbArea[4];      //!deltaPt = Pt(embjet) - Area(embjet) * rho - Pt(embtrack)
 
  private:
   AliAnalysisTaskDeltaPt(const AliAnalysisTaskDeltaPt&);            // not implemented
   AliAnalysisTaskDeltaPt &operator=(const AliAnalysisTaskDeltaPt&); // not implemented
 
-  ClassDef(AliAnalysisTaskDeltaPt, 1) // deltaPt analysis task
+  ClassDef(AliAnalysisTaskDeltaPt, 2) // deltaPt analysis task
 };
 #endif
index 2855b7e..9bee01c 100644 (file)
@@ -33,6 +33,7 @@ AliAnalysisTaskRhoBase::AliAnalysisTaskRhoBase() :
   fCompareRhoScaled(0),
   fHistJetPtvsCent(0),
   fHistJetAreavsCent(0),
+  fHistJetRhovsCent(0),
   fHistNjetvsCent(0),
   fHistJetPtvsNtrack(0),
   fHistJetAreavsNtrack(0),
@@ -73,6 +74,7 @@ AliAnalysisTaskRhoBase::AliAnalysisTaskRhoBase(const char *name, Bool_t histo) :
   fCompareRhoScaled(0),
   fHistJetPtvsCent(0),
   fHistJetAreavsCent(0),
+  fHistJetRhovsCent(0),
   fHistNjetvsCent(0),
   fHistJetPtvsNtrack(0),
   fHistJetAreavsNtrack(0),
@@ -126,6 +128,7 @@ void AliAnalysisTaskRhoBase::UserCreateOutputObjects()
   if (!fJetsName.IsNull()) {
     fHistJetPtvsCent            = new TH2F("JetPtvsCent",           "JetPtvsCent",           101, -1,  100,   fNbins, fMinBinPt, fMaxBinPt);
     fHistJetAreavsCent          = new TH2F("JetAreavsCent",         "JetAreavsCent",         101, -1,  100,   30, 0, fJetRadius * fJetRadius * TMath::Pi() * 3);
+    fHistJetRhovsCent           = new TH2F("fHistJetRhovsCent",     "fHistJetRhovsCent",     101, -1,  100,   fNbins, fMinBinPt, fMaxBinPt*2);
     fHistNjetvsCent             = new TH2F("NjetvsCent",            "NjetvsCent",            101, -1,  100,   150, -0.5, 149.5);
 
     fOutput->Add(fHistJetPtvsCent);
@@ -248,6 +251,7 @@ Bool_t AliAnalysisTaskRhoBase::FillHistograms()
       
       fHistJetPtvsCent->Fill(fCent, jet->Pt());
       fHistJetAreavsCent->Fill(fCent, jet->Area());
+      fHistJetRhovsCent->Fill(fCent, jet->Pt() / jet->Area());
       
       if (fTracks) {
        fHistJetPtvsNtrack->Fill(Ntracks, jet->Pt());
index f217023..74a023d 100644 (file)
@@ -50,6 +50,7 @@ class AliAnalysisTaskRhoBase : public AliAnalysisTaskEmcalJet {
 
   TH2F                  *fHistJetPtvsCent;               //!jet pt vs. centrality
   TH2F                  *fHistJetAreavsCent;             //!jet area vs. centrality
+  TH2F                  *fHistJetRhovsCent;              //!jet pt/area vs. centrality
   TH2F                  *fHistNjetvsCent;                //!no. of jets vs. centrality
   TH2F                  *fHistJetPtvsNtrack;             //!jet pt vs. no. of tracks
   TH2F                  *fHistJetAreavsNtrack;           //!jet area vs. no. of tracks
@@ -73,6 +74,6 @@ class AliAnalysisTaskRhoBase : public AliAnalysisTaskEmcalJet {
   AliAnalysisTaskRhoBase(const AliAnalysisTaskRhoBase&);             // not implemented
   AliAnalysisTaskRhoBase& operator=(const AliAnalysisTaskRhoBase&);  // not implemented
   
-  ClassDef(AliAnalysisTaskRhoBase, 7); // Rho base task
+  ClassDef(AliAnalysisTaskRhoBase, 8); // Rho base task
 };
 #endif
index e774cc9..5b88e7a 100644 (file)
@@ -137,6 +137,7 @@ void AliEmcalJetTask::FindJets()
   // setup fj wrapper
   AliFJWrapper fjw(name, name);
   fjw.SetAreaType(fastjet::active_area_explicit_ghosts);
+  fjw.SetGhostArea(fGhostArea);
   fjw.SetR(fRadius);
   fjw.SetAlgorithm(jalgo);  
   fjw.SetMaxRap(1);
index b65a3ae..6c9a433 100644 (file)
@@ -30,6 +30,7 @@ class AliEmcalJetTask : public AliAnalysisTaskSE {
   void                   SetType(Int_t t)                 { fType          = t     ; }
   void                   SetEtaRange(Double_t emi, Double_t ema) {fEtaMin = emi; fEtaMax = ema; }
   void                   SetPhiRange(Double_t pmi, Double_t pma) {fPhiMin = pmi; fPhiMax = pma; }
+  void                   SetGhostArea(Double_t gharea)    { fGhostArea      = gharea;  }
 
  protected:
   void                   FindJets();
@@ -49,6 +50,7 @@ class AliEmcalJetTask : public AliAnalysisTaskSE {
   Double_t               fEtaMax;                 // maximum eta for constituents (applied before clustering)
   Double_t               fMinJetArea;             // min area to keep jet in output
   Double_t               fMinJetPt;               // min jet pt to keep jet in output
+  Double_t               fGhostArea;              // ghost area
   Bool_t                 fIsInit;                 //!=true if already initialized
   Bool_t                 fIsMcPart;               //!=true if MC particles are given as input
   Bool_t                 fIsEmcPart;              //!=true if emcal particles are given as input (for clusters)
@@ -61,6 +63,6 @@ class AliEmcalJetTask : public AliAnalysisTaskSE {
   AliEmcalJetTask(const AliEmcalJetTask&);            // not implemented
   AliEmcalJetTask &operator=(const AliEmcalJetTask&); // not implemented
 
-  ClassDef(AliEmcalJetTask, 4) // Jet producing task
+  ClassDef(AliEmcalJetTask, 5) // Jet producing task
 };
 #endif
index 8dc65b5..bfb5db7 100644 (file)
@@ -34,6 +34,7 @@ ClassImp(AliAnalysisTaskSAQA)
 AliAnalysisTaskSAQA::AliAnalysisTaskSAQA() : 
   AliAnalysisTaskEmcalJet("AliAnalysisTaskSAQA", kTRUE),
   fCellEnergyCut(0.1),
+  fParticleLevel(kFALSE),
   fNclusters(0),
   fNtracks(0),
   fNjets(0),
@@ -72,6 +73,7 @@ AliAnalysisTaskSAQA::AliAnalysisTaskSAQA() :
 AliAnalysisTaskSAQA::AliAnalysisTaskSAQA(const char *name) : 
   AliAnalysisTaskEmcalJet(name, kTRUE),
   fCellEnergyCut(0.1),
+  fParticleLevel(kFALSE),
   fNclusters(0),
   fNtracks(0),
   fNjets(0),
index aade385..d45ce0e 100644 (file)
@@ -7,7 +7,8 @@ AliEmcalJetTask* AddTaskEmcalJet(
   const Double_t radius      = 0.4,
   const Int_t type           = 0,
   const Double_t minTrPt     = 0.15,
-  const Double_t minClPt     = 0.15
+  const Double_t minClPt     = 0.30,
+  const Double_t ghostArea   = 0.01  
 )
 {  
   // Get the pointer to the existing analysis manager via the static access method.
@@ -94,6 +95,7 @@ AliEmcalJetTask* AddTaskEmcalJet(
   jetTask->SetMinJetPt(minJetPt);
   jetTask->SetRadius(radius);
   jetTask->SetType(type);
+  jetTask->SetGhostArea(ghostArea);
 
   //-------------------------------------------------------
   // Final settings, pass to manager and set the containers
index e37a1d8..5123657 100644 (file)
@@ -3,8 +3,8 @@
 AliJetResponseMaker* AddTaskJetRespPtHard(const char *ntracks            = "Tracks",
                                          const char *nclusters          = "CaloClusters",
                                          const char *njets              = "Jets",
-                                         const char *nmcjets            = "MCJets",
                                          const char *nmctracks          = "MCParticles",
+                                         const char *nmcjets            = "MCJets",
                                          Double_t    jetradius          = 0.2,
                                          Double_t    jetptcut           = 1,
                                          Double_t    jetareacut         = 0.557,
@@ -25,7 +25,7 @@ AliJetResponseMaker* AddTaskJetRespPtHard(const char *ntracks            = "Trac
   AliJetResponseMaker *jetTask = new AliJetResponseMaker[maxPtBin - minPtBin + 1];
 
   for (Int_t i = minPtBin; i <= maxPtBin; i++) {
-    AddTaskJetResponseMaker(ntracks, nclusters, njets, nmcjets, nmctracks, 
+    AddTaskJetResponseMaker(ntracks, nclusters, njets, nmctracks, nmcjets,
                            jetradius, jetptcut, jetareacut, jetBiasTrack, 
                            jetBiasClus, maxDistance, type, i, taskname, jetTask + i - minPtBin);
     jetTask[i - minPtBin].SetDoMatching(domatch);
index 6aa2d82..7643af3 100644 (file)
@@ -34,13 +34,22 @@ AliAnalysisTaskSAJF* AddTaskSAJF(
   // Init the task and do settings
   //-------------------------------------------------------
 
-  TString name(Form("%s_%s_%s_R0%d_",taskname,njets,nrho,(Int_t)floor(jetradius*100+0.5)));
+  TString name(taskname);
+  if (strcmp(njets,"")) {
+    name += "_";
+    name += njets;
+  }
+  if (strcmp(nrho,"")) {
+    name += "_";
+    name += nrho;
+  }
   if (type == AliAnalysisTaskEmcal::kTPC) 
-    name += "TPC";
+    name += "_TPC";
   else if (type == AliAnalysisTaskEmcal::kEMCAL) 
-    name += "EMCAL";
+    name += "_EMCAL";
   else if (type == AliAnalysisTaskEmcal::kUser) 
-    name += "USER";
+    name += "_USER";
+
   AliAnalysisTaskSAJF* jetTask = new AliAnalysisTaskSAJF(name);
   jetTask->SetAnaType(type);
   jetTask->SetTracksName(ntracks);
index 57f5636..82238bf 100644 (file)
@@ -34,19 +34,25 @@ AliAnalysisTaskSAQA* AddTaskSAQA(
   // Init the task and do settings
   //-------------------------------------------------------
   
-  TString name;
-  if (strcmp(njets,"")) {
-    name = Form("%s_%s_%s_%s_R0%d_",taskname,njets,ntracks,nclusters,(Int_t)floor(jetradius*100+0.5));
-    if (type == AliAnalysisTaskEmcal::kTPC) 
-      name += "TPC";
-    else if (type == AliAnalysisTaskEmcal::kEMCAL) 
-      name += "EMCAL";
-    else if (type == AliAnalysisTaskEmcal::kUser) 
-      name += "USER";
+  TString name(taskname);
+  if (strcmp(ntracks,"")) {
+    name += "_";
+    name += ntracks;
+  }
+  if (strcmp(nclusters,"")) {
+    name += "_";
+    name += nclusters;
   }
-  else {
-    name = Form("%s_%s_%s_",taskname,ntracks,nclusters);
+  if (strcmp(njets,"")) {
+    name += "_";
+    name += njets;
   }
+  if (type == AliAnalysisTaskEmcal::kTPC) 
+    name += "_TPC";
+  else if (type == AliAnalysisTaskEmcal::kEMCAL) 
+    name += "_EMCAL";
+  else if (type == AliAnalysisTaskEmcal::kUser) 
+    name += "_USER";
 
   AliAnalysisTaskSAQA* qaTask = new AliAnalysisTaskSAQA(name);
   qaTask->SetTracksName(ntracks);