JToT.cc File Reference

Example program to histogram time over threshold distributions. More...

#include <string>
#include <iostream>
#include <iomanip>
#include "TROOT.h"
#include "TFile.h"
#include "TH1D.h"
#include "JDetector/JPMTParameters.hh"
#include "JDetector/JPMTAnalogueSignalProcessor.hh"
#include "JDetector/JPMTDefaultSimulator.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

Example program to histogram time over threshold distributions.

Author

mdejong

Definition in file JToT.cc.

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 25 of file JToT.cc.

{
  using namespace std;
  using namespace JPP;

  string         outputFile;
  JPMTParameters parameters;
  int            NPE;
  int            numberOfHits;
  double         precision;
  int            debug;

  try {

    JProperties properties = parameters.getProperties();

    JParser<> zap("Example program to histogram time over threshold distributions.");

    zap['o'] = make_field(outputFile)       = "tot.root";
    zap['P'] = make_field(properties)       = JPARSER::initialised();
    zap['N'] = make_field(NPE)              = 1;
    zap['n'] = make_field(numberOfHits)     = 1000000;
    zap['e'] = make_field(precision)        = 0.005;
    zap['d'] = make_field(debug)            = 0;

    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }


  if (debug >= JEEP::debug_t) {
    cout << "PMT parameters:" << endl;
    cout << parameters.getProperties(JEquationParameters("=", "\n", "", "")) << endl;
  }

  const JPMTAnalogueSignalProcessor cpu(parameters);

  const JPMTIdentifier pmt(1,0);
  const JCalibration   calibration;

  const JPMTDefaultSimulator simulator(parameters, pmt);

  JPMTData<JPMTSignal> input;
  JPMTData<JPMTPulse>  output;

  ASSERT(numberOfHits > 0);


  TFile out(outputFile.c_str(), "recreate");

  const double xmin =  -0.5;
  const double xmax = 255.5;
  const double dx   =   1.0;
  const int    nx   = (int) ((xmax - xmin) / dx);

  TH1D h0("h0", NULL, nx, xmin, xmax);
  TH1D h1("h1", NULL, nx, xmin, xmax);
  TH1D h2("h2", NULL, nx, xmin, xmax);

  h1.Sumw2();
  h2.Sumw2();

  DEBUG("    ToT      npe   dP/dnpe dP/dToT" << endl);

  for (int i = 1; i <= h0.GetNbinsX(); ++i) {

    const double x   = h0.GetBinCenter(i);
    const double npe = cpu.getNPE(x);
    const double y   = cpu.getProbability(npe, NPE);
    const double v   = cpu.getDerivative(npe);

    DEBUG("  "
          << FIXED(5,1) << x   << "   "
          << FIXED(7,3) << npe << "   "
          << FIXED(5,3) << y   << "   "
          << FIXED(6,4) << y*v << endl);

    h0.SetBinContent(i, y*v);
  }


  if (numberOfHits > 0) {

    for (int i = 0; i != numberOfHits; ++i) {
      
      const double npe = cpu.getRandomAmplitude(NPE);

      if (cpu.applyThreshold(npe)) {

        const double tot_ns = cpu.getTimeOverThreshold(npe);

        h1.Fill(tot_ns);
      }
    }

    h1.Scale(1.0 / (double) numberOfHits / dx);
  }


  if (numberOfHits > 0) {

    const double t_ns = 25.0;
    
    for (int i = 0; i != numberOfHits; ++i) {
      
      input .clear();
      output.clear();

      input.push_back(JPMTSignal(t_ns, NPE));

      simulator.processHits(pmt, calibration, input, output);

      for (JPMTData<JPMTPulse>::const_iterator hit = output.begin(); hit != output.end(); ++hit) {
        h2.Fill(hit->tot_ns);
      }
    }

    h2.Scale(1.0 / (double) numberOfHits / dx);
  }

  for (int i = 1; i <= h0.GetNbinsX(); ++i) {

    const Double_t x  = h0.GetBinCenter (i);
    const Double_t y0 = h0.GetBinContent(i);
    const Double_t y1 = h1.GetBinContent(i);
    const Double_t y2 = h2.GetBinContent(i);

    DEBUG("[" << setw(3) << i << "]"  << ' ' 
          << FIXED(5,1) << x  << ' '
          << FIXED(6,4) << y0 << ' '
          << FIXED(6,4) << y1 << ' '
          << FIXED(6,4) << y2 << endl);

    ASSERT(fabs(y0 - y1) < precision);
    ASSERT(fabs(y0 - y2) < precision);
  }


  out.Write();
  out.Close();
}