DQA_MUC.h

Go to the documentation of this file.

#ifndef Physics_Analysis_DQA_MUC_H
#define Physics_Analysis_DQA_MUC_H
 
#include "GaudiKernel/AlgFactory.h"
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/NTuple.h"
 
#include "TH1F.h"
#include "TF1.h"
 
#include "MucCalibAlg/MucStructConst.h"
#include "MucCalibAlg/MucMark.h"
#include "MucCalibAlg/MucIdTransform.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