//____________________________________________________________________________
AliT0QADataMakerRec::AliT0QADataMakerRec() :
- AliQADataMakerRec(AliQA::GetDetName(AliQA::kT0), "T0 Quality Assurance Data Maker")
+AliQADataMakerRec(AliQAv1::GetDetName(AliQAv1::kT0),
+ "T0 Quality Assurance Data Maker"),
+ fnEvent(0)
{
// ctor
- fnEvent=0;
for (Int_t i=0; i<6; i++) {
fNumTriggers[i]=0;
fNumTriggersCal[i]=0;
//____________________________________________________________________________
AliT0QADataMakerRec::AliT0QADataMakerRec(const AliT0QADataMakerRec& qadm) :
- AliQADataMakerRec()
+ AliQADataMakerRec(),
+ fnEvent(0)
+
{
//copy ctor
SetName((const char*)qadm.GetName()) ;
return *this;
}
//____________________________________________________________________________
-void AliT0QADataMakerRec::EndOfDetectorCycle(AliQA::TASKINDEX_t task, TObjArray * list)
+void AliT0QADataMakerRec::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** list)
{
//Detector specific actions at end of cycle
// do the QA checking
- AliQAChecker::Instance()->Run(AliQA::kT0, task, list) ;
- if ( task == AliQA::kRAWS ) {
- const Char_t *triggers[6] = {"mean", "vertex","ORA","ORC","central","semi-central"};
- for (Int_t itr=0; itr<6; itr++) {
- GetRawsData(197)->Fill(triggers[itr], fNumTriggersCal[itr]);
- GetRawsData(197)->SetBinContent(itr+1, fNumTriggersCal[itr]);
- }
- GetRawsData(205)->SetOption("COLZ");
- GetRawsData(206)->SetOption("COLZ");
- GetRawsData(207)->SetOption("COLZ");
- GetRawsData(205)->GetXaxis()->SetTitle("#PMT");
- GetRawsData(206)->GetXaxis()->SetTitle("#PMT");
- GetRawsData(205)->GetYaxis()->SetTitle("NeventsReg/Nevents");
- GetRawsData(206)->GetYaxis()->SetTitle("NeventsReg/Nevents");
- GetRawsData(207)->GetXaxis()->SetTitle("#PMT");
- GetRawsData(207)->GetYaxis()->SetTitle("Charge, #channels");
- }
- if ( task == AliQA::kRECPOINTS) {
- GetRecPointsData(0)->SetOption("COLZ");
- GetRecPointsData(1)->SetOption("COLZ");
- GetRecPointsData(0)->GetXaxis()->SetTitle("#PMT");
- GetRecPointsData(1)->GetXaxis()->SetTitle("#PMT");
- GetRecPointsData(0)->GetYaxis()->SetTitle("CFD time");
- GetRecPointsData(1)->GetYaxis()->SetTitle("Charge, #channels");
- }
+ AliQAChecker::Instance()->Run(AliQAv1::kT0, task, list) ;
+ for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
+ if (! AliQAv1::Instance(AliQAv1::GetDetIndex(GetName()))->IsEventSpecieSet(AliRecoParam::ConvertIndex(specie)) )
+ continue ;
+ SetEventSpecie(AliRecoParam::ConvertIndex(specie)) ;
+ if ( task == AliQAv1::kRAWS ) {
+ const Char_t *triggers[6] = {"mean", "vertex","ORA","ORC","central","semi-central"};
+ for (Int_t itr=0; itr<6; itr++) {
+ GetRawsData(197)->Fill(triggers[itr], fNumTriggersCal[itr]);
+ GetRawsData(197)->SetBinContent(itr+1, fNumTriggersCal[itr]);
+ }
+ }
+ }
}
//____________________________________________________________________________
void AliT0QADataMakerRec::InitRaws()
{
// create Raw histograms in Raw subdir
+ const Bool_t expert = kTRUE ;
+ const Bool_t saveCorr = kTRUE ;
+ const Bool_t image = kTRUE ;
+
TString timename, ampname, qtcname, ledname;
TString timeCalname, ampCalname, ledCalname, qtcCalname;
TH1F* fhRefPoint = new TH1F("hRefPoint","Ref Point", 10,1252170, 1252180);
- Add2RawsList( fhRefPoint,0);
+ Add2RawsList( fhRefPoint,0, !expert, image, !saveCorr);
TH1F *fhRawCFD[24]; TH1F * fhRawLEDamp[24];
TH1F *fhRawQTC[24]; TH1F * fhRawLED[24];
ampname += i;
qtcname += i;
ledname += i;
- fhRawCFD[i] = new TH1F(timename.Data(), timename.Data(),500,-250,250);
- Add2RawsList( fhRawCFD[i],i+1);
- fhRawLED[i] = new TH1F(ledname.Data(), ledname.Data(),2000,-1000,1000);
- Add2RawsList( fhRawLED[i],i+24+1);
- fhRawLEDamp[i] = new TH1F(ampname.Data(), ampname.Data(),100,300,600);
- Add2RawsList( fhRawLEDamp[i],i+48+1);
- fhRawQTC[i] = new TH1F(qtcname.Data(), qtcname.Data(),1500,1000,7000);
- Add2RawsList( fhRawQTC[i],i+72+1);
+ fhRawCFD[i] = new TH1F(timename.Data(), Form("%s;Time [ns];Counts", timename.Data()),500,-250,250);
+ Add2RawsList( fhRawCFD[i],i+1, !expert, image, !saveCorr);
+ fhRawLED[i] = new TH1F(ledname.Data(), Form("%s;Time [ns];Counts", ledname.Data()),2000,-1000,1000);
+ Add2RawsList( fhRawLED[i],i+24+1, !expert, image, !saveCorr);
+ fhRawLEDamp[i] = new TH1F(ampname.Data(), Form("%s;Amplitude [ADC counts];Counts", ampname.Data()),100,300,600);
+ Add2RawsList( fhRawLEDamp[i],i+48+1, !expert, image, !saveCorr);
+ fhRawQTC[i] = new TH1F(qtcname.Data(), Form("%s;Charge [??];Counts", qtcname.Data()),1500,1000,7000);
+ Add2RawsList( fhRawQTC[i],i+72+1, !expert, image, !saveCorr);
}
- // const Bool_t saveForCorr = kTRUE;
- const Bool_t saveForCorr = kFALSE;
- TH1F* fhRawTrigger = new TH1F("hRawTrigger"," phys triggers",5,0,5);
- Add2RawsList(fhRawTrigger ,97);
+ 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);
- Add2RawsList( fhRawMean,98);
- TH1F* fhRawVertex = new TH1F("hRawVertex","online vertex signal", 100,0,600);
- Add2RawsList( fhRawVertex,99);
- TH1F* fhRawORA = new TH1F("hRawORA","online OR A", 100,2500,2800);
- Add2RawsList( fhRawORA,100);
- TH1F* fhRawORC = new TH1F("hRawORC","online OR C", 100,2000,2300);
- Add2RawsList( fhRawORC,101);
+ 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);
+ TH1F* fhRawORA = new TH1F("hRawORA","online OR A;??;??", 100,2500,2800);
+ Add2RawsList( fhRawORA,100, !expert, image, !saveCorr);
+ TH1F* fhRawORC = new TH1F("hRawORC","online OR C;??;??", 100,2000,2300);
+ Add2RawsList( fhRawORC,101, !expert, image, !saveCorr);
for (Int_t i=0; i<24; i++)
{
ledCalname += i;
ampCalname += i;
qtcCalname += i;
- fhRawCFDcal[i] = new TH1F(timeCalname.Data(), timeCalname.Data(),2000,-1000,1000);
- Add2RawsList( fhRawCFDcal[i],101+i+1);
- fhRawLEDcal[i] = new TH1F(ledCalname.Data(), ledCalname.Data(),2000,-1000,1000);
- Add2RawsList( fhRawLEDcal[i],101+i+24+1);
- fhRawLEDampcal[i] = new TH1F(ampCalname.Data(), ampCalname.Data(),300,300,600);
- Add2RawsList( fhRawLEDampcal[i],101+i+48+1);
- fhRawQTCcal[i] = new TH1F(qtcCalname.Data(), qtcCalname.Data(),1000,0,7000);
- Add2RawsList( fhRawQTCcal[i],101+i+72+1);
+ fhRawCFDcal[i] = new TH1F(timeCalname.Data(), Form("%s;Time [ns];Counts", timeCalname.Data()),2000,-1000,1000);
+ Add2RawsList( fhRawCFDcal[i],101+i+1, !expert, image, !saveCorr);
+ fhRawLEDcal[i] = new TH1F(ledCalname.Data(), Form("%s;Time [ns];Counts", ledCalname.Data()),2000,-1000,1000);
+ Add2RawsList( fhRawLEDcal[i],101+i+24+1, !expert, image, !saveCorr);
+ fhRawLEDampcal[i] = new TH1F(ampCalname.Data(), Form("%s;Amplitude [ADC counts];Counts", ampCalname.Data()),300,300,600);
+ Add2RawsList( fhRawLEDampcal[i],101+i+48+1, !expert, image, !saveCorr);
+ fhRawQTCcal[i] = new TH1F(qtcCalname.Data(), Form("%s;Charge [??];Counts",qtcCalname.Data()),1000,0,7000);
+ Add2RawsList( fhRawQTCcal[i],101+i+72+1, !expert, image, !saveCorr);
}
TH1F* fhRawTriggerCal = new TH1F("hRawTriggerCal"," laser triggers",6,0,6);
- Add2RawsList(fhRawTriggerCal ,197 ,saveForCorr);
+ Add2RawsList(fhRawTriggerCal ,197 , !expert, !image, saveCorr);
- TH1F* fhRawMeanCal = new TH1F("hRawMeanCal","online mean signal, calibration event",
+ TH1F* fhRawMeanCal = new TH1F("hRawMeanCal","online mean signal, calibration event;??;??",
10000,0,10000);
Add2RawsList( fhRawMeanCal,198);
TH1F* fhRawVertexCal = new TH1F("hRawVertexCal","online vertex signal, calibration even ",
10000,0,10000);
- Add2RawsList( fhRawVertexCal,199);
- TH1F* fhRawORAcal = new TH1F("hRawORAcal","online OR A", 10000,0,10000);
- Add2RawsList( fhRawORAcal,200 );
- TH1F* fhRawORCcal = new TH1F("hRawORCcal","online OR C", 10000,0,10000);
- Add2RawsList( fhRawORCcal,201);
- TH1F* fhMultcal = new TH1F("hMultcal","full mulltiplicity", 10000,0,10000);
- Add2RawsList( fhMultcal,202 );
- TH1F* fhMultScal = new TH1F("hMultScal","full multiplicity with semi-central trigger",
+ Add2RawsList( fhRawVertexCal,199, !expert, image, !saveCorr);
+ TH1F* fhRawORAcal = new TH1F("hRawORAcal","online OR A;??;??", 10000,0,10000);
+ Add2RawsList( fhRawORAcal,200, !expert, image, !saveCorr );
+ TH1F* fhRawORCcal = new TH1F("hRawORCcal","online OR C;??;??", 10000,0,10000);
+ Add2RawsList( fhRawORCcal,201, !expert, image, !saveCorr);
+ TH1F* fhMultcal = new TH1F("hMultcal","full mulltiplicity;Multiplicity;Entries", 10000,0,10000);
+ Add2RawsList( fhMultcal,202, !expert, image, !saveCorr );
+ TH1F* fhMultScal = new TH1F("hMultScal","full multiplicity with semi-central trigger;Multiplicity;Entries",
10000,0,10000);
- Add2RawsList( fhMultScal,203);
- TH1F* fhMultCcal = new TH1F("hMultCcal","full multiplicity with central trigger",
+ Add2RawsList( fhMultScal,203, !expert, image, !saveCorr);
+ TH1F* fhMultCcal = new TH1F("hMultCcal","full multiplicity with central trigger;Multiplicity;Entries",
1000,0,10000);
- Add2RawsList( fhMultCcal,204);
+ Add2RawsList( fhMultCcal,204, !expert, image, !saveCorr);
// TH2F* fhEffCFD = new TH2F("hEffCFD"," CFD time",24, 0 ,24,
// 100,-500,500);
- TH2F* fhEffCFD = new TH2F("hEffCFD"," CFD time",24, 0 ,24, 50, 0,5);
-
- Add2RawsList( fhEffCFD,205, saveForCorr);
- TH2F* fhEffLED = new TH2F("hEffLED","LED time",24, 0 ,24,
+ TH2F* fhEffCFD = new TH2F("hEffCFD"," CFD time;Time [ns];Counts",24, 0 ,24, 50, 0,5);
+ fhEffCFD->SetOption("COLZ");
+ Add2RawsList( fhEffCFD,205, !expert, !image, saveCorr);
+ TH2F* fhEffLED = new TH2F("hEffLED","LED time;Time [ns];Counts",24, 0 ,24,
100, 0, 5);
//100,-500,500);
-
- Add2RawsList( fhEffLED,206, saveForCorr);
- TH2F* fhEffQTC = new TH2F("hEffQTC","QTC amplitude",24, 0 ,24,
- 100,0,7000);
- Add2RawsList( fhEffQTC,207, saveForCorr);
- //yeys guide line
+ fhEffLED->SetOption("COLZ");
+ Add2RawsList( fhEffLED,206, !expert, !image, saveCorr);
+ TH2F* fhEffQTC = new TH2F("hEffQTC","QTC amplitude%s;Amplitude [ADC counts];Counts",24, 0 ,24, 100,0,7000);
+ fhEffQTC->SetOption("COLZ");
+ Add2RawsList( fhEffQTC,207, !expert, !image, saveCorr);
+ //eyes guide line
// TH2F* fhLineQTC = new TH2F("hLineQTC","QTC amplitude boeder",24, 0 ,24,
// 100,0,7000);
-// Add2RawsList( fhLineQTC,208, saveForCorr);
+// Add2RawsList( fhLineQTC,208, !expert, !image, saveCorr);
+ const Char_t *triggers[6] = {"mean", "vertex","ORA","ORC","central","semi-central"};
+ for (Int_t itr=0; itr<6; itr++) {
+ GetRawsData(197)->Fill(triggers[itr], fNumTriggersCal[itr]);
+ GetRawsData(197)->SetBinContent(itr+1, fNumTriggersCal[itr]);
+ }
}
+//____________________________________________________________________________
+void AliT0QADataMakerRec::InitDigits()
+{
+ // create Digits histograms in Digits subdir
+ const Bool_t expert = kTRUE ;
+ const Bool_t image = kTRUE ;
+
+ TH2F * fhDigCFD = new TH2F("fhDigCFD", " CFD digits; #PMT; CFD digits[#channels]",25,-0.5,24.5,100,0,1000);
+ fhDigCFD->SetOption("COLZ");
+ 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);
+ fhDigLEDamp->SetOption("COLZ");
+ Add2DigitsList( fhDigLEDamp,1, !expert, image);
+ TH2F * fhDigQTC = new TH2F("fhDigQTC", " QTC digits; #PMT; QTC amplitude",25,-0.5,24.5,100,100,10000);
+ fhDigQTC->SetOption("COLZ");
+ Add2DigitsList( fhDigQTC,2, !expert, image);}
+
//____________________________________________________________________________
void AliT0QADataMakerRec::InitRecPoints()
{
// create cluster histograms in RecPoint subdir
+ const Bool_t expert = kTRUE ;
+ const Bool_t image = kTRUE ;
-
-
- TH2F* fhRecCFD = new TH2F("hRecCFD"," CFD time",24, 0 ,24,
+ TH2F* fhRecCFD = new TH2F("hRecCFD"," CFD time;Time [ns];Counts",24, 0 ,24,
100,-50,50);
- Add2RecPointsList ( fhRecCFD,0);
+ fhRecCFD->SetOption("COLZ");
+ Add2RecPointsList ( fhRecCFD,0, !expert, image);
- TH2F* fhRecAmpDiff = new TH2F("hRecAmpDiff"," LED-CFD min QTC amplitude",
+ TH2F* fhRecAmpDiff = new TH2F("hRecAmpDiff"," LED-CFD min QTC amplitude;Amplitude [ADC counts];Counts",
24, 0 ,24, 200,-10,10);
- Add2RecPointsList (fhRecAmpDiff, 1);
+ fhRecAmpDiff->SetOption("COLZ");
+ Add2RecPointsList (fhRecAmpDiff, 1, !expert, image);
- TH1F *fhMean = new TH1F("hMean","online - rec mean",1000, -5000, 5000);
- Add2RecPointsList ( fhMean,2);
- }
+ TH1F *fhMean = new TH1F("hMean","online - rec mean;online - rec mean[#channels];",2000, -1000, 1000);
+ Add2RecPointsList ( fhMean,2, !expert, image);
+
+}
//____________________________________________________________________________
void AliT0QADataMakerRec::InitESDs()
{
//create ESDs histograms in ESDs subdir
-
- TH1F *fhESDMean = new TH1F("hESDmean"," ESD mean",100,2400,2500);
- Add2ESDsList(fhESDMean, 0) ;
- TH1F * fhESDVertex = new TH1F("hESDvertex","ESD vertex",100,-50,50);
- Add2ESDsList(fhESDVertex, 1) ;
+ const Bool_t expert = kTRUE ;
+ const Bool_t image = kTRUE ;
+
+ TH1F *fhESDMean = new TH1F("hESDmean"," ESD mean; mean time[%channels]",1000,0,1000);
+ Add2ESDsList(fhESDMean, 0, !expert, image) ;
+ TH1F * fhESDVertex = new TH1F("hESDvertex","ESDvertex; vertex[cm];",82,-30,30);
+ Add2ESDsList(fhESDVertex, 1, !expert, image) ;
}
//____________________________________________________________________________
void AliT0QADataMakerRec::MakeRaws( AliRawReader* rawReader)
{
+
+ // Check id histograms already created for this Event Specie
+ if ( ! GetRawsData(0) )
+ InitRaws() ;
+
rawReader->Reset() ;
//fills QA histos for RAW
Int_t shift=0;
AliT0RawReader *start = new AliT0RawReader(rawReader);
// start->Next();
if (! start->Next())
- AliDebug(1,Form(" no raw data found!!"));
+ AliDebug(AliQAv1::GetQADebugLevel(),Form(" no raw data found!!"));
else
{
fnEvent++;
}
}
-
-
- delete start;
+ delete start;
}
}
+//____________________________________________________________________________
+void AliT0QADataMakerRec::MakeDigits( TTree *digitsTree)
+{
+ //fills QA histos for Digits
+
+ // Check id histograms already created for this Event Specie
+ if ( ! GetDigitsData(0) )
+ InitDigits() ;
+
+ TArrayI *digCFD = new TArrayI(24);
+ TArrayI *digLED = new TArrayI(24);
+ TArrayI *digQT0 = new TArrayI(24);
+ TArrayI *digQT1 = new TArrayI(24);
+ Int_t refpoint=0;
+
+ TBranch *brDigits=digitsTree->GetBranch("T0");
+ AliT0digit *fDigits = new AliT0digit() ;
+ if (brDigits) {
+ brDigits->SetAddress(&fDigits);
+ }else{
+ AliError(Form("EXEC Branch T0 digits not found"));
+ return;
+ }
+ digitsTree->GetEvent(0);
+ digitsTree->GetEntry(0);
+ brDigits->GetEntry(0);
+ fDigits->GetTimeCFD(*digCFD);
+ fDigits->GetTimeLED(*digLED);
+ fDigits->GetQT0(*digQT0);
+ fDigits->GetQT1(*digQT1);
+ refpoint = fDigits->RefPoint();
+ for (Int_t i=0; i<24; i++)
+ {
+ 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(2) -> Fill(i, (digQT1->At(i) - digQT0->At(i)));
+ }
+ }
+
+ delete digCFD;
+ delete digLED;
+ delete digQT0;
+ delete digQT1;
+
+}
+
//____________________________________________________________________________
void AliT0QADataMakerRec::MakeRecPoints(TTree * clustersTree)
{
//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;
}
+ // Check id histograms already created for this Event Specie
+ if ( ! GetRecPointsData(0) )
+ InitRecPoints() ;
+
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
+
+ // Check id histograms already created for this Event Specie
+ if ( ! GetESDsData(0) )
+ InitESDs() ;
GetESDsData(0) -> Fill(esd->GetT0());
GetESDsData(1)-> Fill(esd->GetT0zVertex());
-
+
}