[u/mrichter/AliRoot.git] / PWGHF / hfe / AliHFEtofPIDqa.cxx

/**************************************************************************
* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
*                                                                        *
* Author: The ALICE Off-line Project.                                    *
* Contributors are mentioned in the code where appropriate.              *
*                                                                        *
* Permission to use, copy, modify and distribute this software and its   *
* documentation strictly for non-commercial purposes is hereby granted   *
* without fee, provided that the above copyright notice appears in all   *
* copies and that both the copyright notice and this permission notice   *
* appear in the supporting documentation. The authors make no claims     *
* about the suitability of this software for any purpose. It is          *
* provided "as is" without express or implied warranty.                  *
**************************************************************************/
//
// Class AliHFEtofPIDqa
//
// Monitoring TOF PID in the HFE PID montioring framework. The following
// quantities are monitored:
//   TOF time distance to electron hypothesis (Number of sigmas)
//   TOF time distance to the pion hypothesis (Absolute value)
//   TPC dE/dx (Number of sigmas, as control histogram)
// (Always as function of momentum, particle species and centrality 
// before and after cut)
// More information about the PID monitoring framework can be found in
// AliHFEpidQAmanager.cxx and AliHFEdetPIDqa.cxx
//
// Author:
//
//   Markus Fasel <M.Fasel@gsi.de>
#include <TClass.h>
#include <TArrayD.h>
#include <TH2.h>
#include <THnSparse.h>
#include <TString.h>

#include "AliLog.h"
#include "AliPID.h"
#include "AliPIDResponse.h"
#include "AliVParticle.h"

#include "AliHFEcollection.h"
#include "AliHFEpid.h"
#include "AliHFEpidBase.h"
#include "AliHFEpidQAmanager.h"
#include "AliHFEpidTPC.h"
#include "AliHFEpidTRD.h"
#include "AliHFEpidTOF.h"
#include "AliHFEtools.h"
#include "AliHFEtofPIDqa.h"

ClassImp(AliHFEpidTOF)

//_________________________________________________________
AliHFEtofPIDqa::AliHFEtofPIDqa():
    AliHFEdetPIDqa()
  , fHistos(NULL)
{
  //
  // Dummy constructor
  //
}

//_________________________________________________________
AliHFEtofPIDqa::AliHFEtofPIDqa(const char* name):
    AliHFEdetPIDqa(name, "QA for TOF")
  , fHistos(NULL)
{
  //
  // Default constructor
  //
}

//_________________________________________________________
AliHFEtofPIDqa::AliHFEtofPIDqa(const AliHFEtofPIDqa &o):
    AliHFEdetPIDqa(o)
  , fHistos()
{
  //
  // Copy constructor
  //
  o.Copy(*this);
}

//_________________________________________________________
AliHFEtofPIDqa &AliHFEtofPIDqa::operator=(const AliHFEtofPIDqa &o){
  //
  // Do assignment
  //
  AliHFEdetPIDqa::operator=(o);
  if(&o != this) o.Copy(*this);
  return *this;
}

//_________________________________________________________
AliHFEtofPIDqa::~AliHFEtofPIDqa(){
  //
  // Destructor
  //
  if(fHistos) delete fHistos;
}

//_________________________________________________________
void AliHFEtofPIDqa::Copy(TObject &o) const {
  //
  // Make copy
  //
  AliHFEtofPIDqa &target = dynamic_cast<AliHFEtofPIDqa &>(o);
  if(target.fHistos){
    delete target.fHistos;
    target.fHistos = NULL;
  }
  if(fHistos) target.fHistos = new AliHFEcollection(*fHistos);
}

//_________________________________________________________
Long64_t AliHFEtofPIDqa::Merge(TCollection *coll){
  //
  // Merge with other objects
  //
  if(!coll) return 0;
  if(coll->IsEmpty()) return 1;

  TIter it(coll);
  AliHFEtofPIDqa *refQA = NULL;
  TObject *o = NULL;
  Long64_t count = 0;
  TList listHistos;
  while((o = it())){
    refQA = dynamic_cast<AliHFEtofPIDqa *>(o);
    if(!refQA) continue;

    listHistos.Add(refQA->fHistos);
    count++; 
  }
  fHistos->Merge(&listHistos);
  return count + 1;
}

