GetRawsData(207)->GetXaxis()->SetTitle("#PMT");
GetRawsData(207)->GetYaxis()->SetTitle("Charge, #channels");
}
+ if ( task == AliQAv1::kDIGITSR) {
+ GetDigitsData(0)->SetOption("COLZ");
+ GetDigitsData(1)->SetOption("COLZ");
+ GetDigitsData(2)->SetOption("COLZ");
+
+ }
if ( task == AliQAv1::kRECPOINTS) {
GetRecPointsData(0)->SetOption("COLZ");
GetRecPointsData(1)->SetOption("COLZ");
TH1F* fhRawTrigger = new TH1F("hRawTrigger"," phys triggers;Trigger #;Counts",5,0,5);
Add2RawsList(fhRawTrigger ,97, !expert, image, !saveCorr);
- TH1F* fhRawMean = new TH1F("hRawMean","online mean signal;??;??", 100,2400,2500);
+ TH1F* fhRawMean = new TH1F("hRawMean","online mean signal;;??", 100,2400,2500);
Add2RawsList( fhRawMean,98, !expert, image, !saveCorr);
TH1F* fhRawVertex = new TH1F("hRawVertex","online vertex signal;??;??", 100,0,600);
Add2RawsList( fhRawVertex,99, !expert, image, !saveCorr);
//100,-500,500);
Add2RawsList( fhEffLED,206, !expert, !image, saveCorr);
- TH2F* fhEffQTC = new TH2F("hEffQTC","QTC amplitude%s;Amplitude [ADC counts];Counts",24, 0 ,24,
- 100,0,7000);
+ TH2F* fhEffQTC = new TH2F("hEffQTC","QTC amplitude%s;Amplitude [ADC counts];Counts",24, 0 ,24, 100,0,7000);
Add2RawsList( fhEffQTC,207, !expert, !image, saveCorr);
- //yeys guide line
+ //eyes guide line
// TH2F* fhLineQTC = new TH2F("hLineQTC","QTC amplitude boeder",24, 0 ,24,
// 100,0,7000);
// Add2RawsList( fhLineQTC,208, !expert, !image, saveCorr);
const Bool_t expert = kTRUE ;
const Bool_t image = kTRUE ;
- TH2F * fhDigCFD = new TH2F("fhDigCFD", " CFD digits;something;something else",25,-0.5,24.5,100,100,1000);
- Add2DigitsList( fhDigCFD,0);
- TH2F *fhDigLEDamp = new TH2F("fhDigLEDamp", " LED-CFD digits;something;something else",25,-0.5,24.5,100,100,1000);
+ TH2F * fhDigCFD = new TH2F("fhDigCFD", " CFD digits; #PMT; CFD digits[#channels]",25,-0.5,24.5,100,0,1000);
+ Add2DigitsList( fhDigCFD,0, !expert, image);
+ TH2F *fhDigLEDamp = new TH2F("fhDigLEDamp", " LED-CFD digits; #PMT; LED-CFD amplitude ",25,-0.5,24.5,100,100,1000);
Add2DigitsList( fhDigLEDamp,1, !expert, image);
- TH2F * fhDigQTC = new TH2F("fhDigQTC", " QTC digits;something;something else",25,-0.5,24.5,100,100,1000);
+ TH2F * fhDigQTC = new TH2F("fhDigQTC", " QTC digits; #PMT; QTC amplitude",25,-0.5,24.5,100,100,10000);
Add2DigitsList( fhDigQTC,2, !expert, image);}
//____________________________________________________________________________
24, 0 ,24, 200,-10,10);
Add2RecPointsList (fhRecAmpDiff, 1, !expert, image);
- TH1F *fhMean = new TH1F("hMean","online - rec mean;??;??",1000, -5000, 5000);
+ TH1F *fhMean = new TH1F("hMean","online - rec mean;online - rec mean[#channels];",2000, -1000, 1000);
Add2RecPointsList ( fhMean,2, !expert, image);
}
const Bool_t expert = kTRUE ;
const Bool_t image = kTRUE ;
- TH1F *fhESDMean = new TH1F("hESDmean"," ESD mean;??;??",100,2400,2500);
+ TH1F *fhESDMean = new TH1F("hESDmean"," ESD mean; mean time[%channels]",1000,0,1000);
Add2ESDsList(fhESDMean, 0, !expert, image) ;
- TH1F * fhESDVertex = new TH1F("hESDvertex","ESD vertex;??;??",82,-30,30);
+ TH1F * fhESDVertex = new TH1F("hESDvertex","ESDvertex; vertex[cm];",82,-30,30);
Add2ESDsList(fhESDVertex, 1, !expert, image) ;
}
}
-
-
- delete start;
+ delete start;
}
}
if (digCFD->At(i)>0) {
Int_t cfd=digCFD->At(i)- refpoint;
GetDigitsData(0) ->Fill(i,cfd);
- GetDigitsData(1) -> Fill(i,(digLED->At(i) - digCFD->At(i)));
+ GetDigitsData(1) -> Fill(i, (digLED->At(i) - digCFD->At(i)));
GetDigitsData(2) -> Fill(i, (digQT1->At(i) - digQT0->At(i)));
}
}
//fills QA histos for clusters
AliT0RecPoint* frecpoints= new AliT0RecPoint ();
- if (!frecpoints) {
+ if (!frecpoints) {
AliError("Reconstruct Fill ESD >> no recpoints found");
return;
}
if (brRec) {
brRec->SetAddress(&frecpoints);
}else{
- AliError(Form("EXEC Branch T0 rec not found "));
- return;
+ AliError(Form("EXEC Branch T0 rec not found "));
+ return;
}
brRec->GetEntry(0);
-
+
for ( Int_t i=0; i<24; i++) {
if(i<12)
GetRecPointsData(0) -> Fill(i, frecpoints -> GetTime(i) - frecpoints -> GetTime(0));
GetRecPointsData(1) -> Fill( i, frecpoints -> GetAmp(i) - frecpoints->AmpLED(i));
}
Double_t mmm=frecpoints->GetOnlineMean()- frecpoints->GetMeanTime();
- GetRecPointsData(2) ->Fill(mmm);
- // printf(" AliT0QADataMakerRec: diff mean %f \n",mmm );
- // GetRecPointsData(73) ->Fill(frecpoints->GetMeanTime());
+ GetRecPointsData(2) ->Fill(mmm);
-
-
}
//____________________________________________________________________________
void AliT0QADataMakerRec::MakeESDs(AliESDEvent * esd)
{
//fills QA histos for ESD
-
+
GetESDsData(0) -> Fill(esd->GetT0());
GetESDsData(1)-> Fill(esd->GetT0zVertex());
-
+
}
//---
// Produces the data needed to calculate the quality assurance.
-// All data must be mergeable objects.
-// A. Mastroserio
+// Alla.Maevskaya@cern.ch
+//
//---
// --- ROOT system ---
#include "AliQAChecker.h"
#include "AliT0RawReader.h"
+#include "iostream.h"
+
ClassImp(AliT0QADataMakerSim)
//____________________________________________________________________________
{
// ctor
- /*
- for(Int_t i=0; i<24; i++) {
- fhHitsTime[i]=0x0;
- fhDigCFD[i]=0x0;
- fhDigLEDamp[i]=0x0;
- fhRecCFD[i]=0x0;
- fhRecLEDamp[i]=0x0;
- fhRecQTC[i]=0x0;
- }
- */
// fDetectorDir = fOutput->GetDirectory(GetName()) ;
// if (!fDetectorDir)
// fDetectorDir = fOutput->mkdir(GetName()) ;
AliQADataMakerSim()
{
//copy ctor
- /*
- for(Int_t i=0; i<24; i++) {
- fhHitsTime[i]=0x0;
- fhDigCFD[i]=0x0;
- fhDigLEDamp[i]=0x0;
- fhRecCFD[i]=0x0;
- fhRecLEDamp[i]=0x0;
- fhRecQTC[i]=0x0;
- }
- */
SetName((const char*)qadm.GetName()) ;
SetTitle((const char*)qadm.GetTitle());
}
// Equal operator.
this->~AliT0QADataMakerSim();
new(this) AliT0QADataMakerSim(qadm);
- return *this;
+ return *this;
}
//____________________________________________________________________________
void AliT0QADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** list)
TString timename;
- TH2F *fhHitsTimeA = new TH2F("hHitsTimeA", "Hits Efficiency;Time [ns];Efficiency [%]", 25, 0, 25, 100,12,15 );
+ TH2F *fhHitsTimeA = new TH2F("hHitsTimeA", "Hits Efficiency;#PMT; Time [ns];", 13, 12, 25, 100,12,15 );
Add2HitsList(fhHitsTimeA,0, !expert, image);
- TH2F *fhHitsTimeC = new TH2F("hHitsTimeC", "Hits Efficiency;Time [ns];Efficiency [%]", 25, 0, 25, 100,2,5 );
+ TH2F *fhHitsTimeC = new TH2F("hHitsTimeC", "Hits Efficiency;#PMT; Time [ns];", 13, 0, 13, 100,2,5 );
Add2HitsList(fhHitsTimeC,1, !expert, image);
}
const Bool_t expert = kTRUE ;
const Bool_t image = kTRUE ;
- TH2F * fhDigCFD = new TH2F("fhDigCFD", " CFD digits;something;something else",25,-0.5,24.5,100,100,1000);
+ TH2F * fhDigCFD = new TH2F("fhDigCFD", " CFD digits; #PMT; CFD time [#channel]",25,-0.5,24.5,100,0,1000);
Add2DigitsList( fhDigCFD,0);
- TH2F *fhDigLEDamp = new TH2F("fhDigLEDamp", " LED-CFD digits;something;something else",25,-0.5,24.5,100,100,1000);
+ TH2F *fhDigLEDamp = new TH2F("fhDigLEDamp", " LED-CFD digits; #PMT; amplitude LED-CFD [#channel]",25,-0.5,24.5,100,100,1000);
Add2DigitsList( fhDigLEDamp,1, !expert, image);
- TH2F * fhDigQTC = new TH2F("fhDigQTC", " QTC digits;something;something else",25,-0.5,24.5,100,100,1000);
+ TH2F * fhDigQTC = new TH2F("fhDigQTC", " QTC digits; #PMT; amplitude QTC [#channel]",25,-0.5,24.5,200,500,10000);
Add2DigitsList( fhDigQTC,2, !expert, image);
fDigits->GetQT0(*digQT0);
fDigits->GetQT1(*digQT1);
refpoint = fDigits->RefPoint();
+
for (Int_t i=0; i<24; i++)
{
if (digCFD->At(i)>0) {
GetDigitsData(0) ->Fill(i,cfd);
GetDigitsData(1) -> Fill(i,(digLED->At(i) - digCFD->At(i)));
GetDigitsData(2) -> Fill(i, (digQT1->At(i) - digQT0->At(i)));
+
}
}