Jpp
15.0.3
the software that should make you happy
Main Page
Namespaces
Classes
Files
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Pages
examples
JDetector
JToTProbability.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
10
#include "
JDetector/JPMTParameters.hh
"
11
#include "
JDetector/JPMTAnalogueSignalProcessor.hh
"
12
13
#include "
Jeep/JPrint.hh
"
14
#include "
Jeep/JParser.hh
"
15
#include "
Jeep/JMessage.hh
"
16
17
18
/**
19
* \file
20
*
21
* Example program to histogram time-over-threshold probability distributions.
22
* \author mdejong
23
*/
24
int
main
(
int
argc,
char
**argv)
25
{
26
using namespace
std;
27
using namespace
JPP;
28
29
string
outputFile
;
30
JPMTParameters
parameters
;
31
double
NPE;
32
int
debug
;
33
34
try
{
35
36
JProperties
properties =
parameters
.getProperties();
37
38
JParser<>
zap(
"Example program to histogram time-over-threshold probability distributions."
);
39
40
zap[
'o'
] =
make_field
(
outputFile
) =
"time-over-threshold.root"
;
41
zap[
'P'
] =
make_field
(properties) =
JPARSER::initialised
();
42
zap[
'N'
] =
make_field
(NPE);
43
zap[
'd'
] =
make_field
(
debug
) = 0;
44
45
zap(argc, argv);
46
47
}
48
catch
(
const
exception &error) {
49
FATAL
(error.what() << endl);
50
}
51
52
const
JPMTAnalogueSignalProcessor cpu(
parameters
);
53
54
TFile out(
outputFile
.c_str(),
"recreate"
);
55
56
57
const
double
xmin = 0.5;
58
const
double
xmax = 255.5;
59
const
int
nx = (int) ((xmax - xmin) / 1.0);
60
61
TH1D h0(
"h0"
, NULL, nx, xmin, xmax);
62
63
for
(
int
i = 1; i <= h0.GetNbinsX(); ++i) {
64
65
const
double
x
= h0.GetBinCenter(i);
66
const
double
y =
getTimeOverThresholdProbability
(cpu, x, NPE);
67
68
h0.SetBinContent(i, y);
69
}
70
71
out.Write();
72
out.Close();
73
}
JPARSER::JParser
Utility class to parse command line options.
Definition:
JParser.hh:1500
main
int main(int argc, char *argv[])
Definition:
Main.cc:15
JDETECTOR::getTimeOverThresholdProbability
double getTimeOverThresholdProbability(const JPMTAnalogueSignalProcessor &pmt, const double tot_ns, const double NPE, const double precision=1.0e-4)
Get time-over-threshold probability.
Definition:
JPMTAnalogueSignalProcessor.hh:881
JEEP::JProperties
Utility class to parse parameter values.
Definition:
JProperties.hh:496
parameters
*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
JPARSER::initialised
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Definition:
JParser.hh:66
outputFile
string outputFile
Definition:
JDAQTimesliceSelector.cc:37
x
then chmod x
Definition:
JEvtReweightMupageParameterScan.sh:413
JPrint.hh
I/O formatting auxiliaries.
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition:
JParser.hh:1961
debug
int debug
debug level
Definition:
JSirene.cc:63
JPMTParameters.hh
JMessage.hh
General purpose messaging.
FATAL
#define FATAL(A)
Definition:
JMessage.hh:67
JParser.hh
Utility class to parse command line options.
JPMTAnalogueSignalProcessor.hh
PMT analogue signal processor.
Generated by
1.8.5