Adding copy constructor (Marian)

[u/mrichter/AliRoot.git] / PWG2 / AliAnalysisTaskProtonsQA.cxx

diff --git a/PWG2/AliAnalysisTaskProtonsQA.cxx b/PWG2/AliAnalysisTaskProtonsQA.cxx
index 65c54ca4ef213ccad1eb476761e6de67ad5283ed..18034be7c41f54433e84753cb39cd470de0c55ef 100644 (file)
--- a/PWG2/AliAnalysisTaskProtonsQA.cxx
+++ b/PWG2/AliAnalysisTaskProtonsQA.cxx
@@ -17,7 +17,7 @@
 #include "AliMCEvent.h"
 #include "AliStack.h"
 
-#include "PWG2spectra/SPECTRA/AliProtonQAAnalysis.h"
+#include "AliProtonQAAnalysis.h"
 #include "AliAnalysisTaskProtonsQA.h"
 
 // Analysis task used for the QA of the (anti)proton analysis
@@ -28,7 +28,7 @@ ClassImp(AliAnalysisTaskProtonsQA)
 //________________________________________________________________________ 
 AliAnalysisTaskProtonsQA::AliAnalysisTaskProtonsQA()
   : AliAnalysisTask(), fESD(0), fMC(0),
-    fList0(0), fList1(0), fList2(0), fList3(0), fList4(0),
+    fList0(0), fList1(0), fList2(0), fList3(0), fList4(0), fList5(0),
     fAnalysis(0) {
     //Dummy constructor
 }
@@ -36,7 +36,7 @@ AliAnalysisTaskProtonsQA::AliAnalysisTaskProtonsQA()
 //________________________________________________________________________
 AliAnalysisTaskProtonsQA::AliAnalysisTaskProtonsQA(const char *name) 
 : AliAnalysisTask(name, ""), fESD(0), fMC(0),
-  fList0(0), fList1(0), fList2(0), fList3(0), fList4(0),
+  fList0(0), fList1(0), fList2(0), fList3(0), fList4(0), fList5(0),
   fAnalysis(0) {
   // Constructor
 
@@ -49,6 +49,7 @@ AliAnalysisTaskProtonsQA::AliAnalysisTaskProtonsQA(const char *name)
   DefineOutput(2, TList::Class());
   DefineOutput(3, TList::Class());
   DefineOutput(4, TList::Class());
+  DefineOutput(5, TList::Class());
 }
 
 //________________________________________________________________________
