Jpp - the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Functions
JSlewing.cc File Reference

Example program to histogram time over threshold as a function of number of photo-electrons. More...

#include <string>
#include <iostream>
#include <iomanip>
#include <cmath>
#include "TROOT.h"
#include "TFile.h"
#include "TH1D.h"
#include "JDetector/JPMTAnalogueSignalProcessor.hh"
#include "JTrigger/JGetRiseTime.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 as a function of number of photo-electrons.

Author
mdejong

Definition in file JSlewing.cc.

Function Documentation

int main ( int  argc,
char **  argv 
)

Definition at line 25 of file JSlewing.cc.

26 {
27  using namespace std;
28  using namespace JPP;
29 
30  string outputFile;
31  JPMTParameters parameters;
32  int debug;
33 
34  try {
35 
36  JProperties properties = parameters.getProperties();
37 
38  JParser<> zap("Example program to histogram time over threshold as a function of number of photo-electrons.");
39 
40  zap['o'] = make_field(outputFile) = "histogram.root";
41  zap['P'] = make_field(properties) = JPARSER::initialised();
42  zap['d'] = make_field(debug) = 0;
43 
44  zap(argc, argv);
45  }
46  catch(const exception &error) {
47  FATAL(error.what() << endl);
48  }
49 
50 
51  if (debug >= JEEP::debug_t) {
52  cout << "PMT parameters:" << endl;
53  cout << parameters.getProperties(JEquationParameters("=", "\n", "", "")) << endl;
54  }
55 
56  JPMTAnalogueSignalProcessor cpu(parameters);
57 
58  TFile out(outputFile.c_str(), "recreate");
59 
60  TH1D h0("[model]", NULL, 510, 0.0, 255.0);
61  TH1D h1("[parametrisation]", NULL, 510, 0.0, 255.0);
62 
63  for (int i = 1; i <= h0.GetNbinsX(); ++i) {
64 
65  try {
66 
67  const double x = h0.GetBinCenter(i);
68  const double npe = cpu.getNPE(x);
69  const double y = cpu.getRiseTime(npe);
70 
71  h0.SetBinContent(i, y);
72  h1.SetBinContent(i, getRiseTime(x));
73  }
74  catch(const JValueOutOfRange& error) {}
75  }
76 
77  out.Write();
78  out.Close();
79 }
Utility class to parse command line options.
Definition: JParser.hh:1500
debug
Definition: JMessage.hh:29
Utility class to parse parameter values.
Definition: JProperties.hh:496
*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
Simple data structure to support I/O of equations (see class JLANG::JEquation).
then for HISTOGRAM in h0 h1
Definition: JMatrixNZ.sh:69
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Definition: JParser.hh:66
string outputFile
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition: JParser.hh:1961
static const JGetRiseTime getRiseTime
Function object for rise time evaluation.
int debug
debug level
Definition: JSirene.cc:63
#define FATAL(A)
Definition: JMessage.hh:67