Jpp  15.0.1-rc.2-highQE
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JRoot.cc
Go to the documentation of this file.
1 
2 #include <string>
3 #include <iostream>
4 #include <fstream>
5 #include <iomanip>
6 #include <vector>
7 #include <limits>
8 #include <string.h>
9 
10 #include "TROOT.h"
11 #include "TFile.h"
12 
13 #include "JDAQ/JDAQEventIO.hh"
14 #include "JDAQ/JDAQTimesliceIO.hh"
16 
17 #include "JROOT/JTreeReader.hh"
18 #include "JIO/JStreamIO.hh"
19 #include "JIO/JByteArrayIO.hh"
20 #include "JSupport/JSupport.hh"
22 
23 #include "Jeep/JParser.hh"
24 #include "Jeep/JMessage.hh"
25 
26 
27 namespace JSUPPORT {
28 
29  /**
30  * Get key for given DAQ data type.
31  *
32  * \param type data type
33  * \return map element
34  */
35  template<>
36  template<class T>
38  {
39  return KM3NETDAQ::getDataType<T>();
40  }
41 }
42 
43 
44 /**
45  * \file
46  * Auxiliary program to convert binary data to ROOT formatted data.
47  * \author mdejong
48  */
49 int main(int argc, char **argv)
50 {
51  using namespace std;
52 
53  vector<string> inputFile;
54  string outputFile;
55  Long64_t numberOfEvents;
56  int debug;
57 
58  try {
59 
60  JParser<> zap("Auxiliary program to convert binary data to ROOT formatted data.");
61 
62  zap['f'] = make_field(inputFile);
63  zap['o'] = make_field(outputFile);
64  zap['n'] = make_field(numberOfEvents) = numeric_limits<Long64_t>::max();
65  zap['d'] = make_field(debug) = 1;
66 
67  zap(argc, argv);
68  }
69  catch(const exception &error) {
70  FATAL(error.what() << endl);
71  }
72 
73 
74  cout.tie(&cerr);
75 
76  using namespace KM3NETDAQ;
77  using namespace JPP;
78 
79 
80  JAutoTreeWriter<int> writer;
81 
82  writer.open(outputFile.c_str());
83 
84  if (!writer.is_open()) {
85  FATAL("Error opening file " << outputFile << endl);
86  }
87 
88  writer.insert<JDAQTypes_t>();
89 
90 
91  JDAQPreamble preamble;
92  vector<char> buffer;
93 
94  Long64_t event_count = 0;
95 
96  for (vector<string>::const_iterator i = inputFile.begin(); i != inputFile.end(); ++i) {
97 
98  STATUS("File: " << *i << endl);
99 
100  ifstream in(i->c_str());
101 
102  for (JStreamReader is(in); event_count < numberOfEvents; ++event_count) {
103 
104  STATUS("event: " << setw(10) << event_count << '\r'); DEBUG(endl);
105 
106  if (!(is >> preamble)) {
107  break;
108  }
109 
110  DEBUG("length " << preamble.getLength() << endl);
111  DEBUG("type " << preamble.getDataType() << endl);
112 
113  if (preamble.getLength() < (int) getSizeof<JDAQPreamble>()) {
114 
115  ERROR("File " << *i << " invalid length " << preamble.getLength() << endl);
116 
117  break;
118  }
119 
120  JAutoTreeWriter<int>::iterator i = writer.find(preamble.getDataType());
121 
122  if (i != writer.end()) {
123 
124  buffer.resize(preamble.getLength());
125 
126  memcpy(buffer.data(), static_cast<JDAQAbstractPreamble*>(&preamble), getSizeof<JDAQPreamble>());
127 
128  is.read(buffer.data() + getSizeof<JDAQPreamble>(),
129  preamble.getLength() - getSizeof<JDAQPreamble>());
130 
131  JByteArrayReader bin(buffer.data(), buffer.size());
132 
133  i->second->copy(bin);
134 
135  } else {
136 
137  in.ignore((streamsize) (preamble.getLength() - getSizeof<JDAQPreamble>()));
138 
139  ERROR("Illegal data type, skip" << endl);
140  }
141  }
142  STATUS(endl);
143 
144  in.close();
145  }
146 
147  writer.close();
148 }
Utility class to parse command line options.
Definition: JParser.hh:1500
int main(int argc, char *argv[])
Definition: Main.cc:15
ROOT TTree parameter settings of various packages.
Auxiliary class to copy input data to corresponding TTree.
#define STATUS(A)
Definition: JMessage.hh:63
Auxiliary class for a type holder.
Definition: JType.hh:19
TTree reading for template data type.
string outputFile
is
Definition: JDAQCHSM.chsm:167
Type list.
Definition: JTypeList.hh:22
static JKey_t getKey(JType< T > type)
Get key.
Binary input based on std::istream.
Definition: JStreamIO.hh:24
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition: JParser.hh:1961
Byte array binary input.
Definition: JByteArrayIO.hh:25
#define ERROR(A)
Definition: JMessage.hh:66
map_type::iterator iterator
int debug
debug level
Definition: JSirene.cc:63
Simple data structure for the DAQ preamble required for a correct calculation of the object size for ...
virtual void open(const char *file_name) override
Open file.
General purpose messaging.
#define FATAL(A)
Definition: JMessage.hh:67
Utility class to parse command line options.
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
size_t getSizeof< JDAQPreamble >()
Get size of type.
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2] eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY source JAcoustics sh $DETECTOR_ID CHECK_EXIT_CODE typeset A TRIPODS get_tripods $WORKDIR tripod txt TRIPODS for EMITTER in
Definition: JCanberra.sh:41