Wrapper class to make final fit of muon trajectory. More...

#include <JMuonGandalf.hh>

Inheritance diagram for JRECONSTRUCTION::JMuonGandalf:

Public Types
typedef JRegressor< JLine3Z, JGandalf >	JRegressor_t

typedef JTRIGGER::JHitL0	hit_type

typedef std::vector< hit_type >	buffer_type

Public Member Functions
	JMuonGandalf (const JMuonGandalfParameters_t &parameters, const JModuleRouter &router, const JSummaryRouter &summary, const std::string &pdf_file, const int debug=0)
	Constructor. More...

JEvt	operator() (const KM3NETDAQ::JDAQEvent &event, const JEvt &in)
	Fit function. More...

JEvt	operator() (const buffer_type &data, const JEvt &in)
	Fit function. More...

void	reset ()
	Reset fit parameters. More...

bool	equals (const JMuonGandalfParameters_t &parameters) const
	Equality. More...

	ClassDef (JMuonGandalfParameters_t, 1)

Public Attributes
const JModuleRouter &	router

const JSummaryRouter &	summary

int	debug

double	roadWidth_m

double	R_Hz

size_t	numberOfPrefits

double	TTS_ns

double	E_GeV

double	TMin_ns

double	TMax_ns

double	ZMin_m

double	ZMax_m

bool	reprocess

Static Public Attributes
struct {
}	compare
	Auxiliary data structure for sorting of hits. More...

Detailed Description

Wrapper class to make final fit of muon trajectory.

The JMuonGandalf fit uses one or more start values (usually taken from the output of JMuonSimplex).
All hits of which the PMT position lies within a set road width (JMuonGandalfParameters_t::roadWidth_m) and time is within a set window (JMuonGandalfParameters_t::TMin_ns, JMuonGandalfParameters_t::TMax_ns) around the Cherenkov hypothesis are taken.
In case there are multiple hits from the same PMT is the specified window, the first hit is taken and the other hits are discarded.
The PDF is accordingly evaluated, i.e. the normalised probability for a first hit at the given time of the hit is taken. The normalisation is consistently based on the specified time window.
Note that this hit selection is unbiased with respect to the PDF of a single PMT.

Definition at line 72 of file JMuonGandalf.hh.

Member Typedef Documentation

typedef JRegressor<JLine3Z, JGandalf> JRECONSTRUCTION::JMuonGandalf::JRegressor_t

Definition at line 76 of file JMuonGandalf.hh.

typedef JTRIGGER::JHitL0 JRECONSTRUCTION::JMuonGandalf::hit_type

Definition at line 77 of file JMuonGandalf.hh.

typedef std::vector<hit_type> JRECONSTRUCTION::JMuonGandalf::buffer_type

Definition at line 78 of file JMuonGandalf.hh.

Constructor & Destructor Documentation

JRECONSTRUCTION::JMuonGandalf::JMuonGandalf	(	const JMuonGandalfParameters_t &	parameters,
		const JModuleRouter &	router,
		const JSummaryRouter &	summary,
		const std::string &	pdf_file,
		const int	debug = `0`
	)

inline

Constructor.

Parameters

parameters	parameters
router	module router
summary	summary file router
pdf_file	PDF file
debug	debug

Definition at line 91 of file JMuonGandalf.hh.

                                                             :
       JMuonGandalfParameters_t(parameters),
       JRegressor_t(pdf_file, parameters.TTS_ns),
       router (router),
       summary(summary)
     {
       using namespace JFIT;
 
       if (this->getRmax() < roadWidth_m) {
         roadWidth_m = this->getRmax();
       }
 
       JRegressor_t::debug              = debug;
       JRegressor_t::T_ns.setRange(TMin_ns, TMax_ns);
       JRegressor_t::Vmax_npe           =   10.0;
       JRegressor_t::MAXIMUM_ITERATIONS = 1000;
       
       this->parameters.resize(5);
       
       this->parameters[0] = JLine3Z::pX();
       this->parameters[1] = JLine3Z::pY();
       this->parameters[2] = JLine3Z::pT();
       this->parameters[3] = JLine3Z::pDX();
       this->parameters[4] = JLine3Z::pDY();
     }

Member Function Documentation

JEvt JRECONSTRUCTION::JMuonGandalf::operator()	(	const KM3NETDAQ::JDAQEvent &	event,
		const JEvt &	in
	)

inline

Fit function.

Parameters

event	event
in	start values

Returns: fit results

Definition at line 129 of file JMuonGandalf.hh.

     {
       using namespace std;
       using namespace JFIT;
       using namespace JTRIGGER;
 
       const JBuildL0<hit_type>  buildL0;
 
       buffer_type dataL0;
 
       buildL0(event, router, true, back_inserter(dataL0));
 
       return (*this)(dataL0, in);
     }