//_________________________________________________________
void AliHFEtofPIDqa::Initialize(){
  //
  // Define Histograms
  //

  fHistos = new AliHFEcollection("tofqahistos", "Collection of TOF QA histograms");

  // Make common binning
  const Int_t kPIDbins = AliPID::kSPECIES + 1;
  const Int_t kSteps = 2;
  const Int_t kCentralityBins = 11;
  const Double_t kMinPID = -1;
  const Double_t kMinP = 0.;
  const Double_t kMaxPID = (Double_t)AliPID::kSPECIES;
  const Double_t kMaxP = 20.;
  // Quantities where one can switch between low and high resolution
  Int_t kPbins = fQAmanager->HasHighResolutionHistos() ?  1000 : 100;
  Int_t kSigmaBins = fQAmanager->HasHighResolutionHistos() ? 1400 : 140;

  // 1st histogram: TOF sigmas: (species, p nsigma, step)
  Int_t nBinsSigma[5] = {kPIDbins, kPbins, kSigmaBins, kSteps, kCentralityBins};
  Double_t minSigma[5] = {kMinPID, kMinP, -12., 0, 0};
  Double_t maxSigma[5] = {kMaxPID, kMaxP, 12., 2., 11.};
  fHistos->CreateTHnSparse("tofnSigma", "TOF signal; species; p [GeV/c]; TOF signal [a.u.]; Selection Step; Centrality", 5, nBinsSigma, minSigma, maxSigma);

  // 3rd histogram: TPC sigmas to the electron line: (species, p nsigma, step - only filled if apriori PID information available)
  fHistos->CreateTHnSparse("tofMonitorTPC", "TPC signal; species; p [GeV/c]; TPC signal [a.u.]; Selection Step; Centrality", 5, nBinsSigma, minSigma, maxSigma);

  // 4th histogram: TOF mismatch template
  AliHFEpidTOF *tofpid = dynamic_cast<AliHFEpidTOF *>(fQAmanager->GetDetectorPID(AliHFEpid::kTOFpid));
  if(!tofpid) AliDebug(1, "TOF PID not available");
  if(tofpid && tofpid->IsGenerateTOFmismatch()){
    AliDebug(1, "Prepare histogram for TOF mismatch tracks");
    fHistos->CreateTHnSparse("tofMismatch", "TOF signal; species; p [GeV/c]; TOF signal [a.u.]; Selection Step; Centrality", 5, nBinsSigma, minSigma, maxSigma);
  }
}

//_________________________________________________________
void AliHFEtofPIDqa::ProcessTrack(const AliHFEpidObject *track, AliHFEdetPIDqa::EStep_t step){
  //
  // Fill TPC histograms
  //
  //AliHFEpidObject::AnalysisType_t anatype = track->IsESDanalysis() ? AliHFEpidObject::kESDanalysis : AliHFEpidObject::kAODanalysis;
  //AliHFEpidObject::AnalysisType_t anatype = track->IsESDanalysis() ? AliHFEpidObject::kESDanalysis : AliHFEpidObject::kAODanalysis;

  Float_t centrality = track->GetCentrality();
  Int_t species = track->GetAbInitioPID();
  if(species >= AliPID::kSPECIES) species = -1;

  AliDebug(1, Form("Monitoring particle of type %d for step %d", species, step));
  AliHFEpidTOF *tofpid= dynamic_cast<AliHFEpidTOF *>(fQAmanager->GetDetectorPID(AliHFEpid::kTOFpid));
  
  const AliPIDResponse *pidResponse = tofpid ? tofpid->GetPIDResponse() : NULL;
  Double_t contentSignal[5];
  contentSignal[0] = species;
  contentSignal[1] = track->GetRecTrack()->P();
  contentSignal[2] = pidResponse ? pidResponse->NumberOfSigmasTOF(track->GetRecTrack(), AliPID::kElectron): 0.;
  contentSignal[3] = step;
  contentSignal[4] = centrality;
  fHistos->Fill("tofnSigma", contentSignal);
  contentSignal[2] = pidResponse ? pidResponse->NumberOfSigmasTPC(track->GetRecTrack(), AliPID::kElectron) : 0.;
  fHistos->Fill("tofMonitorTPC", contentSignal);

  if(tofpid->IsGenerateTOFmismatch()){
    AliDebug(1,"Filling TOF mismatch histos");
    TArrayD nsigmismatch(tofpid->GetNmismatchTracks());
    tofpid->GenerateTOFmismatch(track->GetRecTrack(),tofpid->GetNmismatchTracks(), nsigmismatch);
    for(int itrk = 0; itrk < tofpid->GetNmismatchTracks(); itrk++){
      contentSignal[2] = nsigmismatch[itrk];
      fHistos->Fill("tofMismatch",contentSignal);
    }
  }
}

