Jpp
JTransitTimeGenerator.cc
Go to the documentation of this file.
1 
2 #include <string>
3 #include <iostream>
4 #include <iomanip>
5 
6 #include "TROOT.h"
7 #include "TFile.h"
8 #include "TH1D.h"
9 
12 
13 #include "Jeep/JParser.hh"
14 #include "Jeep/JMessage.hh"
15 
16 
17 /**
18  * \file
19  *
20  * Example program to test generation of transit times using PMT simulator.
21  * \author mdejong
22  */
23 int main(int argc, char **argv)
24 {
25  using namespace std;
26  using namespace JPP;
27 
28  string outputFile;
29  int numberOfHits;
30  JPMTParameters parameters;
31  int debug;
32 
33  try {
34 
35  JProperties properties = parameters.getProperties();
36 
37  JParser<> zap("Example program to test generation of transit times using PMT simulator.");
38 
39  zap['o'] = make_field(outputFile) = "pmt.root";
40  zap['n'] = make_field(numberOfHits) = 10000;
41  zap['P'] = make_field(properties) = JPARSER::initialised();
42  zap['d'] = make_field(debug) = 0;
43 
44  zap(argc, argv);
45  }
46  catch(const exception &error) {
47  FATAL(error.what() << endl);
48  }
49 
50 
51  if (debug >= JEEP::debug_t) {
52  cout << "PMT parameters:" << endl;
53  cout << parameters.getProperties(JEquationParameters("=", "\n", "", "")) << endl;
54  }
55 
56  const JPMTIdentifier pmt(1,0);
57  const JCalibration calibration;
58 
59  const JPMTDefaultSimulator simulator(parameters, pmt);
60 
61 
62  TFile out(outputFile.c_str(), "recreate");
63 
64  JPMTData<JPMTSignal> input;
65  JPMTData<JPMTPulse> output;
66 
67 
68  TH1D h1("tts", NULL, 120, -20.0, +100.0);
69 
70  for (int i = 0; i != numberOfHits; ++i) {
71 
72  const double t_ns = 0.0;
73  const int npe = 1;
74  const double W = 1.0 / (double) numberOfHits;
75 
76  input .clear();
77  output.clear();
78 
79  input.push_back(JPMTSignal(t_ns, npe));
80 
81  simulator.processHits(pmt, calibration, input, output);
82 
83  for (JPMTData<JPMTPulse>::const_iterator hit = output.begin(); hit != output.end(); ++hit) {
84  h1.Fill(hit->t_ns, W);
85  }
86  }
87 
88  out.Write();
89  out.Close();
90 }
main
int main(int argc, char **argv)
Definition: JTransitTimeGenerator.cc:23
JPMTSimulator.hh
JMessage.hh
JPARSER::initialised
Empty structure for specification of parser element that is initialised (i.e.
Definition: JParser.hh:63
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
JParser.hh
JLANG::JEquationParameters
Simple data structure to support I/O of equations (see class JLANG::JEquation).
Definition: JEquationParameters.hh:20
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition: JParser.hh:1954
JEEP::JProperties
Utility class to parse parameter values.
Definition: JProperties.hh:496
std
Definition: jaanetDictionary.h:36
JPMTDefaultSimulator.hh
FATAL
#define FATAL(A)
Definition: JMessage.hh:67
outputFile
string outputFile
Definition: JDAQTimesliceSelector.cc:37
JEEP::debug_t
debug
Definition: JMessage.hh:29