JEvt JRECONSTRUCTION::JMuonGandalf::operator()	(	const buffer_type &	data,
		const JEvt &	in
	)

inline

Fit function.

Parameters

data	hit data
in	start values

Returns: fit results

Definition at line 152 of file JMuonGandalf.hh.

     {
       using namespace std;
       using namespace JFIT;
       using namespace JGEOMETRY3D;
 
       JEvt out;
       
       for (JEvt::const_iterator track = in.begin(); track != in.end(); ++track) {
 
         const JRotation3D  R (getDirection(*track));
         const JLine1Z      tz(getPosition (*track).rotate(R), track->getT());
         JRange<double>     Z_m;
       
         if (track->hasW(JSTART_LENGTH_METRES)       && 
             track->getW(JSTART_LENGTH_METRES) > 0.0) {
           Z_m = JZRange(ZMin_m, ZMax_m + track->getW(JSTART_LENGTH_METRES));
         }
 
         const JModel<JLine1Z> match(tz, roadWidth_m, JRegressor_t::T_ns, Z_m);
 
         // hit selection based on start value
 
         vector<JHitW0> buffer;
           
         for (buffer_type::const_iterator i = data.begin(); i != data.end(); ++i) {
             
           const double rate_Hz = summary.getRate(*i);
 
           JHitW0 hit(*i, rate_Hz);
             
           hit.rotate(R);
 
           if (match(hit)) {
             buffer.push_back(hit);
           }
         }
           
         // select first hit
           
         sort(buffer.begin(), buffer.end(), compare);
           
         vector<JHitW0>::iterator __end = unique(buffer.begin(), buffer.end(), equal_to<JDAQPMTIdentifier>());
           
 
         const int NDF = distance(buffer.begin(), __end) - this->parameters.size();
 
         if (NDF > 0) {
 
           // set fit parameters
 
           if (track->getE() > 0.1)
             JRegressor_t::E_GeV   = track->getE();
           else
             JRegressor_t::E_GeV   = this->JMuonGandalfParameters_t::E_GeV;
 
           const double chi2 = (*this)(JLine3Z(tz), buffer.begin(), __end);
             
           JTrack3D tb(this->value);
             
           tb.rotate_back(R);
 
           out.push_back(getFit(JHistory(track->getHistory()).add(JMUONGANDALF), tb, getQuality(chi2), NDF));
 
           // set additional values
 
           out.rbegin()->setW(track->getW());
           out.rbegin()->setW(JGANDALF_BETA0_RAD,            sqrt(this->error.getDX() * this->error.getDX()  +
                                                                  this->error.getDY() * this->error.getDY()));
           out.rbegin()->setW(JGANDALF_BETA1_RAD,            sqrt(this->error.getDX() * this->error.getDY()));
           out.rbegin()->setW(JGANDALF_CHI2,                 chi2);
           out.rbegin()->setW(JGANDALF_NUMBER_OF_HITS,       distance(buffer.begin(), __end));
           out.rbegin()->setW(JGANDALF_LAMBDA,               this->lambda);
           out.rbegin()->setW(JGANDALF_NUMBER_OF_ITERATIONS, this->numberOfIterations);
         }
       }
 
       return out;
     }

void JRECONSTRUCTION::JMuonGandalfParameters_t::reset ( )

inlineinherited

Reset fit parameters.

Definition at line 41 of file JMuonGandalfParameters_t.hh.

     {
       roadWidth_m     = std::numeric_limits<double>::max();
       R_Hz            = 6.0e3;
       numberOfPrefits = 0;
       TTS_ns          = 2;
       E_GeV           = 1.0e3;
       TMin_ns         =  -50.0;     
       TMax_ns         = +450.0;     
       ZMin_m          =    0.0;     
       ZMax_m          =    0.0;     
       reprocess       = false;
     }

bool JRECONSTRUCTION::JMuonGandalfParameters_t::equals ( const JMuonGandalfParameters_t & parameters ) const

inlineinherited

Equality.

Parameters

parameters fit parameters

Returns: true if equals; else false

Definition at line 61 of file JMuonGandalfParameters_t.hh.

     {
       return (this->roadWidth_m      == parameters.roadWidth_m     &&
               this->R_Hz             == parameters.R_Hz            &&
               this->numberOfPrefits  == parameters.numberOfPrefits &&
               this->TTS_ns           == parameters.TTS_ns          &&
               this->E_GeV            == parameters.E_GeV           &&
               this->TMin_ns          == parameters.TMin_ns         &&
               this->TMax_ns          == parameters.TMax_ns         &&
               this->ZMin_m           == parameters.ZMin_m          &&
               this->ZMax_m           == parameters.ZMax_m          &&
               this->reprocess        == parameters.reprocess);
     }