1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
18 #include "AliMUONTrackerQAChecker.h"
20 /// \class AliMUONTrackerQAChecker
22 /// Implementation of AliQACheckerBase for MCH and MTR
24 /// For the moment we only implement the checking of raw data QA for the tracker
25 /// by looking at the occupancy at the bus patch level.
27 /// \author Laurent Aphecetche, Subatech
29 #include "AliCDBManager.h"
30 #include "AliCodeTimer.h"
32 #include "AliMUONQAIndices.h"
33 #include "AliMUONRecoParam.h"
34 #include "AliMpBusPatch.h"
35 #include "AliMpDDLStore.h"
37 #include "Riostream.h"
39 #include "TDirectory.h"
44 #include "TPaveText.h"
45 #include "TVirtualPad.h"
50 ClassImp(AliMUONTrackerQAChecker)
55 //___________________________________________________________________________
65 // Calculates the trimmed (tmean) mean
66 // of a sample (x) and estimates the variance (tvar)
70 // First check input parameters
72 // number of observations
78 if ( alpha < 0 || alpha >= 0.5 )
79 // proportion of observations
80 // to be trimmed at each end of the sorted sample
85 // Input parameters are good. Let's move on.
87 // Insure we use a sample sorted into ascending order.
89 Int_t* indices = new Int_t[n];
91 TMath::Sort(n,x,indices,kFALSE);
93 Double_t* sx = new Double_t[n];
95 for ( Int_t i = 0; i < n; ++i )
97 sx[i] = x[indices[i]];
102 // Number of observations trimmed at each end.
104 Int_t k = TMath::FloorNint(alpha * n);
108 for ( Int_t i = k; i < n - k ; ++i )
113 tmean = sum / ( n - 2 * k );
117 for ( Int_t i = k; i < n - k; ++i )
119 t2 += (sx[i] - tmean) * (sx[i] - tmean);
124 k * (sx[k] - tmean) * (sx[k] - tmean) +
125 k * (sx[n - k - 1] - tmean) * (sx[n - k - 1] - tmean)
128 // get the min and max for the non-rejected values
132 for ( Int_t i = k; i < n-k; ++i )
134 min = TMath::Min(min,sx[i]);
135 max = TMath::Max(max,sx[i]);
143 //___________________________________________________________________________
144 Int_t GetColorFromCheckCode(AliMUONVQAChecker::ECheckCode code)
146 if ( code == AliMUONVQAChecker::kInfo ) return AliMUONVQAChecker::kInfoColor;
147 else if ( code == AliMUONVQAChecker::kWarning ) return AliMUONVQAChecker::kWarningColor;
148 else if ( code == AliMUONVQAChecker::kFatal) return AliMUONVQAChecker::kFatalColor;
149 else return AliMUONVQAChecker::kErrorColor;
152 const char* NOTENOUGHEVENTMESSAGE = "Not enough event to judge. Please wait a bit";
153 const char* NOTIFYEXPERTMESSAGE = "PLEASE NOTIFY EXPERT !";
154 const char* ALLISFINEMESSAGE = "All is fine. Just enjoy.";
156 const int NOTENOUGHEVENTLIMIT = 50;
158 //___________________________________________________________________________
159 void SetupHisto(Int_t neventsseen, Int_t neventsused, const TObjArray& messages, TH1& histo, AliMUONVQAChecker::ECheckCode code, const char* extraopt="")
161 Bool_t allIsFine(kFALSE);
163 if ( code == AliMUONVQAChecker::kInfo )
168 Double_t y1 = 0.99 - (messages.GetLast()+(allIsFine?1:0)+2)*0.075;
170 TPaveText* text = new TPaveText(0.5,y1,0.99,0.99,"NDC");
172 text->AddText(Form("MCH RUN %d - %d events seen - %d events used",AliCDBManager::Instance()->GetRun(),neventsseen,neventsused));
174 TIter next(&messages);
177 while ( ( str = static_cast<TObjString*>(next()) ) )
179 text->AddText(str->String());
184 text->AddText(ALLISFINEMESSAGE);
187 text->SetFillColor(GetColorFromCheckCode(code));
189 Int_t color = GetColorFromCheckCode(code);
191 histo.SetFillStyle(1001);
192 histo.SetFillColor(color);
195 histo.SetOption(sopt.Data());
197 histo.SetBit(TH1::kNoTitle);
199 TList* lstF = histo.GetListOfFunctions();
200 TObject* title = lstF->FindObject("title");
201 if (title) delete lstF->Remove(title);
207 //___________________________________________________________________________
208 void ShowTrueValue(TH1& hrostatusnorm, Int_t v)
210 // For a bar > 100% we put the text inside the bar (as TEXT45 option is putting above which
211 // would go offscreen)
213 Int_t bin = hrostatusnorm.FindBin(1.0*v);
215 Double_t value = hrostatusnorm.GetBinContent(bin);
219 TLatex* l = new TLatex(hrostatusnorm.GetBinCenter(bin),50,Form("%g",value));
220 l->SetTextFont(gStyle->GetTextFont());
223 l->SetTextSize(0.02*hrostatusnorm.GetMarkerSize());
224 hrostatusnorm.GetListOfFunctions()->Add(l);
229 //__________________________________________________________________
230 AliMUONTrackerQAChecker::AliMUONTrackerQAChecker() : AliMUONVQAChecker()
235 //__________________________________________________________________
236 AliMUONTrackerQAChecker::~AliMUONTrackerQAChecker()
241 //__________________________________________________________________
242 AliMUONTrackerQAChecker::AliMUONTrackerQAChecker(const AliMUONTrackerQAChecker& qac) :
243 AliMUONVQAChecker(qac)
248 //______________________________________________________________________________
249 AliMUONVQAChecker::ECheckCode*
250 AliMUONTrackerQAChecker::CheckRecPoints(TObjArray ** list, const AliMUONRecoParam* /*recoParam*/)
253 /// Very binary check for the moment.
255 AliCodeTimerAuto("",0);
257 AliMUONVQAChecker::ECheckCode * rv = new AliMUONVQAChecker::ECheckCode[AliRecoParam::kNSpecies] ;
258 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
259 rv[specie] = AliMUONVQAChecker::kInfo;
262 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
264 TH1* h = AliQAv1::GetData(list,AliMUONQAIndices::kTrackerNumberOfClustersPerDE,AliRecoParam::ConvertIndex(specie));
266 if ( !h ) rv[specie] = AliMUONVQAChecker::kWarning; // only a warning if histo not found, in order not to kill anything because QA failed...
268 else if ( h->GetMean() == 0.0 ) rv[specie] = MarkHisto(*h,AliMUONVQAChecker::kFatal);
273 //______________________________________________________________________________
274 AliMUONVQAChecker::ECheckCode
275 AliMUONTrackerQAChecker::MarkHisto(TH1& histo, AliMUONVQAChecker::ECheckCode value) const
277 /// Mark histo as originator of some QA error/warning
279 if ( value != AliMUONVQAChecker::kInfo )
281 histo.SetBit(AliQAv1::GetQABit());
287 //______________________________________________________________________________
288 AliMUONVQAChecker::ECheckCode*
289 AliMUONTrackerQAChecker::CheckESD(TObjArray ** list, const AliMUONRecoParam* /*recoParam*/)
293 AliCodeTimerAuto("",0);
295 AliMUONVQAChecker::ECheckCode * rv = new AliMUONVQAChecker::ECheckCode[AliRecoParam::kNSpecies] ;
296 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
297 rv[specie] = AliMUONVQAChecker::kInfo;
299 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
301 TH1* h = AliQAv1::GetData(list,AliMUONQAIndices::kESDnTracks,AliRecoParam::ConvertIndex(specie));
303 if (!h) rv[specie] = AliMUONVQAChecker::kWarning;
305 else if ( h->GetMean() == 0.0 ) rv[specie] = MarkHisto(*h,AliMUONVQAChecker::kFatal); // no track -> fatal
307 h = AliQAv1::GetData(list,AliMUONQAIndices::kESDMatchTrig,AliRecoParam::ConvertIndex(specie));
309 if (!h) rv[specie] = AliMUONVQAChecker::kWarning;
311 else if (h->GetMean() == 0.0 ) rv[specie] = MarkHisto(*h,AliMUONVQAChecker::kError); // no trigger matching -> error
316 //______________________________________________________________________________
317 AliMUONVQAChecker::ECheckCode*
318 AliMUONTrackerQAChecker::CheckRaws(TObjArray ** list, const AliMUONRecoParam* recoParam)
322 AliCodeTimerAuto("",0);
324 if (!recoParam) return 0x0;
326 AliMUONVQAChecker::ECheckCode * rv = new AliMUONVQAChecker::ECheckCode[AliRecoParam::kNSpecies] ;
328 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
330 rv[specie] = AliMUONVQAChecker::kInfo;
333 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
335 TH1* hneventsseen = AliQAv1::GetData(list,AliMUONQAIndices::kTrackerNofPhysicsEventsSeen,AliRecoParam::ConvertIndex(specie));
336 TH1* hneventsused = AliQAv1::GetData(list,AliMUONQAIndices::kTrackerNofGoodPhysicsEventsUsed,AliRecoParam::ConvertIndex(specie));
338 if (!hneventsseen || !hneventsused ) continue;
340 Int_t neventsseen = TMath::Nint(hneventsseen->GetBinContent(1));
341 Int_t neventsused = TMath::Nint(hneventsused->GetBinContent(1));
343 AliMUONVQAChecker::ECheckCode c1 = AliMUONVQAChecker::kInfo;
344 AliMUONVQAChecker::ECheckCode c2 = AliMUONVQAChecker::kInfo;
345 AliMUONVQAChecker::ECheckCode c3 = AliMUONVQAChecker::kInfo;
347 TH1* hbp = AliQAv1::GetData(list,AliMUONQAIndices::kTrackerBusPatchOccupancy,AliRecoParam::ConvertIndex(specie));
348 TH1* hbpconfig = AliQAv1::GetData(list,AliMUONQAIndices::kTrackerBusPatchConfig,AliRecoParam::ConvertIndex(specie));
349 TH1* hddl = AliQAv1::GetData(list,AliMUONQAIndices::kTrackerDDLOccupancy,AliRecoParam::ConvertIndex(specie));
351 if ( hbp && hbpconfig && hddl )
353 c1 = BeautifyOccupancyHistograms(*hddl,*hbp,hbpconfig,neventsseen,neventsused,*recoParam);
357 AliError(Form("Could not BeautifyOccupancyHistograms : hddl=%p hbpconfig=%p hbp=%p",hddl,hbpconfig,hbp));
360 TH1* hbuspatchtokenerrors = AliQAv1::GetData(list,AliMUONQAIndices::kTrackerBusPatchTokenLostErrors,AliRecoParam::ConvertIndex(specie));
361 TH1* hrostatus = AliQAv1::GetData(list,AliMUONQAIndices::kTrackerReadoutStatus,AliRecoParam::ConvertIndex(specie));
362 TH1* hrostatusnorm = AliQAv1::GetData(list,AliMUONQAIndices::kTrackerReadoutStatusPerEvent,AliRecoParam::ConvertIndex(specie));
364 if ( hrostatus && hrostatusnorm && hbuspatchtokenerrors )
366 c2 = BeautifyReadoutHistograms(*hrostatus,*hrostatusnorm,*hbuspatchtokenerrors,
367 neventsseen,neventsused,*recoParam);
371 AliError(Form("Could not BeautifyReadoutHistograms : hrostatus=%p hrostatusnorm=%p hbuspatchtokenerrors=%p",
372 hrostatus,hrostatusnorm,hbuspatchtokenerrors));
376 TH1* heventsize = AliQAv1::GetData(list,AliMUONQAIndices::kTrackerDDLEventSize,AliRecoParam::ConvertIndex(specie));
377 TH1* heventsizeperevent = AliQAv1::GetData(list,AliMUONQAIndices::kTrackerDDLEventSizePerEvent,AliRecoParam::ConvertIndex(specie));
379 if ( heventsize && heventsizeperevent )
381 c3 = BeautifyEventsizeHistograms(*heventsize,*heventsizeperevent,
382 neventsseen,neventsused,*recoParam);
386 AliError(Form("Could not BeautifyEventsizeHistograms heventsize=%p heventsizeperevent=%p",heventsize,heventsizeperevent));
389 if ( c1 == AliMUONVQAChecker::kFatal || c2 == AliMUONVQAChecker::kFatal || c3 == AliMUONVQAChecker::kFatal )
391 rv[specie] = AliMUONVQAChecker::kFatal;
393 else if ( c1 == AliMUONVQAChecker::kError || c2 == AliMUONVQAChecker::kError || c3 == AliMUONVQAChecker::kError )
395 rv[specie] = AliMUONVQAChecker::kError;
397 else if ( c1 == AliMUONVQAChecker::kWarning || c2 == AliMUONVQAChecker::kWarning || c3 == AliMUONVQAChecker::kWarning )
399 rv[specie] = AliMUONVQAChecker::kWarning;
403 rv[specie] = AliMUONVQAChecker::kInfo;
410 //____________________________________________________________________________
411 AliMUONVQAChecker::ECheckCode
412 AliMUONTrackerQAChecker::BeautifyOccupancyHistograms(TH1& hddl,
414 const TH1* hbuspatchconfig,
417 const AliMUONRecoParam& recoParam)
419 /// Put labels, limits and so on on the TrackerBusPatchOccupancy histograms
420 /// hbuspatchconfig and hbp must have the same bin definitions
422 if ( hbuspatchconfig )
424 if ( hbp.GetNbinsX() != hbuspatchconfig->GetNbinsX() ||
425 hbp.GetXaxis()->GetXmin() != hbuspatchconfig->GetXaxis()->GetXmin() ||
426 hbp.GetXaxis()->GetXmax() != hbuspatchconfig->GetXaxis()->GetXmax() )
428 AliError("hbp and hbuspatchconfig histograms are not compatible !");
429 return AliMUONVQAChecker::kFatal;
433 hddl.SetStats(kFALSE);
434 hbp.SetXTitle("Absolute Bus Patch Id");
435 hbp.SetYTitle("Occupancy (percent)");
436 hbp.SetStats(kFALSE);
438 Double_t xmin = hbp.GetXaxis()->GetXmin();
439 Double_t xmax = hbp.GetXaxis()->GetXmax();
441 Double_t occMax(0.1); // 0.1% y-limit for the plot
442 Double_t maxToleratedOccupancy(recoParam.BuspatchOccupancyHighLimit()*100.0);
443 Double_t minToleratedOccupancy(recoParam.BuspatchOccupancyLowLimit()*100.0);
444 TLine* line1 = new TLine(xmin,maxToleratedOccupancy,xmax,maxToleratedOccupancy);
445 line1->SetLineColor(1);
446 line1->SetLineWidth(1);
448 TLine* line2 = new TLine(xmin,minToleratedOccupancy,xmax,minToleratedOccupancy);
449 line2->SetLineColor(1);
450 line2->SetLineWidth(1);
452 TList* lstF = hbp.GetListOfFunctions();
454 TObject *stats = lstF->FindObject("stats");
457 while ((obj = lstF->First())) { while(lstF->Remove(obj)) { } delete obj; }
458 if (stats) lstF->Add(stats);
461 hbp.GetListOfFunctions()->Add(line1);
462 hbp.GetListOfFunctions()->Add(line2);
464 TIter next(AliMpDDLStore::Instance()->CreateBusPatchIterator());
465 AliMpBusPatch* bp(0x0);
467 Int_t nBusPatches(0);
468 Int_t nMissingBusPatches(0);
470 while ( ( bp = static_cast<AliMpBusPatch*>(next())) )
474 Int_t bin = hbp.FindBin(bp->GetId());
476 if ( hbp.GetBinContent(bin) <= 0.0 )
478 ++nMissingBusPatches;
486 Int_t nBusPatchesAboveLimit(0);
487 Int_t nBusPatchesBelowLimit(0);
488 Double_t alpha(0.1); // trim 10% of data
489 Double_t tmean(0.0),tvar(0.0);
490 Double_t ymin(0.0),ymax(0.0);
494 Double_t* x = new Double_t[nBusPatches];
496 while ( ( bp = static_cast<AliMpBusPatch*>(next())) )
498 Int_t bin = hbp.FindBin(bp->GetId());
499 if ( hbp.GetBinContent(bin) > 0 )
501 x[n] = hbp.GetBinContent(bin);
504 if ( hbp.GetBinContent(bin) > maxToleratedOccupancy )
506 ++nBusPatchesAboveLimit;
508 if ( hbp.GetBinContent(bin) < minToleratedOccupancy )
510 // check whether this buspatch has a reason to be absent (only valid
511 // if we got the config, otherwise we cannot do the test)
512 if ( hbuspatchconfig && hbuspatchconfig->GetBinContent(bin) > 0 )
514 // should be there, so it's an error
515 ++nBusPatchesBelowLimit;
520 // computed the truncated mean of the occupancy values, in order to get a
521 // reasonable y-range for the histogram (without giant peaks to the roof
522 // for misbehaving buspatches).
523 ok = trim(n,x,alpha,tmean,tvar,ymin,ymax);
534 ymax = TMath::Max(ymax,occMax);
537 hbp.SetMaximum(ymax*1.4);
539 AliMUONVQAChecker::ECheckCode rv(AliMUONVQAChecker::kInfo);
542 messages.SetOwner(kTRUE);
544 if ( neventsseen < NOTENOUGHEVENTLIMIT )
546 messages.Add(new TObjString(NOTENOUGHEVENTMESSAGE));
552 messages.Add(new TObjString("Could not compute truncated mean. Not enough events ?"));
554 if ( neventsused < TMath::Nint(0.1*neventsseen) )
556 messages.Add(new TObjString(Form("We've actually seen %d events, but used only %d !",neventsseen,neventsused)));
557 messages.Add(new TObjString("For a normal physics run, this is highly suspect !"));
558 messages.Add(new TObjString(NOTIFYEXPERTMESSAGE));
559 rv = AliMUONVQAChecker::kFatal;
562 else if (!nBusPatches)
564 messages.Add(new TObjString(Form("Could not get the total number of buspatches (%d). ERROR !!!",
566 messages.Add(new TObjString("Please check with expert !"));
567 rv = AliMUONVQAChecker::kFatal;
571 Float_t missingBusPatchFraction = nMissingBusPatches*100.0/nBusPatches;
572 Float_t aboveLimitFraction = nBusPatchesAboveLimit*100.0/nBusPatches;
573 Float_t belowLimitFraction = nBusPatchesBelowLimit*100.0/nBusPatches;
575 messages.Add(new TObjString(Form("%5.2f %% of missing buspatches (%d out of %d)",missingBusPatchFraction,nMissingBusPatches,nBusPatches)));
576 messages.Add(new TObjString(Form("%5.2f %% bus patches above the %5.2f %% limit",aboveLimitFraction,maxToleratedOccupancy)));
577 messages.Add(new TObjString(Form("%5.2f %% bus patches below the %e %% limit",belowLimitFraction,minToleratedOccupancy)));
578 messages.Add(new TObjString(Form("Bus patch mean occupancy (truncated at %2d %%) is %7.2f %%",(Int_t)(alpha*100),tmean)));
580 if ( aboveLimitFraction > recoParam.FractionOfBuspatchOutsideOccupancyLimit()*100.0 ||
581 belowLimitFraction > recoParam.FractionOfBuspatchOutsideOccupancyLimit()*100.0 )
583 rv = AliMUONVQAChecker::kError;
587 rv = AliMUONVQAChecker::kInfo;
592 SetupHisto(neventsseen,neventsused,messages,hbp,rv);
594 /// Make as well a version for DDL occupancy, that'll be used by the shifter
595 SetupHisto(neventsseen,neventsused,messages,hddl,rv);
598 hddl.SetStats(kFALSE);
603 //______________________________________________________________________________
604 AliMUONVQAChecker::ECheckCode
605 AliMUONTrackerQAChecker::BeautifyReadoutHistograms(TH1& hrostatus,
607 const TH1& hbuspatchtokenerrors,
610 const AliMUONRecoParam& recoParam)
612 /// Normalize and put some text on the readout error histogram
613 /// Note in particular the treatment of tokenlost errors !
615 hrostatusnorm.Reset();
617 AliMUONVQAChecker::ECheckCode rv(AliMUONVQAChecker::kInfo);
620 messages.SetOwner(kTRUE);
622 if ( neventsseen < NOTENOUGHEVENTLIMIT )
624 messages.Add(new TObjString(NOTENOUGHEVENTMESSAGE));
628 hrostatusnorm.Add(&hrostatus,100.0/neventsseen);
629 hrostatusnorm.SetOption("TEXT45"); // note : cannot use HIST option, otherwise the associated function (i.e. the tpavetext) is not drawn...
630 hrostatusnorm.SetMinimum(0.0);
631 hrostatusnorm.SetMaximum(100.0);
633 Double_t tokenLost = hrostatusnorm.GetBinContent(hrostatusnorm.FindBin(1.0*AliMUONQAIndices::kTrackerRawNofTokenLostErrors));
635 if ( tokenLost > recoParam.TokenLostLimit() )
637 rv = AliMUONVQAChecker::kError;
639 messages.Add(new TObjString("There are some token lost errors !"));
640 messages.Add(new TObjString("PLEASE CHECK THE BUSY TIME FOR MUON !"));
641 messages.Add(new TObjString("If above 5 ms please have the MUON expert"));
642 messages.Add(new TObjString("check the following bus patches :"));
644 for ( Int_t i = 1; i <= hbuspatchtokenerrors.GetNbinsX(); ++i )
646 if ( hbuspatchtokenerrors.GetBinContent(i) > 0 )
648 messages.Add(new TObjString(Form("BP %4d",i)));
653 if ( hrostatusnorm.GetBinContent(hrostatusnorm.FindBin(AliMUONQAIndices::kTrackerRawNofEmptyEvents)) > 25.0 )
655 messages.Add(new TObjString("Too many empty events !"));
656 messages.Add(new TObjString(NOTIFYEXPERTMESSAGE));
657 rv = AliMUONVQAChecker::kFatal;
661 SetupHisto(neventsseen,neventsused,messages,hrostatusnorm,rv,"text45");
663 ShowTrueValue(hrostatusnorm,AliMUONQAIndices::kTrackerRawNofTokenLostErrors);
664 ShowTrueValue(hrostatusnorm,AliMUONQAIndices::kTrackerRawNofParityErrors);
665 ShowTrueValue(hrostatusnorm,AliMUONQAIndices::kTrackerRawNofPaddingErrors);
666 ShowTrueValue(hrostatusnorm,AliMUONQAIndices::kTrackerRawNofGlitchErrors);
671 //______________________________________________________________________________
672 AliMUONVQAChecker::ECheckCode
673 AliMUONTrackerQAChecker::BeautifyEventsizeHistograms(TH1& heventsize,
674 TH1& heventsizeperevent,
677 const AliMUONRecoParam& recoParam)
679 /// Normalize and put some text on the event size histogram
681 AliMUONVQAChecker::ECheckCode rv(AliMUONVQAChecker::kInfo);
683 heventsizeperevent.Reset();
686 messages.SetOwner(kTRUE);
688 if ( neventsseen < NOTENOUGHEVENTLIMIT )
690 messages.Add(new TObjString(NOTENOUGHEVENTMESSAGE));
694 heventsizeperevent.Add(&heventsize,1.0/neventsseen/1024.0);
695 heventsizeperevent.SetMinimum(0);
697 Double_t totalEventSizePerEvent = heventsizeperevent.Integral();
702 if ( totalEventSizePerEvent > recoParam.EventSizeHardLimit() )
704 rv = AliMUONVQAChecker::kError;
705 msg = "That is really too high.";
706 action = NOTIFYEXPERTMESSAGE;
708 else if ( totalEventSizePerEvent > recoParam.EventSizeSoftLimit() )
710 msg = "That is a bit high.";
711 action = "Please keep an eye on it.";
712 rv = AliMUONVQAChecker::kWarning;
716 rv = AliMUONVQAChecker::kInfo;
719 messages.Add(new TObjString(Form("<MCH event size> %5.1f KB/event\n",totalEventSizePerEvent)));
722 messages.Add(new TObjString(msg));
723 messages.Add(new TObjString(action));
727 SetupHisto(neventsseen,neventsused,messages,heventsizeperevent,rv);