1
/* This file is part of MAUS: http://micewww.pp.rl.ac.uk:8080/projects/maus
3
* MAUS is free software: you can redistribute it and/or modify
4
* it under the terms of the GNU General Public License as published by
5
* the Free Software Foundation, either version 3 of the License, or
6
* (at your option) any later version.
8
* MAUS is distributed in the hope that it will be useful,
9
* but WITHOUT ANY WARRANTY; without even the implied warranty of
10
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
11
* GNU General Public License for more details.
13
* You should have received a copy of the GNU General Public License
14
* along with MAUS. If not, see <http://www.gnu.org/licenses/>.
18
#ifndef _MAUS_INPUTCPPDATA_FADCDATAPROCESSOR_HH_
19
#define _MAUS_INPUTCPPDATA_FADCDATAPROCESSOR_HH_
25
#include "json/json.h"
27
#include "unpacking/MDprocessor.h"
29
#include "Utils/DAQChannelMap.hh"
30
#include "DataStructure/DAQData.hh"
32
////////////////////////////////////////////////////////////////////////////////
35
* Simple class to hold and arrange the Json document
36
* containing all the daq spill information.
41
class MDarranger : public MDprocessor {
44
virtual ~MDarranger() {}
46
/** Sets the json document to which the MDarranger will write
48
void set_JSON_doc(Json::Value *doc) {_docSpill = doc;}
50
/** Gets the json document where DAQ data is written
52
Json::Value* get_JSON_doc() {return _docSpill;}
54
/** This function sets the pointer to the DAQData object.
55
* \param[in] dd The pointer to the DAQData objevt.
57
static void set_daq_data(MAUS::DAQData *dd) {_daq_data = dd;}
59
/** Gets the json document where DAQ data is written
61
static MAUS::DAQData* get_daq_data() {return _daq_data;}
64
* This function sets the DAQ map.
65
* \param[in] map The DAQ channel map.
67
static void set_DAQ_map(DAQChannelMap* map) {_chMap = map;}
69
string get_equipment_name() {return _equipment;}
72
/** The JSON node to put the data under.
73
* It is created at the beginning of the spill.
75
Json::Value* _docSpill;
78
/**The DAQData object to put the data under.
79
* It is created at the beginning of the spill.
82
static MAUS::DAQData *_daq_data;
84
/** The DAQ channel map.
85
* It is used to group all channels belonging to a given detector.
87
static DAQChannelMap* _chMap;
89
/** The name of the DAQ equipment. **/
93
class ZeroSupressionFilter : public MDarranger {
97
ZeroSupressionFilter() :_zero_suppression(0), _zs_threshold(0) {}
98
virtual ~ZeroSupressionFilter() {}
100
void set_zero_supression(bool zs) { _zero_suppression = zs; }
101
void set_zs_threshold(int zst) { _zs_threshold = zst; }
105
/// If true - do zero suppression
106
bool _zero_suppression;
108
/// Value of the threshold used for zero suppression */
112
////////////////////////////////////////////////////////////////////////////////
115
// do not increment inside following macro! (e.g.: MAX( ++a, ++b );
116
#define MAX(a, b) (((a) > (b)) ? (a) : (b))
118
#define SIGNAL_INTEGRATION_WINDOW 30
120
class fADCDataProcessor : public ZeroSupressionFilter {
123
fADCDataProcessor() {}
124
~fADCDataProcessor() {}
126
/** Return the position of the maximum.
127
* This function returns the number of the sample having maximum amplitude.
129
int get_max_position();
131
/** Return the position of the minimum.
132
* This function returns the number of the sample having minimum amplitude.
134
int get_min_position();
136
int get_arrival_time();
140
/** Return the area of the signal.
141
* This function returns the area of the signal. The pedestal is subtracted.
142
* The integration is done in window starting 10 samples before the maximum
143
* and ending 20 samples after the maximum.
144
* \param[out] pos The number of the sample having maximum amplitude.
146
int get_pos_signal_area(int &pos);
148
int get_neg_signal_area(int &pos);
150
int get_pedestal_area(int &pos);
152
/** Return the measured value for the given sample
153
* \param[in] i The number of the sample.
155
int get_point(int i) const { return _data[i]; }
157
/// Return the vector of samples
158
vector<int> get_data() const { return _data; }
160
/** Return the area of the signal.
161
* This function returns the area of the signal. The pedestal is subtracted.
162
* The integration is done using the whole acquisition window.
167
/// Return the data member _pedestal
168
double get_pedestal() const { return _pedestal; }
170
/** Return the charge of the signal.
171
* This function returns the charge of the signal calculated by the specified
173
* \param[in] Algorithm identifier of the algorithm.
175
int get_charge(int Algorithm = ceaPedMax);
177
enum chargeEstimationAlgorithm {
178
ceaMinMax, /// Simplest algorithm
179
ceaFractionDescriminatorThreshold, /// not implemented
184
/**Return the charge of the Cherenkov Singal (negative pulse)
187
int get_charge_ckov(int Algorithm = ceaPedMin);
189
Json::Value get_samples();
197
/// vector of samples (measurements) */
200
/// Pedestal level of the signal */
added
removed
src/input/InputCppDAQData/fADCDataProcessor.hh
