allow for jet by jet comparisio

[u/mrichter/AliRoot.git] / ANALYSIS / AliAnalysisTaskPIDCombined.cxx

diff --git a/ANALYSIS/AliAnalysisTaskPIDCombined.cxx b/ANALYSIS/AliAnalysisTaskPIDCombined.cxx
index c1ec62804076b332991b783123ccdbd9c91143b4..465b1cec426fe9d785e73ed2fd0be2d70ce0408f 100644 (file)
--- a/ANALYSIS/AliAnalysisTaskPIDCombined.cxx
+++ b/ANALYSIS/AliAnalysisTaskPIDCombined.cxx
@@ -50,10 +50,23 @@ ClassImp(AliAnalysisTaskPIDCombined)
 AliAnalysisTaskPIDCombined::AliAnalysisTaskPIDCombined() :\r
   AliAnalysisTaskSE(),\r
   fHistList(),\r
+  fProbTPCnSigma(),\r
+  fProbTOFnSigma(),\r
+  fProbTPCTOFnSigmaTPC(),\r
+  fProbTPC(),\r
+  fProbTOF(),\r
+  fProbTPCTOF(),\r
+  fPriors(),\r
+  fProbTPCTOFnSigTPCMom(),\r
+  fProbTPCnSigTPCMom(),\r
+  fProbTOFnSigTOFMom(),\r
+  fPriorsUsed(),\r
   fPIDResponse(0x0),\r
   fPIDCombined(0x0),\r
   fTrackCuts(0x0),\r
-  fTrackFilter(0x0)\r
+  fTrackFilter(0x0),\r
+  fDeDx(NULL),
+  fDeDxTuned(NULL)
 {\r
   //\r
   // Constructor\r
@@ -64,10 +77,23 @@ AliAnalysisTaskPIDCombined::AliAnalysisTaskPIDCombined() :
 AliAnalysisTaskPIDCombined::AliAnalysisTaskPIDCombined(const char *name) :\r
   AliAnalysisTaskSE(name),\r
   fHistList(),\r
+  fProbTPCnSigma(),\r
+  fProbTOFnSigma(),\r
+  fProbTPCTOFnSigmaTPC(),\r
+  fProbTPC(),\r
+  fProbTOF(),\r
+  fProbTPCTOF(),\r
+  fPriors(),\r
+  fProbTPCTOFnSigTPCMom(),\r
+  fProbTPCnSigTPCMom(),\r
+  fProbTOFnSigTOFMom(),\r
+  fPriorsUsed(),\r
   fPIDResponse(0x0),\r
   fPIDCombined(0x0),\r
   fTrackCuts(0x0),\r
-  fTrackFilter(0x0)\r
+  fTrackFilter(0x0),\r
+  fDeDx(NULL),
+  fDeDxTuned(NULL)
 {\r
   //\r
   // Constructor\r
@@ -104,13 +130,13 @@ void AliAnalysisTaskPIDCombined::UserCreateOutputObjects()
   fPIDCombined->SetDefaultTPCPriors();\r
   fPIDCombined->SetDetectorMask(AliPIDResponse::kDetTPC+AliPIDResponse::kDetTOF);\r
 \r
-  // no light nuclei\r
-  fPIDCombined->SetSelectedSpecies(AliPID::kSPECIES);\r
+  // no light nuclei - no need to call it, this is default\r
+  //  fPIDCombined->SetSelectedSpecies(AliPID::kSPECIES);\r
 \r
 \r
   fHistList.Add(new TH1D("nEvents","Number of Evnets;Selection",2,0,2));\r
 \r
-  for (Int_t ispec=0; ispec<AliPID::kSPECIES; ++ispec){\r
+  for (Int_t ispec=0; ispec<5; ++ispec){\r
 \r
  \r
     fProbTPC[ispec]=new TH2D(Form("prob%s_mom_TPC",AliPID::ParticleName(ispec)),\r
@@ -199,7 +225,11 @@ void AliAnalysisTaskPIDCombined::UserCreateOutputObjects()
   }\r
 \r
 \r
-\r
+  fDeDx = new TH2D("hDeDx",";p_{TPC};dE/dx (a.u.)",500,0,5,500,0,500);
+  fHistList.Add(fDeDx);
+  fDeDxTuned = new TH2D("hDeDxTuned",";p_{TPC};dE/dx (a.u.)",500,0,5,500,0,500);
+  fHistList.Add(fDeDxTuned);
+
   fHistList.SetOwner();\r
   PostData(1,&fHistList);\r
 \r
@@ -279,6 +309,10 @@ void AliAnalysisTaskPIDCombined::UserExec(Option_t *)
 \r
       }\r
 \r
+      fPIDResponse->GetTPCsignalTunedOnData(track);
+
+      fDeDx->Fill(mom,track->GetTPCsignal());
+      fDeDxTuned->Fill(mom,track->GetTPCsignalTunedOnData());
 \r
     }\r
   }\r