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

◆ JRegressor_t

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

Definition at line 69 of file JMuonSimplex.hh.

◆ hit_type

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

Definition at line 70 of file JMuonSimplex.hh.

◆ buffer_type

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

Definition at line 71 of file JMuonSimplex.hh.

Constructor & Destructor Documentation

◆ JMuonSimplex()

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

◆ operator()() [1/2]

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);
     }

◆ operator()() [2/2]

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;
     }

◆ reset()

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; 
     }

◆ equals()

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);
     }

◆ ClassDef()

JRECONSTRUCTION::JMuonSimplexParameters_t::ClassDef	(	JMuonSimplexParameters_t	,
		2
	)

inherited

Member Data Documentation

◆ router

const JModuleRouter& JRECONSTRUCTION::JMuonSimplex::router

Definition at line 231 of file JMuonSimplex.hh.

◆ sigma_ns

double JRECONSTRUCTION::JMuonSimplexParameters_t::sigma_ns

inherited

time resolution [ns]

Definition at line 75 of file JMuonSimplexParameters_t.hh.

◆ useL0

bool JRECONSTRUCTION::JMuonSimplexParameters_t::useL0

inherited

option for L0 hit use

Definition at line 76 of file JMuonSimplexParameters_t.hh.

◆ numberOfPrefits

size_t JRECONSTRUCTION::JMuonSimplexParameters_t::numberOfPrefits

inherited

number of prefits

Definition at line 77 of file JMuonSimplexParameters_t.hh.

◆ TMaxLocal_ns

double JRECONSTRUCTION::JMuonSimplexParameters_t::TMaxLocal_ns

inherited

time window for local coincidences [ns]

Definition at line 78 of file JMuonSimplexParameters_t.hh.

◆ ctMin

double JRECONSTRUCTION::JMuonSimplexParameters_t::ctMin

inherited

minimal cosine space angle between PMT axes

Definition at line 79 of file JMuonSimplexParameters_t.hh.

◆ roadWidth_m

double JRECONSTRUCTION::JMuonSimplexParameters_t::roadWidth_m

inherited

road width [m]

Definition at line 80 of file JMuonSimplexParameters_t.hh.

◆ TMin_ns

double JRECONSTRUCTION::JMuonSimplexParameters_t::TMin_ns

inherited

minimal time w.r.t. Cherenkov hypothesis [ns]

Definition at line 81 of file JMuonSimplexParameters_t.hh.

◆ TMax_ns

double JRECONSTRUCTION::JMuonSimplexParameters_t::TMax_ns

inherited

maximal time w.r.t. Cherenkov hypothesis [ns]

Definition at line 82 of file JMuonSimplexParameters_t.hh.

◆ NMax

int JRECONSTRUCTION::JMuonSimplexParameters_t::NMax

inherited

maximum number of iterations

Definition at line 83 of file JMuonSimplexParameters_t.hh.

◆ reprocess

bool JRECONSTRUCTION::JMuonSimplexParameters_t::reprocess

inherited

reprocess

Definition at line 84 of file JMuonSimplexParameters_t.hh.

The documentation for this struct was generated from the following file:

JMuonSimplex.hh

Public Types

Public Member Functions

Public Attributes

Detailed Description

Member Typedef Documentation

◆ JRegressor_t

◆ hit_type

◆ buffer_type

Constructor & Destructor Documentation

◆ JMuonSimplex()

Member Function Documentation

◆ operator()() [1/2]

◆ operator()() [2/2]

◆ reset()

◆ equals()

◆ ClassDef()

Member Data Documentation

◆ router

◆ sigma_ns

◆ useL0

◆ numberOfPrefits

◆ TMaxLocal_ns

◆ ctMin

◆ roadWidth_m

◆ TMin_ns

◆ TMax_ns

◆ NMax

◆ reprocess