///////////////////////////////////////////////////////////////////////
/*
+Modified by fbellini on 01/11/2011
+- removed TLines as functions
+- changed shifters plots for 2012 DQM
+
+Modified by fbellini on 01/11/2011
+- added histograms for LTM monitoring
+- fix for coverity
+
Modified by fbellini on 17/10/2011
- fix for memory leak in constructor
- added methods to read histos ranges from config file in DQM
- Added CheckVolumeID() and CheckEquipID() methods
- Updated Raw QA
*/
+#include <iostream>
+#include <fstream>
#include <TClonesArray.h>
#include <TH1F.h>
const Float_t AliTOFQADataMakerRec::fgkNbinsWidthTime=2.44;//width of bins in time plot
Float_t AliTOFQADataMakerRec::fgRangeMinTime=0.0;//range min in time plot
Float_t AliTOFQADataMakerRec::fgRangeMaxTime=620.0; //range max in time plot
-Int_t AliTOFQADataMakerRec::fgCutNmaxFiredMacropad=5;//cut on number of max fired macropad
+Int_t AliTOFQADataMakerRec::fgCutNmaxFiredMacropad=50;//cut on number of max fired macropad
const Int_t AliTOFQADataMakerRec::fgkFiredMacropadLimit=50;//cut on number of max fired macropad
// fLineExpTimeMax = new TLine(250., 0., 250., 0.);
// fLineExpTotMin = new TLine(5., 0., 5., 0.);
// fLineExpTotMax = new TLine(20., 0., 20., 0.);
- for (Int_t sm=0;sm<17;sm++){
+ /*
+ for (Int_t sm=0;sm<17;sm++){
fLineSMid[sm] = new TLine( sm+1, 0., sm+1, 91.);
}
for (Int_t sm=0;sm<22;sm++){
fLineLTMbitId[sm] = new TLine( 0., sm+1, 72. ,sm+1);
}
-
+ */
}
//____________________________________________________________________________
//
SetName((const char*)qadm.GetName()) ;
SetTitle((const char*)qadm.GetTitle());
-
+ /*
for (Int_t sm=0;sm<17;sm++){
fLineSMid[sm]=qadm.fLineSMid[sm];
}
for (Int_t sm=0;sm<22;sm++){
fLineLTMbitId[sm] = qadm.fLineLTMbitId[sm];
}
+ */
}
//__________________________________________________________________
delete fLineExpTotMin;
if (fLineExpTotMax)
delete fLineExpTotMax;
- for (Int_t sm=0;sm<17;sm++){
+ /*
+ for (Int_t sm=0;sm<17;sm++){
if (fLineSMid[sm])
delete fLineSMid[sm];
}
if (fLineLTMbitId[sm])
delete fLineLTMbitId[sm];
}
+ */
}
//----------------------------------------------------------------------------
AliTOFChannelOnlineStatusArray* AliTOFQADataMakerRec::GetCalibData()
TH2F * h25 = new TH2F("hTOFhitsCTTM","Map of hit pads according to CTTM numbering;LTM index;bit index", 72, 0., 72., 23, 0., 23.);
TH2F * h26 = new TH2F("hTOFmacropadCTTM","Map of hit macropads according to CTTM numbering;LTM index; bit index", 72, 0., 72., 23, 0., 23.);
TH2F * h27 = new TH2F("hTOFmacropadDeltaPhiTime","#Deltat vs #Delta#Phi of hit macropads;#Delta#Phi (degrees);#DeltaBX", 18, 0., 180., 20, 0., 20.0);
-
+ TH2I *h28 = new TH2I("hBXVsCttmBit","BX ID in TOF matching window vs trg channel; trg channel; BX", 1728, 0, 1728, 24, 0, 24);
+ TH2F *h29 = new TH2F("hTimeVsCttmBit","TOF raw time vs trg channel; trg channel; raw time (ns)", 1728, 0., 1728., fgNbinsTime, fgRangeMinTime, fgRangeMaxTime);
+
h25->GetYaxis()->SetTickLength(-0.02);
h26->GetYaxis()->SetTickLength(-0.02);
h25->GetYaxis()->SetNdivisions(210);
h25->Sumw2() ;
h26->Sumw2() ;
h27->Sumw2() ;
+ h28->Sumw2() ;
+ h29->Sumw2() ;
+
//add lines for DQM shifter
fLineExpTimeMin = new TLine(200., 0., 200., 0.);
- fLineExpTimeMax = new TLine(250., 0., 250., 0.);
- fLineExpTotMin = new TLine(5., 0., 5., 0.);
- fLineExpTotMax = new TLine(20., 0., 20., 0.);
+ fLineExpTimeMax = new TLine(300., 0., 300., 0.);
+ fLineExpTotMin = new TLine(10., 0., 10., 0.);
+ fLineExpTotMax = new TLine(15., 0., 15., 0.);
fLineExpTimeMin->SetLineColor(kGreen);
fLineExpTimeMin->SetLineWidth(2);
fLineExpTotMax->SetLineColor(kGreen);
fLineExpTotMax->SetLineWidth(2);
- for (Int_t sm=0;sm<17;sm++){
+ /*
+ for (Int_t sm=0;sm<17;sm++){
fLineSMid[sm]->SetLineColor(kMagenta);
fLineSMid[sm]->SetLineWidth(2);
}
-
+ */
h5->GetListOfFunctions()->Add(fLineExpTimeMin);
h5->GetListOfFunctions()->Add(fLineExpTimeMax);
h10->GetListOfFunctions()->Add(fLineExpTotMin);
h10->GetListOfFunctions()->Add(fLineExpTotMax);
-
+ /*
for (Int_t sm=0;sm<17;sm++){
h16->GetListOfFunctions()->Add(fLineSMid[sm]);
h17->GetListOfFunctions()->Add(fLineSMid[sm]);
}
-
for (Int_t sm=0;sm<71;sm++){
fLineLTMid[sm]->SetLineColor(kBlack);
fLineLTMid[sm]->SetLineWidth(1);
h26->GetListOfFunctions()->Add(fLineLTMbitId[sm]);
h25->GetListOfFunctions()->Add(fLineLTMbitId[sm]);
}
-
+ */
TPaveText *phosHoleBox=new TPaveText(13,38,16,53,"b");
phosHoleBox->SetFillStyle(0);
phosHoleBox->SetFillColor(kWhite);
// h10->SetDrawOption("logy");
Add2RawsList(h0, 0, !expert, image, !saveCorr) ;
- Add2RawsList(h1, 1, expert, !image, !saveCorr) ;
- Add2RawsList(h2, 2, expert, !image, !saveCorr) ;
- Add2RawsList(h3, 3, expert, !image, !saveCorr) ;
- Add2RawsList(h4, 4, expert, !image, !saveCorr) ;
+ Add2RawsList(h1, 1, expert, !image, !saveCorr) ;
+ Add2RawsList(h2, 2, expert, !image, !saveCorr) ;
+ Add2RawsList(h3, 3, expert, !image, !saveCorr) ;
+ Add2RawsList(h4, 4, expert, !image, !saveCorr) ;
Add2RawsList(h5, 5, !expert, image, !saveCorr) ;
- Add2RawsList(h6, 6, expert, !image, !saveCorr) ;
- Add2RawsList(h7, 7, expert, !image, !saveCorr) ;
- Add2RawsList(h8, 8, expert, !image, !saveCorr) ;
- Add2RawsList(h9, 9, expert, !image, !saveCorr) ;
+ Add2RawsList(h6, 6, expert, !image, !saveCorr) ;
+ Add2RawsList(h7, 7, expert, !image, !saveCorr) ;
+ Add2RawsList(h8, 8, expert, !image, !saveCorr) ;
+ Add2RawsList(h9, 9, expert, !image, !saveCorr) ;
Add2RawsList(h10, 10, !expert, image, !saveCorr) ;
Add2RawsList(h11, 11, expert, !image, !saveCorr) ;
Add2RawsList(h12, 12, expert, !image, !saveCorr) ;
Add2RawsList(h13, 13, expert, !image, !saveCorr) ;
Add2RawsList(h14, 14, expert, !image, !saveCorr) ;
Add2RawsList(h15, 15, expert, !image, !saveCorr) ;
- Add2RawsList(h16, 16, !expert, image, !saveCorr) ;
- Add2RawsList(h17, 17, !expert, image, !saveCorr) ;
- Add2RawsList(h18, 18, expert, !image, !saveCorr) ;
- Add2RawsList(h19, 19, expert, !image, !saveCorr) ;
- Add2RawsList(h20, 20, expert, !image, !saveCorr) ;
- Add2RawsList(h21, 21, expert, !image, !saveCorr) ;
- Add2RawsList(h22, 22, !expert, image, !saveCorr) ;
- Add2RawsList(h23, 23, !expert, !image, !saveCorr) ;
- Add2RawsList(h24, 24, !expert, !image, !saveCorr) ;
- Add2RawsList(h25, 25, !expert, !image, !saveCorr) ;
- Add2RawsList(h26, 26, !expert, image, !saveCorr) ;
- Add2RawsList(h27, 27, !expert, image, !saveCorr) ;
+ Add2RawsList(h16, 16, !expert, image, !saveCorr) ;
+ Add2RawsList(h17, 17, !expert, image, !saveCorr) ;
+ Add2RawsList(h18, 18, expert, !image, !saveCorr) ;
+ Add2RawsList(h19, 19, expert, !image, !saveCorr) ;
+ Add2RawsList(h20, 20, !expert, image, !saveCorr) ;
+ Add2RawsList(h21, 21, !expert, image, !saveCorr) ;
+ Add2RawsList(h22, 22, expert, !image, !saveCorr) ;
+ Add2RawsList(h23, 23, expert, !image, !saveCorr) ;
+ Add2RawsList(h24, 24, expert, !image, !saveCorr) ;
+ Add2RawsList(h25, 25, expert, !image, !saveCorr) ;
+ Add2RawsList(h26, 26, expert, image, !saveCorr) ;
+ Add2RawsList(h27, 27, expert, !image, !saveCorr) ;
+ Add2RawsList(h28, 28, expert, !image, !saveCorr) ;
+ Add2RawsList(h29, 29, expert, !image, !saveCorr) ;
+
//
ClonePerTrigClass(AliQAv1::kRAWS); // this should be the last line
}
//
// makes data from Raws
//
+ // AliLog::SetClassDebugLevel("AliRawReader",0);
+ // AliLog::SetClassDebugLevel("AliTOFRawStream",0);
+ // AliLog::SetClassDebugLevel("AliTOFDecoderV2",0);
+ AliLog::SetGlobalLogLevel(AliLog::kError);
if (rawReader->GetType()==7) {
Double_t tdc2ns=AliTOFGeometry::TdcBinWidth()*1E-3;//in ns
clonesRawData = (TClonesArray*)fTOFRawStream.GetRawData();
for (Int_t iRawData = 0; iRawData<clonesRawData->GetEntriesFast(); iRawData++) {
AliTOFrawData *tofRawDatum = (AliTOFrawData*)clonesRawData->UncheckedAt(iRawData);
+ Float_t tofRawTime=tofRawDatum->GetTOF()*tdc2ns;
if (tofRawDatum->GetTOF()){
equipmentID[0]=iDDL;
if (equipmentID[2]==1) //crate left, A-side or C-side
FillRawsData(15,equipmentID[0]);
else
- FillRawsData(15,equipmentID[0]-1);
+ FillRawsData(15,equipmentID[0]-1);
+
+ //retrieve CTTM index (Ltm, bit)
+ GetCTTMIndex(equipmentID, indexCTTM);
+
+ //get BX index within TOF-matching window
+ Int_t indexBC=-1;
+ indexBC= TMath::Nint(tofRawTime/24.4);
+
+ Int_t indexCttmChannel=indexCTTM[0]*24+indexCTTM[1];
+ FillRawsData(28,indexCttmChannel,indexBC);
+ FillRawsData(29,indexCttmChannel,tofRawTime);
//fired macropad map (from LTM hits) - only for low multi evts (UPC)
if ((nFiredMacropad<=fgCutNmaxFiredMacropad)){
iFiredMacropad++;
- AliInfo(Form("Event found with %i fired macropads in BCID = %i!",nFiredMacropad,BCID));
-
- GetCTTMIndex(equipmentID, indexCTTM);
+ //AliInfo(Form("Event found with %i fired macropads in BCID = %i!",nFiredMacropad,BCID));
FillRawsData(26,indexCTTM[0],indexCTTM[1]);
Float_t halfSMphi=-999.0;
- Int_t indexBC=-1;
if (indexCTTM[0]<36)
halfSMphi=indexCTTM[0]*10.+5.;
else halfSMphi=(indexCTTM[0]-36)*10.+5.;
macropadPhiTimeUPC[iFiredMacropad][0]=halfSMphi;
- indexBC= TMath::Nint(tofRawDatum->GetTOF()*tdc2ns)/25;
macropadPhiTimeUPC[iFiredMacropad][1]=indexBC;
}
}
ntof[0]++; //counter for tof hits
//fill global spectra for DQM plots
- FillRawsData(5, tofRawDatum->GetTOF()*tdc2ns) ;//in ns
+ FillRawsData(5, tofRawTime) ;//in ns
FillRawsData(10, tofRawDatum->GetTOT()*tot2ns) ;//in ns
- FillRawsData(23, BCID, tofRawDatum->GetTOF()*tdc2ns) ;//in ns
+ FillRawsData(23, BCID, tofRawTime) ;//in ns
//fill side-related spectra for experts plots
Int_t ddlACside=iDDL/36; // 0 or 1
if (volumeID2[0]>4 && volumeID2[0]<14){ //O side
if (ddlPerSm<2){ //A side
ntof[1]++;
- FillRawsData(6, tofRawDatum->GetTOF()*tdc2ns) ;
+ FillRawsData(6, tofRawTime) ;
FillRawsData(11, tofRawDatum->GetTOT()*tot2ns) ;
} else { //C side
ntof[3]++;
- FillRawsData(8, tofRawDatum->GetTOF()*tdc2ns) ;
+ FillRawsData(8, tofRawTime) ;
FillRawsData(13, tofRawDatum->GetTOT()*tot2ns) ;
}
} else {
if (volumeID2[0]<5 || volumeID2[0]>13){ //I side
if (ddlPerSm<2){ //A side
ntof[2]++;
- FillRawsData(7, tofRawDatum->GetTOF()*tdc2ns) ;
+ FillRawsData(7, tofRawTime) ;
FillRawsData(12, tofRawDatum->GetTOT()*tot2ns) ;
} else {//C side
ntof[4]++;
- FillRawsData(9, tofRawDatum->GetTOF()*tdc2ns) ;
+ FillRawsData(9, tofRawTime) ;
FillRawsData(14, tofRawDatum->GetTOT()*tot2ns) ;
}
}
//compute TRM offset
Int_t trm= iDDL*10+(equipmentID[1]-3);
- FillRawsData(20+ddlACside,trm,tofRawDatum->GetTOF()*tdc2ns);
- FillRawsData(22,GetStripIndex(volumeID),tofRawDatum->GetTOF()*tdc2ns) ;
+ FillRawsData(20+ddlACside,trm,tofRawTime);
+ FillRawsData(22,GetStripIndex(volumeID),tofRawTime) ;
Short_t fea = volumeID2[4]/12;
Float_t hitmapx = volumeID2[0] + ((Double_t)(3 - fea) + 0.5) *0.25;
FillRawsData(17,hitmapx,GetStripIndex(volumeID2));
else { //orphans
if (!(fCalibData->GetNoiseStatus(chIndex) == AliTOFChannelOnlineStatusArray::kTOFNoiseBad)
&& (fCalibData->GetHWStatus(chIndex) == AliTOFChannelOnlineStatusArray::kTOFHWOk))
- FillRawsData(19, tofRawDatum->GetTOF()*tdc2ns) ;//in ns
+ FillRawsData(19, tofRawTime) ;//in ns
}//end orphans
}//end volumeID check
}//end equipID check
}
}//end cut on number of fired macropad
} else {
- AliDebug(1,Form("Event of type %d found. Skipping non-physics event for QA.\n", rawReader->GetType()));
+ AliDebug(2,Form("Event of type %d found. Skipping non-physics event for QA.\n", rawReader->GetType()));
}
//fill reference map for DQM shifter only once in a detector cycle
if ( (htmp=(TH1*)arrRW[22]) ) htmp->SetOption("colz");
if ( (htmp=(TH1*)arrRW[23]) ) htmp->SetOption("colz");
if ( (htmp=(TH1*)arrRW[24]) ) htmp->SetOption("colz");
+ if ( (htmp=(TH1*)arrRW[28]) ) htmp->SetOption("colz");
+ if ( (htmp=(TH1*)arrRW[29]) ) htmp->SetOption("colz");
+
if ( (htmp=(TH1*)arrRW[25]) ) {
htmp->SetOption("colz");
htmp->SetTitle(title25.Data());
stripOffset = nStripC+nStripB+nStripA+nStripB;
break;
default:
- AliDebug(1,Form("Wrong plate number in TOF (%d) !",iplate));
+ // AliDebug(2,Form("Wrong plate number in TOF (%d) !",iplate));
break;
};
Int_t zindex=npadZ*(istrip+stripOffset)+(ipadZ+1);
stripOffset = nStripC+nStripB+nStripA+nStripB;
break;
default:
- AliDebug(1,Form("Wrong plate number in TOF (%d) !",iplate));
+ //AliDebug(2,Form("Wrong plate number in TOF (%d) !",iplate));
stripOffset=-1;
break;
};
//
for (Int_t j=0;j<5;j++){
if (volumeID[j]<0) {
- AliDebug(1,Form("Invalid detector volume index for volumeID[%i]",j));
+ //AliDebug(2,Form("Invalid detector volume index for volumeID[%i]",j));
return kFALSE;
}
}
//Checks equipment ID validity
for (Int_t j=0;j<5;j++){
if (equipmentID[j]<0) {
- AliDebug(1,Form("Invalid equipment volume index for equipmentID[%i]",j));
+ // AliDebug(2,Form("Invalid equipment volume index for equipmentID[%i]",j));
return kFALSE;
}
}
return;
}
- FILE * configFile = fopen(filename,"r");
- if (!configFile){
+ std::fstream configFile;
+ configFile.open(filename, std::fstream::in);
+ if (!configFile.is_open()){
AliInfo("Cannot open config file with histograms ranges -> use default values.");
SetDefaultHistogramRange();
return;
}
- if (feof(configFile)){
- AliInfo("Unexpected EOF of config file with histograms ranges -> use default values.");
+ //check file size
+ Int_t begin = configFile.tellg();
+ configFile.seekg(0, std::fstream::end); /* end */
+ Int_t end = configFile.tellg();
+ Int_t size = end - begin;
+ configFile.seekg(0, std::fstream::beg); /* rewind file */
+ if (size <= 0){
+ AliInfo(Form("Unexpected EOF of config file with histograms ranges. File size: %d -> use default values", size));
SetDefaultHistogramRange();
return;
}
Int_t nbinsMulti=0,nbinsTime=0;
Float_t minTime=9999.0, maxTime=-9999.0;
Int_t cutFiredMacropad=0;
+ TString endoflist;
+ while (!configFile.eof()) {
+ configFile >> cutFiredMacropad >> minMulti >> maxMulti >> minTime >> maxTime;
+ configFile >> endoflist;
+ if (endoflist.Contains("end")) break;
+ }
- fscanf(configFile,"%10i %10i %10i %10f %10f", &cutFiredMacropad,&minMulti,&maxMulti,&minTime,&maxTime);
-
//set multiplicity histo ranges
if (minMulti>maxMulti){
- AliInfo("Invalid range for multiplicity histogram set. Changing to defualt values.");
+ AliInfo("Invalid range for multiplicity histogram set. Changing to default values.");
SetDefaultMultiHistogramRange();
} else {
nbinsMulti = maxMulti-minMulti;
SetNbinsMultiplicityHisto(nbinsMulti);
SetMultiplicityHistoRange(minMulti,maxMulti);
- AliInfo(Form("Setting multiplicity histogram ranges to: multMin = %i - multMax = %i - nMultBins = %i",
- fgRangeMinMultiplicity, fgRangeMaxMultiplicity, fgNbinsMultiplicity));
+ //AliInfo(Form("Setting multiplicity histogram ranges to: multMin = %i - multMax = %i - nMultBins = %i", fgRangeMinMultiplicity, fgRangeMaxMultiplicity, fgNbinsMultiplicity));
}
//set time histo ranges
if (minTime>maxTime){
- AliInfo("Invalid range for time histogram set. Changing to defualt values.");
+ AliInfo("Invalid range for time histogram set. Changing to default values.");
SetDefaultTimeHistogramRange();
} else {
nbinsTime = TMath::Nint((maxTime - minTime)/fgkNbinsWidthTime);//ns
maxTime=minTime+nbinsTime*fgkNbinsWidthTime;//ns
SetNbinsTimeHisto(nbinsTime);
SetTimeHistoRange(minTime,maxTime);
- AliInfo(Form("Setting time histogram ranges to: timeMin = %5.2f ns - timeMax = %5.2f ns - nTimeBins = %i",
- fgRangeMinTime, fgRangeMaxTime,fgNbinsTime));
+ //AliInfo(Form("Setting time histogram ranges to: timeMin = %5.2f ns - timeMax = %5.2f ns - nTimeBins = %i", fgRangeMinTime, fgRangeMaxTime,fgNbinsTime));
}
if ((cutFiredMacropad>0)&&(cutFiredMacropad<fgkFiredMacropadLimit)){
SetDefaultCutNmaxFiredMacropad();
} else {
SetCutNmaxFiredMacropad(cutFiredMacropad);
- AliInfo(Form("Setting cut on fired macropad to: = %i",cutFiredMacropad));
+ //AliInfo(Form("Setting cut on fired macropad to: = %i",cutFiredMacropad));
}
-
- fclose(configFile);
+ AliInfo(Form("Setting: multMin = %i - multMax = %i - nMultBins = %i, timeMin = %5.2f ns - timeMax = %5.2f ns - nTimeBins = %i, cutMaxFiredMacropad = %i",
+ fgRangeMinMultiplicity, fgRangeMaxMultiplicity, fgNbinsMultiplicity, fgRangeMinTime, fgRangeMaxTime,fgNbinsTime, cutFiredMacropad));
+ configFile.close();
return;
}
//
// set default histogram ranges (tuned on 2011 pp collisions)
//
- AliInfo("Setting all histogram ranges to default values.");
+ //AliInfo("Setting all histogram ranges to default values.");
SetDefaultMultiHistogramRange();
SetDefaultTimeHistogramRange();
SetDefaultCutNmaxFiredMacropad();
//
SetMultiplicityHistoRange (0, 200);
SetNbinsMultiplicityHisto(200);
- AliInfo("Setting Multiplicity histogram ranges to default values.");
- AliInfo(Form("multMin = %i - multMax = %i - nMultBins = %i",
- fgRangeMinMultiplicity, fgRangeMaxMultiplicity, fgNbinsMultiplicity));
+ //AliInfo("Setting Multiplicity histogram ranges to default values.");
+ //AliInfo(Form("multMin = %i - multMax = %i - nMultBins = %i",
+ // fgRangeMinMultiplicity, fgRangeMaxMultiplicity, fgNbinsMultiplicity));
return;
}
SetNbinsTimeHisto(250);
SetTimeHistoRange (0.0,610.);
- AliInfo("Setting Time histogram ranges to default values:");
- AliInfo(Form("timeMin = %5.2f ns - timeMax = %5.2f ns - nTimeBins = %i",
- fgRangeMinTime, fgRangeMaxTime,fgNbinsTime));
+ // AliInfo("Setting Time histogram ranges to default values:");
+ // AliInfo(Form("timeMin = %5.2f ns - timeMax = %5.2f ns - nTimeBins = %i",
+ // fgRangeMinTime, fgRangeMaxTime,fgNbinsTime));
return;
}
//
// set default cut on fired macropad
//
- SetCutNmaxFiredMacropad(5);
- AliInfo(Form("Setting cut on fired macropad to default values: NfiredMacropad = %i", fgCutNmaxFiredMacropad));
+ SetCutNmaxFiredMacropad(50);
+ // AliInfo(Form("Setting cut on fired macropad to default values: NfiredMacropad = %i", fgCutNmaxFiredMacropad));
return;
}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff