Jpp 20.0.0-195-g190c9e876
the software that should make you happy
Loading...
Searching...
No Matches
Functions
JFitK40.cc File Reference

Auxiliary program to fit PMT parameters from JMergeCalibrateK40.cc output. More...

#include <string>
#include <iostream>
#include <iomanip>
#include <cmath>
#include <set>
#include <vector>
#include <algorithm>
#include <memory>
#include "TROOT.h"
#include "TFile.h"
#include "TH1D.h"
#include "TH2D.h"
#include "km3net-dataformat/online/JDAQ.hh"
#include "km3net-dataformat/online/JDAQHeader.hh"
#include "JTools/JConstants.hh"
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JPMTParametersMap.hh"
#include "JDetector/JStringRouter.hh"
#include "JROOT/JRootFileWriter.hh"
#include "JROOT/JRootToolkit.hh"
#include "JTools/JRange.hh"
#include "JLang/JSTDObjectWriter.hh"
#include "JSupport/JMeta.hh"
#include "JSupport/JSingleFileScanner.hh"
#include "JCalibrate/JCalibrateK40.hh"
#include "JCalibrate/JTDC_t.hh"
#include "JFitK40.hh"
#include "Jeep/JProperties.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

Auxiliary program to fit PMT parameters from JMergeCalibrateK40.cc output.


Note that the contebts of the 2D histograms are converted to standard deviations.

Author
mdejong

Definition in file JFitK40.cc.

Function Documentation

◆ main()

int main ( int argc,
char ** argv )

Definition at line 64 of file JFitK40.cc.

