// do the QA checking
AliQAChecker::Instance()->Run(AliQAv1::kT0, task, list) ;
for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
- SetEventSpecie(specie) ;
+ if (! IsValidEventSpecie(specie, list))
+ 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]);
}
- 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 == AliQAv1::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");
}
}
}
ampname += i;
qtcname += i;
ledname += i;
- fhRawCFD[i] = new TH1F(timename.Data(), timename.Data(),500,-250,250);
+ 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(), ledname.Data(),2000,-1000,1000);
+ 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(), ampname.Data(),100,300,600);
+ 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(), qtcname.Data(),1500,1000,7000);
+ fhRawQTC[i] = new TH1F(qtcname.Data(), Form("%s;Charge [??];Counts", qtcname.Data()),1500,1000,7000);
Add2RawsList( fhRawQTC[i],i+72+1, !expert, image, !saveCorr);
}
- TH1F* fhRawTrigger = new TH1F("hRawTrigger"," phys triggers",5,0,5);
+ 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);
+ 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);
+ 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);
+ 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);
+ 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(), ledCalname.Data(),2000,-1000,1000);
+ 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(), ampCalname.Data(),300,300,600);
+ 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(), qtcCalname.Data(),1000,0,7000);
+ 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 , !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, !expert, image, !saveCorr);
- TH1F* fhRawORAcal = new TH1F("hRawORAcal","online OR A", 10000,0,10000);
+ 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);
+ TH1F* fhRawORCcal = new TH1F("hRawORCcal","online OR C;??;??", 10000,0,10000);
Add2RawsList( fhRawORCcal,201, !expert, image, !saveCorr);
- TH1F* fhMultcal = new TH1F("hMultcal","full mulltiplicity", 10000,0,10000);
+ 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",
+ TH1F* fhMultScal = new TH1F("hMultScal","full multiplicity with semi-central trigger;Multiplicity;Entries",
10000,0,10000);
Add2RawsList( fhMultScal,203, !expert, image, !saveCorr);
- TH1F* fhMultCcal = new TH1F("hMultCcal","full multiplicity with central trigger",
+ TH1F* fhMultCcal = new TH1F("hMultCcal","full multiplicity with central trigger;Multiplicity;Entries",
1000,0,10000);
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);
-
+ 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",24, 0 ,24,
+ TH2F* fhEffLED = new TH2F("hEffLED","LED time;Time [ns];Counts",24, 0 ,24,
100, 0, 5);
//100,-500,500);
-
+ fhEffLED->SetOption("COLZ");
Add2RawsList( fhEffLED,206, !expert, !image, saveCorr);
- TH2F* fhEffQTC = new TH2F("hEffQTC","QTC amplitude",24, 0 ,24,
- 100,0,7000);
+ 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);
- //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 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()
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);
+ 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);
+ fhRecAmpDiff->SetOption("COLZ");
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);
- }
+
+}
//____________________________________________________________________________
void AliT0QADataMakerRec::InitESDs()
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",100,-50,50);
+ TH1F * fhESDVertex = new TH1F("hESDvertex","ESDvertex; vertex[cm];",82,-30,30);
Add2ESDsList(fhESDVertex, 1, !expert, image) ;
//____________________________________________________________________________
void AliT0QADataMakerRec::MakeRaws( AliRawReader* rawReader)
{
+
rawReader->Reset() ;
//fills QA histos for RAW
Int_t shift=0;
}
}
-
-
- delete start;
+ delete start;
}
}
+//____________________________________________________________________________
+void AliT0QADataMakerRec::MakeDigits( TTree *digitsTree)
+{
+ //fills QA histos for Digits
+
+ 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;
}
-
+
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());
-
+
}