X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PWG1%2FTRD%2FAliTRDcheckDET.cxx;fp=PWG1%2FTRD%2FAliTRDcheckDET.cxx;h=d1a2108d5805c6673ccb5cfadcab0f92c0e06eb1;hb=40600e64b92ae52e4a65d50426d7af033e26c920;hp=bfad5a5d712090287e778a61c93614dea18f7183;hpb=277c5bf4c528ef45220585cec52ca36dec2e3f50;p=u%2Fmrichter%2FAliRoot.git diff --git a/PWG1/TRD/AliTRDcheckDET.cxx b/PWG1/TRD/AliTRDcheckDET.cxx index bfad5a5d712..d1a2108d580 100644 --- a/PWG1/TRD/AliTRDcheckDET.cxx +++ b/PWG1/TRD/AliTRDcheckDET.cxx @@ -498,11 +498,11 @@ TObjArray *AliTRDcheckDET::Histos(){ arr->SetOwner(kTRUE); arr->SetName(""); fContainer->AddAt(arr, kPH); if(!(h3 = (TH3F *)gROOT->FindObject("hPHx"))){ - h3 = new TH3F("hPHx", "(x);x[mm];Centrality;Charge[a.u.]", 31, -0.08, 4.88, 100, 0, 1024, AliTRDeventInfo::kCentralityClasses + 1, -1.5, AliTRDeventInfo::kCentralityClasses - 0.5); + h3 = new TH3F("hPHx", "(x);x[mm];Charge[a.u.];Centrality", 31, -0.08, 4.88, 100, 0, 1024, AliTRDeventInfo::kCentralityClasses + 1, -1.5, AliTRDeventInfo::kCentralityClasses - 0.5); } else h3->Reset(); arr->AddAt(h3, 0); if(!(h3 = (TH3F *)gROOT->FindObject("hPHt"))){ - h3 = new TH3F("hPHt", "(t);time[100ns];Centrality;Charge[a.u.]", 31, -0.5, 30.5, 100, 0, 1024, AliTRDeventInfo::kCentralityClasses + 1, -1.5, AliTRDeventInfo::kCentralityClasses - 0.5); + h3 = new TH3F("hPHt", "(t);time[100ns];Charge[a.u.];Centrality", 31, -0.5, 30.5, 100, 0, 1024, AliTRDeventInfo::kCentralityClasses + 1, -1.5, AliTRDeventInfo::kCentralityClasses - 0.5); } else h3->Reset(); arr->AddAt(h3, 1); @@ -1703,6 +1703,9 @@ Bool_t AliTRDcheckDET::MakePlotPulseHeight(){ TH1 *hPhtCent[AliTRDeventInfo::kCentralityClasses]; TH2 *hPtmp; Int_t nHistsCentrality = 0; + TArrayD ptbinning(hPht->GetNbinsX()+1); + for(Int_t ibin = 1; ibin <=hPht->GetNbinsX(); ibin++) ptbinning[ibin-1] = hPht->GetXaxis()->GetBinLowEdge(ibin); + ptbinning[hPht->GetNbinsX()] = hPht->GetXaxis()->GetBinUpEdge(hPht->GetNbinsX()); for(Int_t icent = 0; icent < AliTRDeventInfo::kCentralityClasses; icent++){ hPht->GetZaxis()->SetRange(icent+2,icent+2); hPtmp = dynamic_cast(hPht->Project3D("yx")); @@ -1710,7 +1713,7 @@ Bool_t AliTRDcheckDET::MakePlotPulseHeight(){ hPhtCent[icent] = NULL; continue; } - hPhtCent[icent] = new TH1F(Form("hPhtCent_%d", icent), "Average Ph vs Time", hPtmp->GetNbinsX(), hPtmp->GetXaxis()->GetXbins()->GetArray()); + hPhtCent[icent] = new TH1F(Form("hPhtCent_%d", icent), "Average Ph vs Time", hPtmp->GetNbinsX(), ptbinning.GetArray()); for(Int_t it = 1; it <= hPtmp->GetNbinsX(); it++){ htmp = hPtmp->ProjectionY("htmp", it, it); hPhtCent[icent]->SetBinContent(it, htmp->GetMean());