Updates to the QA framework.
authorcholm <cholm@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 29 Jun 2011 10:40:59 +0000 (10:40 +0000)
committercholm <cholm@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 29 Jun 2011 10:40:59 +0000 (10:40 +0000)
Better structure of code.
Make sure the histograms are properly displayed in the
generated images.
Added # R/O errors to the raw part.

FMD/AliFMDQAChecker.cxx
FMD/AliFMDQAChecker.h
FMD/AliFMDQADataMakerRec.cxx
FMD/AliFMDQADataMakerRec.h

index 11a698a..2da8902 100644 (file)
@@ -34,6 +34,7 @@
 #include <iostream>
 #include <TCanvas.h>
 #include <TPaveText.h>
+#include <TStyle.h>
 
 // --- AliRoot header files ---
 #include "AliLog.h"
@@ -46,6 +47,48 @@ ClassImp(AliFMDQAChecker)
 #if 0
 ; // This is for Emacs! - do not delete
 #endif
+
+//__________________________________________________________________
+Double_t 
+AliFMDQAChecker::CheckOne(AliQAv1::ALITASK_t          what,
+                         AliRecoParam::EventSpecie_t specie, 
+                         TH1*                        hist) const
+{
+  if(what == AliQAv1::kESD) return CheckESD(specie, hist);
+  if(what == AliQAv1::kRAW) return CheckRaw(specie, hist);
+  if(what == AliQAv1::kSIM) return CheckSim(specie, hist);
+  if(what == AliQAv1::kREC) return CheckRec(specie, hist);
+  return 0;
+}
+//__________________________________________________________________
+Double_t 
+AliFMDQAChecker::CheckESD(AliRecoParam::EventSpecie_t /* specie*/, 
+                         TH1*                        hist) const
+{
+  return (hist->GetMean() > 0 ? 1 : 0);
+}
+//__________________________________________________________________
+Double_t 
+AliFMDQAChecker::CheckRaw(AliRecoParam::EventSpecie_t /* specie*/, 
+                         TH1*                        hist) const
+{
+  return (hist->GetMean() > 0 ? 1 : 0);
+}
+//__________________________________________________________________
+Double_t 
+AliFMDQAChecker::CheckSim(AliRecoParam::EventSpecie_t /* specie*/, 
+                         TH1*                        hist) const
+{
+  return (hist->GetMean() > 0 ? 1 : 0);
+}
+//__________________________________________________________________
+Double_t 
+AliFMDQAChecker::CheckRec(AliRecoParam::EventSpecie_t /* specie*/, 
+                         TH1*                        hist) const
+{
+  return (hist->GetMean() > 0 ? 1 : 0);
+}
+
 //__________________________________________________________________
 void AliFMDQAChecker::Check(Double_t*                   rv, 
                            AliQAv1::ALITASK_t          what, 
@@ -73,26 +116,37 @@ void AliFMDQAChecker::Check(Double_t*                   rv,
     
     if(!list[specie]) continue;
     
-    TH1F* hist  = 0;
+    TH1* hist  = 0;
     Int_t nHist = list[specie]->GetEntriesFast();
     for(Int_t i= 0; i< nHist; i++) {
       
-      if (!(hist = static_cast<TH1F*>(list[specie]->At(i)))) continue;
+      if (!(hist = static_cast<TH1*>(list[specie]->At(i)))) continue;
       
-      if(what == AliQAv1::kESD) 
-       rv[specie] += (hist->GetMean() > 0 ? 1 : 0);
-      if(what == AliQAv1::kRAW) 
-       rv[specie] += (hist->GetMean() > 0 ? 1 : 0);
-      if(what == AliQAv1::kSIM) 
-       rv[specie] += (hist->GetMean() > 0 ? 1 : 0);
-      if(what == AliQAv1::kREC) 
-       rv[specie] += (hist->GetMean() > 0 ? 1 : 0);
+      rv[specie] += CheckOne(what, AliRecoParam::ConvertIndex(specie), hist);
     } // for (int i ...)
     // if (count != 0) rv[specie] /= count;
   }
   // return rv;
 }
 
