1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
16 // Class AliHFEitsPIDqa
17 // Monitoring ITS PID in the HFE PID montioring framework. The following
18 // quantities are monitored:
19 // ITS dE/dx (Number of sigmas)
20 // (Always as function of momentum, particle species and centrality
21 // before and after cut)
22 // More information about the PID monitoring framework can be found in
23 // AliHFEpidQAmanager.cxx and AliHFEdetPIDqa.cxx
26 // Yvonne Pachmayer <pachmay@physi.uni-heidelberg.de>
31 #include <THnSparse.h>
34 #include "AliAODTrack.h"
35 #include "AliESDtrack.h"
38 #include "AliPIDResponse.h"
40 #include "AliHFEcollection.h"
41 #include "AliHFEpidBase.h"
42 #include "AliHFEpidQAmanager.h"
43 #include "AliHFEpidITS.h"
44 #include "AliHFEtools.h"
45 #include "AliHFEitsPIDqa.h"
47 ClassImp(AliHFEitsPIDqa)
49 //_________________________________________________________
50 AliHFEitsPIDqa::AliHFEitsPIDqa():
53 , fBrowseCentrality(-1)
60 //_________________________________________________________
61 AliHFEitsPIDqa::AliHFEitsPIDqa(const char* name):
62 AliHFEdetPIDqa(name, "QA for ITS")
64 , fBrowseCentrality(-1)
67 // Default constructor
71 //_________________________________________________________
72 AliHFEitsPIDqa::AliHFEitsPIDqa(const AliHFEitsPIDqa &o):
75 , fBrowseCentrality(o.fBrowseCentrality)
83 //_________________________________________________________
84 AliHFEitsPIDqa &AliHFEitsPIDqa::operator=(const AliHFEitsPIDqa &o){
88 AliHFEdetPIDqa::operator=(o);
89 if(&o != this) o.Copy(*this);
93 //_________________________________________________________
94 AliHFEitsPIDqa::~AliHFEitsPIDqa(){
98 if(fHistos) delete fHistos;
101 //_________________________________________________________
102 void AliHFEitsPIDqa::Copy(TObject &o) const {
106 AliHFEitsPIDqa &target = dynamic_cast<AliHFEitsPIDqa &>(o);
108 delete target.fHistos;
109 target.fHistos = NULL;
111 if(fHistos) target.fHistos = new AliHFEcollection(*fHistos);
112 target.fBrowseCentrality = fBrowseCentrality;
115 //_________________________________________________________
116 Long64_t AliHFEitsPIDqa::Merge(TCollection *coll){
118 // Merge with other objects
121 if(coll->IsEmpty()) return 1;
124 AliHFEitsPIDqa *refQA = NULL;
129 refQA = dynamic_cast<AliHFEitsPIDqa *>(o);
132 listHistos.Add(refQA->fHistos);
135 fHistos->Merge(&listHistos);
139 //_________________________________________________________
140 void AliHFEitsPIDqa::Browse(TBrowser *b){
146 b->Add(fHistos, fHistos->GetName());
148 // Make Projections of the nsigma Spectra and add them to a new Folder
149 TString specnames[AliPID::kSPECIES+1] = {"All", "Electrons", "Muon", "Pions", "Kaon", "Protons"};
150 Int_t specind[AliPID::kSPECIES+1] = {-1, AliPID::kElectron, AliPID::kMuon, AliPID::kPion, AliPID::kKaon, AliPID::kProton};
151 TList *listNsigma = new TList;
152 listNsigma->SetOwner();
155 for(Int_t ispec = 0; ispec < AliPID::kSPECIES+1; ispec++){
156 for(Int_t istep = 0; istep < 2; istep++){
157 hptr = MakeSpectrumNSigma(static_cast<AliHFEdetPIDqa::EStep_t>(istep), specind[ispec], fBrowseCentrality);
158 hptr->SetName(Form("hITSnsigma%s%s%s", specnames[ispec].Data(), istep == 0 ? "Before" : "After", fBrowseCentrality == -1 ? "MinBias" : Form("Cent%d", fBrowseCentrality)));
159 listNsigma->Add(hptr);
163 b->Add(listNsigma, "Projections NSigma");
168 //_________________________________________________________
169 void AliHFEitsPIDqa::Initialize(){
174 fHistos = new AliHFEcollection("itsqahistos", "Collection of ITS QA histograms");
176 // Make common binning
177 const Int_t kNdim = 5;
178 const Int_t kPIDbins = AliPID::kSPECIES + 1;
179 const Int_t kSteps = 2;
180 const Int_t kCentralityBins = 11;
181 const Double_t kMinPID = -1;
182 const Double_t kMinP = 0.;
183 const Double_t kMaxPID = (Double_t)AliPID::kSPECIES;
184 const Double_t kMaxP = 2.;
186 // Quantities where one can switch between low and high resolution
187 Int_t kPbins = fQAmanager->HasHighResolutionHistos() ? 100 : 10;
188 Int_t kSigmaBins = fQAmanager->HasHighResolutionHistos() ? 1400 : 240;
190 // 1st histogram: ITS sigmas: (species, p, nsigma, step, centrality)
191 Int_t nBinsSigma[kNdim] = {kPIDbins, kPbins, kSigmaBins, kSteps, kCentralityBins};
192 Double_t minSigma[kNdim] = {kMinPID, kMinP, -12., 0., 0.};
193 Double_t maxSigma[kNdim] = {kMaxPID, kMaxP, 12., 2., 11.};
194 fHistos->CreateTHnSparse("itsnSigma", "ITS signal; species; p [GeV/c]; ITS signal [a.u.]; Selection Step; Centrality; Eta", kNdim, nBinsSigma, minSigma, maxSigma);
199 //_________________________________________________________
200 void AliHFEitsPIDqa::ProcessTrack(const AliHFEpidObject *track, AliHFEdetPIDqa::EStep_t step){
202 // Fill ITS histograms
204 AliDebug(1, Form("QA started for ITS PID for step %d", (Int_t)step));
205 //AliHFEpidObject::AnalysisType_t anatype = track->IsESDanalysis() ? AliHFEpidObject::kESDanalysis : AliHFEpidObject::kAODanalysis;
206 Float_t centrality = track->GetCentrality();
207 Int_t species = track->GetAbInitioPID();
208 if(species >= AliPID::kSPECIES) species = -1;
210 AliHFEpidITS *itspid = dynamic_cast<AliHFEpidITS *>(fQAmanager->GetDetectorPID(AliHFEpid::kITSpid));
211 Double_t contentSignal[5];
212 contentSignal[0] = species;
213 contentSignal[1] = track->GetRecTrack()->P();
214 contentSignal[2] = itspid ? itspid->GetITSNsigmaCorrected(track->GetRecTrack()) : -100.;
215 contentSignal[3] = step;
216 contentSignal[4] = centrality;
218 fHistos->Fill("itsnSigma", contentSignal);
221 //_________________________________________________________
222 TH2 *AliHFEitsPIDqa::MakeSpectrumNSigma(AliHFEdetPIDqa::EStep_t istep, Int_t species, Int_t centralityClass){
226 THnSparseF *hSignal = dynamic_cast<THnSparseF *>(fHistos->Get("itsnSigma"));
227 if(!hSignal) return NULL;
228 hSignal->GetAxis(3)->SetRange(istep + 1, istep + 1);
229 if(species >= 0 && species < AliPID::kSPECIES)
230 hSignal->GetAxis(0)->SetRange(2 + species, 2 + species);
231 TString hname = Form("hITSsigma%s", istep == AliHFEdetPIDqa::kBeforePID ? "before" : "after"),
232 htitle = Form("ITS dE/dx Spectrum[#sigma] %s selection", istep == AliHFEdetPIDqa::kBeforePID ? "before" : "after");
234 hname += AliPID::ParticleName(species);
235 htitle += Form(" for %ss", AliPID::ParticleName(species));
237 TString centname, centtitle;
238 Bool_t hasCentralityInfo = kTRUE;
239 if(centralityClass > -1){
240 if(hSignal->GetNdimensions() < 5){
241 AliError("Centrality Information not available");
242 centname = centtitle = "MinBias";
243 hasCentralityInfo= kFALSE;
245 // Project centrality classes
247 hSignal->GetAxis(4)->SetRange(centralityClass+1, centralityClass+1);
248 centname = Form("Cent%d", centralityClass);
249 centtitle = Form(" Centrality %d", centralityClass);
252 centname = centtitle = "MinBias";
253 hasCentralityInfo= kFALSE;
258 TH2 *hTmp = hSignal->Projection(2,1);
259 hTmp->SetName(hname.Data());
260 hTmp->SetTitle(htitle.Data());
261 hTmp->SetStats(kFALSE);
262 hTmp->GetXaxis()->SetTitle("p [GeV/c]");
263 hTmp->GetYaxis()->SetTitle("ITS dE/dx - <dE/dx>|_{el} [#sigma]");
264 hSignal->GetAxis(3)->SetRange(0, hSignal->GetAxis(3)->GetNbins());
265 hSignal->GetAxis(0)->SetRange(0, hSignal->GetAxis(0)->GetNbins());
266 if(hasCentralityInfo) hSignal->GetAxis(4)->SetRange(0, hSignal->GetAxis(4)->GetNbins());