Jpp
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Functions
JToT.cc File Reference

Example program to histogram time over threshold distributions. More...

#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

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  for (int i = 0; i != numberOfHits; ++i) {
111 
112  const double npe = cpu.getRandomAmplitude(NPE);
113 
114  if (cpu.applyThreshold(npe)) {
115 
116  const double tot_ns = cpu.getTimeOverThreshold(npe);
117 
118  h1.Fill(tot_ns);
119  }
120  }
121 
122  h1.Scale(1.0 / (double) numberOfHits / dx);
123  }
124 
125 
126  if (numberOfHits > 0) {
127 
128  const double t_ns = 25.0;
129 
130  for (int i = 0; i != numberOfHits; ++i) {
131 
132  input .clear();
133  output.clear();
134 
135  input.push_back(JPMTSignal(t_ns, NPE));
136 
137  simulator.processHits(pmt, calibration, input, output);
138 
139  for (JPMTData<JPMTPulse>::const_iterator hit = output.begin(); hit != output.end(); ++hit) {
140  h2.Fill(hit->tot_ns);
141  }
142  }
143 
144  h2.Scale(1.0 / (double) numberOfHits / dx);
145  }
146 
147  for (int i = 1; i <= h0.GetNbinsX(); ++i) {
148 
149  const Double_t x = h0.GetBinCenter (i);
150  const Double_t y0 = h0.GetBinContent(i);
151  const Double_t y1 = h1.GetBinContent(i);
152  const Double_t y2 = h2.GetBinContent(i);
153 
154  DEBUG("[" << setw(3) << i << "]" << ' '
155  << FIXED(5,1) << x << ' '
156  << FIXED(6,4) << y0 << ' '
157  << FIXED(6,4) << y1 << ' '
158  << FIXED(6,4) << y2 << endl);
159 
160  ASSERT(fabs(y0 - y1) < precision);
161  ASSERT(fabs(y0 - y2) < precision);
162  }
163 
164 
165  out.Write();
166  out.Close();
167 }
Utility class to parse command line options.
Definition: JParser.hh:1410
debug
Definition: JMessage.hh:27
Utility class to parse parameter values.
Definition: JProperties.hh:484
Simple data structure to support I/O of equations (see class JLANG::JEquation).
Empty structure for specification of parser element that is initialised (i.e.
Definition: JParser.hh:64
Auxiliary data structure for floating point format specification.
Definition: JPrint.hh:461
string outputFile
#define ASSERT(A)
Assert macro.
Definition: JMessage.hh:72
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition: JParser.hh:1836
int debug
debug level
Definition: JSirene.cc:59
#define FATAL(A)
Definition: JMessage.hh:65
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:60