@@ -85,50 +86,72 @@ void AliAnalysisTaskProtonsQA::CreateOutputObjects() {
   
   //proton analysis object
   fAnalysis = new AliProtonQAAnalysis();
-  fAnalysis->SetQAOn();
   fAnalysis->SetRunMCAnalysis();
-  //fAnalysis->SetMCProcessId(4);//4: weak decay - 13: hadronic interaction
+  fAnalysis->SetRunEfficiencyAnalysis(kTRUE,kTRUE); //kTRUE,kTRUE for eta-pT efficiencies and if the cuts should be used in the reco and pid efficiencies
+  //fAnalysis->SetMCProcessId(13);//4: weak decay - 13: hadronic interaction
   //fAnalysis->SetMotherParticlePDGCode(3122);//3122: Lambda
 
   //Use of TPConly tracks
-  fAnalysis->SetQAYPtBins(10, -0.5, 0.5, 12, 0.5, 0.9); //TPC only
+  /*fAnalysis->SetQAYPtBins(10, -0.5, 0.5, 12, 0.5, 0.9); //TPC only
   fAnalysis->UseTPCOnly();
-  //fAnalysis->SetMinTPCClusters(100);
-  //fAnalysis->SetMaxChi2PerTPCCluster(2.1);
-  //fAnalysis->SetMaxCov11(1.0);
-  //fAnalysis->SetMaxCov22(1.0);
-  //fAnalysis->SetMaxCov33(0.2);
-  //fAnalysis->SetMaxCov44(0.2);
-  //fAnalysis->SetMaxCov55(0.5);
+  fAnalysis->SetTPCpid();
+  fAnalysis->SetMinTPCClusters(100);
+  fAnalysis->SetMaxChi2PerTPCCluster(2.2);
+  fAnalysis->SetMaxCov11(0.5);
+  fAnalysis->SetMaxCov22(0.5);
+  fAnalysis->SetMaxCov33(0.5);
+  fAnalysis->SetMaxCov44(0.5);
+  fAnalysis->SetMaxCov55(0.5);
   //fAnalysis->SetMaxSigmaToVertexTPC(2.0);
-  //fAnalysis->SetMaxDCAXY(1.0);
-  //fAnalysis->SetMaxDCAZ(1.0);
-  //fAnalysis->SetTPCRefit();
+  fAnalysis->SetMaxDCAXYTPC(1.5);
+  fAnalysis->SetMaxDCAZTPC(1.5);*/
+
+  //Use of HybridTPC tracks
+  fAnalysis->SetQAYPtBins(20, -1.0, 1.0, 26, 0.2, 1.5); //HybridTPC
+  fAnalysis->UseHybridTPC();
   fAnalysis->SetTPCpid();
+  fAnalysis->SetMinTPCClusters(110);
+  fAnalysis->SetMaxChi2PerTPCCluster(2.2);
+  fAnalysis->SetMaxCov11(0.5);
+  fAnalysis->SetMaxCov22(0.5);
+  fAnalysis->SetMaxCov33(0.5);
+  fAnalysis->SetMaxCov44(0.5);
+  fAnalysis->SetMaxCov55(0.5);
+  fAnalysis->SetMaxSigmaToVertex(2.0);
+  //fAnalysis->SetMaxDCAXY(1.5);
+  //fAnalysis->SetMaxDCAZ(1.5);
+  //fAnalysis->SetPointOnITSLayer6();
+  //fAnalysis->SetPointOnITSLayer5();
+  //fAnalysis->SetPointOnITSLayer4();
+  //fAnalysis->SetPointOnITSLayer3();
+  fAnalysis->SetPointOnITSLayer2();
+  fAnalysis->SetPointOnITSLayer1();
+  fAnalysis->SetMinITSClusters(5);
 
   //Combined tracking
+  /*fAnalysis->SetQAYPtBins(10, -0.5, 0.5, 12, 0.5, 0.9); //HybridTPC
   //fAnalysis->SetQAYPtBins(20, -1.0, 1.0, 27, 0.4, 3.1); //combined tracking
-  //fAnalysis->SetMinTPCClusters(50);
-  //fAnalysis->SetMaxChi2PerTPCCluster(3.5);
-  //fAnalysis->SetMaxCov11(2.0);
-  //fAnalysis->SetMaxCov22(2.0);
-  //fAnalysis->SetMaxCov33(0.5);
-  //fAnalysis->SetMaxCov44(0.5);
-  //fAnalysis->SetMaxCov55(2.0);
-  //fAnalysis->SetMaxSigmaToVertex(2.0);
-  //fAnalysis->SetTPCRefit();
-  //ITS related cuts - to be used in the case of the analysis of global tracks
-  //fAnalysis->SetPointOnITSLayer1();
-  //fAnalysis->SetPointOnITSLayer2();
+  fAnalysis->SetMinTPCClusters(110);
+  fAnalysis->SetMaxChi2PerTPCCluster(2.2);
+  fAnalysis->SetMaxCov11(0.5);
+  fAnalysis->SetMaxCov22(0.5);
+  fAnalysis->SetMaxCov33(0.5);
+  fAnalysis->SetMaxCov44(0.5);
+  fAnalysis->SetMaxCov55(0.5);
+  fAnalysis->SetMaxSigmaToVertex(2.0);
+  //fAnalysis->SetMaxDCAXY(2.0);
+  //fAnalysis->SetMaxDCAZ(2.0);
+  fAnalysis->SetTPCRefit();
+  fAnalysis->SetPointOnITSLayer1();
+  fAnalysis->SetPointOnITSLayer2();
   //fAnalysis->SetPointOnITSLayer3();
   //fAnalysis->SetPointOnITSLayer4();
   //fAnalysis->SetPointOnITSLayer5();
   //fAnalysis->SetPointOnITSLayer6();
-  //fAnalysis->SetMinITSClusters(1);
-  //fAnalysis->SetITSRefit();
-  //fAnalysis->SetESDpid();
+  fAnalysis->SetMinITSClusters(1);
+  fAnalysis->SetITSRefit();
+  fAnalysis->SetESDpid();*/
 
-  fAnalysis->InitQA();
   fAnalysis->SetPriorProbabilities(partFrac);
 
   fList0 = new TList();
@@ -145,6 +168,9 @@ void AliAnalysisTaskProtonsQA::CreateOutputObjects() {
 
   fList4 = new TList();
   fList4 = fAnalysis->GetAcceptedDCAList();
+
+  fList5 = new TList();
+  fList5 = fAnalysis->GetEfficiencyQAList();
 }
 
 //________________________________________________________________________
@@ -168,8 +194,9 @@ void AliAnalysisTaskProtonsQA::Exec(Option_t *) {
     return;
   }
   
-  fAnalysis->RunQA(stack, fESD);
+  fAnalysis->RunQAAnalysis(stack, fESD);
   fAnalysis->RunMCAnalysis(stack);
+  fAnalysis->RunEfficiencyAnalysis(stack, fESD);
 
   // Post output data.
   PostData(0, fList0);
@@ -177,6 +204,7 @@ void AliAnalysisTaskProtonsQA::Exec(Option_t *) {
   PostData(2, fList2);
   PostData(3, fList3);
   PostData(4, fList4);
+  PostData(5, fList5);
 }      
 
 //________________________________________________________________________