+namespace {
+  Int_t CheckForLog(TAxis*       axis,
+                   TVirtualPad* pad, 
+                   Int_t        xyz)
+  {
+    Int_t ret = 0;
+    TString t(axis->GetTitle());
+    if (!t.Contains("[log]", TString::kIgnoreCase)) return 0;
+    t.ReplaceAll("[log]", "");
+    switch (xyz) { 
+    case 1: pad->SetLogx(); ret |= 0x1; break;
+    case 2: pad->SetLogy(); ret |= 0x2; break;
+    case 3: pad->SetLogz(); ret |= 0x4; break;
+    }
+    axis->SetTitle(t);
+    return ret;
+  }
+}
 
 //____________________________________________________________________________ 
 void 
@@ -112,27 +166,36 @@ AliFMDQAChecker::MakeImage(TObjArray** list,
   Double_t max     = 0;
   Double_t min     = 10000;
   for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
-    if (! AliQAv1::Instance(AliQAv1::GetDetIndex(GetName()))
-       ->IsEventSpecieSet(AliRecoParam::ConvertIndex(specie))) 
+    AliRecoParam::EventSpecie_t spe = AliRecoParam::ConvertIndex(specie);
+    if (!AliQAv1::Instance(AliQAv1::GetDetIndex(GetName()))
+       ->IsEventSpecieSet(spe)) 
       continue;
     // if (!AliQAv1::Instance()->IsEventSpecieSet(specie)) continue ;
     
-    if(!list[specie]) continue;
+    if(!list[specie] || list[specie]->GetEntriesFast() == 0) continue;
 
-    TH1F* hist  = 0;
+    TH1* hist  = 0;
     Int_t nHist = list[specie]->GetEntriesFast();
     for(Int_t i= 0; i< nHist; i++) {
       hist = static_cast<TH1F*>(list[specie]->At(i));
       if (hist && hist->TestBit(AliQAv1::GetImageBit())) {
         nImages++; 
-       max = TMath::Max(max, hist->GetMaximum());
-       min = TMath::Min(min, hist->GetMinimum());
+       TString name(hist->GetName());
+       if (name.Contains("readouterrors", TString::kIgnoreCase)) continue;
+       Double_t hMax = hist->GetMaximum(); 
+       // hist->GetBinContent(hist->GetMaximumBin());
+       Double_t hMin = hist->GetMinimum();
+       // hist->GetBinContent(hist->GetMinimumBin());
+       max = TMath::Max(max, hMax);
+       min = TMath::Min(min, hMin);
+       AliInfo(Form("Histogram %30s min=%f, max=%f (min=%f,max=%f)",
+                    name.Data(), hMin, hMax, min, max));
       }
     }
     break ; 
   }
   min = TMath::Max(0.1, min);
