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

#include <JMuonSimplex.hh>

Inheritance diagram for JRECONSTRUCTION::JMuonSimplex:

Public Types
typedef JRegressor< JLine3Z, JSimplex >	JRegressor_t

typedef JTRIGGER::JHitR1	hit_type

typedef std::vector< hit_type >	buffer_type

Public Member Functions
	JMuonSimplex (const JMuonSimplexParameters_t &parameters, const JModuleRouter &router, const int debug=0)
	Constructor. More...

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

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

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

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

	ClassDef (JMuonSimplexParameters_t, 2)

Public Attributes
const JModuleRouter &	router

double	sigma_ns
	time resolution [ns] More...

bool	useL0
	option for L0 hit use More...

size_t	numberOfPrefits
	number of prefits More...

double	TMaxLocal_ns
	time window for local coincidences [ns] More...

double	ctMin
	minimal cosine space angle between PMT axes More...

double	roadWidth_m
	road width [m] More...

double	TMin_ns
	minimal time w.r.t. Cherenkov hypothesis [ns] More...

double	TMax_ns
	maximal time w.r.t. Cherenkov hypothesis [ns] More...

int	NMax
	maximum number of iterations More...

bool	reprocess
	reprocess More...

Detailed Description

Wrapper class to make intermediate fit of muon trajectory.

The JMuonSimplex fit uses one or more start values (usually taken from the output of JMuonPrefit) and produces new start values for subsequent fits (usually JMuonGandalf).
All hits of which the PMT position lies within a set road width (JMuonSimplexParameters_t::roadWidth_m) and time is within a set window (JMuonSimplexParameters_t::TMin_ns, JMuonSimplexParameters_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 chi-squared is based on an M-estimator of the time residuals.

Definition at line 65 of file JMuonSimplex.hh.

Member Typedef Documentation

typedef JRegressor<JLine3Z, JSimplex> JRECONSTRUCTION::JMuonSimplex::JRegressor_t

Definition at line 69 of file JMuonSimplex.hh.

typedef JTRIGGER::JHitR1 JRECONSTRUCTION::JMuonSimplex::hit_type

Definition at line 70 of file JMuonSimplex.hh.

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

Definition at line 71 of file JMuonSimplex.hh.

Constructor & Destructor Documentation

JRECONSTRUCTION::JMuonSimplex::JMuonSimplex	(	const JMuonSimplexParameters_t &	parameters,
		const JModuleRouter &	router,
		const int	debug = `0`
	)

inline

Constructor.

Parameters

parameters	parameters
router	module router
debug	debug

Definition at line 82 of file JMuonSimplex.hh.

                                                             :
       JMuonSimplexParameters_t(parameters),
       JRegressor_t(parameters.sigma_ns),
       router(router)
     {
       using namespace JFIT;
 
       this->estimator.reset(new JMEstimatorLorentzian());
 
       JRegressor_t::debug              = debug;
       JRegressor_t::MAXIMUM_ITERATIONS = NMax;
     }

Member Function Documentation

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

inline

Fit function.

Parameters

event	event
in	start values

Returns: fit results

Definition at line 105 of file JMuonSimplex.hh.

     {
       using namespace std;
       using namespace JFIT;
       using namespace JTRIGGER;
 
       const JBuildL0<hit_type> buildL0;
       const JBuildL2<hit_type> buildL2(JL2Parameters(2, TMaxLocal_ns, ctMin));
 
       buffer_type dataL0;
       buffer_type dataL1;
 
       const KM3NETDAQ::JDAQTimeslice timeslice(event, true);
       JSuperFrame2D<JHit>            buffer;
 
       for (JDAQTimeslice::const_iterator i = timeslice.begin(); i != timeslice.end(); ++i) {
 
         if (router.hasModule(i->getModuleID())) {
 
           buffer(*i, router.getModule(i->getModuleID()));
 
           if (useL0) {
             buildL0(buffer, back_inserter(dataL0));
           }
           {
             buildL2(buffer, back_inserter(dataL1));
           }
         }
       }
 
       return (*this)(dataL0, dataL1, in);
     }

JEvt JRECONSTRUCTION::JMuonSimplex::operator()	(	const buffer_type &	dataL0,
		const buffer_type &	dataL1,
		const JEvt &	in
	)

inline

Fit function.

Parameters

dataL0	L0 hit data
dataL1	L1 hit data
in	start values

Returns: fit results

Definition at line 147 of file JMuonSimplex.hh.

     {
       using namespace std;
       using namespace JFIT;
       using namespace JGEOMETRY3D;
       using namespace JLANG;
 
       JEvt out;
 
       buffer_type data;
 
       data.reserve(dataL0.size() + 
                    dataL1.size());
       
 
       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());
         const JModel<JLine1Z> match(tz, roadWidth_m, JTimeRange(TMin_ns, TMax_ns));
             
         data.clear();
 
         for (buffer_type::const_iterator i = dataL1.begin(); i != dataL1.end(); ++i) {
 
           hit_type hit(*i);
               
           hit.rotate(R);
 
           if (match(hit)) {
             data.push_back(hit);
           }
         }
 
         if (useL0) {
 
           buffer_type::iterator __end = data.end();     
 
           for (buffer_type::const_iterator i = dataL0.begin(); i != dataL0.end(); ++i) {
               
             if (find_if(data.begin(), __end, make_predicate(&hit_type::getModuleID, i->getModuleID())) == __end) {
 
               hit_type hit(*i);
               
               hit.rotate(R);
               
               if (match(hit)) {
                 data.push_back(hit);
               }
             }
           }
         }
 
 
         this->step.resize(5);
         
         this->step[0] = JLine3Z(JLine1Z(JVector3D(0.5, 0.0, 0.0), 0.0));
         this->step[1] = JLine3Z(JLine1Z(JVector3D(0.0, 0.5, 0.0), 0.0));
         this->step[2] = JLine3Z(JLine1Z(JVector3D(0.0, 0.0, 0.0), 1.0));
         this->step[3] = JLine3Z(JLine1Z(JVector3D(), 0.0), JVersor3Z(0.005, 0.0));
         this->step[4] = JLine3Z(JLine1Z(JVector3D(), 0.0), JVersor3Z(0.0, 0.005));
         
         const int NDF = getCount(data.begin(), data.end()) - this->step.size();
 
         if (NDF > 0) {
 
           const double chi2 = (*this)(JLine3Z(tz), data.begin(), data.end());
               
           JTrack3D tb(this->value);
               
           tb.rotate_back(R);
 
           out.push_back(getFit(JHistory(track->getHistory()).add(JMUONSIMPLEX), tb, getQuality(chi2, NDF), NDF)); 
 
           out.rbegin()->setW(track->getW());
         }
       }
 
       return out;
     }