65{
66 using namespace std;
67 using namespace JPP;
68 using namespace KM3NETDAQ;
69
70 typedef JRange<double> JRange_t;
71
72 JK40Parameters K40 = JK40Parameters::getInstance();
73
74 string inputFile;
75 string outputFile;
76 string detectorFile;
77 string pmtFile;
78 JTDC_t TDC;
79 bool reverse;
80 bool overwriteDetector;
81 JCounter writeFits;
82 bool fitAngle;
83 bool fitNoise;
84 bool fitModel;
85 bool fixQE;
86 JRange_t X;
87 JRange_t Y;
88 int debug;
89
90 try {
91
92 JProperties properties;
93
94 properties.insert(gmake_property(MINIMAL_RATE_HZ));
95 properties.insert(gmake_property(STDEV));
96 properties.insert(gmake_property(MAXIMAL_COUNTS));
97 properties.insert(gmake_property(QE_MIN));
98 properties.insert(gmake_property(T0_NS));
99 properties.insert(gmake_property(K40.R));
100 properties.insert(gmake_property(K40.p1));
101 properties.insert(gmake_property(K40.p2));
102 properties.insert(gmake_property(K40.p3));
103 properties.insert(gmake_property(K40.p4));
104 properties.insert(gmake_property(TEROSTAT_R1));
105 properties.insert(gmake_property(TEROSTAT_DZ));
106 properties.insert(gmake_property(BELL_SHAPE));
107 properties.insert(gmake_property(MESTIMATOR));
108
109 JParser<> zap("Auxiliary program to fit PMT parameters from JMergeCalibrateK40 output.");
110
111 zap['@'] = make_field(properties) = JPARSER::initialised();
112 zap['f'] = make_field(inputFile, "input file (output from JMergeCalibrateK40).");
113 zap['o'] = make_field(outputFile, "output file.") = "fit.root";
114 zap['a'] = make_field(detectorFile, "detector file.");
115 zap['P'] = make_field(pmtFile, "specify PMT file name that can be given to JTriggerEfficiency.") = "";
116 zap['!'] = make_field(TDC,
117 "Fix time offset(s) of PMT(s) of certain module(s), e.g."
118 "\n-! \"808969848 0 808982077 23\" will fix time offsets of PMT 0 of module 808969848 and of PMT 23 of module 808982077."
119 "\nSame PMT offset can be fixed for all optical modules, e.g."
120 "\n-! \"-1 0 -1 23\" will fix time offsets of PMTs 0 and 23 of all optical modules.") = JPARSER::initialised();
121 zap['r'] = make_field(reverse, "reverse TDC constraints due to option -! <TDC>.");
122 zap['A'] = make_field(overwriteDetector, "overwrite detector file provided through '-a' with fitted time offsets.");
123 zap['w'] = make_field(writeFits, "write fit results to ROOT file; -ww also write fitted TTS to PMT file.");
124 zap['D'] = make_field(fitAngle, "fit angular distribution; fix normalisation.");
125 zap['B'] = make_field(fitNoise, "fit background.");
126 zap['M'] = make_field(fitModel, "fit angular distribution as well as normalisation; fix PMT QEs = 1.0.");
127 zap['Q'] = make_field(fixQE, "fix PMT QEs = 1.0.");
128 zap['x'] = make_field(X, "fit range (PMT pairs).") = JRange_t();
129 zap['y'] = make_field(Y, "fit range (time residual [ns]).") = JRange_t();
130 zap['d'] = make_field(debug, "debug.") = 1;
131
132 zap(argc, argv);
133 }
134 catch(const exception &error) {
135 FATAL(error.what() << endl);
136 }
137
138
139 if ((fitModel ? 1 : 0) +
140 (fitAngle ? 1 : 0) +
141 (fitNoise ? 1 : 0) +
142 (fixQE ? 1 : 0) > 1) {
143 FATAL("Use either option -M, -D, -B or -Q" << endl);
144 }
145
146 const int option = (fitModel ? FIT_MODEL_t :
147 fitAngle ? FIT_PMTS_AND_ANGULAR_DEPENDENCE_t :
148 fitNoise ? FIT_PMTS_AND_BACKGROUND_t :
149 fixQE ? FIT_PMTS_QE_FIXED_t :
150 FIT_PMTS_t);
151
152 if (reverse) {
153 TDC.reverse();
154 }
155
156 for (JTDC_t::const_iterator i = TDC.begin(); i != TDC.end(); ++i) {
157 DEBUG("PMT " << setw(10) << i->first << ' ' << setw(2) << i->second << " constrain t0." << endl);
158 }
159
160 try {
161 TDC.is_valid(true);
162 }
163 catch(const exception &error) {
164 FATAL(error.what() << endl);
165 }
166
167 JDetector detector;
168
169 try {
170 load(detectorFile, detector);
171 }
172 catch(const JException& error) {
173 FATAL(error);
174 }
175
176
177 JPMTParametersMap parameters;
178
179 if (pmtFile != "") {
180 try {
181 parameters.load(pmtFile.c_str());
182 }
183 catch(const exception& error) {}
184 }
185
186 parameters.comment.add(JMeta(argc, argv));
187
188
189 TFile* in = TFile::Open(inputFile.c_str(), "exist");
190
191 if (in == NULL || !in->IsOpen()) {
192 FATAL("File: " << inputFile << " not opened." << endl);
193 }
194
195
196 TFile out(outputFile.c_str(), "recreate");
197
198 TH1D h0("chi2", NULL, 500, 0.0, 5.0);
199 TH1D hn("hn", NULL, 501, -0.5, 500.0);
200 TH1D hr("rate", NULL, 500, 0.0, 25.0);
201 TH1D h1("p1", NULL, 500, -5.0, +5.0);
202 TH1D h2("p2", NULL, 500, -5.0, +5.0);
203 TH1D h3("p3", NULL, 500, -5.0, +5.0);
204 TH1D h4("p4", NULL, 500, -5.0, +5.0);
205 TH1D hc("cc", NULL, 500, -0.1, +0.1);
206 TH1D hb("bc", NULL, 500, -0.1, +0.1);
207
208 const floor_range range = getRangeOfFloors(detector);
209 const JStringRouter string(detector);
210
211 TH2D H2("detector", NULL,
212 string.size() + 0, -0.5, string.size() - 0.5,
213 range.getUpperLimit(), 1 - 0.5, range.getUpperLimit() + 0.5);
214
215 for (Int_t i = 1; i <= H2.GetXaxis()->GetNbins(); ++i) {
216 H2.GetXaxis()->SetBinLabel(i, MAKE_CSTRING(string.at(i-1)));
217 }
218 for (Int_t i = 1; i <= H2.GetYaxis()->GetNbins(); ++i) {
219 H2.GetYaxis()->SetBinLabel(i, MAKE_CSTRING(i));
220 }
221
222 TH2D* HN = (TH2D*) H2.Clone("iterations");
223
224 JFit fit(MESTIMATOR, debug);
225
226 for (JDetector::iterator module = detector.begin(); module != detector.end(); ++module) {
227
228 if (module->getFloor() == 0) {
229 continue;
230 }
231
232 const JTDC_t::range_type range = TDC.equal_range(module->getID());
233
234 NOTICE("Module " << setw(10) << module->getID() << ' ' << getLabel(module->getLocation()) << " !" << distance(range.first, range.second) << endl);
235
236 TH2D* h2d = (TH2D*) in->Get(MAKE_CSTRING(module->getID() << _2R));
237
238 if (h2d == NULL || h2d->GetEntries() == 0) {
239
240 NOTICE("No data for module " << module->getID() << " -> set QEs to 0." << endl);
241
242 for (int pmt = 0; pmt != NUMBER_OF_PMTS; ++pmt) {
243 parameters[JPMTIdentifier(module->getID(), pmt)].QE = 0.0;
244 }
245
246 continue;
247 }
248
249 JModel model(*module, K40, range, option);
250
251 data_type data; // input data
252
253 vector<size_t> count[2] = {
254 vector<size_t>(NUMBER_OF_PMTS, 0), // counter (exclude self)
255 vector<size_t>(NUMBER_OF_PMTS, 1) // counter (include self)
256 };
257
258 for (int ix = 1; ix <= h2d->GetXaxis()->GetNbins(); ++ix) {
259
260 const pair_type pair = model.getPair(ix - 1);
261
262 auto& buffer = data[pair]; // storage
263
264 double V = 0.0; // integrated value
265 double W = 0.0; // integrated error
266
267 for (int iy = 1; iy <= h2d->GetYaxis()->GetNbins(); ++iy) {
268
269 const double x = h2d->GetXaxis()->GetBinCenter(ix);
270 const double y = h2d->GetYaxis()->GetBinCenter(iy);
271
272 if (X(x) && Y(y)) {
273
274 double value = h2d->GetBinContent(ix,iy);
275 double error = h2d->GetBinError (ix,iy);
276
277 buffer.push_back(rate_type(y, value, error));
278
279 double width = h2d->GetYaxis()->GetBinWidth(iy);
280
281 value *= width;
282 error *= width;
283
284 V += value;
285 W += error * error;
286 }
287 }
288
289 W = sqrt(W);
290
291 if (V <= 0.0 - STDEV*W) {
292 count[0][pair.first] += 1;
293 count[0][pair.second] += 1;
294 }
295
296 if (V <= MINIMAL_RATE_HZ + STDEV*W) {
297 count[1][pair.first] += 1;
298 count[1][pair.second] += 1;
299 }
300 }
301
302 for (int pmt = 0; pmt != NUMBER_OF_PMTS; ++pmt) {
303
304 if (count[0][pmt] >= MAXIMAL_COUNTS) { // too many paired rates negative
305
306 WARNING("PMT " << setw(10) << module->getID() << '.' << FILL(2,'0') << pmt << FILL() << " some rates negative -> fit background" << endl);
307
308 if (fit.value.parameters[pmt].status) {
309 model.parameters[pmt].bg.set();
310 }
311 }
312
313 if (count[1][pmt] == NUMBER_OF_PMTS) { // all paired rates too low
314
315 WARNING("PMT " << setw(10) << module->getID() << '.' << FILL(2,'0') << pmt << FILL() << " all rates to low -> disable" << endl);
316
317 model.parameters[pmt].disable();
318 }
319 }
320
321 DEBUG("Start value:" << endl << model << endl);
322
323 try {
324
325 fit.value = model; // start value
326
327 auto result = fit(data);
328
329 if (result.ndf <= 0) {
330
331 ERROR("Fit result " << setw(10) << module->getID() << " NDF " << setw(5) << result.ndf << " -> skip" << endl);
332
333 continue;
334 }
335
336 bool refit = false;
337
338 for (int pmt = 0; pmt != NUMBER_OF_PMTS; ++pmt) {
339
340 if (fit.value.parameters[pmt].status) {
341
342 if (fit.value.parameters[pmt].QE() <= QE_MIN ||
343 fit.value.parameters[pmt].QE() <= 0.0 + STDEV * fit.error.parameters[pmt].QE()) {
344
345 WARNING("PMT " << setw(10) << module->getID() << '.' << FILL(2,'0') << pmt << FILL() << ' '
346 << "QE = "
347 << FIXED(5,3) << fit.value.parameters[pmt].QE() << " +/- "
348 << FIXED(5,3) << fit.error.parameters[pmt].QE() << " "
349 << " -> disable" << (!refit ? " and refit" : "") << endl);
350
351 fit.value.parameters[pmt].disable();
352
353 refit = true;
354 }
355
356 if (fit.value.parameters[pmt].t0.atLimit(T0_NS)) {
357
358 WARNING("PMT " << setw(10) << module->getID() << '.' << FILL(2,'0') << pmt << FILL() << ' '
359 << "t0 at limit "
360 << FIXED(5,3) << fit.value.parameters[pmt].t0() << " +/- "
361 << FIXED(5,3) << fit.error.parameters[pmt].t0() << " "
362 << " -> disable" << (!refit ? " and refit" : "") << endl);
363
364 fit.value.parameters[pmt].disable();
365
366 refit = true;
367 }
368 }
369 }
370
371 if (refit) {
372
373 for (int pmt = 0; pmt != NUMBER_OF_PMTS; ++pmt) {
374 if (fit.value.parameters[pmt].status) {
375 fit.value.parameters[pmt].set(JPMTParameters_t::getInstance());
376 }
377 }
378
379 refit = false;
380 result = fit(data);
381 }
382
383 NOTICE("Fit result " << setw(10) << module->getID() << " chi2 / NDF " << FIXED(10,2) << result.chi2 << " / " << setw(5) << result.ndf << ' ' << setw(5) << fit.numberOfIterations << endl);
384
385 if (fitModel && debug < debug_t) {
386 fit.value.print(cout);
387 }
388
389 DEBUG(fit.value);
390
391 if (writeFits) {
392
393 h0.Fill(result.chi2 / result.ndf);
394 hn.Fill(fit.numberOfIterations);
395 hr.Fill(fit.value.R );
396 h1.Fill(fit.value.p1);
397 h2.Fill(fit.value.p2);
398 h3.Fill(fit.value.p3);
399 h4.Fill(fit.value.p4);
400 hc.Fill(fit.value.cc);
401 hb.Fill(fit.value.bc);
402
403 TH1D h1t(MAKE_CSTRING(module->getID() << ".1t0"), NULL, NUMBER_OF_PMTS, -0.5, NUMBER_OF_PMTS - 0.5);
404 TH1D h1s(MAKE_CSTRING(module->getID() << ".1TTS"), NULL, NUMBER_OF_PMTS, -0.5, NUMBER_OF_PMTS - 0.5);
405 TH1D h1q(MAKE_CSTRING(module->getID() << ".1QE"), NULL, NUMBER_OF_PMTS, -0.5, NUMBER_OF_PMTS - 0.5);
406
407 for (int pmt = 0; pmt != NUMBER_OF_PMTS; ++pmt) {
408 h1t.SetBinContent(pmt + 1, fit.value.parameters[pmt].t0 ());
409 h1t.SetBinError (pmt + 1, fit.error.parameters[pmt].t0 () + numeric_limits<double>::epsilon());
410 h1s.SetBinContent(pmt + 1, fit.value.parameters[pmt].TTS());
411 h1s.SetBinError (pmt + 1, fit.error.parameters[pmt].TTS() + numeric_limits<double>::epsilon());
412 h1q.SetBinContent(pmt + 1, fit.value.parameters[pmt].QE ());
413 h1q.SetBinError (pmt + 1, fit.error.parameters[pmt].QE () + numeric_limits<double>::epsilon());
414 }
415
416 out << h1t << h1s << h1q;
417
418 for (int ix = 1; ix <= h2d->GetXaxis()->GetNbins(); ++ix) {
419
420 const pair_type pair = fit.value.getPair(ix - 1);
421
422 for (int iy = 1; iy <= h2d->GetYaxis()->GetNbins(); ++iy) {
423
424 const double dt_ns = h2d->GetYaxis()->GetBinCenter(iy);
425
426 h2d->SetBinContent(ix, iy, fit.value.getValue(pair, dt_ns));
427 h2d->SetBinError (ix, iy, 0.0);
428 }
429 }
430
431 h2d->SetName(MAKE_CSTRING(module->getID() << _2F));
432 h2d->Write();
433
434 const double x = string.getIndex(module->getString());
435 const double y = module->getFloor();
436
437 H2 .Fill(x, y, result.chi2 / result.ndf);
438 HN->Fill(x, y, fit.numberOfIterations);
439 }
440
441 const double t0 = (fit.value.hasFixedTimeOffset() ? fit.value.getFixedTimeOffset() : 0.0);
442
443 for (int pmt = 0; pmt != NUMBER_OF_PMTS; ++pmt) {
444
445 JPMTParameters& data = parameters[JPMTIdentifier(module->getID(), pmt)];
446
447 const double R = getSurvivalProbability(data) / getSurvivalProbability(JPMTParameters());
448
449 if (R > 0.0)
450 data.QE = fit.value.parameters[pmt].QE / R;
451 else
452 data.QE = 0.0;
453
454 if (writeFits > 1) {
455 data.TTS_ns = fit.value.parameters[pmt].TTS();
456 }
457
458 module->getPMT(pmt).addT0(fit.value.parameters[pmt].t0.get() - t0);
459 }
460 }
461 catch(const exception& error) {
462
463 ERROR("Module " << setw(10) << module->getID() << ' ' << error.what() << " -> set QEs to 0." << endl);
464
465 for (int pmt = 0; pmt != NUMBER_OF_PMTS; ++pmt) {
466 parameters[JPMTIdentifier(module->getID(), pmt)].QE = 0.0;
467 }
468 }
469 }
470
471
472 vector<JMeta> meta(1, JMeta(argc, argv));
473
474 {
475 JSingleFileScanner<JMeta> reader(inputFile);
476 JSTDObjectWriter <JMeta> writer(meta);
477
478 writer << reader;
479 }
480
481 for (vector<JMeta>::const_reverse_iterator i = meta.rbegin(); i != meta.rend(); ++i) {
482 parameters.comment.add(*i);
483 detector .comment.add(*i);
484 }
485
486 if (overwriteDetector) {
487
488 NOTICE("Store calibration data on file " << detectorFile << endl);
489
490 store(detectorFile, detector);
491 }
492
493 if (pmtFile != "") {
494 parameters.store(pmtFile.c_str());
495 }
496
497
498 for (vector<JMeta>::const_iterator i = meta.begin(); i != meta.end(); ++i) {
499 putObject(out,*i);
500 }
501
502 for (JRootFileReader<JDAQHeader> in(inputFile.c_str()); in.hasNext(); ) {
503 putObject(out, *in.next());
504 }
505
506 if (writeFits) {
507 out << h0 << hn << hr << h1 << h2 << h3 << h4 << hc << hb << H2 << *HN;
508 }
509
510 out.Close();
511
512 return 0;
513}
outputFile
string outputFile
Definition JDAQTimesliceSelector.cc:37
DEBUG
#define DEBUG(A)
Message macros.
Definition JMessage.hh:62
NOTICE
#define NOTICE(A)
Definition JMessage.hh:64
FATAL
#define FATAL(A)
Definition JMessage.hh:67
debug
int debug
debug level
Definition JSirene.cc:74
WARNING
#define WARNING(A)
Definition JMessage.hh:65
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition JParser.hh:2142
MAKE_CSTRING
#define MAKE_CSTRING(A)
Make C-string.
Definition JPrint.hh:72
gmake_property
#define gmake_property(A)
macros to convert (template) parameter to JPropertiesElement object
Definition JProperties.hh:1139
distance
std::vector< T >::difference_type distance(typename std::vector< T >::const_iterator first, typename PhysicsEvent::const_iterator< T > second)
Specialisation of STL distance.
Definition build/antares-dataformat/PhysicsEvent.hh:434
JDETECTOR::JDetector
Detector data structure.
Definition JDetector.hh:96
JDETECTOR::JPMTParametersMap
Auxiliary class for map of PMT parameters.
Definition JPMTParametersMap.hh:93
JDETECTOR::JPMTParameters
Data structure for PMT parameters.
Definition JPMTParameters.hh:30
JEEP::JProperties
Utility class to parse parameter values.
Definition JProperties.hh:503
JLANG::JException
General exception.
Definition JException.hh:24
JPARSER::JParser
Utility class to parse command line options.
Definition JParser.hh:1698
JSUPPORT::JSingleFileScanner
Object reading from a list of files.
Definition JSingleFileScanner.hh:172
JTOOLS::JRange
Range of values.
Definition JRange.hh:42
JTOOLS::JRange::getUpperLimit
T getUpperLimit() const
Get upper limit.
Definition JRange.hh:213
JCALIBRATE::TEROSTAT_R1
static double TEROSTAT_R1
scaling factor
Definition JFitK40.hh:877
JCALIBRATE::_2F
static const char *const _2F
Name extension for 2F rate fitted.
Definition JCalibrateK40.hh:37
JCALIBRATE::FIT_PMTS_QE_FIXED_t
@ FIT_PMTS_QE_FIXED_t
fit parameters of PMTs with QE fixed
Definition JFitK40.hh:56
JCALIBRATE::FIT_PMTS_AND_ANGULAR_DEPENDENCE_t
@ FIT_PMTS_AND_ANGULAR_DEPENDENCE_t
fit parameters of PMTs and angular dependence of K40 rate
Definition JFitK40.hh:54
JCALIBRATE::FIT_MODEL_t
@ FIT_MODEL_t
fit parameters of K40 rate and TTSs of PMTs
Definition JFitK40.hh:57
JCALIBRATE::FIT_PMTS_AND_BACKGROUND_t
@ FIT_PMTS_AND_BACKGROUND_t
fit parameters of PMTs and background
Definition JFitK40.hh:55
JCALIBRATE::FIT_PMTS_t
@ FIT_PMTS_t
fit parameters of PMTs
Definition JFitK40.hh:53
JCALIBRATE::_2R
static const char *const _2R
Name extension for 2D rate measured.
Definition JCalibrateK40.hh:36
JCALIBRATE::TEROSTAT_DZ
static double TEROSTAT_DZ
maximal PMT inclination
Definition JFitK40.hh:876
JCALIBRATE::BELL_SHAPE
static double BELL_SHAPE
Bell shape.
Definition JFitK40.hh:878
JDETECTOR::getLabel
std::string getLabel(const JLocation &location)
Get module label for monitoring and other applications.
Definition JLocation.hh:247
JDETECTOR::getRangeOfFloors
floor_range getRangeOfFloors(const JDetector &detector)
Get range of floors.
Definition JDetectorToolkit.hh:295
JDETECTOR::load
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
Definition JDetectorToolkit.hh:520
JDETECTOR::store
void store(const std::string &file_name, const JDetector &detector)
Store detector to output file.
Definition JDetectorToolkit.hh:617
JDETECTOR::getSurvivalProbability
double getSurvivalProbability(const JPMTParameters &parameters)
Get model dependent probability that a one photo-electron hit survives the simulation of the PMT assu...
Definition JPMTParametersToolkit.hh:32
JEEP::debug_t
@ debug_t
debug
Definition JMessage.hh:29
JFIT::model
void model(JModel_t &value)
Auxiliary function to constrain model during fit.
Definition JGandalf.hh:57
JPP
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Definition JAAnetToolkit.hh:48
JROOT::putObject
bool putObject(TDirectory &dir, const TObject &object)
Write object to ROOT directory.
Definition JRootFileWriter.hh:41
JTOOLS::result
return result
Definition JPolint.hh:862
KM3NETDAQ
KM3NeT DAQ data structures and auxiliaries.
Definition DataQueue.cc:39
KM3NETDAQ::NUMBER_OF_PMTS
static const int NUMBER_OF_PMTS
Total number of PMTs in module.
Definition JDAQ.hh:26
FILL
Auxiliary data structure for sequence of same character.
Definition JManip.hh:330
FIXED
Auxiliary data structure for floating point format specification.
Definition JManip.hh:448
JAANET::JRange_t
Type definition of range.
Definition JHead.hh:43
JAANET::K40
Livetime of noise data.
Definition JHead.hh:1062
JAANET::detector
Detector file.
Definition JHead.hh:227
JACOUSTICS::JFit
Acoustic single fit.
Definition JAcoustics/JEvt.hh:36
JACOUSTICS::JModel
Model for fit to acoustics data.
Definition JAcoustics/JModel.hh:480
JCALIBRATE::JK40Parameters
Fit parameters for two-fold coincidence rate due to K40.
Definition JFitK40.hh:729
JCALIBRATE::JK40Parameters::getInstance
static const JK40Parameters & getInstance()
Get default values.
Definition JFitK40.hh:745
JCALIBRATE::JPMTParameters_t::getInstance
static const JPMTParameters_t & getInstance()
Get default values.
Definition JFitK40.hh:475
JCALIBRATE::JTDC_t
Auxiliary class for TDC constraints.
Definition JTDC_t.hh:39
JCALIBRATE::JTDC_t::equal_range
range_type equal_range(const int id)
Get range of constraints for given module.
Definition JTDC_t.hh:101
JCALIBRATE::JTDC_t::reverse
void reverse()
Reverse constraints.
Definition JTDC_t.hh:137
JCALIBRATE::JTDC_t::is_valid
bool is_valid(const bool option=false) const
Check validity of TDC constrants.
Definition JTDC_t.hh:171
JCALIBRATE::data_type
Data structure for measured coincidence rates of all pairs of PMTs in optical module.
Definition JFitK40.hh:103
JCALIBRATE::rate_type
Data structure for measured coincidence rate of pair of PMTs.
Definition JFitK40.hh:66
JDETECTOR::JStringRouter
Router for mapping of string identifier to index.
Definition JStringRouter.hh:23
JEEP::JComment::add
JComment & add(const std::string &comment)
Add comment.
Definition JComment.hh:100
JLANG::JObjectStreamIO::store
void store(const char *file_name) const
Store to output file.
Definition JObjectStreamIO.hh:41
JLANG::JObjectStreamIO::load
void load(const char *file_name)
Load from input file.
Definition JObjectStreamIO.hh:30
JPARSER::initialised
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Definition JParser.hh:68
JSUPPORT::JMeta
Auxiliary class for ROOT I/O of application specific meta data.
Definition JMeta.hh:72
JTOOLS::JCombinatorics::pair_type
Data structure for a pair of indices.
Definition JCombinatorics.hh:29
Generated by doxygen 1.12.0