25 int main(
int argc,
char **argv)
31 JPMTParameters parameters;
36 JProperties properties = parameters.getProperties();
38 JParser<> zap(
"Example program to histogram time over threshold as a function of number of photo-electrons.");
46 catch(
const exception &error) {
47 FATAL(error.what() << endl);
52 cout <<
"PMT parameters:" << endl;
56 JPMTAnalogueSignalProcessor cpu(parameters);
60 TH1D h0(
"[model]", NULL, 510, 0.0, 255.0);
61 TH1D h1(
"[parametrisation]", NULL, 510, 0.0, 255.0);
63 for (
int i = 1; i <= h0.GetNbinsX(); ++i) {
67 const double x = h0.GetBinCenter(i);
68 const double npe = cpu.getNPE(x);
69 const double y = cpu.getRiseTime(npe);
71 h0.SetBinContent(i, y);
74 catch(
const JValueOutOfRange& error) {}