//_________________________________________________________
TH2 *AliHFEtofPIDqa::MakeTPCspectrumNsigma(AliHFEdetPIDqa::EStep_t step, Int_t species, Int_t centralityClass){
  //
  // Get the TPC control histogram for the TRD selection step (either before or after PID)
  //
  THnSparseF *histo = dynamic_cast<THnSparseF *>(fHistos->Get("tofMonitorTPC"));
  if(!histo){
    AliError("QA histogram monitoring TPC nSigma not available");
    return NULL;
  }
  if(species > -1 && species < AliPID::kSPECIES){
    // cut on species (if available)
    histo->GetAxis(0)->SetRange(species + 2, species + 2); // undef + underflow
  }
  TString centname, centtitle;
  Bool_t hasCentralityInfo = kTRUE;
  if(centralityClass > -1){
    if(histo->GetNdimensions() < 5){
      AliError("Centrality Information not available");
      centname = centtitle = "MinBias";
      hasCentralityInfo = kFALSE;
    } else {
      // Project centrality classes
      // -1 is Min. Bias
      histo->GetAxis(4)->SetRange(centralityClass+1, centralityClass+1);
      centname = Form("Cent%d", centralityClass);
      centtitle = Form("Centrality %d", centralityClass);
    }
  } else {
    histo->GetAxis(4)->SetRange(1, histo->GetAxis(4)->GetNbins()-1);
    centname = centtitle = "MinBias";
    hasCentralityInfo = kTRUE;
  }
  histo->GetAxis(3)->SetRange(step + 1, step + 1); 

  TH2 *hSpec = histo->Projection(2, 1);
  // construct title and name
  TString stepname = step == AliHFEdetPIDqa::kBeforePID ? "before" : "after";
  TString speciesname = species > -1 && species < AliPID::kSPECIES ? AliPID::ParticleName(species) : "all Particles";
  TString specID = species > -1 && species < AliPID::kSPECIES ? AliPID::ParticleName(species) : "unid";
  TString histname = Form("hSigmaTPC%s%s%s", specID.Data(), stepname.Data(), centname.Data());
  TString histtitle = Form("TPC Sigma for %s %s PID %s", speciesname.Data(), stepname.Data(), centtitle.Data());
  hSpec->SetName(histname.Data());
  hSpec->SetTitle(histtitle.Data());

  // Unset range on the original histogram
  histo->GetAxis(0)->SetRange(0, histo->GetAxis(0)->GetNbins());
  histo->GetAxis(3)->SetRange(0, histo->GetAxis(3)->GetNbins());
  if(hasCentralityInfo)histo->GetAxis(4)->SetRange(0, histo->GetAxis(4)->GetNbins());
  return hSpec; 
}