-  max = TMath::Min(1.0, max);
+  max = TMath::Max(1.0, max);
 
   if (nImages == 0) {
     AliDebug(AliQAv1::GetQADebugLevel(), 
@@ -144,16 +207,18 @@ AliFMDQAChecker::MakeImage(TObjArray** list,
   AliDebug(AliQAv1::GetQADebugLevel(), 
           Form("%d histograms will be plotted for %s %s\n", 
                nImages, GetName(), AliQAv1::GetTaskName(task).Data()));  
-
+  gStyle->SetOptStat(0);
   for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
-    if (!AliQAv1::Instance()->IsEventSpecieSet(specie)) continue ;
+    if (!AliQAv1::Instance(AliQAv1::GetDetIndex(GetName()))
+       ->IsEventSpecieSet(specie)) continue;
     
-    if(!list[specie]) continue;
+    if(!list[specie] || list[specie]->GetEntries() <= 0) continue;
 
     const Char_t * title = Form("QA_%s_%s_%s", GetName(), 
                                AliQAv1::GetTaskName(task).Data(), 
                                AliRecoParam::GetEventSpecieName(specie)); 
-    if (!fImage[specie]) fImage[specie] = new TCanvas(title, title) ;
+    if (!fImage[specie]) 
+      fImage[specie] = new TCanvas(title, title) ;
     fImage[specie]->Clear() ; 
     fImage[specie]->SetTitle(title) ; 
     fImage[specie]->cd() ; 
@@ -161,10 +226,11 @@ AliFMDQAChecker::MakeImage(TObjArray** list,
     TPaveText someText(0.015, 0.015, 0.98, 0.98) ;
     someText.AddText(title) ;
     someText.Draw() ; 
-    fImage[specie]->Print(Form("%s%s%d.%s", AliQAv1::GetImageFileName(), 
-                               AliQAv1::GetModeName(mode), 
-                               AliQAChecker::Instance()->GetRunNumber(), 
-                               AliQAv1::GetImageFileFormat()), "ps") ; 
+    TString outName(Form("%s%s%d.%s", AliQAv1::GetImageFileName(), 
+                        AliQAv1::GetModeName(mode), 
+                        AliQAChecker::Instance()->GetRunNumber(), 
+                        AliQAv1::GetImageFileFormat()));
+    fImage[specie]->Print(outName, "ps") ; 
     fImage[specie]->Clear(); 
 
     Int_t nx = int(nImages + .5) / 2;
@@ -172,24 +238,61 @@ AliFMDQAChecker::MakeImage(TObjArray** list,
     fImage[specie]->Divide(nx, ny, 0, 0);
     
     
-    TH1F* hist  = 0;
+    TH1*  hist  = 0;
     Int_t nHist = list[specie]->GetEntriesFast();
+    // list[specie]->Print();
     Int_t j     = 0;
     for (Int_t i = 0; i < nHist; i++) { 
-      hist = static_cast<TH1F*>(list[specie]->At(i));
-      if (!(hist && hist->TestBit(AliQAv1::GetImageBit()))) continue;
+      hist = static_cast<TH1*>(list[specie]->At(i));
+      if (!hist || !hist->TestBit(AliQAv1::GetImageBit())) continue;
       
       TVirtualPad* pad = fImage[specie]->cd(++j);
-      pad->SetLogy();
-      hist->SetMinimum(min);
-      hist->SetMaximum(max);
-      hist->DrawCopy();
+      pad->SetRightMargin(0.001);
+      Int_t logOpts = 0;
+      logOpts |= CheckForLog(hist->GetXaxis(), pad, 1);
+      logOpts |= CheckForLog(hist->GetYaxis(), pad, 2);
+      logOpts |= CheckForLog(hist->GetZaxis(), pad, 3);
+      if (hist->GetEntries() > 0)
+       AliInfo(Form("Drawing %30s (min=%f, max=%f) in pad %d",
+                    hist->GetName(), min, max, j));
+
+      TString opt;
+      TString name(hist->GetName());
+      if (name.Contains("readouterrors", TString::kIgnoreCase)) {
+       pad->SetRightMargin(0.15);
+       pad->SetBottomMargin(0.10);
+       pad->SetTopMargin(0.02);
+       opt="COLZ";
+      }
+      else {
+       hist->SetMinimum(min);
+       hist->SetMaximum(max);
+
+      }
+      hist->DrawCopy(opt);
+
+      if (name.Contains("readouterrors", TString::kIgnoreCase)) 
+       continue;
+
+      TPad* insert = new TPad("insert", "Zoom", 
+                             .4,.4, .99, .95, 0, 0, 0);
+      insert->SetTopMargin(0);
+      insert->SetRightMargin(0.001);
+      insert->SetFillColor(0);
+      insert->SetBorderSize(1);
+      insert->SetBorderMode(0);
+      insert->Draw();
+      insert->cd();
+      if (logOpts & 0x1) insert->SetLogx();
+      if (logOpts & 0x2) insert->SetLogy();
+      if (logOpts & 0x4) insert->SetLogz();
+      hist->GetXaxis()->SetRangeUser(hist->GetXaxis()->GetXmin(), 
+                                    hist->GetXaxis()->GetXmax()/8);
+      hist->DrawCopy(opt);
+      pad->cd();
     }
 
-    fImage[specie]->Print(Form("%s%s%d.%s", AliQAv1::GetImageFileName(), 
-                              AliQAv1::GetModeName(mode), 
-                              AliQAChecker::Instance()->GetRunNumber(), 
-                              AliQAv1::GetImageFileFormat()), "ps"); 
+    fImage[specie]->Print(outName, "ps");
   }
 }
 
index 21f2bcb..a5c32bd 100644 (file)
@@ -45,7 +45,27 @@ public:
                  AliQAv1::TASKINDEX_t task, 
                  AliQAv1::MODE_t mode);
 
-private:
+protected:
+  /** 
+   * Check one histogram 
+   * 
+   * @param specie 
+   * @param hist 
+   * 
+   * @return 
+   */
+  Double_t CheckOne(AliQAv1::ALITASK_t          what,
+                   AliRecoParam::EventSpecie_t specie, 
+                   TH1*                        hist) const;
+  Double_t CheckRaw(AliRecoParam::EventSpecie_t specie, 
+                   TH1*                        hist) const;
+  Double_t CheckSim(AliRecoParam::EventSpecie_t specie, 
+                   TH1*                        hist) const;
+  Double_t CheckESD(AliRecoParam::EventSpecie_t specie, 
+                   TH1*                        hist) const;
+  Double_t CheckRec(AliRecoParam::EventSpecie_t specie, 
+                   TH1*                        hist) const;
+
   ClassDef(AliFMDQAChecker,0)  // Yves? what to do? 
 };
 
index c528da6..3f63e8e 100644 (file)
@@ -18,6 +18,7 @@
 #include <TFile.h> 
 #include <TH1F.h> 
 #include <TH1I.h> 
+#include <TH2I.h> 
 
 // --- AliRoot header files ---
 #include "AliESDEvent.h"
@@ -33,6 +34,9 @@
 #include "AliFMDAltroMapping.h"
 #include "AliFMDDebug.h"
 
+namespace {
+  Int_t colors[3] = {kRed,kGreen,kBlue};
+}
 //_____________________________________________________________________
 // This is the class that collects the QA data for the FMD during
 // reconstruction.  
@@ -112,16 +116,69 @@ AliFMDQADataMakerRec::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task,
 }
 
 //_____________________________________________________________________ 
