ConvRawData.h

Go to the documentation of this file.

#ifndef CONVRAWDATA_H_
#define CONVRAWDATA_H_
 
#include <TClonesArray.h>
#include <TFile.h>
#include <TMap.h>
#include <TTree.h>
#include "FairTask.h"           // for FairTask, InitStatus
#include "FairEventHeader.h"    // for FairEventHeader
#include "SNDLHCEventHeader.h"  // for EventHeader
#include "Scifi.h"              // for Scifi detector
#include "sndScifiHit.h"    // for SciFi Hit
#include "MuFilterHit.h"    // for Muon Filter Hit
 
#include <iostream>
#include <tuple>
#include <map>
#include <vector>
 
class ConvRawData : public FairTask
{
  public:
    ConvRawData();
 
    ~ConvRawData();
 
    virtual InitStatus Init();
 
    virtual void Exec(Option_t* opt);
 
    void UpdateInput(int n);
 
    private:
      void StartTimeofRun(std::string path);
      void DetMapping(std::string path);
      void checkBoardMapping(std::string path);
      void debugMapping(std::string board, int tofpetID, int tofpetChannel);
      double qdc_calibration( int board_id, int tofpet_id, int channel, int tac, 
                              uint16_t v_coarse, uint16_t v_fine, uint16_t tf);
      double qdc_chi2( int board_id, int tofpet_id, int channel, int tac, int TDC);
      double qdc_sat( int board_id, int tofpet_id, int channel, int tac, uint16_t v_fine);
      double time_calibration( int board_id, int tofpet_id, int channel, int tac,
                               int64_t t_coarse, uint16_t t_fine, int TDC);
      std::tuple<double, double, double, double> comb_calibration( int board_id, int tofpet_id, 
                                                              int channel, int tac,
                                                              uint16_t v_coarse, uint16_t v_fine,
                                                              int64_t t_coarse, uint16_t t_fine,
                                                              double GQDC , int TDC);
                                                              // max gain QDC = 3.6
      std::map<double, std::pair<double, double> > calibrationReport();
 
      int channel_func( int tofpet_id, int tofpet_channel, int position);
      void read_csv(std::string path);
      void Process0();
      void Process1();
    
      std::map<int, MuFilterHit* > digiMuFilterStore{};
      std::map<int, sndScifiHit* > digiSciFiStore{}; 
      std::map<std::vector<int>, std::map<std::string, double> > X_qdc{};
      std::map<std::vector<int>, std::map<std::string, double> > X_tdc{};
      std::map<std::string, std::map<std::string, std::map<std::string, int>> > boardMaps{};
      std::map<int, std::map<int, int> > MufiSystem{}; // <board_id_mu, <slot(tofpetID), s>
      std::map<int, std::string > slots = { {0,"A"}, {1,"A"}, {2,"B"}, {3,"B"},
                                  {4,"C"}, {5,"C"}, {6,"D"}, {7,"D"} };
      std::map<int, std::map<int, int> > TofpetMap{};
      std::map<std::string, std::map<std::string, std::map<std::string, std::string>> > boardMapsMu{};
      std::map<std::string, std::vector<int> > offMap{};// name is key, vector is first bar, number of sipm channels / bar and direction
      // For time monitoring
      std::map<std::string, double> counters = { {"N",0}, {"event",0}, {"qdc",0}, {"tdc",0}, {"chi2",0},
                                       {"make",0}, {"storage",0}, {"createScifi",0}, {"createMufi",0} };
  
      Scifi* ScifiDet;
      TFile* fOut;
      TTree* fEventTree;
      // Board_x data
      std::map<std::string, TTree*> boards{};
      int frunNumber;
      int fnStart, fnEvents;
      int fheartBeat;
      int debug, stop, makeCalibration, local, newFormat;
      std::string fpathCalib, fpathJSON; 
      int eventNumber;
      double chi2Max, saturationLimit;
      double runStartUTC;
      TMap* FSmap;
    
      SNDLHCEventHeader* fSNDLHCEventHeader;
      FairEventHeader* fEventHeader;
      TClonesArray* fDigiSciFi;
      TClonesArray* fDigiMuFilter;
    
      ConvRawData(const ConvRawData&);
      ConvRawData& operator=(const ConvRawData&);
 
      ClassDef(ConvRawData, 3);
};
#endif /* CONVRAWDATA_H_ */