Initialisation.

[u/mrichter/AliRoot.git] / MUON / AliMUONTriggerEfficiencyCells.h

diff --git a/MUON/AliMUONTriggerEfficiencyCells.h b/MUON/AliMUONTriggerEfficiencyCells.h
index 72a9bcea9ee87e400e7e66cee7ce94e8f37a7983..03816e30b846d1f195e4c51009672a8902369833 100755 (executable)
--- a/MUON/AliMUONTriggerEfficiencyCells.h
+++ b/MUON/AliMUONTriggerEfficiencyCells.h
@@ -11,50 +11,84 @@
 #define ALIMUONTRIGGEREFFICIENCYCELLS_H
 
 #include "TObject.h"
+#include "TArrayF.h"
 #include "TArrayI.h"
 #include "TVector2.h"
-#include "TString.h"
+#include "TMatrix.h"
+#include "TH1F.h"
+#include "TList.h"
 
 class AliMUONTriggerEfficiencyCells : public TObject
 {
 public:
   AliMUONTriggerEfficiencyCells();
-  AliMUONTriggerEfficiencyCells(const char* filename);
+  AliMUONTriggerEfficiencyCells(const Char_t* filename);
+  AliMUONTriggerEfficiencyCells(TList *countHistoList, TList *noCountHistoList);
+
+  AliMUONTriggerEfficiencyCells(const AliMUONTriggerEfficiencyCells& other); // copy constructor
+  AliMUONTriggerEfficiencyCells& operator=(const AliMUONTriggerEfficiencyCells& other); // assignment operator
 
   virtual ~AliMUONTriggerEfficiencyCells();
 
-  void GetCellEfficiency(Int_t detElemId, Float_t x, Float_t y, Float_t &eff1, Float_t &eff2);
-  void GetCellEfficiency(Int_t detElemId, Int_t localBoard, Float_t &eff1, Float_t &eff2);
+  void GetCellEfficiency(Int_t detElemId, Float_t x, Float_t y, Float_t &eff1, Float_t &eff2) const;
+  void GetCellEfficiency(Int_t detElemId, Int_t localBoard, Float_t &eff1, Float_t &eff2) const;
     
-  void IsTriggered(Int_t detElemId, Float_t x, Float_t y, Bool_t &trig1, Bool_t &trig2);
-  void IsTriggered(Int_t detElemId, Int_t localBoard, Bool_t &trig1, Bool_t &trig2);
+  void IsTriggered(Int_t detElemId, Float_t x, Float_t y, Bool_t &trig1, Bool_t &trig2) const;
+  void IsTriggered(Int_t detElemId, Int_t localBoard, Bool_t &trig1, Bool_t &trig2) const;
 
-  TVector2 ChangeReferenceFrame(Float_t x, Float_t y, Float_t x0, Float_t y0);
+  void DisplayEfficiency(Bool_t perSlat=kFALSE, const Char_t* geoFilename="geometry.root");
+  Bool_t SumRunEfficiency(const AliMUONTriggerEfficiencyCells &other);
 
+  /// Set the list of fired strips
+  void SetFiredStrips(TList *firedStrips){fFiredStrips = firedStrips;}
+  void CheckFiredStrips(const Char_t *geoFilename="geometry.root");
+                                // Check for strips with lower counts than others:
+                                // syntomatic of possible read-out problems in boards
   void Reset();
     
 protected:
-    TArrayI CellByCoord(Int_t detElemId, Float_t x, Float_t y);
-    void ReadFile(const char* filename="$ALICE_ROOT/MUON/data/efficiencyCells.dat");
+    TArrayI CellByCoord(Int_t detElemId, Float_t x, Float_t y) const;
+    TVector2 ChangeReferenceFrame(Float_t x, Float_t y, Float_t x0, Float_t y0);
+    void ReadFile(const Char_t* filename="$ALICE_ROOT/MUON/data/efficiencyCells.dat");
+    void CalculateEfficiency(Int_t trigger44, Int_t trigger34,
+                            Float_t &efficiency, Float_t &error,
+                            Bool_t failuresAsInput);
+
 
 private:
-    Int_t FindChamberIndex(Int_t detElemId);
-    Int_t FindSlatIndex(Int_t detElemId);
+    void CheckConstants() const;
+    Int_t FindChamberIndex(Int_t detElemId) const;
+    Int_t FindSlatIndex(Int_t detElemId) const;
     void ReadFileXY(ifstream &file);
     void ReadFileBoards(ifstream &file);
+    void ReadHistoBoards(const Char_t* filename="MUON.TriggerEfficiencyMap.root");
+    void InitHistos();
+    void FillHistosFromList();
+    Bool_t GetListsForCheck(const Char_t* geoFilename="geometry.root");
+    
+    static const Int_t fgkNcells=80;   ///< Number of cells
+    static const Int_t fgkNcathodes=2; ///<Number of cathodes
+    static const Int_t fgkNchambers=4; ///<Number of chambers
+    static const Int_t fgkNplanes=8;   ///<Number of planes
+    static const Int_t fgkNslats=18;   ///<Number of slats
+
     
-    static const Int_t fgkNofCells=80; ///< Number of cells
+    TMatrixF fCellContent[fgkNplanes][fgkNslats]; ///< the cells content
+    TArrayF  fCellSize[fgkNplanes];    ///< the size of the cells
+    TArrayI  fCellNumber[fgkNplanes];  ///< id of the cells
     
-    /// the cells content
-    Float_t fCellContent[4][18][2][fgkNofCells][fgkNofCells]; //[trig. chambers][RPCs][cathode][cellsX][cellsY]
+    TH1F *fBoardEfficiency[fgkNplanes];///< the boards content
+    TH1F *fSlatEfficiency[fgkNplanes];///< the slats content
 
-    Float_t fCellSize[4][18][2]; ///< the size of the cells
-    Int_t fCellNumber[4][18][2]; ///< id of the cells
+    TList *fCountHistoList; ///<list of efficiency numerators
+    TList *fNoCountHistoList; ///<list of efficiency denominators
+    TList *fFiredStrips; ///<list of fired strips for efficiency check
 
-    static const Int_t fgkNofBoards=234; ///< Number of boards
-    /// the boards content
-    Float_t fBoardContent[4][2][fgkNofBoards]; //[trig. chambers][RPCs][cathode][board]
-    
-    ClassDef(AliMUONTriggerEfficiencyCells,2) // Trigger efficiency store
+    TList *fDisplayHistoList; //!< list of efficiency histograms for display
+    TList *fBoardLabelList; //!< list of board labels for display
+    TList *fFiredFitHistoList; //!< list of fired strips for checks
+    TList *fFiredDisplayHistoList; //!< list of fired strips for display
+
+    ClassDef(AliMUONTriggerEfficiencyCells,5) // Trigger efficiency store
 };
 #endif