+TH1* AliFMDQADataMakerRec::MakeADCHist(UShort_t d, Char_t r, Short_t b)
+{
+  TString name("adc"); 
+  TString title("ADC counts");
+  Int_t   color = kRed+1; 
+  if (d > 0) { 
+    name.Append(Form("FMD%d%c", d, r));
+    title.Append(Form(" in FMD%d%c", d, r));
+    color = colors[d-1]+3+(r == 'I' || r == 'i' ? 0 : 1);
+    if (b >= 0) {
+      name.Append(Form("_0x%02x", b));
+      title.Append(Form("[0x%02x]", b));
+    }
+  }
+  TH1* hist = new TH1F(name, title,1024,0,1024);
+  hist->SetXTitle("Amplitude [ADC counts]");
+  hist->SetYTitle("Events [log]");
+  hist->SetFillStyle(3001);
+  hist->SetFillColor(color);
+  hist->SetLineColor(color);
+  hist->SetMarkerColor(color);
+  hist->GetXaxis()->SetNdivisions(408,false);
+  // hist->SetStats(0);
+
+  return hist;
+}
+//_____________________________________________________________________ 
+TH1* AliFMDQADataMakerRec::MakeELossHist(UShort_t d, Char_t r, Short_t b)
+{
+  TString name("eloss"); 
+  TString title("Energy loss");
+  Int_t   color = kBlue+1;
+  if (d > 0) { 
+    name.Append(Form("FMD%d%c", d, r));
+    title.Append(Form(" in FMD%d%c", d, r));
+    color = colors[d-1]+3+(r == 'I' || r == 'i' ? 0 : 1);
+    if (b >= 0) {
+      name.Append(Form("_0x%02x", b));
+      title.Append(Form("[0x%02x]", b));
+    }
+  }
+  TH1* hist = new TH1F(name, title,300,0, 15);
+  hist->SetXTitle("#Delta E/#Delta_{mip}");
+  hist->SetYTitle("Events [log]");
+  hist->SetFillStyle(3001);
+  hist->SetFillColor(color);
+  hist->SetLineColor(color);
+  hist->SetMarkerColor(color);
+  // hist->SetStats(0);
+
+  return hist;
+}
+
+
+//_____________________________________________________________________ 
 void AliFMDQADataMakerRec::InitESDs()
 {
   // create Digits histograms in Digits subdir
   const Bool_t expert   = kTRUE ; 
   const Bool_t image    = kTRUE ; 
   
-  TH1F* hEnergyOfESD = new TH1F("hEnergyOfESD","Energy distribution",100,0,3);
-  hEnergyOfESD->SetXTitle("Edep/Emip");
-  hEnergyOfESD->SetYTitle("Counts");
-  Add2ESDsList(hEnergyOfESD, 0, !expert, image);
+  TH1* hist = MakeELossHist();
+  Add2ESDsList(hist, 0, !expert, image);
     
 }
 