//_________________________________________________________
TH2 *AliHFEtofPIDqa::MakeSpectrumNSigma(AliHFEdetPIDqa::EStep_t istep, Int_t species, Int_t centralityClass){
  //
  // Plot the Spectrum
  //
  THnSparseF *hSignal = dynamic_cast<THnSparseF *>(fHistos->Get("tofnSigma"));
  if(!hSignal) return NULL;
  hSignal->GetAxis(3)->SetRange(istep + 1, istep + 1);
  if(species >= 0 && species < AliPID::kSPECIES)
    hSignal->GetAxis(0)->SetRange(2 + species, 2 + species);
  TString centname, centtitle;
  Bool_t hasCentralityInfo = kTRUE;
  if(centralityClass > -1){
    if(hSignal->GetNdimensions() < 5){
      AliError("Centrality Information not available");
      centname = centtitle = "MinBias";
      hasCentralityInfo = kFALSE;
    } else {
      // Project centrality classes
      // -1 is Min. Bias
      AliInfo(Form("Centrality Class: %d", centralityClass));
      hSignal->GetAxis(4)->SetRange(centralityClass+1, centralityClass+1);
      centname = Form("Cent%d", centralityClass);
      centtitle = Form("Centrality %d", centralityClass);
    }
  } else {
    hSignal->GetAxis(4)->SetRange(1, hSignal->GetAxis(4)->GetNbins()-1);
    centname = centtitle = "MinBias";
    hasCentralityInfo = kTRUE;
  }
  TH2 *hTmp = hSignal->Projection(2,1);
  TString hname = Form("hTOFsigma%s%s", istep == AliHFEdetPIDqa::kBeforePID ? "before" : "after", centname.Data()), 
          htitle = Form("Time-of-flight spectrum[#sigma] %s selection %s", istep == AliHFEdetPIDqa::kBeforePID ? "before" : "after", centtitle.Data());
  if(species > -1){
    hname += AliPID::ParticleName(species);
    htitle += Form(" for %ss", AliPID::ParticleName(species));
  }
  hTmp->SetName(hname.Data());
  hTmp->SetTitle(htitle.Data());
  hTmp->SetStats(kFALSE);
  hTmp->GetXaxis()->SetTitle("p [GeV/c]");
  hTmp->GetYaxis()->SetTitle("TOF time|_{el} - expected time|_{el} [#sigma]");
  hSignal->GetAxis(3)->SetRange(0, hSignal->GetAxis(3)->GetNbins());
  hSignal->GetAxis(0)->SetRange(0, hSignal->GetAxis(0)->GetNbins());
  if(hasCentralityInfo)hSignal->GetAxis(4)->SetRange(0, hSignal->GetAxis(4)->GetNbins());
  return hTmp;
}

//_________________________________________________________
TH1 *AliHFEtofPIDqa::GetHistogram(const char *name){
  // 
  // Get the histogram
  //
  if(!fHistos) return NULL;
  TObject *histo = fHistos->Get(name);
  if(!histo->InheritsFrom("TH1")) return NULL;
  return dynamic_cast<TH1 *>(histo);
}
Commit	Line	Data
3a72645a	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	// Class AliHFEtofPIDqa
	17	//
	18	// Monitoring TOF PID in the HFE PID montioring framework. The following
	19	// quantities are monitored:
	20	// TOF time distance to electron hypothesis (Number of sigmas)
	21	// TOF time distance to the pion hypothesis (Absolute value)
	22	// TPC dE/dx (Number of sigmas, as control histogram)
	23	// (Always as function of momentum, particle species and centrality
	24	// before and after cut)
	25	// More information about the PID monitoring framework can be found in
	26	// AliHFEpidQAmanager.cxx and AliHFEdetPIDqa.cxx
	27	//
	28	// Author:
	29	//
	30	// Markus Fasel <M.Fasel@gsi.de>
	31	#include <TClass.h>
afb48e1d	32	#include <TArrayD.h>
3a72645a	33	#include <TH2.h>
	34	#include <THnSparse.h>
	35	#include <TString.h>
	36
