Jpp
Functions
JRemovePMT.cc File Reference
#include <string>
#include <iostream>
#include <iomanip>
#include <set>
#include "JDAQ/JDAQ.hh"
#include "JDAQ/JDAQPMTIdentifier.hh"
#include "JDAQ/JDAQTimeslice.hh"
#include "JDAQ/JDAQSummaryslice.hh"
#include "JDAQ/JDAQEvent.hh"
#include "JTrigger/JTriggerParameters.hh"
#include "JSupport/JMultipleFileScanner.hh"
#include "JSupport/JFileRecorder.hh"
#include "JSupport/JSupport.hh"
#include "JLang/JTypeSelector.hh"
#include "JLang/JObjectMultiplexer.hh"
#include "JLang/JObjectDemultiplexer.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"

Go to the source code of this file.

Functions

int main (int argc, char **argv)
 

Detailed Description

Example program to remove PMT(s) from data (and set corresponding rate to 0).

Author
mdejong

Definition in file JRemovePMT.cc.

Function Documentation

◆ main()

int main ( int  argc,
char **  argv 
)

Definition at line 60 of file JRemovePMT.cc.

61 {
62  using namespace std;
63  using namespace JPP;
64  using namespace KM3NETDAQ;
65 
66  JMultipleFileScanner<JDAQTimesliceTypes_t> inputFile;
67  JLimit_t& numberOfEvents = inputFile.getLimit();
68  JFileRecorder<JTriggerTypes_t> outputFile;
70  int debug;
71 
72  try {
73 
74  JParser<> zap("Example program to remove PMT(s) from data (and set corresponding rate to 0).");
75 
76  zap['f'] = make_field(inputFile);
77  zap['o'] = make_field(outputFile) = "abc.root";
78  zap['n'] = make_field(numberOfEvents) = JLimit::max();
79  zap['P'] = make_field(PMT) = JPARSER::initialised();
80  zap['d'] = make_field(debug) = 1;
81 
82  zap(argc, argv);
83  }
84  catch(const exception& error) {
85  FATAL(error.what() << endl);
86  }
87 
88 
89  cout.tie(&cerr);
90 
91  outputFile.open();
92 
93  {
94  JObjectMultiplexer<JDAQTimesliceTypes_t, JDAQTimeslice> in(inputFile);
95 
96  // put base class JDAQTimeslice at end of type list
97 
99 
100  JObjectDemultiplexer<JDAQTimeslice, daq_timeslicetypes_t> out(outputFile);
101 
102  for (counter_type counter = 0; in.hasNext() && counter != inputFile.getLimit(); ++counter) {
103 
104  STATUS("event: " << setw(10) << counter << '\r'); DEBUG(endl);
105 
106  JDAQTimeslice* timeslice = in.next();
107 
108  for (JDAQTimeslice::iterator frame = timeslice->begin(); frame != timeslice->end(); ++frame) {
109 
110  bool rm = false;
111 
112  for (int pmt = 0; pmt != NUMBER_OF_PMTS; ++pmt) {
113 
114  if (PMT.find(JDAQPMTIdentifier(frame->getModuleID(), pmt)) != PMT.end() ||
115  PMT.find(JDAQPMTIdentifier(-1, pmt)) != PMT.end()) {
116  rm = true;
117  }
118  }
119 
120  if (rm) {
121 
123 
124  for (JDAQSuperFrame::const_iterator hit = frame->begin(); hit != frame->end(); ++hit) {
125  buffer[hit->getPMT()].push_back(*hit);
126  }
127 
128  vector<JDAQHit> data;
129 
130  for (int pmt = 0; pmt != NUMBER_OF_PMTS; ++pmt) {
131 
132  if (PMT.find(JDAQPMTIdentifier(frame->getModuleID(), pmt)) == PMT.end() &&
133  PMT.find(JDAQPMTIdentifier(-1, pmt)) == PMT.end()) {
134  copy(buffer[pmt].begin(), buffer[pmt].end(), back_inserter(data));
135  }
136  }
137 
138  sort(data.begin(), data.end());
139 
140  frame->clear();
141  frame->add(data.size(), data.data());
142  }
143  }
144 
145  out.put(*timeslice);
146  }
147  STATUS(endl);
148  }
149  {
150  JMultipleFileScanner<JDAQEvent> in(inputFile);
151 
152  while (in.hasNext()) {
153 
154  STATUS("event: " << setw(10) << in.getCounter() << '\r'); DEBUG(endl);
155 
156  JDAQEvent* event = in.next();
157 
158  remove(event->getHits<JDAQTriggeredHit>(), PMT);
159  remove(event->getHits<JDAQSnapshotHit> (), PMT);
160 
161  outputFile.put(*event);
162  }
163  STATUS(endl);
164  }
165  {
166  JMultipleFileScanner<JDAQSummaryslice> in(inputFile);
167 
168  while (in.hasNext()) {
169 
170  STATUS("event: " << setw(10) << in.getCounter() << '\r'); DEBUG(endl);
171 
172  JDAQSummaryslice* summaryslice = in.next();
173 
174  for (JDAQSummaryslice::iterator frame = summaryslice->begin(); frame != summaryslice->end(); ++frame) {
175 
176  for (int pmt = 0; pmt != NUMBER_OF_PMTS; ++pmt) {
177 
178  if (PMT.find(JDAQPMTIdentifier(frame->getModuleID(), pmt)) != PMT.end() ||
179  PMT.find(JDAQPMTIdentifier(-1, pmt)) != PMT.end()) {
180  (*frame)[pmt].setValue(0.0);
181  }
182  }
183  }
184 
185  outputFile.put(*summaryslice);
186  }
187  STATUS(endl);
188  }
189 
190  JMultipleFileScanner<JTriggerParameters> io(inputFile);
191 
192  io >> outputFile;
193 
194  outputFile.close();
195 }
KM3NETDAQ::JDAQEvent
DAQ Event.
Definition: JDAQEvent.hh:34
JPARSER::initialised
Empty structure for specification of parser element that is initialised (i.e.
Definition: JParser.hh:63
std::vector
Definition: JSTDTypes.hh:12
JLANG::JAppend
Append to type list.
Definition: JTypeList.hh:61
JSUPPORT::JLimit_t
JLimit JLimit_t
Type definition of limit.
Definition: JLimit.hh:215
KM3NETDAQ::NUMBER_OF_PMTS
static const int NUMBER_OF_PMTS
Total number of PMTs in module.
Definition: JDAQ.hh:26
KM3NETDAQ::JDAQTimeslice
Data time slice.
Definition: JDAQTimeslice.hh:36
JPARSER::JParser
Utility class to parse command line options.
Definition: JParser.hh:1493
PMT
JDAQPMTIdentifier PMT
Command line options.
Definition: JDAQTimesliceSelector.cc:35
JAANET::copy
void copy(const Head &from, JHead &to)
Copy header from from to to.
Definition: JHead.cc:152
std::set
Definition: JSTDTypes.hh:13
JPP
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Definition: JAAnetToolkit.hh:37
KM3NETDAQ::JDAQTriggeredHit
DAQ triggered hit.
Definition: JDAQTriggeredHit.hh:25
KM3NETDAQ::JDAQSuperFrame::const_iterator
JDAQFrame::const_iterator const_iterator
Definition: JDAQSuperFrame.hh:34
debug
int debug
debug level
Definition: JSirene.cc:59
KM3NETDAQ::JDAQTriggerCounter::next
JTriggerCounter_t next()
Increment trigger counter.
Definition: JDAQTriggerCounter.hh:120
JSUPPORT::JLimit::getLimit
const JLimit & getLimit() const
Get limit.
Definition: JLimit.hh:73
KM3NETDAQ::JDAQSummaryslice
Data summary slice.
Definition: JDAQSummaryslice.hh:31
KM3NETDAQ::JDAQKeyHit
DAQ key hit.
Definition: JDAQKeyHit.hh:24
STATUS
#define STATUS(A)
Definition: JMessage.hh:63
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition: JParser.hh:1954
DEBUG
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
std
Definition: jaanetDictionary.h:36
KM3NETDAQ
KM3NeT DAQ data structures and auxiliaries.
Definition: DataQueue.cc:39
FATAL
#define FATAL(A)
Definition: JMessage.hh:67
outputFile
string outputFile
Definition: JDAQTimesliceSelector.cc:37
JROOT::counter_type
Long64_t counter_type
Type definition for counter.
Definition: JCounter.hh:24