Jpp
Functions
JToT.cc File Reference
#include <string>
#include <iostream>
#include <iomanip>
#include "TROOT.h"
#include "TFile.h"
#include "TH1D.h"
#include "JDetector/JPMTParameters.hh"
#include "JDetector/JPMTAnalogueSignalProcessor.hh"
#include "JDetector/JPMTDefaultSimulator.hh"
#include "Jeep/JPrint.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 histogram time over threshold distributions.

Author
mdejong

Definition in file JToT.cc.

Function Documentation

◆ main()

int main ( int  argc,
char **  argv 
)

Definition at line 25 of file JToT.cc.

26 {
27  using namespace std;
28  using namespace JPP;
29 
30  string outputFile;
31  JPMTParameters parameters;
32  int NPE;
33  int numberOfHits;
34  double precision;
35  int debug;
36 
37  try {
38 
39  JProperties properties = parameters.getProperties();
40 
41  JParser<> zap("Example program to histogram time over threshold distributions.");
42 
43  zap['o'] = make_field(outputFile) = "tot.root";
44  zap['P'] = make_field(properties) = JPARSER::initialised();
45  zap['N'] = make_field(NPE) = 1;
46  zap['n'] = make_field(numberOfHits) = 1000000;
47  zap['e'] = make_field(precision) = 0.005;
48  zap['d'] = make_field(debug) = 0;
49 
50  zap(argc, argv);
51  }
52  catch(const exception &error) {
53  FATAL(error.what() << endl);
54  }
55 
56 
57  if (debug >= JEEP::debug_t) {
58  cout << "PMT parameters:" << endl;
59  cout << parameters.getProperties(JEquationParameters("=", "\n", "", "")) << endl;
60  }
61 
62  const JPMTAnalogueSignalProcessor cpu(parameters);
63 
64  const JPMTIdentifier pmt(1,0);
65  const JCalibration calibration;
66 
67  const JPMTDefaultSimulator simulator(parameters, pmt);
68 
69  JPMTData<JPMTSignal> input;
70  JPMTData<JPMTPulse> output;
71 
72  ASSERT(numberOfHits > 0);
73 
74 
75  TFile out(outputFile.c_str(), "recreate");
76 
77  const double xmin = -0.5;
78  const double xmax = 255.5;
79  const double dx = 1.0;
80  const int nx = (int) ((xmax - xmin) / dx);
81 
82  TH1D h0("h0", NULL, nx, xmin, xmax);
83  TH1D h1("h1", NULL, nx, xmin, xmax);
84  TH1D h2("h2", NULL, nx, xmin, xmax);
85 
86  h1.Sumw2();
87  h2.Sumw2();
88 
89  DEBUG(" ToT npe dP/dnpe dP/dToT" << endl);
90 
91  for (int i = 1; i <= h0.GetNbinsX(); ++i) {
92 
93  const double x = h0.GetBinCenter(i);
94  const double npe = cpu.getNPE(x);
95  const double y = cpu.getProbability(npe, NPE);
96  const double v = cpu.getDerivative(npe);
97 
98  DEBUG(" "
99  << FIXED(5,1) << x << " "
100  << FIXED(7,3) << npe << " "
101  << FIXED(5,3) << y << " "
102  << FIXED(6,4) << y*v << endl);
103 
104  h0.SetBinContent(i, y*v);
105  }
106 
107 
108  if (numberOfHits > 0) {
109 
110  int number_of_hits = 0;
111 
112  for (int i = 0; i != numberOfHits; ++i) {
113 
114  const double npe = cpu.getRandomAmplitude(NPE);
115 
116  if (cpu.applyThreshold(npe)) {
117 
118  const double tot_ns = cpu.getTimeOverThreshold(npe);
119 
120  ++number_of_hits;
121 
122  h1.Fill(tot_ns);
123  }
124  }
125 
126  h1.Scale(1.0 / (double) number_of_hits / dx);
127  }
128 
129 
130  if (numberOfHits > 0) {
131 
132  int number_of_hits = 0;
133 
134  const double t_ns = 25.0;
135 
136  for (int i = 0; i != numberOfHits; ++i) {
137 
138  input .clear();
139  output.clear();
140 
141  input.push_back(JPMTSignal(t_ns, NPE));
142 
143  simulator.processHits(pmt, calibration, input, output);
144 
145  for (JPMTData<JPMTPulse>::const_iterator hit = output.begin(); hit != output.end(); ++hit) {
146 
147  ++number_of_hits;
148 
149  h2.Fill(hit->tot_ns);
150  }
151  }
152 
153  h2.Scale(1.0 / (double) number_of_hits / dx);
154  }
155 
156  out.Write();
157  out.Close();
158 
159  for (int i = 1; i <= h0.GetNbinsX(); ++i) {
160 
161  const Double_t x = h0.GetBinCenter (i);
162  const Double_t y0 = h0.GetBinContent(i);
163  const Double_t y1 = h1.GetBinContent(i);
164  const Double_t y2 = h2.GetBinContent(i);
165 
166  DEBUG("[" << setw(3) << i << "]" << ' '
167  << FIXED(5,1) << x << ' '
168  << FIXED(6,4) << y0 << ' '
169  << FIXED(6,4) << y1 << ' '
170  << FIXED(6,4) << y2 << endl);
171 
172  ASSERT(fabs(y0 - y1) < precision);
173  ASSERT(fabs(y0 - y2) < precision);
174  }
175 }
FIXED
Auxiliary data structure for floating point format specification.
Definition: JPrint.hh:481
ASSERT
#define ASSERT(A,...)
Assert macro.
Definition: JMessage.hh:90
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
JTOOLS::v
data_type v[N+1][M+1]
Definition: JPolint.hh:707
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
DEBUG
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
std
Definition: jaanetDictionary.h:36
FATAL
#define FATAL(A)
Definition: JMessage.hh:67
outputFile
string outputFile
Definition: JDAQTimesliceSelector.cc:37
JEEP::debug_t
debug
Definition: JMessage.hh:29