/* This file is part of MAUS: http://micewww.pp.rl.ac.uk:8080/projects/maus * * MAUS is free software: you can redistribute it and/or modify * it under the terms of the GNU General Public License as published by * the Free Software Foundation, either version 3 of the License, or * (at your option) any later version. * * MAUS is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * GNU General Public License for more details. * * You should have received a copy of the GNU General Public License * along with MAUS. If not, see . * */ /** @class RealDataDigitization * D... * */ #ifndef REALDATADIGITIZATION_HH #define REALDATADIGITIZATION_HH // C headers #include #include #include #include #include #include // C++ headers #include #include #include #include #include #include #include #include // G4MICE from commonCpp #include "Interface/dataCards.hh" #include "Interface/MICEEvent.hh" // from old file #include "Config/MiceModule.hh" #include "src/common_cpp/DataStructure/SciFiDigit.hh" #include "src/common_cpp/DataStructure/SciFiEvent.hh" #include "src/common_cpp/DataStructure/ReconEvent.hh" #include "src/common_cpp/DataStructure/Spill.hh" namespace MAUS { class RealDataDigitization { public: typedef std::vector IntChannelArray; RealDataDigitization(); ~RealDataDigitization(); void initialise(double npe_cut, const std::string& map_file, const std::string& calib_file); /** @brief Processes a spill from DAQ * @params spill A SciFiSpill to be filled */ void process(Spill *spill); /** @brief Reads the VLSB data * @params input_event The JSON input vlsb event to be processed. * @params event The SciFiEvent to be filled. * @params TrackerDaq The TrackerDaq object to be created for T0. * @params TrackerDaq The TrackerDaq object to be created for T1. */ std::vector process_VLSB(int SpillNum, TrackerDaq* td); /** @brief Reads the VLSB data in the cosmic run format * @params input_event The JSON input vlsb event to be processed. * @params event The SciFiEvent to be filled. * @params TrackerDaq The TrackerDaq object to be created for T0. * @params TrackerDaq The TrackerDaq object to be created for T1. */ void process_VLSB_c(Json::Value input_event, SciFiEvent* event, TrackerDaq *tracker0daq_event, TrackerDaq *tracker1daq_event); /** @brief Reads in the calibration. */ bool load_calibration(std::string filename); /** @brief Loads the mapping. */ bool load_mapping(std::string file); /** @brief Converts read-out map into SciFi Fibre map */ bool get_StatPlaneChannel(int &board, int &bank, int &chan_ro, int &tracker, int &station, int &plane, int &channel, int &extWG, int &inWG, int &WGfib) const; /** @brief Reads the bad channel list from file. */ bool load_bad_channels(); /** @brief Returns value depends on the goodness of the channel. */ bool is_good_channel(int bank, int chan_ro) const; private: static const int _number_channels = 128; static const int _number_banks = 64; static const int _number_boards = 16; static const int _total_number_channels = 6403; double _npe_cut; // static const double _min = 0.000000001; /// Arrays containing calibration values for every channel in the 4 banks of the 16 boards. Json::Value _calibration[_number_banks][_number_channels]; /// This is an array storing the goodness of each channel. bool _good_chan[_number_banks][_number_channels]; /// This is for the mapping storage. IntChannelArray _board; IntChannelArray _bank; IntChannelArray _chan_ro; IntChannelArray _tracker; IntChannelArray _station; IntChannelArray _view; IntChannelArray _fibre; IntChannelArray _extWG; IntChannelArray _inWG; IntChannelArray _WGfib; }; // Don't forget this trailing colon!!!! } // ~namespace MAUS #endif