DQA_MUC.h

Go to the documentation of this file.

#ifndef Physics_Analysis_DQA_MUC_H
#define Physics_Analysis_DQA_MUC_H
 
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/NTuple.h"
#include "TF1.h"
#include "TH1F.h"
 
#include "MucCalib/MucIdTransform.h"
#include "MucCalib/MucMark.h"
#include "MucCalib/MucStructConst.h"
 
class ITHistSvc;
const int  TAGN            = 2; // 0-all; 1-dimu
const int  LVLN            = 3; // 0-layer; 1-box; 2-strip
const char ENAME[TAGN][10] = { "All", "Dimu" };
const char HTYPE[2][10]    = { "Num", "Den" };
const char LNAME[LVLN][10] = { "Layer", "Box", "Strip" };
 
class DQA_MUC : public Algorithm {
 
public:
  DQA_MUC( const std::string& name, ISvcLocator* pSvcLocator );
  StatusCode initialize();
  StatusCode execute();
  StatusCode finalize();
 
private:
  StatusCode FillHistograms( bool isDimu );
 
private:
  // Interface
  ITHistSvc* m_thsvc;
  long       m_run, m_event;
  bool       m_dstFileOnly;
  int        m_effWindow;
  // Muc variables
  // Digi and hit container for total events
  // 0-digi number, 1-rec track number, 2-efficiency hit number,
  // 3-incidental and noise hit number, 4-rec/fit hit number
  long m_recordAll[PART_MAX][SEGMENT_MAX][LAYER_MAX][STRIP_INBOX_MAX][5];
  long m_recordDimu[PART_MAX][SEGMENT_MAX][LAYER_MAX][STRIP_INBOX_MAX][5];
 
  // Digi or hit mark collection for one event
  MucMark* m_ptrMucMark;
  mark_col m_digiCol;   // container of digis in an event
  mark_col m_expHitCol; // container of fitting hits in all rec tracks in an  event
  mark_col m_effHitCol; // container of fired hits in all rec tracks in an event
  mark_col m_nosHitCol; // container of noise hits in all rec tracks in an event
 
  mark_col m_segDigiCol[PART_MAX][SEGMENT_MAX]; // container of digis in segment of part
  vector<mark_col> m_clusterCol; // container of clusters in all digits of an event
  vector<float>    m_lineResCol; // container of residual by line fit
 
  // Muc id transformer
  MucIdTransform* m_ptrIdTr;
 
  // Spacial resolution
  TH1F* m_hBrResDist[B_LAY_NUM][TAGN]; // barrel spacial resolution
  TH1F* m_hEcResDist[E_LAY_NUM][TAGN]; // endcap spacial resolution
  //  TH1F* m_hBrRes[TAGN];                 // barrel spacial resolution vs layer
  //  TH1F* m_hEcRes[TAGN];                 // endcap spacial resolution vs layer
 
  // Efficiency and noise ratio
  TH1F* m_hEff[LVLN][TAGN][2];
  TH1F* m_hNosRatio[TAGN][2];
 
  // Costheta and phi of MUC track position on 1st layer
  TH1F* m_hCostheta[TAGN];
  TH1F* m_hPhi[TAGN];
};
 
#endif