void JRECONSTRUCTION::JMuonSimplexParameters_t::reset ( )

inlineinherited

Reset fit parameters.

Definition at line 39 of file JMuonSimplexParameters_t.hh.

     {
       sigma_ns        = 3;
       useL0           = true;
       numberOfPrefits = 0;
       TMaxLocal_ns    = 15.0;
       ctMin           = 0.0;
       roadWidth_m     = 200.0;
       TMin_ns         = -50.0;
       TMax_ns         = +50.0;
       NMax            = 1000;
       reprocess       = false; 
     }

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

inlineinherited

Equality.

Parameters

parameters fit parameters

Returns: true if equals; else false

Definition at line 59 of file JMuonSimplexParameters_t.hh.

     {
       return (this->sigma_ns         == parameters.sigma_ns         &&
               this->useL0            == parameters.useL0            &&
               this->numberOfPrefits  == parameters.numberOfPrefits  &&
               this->TMaxLocal_ns     == parameters.TMaxLocal_ns     &&
               this->ctMin            == parameters.ctMin            &&
               this->roadWidth_m      == parameters.roadWidth_m      &&
               this->TMin_ns          == parameters.TMin_ns          &&
               this->TMax_ns          == parameters.TMax_ns          &&
               this->NMax             == parameters.NMax             &&
               this->reprocess        == parameters.reprocess);
     }