3a72645a	37	#include "AliLog.h"
3a72645a	38	#include "AliPID.h"
8c1c76e9	39	#include "AliPIDResponse.h"
6555e2ad	40	#include "AliVParticle.h"
3a72645a	41
3a72645a	42	#include "AliHFEcollection.h"
6555e2ad	43	#include "AliHFEpid.h"
3a72645a	44	#include "AliHFEpidBase.h"
6555e2ad	45	#include "AliHFEpidQAmanager.h"
6555e2ad	46	#include "AliHFEpidTPC.h"
cedf0381	47	#include "AliHFEpidTRD.h"
6555e2ad	48	#include "AliHFEpidTOF.h"
3a72645a	49	#include "AliHFEtools.h"
	50	#include "AliHFEtofPIDqa.h"
	51
6555e2ad	52	ClassImp(AliHFEpidTOF)
6555e2ad	53
3a72645a	54	//_________________________________________________________
	55	AliHFEtofPIDqa::AliHFEtofPIDqa():
	56	AliHFEdetPIDqa()
	57	, fHistos(NULL)
	58	{
	59	//
	60	// Dummy constructor
	61	//
	62	}
	63
	64	//_________________________________________________________
	65	AliHFEtofPIDqa::AliHFEtofPIDqa(const char* name):
6555e2ad	66	AliHFEdetPIDqa(name, "QA for TOF")
3a72645a	67	, fHistos(NULL)
	68	{
	69	//
	70	// Default constructor
	71	//
	72	}
	73
	74	//_________________________________________________________
	75	AliHFEtofPIDqa::AliHFEtofPIDqa(const AliHFEtofPIDqa &o):
	76	AliHFEdetPIDqa(o)
	77	, fHistos()
	78	{
	79	//
	80	// Copy constructor
	81	//
	82	o.Copy(*this);
	83	}
	84
	85	//_________________________________________________________
	86	AliHFEtofPIDqa &AliHFEtofPIDqa::operator=(const AliHFEtofPIDqa &o){
	87	//
	88	// Do assignment
	89	//
	90	AliHFEdetPIDqa::operator=(o);
	91	if(&o != this) o.Copy(*this);
	92	return *this;
	93	}
	94
	95	//_________________________________________________________
	96	AliHFEtofPIDqa::~AliHFEtofPIDqa(){
	97	//
	98	// Destructor
	99	//
	100	if(fHistos) delete fHistos;
	101	}
	102
	103	//_________________________________________________________
	104	void AliHFEtofPIDqa::Copy(TObject &o) const {
	105	//
	106	// Make copy
	107	//
	108	AliHFEtofPIDqa &target = dynamic_cast<AliHFEtofPIDqa &>(o);
	109	if(target.fHistos){
	110	delete target.fHistos;
	111	target.fHistos = NULL;
	112	}
	113	if(fHistos) target.fHistos = new AliHFEcollection(*fHistos);
	114	}
	115
	116	//_________________________________________________________
	117	Long64_t AliHFEtofPIDqa::Merge(TCollection *coll){
	118	//
	119	// Merge with other objects
	120	//
	121	if(!coll) return 0;
	122	if(coll->IsEmpty()) return 1;
	123
	124	TIter it(coll);
	125	AliHFEtofPIDqa *refQA = NULL;
	126	TObject *o = NULL;
	127	Long64_t count = 0;
	128	TList listHistos;
	129	while((o = it())){
	130	refQA = dynamic_cast<AliHFEtofPIDqa *>(o);
131	if(!refQA) continue;
132
133	listHistos.Add(refQA->fHistos);
134	count++;
135	}
136	fHistos->Merge(&listHistos);
137	return count + 1;
138	}
139
140	//_________________________________________________________
141	void AliHFEtofPIDqa::Initialize(){
142	//
143	// Define Histograms
144	//
145
146	fHistos = new AliHFEcollection("tofqahistos", "Collection of TOF QA histograms");
147
148	// Make common binning
149	const Int_t kPIDbins = AliPID::kSPECIES + 1;
3a72645a	150	const Int_t kSteps = 2;
c2690925	151	const Int_t kCentralityBins = 11;
3a72645a	152	const Double_t kMinPID = -1;
	153	const Double_t kMinP = 0.;
	154	const Double_t kMaxPID = (Double_t)AliPID::kSPECIES;
	155	const Double_t kMaxP = 20.;
bf892a6a	156	// Quantities where one can switch between low and high resolution
bf892a6a	157	Int_t kPbins = fQAmanager->HasHighResolutionHistos() ? 1000 : 100;
cedf0381	158	Int_t kSigmaBins = fQAmanager->HasHighResolutionHistos() ? 1400 : 140;
bf892a6a	159
3a72645a	160	// 1st histogram: TOF sigmas: (species, p nsigma, step)
c2690925	161	Int_t nBinsSigma[5] = {kPIDbins, kPbins, kSigmaBins, kSteps, kCentralityBins};
	162	Double_t minSigma[5] = {kMinPID, kMinP, -12., 0, 0};
	163	Double_t maxSigma[5] = {kMaxPID, kMaxP, 12., 2., 11.};
	164	fHistos->CreateTHnSparse("tofnSigma", "TOF signal; species; p [GeV/c]; TOF signal [a.u.]; Selection Step; Centrality", 5, nBinsSigma, minSigma, maxSigma);
3a72645a	165
3a72645a	166	// 3rd histogram: TPC sigmas to the electron line: (species, p nsigma, step - only filled if apriori PID information available)
c2690925	167	fHistos->CreateTHnSparse("tofMonitorTPC", "TPC signal; species; p [GeV/c]; TPC signal [a.u.]; Selection Step; Centrality", 5, nBinsSigma, minSigma, maxSigma);
afb48e1d	168
	169	// 4th histogram: TOF mismatch template
	170	AliHFEpidTOF tofpid = dynamic_cast<AliHFEpidTOF >(fQAmanager->GetDetectorPID(AliHFEpid::kTOFpid));
	171	if(!tofpid) AliDebug(1, "TOF PID not available");
	172	if(tofpid && tofpid->IsGenerateTOFmismatch()){
	173	AliDebug(1, "Prepare histogram for TOF mismatch tracks");
	174	fHistos->CreateTHnSparse("tofMismatch", "TOF signal; species; p [GeV/c]; TOF signal [a.u.]; Selection Step; Centrality", 5, nBinsSigma, minSigma, maxSigma);
	175	}
3a72645a	176	}
	177
	178	//_________________________________________________________
