+/**************************************************************************
+ * 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. *
+ **************************************************************************/
+
/* $Id$ */
#include "AlidNdEtaCorrection.h"
+#include <AliLog.h>
#include <TCanvas.h>
+#include <TH3F.h>
#include <TH2F.h>
+#include <TH1D.h>
+#include <TDirectory.h>
+#include <AliCorrection.h>
+#include <AliCorrectionMatrix2D.h>
+#include <AliCorrectionMatrix3D.h>
//____________________________________________________________________
ClassImp(AlidNdEtaCorrection)
//____________________________________________________________________
-AlidNdEtaCorrection::AlidNdEtaCorrection(Char_t* name)
- : TNamed(name, name),
- fNEvents(0),
- fNTriggeredEvents(0)
-{
+AlidNdEtaCorrection::AlidNdEtaCorrection()
+ : TNamed(),
+ fTrack2ParticleCorrection(0),
+ fVertexRecoCorrection(0),
+ fTriggerBiasCorrectionMBToINEL(0),
+ fTriggerBiasCorrectionMBToNSD(0),
+ fTriggerBiasCorrectionMBToND(0)
+{
+ // default constructor
+}
+
+//____________________________________________________________________
+AlidNdEtaCorrection::AlidNdEtaCorrection(const Char_t* name, const Char_t* title, AliPWG0Helper::AnalysisMode analysis)
+ : TNamed(name, title),
+ fTrack2ParticleCorrection(0),
+ fVertexRecoCorrection(0),
+ fTriggerBiasCorrectionMBToINEL(0),
+ fTriggerBiasCorrectionMBToNSD(0),
+ fTriggerBiasCorrectionMBToND(0)
+{
+ //
// constructor
//
- fTrack2ParticleCorrection = new AliCorrectionMatrix3D("nTrackToNPart", "nTrackToNPart",80,-20,20,120,-6,6, 100, 0, 10);
+ fTrack2ParticleCorrection = new AliCorrection("Track2Particle", "Track2Particle", analysis);
+ fVertexRecoCorrection = new AliCorrection("VertexReconstruction", "VertexReconstruction", analysis);
- Float_t binLimitsN[] = {-0.5, 0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5,
- 10.5, 12.5, 14.5, 16.5, 18.5, 20.5, 25.5, 30.5, 40.5, 50.5, 100.5, 300.5};
- Float_t binLimitsVtx[] = {-20,-15,-10,-6,-3,0,3,6,10,15,20};
-
- fVertexRecoCorrection = new AliCorrectionMatrix2D("vtxReco", "vtxReco",10,binLimitsVtx ,22,binLimitsN);
+ fTriggerBiasCorrectionMBToINEL = new AliCorrection("TriggerBias_MBToINEL", "TriggerBias_MBToINEL", analysis);
+ fTriggerBiasCorrectionMBToNSD = new AliCorrection("TriggerBias_MBToNSD", "TriggerBias_MBToNSD", analysis);
+ fTriggerBiasCorrectionMBToND = new AliCorrection("TriggerBias_MBToND", "TriggerBias_MBToND", analysis);
+}
+
+//____________________________________________________________________
+AlidNdEtaCorrection::~AlidNdEtaCorrection()
+{
+ // destructor
+
+ if (fTrack2ParticleCorrection) {
+ delete fTrack2ParticleCorrection;
+ fTrack2ParticleCorrection = 0;
+ }
+
+ if (fVertexRecoCorrection) {
+ delete fVertexRecoCorrection;
+ fVertexRecoCorrection = 0;
+ }
- fTriggerBiasCorrection = new AliCorrectionMatrix2D("triggerBias", "triggerBias",120,-6,6,100, 0, 10);
+ if (fTriggerBiasCorrectionMBToINEL) {
+ delete fTriggerBiasCorrectionMBToINEL;
+ fTriggerBiasCorrectionMBToINEL = 0;
+ }
- fTrack2ParticleCorrection ->SetAxisTitles("vtx z [cm]", "#eta", "p_{T}");
- fVertexRecoCorrection ->SetAxisTitles("vtx z [cm]", "n particles/tracks/tracklets?");
+ if (fTriggerBiasCorrectionMBToNSD) {
+ delete fTriggerBiasCorrectionMBToNSD;
+ fTriggerBiasCorrectionMBToNSD = 0;
+ }
- fTriggerBiasCorrection ->SetAxisTitles("#eta", "p_{T} [GeV/c]");
+ if (fTriggerBiasCorrectionMBToND) {
+ delete fTriggerBiasCorrectionMBToND;
+ fTriggerBiasCorrectionMBToND = 0;
+ }
}
//____________________________________________________________________
//
// divide the histograms in the AliCorrectionMatrix2D objects to get the corrections
-
fTrack2ParticleCorrection->Divide();
-
fVertexRecoCorrection->Divide();
-
- fTriggerBiasCorrection->GetMeasuredHistogram()->Scale(Double_t(fNTriggeredEvents)/Double_t(fNEvents));
- fTriggerBiasCorrection->Divide();
-
+ fTriggerBiasCorrectionMBToINEL->Divide();
+ fTriggerBiasCorrectionMBToNSD->Divide();
+ fTriggerBiasCorrectionMBToND->Divide();
}
//____________________________________________________________________
-Long64_t
-AlidNdEtaCorrection::Merge(TCollection* list) {
+Long64_t AlidNdEtaCorrection::Merge(TCollection* list)
+{
// Merge a list of dNdEtaCorrection objects with this (needed for
- // PROOF).
+ // PROOF).
// Returns the number of merged objects (including this).
if (!list)
return 0;
-
+
if (list->IsEmpty())
return 1;
TObject* obj;
// collections of measured and generated histograms
- TList* collectionNtrackToNparticle = new TList;
- TList* collectionVertexReco = new TList;
- TList* collectionTriggerBias = new TList;
+ TList* collectionNtrackToNparticle = new TList;
+ TList* collectionVertexReco = new TList;
+ TList* collectionTriggerBiasMBToINEL = new TList;
+ TList* collectionTriggerBiasMBToNSD = new TList;
+ TList* collectionTriggerBiasMBToND = new TList;
Int_t count = 0;
while ((obj = iter->Next())) {
-
+
AlidNdEtaCorrection* entry = dynamic_cast<AlidNdEtaCorrection*> (obj);
- if (entry == 0)
+ if (entry == 0)
continue;
- collectionNtrackToNparticle ->Add(entry->GetTrack2ParticleCorrection());
- collectionVertexReco ->Add(entry->GetVertexRecoCorrection());
- collectionTriggerBias ->Add(entry->GetTriggerBiasCorrection());
+ collectionNtrackToNparticle ->Add(entry->fTrack2ParticleCorrection);
+ collectionVertexReco ->Add(entry->fVertexRecoCorrection);
+ collectionTriggerBiasMBToINEL->Add(entry->fTriggerBiasCorrectionMBToINEL);
+ collectionTriggerBiasMBToNSD ->Add(entry->fTriggerBiasCorrectionMBToNSD);
+ collectionTriggerBiasMBToND ->Add(entry->fTriggerBiasCorrectionMBToND);
count++;
}
- fTrack2ParticleCorrection ->Merge(collectionNtrackToNparticle);
- fVertexRecoCorrection ->Merge(collectionVertexReco);
- fTriggerBiasCorrection ->Merge(collectionTriggerBias);
-
+ fTrack2ParticleCorrection ->Merge(collectionNtrackToNparticle);
+ fVertexRecoCorrection ->Merge(collectionVertexReco);
+ fTriggerBiasCorrectionMBToINEL ->Merge(collectionTriggerBiasMBToINEL);
+ fTriggerBiasCorrectionMBToNSD ->Merge(collectionTriggerBiasMBToNSD);
+ fTriggerBiasCorrectionMBToND ->Merge(collectionTriggerBiasMBToND);
+
delete collectionNtrackToNparticle;
delete collectionVertexReco;
- delete collectionTriggerBias;
+ delete collectionTriggerBiasMBToINEL;
+ delete collectionTriggerBiasMBToNSD;
+ delete collectionTriggerBiasMBToND;
return count+1;
}
+//____________________________________________________________________
+void AlidNdEtaCorrection::Add(AlidNdEtaCorrection* aCorrectionsToAdd, Float_t c) {
+ //
+ // adds the measured and generated of aCorrectionsToAdd to measured and generated
+ // of all corrections in this
+
+ fTrack2ParticleCorrection ->Add(aCorrectionsToAdd->GetTrack2ParticleCorrection() ,c);
+ fVertexRecoCorrection ->Add(aCorrectionsToAdd->GetVertexRecoCorrection() ,c);
+ fTriggerBiasCorrectionMBToINEL ->Add(aCorrectionsToAdd->GetTriggerBiasCorrectionINEL(),c);
+ fTriggerBiasCorrectionMBToNSD ->Add(aCorrectionsToAdd->GetTriggerBiasCorrectionNSD() ,c);
+ fTriggerBiasCorrectionMBToND ->Add(aCorrectionsToAdd->GetTriggerBiasCorrectionND() ,c);
+
+}
+
+//____________________________________________________________________
+void AlidNdEtaCorrection::Reset(void) {
+ //
+ // reset all corrections
+ //
+
+ fTrack2ParticleCorrection ->Reset();
+ fVertexRecoCorrection ->Reset();
+ fTriggerBiasCorrectionMBToINEL ->Reset();
+ fTriggerBiasCorrectionMBToNSD ->Reset();
+ fTriggerBiasCorrectionMBToND ->Reset();
+
+}
+
+
//____________________________________________________________________
-Bool_t
-AlidNdEtaCorrection::LoadHistograms(Char_t* fileName, Char_t* dir) {
+Bool_t AlidNdEtaCorrection::LoadHistograms(const Char_t* dir)
+{
//
// loads the histograms
+ // if dir is empty a directory with the name of this object is taken (like in SaveHistogram)
//
- fTrack2ParticleCorrection ->LoadHistograms(fileName, dir);
- fVertexRecoCorrection ->LoadHistograms(fileName, dir);
- fTriggerBiasCorrection ->LoadHistograms(fileName, dir);
-
+ if (!dir)
+ dir = GetName();
+
+ if (!gDirectory->cd(dir))
+ return kFALSE;
+
+ fTrack2ParticleCorrection ->LoadHistograms();
+ fVertexRecoCorrection ->LoadHistograms();
+ fTriggerBiasCorrectionMBToINEL ->LoadHistograms();
+ fTriggerBiasCorrectionMBToNSD ->LoadHistograms();
+ fTriggerBiasCorrectionMBToND ->LoadHistograms();
+
+ gDirectory->cd("..");
+
return kTRUE;
}
-
//____________________________________________________________________
-void
-AlidNdEtaCorrection::SaveHistograms() {
+void AlidNdEtaCorrection::SaveHistograms()
+{
//
// save the histograms
//
gDirectory->mkdir(fName.Data());
gDirectory->cd(fName.Data());
- fTrack2ParticleCorrection ->SaveHistograms();
- fVertexRecoCorrection ->SaveHistograms();
- fTriggerBiasCorrection ->SaveHistograms();
+ fTrack2ParticleCorrection ->SaveHistograms();
+ fVertexRecoCorrection ->SaveHistograms();
+ fTriggerBiasCorrectionMBToINEL->SaveHistograms();
+ fTriggerBiasCorrectionMBToNSD ->SaveHistograms();
+ fTriggerBiasCorrectionMBToND ->SaveHistograms();
- gDirectory->cd("../");
+ gDirectory->cd("..");
}
//____________________________________________________________________
void AlidNdEtaCorrection::DrawHistograms()
{
//
- // call the draw histogram method of the two AliCorrectionMatrix2D objects
+ // call the draw histogram method of the corrections
+ //
+
+ fTrack2ParticleCorrection ->DrawHistograms();
+ fVertexRecoCorrection ->DrawHistograms();
+ fTriggerBiasCorrectionMBToINEL->DrawHistograms();
+ fTriggerBiasCorrectionMBToNSD ->DrawHistograms();
+ fTriggerBiasCorrectionMBToND ->DrawHistograms();
+}
+
+//____________________________________________________________________
+void AlidNdEtaCorrection::DrawOverview(const char* canvasName)
+{
+ //
+ // call the DrawOverview histogram method of the corrections
+ //
+
+ fTrack2ParticleCorrection ->DrawOverview(canvasName);
+ fVertexRecoCorrection ->DrawOverview(canvasName);
+ fTriggerBiasCorrectionMBToINEL->DrawOverview(canvasName);
+ fTriggerBiasCorrectionMBToNSD ->DrawOverview(canvasName);
+ fTriggerBiasCorrectionMBToND ->DrawOverview(canvasName);
+}
+
+//____________________________________________________________________
+void AlidNdEtaCorrection::FillMCParticle(Float_t vtx, Float_t eta, Float_t pt, Bool_t trigger, Bool_t vertex, Int_t processType)
+{
+ // fills a particle in the corrections
+ // it is filled in generated or measured depending of the flags
- fTrack2ParticleCorrection ->DrawHistograms();
- fVertexRecoCorrection ->DrawHistograms();
- fTriggerBiasCorrection ->DrawHistograms();
+ fTriggerBiasCorrectionMBToINEL->GetTrackCorrection()->FillGene(vtx, eta, pt);
+
+ if (processType != 92 && processType != 93)
+ fTriggerBiasCorrectionMBToNSD->GetTrackCorrection()->FillGene(vtx, eta, pt);
+
+ if (processType!=92 && processType!=93 && processType!=94)
+ fTriggerBiasCorrectionMBToND->GetTrackCorrection()->FillGene(vtx, eta, pt);
+
+ if (!trigger)
+ return;
+
+ fTriggerBiasCorrectionMBToINEL->GetTrackCorrection()->FillMeas(vtx, eta, pt);
+ fTriggerBiasCorrectionMBToNSD->GetTrackCorrection()->FillMeas(vtx, eta, pt);
+ fTriggerBiasCorrectionMBToND->GetTrackCorrection()->FillMeas(vtx, eta, pt);
+ fVertexRecoCorrection->GetTrackCorrection()->FillGene(vtx, eta, pt);
+
+ if (!vertex)
+ return;
+
+ fVertexRecoCorrection->GetTrackCorrection()->FillMeas(vtx, eta, pt);
+ fTrack2ParticleCorrection->GetTrackCorrection()->FillGene(vtx, eta, pt);
+}
+
+//____________________________________________________________________
+void AlidNdEtaCorrection::FillTrackedParticle(Float_t vtx, Float_t eta, Float_t pt)
+{
+ // fills a tracked particle in the corrections
+
+ fTrack2ParticleCorrection->GetTrackCorrection()->FillMeas(vtx, eta, pt);
+}
+
+//____________________________________________________________________
+void AlidNdEtaCorrection::FillEvent(Float_t vtx, Float_t n, Bool_t trigger, Bool_t vertex, Int_t processType)
+{
+ // fills an event int he correction
+ // it is filled in generated or measured depending of the flags
+
+ fTriggerBiasCorrectionMBToINEL->GetEventCorrection()->FillGene(vtx, n);
+
+ if (processType != 92 && processType != 93)
+ fTriggerBiasCorrectionMBToNSD->GetEventCorrection()->FillGene(vtx, n);
+
+ if (processType!=92 && processType!=93 && processType!=94)
+ fTriggerBiasCorrectionMBToND->GetEventCorrection()->FillGene(vtx, n);
+
+ if (!trigger)
+ return;
+
+ fTriggerBiasCorrectionMBToINEL->GetEventCorrection()->FillMeas(vtx, n);
+ fTriggerBiasCorrectionMBToNSD->GetEventCorrection()->FillMeas(vtx, n);
+ fTriggerBiasCorrectionMBToND->GetEventCorrection()->FillMeas(vtx, n);
+ fVertexRecoCorrection->GetEventCorrection()->FillGene(vtx, n);
+
+ if (!vertex)
+ return;
+
+ fVertexRecoCorrection->GetEventCorrection()->FillMeas(vtx, n);
+}
+
+//____________________________________________________________________
+Float_t AlidNdEtaCorrection::GetMeasuredFraction(CorrectionType correctionType, Float_t ptCutOff, Float_t eta, Bool_t debug)
+{
+ // calculates the fraction of particles measured (some are missed due to the pt cut off)
+ //
+ // uses the generated particle histogram from the correction passed, e.g. pass GetTrack2ParticleCorrection()
+
+ if (!GetCorrection(correctionType))
+ return -1;
+
+ const TH3F* generated = GetCorrection(correctionType)->GetTrackCorrection()->GetGeneratedHistogram();
+
+ // find eta borders, if eta is negative assume -0.8 ... 0.8
+ Int_t etaBegin = 0;
+ Int_t etaEnd = 0;
+ if (eta < -99)
+ {
+ etaBegin = generated->GetYaxis()->FindBin(-0.8);
+ etaEnd = generated->GetYaxis()->FindBin(0.8);
+ }
+ else
+ {
+ etaBegin = generated->GetYaxis()->FindBin(eta);
+ etaEnd = etaBegin;
+ }
+
+ //Printf("AlidNdEtaCorrection::GetMeasuredFraction: Using vtx range of +- 10 cm");
+ Int_t vertexBegin = generated->GetXaxis()->FindBin(-9.99);
+ Int_t vertexEnd = generated->GetXaxis()->FindBin(9.99);
+
+ TH1D* ptProj = dynamic_cast<TH1D*> (generated->ProjectionZ(Form("%s_pt", generated->GetName()), vertexBegin, vertexEnd, etaBegin, etaEnd));
+ //printf("GetMeasuredFraction: bin range %d %d %d %d\n", vertexBegin, vertexEnd, etaBegin, etaEnd);
+ ptProj->GetXaxis()->SetTitle(generated->GetZaxis()->GetTitle());
+
+ Int_t ptBin = ptProj->FindBin(ptCutOff);
+ //printf("GetMeasuredFraction: bin range %d %d\n", ptBin, ptProj->GetNbinsX());
+ Float_t abovePtCut = ptProj->Integral(ptBin, ptProj->GetNbinsX()+1);
+ Float_t all = ptProj->Integral(1, ptProj->GetNbinsX()+1);
+
+ if (all == 0)
+ return -1;
+
+ Float_t fraction = abovePtCut / all;
+
+ //printf("GetMeasuredFraction: all %f above %f fraction %f\n", all, abovePtCut, fraction);
+
+ if (debug)
+ {
+ new TCanvas;
+ ptProj->Draw();
+ }
+ else
+ delete ptProj;
+
+ if (debug)
+ printf("AlidNdEtaCorrection::GetMeasuredFraction: pt cut off = %f, eta = %f, => fraction = %f\n", ptCutOff, eta, fraction);
+
+ return fraction;
+}
+
+//____________________________________________________________________
+TH1* AlidNdEtaCorrection::GetMeasuredEventFraction(CorrectionType correctionType, Int_t multCut)
+{
+ // calculates the fraction of events above multCut (but including it)
+ //
+ // uses the generated event histogram from the correction passed, e.g. pass GetTrack2ParticleCorrection()
+
+ if (!GetCorrection(correctionType))
+ return 0;
+
+ const TH2F* generated = GetCorrection(correctionType)->GetEventCorrection()->GetGeneratedHistogram();
+
+ TH1* allEvents = generated->ProjectionX(Form("%s_all", generated->GetName()), 1, generated->GetNbinsY());
+ TH1* aboveEvents = generated->ProjectionX(Form("%s_above", generated->GetName()), generated->GetYaxis()->FindBin(multCut), generated->GetNbinsY());
+
+ aboveEvents->Divide(aboveEvents, allEvents, 1, 1, "B");
+
+ return aboveEvents;
+}
+
+//____________________________________________________________________
+void AlidNdEtaCorrection::ReduceInformation()
+{
+ // this function deletes the measured and generated histograms from the corrections to reduce the amount of data
+ // in memory
+
+ // these are needed for GetMeasuredFraction(): fTrack2ParticleCorrection->ReduceInformation();
+ fVertexRecoCorrection ->ReduceInformation();
+ fTriggerBiasCorrectionMBToINEL ->ReduceInformation();
+ fTriggerBiasCorrectionMBToNSD ->ReduceInformation();
+ fTriggerBiasCorrectionMBToND ->ReduceInformation();
+}
+
+//____________________________________________________________________
+AliCorrection* AlidNdEtaCorrection::GetCorrection(CorrectionType correctionType)
+{
+ // returns the given correction
+
+ switch (correctionType)
+ {
+ case kNone : return 0;
+ case kTrack2Particle : return fTrack2ParticleCorrection;
+ case kVertexReco : return fVertexRecoCorrection;
+ case kINEL : return fTriggerBiasCorrectionMBToINEL;
+ case kNSD : return fTriggerBiasCorrectionMBToNSD;
+ case kND : return fTriggerBiasCorrectionMBToND;
+ }
+ return 0;
}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff