Jpp test-rotations-new
the software that should make you happy
Loading...
Searching...
No Matches
JCDFTable1D.cc
Go to the documentation of this file.
1
2#include <string>
3#include <iostream>
4#include <fstream>
5#include <iomanip>
6#include <vector>
7
8#include "TROOT.h"
9#include "TFile.h"
10#include "TH1D.h"
11
14#include "JPhysics/JCDFTable.hh"
15#include "JPhysics/JPDFTypes.hh"
18#include "Jeep/JProperties.hh"
19#include "Jeep/JPrint.hh"
20#include "Jeep/JParser.hh"
21#include "Jeep/JMessage.hh"
22
23
24/**
25 * \file
26 *
27 * Program to plot CDF of Cherenkov light from muon using interpolation tables.
28 * \author mdejong
29 */
30int main(int argc, char **argv)
31{
32 using namespace std;
33 using namespace JPP;
34
35 string inputFile;
36 string outputFile;
37 int numberOfBins = 200; //!< number of bins for average CDF integral of optical module
38 double safetyFactor = 1.7; //!< safety factor for average CDF integral of optical module
39 int debug;
40
41 try {
42
43 JProperties properties;
44
45 properties.insert(gmake_property(numberOfBins));
46 properties.insert(gmake_property(safetyFactor));
47
48 JParser<> zap("Program to plot PDF of Cherenkov light from muon using interpolation tables.");
49
50 zap['@'] = make_field(properties) = JPARSER::initialised();
51 zap['f'] = make_field(inputFile);
52 zap['o'] = make_field(outputFile) = "cdf.root";
53 zap['d'] = make_field(debug) = 0;
54
55 zap(argc, argv);
56 }
57 catch(const exception &error) {
58 FATAL(error.what() << endl);
59 }
60
61
62 typedef JHermiteSplineFunction1D_t JFunction1D_t;
65 JPolint1FunctionalGridMap>::maplist JMapList_t;
67
69
70 JCDF_t function;
71
72 try {
73
74 NOTICE("loading input from file " << inputFile << "... " << flush);
75
76 function.load(inputFile.c_str());
77
78 NOTICE("OK" << endl);
79 }
80 catch(const JException& error) {
81 FATAL(error.what() << endl);
82 }
83
84 JCDF1D_t integral(function, numberOfBins, safetyFactor);
85
86
87 TFile out(outputFile.c_str(), "recreate");
88
89 TH1D h0("h0", NULL,
90 1000, function.intensity.getXmin(), function.intensity.getXmax());
91
92 for (int i = 1; i <= h0.GetNbinsX(); ++i) {
93
94 const double R = h0.GetBinCenter(i);
95 const double y = integral.getNPE(R);
96
97 h0.SetBinContent(i, y);
98 }
99
100 out.Write();
101 out.Close();
102}
int main(int argc, char **argv)
string outputFile
Various implementations of functional maps.
General purpose messaging.
#define NOTICE(A)
Definition JMessage.hh:64
#define FATAL(A)
Definition JMessage.hh:67
int debug
debug level
Definition JSirene.cc:72
Auxiliary methods for PDF calculations.
Numbering scheme for PDF types.
Utility class to parse command line options.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition JParser.hh:2142
I/O formatting auxiliaries.
Utility class to parse parameter values.
#define gmake_property(A)
macros to convert (template) parameter to JPropertiesElement object
int numberOfBins
number of bins for average CDF integral of optical module
Definition JSirene.cc:73
double safetyFactor
safety factor for average CDF integral of optical module
Definition JSirene.cc:74
Utility class to parse parameter values.
General exception.
Definition JException.hh:24
Utility class to parse command line options.
Definition JParser.hh:1698
Custom class for CDF table in 1 dimension.
Multi-dimensional CDF table for arrival time of Cherenkov light.
Definition JCDFTable.hh:58
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Definition JParser.hh:68
Type definition of a spline interpolation method based on a JCollection with double result type.
Auxiliary class for recursive map list generation.
Definition JMapList.hh:109
Type definition of a 1st degree polynomial interpolation based on a JGridMap implementation.
Type definition of a 1st degree polynomial interpolation based on a JMap implementation.