MuonSimpleAlg.h

Go to the documentation of this file.

#ifndef MUONSIMPLEALG_H
#define MUONSIMPLEALG_H 1

#include "DybAlg/DybAlgorithm.h"
#include "Event/UserDataHeader.h"
#include "Event/ReadoutHeader.h"
#include "Event/GenHeader.h"
#include "Event/CalibReadoutHeader.h"
#include "DataSvc/ICableSvc.h"
#include "DataSvc/ICalibDataSvc.h"

#include "DetHelpers/ICoordSysSvc.h"
#include "DetDesc/DetectorElement.h"
#include <vector>
#include <map>
using namespace std;
using namespace DayaBay;

namespace DayaBay{
  class HeaderObject;
  class RecMuon;
}

class MuonSimpleAlg : public DybAlgorithm<DayaBay::UserDataHeader>
{
 public:  
  MuonSimpleAlg(const std::string& name, ISvcLocator* pSvcLocator);
  
  virtual ~MuonSimpleAlg();
  
  virtual StatusCode initialize();
  virtual StatusCode execute();
  virtual StatusCode finalize();

 private:
  //define PMT info
  struct PmtProp{
    Gaudi::XYZPoint Position;
    double Charge;
    double FirstHit;
    //double ExpectTime;
    //int Alive;
    //int Used;
    
  };

  map<unsigned int/*PmtId*/,PmtProp> m_PmtPropMap;
  int           m_PmtNumber;
  map<unsigned int,int> m_PmtId;//transform pmtId from unsigned int to int (0,1,....)  
  // CableSvc
  ICableSvc* m_cableSvc;
  
  // Pmt calibration data Service
  ICalibDataSvc* m_calibDataSvc;

  // find the detector which contain a global point
  ICoordSysSvc* mCss;
  
  string m_site;
  IDetectorElement *m_de;
  float m_timeCutLow;
  float m_timeCutHigh;
};

#endif // MUONSIMPLEALG_H