PWG4/JetTasks/AliPWG4HighPtQATPConly.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16 //-----------------------------------------------------------------------
  17 // This class compares the global reconstruction with the TPConly
  18 // reconstruction
  19 // Momentum resolution is stored as function of track cuts and pt.
  20 // Output: Histograms for different set of cuts
  21 //-----------------------------------------------------------------------
  22 // Author : Marta Verweij - UU
  23 //-----------------------------------------------------------------------
  24
  25 #ifndef ALIPWG4HighPtQATPCONLY_CXX
  26 #define ALIPWG4HighPtQATPCONLY_CXX
  27
  28 #include "AliPWG4HighPtQATPConly.h"
  29
  30 #include "TVector3.h"
  31 #include <iostream>
  32 #include "TH1.h"
  33 #include "TH2.h"
  34 #include "TH3.h"
  35 #include "TList.h"
  36 #include "TChain.h"
  37 #include "TH3F.h"
  38
  39 #include "AliAnalysisManager.h"
  40 #include "AliESDInputHandler.h"
  41 #include "AliESDtrack.h"
  42 #include "AliESDtrackCuts.h"
  43 #include "AliExternalTrackParam.h"
  44 #include "AliLog.h"
  45 #include "AliAnalysisHelperJetTasks.h"
  46
  47 using namespace std; //required for resolving the 'cout' symbol
  48
  49 ClassImp(AliPWG4HighPtQATPConly)
  50
  51 AliPWG4HighPtQATPConly::AliPWG4HighPtQATPConly(): AliAnalysisTask("AliPWG4HighPtQATPConly", ""),
  52   fESD(0),
  53   fTrackCuts(0),
  54   fTrackCutsITS(0),
  55   fTrigger(0),
  56   fNEventAll(0),
  57   fNEventSel(0),
  58   fPtAll(0),
  59   fPtSel(0),
  60   fPtAllminPtTPCvsPtAll(0),
  61   fPtAllminPtTPCvsPtAllNPointTPC(0),
  62   fPtAllminPtTPCvsPtAllDCAR(0),
  63   fPtAllminPtTPCvsPtAllDCAZ(0),
  64   fPtAllminPtTPCvsPtAllPhi(0),
  65   fPtAllminPtTPCvsPtAllNPointITS(0),
  66   fPtAllminPtTPCvsPtAllNSigmaToVertex(0),
  67   fPtAllminPtTPCvsPtAllChi2C(0),
  68   fPtAllminPtTPCvsPtAllRel1PtUncertainty(0),
  69   fHistList(0),
  70   fPtAllTPC(0),
  71   fPtSelTPC(0),
  72   fPtSelTPCITS(0),
  73   fHistListTPC(0),
  74   fPtSelITS(0),
  75   fPtITSminPtTPCvsPtITS(0),
  76   fPtITSminPtTPCvsPtITSNPointTPC(0),
  77   fPtITSminPtTPCvsPtITSDCAR(0),
  78   fPtITSminPtTPCvsPtITSDCAZ(0),
  79   fPtITSminPtTPCvsPtITSPhi(0),
  80   fPtITSminPtTPCvsPtITSNPointITS(0),
  81   fPtITSminPtTPCvsPtITSNSigmaToVertex(0),
  82   fPtITSminPtTPCvsPtITSChi2C(0),
  83   fPtITSminPtTPCvsPtITSRel1PtUncertainty(0),
  84   fHistListITS(0)
  85 {
  86
  87 }
  88 //________________________________________________________________________
  89 AliPWG4HighPtQATPConly::AliPWG4HighPtQATPConly(const char *name):
  90   AliAnalysisTask(name, ""),
  91   fESD(0),
  92   fTrackCuts(),
  93   fTrackCutsITS(),
  94   fTrigger(0),
  95   fNEventAll(0),
  96   fNEventSel(0),
  97   fPtAll(0),
  98   fPtSel(0),
  99   fPtAllminPtTPCvsPtAll(0),
 100   fPtAllminPtTPCvsPtAllNPointTPC(0),
 101   fPtAllminPtTPCvsPtAllDCAR(0),
 102   fPtAllminPtTPCvsPtAllDCAZ(0),
 103   fPtAllminPtTPCvsPtAllPhi(0),
 104   fPtAllminPtTPCvsPtAllNPointITS(0),
 105   fPtAllminPtTPCvsPtAllNSigmaToVertex(0),
 106   fPtAllminPtTPCvsPtAllChi2C(0),
 107   fPtAllminPtTPCvsPtAllRel1PtUncertainty(0),
 108   fHistList(0),
 109   fPtAllTPC(0),
 110   fPtSelTPC(0),
 111   fPtSelTPCITS(0),
 112   fHistListTPC(0),
 113   fPtSelITS(0),
 114   fPtITSminPtTPCvsPtITS(0),
 115   fPtITSminPtTPCvsPtITSNPointTPC(0),
 116   fPtITSminPtTPCvsPtITSDCAR(0),
 117   fPtITSminPtTPCvsPtITSDCAZ(0),
 118   fPtITSminPtTPCvsPtITSPhi(0),
 119   fPtITSminPtTPCvsPtITSNPointITS(0),
 120   fPtITSminPtTPCvsPtITSNSigmaToVertex(0),
 121   fPtITSminPtTPCvsPtITSChi2C(0),
 122   fPtITSminPtTPCvsPtITSRel1PtUncertainty(0),
 123   fHistListITS(0)
 124 {
 125   //
 126   // Constructor. Initialization of Inputs and Outputs
 127   //
 128   Info("AliPWG4HighPtQATPConly","Calling Constructor");
 129   // Input slot #0 works with a TChain ESD
 130   DefineInput(0, TChain::Class());
 131   // Output slot #0 writes into a TList
 132   DefineOutput(0, TList::Class());
 133   // Output slot #1 writes into a TList
 134   DefineOutput(1, TList::Class());
 135   // Output slot #2 writes into a TList
 136   DefineOutput(2, TList::Class());
 137 }
 138
 139 //________________________________________________________________________
 140 void AliPWG4HighPtQATPConly::ConnectInputData(Option_t *)
 141 {
 142   // Connect ESD here
 143   // Called once
 144   AliDebug(2,Form(">> AliPWG4HighPtSpectra::ConnectInputData \n"));
 145   TTree* tree = dynamic_cast<TTree*> (GetInputData(0));
 146   if (!tree) {
 147     Printf("ERROR: Could not read chain from input slot 0");
 148   } else {
 149
 150     AliESDInputHandler *esdH = dynamic_cast<AliESDInputHandler*> (AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());
 151
 152     if (!esdH) {
 153       AliDebug(2,Form("ERROR: Could not get ESDInputHandler"));
 154     } else {
 155       fESD = esdH->GetEvent();
 156     }
 157   }
 158 }
 159
 160 //________________________________________________________________________
 161 void AliPWG4HighPtQATPConly::CreateOutputObjects() {
 162   //Create output objects
 163   AliDebug(2,Form(">> AliPWG4HighPtQATPConly::CreateOutputObjects \n"));
 164
 165   Bool_t oldStatus = TH1::AddDirectoryStatus();
 166   TH1::AddDirectory(kFALSE);
 167
 168   OpenFile(0);
 169   fHistList = new TList();
 170   OpenFile(1);
 171   fHistListTPC = new TList();
 172   OpenFile(2);
 173   fHistListITS = new TList();
 174
 175   Int_t fgkNPhiBins=18;
 176   Float_t kMinPhi = 0.;
 177   Float_t kMaxPhi = 2.*TMath::Pi();
 178
 179   Int_t fgkNPtBins=98;
 180   Float_t fgkPtMin=2.;
 181   Float_t fgkPtMax=100.;
 182   Int_t fgkResPtBins=80;
 183
 184   fNEventAll = new TH1F("fNEventAll","NEventAll",1,-0.5,0.5);
 185   fHistList->Add(fNEventAll);
 186   fNEventSel = new TH1F("fNEventSel","NEvent Selected for analysis",1,-0.5,0.5);
 187   fHistList->Add(fNEventSel);
 188   fPtAll = new TH1F("fPtAll","PtAll",fgkNPtBins, fgkPtMin, fgkPtMax);
 189   fHistList->Add(fPtAll);
 190   fPtSel = new TH1F("fPtSel","PtSel",fgkNPtBins, fgkPtMin, fgkPtMax);
 191   fHistList->Add(fPtSel);
 192
 193   fPtAllminPtTPCvsPtAll = new TH2F("fPtAllminPtTPCvsPtAll","PtAllminPtTPCvsPtAll",fgkNPtBins, fgkPtMin,fgkPtMax,fgkResPtBins,-1.,1.);
 194   fPtAllminPtTPCvsPtAll->SetXTitle("p_{t}^{All}");
 195   fPtAllminPtTPCvsPtAll->SetYTitle("(1/p_{t}^{All}-1/p_{t}^{TPC})/(1/p_{t}^{All})");
 196   fHistList->Add(fPtAllminPtTPCvsPtAll);
 197
 198   fPtAllminPtTPCvsPtAllNPointTPC = new TH3F("fPtAllminPtTPCvsPtAllNPointTPC","PtAllminPtTPCvsPtAllNPointTPC",fgkNPtBins, fgkPtMin,fgkPtMax,fgkResPtBins,-1.,1.,160,0.5,160.5);
 199   fPtAllminPtTPCvsPtAllNPointTPC->SetXTitle("p_{t}^{All}");
 200   fPtAllminPtTPCvsPtAllNPointTPC->SetYTitle("(1/p_{t}^{All}-1/p_{t}^{TPC})/(1/p_{t}^{All})");
 201   fPtAllminPtTPCvsPtAllNPointTPC->SetZTitle("N_{point,TPC}");
 202   fHistList->Add(fPtAllminPtTPCvsPtAllNPointTPC);
 203
 204   fPtAllminPtTPCvsPtAllDCAR = new TH3F("fPtAllminPtTPCvsPtAllDCAR","PtAllminPtTPCvsPtAllDCAR",fgkNPtBins, fgkPtMin,fgkPtMax,fgkResPtBins,-1.,1.,20,-1.,1.);
 205   fPtAllminPtTPCvsPtAllDCAR->SetXTitle("p_{t}^{All}");
 206   fPtAllminPtTPCvsPtAllDCAR->SetYTitle("(1/p_{t}^{All}-1/p_{t}^{TPC})/(1/p_{t}^{All})");
 207   fPtAllminPtTPCvsPtAllDCAR->SetZTitle("DCA_{R}");
 208   fHistList->Add(fPtAllminPtTPCvsPtAllDCAR);
 209
 210   fPtAllminPtTPCvsPtAllDCAZ = new TH3F("fPtAllminPtTPCvsPtAllDCAZ","PtAllminPtTPCvsPtAllDCAZ",fgkNPtBins, fgkPtMin,fgkPtMax,fgkResPtBins,-1.,1.,40,-2.,2.);
 211   fPtAllminPtTPCvsPtAllDCAZ->SetXTitle("p_{t}^{All}");
 212   fPtAllminPtTPCvsPtAllDCAZ->SetYTitle("(1/p_{t}^{All}-1/p_{t}^{TPC})/(1/p_{t}^{All})");
 213   fPtAllminPtTPCvsPtAllDCAZ->SetZTitle("DCA_{Z}");
 214   fHistList->Add(fPtAllminPtTPCvsPtAllDCAZ);
 215
 216   fPtAllminPtTPCvsPtAllPhi = new TH3F("fPtAllminPtTPCvsPtAllPhi","PtAllminPtTPCvsPtAllPhi",fgkNPtBins, fgkPtMin,fgkPtMax,fgkResPtBins,-1.,1.,fgkNPhiBins,kMinPhi,kMaxPhi);
 217   fPtAllminPtTPCvsPtAllPhi->SetXTitle("p_{t}^{All}");
 218   fPtAllminPtTPCvsPtAllPhi->SetYTitle("(1/p_{t}^{All}-1/p_{t}^{TPC})/(1/p_{t}^{All})");
 219   fPtAllminPtTPCvsPtAllPhi->SetZTitle("#phi");
 220   fHistList->Add(fPtAllminPtTPCvsPtAllPhi);
 221
 222   fPtAllminPtTPCvsPtAllNPointITS = new TH3F("fPtAllminPtTPCvsPtAllNPointITS","PtAllminPtTPCvsPtAllNPointITS",fgkNPtBins, fgkPtMin,fgkPtMax,fgkResPtBins,-1.,1.,9,-0.5,8.5);
 223   fPtAllminPtTPCvsPtAllNPointITS->SetXTitle("p_{t}^{All}");
 224   fPtAllminPtTPCvsPtAllNPointITS->SetYTitle("(1/p_{t}^{All}-1/p_{t}^{TPC})/(1/p_{t}^{All})");
 225   fPtAllminPtTPCvsPtAllNPointITS->SetZTitle("N_{point,ITS}}");
 226   fHistList->Add(fPtAllminPtTPCvsPtAllNPointITS);
 227
 228   fPtAllminPtTPCvsPtAllNSigmaToVertex = new TH3F("fPtAllminPtTPCvsPtAllNSigmaToVertex","PtAllminPtTPCvsPtAllNSigmaToVertex",fgkNPtBins, fgkPtMin,fgkPtMax,fgkResPtBins,-1.,1.,40,0.,8.);
 229   fPtAllminPtTPCvsPtAllNSigmaToVertex->SetXTitle("p_{t}^{All}");
 230   fPtAllminPtTPCvsPtAllNSigmaToVertex->SetYTitle("(1/p_{t}^{All}-1/p_{t}^{TPC})/(1/p_{t}^{All})");
 231   fPtAllminPtTPCvsPtAllNSigmaToVertex->SetZTitle("N#sigma to vertex");
 232   fHistList->Add(fPtAllminPtTPCvsPtAllNSigmaToVertex);
 233
 234   fPtAllminPtTPCvsPtAllChi2C = new TH3F("fPtAllminPtTPCvsPtAllChi2C","PtAllminPtTPCvsPtAllChi2C",fgkNPtBins, fgkPtMin,fgkPtMax,fgkResPtBins,-1.,1.,20,0.,10.);
 235   fPtAllminPtTPCvsPtAllChi2C->SetXTitle("p_{t}^{All}");
 236   fPtAllminPtTPCvsPtAllChi2C->SetYTitle("(1/p_{t}^{All}-1/p_{t}^{TPC})/(1/p_{t}^{All})");
 237   fPtAllminPtTPCvsPtAllChi2C->SetZTitle("Constrained #chi^{2}");
 238   fHistList->Add(fPtAllminPtTPCvsPtAllChi2C);
 239
 240   fPtAllminPtTPCvsPtAllRel1PtUncertainty = new TH3F("fPtAllminPtTPCvsPtAllRel1PtUncertainty","PtAllminPtTPCvsPtAllRel1PtUncertainty",fgkNPtBins, fgkPtMin,fgkPtMax,fgkResPtBins,-1.,1.,30,0.,0.3);
 241   fPtAllminPtTPCvsPtAllRel1PtUncertainty->SetXTitle("p_{t}^{All}");
 242   fPtAllminPtTPCvsPtAllRel1PtUncertainty->SetYTitle("(1/p_{t}^{All}-1/p_{t}^{TPC})/(1/p_{t}^{All})");
 243   fPtAllminPtTPCvsPtAllRel1PtUncertainty->SetZTitle("Rel1PtUncertainty");
 244   fHistList->Add(fPtAllminPtTPCvsPtAllRel1PtUncertainty);
 245
 246   //ITSrefit
 247   fPtSelITS = new TH1F("fPtSelITSrefit","PtSel",fgkNPtBins, fgkPtMin, fgkPtMax);
 248   fHistListITS->Add(fPtSelITS);
 249
 250   fPtITSminPtTPCvsPtITS = new TH2F("fPtITSminPtTPCvsPtITS","PtITSminPtTPCvsPtITS",fgkNPtBins, fgkPtMin,fgkPtMax,fgkResPtBins,-1.,1.);
 251   fPtITSminPtTPCvsPtITS->SetXTitle("p_{t}^{ITS}");
 252   fPtITSminPtTPCvsPtITS->SetYTitle("(1/p_{t}^{ITS}-1/p_{t}^{TPC})/(1/p_{t}^{ITS})");
 253   fHistListITS->Add(fPtITSminPtTPCvsPtITS);
 254
 255   fPtITSminPtTPCvsPtITSNPointTPC = new TH3F("fPtITSminPtTPCvsPtITSNPointTPC","PtITSminPtTPCvsPtITSNPointTPC",fgkNPtBins, fgkPtMin,fgkPtMax,fgkResPtBins,-1.,1.,160,0.5,160.5);
 256   fPtITSminPtTPCvsPtITSNPointTPC->SetXTitle("p_{t}^{ITSrefit}");
 257   fPtITSminPtTPCvsPtITSNPointTPC->SetYTitle("(1/p_{t}^{ITSrefit}-1/p_{t}^{TPC})/(1/p_{t}^{ITSrefit})");
 258   fPtITSminPtTPCvsPtITSNPointTPC->SetZTitle("N_{point,TPC}");
 259   fHistListITS->Add(fPtITSminPtTPCvsPtITSNPointTPC);
 260
 261   fPtITSminPtTPCvsPtITSDCAR = new TH3F("fPtITSminPtTPCvsPtITSDCAR","PtITSminPtTPCvsPtITSDCAR",fgkNPtBins, fgkPtMin,fgkPtMax,fgkResPtBins,-1.,1.,20,-1.,1.);
 262   fPtITSminPtTPCvsPtITSDCAR->SetXTitle("p_{t}^{ITSrefit}");
 263   fPtITSminPtTPCvsPtITSDCAR->SetYTitle("(1/p_{t}^{ITSrefit}-1/p_{t}^{TPC})/(1/p_{t}^{ITSrefit})");
 264   fPtITSminPtTPCvsPtITSDCAR->SetZTitle("DCA_{R}");
 265   fHistListITS->Add(fPtITSminPtTPCvsPtITSDCAR);
 266
 267   fPtITSminPtTPCvsPtITSDCAZ = new TH3F("fPtITSminPtTPCvsPtITSDCAZ","PtITSminPtTPCvsPtITSDCAZ",fgkNPtBins, fgkPtMin,fgkPtMax,fgkResPtBins,-1.,1.,40,-2.,2.);
 268   fPtITSminPtTPCvsPtITSDCAZ->SetXTitle("p_{t}^{ITSrefit}");
 269   fPtITSminPtTPCvsPtITSDCAZ->SetYTitle("(1/p_{t}^{ITSrefit}-1/p_{t}^{TPC})/(1/p_{t}^{ITSrefit})");
 270   fPtITSminPtTPCvsPtITSDCAZ->SetZTitle("DCA_{Z}");
 271   fHistListITS->Add(fPtITSminPtTPCvsPtITSDCAZ);
 272
 273   fPtITSminPtTPCvsPtITSPhi = new TH3F("fPtITSminPtTPCvsPtITSPhi","PtITSminPtTPCvsPtITSPhi",fgkNPtBins, fgkPtMin,fgkPtMax,fgkResPtBins,-1.,1.,fgkNPhiBins,kMinPhi,kMaxPhi);
 274   fPtITSminPtTPCvsPtITSPhi->SetXTitle("p_{t}^{ITSrefit}");
 275   fPtITSminPtTPCvsPtITSPhi->SetYTitle("(1/p_{t}^{ITSrefit}-1/p_{t}^{TPC})/(1/p_{t}^{ITSrefit})");
 276   fPtITSminPtTPCvsPtITSPhi->SetZTitle("#phi");
 277   fHistListITS->Add(fPtITSminPtTPCvsPtITSPhi);
 278
 279   fPtITSminPtTPCvsPtITSNPointITS = new TH3F("fPtITSminPtTPCvsPtITSNPointITS","PtITSminPtTPCvsPtITSNPointITS",fgkNPtBins, fgkPtMin,fgkPtMax,fgkResPtBins,-1.,1.,9,-0.5,8.5);
 280   fPtITSminPtTPCvsPtITSNPointITS->SetXTitle("p_{t}^{ITSrefit}");
 281   fPtITSminPtTPCvsPtITSNPointITS->SetYTitle("(1/p_{t}^{ITSrefit}-1/p_{t}^{TPC})/(1/p_{t}^{ITSrefit})");
 282   fPtITSminPtTPCvsPtITSNPointITS->SetZTitle("N_{point,ITS}}");
 283   fHistListITS->Add(fPtITSminPtTPCvsPtITSNPointITS);
 284
 285   fPtITSminPtTPCvsPtITSNSigmaToVertex = new TH3F("fPtITSminPtTPCvsPtITSNSigmaToVertex","PtITSminPtTPCvsPtITSNSigmaToVertex",fgkNPtBins, fgkPtMin,fgkPtMax,fgkResPtBins,-1.,1.,40,0.,8.);
 286   fPtITSminPtTPCvsPtITSNSigmaToVertex->SetXTitle("p_{t}^{ITS}");
 287   fPtITSminPtTPCvsPtITSNSigmaToVertex->SetYTitle("(1/p_{t}^{ITS}-1/p_{t}^{TPC})/(1/p_{t}^{ITS})");
 288   fPtITSminPtTPCvsPtITSNSigmaToVertex->SetZTitle("N#sigma to vertex");
 289   fHistListITS->Add(fPtITSminPtTPCvsPtITSNSigmaToVertex);
 290
 291   fPtITSminPtTPCvsPtITSChi2C = new TH3F("fPtITSminPtTPCvsPtITSChi2C","PtITSminPtTPCvsPtITSChi2C",fgkNPtBins, fgkPtMin,fgkPtMax,fgkResPtBins,-1.,1.,20,0.,10.);
 292   fPtITSminPtTPCvsPtITSChi2C->SetXTitle("p_{t}^{ITS}");
 293   fPtITSminPtTPCvsPtITSChi2C->SetYTitle("(1/p_{t}^{ITS}-1/p_{t}^{TPC})/(1/p_{t}^{ITS})");
 294   fPtITSminPtTPCvsPtITSChi2C->SetZTitle("Constrained #chi^{2}");
 295   fHistListITS->Add(fPtITSminPtTPCvsPtITSChi2C);
 296
 297   fPtITSminPtTPCvsPtITSRel1PtUncertainty = new TH3F("fPtITSminPtTPCvsPtITSRel1PtUncertainty","PtITSminPtTPCvsPtITSRel1PtUncertainty",fgkNPtBins, fgkPtMin,fgkPtMax,fgkResPtBins,-1.,1.,30,0.,0.3);
 298   fPtITSminPtTPCvsPtITSRel1PtUncertainty->SetXTitle("p_{t}^{ITS}");
 299   fPtITSminPtTPCvsPtITSRel1PtUncertainty->SetYTitle("(1/p_{t}^{ITS}-1/p_{t}^{TPC})/(1/p_{t}^{ITS})");
 300   fPtITSminPtTPCvsPtITSRel1PtUncertainty->SetZTitle("Rel1PtUncertainty");
 301   fHistListITS->Add(fPtITSminPtTPCvsPtITSRel1PtUncertainty);
 302
 303   fPtAllTPC = new TH1F("fPtAllTPC","PtAll",fgkNPtBins, fgkPtMin, fgkPtMax);
 304   fHistListTPC->Add(fPtAllTPC);
 305   fPtSelTPC = new TH1F("fPtSelTPC","PtSel",fgkNPtBins, fgkPtMin, fgkPtMax);
 306   fHistListTPC->Add(fPtSelTPC);
 307   fPtSelTPCITS = new TH1F("fPtSelTPCITS","PtSel",fgkNPtBins, fgkPtMin, fgkPtMax);
 308   fHistListTPC->Add(fPtSelTPCITS);
 309
 310   TH1::AddDirectory(oldStatus);
 311
 312 }
 313 //________________________________________________________________________
 314 void AliPWG4HighPtQATPConly::Exec(Option_t *) {
 315   // Main loop
 316   // Called for each event
 317   AliDebug(2,Form(">> AliPWG4HighPtQATPConly::Exec \n"));
 318
 319   // All events without selection
 320   cout << "Fill fNEventAll" << endl;
 321   fNEventAll->Fill(0.);
 322
 323   if (!fESD) {
 324     AliDebug(2,Form("ERROR: fESD not available"));
 325     return;
 326   }
 327
 328   //Trigger selection
 329   AliAnalysisHelperJetTasks::Trigger trig;
 330   trig = (const enum AliAnalysisHelperJetTasks::Trigger)fTrigger;
 331   if (AliAnalysisHelperJetTasks::IsTriggerFired(fESD,trig)){
 332     AliDebug(2,Form(" Trigger Selection: event ACCEPTED ... "));
 333   }else{
 334     AliDebug(2,Form(" Trigger Selection: event REJECTED ... "));
 335     return;
 336   }
 337
 338 //  if(!fESD->IsTriggerClassFired("CINT1B-ABCE-NOPF-ALL") || !fESD->IsTriggerClassFired("CSMBB-ABCE-NOPF-ALL")) return;
 339
 340   const AliESDVertex *vtx = fESD->GetPrimaryVertexTracks();
 341   // Need vertex cut
 342   if (vtx->GetNContributors() < 2)
 343     return;
 344
 345   AliDebug(2,Form("Vertex title %s, status %d, nCont %d\n",vtx->GetTitle(), vtx->GetStatus(), vtx->GetNContributors()));
 346   double primVtx[3];
 347   vtx->GetXYZ(primVtx);
 348   //  printf("primVtx: %g  %g  %g \n",primVtx[0],primVtx[1],primVtx[2]);
 349   if(primVtx[0]>1. || primVtx[1]>1. || primVtx[2]>10.) return;
 350
 351   if(!fESD->GetNumberOfTracks() || fESD->GetNumberOfTracks()<2) return;
 352   Int_t nTracks = fESD->GetNumberOfTracks();
 353   AliDebug(2,Form("nTracks %d\n", nTracks));
 354
 355   if(!fTrackCuts) return;
 356
 357   // Selected events for analysis
 358   fNEventSel->Fill(0.);
 359
 360   for (Int_t iTrack = 0; iTrack < nTracks; iTrack++) {
 361
 362     AliESDtrack *track = fESD->GetTrack(iTrack);
 363     AliExternalTrackParam *trackTPC = (AliExternalTrackParam *)track->GetTPCInnerParam();
 364     //    AliESDtrack *trackTPConly = fTrackCuts->GetTPCOnlyTrack(fESD,iTrack);
 365     if(!track) continue;
 366     if(!trackTPC) continue;
 367 //     if(!trackTPConly) continue;
 368     Float_t pt = track->Pt();
 369     Float_t ptTPC = trackTPC->Pt();
 370     Float_t phi = track->Phi();
 371     Float_t dca2D, dcaZ;
 372     track->GetImpactParameters(dca2D,dcaZ);
 373     // Float_t dca2DTPC, dcaZTPC;
 374     //track->GetImpactParametersTPC(dca2DTPC,dcaZTPC);
 375     UChar_t itsMap = track->GetITSClusterMap();
 376     Int_t nPointITS = 0;
 377     for (Int_t i=0; i < 6; i++) {
 378       if (itsMap & (1 << i))
 379         nPointITS ++;
 380     }
 381     double mom[3];
 382     track->GetPxPyPz(mom);
 383     Float_t nSigmaToVertex = fTrackCuts->GetSigmaToVertex(track);// Calculates the number of sigma to the vertex for a track.
 384     Float_t chi2C = track->GetConstrainedChi2();
 385     Float_t relUncertainty1Pt = TMath::Sqrt(TMath::Abs(track->GetSigma1Pt2()))*pt;
 386
 387     fPtAll->Fill(pt);
 388     fPtAllTPC->Fill(ptTPC);
 389
 390     if (fTrackCuts->AcceptTrack(track)) {
 391
 392       fPtSel->Fill(pt);
 393
 394       fPtSelTPC->Fill(ptTPC);
 395       fPtAllminPtTPCvsPtAll->Fill(pt,(1./pt-1./ptTPC)/(1./pt) );
 396       fPtAllminPtTPCvsPtAllNPointTPC->Fill(pt,(1./pt-1./ptTPC)/(1./pt),track->GetTPCNcls());
 397       fPtAllminPtTPCvsPtAllDCAR->Fill(pt,(1./pt-1./ptTPC)/(1./pt),dca2D);
 398       fPtAllminPtTPCvsPtAllDCAZ->Fill(pt,(1./pt-1./ptTPC)/(1./pt),dcaZ);
 399       fPtAllminPtTPCvsPtAllPhi->Fill(pt,(1./pt-1./ptTPC)/(1./pt),phi);
 400       fPtAllminPtTPCvsPtAllNPointITS->Fill(pt,(1./pt-1./ptTPC)/(1./pt),nPointITS);
 401       fPtAllminPtTPCvsPtAllNSigmaToVertex->Fill(pt,(1./pt-1./ptTPC)/(1./pt),nSigmaToVertex);
 402       fPtAllminPtTPCvsPtAllChi2C->Fill(pt,(1./pt-1./ptTPC)/(1./pt),chi2C);
 403       fPtAllminPtTPCvsPtAllRel1PtUncertainty->Fill(pt,(1./pt-1./ptTPC)/(1./pt),relUncertainty1Pt);
 404     }//fTrackCuts selection
 405
 406
 407     //ITSrefit selection
 408     if (fTrackCutsITS->AcceptTrack(track)) {
 409
 410       fPtSelITS->Fill(pt);
 411       fPtSelTPCITS->Fill(ptTPC);
 412       fPtITSminPtTPCvsPtITS->Fill(pt,(1./pt-1./ptTPC)/(1./pt) );
 413       fPtITSminPtTPCvsPtITSNPointTPC->Fill(pt,(1./pt-1./ptTPC)/(1./pt),track->GetTPCNcls());
 414       fPtITSminPtTPCvsPtITSDCAR->Fill(pt,(1./pt-1./ptTPC)/(1./pt),dca2D);
 415       fPtITSminPtTPCvsPtITSDCAZ->Fill(pt,(1./pt-1./ptTPC)/(1./pt),dcaZ);
 416       fPtITSminPtTPCvsPtITSPhi->Fill(pt,(pt-ptTPC)/(pt),phi);
 417       fPtITSminPtTPCvsPtITSNPointITS->Fill(pt,(pt-ptTPC)/(pt),nPointITS);
 418       fPtITSminPtTPCvsPtITSNSigmaToVertex->Fill(pt,(1./pt-1./ptTPC)/(1./pt),nSigmaToVertex);
 419       fPtITSminPtTPCvsPtITSChi2C->Fill(pt,(1./pt-1./ptTPC)/(1./pt),chi2C);
 420       fPtITSminPtTPCvsPtITSRel1PtUncertainty->Fill(pt,(1./pt-1./ptTPC)/(1./pt),relUncertainty1Pt);
 421     }//fTrackCutsITS loop
 422
 423   }//ESD track loop
 424
 425   // Post output data
 426   PostData(0, fHistList);
 427   PostData(1, fHistListTPC);
 428   PostData(2, fHistListITS);
 429
 430 }
 431 //________________________________________________________________________
 432 void AliPWG4HighPtQATPConly::Terminate(Option_t *)
 433 {
 434
 435 }
 436
 437 #endif