@@ -132,10 +189,8 @@ void AliFMDQADataMakerRec::InitDigits()
   const Bool_t expert   = kTRUE ; 
   const Bool_t image    = kTRUE ; 
   
-  TH1I* hADCCounts = new TH1I("hADCCounts",
-                             "Dist of ADC counts;ADC counts;Counts",
-                             1024,0,1024);
-  Add2DigitsList(hADCCounts, 0, !expert, image);
+  TH1* hist = MakeADCHist();
+  Add2DigitsList(hist, 0, !expert, image);
 }
 
 
@@ -146,11 +201,8 @@ void AliFMDQADataMakerRec::InitRecPoints()
   const Bool_t expert   = kTRUE ; 
   const Bool_t image    = kTRUE ; 
 
-  TH1F* hEnergyOfRecpoints = new TH1F("hEnergyOfRecpoints",
-                                     "Energy Distribution",100,0,3);
-  hEnergyOfRecpoints->SetXTitle("Edep/Emip");
-  hEnergyOfRecpoints->SetYTitle("Counts");
-  Add2RecPointsList(hEnergyOfRecpoints,0, !expert, image);
+  TH1* hist = MakeELossHist();
+  Add2RecPointsList(hist,0, !expert, image);
 }
 
 //_____________________________________________________________________ 
@@ -160,42 +212,38 @@ void AliFMDQADataMakerRec::InitRaws()
   const Bool_t expert   = kTRUE ; 
   const Bool_t saveCorr = kTRUE ; 
   const Bool_t image    = kTRUE ; 
