JPulseFinder.cc File Reference

#include <iostream>
#include <iomanip>
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JModuleLocation.hh"
#include "TH1D.h"
#include "TH2D.h"
#include "JPulseFinder_IO.hh"
#include "Physics.hh"
#include "Detector.hh"
#include "PMT_choice.hh"
#include "JPulseFinder_utils.hh"

Go to the source code of this file.

Functions
int	main (int argc, char **argv)

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 34 of file JPulseFinder.cc.

                                {

  //------------------------------------------------------------
  //  READ COMMAND-LINE OPTIONS AND DECLARE SOME VARIABLES
  //------------------------------------------------------------

  IO options ;

  read_user_options(options, argc , argv);

  vector < int > top_pmts = setTopPMTs(options.up_pmts) ;

  vector < int > bottom_pmts = setBottomPMTs(options.down_pmts) ;

  int max_d = options.number_neighbors ;

  const double cwater = c_water(470 , 3000);

  //------------------------------------------------------------
  //  LOAD DETECTOR FROM DETECTOR FILE
  //------------------------------------------------------------

  JDetector detector = loadDetector(options.detector_file) ;

  //------------------------------------------------------------
  //  ALLOCATE HISTOGRAMS
  //------------------------------------------------------------

  cout << "allocating histograms" << endl;

  vector < vector < TH2D* > > histograms = allocate_histograms(detector , options , cwater, top_pmts , bottom_pmts , max_d) ;

  //------------------------------------------------------------
  //  READ TIMESLICES AND FILL HISTOGRAMS
  //------------------------------------------------------------

  cout << "reading time slices" << endl;

  fill_histograms (histograms , options , detector, top_pmts , bottom_pmts , max_d);

  //------------------------------------------------------------
  //  PROJECT 2D HISTOGRAMS TO OBTAIN TIME AND TOT DISTRIBUTIONS
  //------------------------------------------------------------

  vector < vector < TH1D* > > hit_time = project_time(histograms) ;

  vector < vector < TH1D* > > hit_ToT = project_ToT(histograms) ; 
  

  //------------------------------------------------------------
  //  WRITE OUTPUT
  //------------------------------------------------------------

  cout << "writing output " << endl;

  write_output(options , histograms , hit_ToT , hit_time) ;

  //------------------------------------------------------------
  // DONE
  //------------------------------------------------------------

  cout << endl
       << "---------- Done! ----------" 
       << endl ;
  cout << "Output in '" << options.ofname << "'." << endl ;
}