From d38f1e8800ce223867f77b96404baf91a9287675 Mon Sep 17 00:00:00 2001
From: jotwinow <jotwinow@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Tue, 14 Jul 2009 17:47:06 +0000
Subject: [PATCH] QA for the on-the-fly V0

---
 PWG1/AliAnalysisTaskV0QA.h | 177 +++++++++++++++++++++++++++++++++++++
 1 file changed, 177 insertions(+)
 create mode 100644 PWG1/AliAnalysisTaskV0QA.h

diff --git a/PWG1/AliAnalysisTaskV0QA.h b/PWG1/AliAnalysisTaskV0QA.h
new file mode 100644
index 00000000000..0005510f304
--- /dev/null
+++ b/PWG1/AliAnalysisTaskV0QA.h
@@ -0,0 +1,177 @@
+#include "TH1.h"
+#include "THnSparse.h"
+#include "TList.h"
+#include "AliESDEvent.h"
+#include "AliESD.h"
+#include "TLorentzVector.h"
+
+#include "AliAnalysisTask.h"
+#include "AliAnalysisManager.h"
+#include "AliAnalysisDataContainer.h"
+#include "AliESDtrackCuts.h"
+class AliESDtrackCuts;
+class TChain;
+
+class AliAnalysisTaskV0QA : public AliAnalysisTask {
+ public:
+  //  AliAnalysisTaskV0QA() : AliAnalysisTask(), fESD(0), fChain(0) {}
+  AliAnalysisTaskV0QA(const char *name);
+  virtual ~AliAnalysisTaskV0QA();
+  
+  virtual void   ConnectInputData(Option_t *);
+  virtual void   CreateOutputObjects();
+  virtual void   Exec(Option_t *option);
+  virtual void   Terminate(Option_t *);
+
+  void InspectListOfChargedParticles();
+  void InspectListOfV0s();
+  void FillHnSparseGamma();
+  void FillHnSparseK0();
+  void FillHnSparseL();
+  void FillHnSparseAL();
+
+  //  void getPID(AliESDtrack *esdTrack, Stat_t &fpid, Stat_t &fweight);
+
+
+  //  AliESDtrackCuts* fEsdTrackCuts;           // Object containing the parameters of the esd track cuts
+  //  void SetESDtrackCuts();
+
+
+
+
+  Int_t GetTPCReference(Int_t label);
+
+
+
+
+ private:
+
+
+
+  AliESDEvent *fESD; //ESD object
+  AliStack * stack;
+  AliMCEventHandler* mctruth; 
+ 
+  TChain * fChain;
+  TList * fOutputContainer; // ! output data container
+
+  THnSparse   *fSparseV0;
+  THnSparse   *fSparseK0;
+  THnSparse   *fSparseL;
+  THnSparse   *fSparseAL;
+
+  //////////////////////////////////
+
+  Int_t nEv;
+
+  Int_t nConvGamGeant;
+
+  Int_t * gConvGamGeantIndex;
+  Int_t * eNegConvGamGeantIndex;
+  Int_t * ePosConvGamGeantIndex;
+  Float_t * eNegConvGamGeantLength;
+  Float_t * ePosConvGamGeantLength;
+
+
+  Int_t * eNegConvGamSingleRecIndex;
+  Int_t * ePosConvGamSingleRecIndex;
+  Int_t * eNegConvGamV0RecIndex;
+  Int_t * ePosConvGamV0RecIndex;
+  Int_t * ConvGamV0RecIndexPos;
+  Int_t * ConvGamV0RecIndexNeg;
+
+
+  Int_t gDim;
+  // Lambda
+  Int_t nDecayLGeant;
+  Int_t * lDecayLGeantIndex;
+  Int_t * piNegDecayLGeantIndex;
+  Int_t * pPosDecayLGeantIndex;
+  Float_t * piNegDecayLGeantLength;
+  Float_t * pPosDecayLGeantLength;
+
+  Int_t * piNegDecayLSingleRecIndex;
+  Int_t * pPosDecayLSingleRecIndex;
+  Int_t * piNegDecayLV0RecIndex;
+  Int_t * pPosDecayLV0RecIndex;
+  Int_t * DecayLV0RecIndexPos;
+  Int_t * DecayLV0RecIndexNeg;
+
+
+  // AntiLambda
+  Int_t nDecayALGeant;
+  Int_t * alDecayALGeantIndex;
+  Int_t * piPosDecayALGeantIndex;
+  Int_t * apNegDecayALGeantIndex;
+  Float_t * piPosDecayALGeantLength;
+  Float_t * apNegDecayALGeantLength;
+
+  Int_t * piPosDecayALSingleRecIndex;
+  Int_t * apNegDecayALSingleRecIndex;
+  Int_t * piPosDecayALV0RecIndex;
+  Int_t * apNegDecayALV0RecIndex;
+  Int_t * DecayALV0RecIndexPos;
+  Int_t * DecayALV0RecIndexNeg;
+
+
+  // K0S
+  Int_t nDecayK0Geant;
+  Int_t * K0DecayK0GeantIndex;
+  Int_t * piNegDecayK0GeantIndex;
+  Int_t * piPosDecayK0GeantIndex;
+  Float_t * piNegDecayK0GeantLength;
+  Float_t * piPosDecayK0GeantLength;
+
+  Int_t * piNegDecayK0SingleRecIndex;
+  Int_t * piPosDecayK0SingleRecIndex;
+  Int_t * piNegDecayK0V0RecIndex;
+  Int_t * piPosDecayK0V0RecIndex;
+  Int_t * DecayK0V0RecIndexPos;
+  Int_t * DecayK0V0RecIndexNeg;
+
+  Int_t piPosK0Index;
+  Int_t piNegK0Index;
+
+  
+  Int_t nTracksPrim;
+
+
+  Int_t tpcRefit;
+  Int_t itsRefit;
+  Int_t trdRefit;
+  Int_t trdOut;
+
+
+
+
+  Double_t probabilityPos[AliPID::kSPECIES];
+  Double_t probabilityNeg[AliPID::kSPECIES];
+
+  Double_t * fValueL;
+  Double_t * fValueAL;
+  Double_t * fValueK0;
+  Double_t * fValueV0;
+
+  Double_t * xminV0;
+  Double_t * xmaxV0;
+  Int_t    * binsV0;
+  Int_t    fDim;
+
+
+
+  TObjArray* fRefTPC;
+  int fLabelsTPC[100000];
+
+
+  TClonesArray *clRefsN;
+  TClonesArray *clRefsP;
+
+  // MC variables
+
+  AliAnalysisTaskV0QA (const AliAnalysisTaskV0QA & v0QA );
+  AliAnalysisTaskV0QA & operator=(const AliAnalysisTaskV0QA & v0QA);
+
+
+
+  ClassDef(AliAnalysisTaskV0QA, 1); // example of analysis
+};
-- 
2.43.0