6555e2ad	179	void AliHFEtofPIDqa::ProcessTrack(const AliHFEpidObject *track, AliHFEdetPIDqa::EStep_t step){
3a72645a	180	//
	181	// Fill TPC histograms
	182	//
8c1c76e9	183	//AliHFEpidObject::AnalysisType_t anatype = track->IsESDanalysis() ? AliHFEpidObject::kESDanalysis : AliHFEpidObject::kAODanalysis;
afb48e1d	184	//AliHFEpidObject::AnalysisType_t anatype = track->IsESDanalysis() ? AliHFEpidObject::kESDanalysis : AliHFEpidObject::kAODanalysis;
cedf0381	185
c2690925	186	Float_t centrality = track->GetCentrality();
3a72645a	187	Int_t species = track->GetAbInitioPID();
3a72645a	188	if(species >= AliPID::kSPECIES) species = -1;
3a72645a	189
3a72645a	190	AliDebug(1, Form("Monitoring particle of type %d for step %d", species, step));
6555e2ad	191	AliHFEpidTOF tofpid= dynamic_cast<AliHFEpidTOF >(fQAmanager->GetDetectorPID(AliHFEpid::kTOFpid));
3a72645a	192
8c1c76e9	193	const AliPIDResponse *pidResponse = tofpid ? tofpid->GetPIDResponse() : NULL;
c2690925	194	Double_t contentSignal[5];
3a72645a	195	contentSignal[0] = species;
6555e2ad	196	contentSignal[1] = track->GetRecTrack()->P();
8c1c76e9	197	contentSignal[2] = pidResponse ? pidResponse->NumberOfSigmasTOF(track->GetRecTrack(), AliPID::kElectron): 0.;
3a72645a	198	contentSignal[3] = step;
c2690925	199	contentSignal[4] = centrality;
bf892a6a	200	fHistos->Fill("tofnSigma", contentSignal);
cedf0381	201	contentSignal[2] = pidResponse ? pidResponse->NumberOfSigmasTPC(track->GetRecTrack(), AliPID::kElectron) : 0.;
afb48e1d	202	fHistos->Fill("tofMonitorTPC", contentSignal);
	203
	204	if(tofpid->IsGenerateTOFmismatch()){
	205	AliDebug(1,"Filling TOF mismatch histos");
	206	TArrayD nsigmismatch(tofpid->GetNmismatchTracks());
	207	tofpid->GenerateTOFmismatch(track->GetRecTrack(),tofpid->GetNmismatchTracks(), nsigmismatch);
	208	for(int itrk = 0; itrk < tofpid->GetNmismatchTracks(); itrk++){
	209	contentSignal[2] = nsigmismatch[itrk];
	210	fHistos->Fill("tofMismatch",contentSignal);
	211	}
	212	}
cedf0381	213	}
	214
	215	//_________________________________________________________
	216	TH2 *AliHFEtofPIDqa::MakeTPCspectrumNsigma(AliHFEdetPIDqa::EStep_t step, Int_t species, Int_t centralityClass){
	217	//
	218	// Get the TPC control histogram for the TRD selection step (either before or after PID)
	219	//
	220	THnSparseF histo = dynamic_cast<THnSparseF >(fHistos->Get("tofMonitorTPC"));
	221	if(!histo){
	222	AliError("QA histogram monitoring TPC nSigma not available");
	223	return NULL;
	224	}
	225	if(species > -1 && species < AliPID::kSPECIES){
	226	// cut on species (if available)
	227	histo->GetAxis(0)->SetRange(species + 2, species + 2); // undef + underflow
	228	}
	229	TString centname, centtitle;
	230	Bool_t hasCentralityInfo = kTRUE;
	231	if(centralityClass > -1){
	232	if(histo->GetNdimensions() < 5){
	233	AliError("Centrality Information not available");
	234	centname = centtitle = "MinBias";
	235	hasCentralityInfo = kFALSE;
	236	} else {
	237	// Project centrality classes
	238	// -1 is Min. Bias
	239	histo->GetAxis(4)->SetRange(centralityClass+1, centralityClass+1);
	240	centname = Form("Cent%d", centralityClass);
	241	centtitle = Form("Centrality %d", centralityClass);
	242	}
	243	} else {
	244	histo->GetAxis(4)->SetRange(1, histo->GetAxis(4)->GetNbins()-1);
	245	centname = centtitle = "MinBias";
	246	hasCentralityInfo = kTRUE;
3a72645a	247	}
cedf0381	248	histo->GetAxis(3)->SetRange(step + 1, step + 1);
	249
	250	TH2 *hSpec = histo->Projection(2, 1);
	251	// construct title and name
	252	TString stepname = step == AliHFEdetPIDqa::kBeforePID ? "before" : "after";
	253	TString speciesname = species > -1 && species < AliPID::kSPECIES ? AliPID::ParticleName(species) : "all Particles";
	254	TString specID = species > -1 && species < AliPID::kSPECIES ? AliPID::ParticleName(species) : "unid";
	255	TString histname = Form("hSigmaTPC%s%s%s", specID.Data(), stepname.Data(), centname.Data());
	256	TString histtitle = Form("TPC Sigma for %s %s PID %s", speciesname.Data(), stepname.Data(), centtitle.Data());
	257	hSpec->SetName(histname.Data());
	258	hSpec->SetTitle(histtitle.Data());
	259
	260	// Unset range on the original histogram
	261	histo->GetAxis(0)->SetRange(0, histo->GetAxis(0)->GetNbins());
	262	histo->GetAxis(3)->SetRange(0, histo->GetAxis(3)->GetNbins());
	263	if(hasCentralityInfo)histo->GetAxis(4)->SetRange(0, histo->GetAxis(4)->GetNbins());
	264	return hSpec;
3a72645a	265	}
3a72645a	266
3a72645a	267	//_________________________________________________________
4437a0d2	268	TH2 *AliHFEtofPIDqa::MakeSpectrumNSigma(AliHFEdetPIDqa::EStep_t istep, Int_t species, Int_t centralityClass){
3a72645a	269	//
	270	// Plot the Spectrum
	271	//
	272	THnSparseF hSignal = dynamic_cast<THnSparseF >(fHistos->Get("tofnSigma"));
bf892a6a	273	if(!hSignal) return NULL;
3a72645a	274	hSignal->GetAxis(3)->SetRange(istep + 1, istep + 1);
e3ae862b	275	if(species >= 0 && species < AliPID::kSPECIES)
3a72645a	276	hSignal->GetAxis(0)->SetRange(2 + species, 2 + species);
4437a0d2	277	TString centname, centtitle;
	278	Bool_t hasCentralityInfo = kTRUE;
	279	if(centralityClass > -1){
	280	if(hSignal->GetNdimensions() < 5){
	281	AliError("Centrality Information not available");
	282	centname = centtitle = "MinBias";
	283	hasCentralityInfo = kFALSE;
	284	} else {
	285	// Project centrality classes
	286	// -1 is Min. Bias
	287	AliInfo(Form("Centrality Class: %d", centralityClass));
	288	hSignal->GetAxis(4)->SetRange(centralityClass+1, centralityClass+1);
	289	centname = Form("Cent%d", centralityClass);
	290	centtitle = Form("Centrality %d", centralityClass);
	291	}
	292	} else {
	293	hSignal->GetAxis(4)->SetRange(1, hSignal->GetAxis(4)->GetNbins()-1);
	294	centname = centtitle = "MinBias";
	295	hasCentralityInfo = kTRUE;
	296	}
3a72645a	297	TH2 *hTmp = hSignal->Projection(2,1);
4437a0d2	298	TString hname = Form("hTOFsigma%s%s", istep == AliHFEdetPIDqa::kBeforePID ? "before" : "after", centname.Data()),
4437a0d2	299	htitle = Form("Time-of-flight spectrum[#sigma] %s selection %s", istep == AliHFEdetPIDqa::kBeforePID ? "before" : "after", centtitle.Data());
3a72645a	300	if(species > -1){
bf892a6a	301	hname += AliPID::ParticleName(species);
bf892a6a	302	htitle += Form(" for %ss", AliPID::ParticleName(species));
3a72645a	303	}
bf892a6a	304	hTmp->SetName(hname.Data());
bf892a6a	305	hTmp->SetTitle(htitle.Data());
3a72645a	306	hTmp->SetStats(kFALSE);
	307	hTmp->GetXaxis()->SetTitle("p [GeV/c]");
	308	hTmp->GetYaxis()->SetTitle("TOF time\|_{el} - expected time\|_{el} [#sigma]");
	309	hSignal->GetAxis(3)->SetRange(0, hSignal->GetAxis(3)->GetNbins());
	310	hSignal->GetAxis(0)->SetRange(0, hSignal->GetAxis(0)->GetNbins());
4437a0d2	311	if(hasCentralityInfo)hSignal->GetAxis(4)->SetRange(0, hSignal->GetAxis(4)->GetNbins());
3a72645a	312	return hTmp;
	313	}
	314
	315	//_________________________________________________________
	316	TH1 AliHFEtofPIDqa::GetHistogram(const char name){
	317	//
	318	// Get the histogram
	319	//
	320	if(!fHistos) return NULL;
	321	TObject *histo = fHistos->Get(name);
	322	if(!histo->InheritsFrom("TH1")) return NULL;
	323	return dynamic_cast<TH1 *>(histo);
	324	}
	325