-  Int_t colors[3] = {kRed,kGreen,kBlue};
-  TH1I* hADCCounts;
-  for(Int_t det = 1; det<=3; det++) {
-    Int_t firstring = (det==1 ? 1 : 0);
-    for(Int_t iring = firstring;iring<=1;iring++) {
-      Char_t ring = (iring == 1 ? 'I' : 'O');
-      hADCCounts = 
-       new TH1I(Form("hADCCounts_FMD%d%c", det, ring), 
-                Form("FMD%d%c ADC counts;Amplitude [ADC counts];Counts",
-                     det,ring),1024,0,1023);
-      hADCCounts->SetFillStyle(3001);
-      hADCCounts->SetFillColor(colors[det-1]+2+iring);
-      hADCCounts->SetDrawOption("LOGY");
-
-      Int_t index1 = GetHalfringIndex(det, ring, 0,1);
-      Add2RawsList(hADCCounts, index1, !expert, image, !saveCorr);
+  TH2I* hErrors = new TH2I("readoutErrors", "Read out errors", 3, .5, 3.5,
+                          160, 0, 160); 
+  hErrors->GetXaxis()->SetBinLabel(1, "FMD1");
+  hErrors->GetXaxis()->SetBinLabel(2, "FMD2");
+  hErrors->GetXaxis()->SetBinLabel(3, "FMD3");
+  hErrors->SetYTitle("# errors");
+  hErrors->SetZTitle("Events [log]");
+  Add2RawsList(hErrors, 1, !expert, image, !saveCorr);
+  //AliInfo(Form("Adding %30s to raw list @ %2d", hErrors->GetName(), 1));
+
+  TH1* hist;
+  Int_t idx = 0;
+  for(UShort_t d = 1; d<=3; d++) {
+    UShort_t nR = (d == 1 ? 1 : 2); 
+    for(UShort_t q = 0; q < nR; q++) {
+      Char_t r = (q == 1 ? 'O' : 'I');
+      hist = MakeADCHist(d, r, -1);
+
+      Int_t index1 = GetHalfringIndex(d, r, 0, 1);
+      idx          = TMath::Max(index1, idx);
+      Add2RawsList(hist, index1, !expert, image, !saveCorr);
+      //AliInfo(Form("Adding %30s to raw list @ %2d", hist->GetName(), index1));
       
-      for(Int_t b = 0; b<=1;b++) {
-       
+      for(UShort_t b = 0; b <= 1; b++) {
        //Hexadecimal board numbers 0x0, 0x1, 0x10, 0x11;
-       UInt_t board = (iring == 1 ? 0 : 1);
-       board = board + b*16;
-       
-       
-       hADCCounts      = 
-         new TH1I(Form("hADCCounts_FMD%d%c_0x%02x", det, ring, board), 
-                  Form("FMD%d%c_0x%x ADC counts;Amplitude [ADC counts];Counts",
-                       det,ring, board),1024,0,1023);
-       hADCCounts->SetXTitle("ADC counts");
-       hADCCounts->SetYTitle("");
-       hADCCounts->SetFillStyle(3001);
-       hADCCounts->SetFillColor(colors[det-1]+2+iring);
-       hADCCounts->SetDrawOption("LOGY");
-       Int_t index2 = GetHalfringIndex(det, ring, board/16,0);
-       Add2RawsList(hADCCounts, index2, expert, !image, !saveCorr);
+       UShort_t board = (q == 1 ? 0 : 1) + b*16;
 
+       hist = MakeADCHist(d, r, board);
+       Int_t index2 = GetHalfringIndex(d, r, board/16,0);
+       idx          = TMath::Max(index2, idx);
+       Add2RawsList(hist, index2, expert, !image, !saveCorr);
+       //AliInfo(Form("Adding %30s to raw list @ %2d",hist->GetName(),index2));
       }
     }
   }
@@ -340,11 +388,15 @@ void AliFMDQADataMakerRec::MakeRaws(AliRawReader* rawReader)
     Short_t           board = pars->GetAltroMap()->Sector2Board(ring, sec);
     
     Int_t index1 = GetHalfringIndex(det, ring, 0, 1);
-    GetRawsData(index1)->Fill(digit->Counts());
     Int_t index2 = GetHalfringIndex(det, ring, board/16,0);
-    GetRawsData(index2)->Fill(digit->Counts());
-    
+    // AliInfo(Form("FMD%d%c[0x%02x]->%2d,%2d", 
+    //              det, ring, board, index1, index2));
+    GetRawsData(index1)->Fill(digit->Counts());
+    GetRawsData(index2)->Fill(digit->Counts()); 
   }
+  GetRawsData(1)->Fill(1, fmdReader.GetNErrors(0));
+  GetRawsData(1)->Fill(2, fmdReader.GetNErrors(1));
+  GetRawsData(1)->Fill(3, fmdReader.GetNErrors(2));
 }
 
 //_____________________________________________________________________
@@ -399,10 +451,16 @@ Int_t AliFMDQADataMakerRec::GetHalfringIndex(UShort_t det,
   // Return:
   //    Half ring index
   //  
-  UShort_t iring  =  (ring == 'I' ? 1 : 0);
-  Int_t index = ( ((det-1) << 3) | (iring << 2) | (board << 1) | 
-                 (monitor << 0));
-  
+  UShort_t iring = (ring == 'I' || ring == 'i' ? 1 : 0);
+  Int_t    index = ((((det-1) & 0x3) << 3) | 
+                   ((iring   & 0x1) << 2) |
+                   ((board   & 0x1) << 1) | 
+                   ((monitor & 0x1) << 0));
+#if 0
+  AliInfo(Form("d=%d, r=%c, b=%d, m=%d -> (%d<<3)|(%d<<2)|(%d<<1)|(%d<<0)=%2d",
+              det, ring, board, monitor, (det-1) & 0x3, iring & 0x1, 
+              board & 0x1, monitor & 0x1, index));
+#endif
   return index-2;
 }
 
index 8ab088f..d775b67 100644 (file)
@@ -44,6 +44,9 @@ public:
    */
   virtual ~AliFMDQADataMakerRec();
 private:
+  static TH1* MakeADCHist(UShort_t d=0, Char_t r='\0', Short_t b=-1);
+  static TH1* MakeELossHist(UShort_t d=0, Char_t r='\0', Short_t b=-1);
+  
   /** 
    * Called at end of monitor cycle 
    * 
@@ -114,6 +117,7 @@ private:
                         UShort_t board, UShort_t monitor = 0) const;
   ClassDef(AliFMDQADataMakerRec,0)  // description 
   TClonesArray fRecPointsArray; // Rec points
+
 };
 
 #endif // AliFMDQADataMakerRec_H