ARCA muon reconstruction. More...

#include <JARCAMuonReconstruction.hh>

Classes
struct	input_type
	Input data type. More...

struct	JParameters_t
	Parameters. More...

struct	JStorage_t
	Storage for PDFs. More...

Public Member Functions
	JARCAMuonReconstruction (const JParameters_t &parameters, const JStorage_t &storage, const JK40Rates &rates_Hz, const JPMTParametersMap &pmtParameters, const JEnergyCorrection &correct, const int debug)
	Constructor.

input_type	getInput (const JModuleRouter &router, const JSummaryRouter &summary, const JDAQEvent &event, const coverage_type &coverage) const
	Get input data.

JEvt	operator() (input_type &input)
	Fit function.

Public Attributes
JMuonPrefit	prefit

JMuonSimplex	simplex

JMuonGandalf	gandalf

JMuonStart	start

JMuonEnergy	energy

JMuonFeatures	features

std::vector< JLANG::JClock >	clock

Detailed Description

ARCA muon reconstruction.

Definition at line 41 of file JARCAMuonReconstruction.hh.

Constructor & Destructor Documentation

◆ JARCAMuonReconstruction()

JRECONSTRUCTION::JARCAMuonReconstruction::JARCAMuonReconstruction	(	const JParameters_t &	parameters,
		const JStorage_t &	storage,
		const JK40Rates &	rates_Hz,
		const JPMTParametersMap &	pmtParameters,
		const JEnergyCorrection &	correct,
		const int	debug )

inline

Constructor.

Parameters

parameters	parameters
storage	storage
rates_Hz	K40 rates [Hz]
pmtParameters	PMT parameters
correct	energy correction
debug	debug

Definition at line 137 of file JARCAMuonReconstruction.hh.

                                                            :
      prefit  (parameters.prefit,   debug),
      simplex (parameters.simplex,  debug),
      gandalf (parameters.gandalf,  storage.gandalf,  pmtParameters, debug),
      start   (parameters.start,    storage.start,    rates_Hz, debug),
      energy  (parameters.energy,   storage.energy,   pmtParameters, correct,  debug),
      features(parameters.features, storage.features, debug)
    {
      clock.resize(6);
 
      clock[0].setTitle("prefit");
      clock[1].setTitle("simplex");
      clock[2].setTitle("gandalf");
      clock[3].setTitle("start");
      clock[4].setTitle("energy");
      clock[5].setTitle("features");
    }

Member Function Documentation

◆ getInput()

input_type JRECONSTRUCTION::JARCAMuonReconstruction::getInput	(	const JModuleRouter &	router,
		const JSummaryRouter &	summary,
		const JDAQEvent &	event,
		const coverage_type &	coverage ) const

inline

Get input data.

Parameters

router	module router
summary	summary data
event	event
coverage	coverage

Returns: input data

Definition at line 170 of file JARCAMuonReconstruction.hh.

    {
      input_type input(event.getDAQEventHeader());
 
      input.prefit   = prefit  .getInput(router, event, coverage);
      input.simplex  = simplex .getInput(router, event, JEvt(), coverage);
      input.gandalf  = gandalf .getInput(router, summary, event, JEvt(), coverage);
      input.start    = start   .getInput(router, summary, event, JEvt(), coverage);
      input.energy   = energy  .getInput(router, summary, event, JEvt(), coverage);
      input.features = features.getInput(router, summary, event, JEvt(), coverage);
 
      return input;
    }

◆ operator()()

JEvt JRECONSTRUCTION::JARCAMuonReconstruction::operator() ( input_type & input )

inline

Fit function.

Parameters

input input data

Returns: fit results

Definition at line 194 of file JARCAMuonReconstruction.hh.

    {
      std::vector<int> apps;
 
      {
        JLANG::JClock::sentry sentry(clock[0]);
 
        input.simplex.in  = prefit(input.prefit);
 
        apps.push_back(JMUONPREFIT);
      }
      {
        JLANG::JClock::sentry sentry(clock[1]);
 
        input.gandalf.in  = simplex(input.simplex);
 
        apps.push_back(JMUONSIMPLEX);
      }
      {
        JLANG::JClock::sentry sentry(clock[2]);
 
        input.start.in    = gandalf(input.gandalf);
 
        apps.push_back(JMUONGANDALF);
      }
      {
        JLANG::JClock::sentry sentry(clock[3]);
 
        input.energy.in   = start(input.start);
 
        apps.push_back(JMUONSTART);
      }
      {
        JLANG::JClock::sentry sentry(clock[4]);
 
        input.features.in = energy(input.energy);
 
        apps.push_back(JMUONENERGY);
      }
 
      JLANG::JClock::sentry sentry(clock[5]);
 
      JEvt out = features(input.features);
 
      apps.push_back(JMUONFEATURES);
 
      for (JFIT::JEvt::iterator i = out.begin(); i != out.end(); ++i) {
        i->setStatus(i->getHistory().getStatus(apps) ? COMPLETE_CHAIN : INCOMPLETE_CHAIN);
      }
 
      return out;
    }