DetectorDataStructures.h

Go to the documentation of this file.

#ifndef _DETECTOR_DATA_STRUCTURES_H_
#define _DETECTOR_DATA_STRUCTURES_H_
#include <vector>
#include <cstdint>
using std::vector;
 
#include "MDPP16DataStructures.h"
#include "ConfigClasses.h"
 
 
class MDPPDetBase
{
public:
  MDPPDetBase();
  virtual void defineClass(MDPP16Event& a_ev,
                           int          a_detID,
                           unsigned int a_ADC,
                           unsigned int a_TDC,
                           int          a_trigInd = 0,
                           int          a_rfInd = 1
                         );
  //this constructor works similarly to the 
  //one above but uses a specified rf array potentiall from another board 
  virtual void defineClass(MDPP16Event&      a_ev,
                           int               a_detID,
                           unsigned int      a_ADC,
                           unsigned int      a_TDC,
                           vector<uint16_t>& a_trigArray,
                           int               a_trigInd = 0
                         );
 
  //getters for internal class member data
  uint16_t getADC() const;
  uint16_t getTDC() const;
  uint32_t getendEventWord() const;
  uint16_t getextTime() const;
  uint8_t getdetNumber() const;
  uint8_t getmodID() const;
  uint16_t getRF() const;
  vector<uint16_t> gettdcTriggers() const;
  //setters for interanl class member data
  void setADC(uint16_t a_ADC);
  void setTDC(uint16_t a_TDC);
  void setendEventWord(uint32_t a_endEventWord);
  void setextTime(uint16_t a_extTime);
  void setdetNumber(uint8_t a_detNumber);
  void setmodID(uint8_t a_modID);
  void settdcTrigger(vector<uint16_t> a_tdcTrigger);
 
  double getGlobalTriggerTime();
  double getGlobalHitTime();
 
  uint16_t m_ADC;
  uint16_t m_TDC;
  uint32_t m_endEventWord;
  uint16_t m_extTime;
  uint16_t m_detNumber;
  uint8_t m_modID;
  uint16_t m_nearestRF;
  vector<uint16_t> m_tdcTrigger;
 
  std::ostream& print( std::ostream& a_stream );
};
 
 
class MDPPCloverEvent 
{
public:
  MDPPCloverEvent();
 
  MDPPCloverEvent(const std::vector<MDPP16Event>& a_daqEvs, 
                  const CloverInfo& a_config,
                  int   a_trigInd = 0,
                  int   a_rfInd   = 1
                 );
  MDPPCloverEvent( MDPP16Event* a_daqEvs, 
                  int                a_num,
                  const CloverInfo& a_config,
                  vector<uint16_t>& a_trigArray,
                  int   a_trigInd = 0
                 );
  void defineClass( MDPP16Event* a_daqEvs, 
                  int                a_num,
                  const CloverInfo& a_config,
                  vector<uint16_t>& a_trigArray,
                  int   a_trigInd = 0
                 );
  //basic getters for internal data members
  unsigned int getcloverID() const;
  vector<uint16_t> getADC() const;
  vector<uint16_t> getTDC() const;
  vector<uint16_t> gettdcBGO() const;
  vector<uint16_t> getadcBGO() const;
  vector<uint8_t> getleafNumber() const;
  vector<uint8_t> getbgoNumber() const;
  vector<uint16_t> gettdcTriggers() const;
  vector<uint16_t> getRFs() const;
  uint32_t getendEventWord() const;
  uint16_t getextTime() const;
  //basic setters for internal class member data, constructors should be 
  //preferred
  //basic getters for internal data members
  void setcloverID(uint8_t a_cloverID);
  void setADC(vector<uint16_t> a_ADC);
  void setTDC(vector<uint16_t> a_TDC);
  void settdcBGO(vector<uint16_t> a_tdcBGO);
  void setadcBGO(vector<uint16_t> a_adcBGO);
  void setleafNumber(vector<uint8_t> a_leafNumber);
  void setbgoNumber(vector<uint8_t> a_bgoNumber);
  void settdcTrigger(vector<uint16_t> a_tdcTrigger);
  void setendEventWord(uint32_t a_endEventWord);
  void setextTime(uint16_t a_extTime);
 
 
  //not sure why this class is entirely transparent to access at this point
  uint8_t m_cloverID;
  vector<uint16_t> m_ADC;
  vector<uint16_t> m_TDC;
  vector<uint16_t> m_tdcBGO;
  vector<uint16_t> m_adcBGO;
  vector<uint8_t> m_leafNumber;
  vector<uint8_t> m_bgoNumber;
  vector<uint16_t> m_nearestRF;
  vector<uint16_t> m_tdcTrigger;
  uint32_t m_endEventWord;
  uint16_t m_extTime;
 
 
  //more complex getters
  double getGlobalTriggerTime();
  double getGlobalHitTime(int a_hitNum);
  int getFirstHitLeaf(const CloverInfo& a_config) const;
  double getEnergy(int a_hitNum,const CloverInfo& a_config) const;
  double getTotalED(const CloverInfo& a_config,int& a_multiplicity) const;
  double getTotalED(const CloverInfo& a_config,
                                 int& a_multiplicity,
                                 int& a_minTimeHit,
                                 std::vector<int>& a_rejHits
                                 ) const;
 
