From 15b086e2e6aaf9c6d71fafb84881112315af7769 Mon Sep 17 00:00:00 2001 From: morsch Date: Thu, 6 Jan 2011 16:34:35 +0000 Subject: [PATCH] Additional histograms Marek Chojnacki --- PWG1/ITS/AliAnalysisTaskdEdxSSDQA.cxx | 13 +++++++++---- PWG1/ITS/AliAnalysisTaskdEdxSSDQA.h | 7 ++++--- 2 files changed, 13 insertions(+), 7 deletions(-) diff --git a/PWG1/ITS/AliAnalysisTaskdEdxSSDQA.cxx b/PWG1/ITS/AliAnalysisTaskdEdxSSDQA.cxx index bcd5f49e194..62a5eab4ace 100644 --- a/PWG1/ITS/AliAnalysisTaskdEdxSSDQA.cxx +++ b/PWG1/ITS/AliAnalysisTaskdEdxSSDQA.cxx @@ -28,7 +28,7 @@ ClassImp(AliAnalysisTaskdEdxSSDQA) //________________________________________________________________________ AliAnalysisTaskdEdxSSDQA::AliAnalysisTaskdEdxSSDQA(const char *name) -: AliAnalysisTaskSE(name), fHist1(0),fHist2(0),fListOfHistos(0),fPcut(0.0) +: AliAnalysisTaskSE(name), fHist1(0),fHist2(0),fHist3(0),fListOfHistos(0),fPcut(0.0) { // Constructor // Define input and output slots here @@ -47,6 +47,8 @@ void AliAnalysisTaskdEdxSSDQA::UserCreateOutputObjects() fListOfHistos->Add(fHist1); fHist2=new TH2F("QACharge","QACharge;Module;Q",1698,-0.5,1697.5,150,0,300); fListOfHistos->Add(fHist2); + fHist3=new TH2F("ChargeRatiovCharge","CRvQ;Q;CR",200,0,2000,100,-1.0,1.0); + fListOfHistos->Add(fHist3); } //________________________________________________________________________ void AliAnalysisTaskdEdxSSDQA::LocalInit() @@ -85,7 +87,8 @@ void AliAnalysisTaskdEdxSSDQA::UserExec(Option_t *) AliTrackPoint point; Int_t nlayer=0; Int_t id=0; - Float_t chargeratio=0; + Float_t chargeratio=0.0; + Float_t charge=0.0; for(int itr=0;itrGetTrack(itr); @@ -116,7 +119,8 @@ void AliAnalysisTaskdEdxSSDQA::UserExec(Option_t *) { if(point.GetCharge()>0.0&&point.IsExtra()==kFALSE)//do not use additional clusters { - chargeratio=point.GetChargeRatio(); + chargeratio=point.GetChargeRatio(); + charge=point.GetCharge(); if(nlayer==5&&tmpQESD[2]>0.0) { fHist1->Fill(id,chargeratio); @@ -128,7 +132,8 @@ void AliAnalysisTaskdEdxSSDQA::UserExec(Option_t *) fHist1->Fill(id+748,chargeratio); if(track->GetP()>fPcut) fHist2->Fill(id+748,tmpQESD[3]); - } + } + fHist3->Fill(charge,chargeratio); } } } diff --git a/PWG1/ITS/AliAnalysisTaskdEdxSSDQA.h b/PWG1/ITS/AliAnalysisTaskdEdxSSDQA.h index 404d088046a..d8ca925a3de 100644 --- a/PWG1/ITS/AliAnalysisTaskdEdxSSDQA.h +++ b/PWG1/ITS/AliAnalysisTaskdEdxSSDQA.h @@ -22,8 +22,9 @@ class AliAnalysisTaskdEdxSSDQA : public AliAnalysisTaskSE { Float_t GetPcut() const{return fPcut;} private: - TH2F* fHist1; // signal differenc between ITSsignal stored in ESD and calcuated after using AliTrackPoint - TH2F* fHist2; // landau distributions for each module + TH2F* fHist1; // CR for each module + TH2F* fHist2; // landau distributions for each module + TH2F* fHist3; // CR as function of Charge for the AliTrackPoint TList* fListOfHistos; // output list Float_t fPcut; // Momentum cut @@ -31,7 +32,7 @@ class AliAnalysisTaskdEdxSSDQA : public AliAnalysisTaskSE { AliAnalysisTaskdEdxSSDQA(const AliAnalysisTaskdEdxSSDQA&); // not implemented AliAnalysisTaskdEdxSSDQA& operator=(const AliAnalysisTaskdEdxSSDQA&); // not implemented - ClassDef(AliAnalysisTaskdEdxSSDQA, 1); // example of analysis + ClassDef(AliAnalysisTaskdEdxSSDQA, 2); // example of analysis }; #endif -- 2.43.0