// --- ROOT system ---
#include <TClonesArray.h>
#include <TFile.h>
-#include <TH1D.h>
-#include <TH2D.h>
-#include <TH3D.h>
+#include <TH1F.h>
+#include <TH2F.h>
+#include <TH3F.h>
#include <TProfile.h>
#include <TF1.h>
#include <TCanvas.h>
#include "AliTRDcluster.h"
#include "AliTRDQADataMakerRec.h"
#include "AliTRDgeometry.h"
-//#include "AliTRDdataArrayI.h"
-#include "AliTRDrawStreamOld.h"
+#include "AliTRDrawStream.h"
#include "AliTRDdigitsManager.h"
#include "AliTRDSignalIndex.h"
//
//TStopwatch watch;
//watch.Start();
- /**/
+ ResetEventTrigClasses();
AliDebug(AliQAv1::GetQADebugLevel(), "End of TRD cycle");
-
-
- if (task == AliQAv1::kRAWS) {
-
- // loop over event types
- for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
-
- if (! IsValidEventSpecie(specie, list)) continue;
-
-
- TH2D *mnCls = (TH2D*)list[specie]->At(1);
- TH2D *mClsDet = (TH2D*)list[specie]->At(2);
-
- mClsDet->Divide(mnCls);
- }
- }
-
-
- if (task == AliQAv1::kRECPOINTS) {
-
- TH1D * hist = new TH1D("fitHist", "", 200, -0.5, 199.5);
-
- // loop over event types
- for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
- if (! IsValidEventSpecie(specie, list))
- continue ;
- //list[specie]->Print();
-
- // fill detector map;
- for(Int_t i = 0 ; i < 540 ; i++) {
- Double_t v = ((TH1D*)list[specie]->At(0))->GetBinContent(i+1);
- Int_t sm = i/30;
- Int_t det = i%30;
- TH2D *detMap = (TH2D*)list[specie]->At(87);
- Int_t bin = detMap->FindBin(sm, det);
- detMap->SetBinContent(bin, v);
- }
-
- // Rec points full chambers
- for (Int_t i = 0 ; i < 540 ; i++) {
- //AliDebug(AliQAv1::GetQADebugLevel(), Form("I = %d", i));
- //TH1D *h = ((TH2D*)list[specie]->At(1))->ProjectionY(Form("qaTRD_recPoints_amp_%d",i), i+1, i+1);
- hist->Reset();
-
- // project TH2D into TH1D
- for(Int_t b = 1 ; b < hist->GetXaxis()->GetNbins()-1 ; b++) {
- Double_t xvalue = hist->GetBinCenter(b);
- Int_t bin = ((TH2D*)list[specie]->At(1))->FindBin(i,xvalue);
- Double_t value = ((TH2D*)list[specie]->At(1))->GetBinContent(bin);
- hist->SetBinContent(b, value);
- }
- //AliDebug(AliQAv1::GetQADebugLevel(), Form("Sum = %d %f\n", i, hist->GetSum()));
- if (hist->GetSum() < 100)
- continue; // not enougth data in a chamber
-
- hist->Fit("landau", "q0", "goff", 10, 180);
- TF1 *fit = hist->GetFunction("landau");
- ((TH1D*)list[specie]->At(12))->Fill(fit->GetParameter(1));
- ((TH1D*)list[specie]->At(13))->Fill(fit->GetParameter(2));
+ //
+ for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) { // loop over event types
+ //
+ if (!IsValidEventSpecie(specie, list)) continue;
+ SetEventSpecie(AliRecoParam::ConvertIndex(specie));
+ //
+ for (int itc=-1;itc<GetNTrigClasses();itc++) { // RS: loop over eventual clones per trigger class
+ //
+ if (task == AliQAv1::kRAWS) {
+ TH2 *mnCls = (TH2*)GetRawsData(1, itc);
+ TH2 *mClsDet = (TH2*)GetRawsData(2, itc);
+ if (mnCls && mClsDet) mClsDet->Divide(mnCls);
}
-
-
- // time-bin by time-bin sm by sm
- for(Int_t i=0; i<18; i++) { // loop over super-modules
- for(Int_t j=0; j<kTimeBin; j++) { // loop over time bins
-
- hist->Reset();
- for(Int_t b = 1 ; b < hist->GetXaxis()->GetNbins()-1 ; b++) {
- Double_t xvalue = hist->GetBinCenter(b);
- Double_t svalue = 0.0;
- for(Int_t det = i*30 ; det < (i+1)*30 ; det++) { // loop over detectors
- Int_t bin = ((TH3D*)list[specie]->At(10))->FindBin(det,j,xvalue);
- Double_t value = ((TH3D*)list[specie]->At(10))->GetBinContent(bin);
- svalue += value;
- }
- //AliDebug(AliQAv1::GetQADebugLevel(), Form("v = %f\n", value));
- hist->SetBinContent(b, svalue);
- }
-
- if (hist->GetSum() < 100)
- continue;
-
- hist->Fit("landau", "q0", "goff", 10, 180);
- TF1 *fit = hist->GetFunction("landau");
-
- TH1D *h1 = (TH1D*)list[specie]->At(14+18+i);
- h1->SetMarkerStyle(20);
- Int_t bin = h1->FindBin(j);
- // printf("%d %d %d\n", det, j, bin);
-
- Double_t value = TMath::Abs(fit->GetParameter(1));
- Double_t error = TMath::Abs(fit->GetParError(1));
-
- if (value/error < 3) continue; // insuficient statistics
-
- h1->SetBinContent(bin, value);
- h1->SetBinError(bin, error);
+ //
+ if (task == AliQAv1::kRECPOINTS) {
+ //
+ TH1F * hist = new TH1F("fitHist", "", 200, -0.5, 199.5);
+ TObjArray& arrRP = *GetRecPointsDataOfTrigClass(itc); // RS Histos matching to trigger class
+ // fill detector map;
+ TH1* h0 = (TH1*) arrRP[0];
+ TH1* detMap = (TH1*) arrRP[87];
+ if (h0 && detMap) {
+ for(Int_t i = 0 ; i < 540 ; i++) {
+ Double_t v = h0->GetBinContent(i+1);
+ Int_t sm = i/30;
+ Int_t det = i%30;
+ Int_t bin = detMap->FindBin(sm, det);
+ detMap->SetBinContent(bin, v);
+ }
}
- }
-
- // for numerical convergence
- TF1 *form = new TF1("formLandau", "landau", 0, 200);
-
- // time-bin by time-bin chamber by chamber
- for (Int_t i=0; i<540; i++) {
- for(Int_t j=0; j<kTimeBin; j++) {
-
- hist->Reset();
- for(Int_t b = 1 ; b < hist->GetXaxis()->GetNbins()-1 ; b++) {
- Double_t xvalue = hist->GetBinCenter(b);
- Int_t bin = ((TH3D*)list[specie]->At(10))->FindBin(i,j,xvalue);
- Double_t value = ((TH3D*)list[specie]->At(10))->GetBinContent(bin);
- //AliDebug(AliQAv1::GetQADebugLevel(), Form("v = %f\n", value));
- hist->SetBinContent(b, value);
+ //
+ // Rec points full chambers
+ TH2* h2tmp = (TH2*) arrRP[1];
+ TH1* h12 = (TH1*) arrRP[12];
+ TH1* h13 = (TH1*) arrRP[13];
+ if (h2tmp && h12 && h13) {
+ for (Int_t i = 0 ; i < 540 ; i++) {
+ //AliDebug(AliQAv1::GetQADebugLevel(), Form("I = %d", i));
+ hist->Reset();
+ // project TH2F into TH1F
+ for(Int_t b = 1 ; b < hist->GetXaxis()->GetNbins()-1 ; b++) {
+ Double_t xvalue = hist->GetBinCenter(b);
+ Int_t bin = h2tmp->FindBin(i,xvalue);
+ Double_t value = h2tmp->GetBinContent(bin);
+ hist->SetBinContent(b, value);
+ }
+ //AliDebug(AliQAv1::GetQADebugLevel(), Form("Sum = %d %f\n", i, hist->GetSum()));
+ if (hist->GetSum() < 100) continue; // not enougth data in a chamber
+ //
+ hist->Fit("landau", "q0", "goff", 10, 180);
+ TF1 *fit = hist->GetFunction("landau");
+ h12->Fill(fit->GetParameter(1));
+ h13->Fill(fit->GetParameter(2));
}
-
- if (hist->GetSum() < 100)
- continue;
-
- form->SetParameters(1000, 60, 20);
- hist->Fit(form, "q0", "goff", 20, 180);
-
- Int_t sm = i/30;
- Int_t det = i%30;
- TH2D *h2 = (TH2D*)list[specie]->At(14+sm);
- Int_t bin = h2->FindBin(det,j);
- // printf("%d %d %d\n", det, j, bin);
-
- Double_t value = TMath::Abs(form->GetParameter(1));
- Double_t error = TMath::Abs(form->GetParError(1));
-
- if (value/error < 3) continue;
-
- h2->SetBinContent(bin, value);
- h2->SetBinError(bin, error);
}
- }
- }
- if (hist)
- delete hist;
- }
- //////////////////////////
- // const Int_t knbits = 6;
- // const char *suf[knbits] = {"TPCi", "TPCo", "TPCz", "TRDo", "TRDr", "TRDz"};
- //const char *sufRatio[4] = {"TRDrTRDo", "TRDoTPCo", "TRDrTPCo", "TRDzTPCo"};
-
- if (task == AliQAv1::kESDS) {
-
- const Int_t knRatio = 4;
- const Int_t kN[knRatio] = {4,3,4,5};
- const Int_t kD[knRatio] = {3,1,1,3};
-
- // create ratios
- for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
- if (! IsValidEventSpecie(specie, list))
- continue ;
- for(Int_t type = 0 ; type < 2 ; type++) {
- for(Int_t i = 0 ; i < knRatio ; i++) {
- TH1D *ratio = (TH1D*)list[specie]->At(19 + 2*i + type);
- TH1D *histN = (TH1D*)list[specie]->At(3 + 2*kN[i] + type);
- TH1D *histD = (TH1D*)list[specie]->At(3 + 2*kD[i] + type);
- BuildRatio(ratio, histN, histD);
- //ratio->Reset();
- //ratio->Add(histN);
- //ratio->Divide(histD);
- }
- }
- // ratio for the fraction of electrons per stack
- TH1D *histN = (TH1D*)list[specie]->At(33);
- TH1D *histD = (TH1D*)list[specie]->At(32);
- TH1D *ratio = (TH1D*)list[specie]->At(34);
- BuildRatio(ratio, histN, histD);
- }
- }
+ //
+ // time-bin by time-bin sm by sm
+ TH3* h3tmp = (TH3*) arrRP[10];
+ if (h3tmp) {
+ for(Int_t i=0; i<18; i++) { // loop over super-modules
+ for(Int_t j=0; j<kTimeBin; j++) { // loop over time bins
+ hist->Reset();
+ for(Int_t b = 1 ; b < hist->GetXaxis()->GetNbins()-1 ; b++) {
+ Double_t xvalue = hist->GetBinCenter(b);
+ Double_t svalue = 0.0;
+ for(Int_t det = i*30 ; det < (i+1)*30 ; det++) { // loop over detectors
+ Int_t bin = h3tmp->FindBin(det,j,xvalue);
+ Double_t value = h3tmp->GetBinContent(bin);
+ svalue += value;
+ }
+ //AliDebug(AliQAv1::GetQADebugLevel(), Form("v = %f\n", value));
+ hist->SetBinContent(b, svalue);
+ }
+ //
+ if (hist->GetSum() < 100) continue;
+ //
+ hist->Fit("landau", "q0", "goff", 10, 180);
+ TF1 *fit = hist->GetFunction("landau");
+ //
+ TH1 *hi = (TH1*)arrRP[14+18+i];
+ if (!hi) continue;
+ hi->SetMarkerStyle(20);
+ Int_t bin = hi->FindBin(j);
+ // printf("%d %d %d\n", det, j, bin);
+ //
+ Double_t value = TMath::Abs(fit->GetParameter(1));
+ Double_t error = TMath::Abs(fit->GetParError(1));
+ //
+ if (value/error < 3) continue; // insuficient statistics
+ //
+ hi->SetBinContent(bin, value);
+ hi->SetBinError(bin, error);
+ } // j
+ } // i
+ //
+ // for numerical convergence
+ TF1 *form = new TF1("formLandau", "landau", 0, 200);
+ //
+ // time-bin by time-bin chamber by chamber
+ for (Int_t i=0; i<540; i++) {
+ for(Int_t j=0; j<kTimeBin; j++) {
+ //
+ hist->Reset();
+ for(Int_t b = 1 ; b < hist->GetXaxis()->GetNbins()-1 ; b++) {
+ Double_t xvalue = hist->GetBinCenter(b);
+ Int_t bin = h3tmp->FindBin(i,j,xvalue);
+ Double_t value = h3tmp->GetBinContent(bin);
+ //AliDebug(AliQAv1::GetQADebugLevel(), Form("v = %f\n", value));
+ hist->SetBinContent(b, value);
+ }
+ //
+ if (hist->GetSum() < 100) continue;
+
+ form->SetParameters(1000, 60, 20);
+ hist->Fit(form, "q0", "goff", 20, 180);
+
+ Int_t sm = i/30;
+ Int_t det = i%30;
+ TH1 *hi = (TH1*)arrRP[14+sm];
+ if (!hi) continue;
+ Int_t bin = hi->FindBin(det,j);
+ // printf("%d %d %d\n", det, j, bin);
+ //
+ Double_t value = TMath::Abs(form->GetParameter(1));
+ Double_t error = TMath::Abs(form->GetParError(1));
+ //
+ if (value/error < 3) continue;
+ //
+ hi->SetBinContent(bin, value);
+ hi->SetBinError(bin, error);
+ } // j
+ } // i
+ } // h3tmp
+ if (hist) delete hist;
+ } // RECPOINTS
+ //
+ //////////////////////////
+ // const Int_t knbits = 6;
+ // const char *suf[knbits] = {"TPCi", "TPCo", "TPCz", "TRDo", "TRDr", "TRDz"};
+ //const char *sufRatio[4] = {"TRDrTRDo", "TRDoTPCo", "TRDrTPCo", "TRDzTPCo"};
+
+ if (task == AliQAv1::kESDS) {
+ //
+ const Int_t knRatio = 4;
+ const Int_t kN[knRatio] = {4,3,4,5};
+ const Int_t kD[knRatio] = {3,1,1,3};
+ //
+ TObjArray& arrES = *GetESDsDataOfTrigClass(itc); // RS Histos matching to trigger class
+ // create ratios
+ for(Int_t type = 0 ; type < 2 ; type++) {
+ for(Int_t i = 0 ; i < knRatio ; i++) {
+ TH1 *ratio = (TH1*)arrES[19 + 2*i + type];
+ TH1 *histN = (TH1*)arrES[3 + 2*kN[i] + type];
+ TH1 *histD = (TH1*)arrES[3 + 2*kD[i] + type];
+ if ( ! (ratio && histN && histD) ) continue;
+ BuildRatio(ratio, histN, histD);
+ //ratio->Reset();
+ //ratio->Add(histN);
+ //ratio->Divide(histD);
+ }
+ }
+ // ratio for the fraction of electrons per stack
+ TH1 *histN = (TH1*)arrES[33];
+ TH1 *histD = (TH1*)arrES[32];
+ TH1 *ratio = (TH1*)arrES[34];
+ BuildRatio(ratio, histN, histD);
+ } // ESDs
+ } // RS: loop over eventual clones per trigger class
+ } // loop over species
// call the checker
AliQAChecker::Instance()->Run(AliQAv1::kTRD, task, list) ;
-
+ //
}
//____________________________________________________________________________
TH1 *hist[kNhist];
Int_t histoCounter = -1 ;
- hist[++histoCounter] = new TH1D("qaTRD_esd_ntracks", "TRD esd ntracks;Number of tracks;Counts", 300, -0.5, 299.5);
- hist[++histoCounter] = new TH1D("qaTRD_esd_sector", "TRD esd sector;Sector;Counts", 18, -0.5, 17.7);
- hist[++histoCounter] = new TH1D("qaTRD_esd_bits", "TRD esd bits;Bits;Counts", 64, -0.5, 63.5);
+ hist[++histoCounter] = new TH1F("qaTRD_esd_ntracks", "TRD esd ntracks;Number of tracks;Counts", 300, -0.5, 299.5);
+ hist[++histoCounter] = new TH1F("qaTRD_esd_sector", "TRD esd sector;Sector;Counts", 18, -0.5, 17.7);
+ hist[++histoCounter] = new TH1F("qaTRD_esd_bits", "TRD esd bits;Bits;Counts", 64, -0.5, 63.5);
const Int_t knbits = 6;
const char *suf[knbits] = {"TPCi", "TPCo", "TPCz", "TRDo", "TRDr", "TRDz"};
// histo = 3
for(Int_t i=0; i<knbits; i++) {
- hist[++histoCounter] = new TH1D(Form("qaTRD_esd_pt%s",suf[i]), Form("qaTRD_esd_pt%s;p_{T} (GeV/c);Counts",suf[i]), 100, 0, 10);
- hist[++histoCounter] = new TH1D(Form("qaTRD_esd_trdz%s", suf[i]), ";z (cm)", 200, -400, 400);
+ hist[++histoCounter] = new TH1F(Form("qaTRD_esd_pt%s",suf[i]), Form("qaTRD_esd_pt%s;p_{T} (GeV/c);Counts",suf[i]), 100, 0, 10);
+ hist[++histoCounter] = new TH1F(Form("qaTRD_esd_trdz%s", suf[i]), ";z (cm)", 200, -400, 400);
}
- hist[++histoCounter] = new TH1D("qaTRD_esd_clsTRDo", "TRDo;number of clusters;Counts", 180, -0.5, 179.5);;
- hist[++histoCounter] = new TH1D("qaTRD_esd_clsTRDr", "TRDr;number of clusters;Counts", 180, -0.5, 179.5);;
- hist[++histoCounter] = new TH1D("qaTRD_esd_clsTRDz", "TRDz;number of clusters;Counts", 180, -0.5, 179.5);;
- //hist[++histoCounter] = new TH1D("qaTRD_esd_clsRatio", ";cluster ratio", 100, 0., 1.3);;
+ hist[++histoCounter] = new TH1F("qaTRD_esd_clsTRDo", "TRDo;number of clusters;Counts", 180, -0.5, 179.5);;
+ hist[++histoCounter] = new TH1F("qaTRD_esd_clsTRDr", "TRDr;number of clusters;Counts", 180, -0.5, 179.5);;
+ hist[++histoCounter] = new TH1F("qaTRD_esd_clsTRDz", "TRDz;number of clusters;Counts", 180, -0.5, 179.5);;
+ //hist[++histoCounter] = new TH1F("qaTRD_esd_clsRatio", ";cluster ratio", 100, 0., 1.3);;
- hist[++histoCounter] = new TH2D("qaTRD_esd_sigMom", "TRD esd sig Mom;momentum (GeV/c);signal", 100, 0, 5, 200, 0, 1e3);
+ hist[++histoCounter] = new TH2F("qaTRD_esd_sigMom", "TRD esd sig Mom;momentum (GeV/c);signal", 100, 0, 5, 200, 0, 1e3);
// end of cycle plots (hist 19)
const char *sufRatio[4] = {"TRDrTRDo", "TRDoTPCo", "TRDrTPCo", "TRDzTPCo"};
for(int i=0; i<4; i++) {
- hist[++histoCounter] = new TH1D(Form("qaTRD_esd_pt%s",sufRatio[i]),
+ hist[++histoCounter] = new TH1F(Form("qaTRD_esd_pt%s",sufRatio[i]),
Form("Efficiency in Pt %s;p_{T};eff", sufRatio[i]),
100, 0, 10);
- hist[++histoCounter] = new TH1D(Form("qaTRD_esd_trdz%s",sufRatio[i]),
+ hist[++histoCounter] = new TH1F(Form("qaTRD_esd_trdz%s",sufRatio[i]),
Form("Efficiency in Z %s;z (cm);eff", sufRatio[i]),
200, -400, 400);
}
// 27 - 31
- hist[27] = new TH1D("qaTRD_esd_quality", "TRD esd quality;quality;Counts", 120, 0, 12);
- hist[28] = new TH1D("qaTRD_esd_budget", "TRD esd budget;NN;Counts", 110, -1000, 100);
- hist[29] = new TH1D("qaTRD_esd_chi2", "TRD esd chi2;chi2;Counts", 200, 0, 100);
- hist[30] = new TH1D("qaTRD_esd_timeBin", "TRD esd timeBin;time bin;Counts", 7, -0.5, 6.5);
- hist[31] = new TH1D("qaTRD_esd_pidQuality", "pid Quality;quality;;Counts", 7, -0.5, 6.5);
+ hist[27] = new TH1F("qaTRD_esd_quality", "TRD esd quality;quality;Counts", 120, 0, 12);
+ hist[28] = new TH1F("qaTRD_esd_budget", "TRD esd budget;NN;Counts", 110, -1000, 100);
+ hist[29] = new TH1F("qaTRD_esd_chi2", "TRD esd chi2;chi2;Counts", 200, 0, 100);
+ hist[30] = new TH1F("qaTRD_esd_timeBin", "TRD esd timeBin;time bin;Counts", 7, -0.5, 6.5);
+ hist[31] = new TH1F("qaTRD_esd_pidQuality", "pid Quality;quality;;Counts", 7, -0.5, 6.5);
// stack by stack electron identyfication
- hist[32] = new TH1D("qaTRD_esd_tracksStack", "number of all tracks;stack;Counts", 90, -0.5, 89.5);
- hist[33] = new TH1D("qaTRD_esd_electronStack", "number of electron tracks;stack;Counts", 90, -0.5, 89.5);
- hist[34] = new TH1D("qaTRD_esd_elRatioStack", "fraction of electron tracks;stack;Counts", 90, -0.5, 89.5);
- hist[35] = new TH1D("qaTRD_esd_thetaOut", ";tan(theta);", 100, -1, 1);
+ hist[32] = new TH1F("qaTRD_esd_tracksStack", "number of all tracks;stack;Counts", 90, -0.5, 89.5);
+ hist[33] = new TH1F("qaTRD_esd_electronStack", "number of electron tracks;stack;Counts", 90, -0.5, 89.5);
+ hist[34] = new TH1F("qaTRD_esd_elRatioStack", "fraction of electron tracks;stack;Counts", 90, -0.5, 89.5);
+ hist[35] = new TH1F("qaTRD_esd_thetaOut", ";tan(theta);", 100, -1, 1);
const char *partType[5] = {"Electron", "Muon", "Pion", "Kaon", "Proton"};
for(Int_t i=0; i<AliPID::kSPECIES; i++)
- hist[36+i] = new TH1D(Form("qaTRD_esd_pid%d",i),
+ hist[36+i] = new TH1F(Form("qaTRD_esd_pid%d",i),
Form("%s;probability;Counts",partType[i]), 100, 0, 1);
// dE/dX vs momentum in three regions
for(Int_t i=0; i<4; i++) {
- hist[41+i] = new TH2D(Form("qaTRD_esd_signalPzone_%d",i),
+ hist[41+i] = new TH2F(Form("qaTRD_esd_signalPzone_%d",i),
Form("%s;momentum (GeV/c);signal (a.u.)", zoneName[i]),
nscalex, scalex, nscaley, scaley);
}
//hist[i]->Sumw2();
Add2ESDsList(hist[i], i, !expert, image);
}
-
+ //
+ ClonePerTrigClass(AliQAv1::kESDS); // this should be the last line
}
//____________________________________________________________________________
const Int_t kNhist = 14 + 4 * 18 + 2 + 9;// + 540;
TH1 *hist[kNhist];
- hist[0] = new TH1D("qaTRD_recPoints_det", "RRD recPoints det;Detector ID of the cluster;Counts", 540, -0.5, 539.5);
- hist[1] = new TH2D("qaTRD_recPoints_amp", "TRD recPoints amp;Amplitude;??", 540, -0.5, 539, 200, -0.5, 199.5);
- hist[2] = new TH1D("qaTRD_recPoints_npad", "TRD recPoints npad;Number of Pads;Counts", 12, -0.5, 11.5);
+ hist[0] = new TH1F("qaTRD_recPoints_det", "RRD recPoints det;Detector ID of the cluster;Counts", 540, -0.5, 539.5);
+ hist[1] = new TH2F("qaTRD_recPoints_amp", "TRD recPoints amp;Amplitude;??", 540, -0.5, 539, 200, -0.5, 199.5);
+ hist[2] = new TH1F("qaTRD_recPoints_npad", "TRD recPoints npad;Number of Pads;Counts", 12, -0.5, 11.5);
- hist[3] = new TH1D("qaTRD_recPoints_dist2", "TRD recPoints dist2;residuals [2pad];Counts", 100, -1, 1);
- hist[4] = new TH1D("qaTRD_recPoints_dist3", "TRD recPoints dist3;residuals [3pad];Counts", 100, -1, 1);
- hist[5] = new TH1D("qaTRD_recPoints_dist4", "TRD recPoints dist4;residuals [4pad];Counts", 100, -1, 1);
- hist[6] = new TH1D("qaTRD_recPoints_dist5", "TRD recPoints dist5;residuals [5pad];Counts", 100, -1, 1);
+ hist[3] = new TH1F("qaTRD_recPoints_dist2", "TRD recPoints dist2;residuals [2pad];Counts", 100, -1, 1);
+ hist[4] = new TH1F("qaTRD_recPoints_dist3", "TRD recPoints dist3;residuals [3pad];Counts", 100, -1, 1);
+ hist[5] = new TH1F("qaTRD_recPoints_dist4", "TRD recPoints dist4;residuals [4pad];Counts", 100, -1, 1);
+ hist[6] = new TH1F("qaTRD_recPoints_dist5", "TRD recPoints dist5;residuals [5pad];Counts", 100, -1, 1);
- hist[7] = new TH2D("qaTRD_recPoints_rowCol", "TRDrecPointsrowCol;row;col", 16, -0.5, 15.5, 145, -0.5, 144.5);
- hist[8] = new TH1D("qaTRD_recPoints_time", "TRDrecPoints time;time bin;Counts", kTimeBin, -0.5, kTimeBin-0.5);
- hist[9] = new TH1D("qaTRD_recPoints_nCls", "TRD recPoints nCls;number of clusters;Counts", 500, -0.5, 499.5);
+ hist[7] = new TH2F("qaTRD_recPoints_rowCol", "TRDrecPointsrowCol;row;col", 16, -0.5, 15.5, 145, -0.5, 144.5);
+ hist[8] = new TH1F("qaTRD_recPoints_time", "TRDrecPoints time;time bin;Counts", kTimeBin, -0.5, kTimeBin-0.5);
+ hist[9] = new TH1F("qaTRD_recPoints_nCls", "TRD recPoints nCls;number of clusters;Counts", 500, -0.5, 499.5);
- hist[10] = new TH3D("qaTRD_recPoints_sigTime", "TRD recPoints sigTime;chamber;time bin;signal",
+ hist[10] = new TH3F("qaTRD_recPoints_sigTime", "TRD recPoints sigTime;chamber;time bin;signal",
540, -0.5, 539.5, kTimeBin, -0.5, kTimeBin-0.5, 200, -0.5, 199.5);
hist[11] = new TProfile("qaTRD_recPoints_prf", "TRD recPoints prf;distance;center of gravity;Counts"
, 120, -0.6, 0.6, -1.2, 1.2, "");
- hist[12] = new TH1D("qaTRD_recPoints_ampMPV", "TRD recPoints ampMPV;amplitude MPV;Counts", 150, 0, 150);
- hist[13] = new TH1D("qaTRD_recPoints_ampSigma", "TRD recPoints ampSigma;amplitude Sigma;Counts", 200, 0, 200);
+ hist[12] = new TH1F("qaTRD_recPoints_ampMPV", "TRD recPoints ampMPV;amplitude MPV;Counts", 150, 0, 150);
+ hist[13] = new TH1F("qaTRD_recPoints_ampSigma", "TRD recPoints ampSigma;amplitude Sigma;Counts", 200, 0, 200);
// chamber by chamber
for(Int_t i=0; i<18; i++) {
- hist[14+i] = new TH2D(Form("qaTRD_recPoints_sigTime_sm%d",i), Form("sm%d;det;time bin",i),
+ hist[14+i] = new TH2F(Form("qaTRD_recPoints_sigTime_sm%d",i), Form("sm%d;det;time bin",i),
30, -0.5, 29.5, kTimeBin, -0.5, kTimeBin-0.5);
hist[14+i]->SetMinimum(0);
hist[14+i]->SetMaximum(150);
// time bin by time bin sm-by-sm
for(Int_t i=0; i<18; i++) {
- hist[14+18+i] = new TH1D(Form("qaTRD_recPoints_sigTimeShape_sm%d", i),
+ hist[14+18+i] = new TH1F(Form("qaTRD_recPoints_sigTimeShape_sm%d", i),
Form("sm%d;time bin;signal",i),
kTimeBin, -0.5, kTimeBin-0.5);
// str = 50
for(Int_t i=0; i<18; i++) {
- hist[50+i] = new TH1D(Form("qaTRD_recPoints_nCls_sm%d",i),
+ hist[50+i] = new TH1F(Form("qaTRD_recPoints_nCls_sm%d",i),
Form("sm%d;time bin;number of clusters",i),
kTimeBin, -0.5, kTimeBin-0.5);
}
// str = 68
for(Int_t i=0; i<18; i++) {
- hist[68+i] = new TH1D(Form("qaTRD_recPoints_totalCharge_sm%d", i),
+ hist[68+i] = new TH1F(Form("qaTRD_recPoints_totalCharge_sm%d", i),
Form("sm%d;time bin;total charge", i),
kTimeBin, -0.5, kTimeBin-0.5);
}
- hist[86] = new TH1D("qaTRD_recPoints_signal", "TRD recPoints signal;amplitude;Counts", 400, -0.5, 399.5);
- hist[87] = new TH2D("qaTRD_recPoints_detMap", "TRD recPoints detMap;sm;chamber;Counts", 18, -0.5, 17.5, 30, -0.5, 29.5);
+ hist[86] = new TH1F("qaTRD_recPoints_signal", "TRD recPoints signal;amplitude;Counts", 400, -0.5, 399.5);
+ hist[87] = new TH2F("qaTRD_recPoints_detMap", "TRD recPoints detMap;sm;chamber;Counts", 18, -0.5, 17.5, 30, -0.5, 29.5);
// amplitude as a function of the pad size
for(Int_t i=0; i<9; i++) {
- hist[88+i] = new TH1D(Form("qaTRD_recPoints_signalNpad_%d", i+2), Form("qaTRD_recPoints_signalNpad_%d;amplitude, ADC", i+2), 400, -0.5, 399.5);
+ hist[88+i] = new TH1F(Form("qaTRD_recPoints_signalNpad_%d", i+2), Form("qaTRD_recPoints_signalNpad_%d;amplitude, ADC", i+2), 400, -0.5, 399.5);
}
// one 2D histogram per chamber
// for(Int_t i=0; i<540; i++) {
- // hist[88+i] = new TH2D(Form("qaTRD_recPoints_map%d", i), ";col;row", 16, -0.5, 15.5, 144, -0.5, 143.5);
+ // hist[88+i] = new TH2F(Form("qaTRD_recPoints_map%d", i), ";col;row", 16, -0.5, 15.5, 144, -0.5, 143.5);
//}
//hist[i]->Sumw2();
Add2RecPointsList(hist[i], i, !expert, image);
}
+ //
+ ClonePerTrigClass(AliQAv1::kRECPOINTS); // this should be the last line
}
//____________________________________________________________________________
const Int_t kNhist = 7;
TH1 *hist[kNhist];
- hist[0] = new TH2D("qaTRD_raws_nADC","number of ADC channels;sector;detector", 18, -0.5, 17.5, 30, -0.5, 29.5);
- hist[1] = new TH2D("qaTRD_raws_nCls", "number of clusters;sector;detector", 18, -0.5, 17.5, 30, -0.5, 29.5);
- hist[2] = new TH2D("qaTRD_raws_meanSig", "mean signal;sector;detector", 18, -0.5, 17.5, 30, -0.5, 29.5);
+ hist[0] = new TH2F("qaTRD_raws_nADC","number of ADC channels;sector;detector", 18, -0.5, 17.5, 30, -0.5, 29.5);
+ hist[1] = new TH2F("qaTRD_raws_nCls", "number of clusters;sector;detector", 18, -0.5, 17.5, 30, -0.5, 29.5);
+ hist[2] = new TH2F("qaTRD_raws_meanSig", "mean signal;sector;detector", 18, -0.5, 17.5, 30, -0.5, 29.5);
- hist[3] = new TH1D("qaTRD_raws_ADC", "ADC amplitude;ADC counts", 100, -0.5, 99.5);
- hist[4] = new TH1D("qaTRD_raws_Cls", "Cluster amplitude; ADC counts", 100, -0.5, 199.5);
- hist[5] = new TH2D("qaTRD_raws_ClsTb", "Clusters vs Time Bin;time bin;amoplitude", 30, -0.5, 29.5, 200, -0.5, 199.5);
+ hist[3] = new TH1F("qaTRD_raws_ADC", "ADC amplitude;ADC counts", 100, -0.5, 99.5);
+ hist[4] = new TH1F("qaTRD_raws_Cls", "Cluster amplitude; ADC counts", 100, -0.5, 199.5);
+ hist[5] = new TH2F("qaTRD_raws_ClsTb", "Clusters vs Time Bin;time bin;amoplitude", 30, -0.5, 29.5, 200, -0.5, 199.5);
- hist[6] = new TH2D("qaTRD_raws_ClsAmpDet", ";detector;amplitude", 540, -0.5, 539.5, 100, 0, 200);
+ hist[6] = new TH2F("qaTRD_raws_ClsAmpDet", ";detector;amplitude", 540, -0.5, 539.5, 100, 0, 200);
/*
- hist[0] = new TH2D("qaTRD_raws_DataVolume", ";Sector;Data Volume, kB", 18, -0.5, 17.5, 100, 0, 30);
- hist[1] = new TH2D("qaTRD_raws_HC", "Data Headers;Sector;HC", 18, -0.5, 17.5, 60, -0.5, 59.5);
- hist[2] = new TH2D("qaTRD_raws_LME", "Link Monitor Error;Sector;HC", 18, -0.5, 17.5, 60, -0.5, 59.5);
+ hist[0] = new TH2F("qaTRD_raws_DataVolume", ";Sector;Data Volume, kB", 18, -0.5, 17.5, 100, 0, 30);
+ hist[1] = new TH2F("qaTRD_raws_HC", "Data Headers;Sector;HC", 18, -0.5, 17.5, 60, -0.5, 59.5);
+ hist[2] = new TH2F("qaTRD_raws_LME", "Link Monitor Error;Sector;HC", 18, -0.5, 17.5, 60, -0.5, 59.5);
- hist[3] = new TH1D("qaTRD_rawd_cls", "Clusters amplitude;ADC counts", 100, 0, 200);
- hist[4] = new TH2D("qaTRD_raws_clsTB", "amplitude - time bins;time bin;amplitude"
+ hist[3] = new TH1F("qaTRD_rawd_cls", "Clusters amplitude;ADC counts", 100, 0, 200);
+ hist[4] = new TH2F("qaTRD_raws_clsTB", "amplitude - time bins;time bin;amplitude"
30, -0.5, 29.5, 100, 0, 200);
- hist[5] = new TH2D("qaTRD_raws_clsSec", "amplitude in sectors;Sector;amplitude, ADCs"
+ hist[5] = new TH2F("qaTRD_raws_clsSec", "amplitude in sectors;Sector;amplitude, ADCs"
18, -0.5, 17.5);
*/
//hist[i]->Sumw2();
Add2RawsList(hist[i], i, !expert, image, !saveCorr);
}
-
+ //
+ ClonePerTrigClass(AliQAv1::kRAWS); // this should be the last line
}
//____________________________________________________________________________
//
Int_t nTracks = esd->GetNumberOfTracks();
- GetESDsData(0)->Fill(nTracks);
+ FillESDsData(0,nTracks);
// track loop
for (Int_t iTrack = 0; iTrack<nTracks; iTrack++) {
UInt_t u = 1;
UInt_t status = track->GetStatus();
for(Int_t bit=0; bit<32; bit++)
- if (u<<bit & status) GetESDsData(2)->Fill(bit);
+ if (u<<bit & status) FillESDsData(2,bit);
const Int_t knbits = 6;
Int_t bit[6] = {0,0,0,0,0,0};
for(Int_t b=0; b<knbits; b++) {
if (bit[b]) {
- GetESDsData(2*b+3)->Fill(pt);
- GetESDsData(2*b+4)->Fill(extZ);
+ FillESDsData(2*b+3,pt);
+ FillESDsData(2*b+4,extZ);
}
}
// clusters
for(Int_t b=0; b<3; b++)
- if (bit[3+b]) GetESDsData(b+15)->Fill(track->GetTRDncls0());
+ if (bit[3+b]) FillESDsData(b+15,track->GetTRDncls0());
// refitted only
if (!bit[4]) continue;
//fBudget->Fill(track->GetTRDBudget());
//fSignal->Fill(track->GetTRDsignal());
- GetESDsData(1)->Fill(sector);
- GetESDsData(18)->Fill(track->GetP(), track->GetTRDsignal());
+ FillESDsData(1,sector);
+ FillESDsData(18,track->GetP(), track->GetTRDsignal());
- GetESDsData(27)->Fill(track->GetTRDQuality());
- GetESDsData(28)->Fill(track->GetTRDBudget());
- GetESDsData(29)->Fill(track->GetTRDchi2());
- GetESDsData(30)->Fill(track->GetTRDTimBin(0));
- GetESDsData(31)->Fill(track->GetTRDntrackletsPID());
+ FillESDsData(27,track->GetTRDQuality());
+ FillESDsData(28,track->GetTRDBudget());
+ FillESDsData(29,track->GetTRDchi2());
+ FillESDsData(30,track->GetTRDTimBin(0));
+ FillESDsData(31,track->GetTRDntrackletsPID());
// dedx
for(Int_t j=0; j<6; j++) {
dedx += track->GetTRDslice(j, k-1);
}
- GetESDsData(41+k)->Fill(paramOut->GetP(), dedx/6.);
+ FillESDsData(41+k,paramOut->GetP(), dedx/6.);
}
// probabilities
if (status & AliESDtrack::kTRDpid) {
for(Int_t k=0; k<AliPID::kSPECIES; ++k)
- GetESDsData(36+k)->Fill(track->GetTRDpid(k));
+ FillESDsData(36+k,track->GetTRDpid(k));
}
// probabilities uniformity
if (track->GetTRDntrackletsPID() < 6) continue;
- GetESDsData(35)->Fill(paramOut->GetZ()/paramOut->GetX());
+ FillESDsData(35,paramOut->GetZ()/paramOut->GetX());
Int_t idx = 5 * sector + stack;
- GetESDsData(32)->Fill(idx); // all tracks
+ FillESDsData(32,idx); // all tracks
if (track->GetTRDpid(AliPID::kElectron) > 0.9)
- GetESDsData(33)->Fill(idx); // electrons only
+ FillESDsData(33,idx); // electrons only
/*
- hist[27] = new TH1D("qaTRD_esd_quality", ";quality", 120, 0, 12);
- hist[28] = new TH1D("qaTRD_esd_budget", ";NN", 110, -1000, 100);
- hist[29] = new TH1D("qaTRD_esd_chi2", ";chi2", 300, 0, 100);
- hist[30] = new TH1D("qaTRD_esd_timeBin", 7, -0.5, 6.5);
- hist[31] = new TH1D("qaTRD_esd_pidQuality", 7, -0.5, 6.5);
+ hist[27] = new TH1F("qaTRD_esd_quality", ";quality", 120, 0, 12);
+ hist[28] = new TH1F("qaTRD_esd_budget", ";NN", 110, -1000, 100);
+ hist[29] = new TH1F("qaTRD_esd_chi2", ";chi2", 300, 0, 100);
+ hist[30] = new TH1F("qaTRD_esd_timeBin", 7, -0.5, 6.5);
+ hist[31] = new TH1F("qaTRD_esd_pidQuality", 7, -0.5, 6.5);
*/
/*
*/
}
-
+ //
+ IncEvCountCycleESDs();
+ IncEvCountTotalESDs();
+ //
}
//______________________________________________________________________________
AliInfo("Making QA for Raws");
- // take histograms
-
- TH2D *mnADC = (TH2D*)GetRawsData(0);
- TH2D *mnCls = (TH2D*)GetRawsData(1);
- TH2D *mClsDet = (TH2D*)GetRawsData(2);
-
- TH1D *mADC = (TH1D*)GetRawsData(3);
- TH1D *mCls = (TH1D*)GetRawsData(4);
- TH2D *mClsTb = (TH2D*)GetRawsData(5);
-
- TH2D *mClsDetAmp = (TH2D*)GetRawsData(6);
+ // take histograms (RS: arrays of histos)
+ static TObjArray mnADC,mnCls,mClsDet,mADC,mCls,mClsTb,mClsDetAmp;
+ GetMatchingRawsData(0,&mnADC);
+ GetMatchingRawsData(1,&mnCls);
+ GetMatchingRawsData(2,&mClsDet);
+ GetMatchingRawsData(3,&mADC);
+ GetMatchingRawsData(4,&mCls);
+ GetMatchingRawsData(5,&mClsTb);
+ GetMatchingRawsData(6,&mClsDetAmp);
const Int_t baseline = 10;
rawReader->SelectEquipment(0, 1024, 1041);
rawReader->Select("TRD");
- AliTRDrawStreamBase::SetRawStreamVersion("FAST");
- AliTRDrawStreamBase *data = AliTRDrawStreamBase::GetRawStream(rawReader);
- data->SetSharedPadReadout(kFALSE);
+ AliTRDrawStream *data = new AliTRDrawStream(rawReader);
// build data manager
AliTRDdigitsManager *digitsManager;
while(indexes->NextRCIndex(row, col)) {
- mnADC->Fill(det/30, det%30);
+ for (int ih=mnADC.GetEntriesFast();ih--;) ((TH2*)mnADC.UncheckedAt(ih))->Fill(det/30, det%30);
for(Int_t tb = 0; tb < digits->GetNtime(); tb++) {
Int_t value = digits->GetData(row, col, tb);
- mADC->Fill(value);
+ for (int ih=mADC.GetEntriesFast();ih--;) ((TH1*)mADC.UncheckedAt(ih))->Fill(value);
// simple clusterizer
if (col < 1 || col > digits->GetNcol()-2) continue;
Int_t valueTot = value + valueL + valueR;
if (valueTot < 0) continue;
- mCls->Fill(valueTot);
- mClsTb->Fill(tb, valueTot);
- mClsDetAmp->Fill(det, valueTot);
+ for (int ih=mCls.GetEntriesFast();ih--;) ((TH1*)mCls.UncheckedAt(ih))->Fill(valueTot);
+ for (int ih=mClsTb.GetEntriesFast();ih--;) ((TH2*)mClsTb.UncheckedAt(ih))->Fill(tb, valueTot);
+ for (int ih=mClsDetAmp.GetEntriesFast();ih--;) ((TH2*)mClsDetAmp.UncheckedAt(ih))->Fill(det, valueTot);
if (valueTot < 200) {
- mnCls->Fill(det/30, det%30);
- mClsDet->Fill(det/30, det%30, valueTot);
+ for (int ih=mnCls.GetEntriesFast();ih--;) ((TH2*)mnCls.UncheckedAt(ih))->Fill(det/30, det%30);
+ for (int ih=mClsDet.GetEntriesFast();ih--;) ((TH2*)mClsDet.UncheckedAt(ih))->Fill(det/30, det%30, valueTot);
}
}
delete digitsManager;
delete data;
+ //
+ IncEvCountCycleRaws();
+ IncEvCountTotalRaws();
+ //
}
//____________________________________________________________________________
*/
/////
-
+ TObjArray *hists3D = GetMatchingRecPointsData(10); //RS no alias for 3d histo filling, to directly
+ //
for (Int_t iEntry = 0; iEntry < nEntries; iEntry++) {
// Import the tree
Int_t nPads = c->GetNPads();
nDet[iDet]++;
- GetRecPointsData(0)->Fill(iDet);
- GetRecPointsData(86)->Fill(c->GetQ());
- GetRecPointsData(1)->Fill(iDet, c->GetQ());
- GetRecPointsData(2)->Fill(nPads);
+ FillRecPointsData(0,iDet);
+ FillRecPointsData(86,c->GetQ());
+ FillRecPointsData(1,iDet, c->GetQ());
+ FillRecPointsData(2,nPads);
if (nPads < 6)
- GetRecPointsData(1+c->GetNPads())->Fill(c->GetCenter());
+ FillRecPointsData(1+c->GetNPads(),c->GetCenter());
if (nPads < 10)
- GetRecPointsData(88+nPads-2)->Fill(c->GetQ());
- else GetRecPointsData(96)->Fill(c->GetQ());
+ FillRecPointsData(88+nPads-2,c->GetQ());
+ else FillRecPointsData(96,c->GetQ());
//if (c->GetPadTime() < 5)
- ((TH2D*)GetRecPointsData(7))->Fill(c->GetPadRow(), c->GetPadCol());
- GetRecPointsData(8)->Fill(c->GetPadTime());
-
- ((TH3D*)GetRecPointsData(10))->Fill(iDet, c->GetPadTime(), c->GetQ());
+ FillRecPointsData(7,c->GetPadRow(), c->GetPadCol());
+ FillRecPointsData(8,c->GetPadTime());
+ //
+ if (hists3D) {
+ for (int ih=hists3D->GetEntriesFast();ih--;) {
+ TH3F * ahist = dynamic_cast<TH3F*>(hists3D->At(ih));
+ if (ahist) ahist->Fill(iDet, c->GetPadTime(), c->GetQ());
+ }
+ }
Int_t iSM = iDet / 30;
- GetRecPointsData(50+iSM)->Fill(c->GetPadTime());
- GetRecPointsData(68+iSM)->Fill(c->GetPadTime(), c->GetQ());
+ FillRecPointsData(50+iSM,c->GetPadTime());
+ FillRecPointsData(68+iSM,c->GetPadTime(), c->GetQ());
// total charge sm / det / timeBin
- //((TH2D*)GetRecPointsData(14+iSM))->Fill(iDet-iSM*30, c->GetPadTime(), c->GetQ());
+ // FillRecPointsData(14+iSM,iDet-iSM*30, c->GetPadTime(), c->GetQ());
// PRF for 2pad
if (sig[0] == 0 && sig[1] == 0 && sig[4] == 0 && sig[5] == 0 && sig[6] == 0)
frac = -1. * sig[2] / (sig[2] + sig[3]);
- if (frac > -10) ((TProfile*)GetRecPointsData(11))->Fill(c->GetCenter(), frac);
+ if (frac > -10) FillRecPointsData(11,c->GetCenter(), frac);
//}
}
/*
for(Int_t i=0; i<540; i++)
- if (nDet[i] > 0) GetRecPointsData(9)->Fill(nDet[i]);
+ if (nDet[i] > 0) FillRecPointsData(9,nDet[i]);
*/
- GetRecPointsData(9)->Fill(nCls);
+ FillRecPointsData(9,nCls);
delete clusterArray;
-
+ //
+ IncEvCountCycleRecPoints();
+ IncEvCountTotalRecPoints();
+ //
}
//____________________________________________________________________________
return !!obj;
}
//__________________________________________________________________________
-void AliTRDQADataMakerRec::BuildRatio(TH1D *ratio, TH1D *histN, TH1D*histD) {
+void AliTRDQADataMakerRec::BuildRatio(TH1 *ratio, TH1 *histN, TH1* histD) {
//
// Calculate the ratio of two histograms
// error are calculated assuming the histos have the same counts
}
//__________________________________________________________________________
-Int_t AliTRDQADataMakerRec::FillBits(TH1D *hist, Int_t code, Int_t offset) {
+Int_t AliTRDQADataMakerRec::FillBits(TH1F *hist, Int_t code, Int_t offset) {
Int_t nb = 0;
UInt_t test = 1;