Jpp
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
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  const JStatus status;
59 
60  const JPMTDefaultSimulator simulator(parameters, pmt);
61 
62 
63  TFile out(outputFile.c_str(), "recreate");
64 
65  JPMTData<JPMTSignal> input;
66  JPMTData<JPMTPulse> output;
67 
68 
69  TH1D h1("tts", NULL, 120, -20.0, +100.0);
70 
71  for (int i = 0; i != numberOfHits; ++i) {
72 
73  const double t_ns = 0.0;
74  const int npe = 1;
75  const double W = 1.0 / (double) numberOfHits;
76 
77  input .clear();
78  output.clear();
79 
80  input.push_back(JPMTSignal(t_ns, npe));
81 
82  simulator.processHits(pmt, calibration, status, input, output);
83 
84  for (JPMTData<JPMTPulse>::const_iterator hit = output.begin(); hit != output.end(); ++hit) {
85  h1.Fill(hit->t_ns, W);
86  }
87  }
88 
89  out.Write();
90  out.Close();
91 }
Utility class to parse command line options.
Definition: JParser.hh:1500
debug
Definition: JMessage.hh:29
Utility class to parse parameter values.
Definition: JProperties.hh:496
*fatal Wrong number of arguments esac JCookie sh typeset Z DETECTOR typeset Z SOURCE_RUN typeset Z TARGET_RUN set_variable PARAMETERS_FILE $WORKDIR parameters
Definition: diff-Tuna.sh:38
Simple data structure to support I/O of equations (see class JLANG::JEquation).
then for HISTOGRAM in h0 h1
Definition: JMatrixNZ.sh:69
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Definition: JParser.hh:66
string outputFile
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition: JParser.hh:1961
int debug
debug level
Definition: JSirene.cc:63
General purpose messaging.
#define FATAL(A)
Definition: JMessage.hh:67
Utility class to parse command line options.
int main(int argc, char *argv[])
Definition: Main.cpp:15