Jpp
Functions
JAAnetPreprocessor.cc File Reference
#include "evt/Head.hh"
#include "evt/Evt.hh"
#include "JAAnet/JAAnetToolkit.hh"
#include "JSupport/JMultipleFileScanner.hh"
#include "JSupport/JFileRecorder.hh"
#include "JSupport/JMonteCarloFileSupportkit.hh"
#include "JSupport/JSupport.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
#include "JSupport/JMeta.hh"

Go to the source code of this file.

Functions

int main (int argc, char **argv)
 

Detailed Description

Auxiliary program to preprocess an aanet file.
Currently it provide support for the suppression of empty events from a file, but implementation of other options is encouraged.

Author
mlincetto

Definition in file JAAnetPreprocessor.cc.

Function Documentation

◆ main()

int main ( int  argc,
char **  argv 
)

Definition at line 23 of file JAAnetPreprocessor.cc.

24 {
25  using namespace std;
26  using namespace JPP;
27 
28  JMultipleFileScanner<JAAnetTypes_t> inputFile;
29  JFileRecorder <JAAnetTypes_t> outputFile;
30  JLimit_t& numberOfEvents = inputFile.getLimit();
31  int debug;
32  bool filterEmptyMCH;
33  double timeOffset_ms;
34 
35  try {
36 
37  JParser<> zap("Auxiliary program to preprocess aanet files.");
38 
39  zap['f'] = make_field(inputFile);
40  zap['o'] = make_field(outputFile);
41  zap['n'] = make_field(numberOfEvents) = JLimit::max();
42  zap['d'] = make_field(debug) = 1;
43  zap['H'] = make_field(filterEmptyMCH, "Filter events without hits.");
44  zap['T'] = make_field(timeOffset_ms, "Add time offset [ms] to events.") = 0;
45 
46  zap(argc, argv);
47  }
48  catch(const exception& error) {
49  FATAL(error.what() << endl);
50  }
51 
52  outputFile.open();
53 
54  // put header
55  Head header = inputFile.getHeader();
56 
57  JMultipleFileScanner<Evt> in = inputFile;
58 
59  outputFile.put(header);
60 
61  while (in.hasNext()) {
62 
63  Evt* cur = in.next();
64 
65  if (filterEmptyMCH and cur->mc_hits.empty()) { continue; }
66 
67  if (timeOffset_ms != 0) { add_time_offset (*cur, 1.0e6 * timeOffset_ms); }
68 
69  outputFile.put(*cur);
70 
71  } // while
72 
73  outputFile.close();
74 }
JSUPPORT::JLimit_t
JLimit JLimit_t
Type definition of limit.
Definition: JLimit.hh:215
JPARSER::JParser
Utility class to parse command line options.
Definition: JParser.hh:1493
JPP
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Definition: JAAnetToolkit.hh:37
debug
int debug
debug level
Definition: JSirene.cc:59
JSUPPORT::JLimit::getLimit
const JLimit & getLimit() const
Get limit.
Definition: JLimit.hh:73
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition: JParser.hh:1954
std
Definition: jaanetDictionary.h:36
JAANET::add_time_offset
void add_time_offset(Evt &evt, const double tOff)
Add time offset to mc event; according to current definition, the absolute time of the event is defin...
Definition: JAAnetToolkit.hh:653
FATAL
#define FATAL(A)
Definition: JMessage.hh:67
outputFile
string outputFile
Definition: JDAQTimesliceSelector.cc:37