JRECONSTRUCTION::JMuonSimplex Struct Reference

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

#include <JMuonSimplex.hh>

Inheritance diagram for JRECONSTRUCTION::JMuonSimplex:

Classes
struct	input_type
	Input data type. More...

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 int debug=0)
	Constructor.
input_type	getInput (const JModuleRouter &router, const JDAQEvent &event, const JEvt &in, const coverage_type &coverage) const
	Get input data.
JEvt	operator() (const input_type &input)
	Fit function.
void	reset ()
	Reset fit parameters.
bool	equals (const JMuonSimplexParameters_t &parameters) const
	Equality.
	ClassDef (JMuonSimplexParameters_t, 3)

Public Attributes
double	sigma_ns
	time resolution [ns]
bool	useL0
	option for L0 hit use
size_t	numberOfPrefits
	number of prefits
double	TMaxLocal_ns
	time window for local coincidences [ns]
double	ctMin
	minimal cosine space angle between PMT axes
double	roadWidth_m
	road width [m]
double	TMin_ns
	minimal time w.r.t. Cherenkov hypothesis [ns]
double	TMax_ns
	maximal time w.r.t. Cherenkov hypothesis [ns]
int	NMax
	maximum number of iterations

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 70 of file JMuonSimplex.hh.

Member Typedef Documentation

JRegressor_t

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

Definition at line 74 of file JMuonSimplex.hh.

hit_type

JTRIGGER::JHitR1 JRECONSTRUCTION::JMuonSimplex::hit_type

Definition at line 75 of file JMuonSimplex.hh.

buffer_type

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

Definition at line 76 of file JMuonSimplex.hh.

Constructor & Destructor Documentation

JMuonSimplex()

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

inline

Constructor.

Parameters

parameters	parameters
debug	debug

Definition at line 121 of file JMuonSimplex.hh.

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

Member Function Documentation

getInput()

input_type JRECONSTRUCTION::JMuonSimplex::getInput ( const JModuleRouter & router, const JDAQEvent & event, const JEvt & in, const coverage_type & coverage ) const

inline

Get input data.

Parameters

router	module router
event	event
in	start values
coverage	coverage

Returns

input data

Definition at line 144 of file JMuonSimplex.hh.

    {
      using namespace std;
      using namespace JTRIGGER;
      using namespace KM3NETDAQ;
 
      const JBuildL0<hit_type> buildL0;
      const JBuildL2<hit_type> buildL2(JL2Parameters(2, TMaxLocal_ns, ctMin));
 
      input_type   input(event.getDAQEventHeader(), in, coverage);
 
      buffer_type& dataL0 = input.dataL0;
      buffer_type& dataL1 = input.dataL1;
 
      const 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 input;
    }

operator()()

JEvt JRECONSTRUCTION::JMuonSimplex::operator() ( const input_type & input )

inline

Fit function.

Parameters

input

input data

Returns

fit results

Definition at line 189 of file JMuonSimplex.hh.

    {
      using namespace std;
      using namespace JFIT;
      using namespace JGEOMETRY3D;
      using namespace JLANG;
 
      JEvent event(JMUONSIMPLEX);
 
      JEvt out;
 
      const buffer_type& dataL0 = input.dataL0;
      const buffer_type& dataL1 = input.dataL1;
 
      // select start values
 
      JEvt in = input.in;
 
      in.select(numberOfPrefits, qualitySorter);
 
      if (!in.empty()) {
        in.select(JHistory::is_event(in.begin()->getHistory()));
      }
 
      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(), event()), tb, getQuality(chi2, NDF), NDF)); 
 
          // set additional values
 
          out.rbegin()->setW(track->getW());
          out.rbegin()->setW(JPP_COVERAGE_ORIENTATION, input.coverage.orientation);
          out.rbegin()->setW(JPP_COVERAGE_POSITION,    input.coverage.position);
        }
      }
 
      // apply default sorter
 
      sort(out.begin(), out.end(), qualitySorter);
 
      copy(input.in.begin(), input.in.end(), back_inserter(out));
 
      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;
    }

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

ClassDef()

JRECONSTRUCTION::JMuonSimplexParameters_t::ClassDef ( JMuonSimplexParameters_t , 3  )

inherited

Member Data Documentation

sigma_ns

double JRECONSTRUCTION::JMuonSimplexParameters_t::sigma_ns

inherited

time resolution [ns]

Definition at line 73 of file JMuonSimplexParameters_t.hh.

useL0

bool JRECONSTRUCTION::JMuonSimplexParameters_t::useL0

inherited

option for L0 hit use

Definition at line 74 of file JMuonSimplexParameters_t.hh.

numberOfPrefits

size_t JRECONSTRUCTION::JMuonSimplexParameters_t::numberOfPrefits

inherited

number of prefits

Definition at line 75 of file JMuonSimplexParameters_t.hh.

TMaxLocal_ns

double JRECONSTRUCTION::JMuonSimplexParameters_t::TMaxLocal_ns

inherited

time window for local coincidences [ns]

Definition at line 76 of file JMuonSimplexParameters_t.hh.

ctMin

double JRECONSTRUCTION::JMuonSimplexParameters_t::ctMin

inherited

minimal cosine space angle between PMT axes

Definition at line 77 of file JMuonSimplexParameters_t.hh.

roadWidth_m

double JRECONSTRUCTION::JMuonSimplexParameters_t::roadWidth_m

inherited

road width [m]

Definition at line 78 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 79 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 80 of file JMuonSimplexParameters_t.hh.

NMax

int JRECONSTRUCTION::JMuonSimplexParameters_t::NMax

inherited

maximum number of iterations

Definition at line 81 of file JMuonSimplexParameters_t.hh.

---

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

• JMuonSimplex.hh