  //tests if a bgo is in coincidence, returns true if it is
  bool   rejectBGO(const CloverInfo& a_config) const;
  vector<uint16_t> getTDC(int a_leafNumber) const;
  int getMinTDCValue(int& a_leafID);
 
  //this function returns the average time for the leafs minus the average time
  //for the bgo's and is generically usefull for generating a calibration
  //routine for them. It returns an abnormally large negative value if there is
  //no bgo present
  double getBGODtTDC();
 
  bool constructed();
  bool hasGamma();
  //this function returns true if tdc information for the bgo is present 
  //in the event
  bool hasBGO();
  std::ostream& print( std::ostream& a_stream = std::cout);
};
 
class MDPPScintEvent : public MDPPDetBase
{
public:
  MDPPScintEvent();
  void defineClass(MDPP16Event& a_ev,
                   int a_detID,
                   unsigned int a_ADC,
                   unsigned int a_pkADC,
                   unsigned int a_TDC
                 );
  void defineClass(MDPP16Event& a_ev,
                   int a_detID,
                   unsigned int a_ADC,
                   unsigned int a_pkADC,
                   unsigned int a_TDC,
                   std::vector<uint16_t> a_rfVector
                 );
 
  // MDPPScintEvent(const MDPP16Event& a_daqEv,
  //                 const ScintConfig& a_config,
  //                 int a_eventNumber);
  double getTTPSD() const;
  void setTTPSD(unsigned int a_ADC, unsigned int a_peakADC);
 
  float m_TTPSD;
};
 
 
class MDPPHPGeEvent : public MDPPDetBase
{
public:
  MDPPHPGeEvent();
 
};
 
class MDPPInorgScintEvent : public MDPPDetBase
{
public:
  MDPPInorgScintEvent();
};
 
class DETHelperFunctions
{
public:
  static std::vector<double> getDeltaT(const MDPPCloverEvent& a_clovEvent,
                                       const MDPPScintEvent&  a_scintEvent,
                                       int                    a_leafNumber
                                      );
  static double getDeltaT(const MDPPHPGeEvent& a_hpgeEvent,
                                       const MDPPScintEvent&  a_scintEvent
                                      );
  static double getDeltaT(const MDPPInorgScintEvent& a_inorgEvent,
                                       const MDPPScintEvent&  a_scintEvent
                                      );
  static double getDeltaT(const MDPPScintEvent& a_scintEvent1,
                                       const MDPPScintEvent&  a_scintEvent2
                                      );
  static double getNTOF(const MDPPCloverEvent& a_clovEvent,
                        const MDPPScintEvent&  a_scintEvent,
                        const CloverInfo&      a_clovInfo,
                        const OrgScintInfo&    a_scintInfo
                        );
 
  static std::vector<double> getRFDeltaT(const MDPPCloverEvent& a_clovEvent,
                                         const std::vector<uint16_t>&,
                                         const CloverInfo&      a_clovInfo
                                         );
};
#endif