QA plots for q vector

diff --git a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliAnalysisTaskV2AllChAOD.cxx b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliAnalysisTaskV2AllChAOD.cxx
index 487514750102a58865551339713b580cc7365cd3..dbcb71f447ff9a275ae4ae7b7514f98ab07693ff 100644 (file)
--- a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliAnalysisTaskV2AllChAOD.cxx
+++ b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliAnalysisTaskV2AllChAOD.cxx
@@ -76,6 +76,9 @@ AliAnalysisTaskV2AllChAOD::AliAnalysisTaskV2AllChAOD(const char *name) : AliAnal
   fMinTPCNcls(70),
   fFillTHn(kFALSE),
   fCentrality(0),
+  fQvector(0),
+  fQvector_lq(0),
+  fQvector_sq(0),
   fResSP(0),
   fResSP_vs_Cent(0),
   f2partCumQA_vs_Cent(0),
@@ -224,6 +227,15 @@ void AliAnalysisTaskV2AllChAOD::UserCreateOutputObjects()
   fCentrality = new TH1D("fCentrality", "centrality distribution; centrality", 200, 0., 100);
   fOutput->Add(fCentrality);
   
+  fQvector = new TH1D("fQvector", "q-vector distribution; q-vector", fnQvecBins, 0., fQvecUpperLim);
+  fOutput->Add(fQvector);
+  
+  fQvector_lq = new TH1D("fQvector_lq", "q-vector distribution; q-vector", fnQvecBins, 0., fQvecUpperLim);
+  fOutput_lq->Add(fQvector_lq);
+  
+  fQvector_sq = new TH1D("fQvector_sq", "q-vector distribution; q-vector", fnQvecBins, 0., fQvecUpperLim);
+  fOutput_sq->Add(fQvector_sq);
+  
   // binning common to all the THn
   //change it according to your needs + move it to global variables -> setter/getter
 //   Double_t ptBins[] = {0., 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.2, 1.4, 1.6, 1.8, 2.0, 2.4, 2.8, 3.2, 3.6, 4.0, 4.5, 5.0, 5.5, 6.0, 7.0, 8.0, 9.0, 10.0, 12.0, 14.0, 16.0, 20.0};
@@ -458,7 +470,10 @@ void AliAnalysisTaskV2AllChAOD::UserExec(Option_t *)
   if(fIsMC && fQvecGen) Qvec = fEventCuts->GetQvecPercentileMC(fVZEROside, fQgenType);
   else Qvec = fEventCuts->GetQvecPercentile(fVZEROside);
 
-  
+  fQvector->Fill(Qvec);
+  if (Qvec > fCutLargeQperc && Qvec < 100.) fQvector_lq->Fill(Qvec);
+  if (Qvec > 0. && Qvec < fCutSmallQperc) fQvector_sq->Fill(Qvec);
+
   Double_t Cent=(fDoCentrSystCentrality)?1.01*fEventCuts->GetCent():fEventCuts->GetCent();
   fCentrality->Fill(Cent);
   

diff --git a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliAnalysisTaskV2AllChAOD.h b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliAnalysisTaskV2AllChAOD.h
index 05b3f519ac1e59529ef15d859f8def8c4e800d70..978dc6666b1c5709922b6d32ffedc76657855394 100644 (file)
--- a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliAnalysisTaskV2AllChAOD.h
+++ b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliAnalysisTaskV2AllChAOD.h
@@ -49,6 +49,9 @@ class AliAnalysisTaskV2AllChAOD : public AliAnalysisTaskSE
     fMinTPCNcls(70),
     fFillTHn(kTRUE),
     fCentrality(0),
+    fQvector(0),
+    fQvector_lq(0),
+    fQvector_sq(0),
     fResSP(0),
     fResSP_vs_Cent(0),
     f2partCumQA_vs_Cent(0),
@@ -175,7 +178,10 @@ class AliAnalysisTaskV2AllChAOD : public AliAnalysisTaskSE
   Bool_t fFillTHn;
   
   TH1D * fCentrality;
-  
+  TH1D * fQvector;
+  TH1D * fQvector_lq;
+  TH1D * fQvector_sq;
+
   //output object
   TProfile*     fResSP;             //! resolution
   TProfile*     fResSP_vs_Cent;
@@ -260,7 +266,7 @@ class AliAnalysisTaskV2AllChAOD : public AliAnalysisTaskSE
   AliAnalysisTaskV2AllChAOD(const AliAnalysisTaskV2AllChAOD&);
   AliAnalysisTaskV2AllChAOD& operator=(const AliAnalysisTaskV2AllChAOD&);
   
-  ClassDef(AliAnalysisTaskV2AllChAOD, 12);
+  ClassDef(AliAnalysisTaskV2AllChAOD, 13);
 };
 
 #endif