]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWG1/AliAnalysisTaskQASym.h
removing obsolete cout
[u/mrichter/AliRoot.git] / PWG1 / AliAnalysisTaskQASym.h
index 3958df30341f3f5b392306526da5b40044277dd9..c80a2d03ca2407bea1c9f7b0bcb3e80734882e2f 100644 (file)
@@ -4,6 +4,7 @@
 
 class TH1F;
 class TH2F;
+class TH3F;
 class TList;
 class TNtuple;
 
@@ -24,17 +25,18 @@ class AliAnalysisTaskQASym : public AliAnalysisTaskSE {
   virtual void   UserCreateOutputObjects();
   virtual void   UserExec(Option_t *option);
   virtual void   Terminate(Option_t *);
-
+  virtual void   SetTrackType(Int_t type) {fTrackType = type;}  
   
   
   virtual void   SetCuts(AliESDtrackCuts* cuts)
      {fCuts = cuts;}
 
   virtual void   SetFieldOn(Bool_t b = kTRUE){fFieldOn = b;} 
+
   
  private:
 
+  Int_t       fTrackType;
   Bool_t      fFieldOn;
 
   TList       *fHists;          // List of histos
@@ -46,102 +48,146 @@ class AliAnalysisTaskQASym : public AliAnalysisTaskSE {
   TH1F        *fEtaPt;          // eta over pt 
   TH1F        *fQPt;            // charge over pt 
   TH1F        *fDca;            // distance of closest approach
-  TH1F        *fqPtRec[7];      // charge over pt divided for ITS layer cases
   TH1F        *fqRec;           // reconstrcuted charge
   TH1F        *fsigmaPt;        // sigma_pT
-  TH2F        *fDcaSigmaPos[7]; // distance of closest approach against sigma_pT
-  TH2F        *fDcaSigmaNeg[7]; // distance of closest approach against sigma_pT
-
-
-  //positive und negative particles
-  TH1F        *fRecPtPos;      // pt of pos partice
-  TH1F        *fRecPtNeg;      // pt of neg particle
-  TH1F        *fRecPhiPos;     // phi of pos. particle
-  TH1F        *fRecPhiNeg;     // phi of neg. particle
-  TH1F        *fRecEtaPos;     // eta of neg. particle
-  TH1F        *fRecEtaNeg;     // eta of neg. particle
-  TH1F        *fRecEtaPtPos;   // eta over pt of neg. particle
-  TH1F        *fRecEtaPtNeg;   // eta over pt of neg. particle
-  TH1F        *fRecDcaPos;     // distance of closest approach of neg. particle
-  TH1F        *fRecDcaNeg;     // distance of closest of neg. particle
-  TH1F        *fRecDcaNegInv;  // invers dca of neg. particle
-  TH1F        *fRecDPos;       // impact parameter of neg. particle
-  TH1F        *fRecDNeg;       // impact parameter of neg. particle
-
-
 
+   //positive und negative tracks
+  TH1F        *fRecPtPos;      // pt of pos tracks
+  TH1F        *fRecPtNeg;      // pt of neg tracks
+  TH1F        *fRecPhiPos;     // phi of pos. tracks
+  TH1F        *fRecPhiNeg;     // phi of neg. tracks
+  TH1F        *fRecEtaPos;     // eta of neg. tracks
+  TH1F        *fRecEtaNeg;     // eta of neg. tracks
+  TH1F        *fRecEtaPtPos;   // eta over pt of neg. tracks
+  TH1F        *fRecEtaPtNeg;   // eta over pt of neg. tracks
+  TH1F        *fRecDcaPos;     // distance of closest approach of neg. tracks
+  TH1F        *fRecDcaNeg;     // distance of closest of neg. tracks
+  TH1F        *fRecDcaNegInv;  // invers dca of neg. tracks
+  TH1F        *fRecDPos;       // impact parameter of neg. tracks
+  TH1F        *fRecDNeg;       // impact parameter of neg. tracks
 
   // two sides of TPC -> Eta/Theta
-  TH1F        *fRecQPtPosEta; 
-  TH1F        *fRecQPtNegEta;
-  TH1F        *fRecPtPosEta;  
-  TH1F        *fRecPtNegEta; 
-  TH1F        *fRecPhiPosEta;
-  TH1F        *fRecPhiNegEta;
-  TH1F        *fRecDcaPosEta;
-  TH1F        *fRecDcaNegEta;
-  TH1F        *fRecDPosEta;
-  TH1F        *fRecDNegEta;
-
-
-  // sectors of TPC (with pt>xGeV?), TODO: extent to TPC standalone tracks
-  TH1F        *fRecPtTpcSector[18];
-  TH1F        *fRecEtaTpcSector[18];
-  TH1F        *fRecQPtTpcSector[18];
-  TH1F        *fRecEtaPtTpcSector[18];
-  TH1F        *fSignedDcaTpcSector[18];
-
-
-  // 7 different case of hit in ITS ladders
-  TH1F        *fSignDcaPos[7];
-  TH1F        *fSignDcaNeg[7];
-  TH1F        *fSignDcaNegInv[7];
-  TH1F        *fPtSigmaPos[7];
-  TH1F        *fPtSigmaNeg[7];
-  TH1F        *fRecPtPosLadder[7];
-  TH1F        *fRecPtNegLadder[7];
-  TH1F        *fRecPhiPosLadder[7];
-  TH1F        *fRecPhiNegLadder[7];
-  TH1F        *fRecEtaPosLadder[7];
-  TH1F        *fRecEtaNegLadder[7];
-
-
-  // 2D: all measures as function of z of vertex
-  TH2F        *fRecPtPosVz;
-  TH2F        *fRecPtNegVz;
-  TH2F        *fRecEtaPosVz;
-  TH2F        *fRecEtaNegVz;
-  TH2F        *fRecPhiPosVz;
-  TH2F        *fRecPhiNegVz;
-  TH2F        *fSignedDcaPosVz;
-  TH2F        *fSignedDcaNegVz;
-  TH2F        *fRecQPtPosEtaVz;
-  TH2F        *fRecQPtNegEtaVz;
-  TH2F        *fRecEtaPtPosVz;
-  TH2F        *fRecEtaPtNegVz;
+  TH1F        *fRecQPtPosEta;   //charge/pT for pos. eta
+  TH1F        *fRecQPtNegEta;   //charge/pT for neg. eta
+  TH1F        *fRecPtPosEta;    //pT        for pos. eta
+  TH1F        *fRecPtNegEta;    //pT        for neg. eta
+  TH1F        *fRecPhiPosEta;   //phi       for pos. eta
+  TH1F        *fRecPhiNegEta;   //phi       for neg. eta
+  TH1F        *fRecDcaPosEta;   //dca       for pos. eta
+  TH1F        *fRecDcaNegEta;   //dca       for neg. eta
+  TH1F        *fRecDPosEta;     //d         for pos. eta
+  TH1F        *fRecDNegEta;     //d         for neg. eta
+
+  // 2D: all measures as function of z of first trackpoint 
+  TH2F        *fRecPtPosVz;     //pt-zfirst of pos tracks
+  TH2F        *fRecPtNegVz;      //pt-zfirst of neg tracks
+  TH2F        *fRecEtaPosVz;    //eta-zfirst of pos tracks
+  TH2F        *fRecEtaNegVz;    //eta-zfirst of neg tracks
+  TH2F        *fRecPhiPosVz;    //phi-zfirst of pos tracks
+  TH2F        *fRecPhiNegVz;    //phi-zfirst of neg tracks
+  TH2F        *fSignedDcaPosVz; //dca-zfirst of pos tracks
+  TH2F        *fSignedDcaNegVz; //dca-zfirst of neg tracks
+  TH2F        *fRecQPtPosEtaVz; //charge/pT-zfirst of pos tracks
+  TH2F        *fRecQPtNegEtaVz; //charge/pT-zfirst of neg tracks
+  TH2F        *fRecEtaPtPosVz;  //eta/pT-zfirst of pos tracks
+  TH2F        *fRecEtaPtNegVz;  //eta/pT-zfirst of neg tracks
 
 
   //high
-  TH1F * fDeltaPhiAll;
-  TH2F * fDeltaPhiLeading;
-  TH1F * fDiffDcaD;
+  TH1F * fDeltaPhiAll;         // phiLeaingTracks-phiOthers
+  TH2F * fDeltaPhiLeading;     // phiLeaingTracks-phiOthers vs. phiLeading
+  TH1F * fDiffDcaD;            // d-dca
 
   //sim
-  TH1F * fPhiRec;
-  TH1F * fThetaRec;
-  TH1F * fNumber;
-  TH1F * fVx;
-  TH1F * fVy;
-  TH1F * fVz;
+  TH1F * fPhiRec;              //phi
+  TH1F * fThetaRec;            //theta
+  TH1F * fNumber;              //Number of tracks per event
+  TH1F * fVx;                  // x of first track point
+  TH1F * fVy;                  // y of first track point
+  TH1F * fVz;                  // z of first track point
   TNtuple * test;
 
+  //new
+  TH2F        *fRecDcaPosPhi;     //dca-phi for pos.
+  TH2F        *fRecDcaNegPhi;     //dca-phi for neg.
+  TH2F        *fRecPtPosPhi;      //pt-phi for pos.
+  TH2F        *fRecPtNegPhi;      //pt-phi for neg.
+  TH2F        *fRecEtaPosPhi;     //eta-phi for pos.
+  TH2F        *fRecEtaNegPhi;     //eta-phi for neg.
+  TH2F        *fRecQPtPhi;        //charge/pt-phi
+  TH2F        *fRecEtaPtPosPhi;   //eta/pt-phi for neg.
+  TH2F        *fRecEtaPtNegPhi;   //eta/pt-phi for pos.
+
+  TH1F        *fRecPtPosEtaPos;   //pt for pos tracks and pos eta
+  TH1F        *fRecPtNegEtaPos;   //pt for neg tracks and pos eta
+  TH1F        *fRecPtPosEtaNeg;   //pt for pos tracks and neg eta
+  TH1F        *fRecPtNegEtaNeg;   //pt for neg tracks and neg eta
+
+  TH1F        *fRec1PtPosEtaPos;   //1/pt for pos tracks and pos eta
+  TH1F        *fRec1PtNegEtaPos;   //1/pt for neg tracks and pos eta
+  TH1F        *fRec1PtPosEtaNeg;   //1/pt for pos tracks and neg eta
+  TH1F        *fRec1PtNegEtaNeg;   //1/pt for neg tracks and neg eta
+
+  TH1F        *fRecPhiPosEtaPos;   //phi for pos tracks and pos eta
+  TH1F        *fRecPhiNegEtaPos;   //phi for neg tracks and pos eta
+  TH1F        *fRecPhiPosEtaNeg;   //phi for pos tracks and neg eta
+  TH1F        *fRecPhiNegEtaNeg;   //phi for neg tracks and neg eta
+
+  TH2F        *fRecDcaPosPhiEtaPos;  //dca-phi for pos tracks and pos eta
+  TH2F        *fRecDcaNegPhiEtaPos;  //dca-phi for neg tracks and pos eta
+  TH2F        *fRecDcaPosPhiEtaNeg;  //dca-phi for pos tracks and neg eta
+  TH2F        *fRecDcaNegPhiEtaNeg;  //dca-phi for neg tracks and neg eta
+
+  TH2F        *fRecDcaPosPtEtaPos;  //dca-pt for pos tracks and pos eta
+  TH2F        *fRecDcaNegPtEtaPos;  //dca-pt for neg tracks and pos eta
+  TH2F        *fRecDcaPosPtEtaNeg;  //dca-pt for pos tracks and neg eta
+  TH2F        *fRecDcaNegPtEtaNeg;  //dca-pt for neg tracks and neg eta
+
+  TH2F        *fRecPtPosPhiEtaPos;  //pt-phi for pos tracks and pos eta
+  TH2F        *fRecPtNegPhiEtaPos;  //pt-phi for neg tracks and pos eta 
+  TH2F        *fRecPtPosPhiEtaNeg;  //pt-phi for pos tracks and neg eta
+  TH2F        *fRecPtNegPhiEtaNeg;  //pt-phi for neg tracks and neg eta
+
+  TH3F        *fRecDcaPhiPtPosEtaPos; //dca-pt-phi for pos tracks and pos eta
+  TH3F        *fRecDcaPhiPtNegEtaPos; //dca-pt-phi for neg tracks and pos eta
+  TH3F        *fRecDcaPhiPtPosEtaNeg; //dca-pt-phi for pos tracks and neg eta
+  TH3F        *fRecDcaPhiPtNegEtaNeg; //dca-pt-phi for neg tracks and neg eta
+
+  TH2F        *fEtavPt;
+  TH2F        *fCompareTPCparam;
+
+  Double_t  sdca;
+  Float_t xy, z, xvertexcor, yvertexcor;
+  AliESDtrackCuts* fCuts;                      // List of cuts
 
+  // sectors of TPC 
+  TH1F        *fRecPtTpcSector[18];     //pt for TPC sectors
+  TH1F        *fRecEtaTpcSector[18];    //eta for TPC sectors
+  TH1F        *fRecQPtTpcSector[18];    //charge/pt for TPC sectors
+  TH1F        *fRecEtaPtTpcSector[18];  //eta/pt for TPC sectors
+  TH1F        *fSignedDcaTpcSector[18]; //dca for TPC sectors
+
+
+  // 7 different case of hit in ITS ladders
+  TH1F        *fRecPtPosLadder[7];  //pt for pos tracks
+  TH1F        *fRecPtNegLadder[7];  //pt for neg tracks
+  TH1F        *fRecPhiPosLadder[7]; //phi for pos tracks
+  TH1F        *fRecPhiNegLadder[7]; //phi for neg tracks
+  TH1F        *fRecEtaPosLadder[7]; //eta for pos tracks
+  TH1F        *fRecEtaNegLadder[7]; //eta for neg tracks
+  TH1F        *fSignDcaPos[7];      //dca for pos tracks
+  TH1F        *fSignDcaNeg[7];      //dca for neg tracks
+  TH1F        *fSignDcaNegInv[7];   //-dca for neg tracks
+  TH1F        *fPtSigmaPos[7];      //sigma_pT for pos tracks
+  TH1F        *fPtSigmaNeg[7];      //sigma_pT for neg tracks
+  TH1F        *fqPtRec[7];          // charge/pt 
+  TH2F        *fDcaSigmaPos[7];     // dca - sigma_pT for pos tracks
+  TH2F        *fDcaSigmaNeg[7];     // dca - sigma_pT for neg tracks
+
+  
+  
 
-  Double_t  sdca_tr;
-  Float_t xy, z;
 
-  AliESDtrackCuts* fCuts;                      // List of cuts
   AliAnalysisTaskQASym(const AliAnalysisTaskQASym&); // not implemented
   AliAnalysisTaskQASym& operator=(const AliAnalysisTaskQASym&); // not implemented