GenesisEventStructures.h

Go to the documentation of this file.

 
#ifndef GENESIS_DATA_STRUCTURES
#define GENESIS_DATA_STRUCTURES
 
//project includes 
#include "DetectorDataStructures.h"
#include "MDPP16DataStructures.h"
#include "MesyDAQPostProcessing.h" //needed for static rf filter function
#include "NSDPhysicsEvents.h"
#include "GENESISConfig.h"
//ROOT includes
#include "TH1.h"
//c++includes
#include <cstdint>
#include <vector>
#include <utility>
#include <deque>
using std::vector;
#include <map>
using std::map;
using std::pair;
 
class GenesisEvent
{
  public:
    GenesisEvent();
 
    // /// add leaf/BGO hit to clover event
    // void addMDPPCloverEvent(unsigned int a_ADC, unsigned int a_TDC, unsigned int a_chNum, 
    //                         unsigned int a_endWord, unsigned int a_extTimeStamp, GenConfig* a_detConfig);
 
    // /// add scintillator hit m_scintEvents
    // void addMDPPScintEvent(unsigned int a_total, unsigned int a_TDC, unsigned int a_peak, unsigned int a_chNum, 
    //                         unsigned int a_endWord, unsigned int a_extTimeStamp, GenConfig* a_detConfig);
 
    void addMDPPCloverEvent(unsigned int a_ADC,
                            unsigned int a_TDC,
                            unsigned int a_chNum,
                            unsigned int a_endWord,
                            unsigned int a_extTimeStamp,
                            int          a_cloverID,
                            bool         a_isBGO
                            );
 
    void addMDPPScintEvent(unsigned int     a_total,
                           unsigned int     a_TDC,
                           unsigned int     a_peak,
                           unsigned int     a_chNum,
                           unsigned int     a_endWord,
                           unsigned int     a_extTimeStamp,
                           int              a_modID,
                           vector<uint16_t> a_triggers
                           );
void addMDPP16Evs(MDPP16Event* a_start,
                                int a_num,
                                ExperimentConfig& a_config,
                                time_t fileStartSinceEpoch,
                                bool a_saveBGOSingles=false,
                                int a_rfTrigIDs=1
                                );
 
    void addMDPPCloverEvent(MDPPCloverEvent* a_cloverEvent);
    void addMDPPScintEvent(MDPPScintEvent* a_scintEvent);
    void addMDPPCloverEvent(MDPPCloverEvent& a_cloverEvent);
    void addMDPPScintEvent(MDPPScintEvent& a_scintEvent);
    void addMDPPHPGeEvent(MDPPHPGeEvent& a_event);
    void addMDPPInorgScintEvent(MDPPInorgScintEvent& a_event);
    void addBGOData(MDPPCloverEvent* a_cloverEvent, int a_cloverID);
    void addBGOData(MDPPCloverEvent& a_cloverEvent, int a_cloverID);
    void addRFTimes(int a_moduleID, vector<uint16_t>& a_RFTimes);
    bool setCoarseTime(unsigned int a_endWord, unsigned int a_extTimeStamp);
 
    void setWallClockTime(uint32_t a_t);
 
 
    MDPPCloverEvent getCloverEvent(int a_cloverID) const;
    map<int , MDPPCloverEvent> getCloverEvents() const;
    vector<MDPPScintEvent> getScintEvent() const;
    MDPPHPGeEvent getHPGeEvent(int a_detID) const;
    MDPPInorgScintEvent getInorgScintEvent(int a_detID) const;
    unsigned int getCoarseTime() const;
    double getMinTime() const;
    vector<int> getMinTimeChannel() const;
    int getNHits() const;
    std::vector<DetType> getDetTypes() const;
    std::map<DetType,int> getDetsMultiplicity() const;
    vector<unsigned int> getLeafChannelList() const;
    vector<unsigned int> getScintChannelList() const;
    vector<int> getInorgScintIDList() const;
    bool hasBGO(int modID) const;
    vector<int> getCloverIDList() const;
    vector<int> getCloverIDListLeaf() const;
    vector<int> getHPGeIDList() const;
    uint32_t getWallClockTime() const;
 
    //functions below here modify the class in some way that is not just
    //adding data
 
    void reduceStoredRFs(int a_moduleID, bool a_isClover=0);
 
    void sortEvents();
 
 
    bool isCoin();
 
    void reset();
 
    bool isEmpty();
 
    std::ostream&  print(std::ostream& a_stream = std::cout);
  private:
    //private for making sure build is currently independent of old RF methods
    vector<uint16_t> getRFTimes(int a_moduleID) const;
    bool m_hasGoodRF;
    unsigned int m_coarseTime;
 
    uint32_t m_wallClockTime;
    map<int , MDPPCloverEvent> m_cloverEvents;
    vector<MDPPScintEvent> m_scintEvents;
    vector<MDPPHPGeEvent> m_HPGeEvents;
    vector<MDPPInorgScintEvent> m_inorgScintEvents;
 
    //stores the trigger array for rf times for a given module
    map<int, vector<uint16_t>> m_RFTimes;
 
 
};
class GenesisCalibratedEvent : public NSDPhysicsEvents::CalibratedCoinEvent
{
public:
  GenesisCalibratedEvent();
 
  void defineClass(const GenesisEvent* a_ev,
                   const GenesisConfig& a_conf,
                   int a_detOption=0
                   );
 
 
 
 
  //this function takes unpacked data from a GenesisEvent and 
  //builds out the relevant information for calibrated clover event
  //and adds it to the class
  void addCloverEvent(int                 a_clovID,
                      const GenesisEvent* a_ev,
                      const GenesisConfig& a_conf
                      );
  void addOrgScintEvent(const GenesisEvent* a_hit,
                        const GenesisConfig& a_conf
                        );
 
};
 
 
 
 
#endif