MCoincidence.hh

Go to the documentation of this file.

#ifndef _M_COINCIDENCE_HH_
#define _M_COINCIDENCE_HH_
 
#include "QModule.hh"
#include <map>
#include "QRunData.hh"
 
class MCoincidence : public Diana::QModule {
  
public:
  
  void Init(Diana::QEvent& ev);
  
  void Do(Diana::QEvent& ev, const Diana::QEventList& neighbours);
  
  void Done();
  
  double GetTime(const Diana::QEvent& ev);
private:
  enum TimeVariable {
    TV_RAW,
    TV_OF,
  };
  
  void GetEventEnergy(const Diana::QEvent &ev,double &energy,bool &valid) const;
  double GetDistanceMm(const Diana::QEvent &ev, const Diana::QEvent* coincident);
  
  int fDataset;
  bool fUseRunningWindow;
  bool fSyncTime;
  std::string fSyncType;
  bool fOkInvalidTotalEnergy;
  std::string fSyncSource;
  std::string fSyncTowerSource;
  double fCoincidenceWindow;
  std::string fEnergyLabel;
  std::vector<std::string> fPSALabels;
  TimeVariable fTimeVariable;
  QRunData fRunData;
  int fThisRun;
  struct ChannelInfo {
    bool BlackSheep;
    double DeltaT;
  };
  std::map<int,ChannelInfo> fSyncTimeDifference;
  bool fStorePosition;
  std::set<int> fBlackList;
  std::string fOFLabel; //label for the optimal filter (was COF in DIANA-0, different for DIANA)
  double fCoincidenceRadius;
  std::map<int,double> fChannelPositionX;
  std::map<int,double> fChannelPositionY;
  std::map<int,double> fChannelPositionZ;
};
 
#endif