JRECONSTRUCTION::JShowerPointSimplex Class Reference

class to handle the second position fit of the shower reconstruction, mainly dedicated for ORCA More...

#include <JShowerPointSimplex.hh>

Inheritance diagram for JRECONSTRUCTION::JShowerPointSimplex:

Classes
struct	input_type
	Input data type. More...

Public Types
typedef JRegressor< JPoint4D, JGandalf >	JRegressor_t
typedef JTRIGGER::JHitR1	hit_type
typedef std::vector< hit_type >	buffer_type

Public Member Functions
	JShowerPointSimplex (const JShowerPointSimplexParameters_t &parameters, const int debug=0)
	Parameterized constructor.
input_type	getInput (const JModuleRouter &router, const KM3NETDAQ::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 JShowerPointSimplexParameters_t &parameters) const
	Equality.
	ClassDef (JShowerPointSimplexParameters_t, 3)

Public Attributes
size_t	numberOfPrefits
	number of prefits
size_t	numberOfPostfits
	number of postfits
double	Qwatershed
	quality threshold
double	sigma_ns
	time resolution [ns]
double	TMaxLocal_ns
	time window for local coincidences [ns]
double	TMax_ns
	maximum time for local coincidences [ns]
double	TMin_ns
	minimum time for local coincidences [ns]
double	DMax_m
	maximal distance to optical module [m]
double	ctMin
	minimal cosine space angle between PMT axes
int	mestimator
	m-estimator (see JFIT::JMEstimator_t)
int	NMax
	maximum number of iterations

Detailed Description

class to handle the second position fit of the shower reconstruction, mainly dedicated for ORCA

Definition at line 73 of file JShowerPointSimplex.hh.

Member Typedef Documentation

JRegressor_t

JRegressor<JPoint4D, JGandalf> JRECONSTRUCTION::JShowerPointSimplex::JRegressor_t

Definition at line 79 of file JShowerPointSimplex.hh.

hit_type

JTRIGGER::JHitR1 JRECONSTRUCTION::JShowerPointSimplex::hit_type

Definition at line 80 of file JShowerPointSimplex.hh.

buffer_type

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

Definition at line 81 of file JShowerPointSimplex.hh.

Constructor & Destructor Documentation

JShowerPointSimplex()

JRECONSTRUCTION::JShowerPointSimplex::JShowerPointSimplex ( const JShowerPointSimplexParameters_t & parameters, const int debug = 0 )

inline

Parameterized constructor.

Parameters

parameters	struct that holds default-optimized parameters for the reconstruction.
debug	debug

Definition at line 125 of file JShowerPointSimplex.hh.

                                                                           :
      JShowerPointSimplexParameters_t(parameters),
      JRegressor_t(parameters.sigma_ns)
    {
      using namespace JPP;
      
      JRegressor_t::debug  = debug;
      JRegressor_t::MAXIMUM_ITERATIONS = NMax;
      JRegressor_t::EPSILON = 1e-4;
 
      this->estimator.reset(getMEstimator(parameters.mestimator));
      this->parameters.resize(4);
      this->parameters[0] = JPoint4D::pX();
      this->parameters[1] = JPoint4D::pY();
      this->parameters[2] = JPoint4D::pZ();
      this->parameters[3] = JPoint4D::pT();
    }

Member Function Documentation

getInput()

input_type JRECONSTRUCTION::JShowerPointSimplex::getInput ( const JModuleRouter & router, const KM3NETDAQ::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 152 of file JShowerPointSimplex.hh.

    {
      using namespace std;
      using namespace KM3NETDAQ;
      using namespace JTRIGGER; 
      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;
  
      buildL0(JDAQTimeslice(event, true), router, back_inserter(dataL0));
      buildL2(JDAQTimeslice(event, false), router, back_inserter(dataL1));  // false = triggered  
 
      return input;
    }

operator()()

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

inline

Fit function.

Parameters

input

input data

Returns

fit results

Definition at line 177 of file JShowerPointSimplex.hh.

    {
      using namespace std;
      using namespace JPP;
 
      JEvent event(JSHOWERPOINTSIMPLEX);
 
      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()));
      }
      
      for (JEvt::const_iterator shower = in.begin(); shower != in.end(); ++shower) {
 
        JPoint4D vx(getPosition(*shower), shower->getT());
 
        JFIT::JModel<JPoint4D> match(vx, DMax_m, JTimeRange(TMin_ns, TMax_ns));
 
        buffer_type data;
 
        data.reserve(dataL0.size() + dataL1.size());
 
        for (buffer_type::const_iterator i = dataL1.begin(); i != dataL1.end(); ++i) {
          
          if (match(*i)) {
            data.push_back(*i);
          }
        }
 
        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) {
            if (match(*i)) {
              data.push_back(*i);
            }
          }
        }
        
        const int NDF = getCount(data.begin(), data.end()) - this->parameters.size();
        
        if (NDF > 0) {
 
          double chi2 = (*this)(vx, data.begin(), data.end());
          
          JShower3D result(JVertex3D(this->value.getPosition(), this->value.getT()), JDirection3D());
 
          out.push_back(getFit(JHistory(shower->getHistory(), event()), result, getQuality(chi2, NDF), NDF, shower->getE()));
 
          // set additional values
 
          out.rbegin()->setW(JPP_COVERAGE_ORIENTATION, input.coverage.orientation);
          out.rbegin()->setW(JPP_COVERAGE_POSITION,    input.coverage.position);
        }
      }
      if (numberOfPostfits > 0) {
 
        JEvt::iterator __end = out.end();
 
        if (Qwatershed > 0) {                         // sort distinct maxima
 
          __end = out.begin();                        // output
 
          for (JEvt::iterator p1 = out.begin(); p1 != out.end() && (size_t) distance(out.begin(), __end) < numberOfPostfits; ) {
 
            JEvt::iterator p2 = p1;                   // p2 -> best-of-the-rest
 
            for (JEvt::iterator i = p1; i != out.end(); ++i) {
              if (qualitySorter(*i, *p2)) {
                p2 = i;
              }
            }
 
            swap(*p2, *p1);                           // put best-of-the-rest at head
 
            p2 = p1++;                                // p2 = head; p1 = next to head
 
            const JPoint point(*p2);             // reference position
 
            sort(p1, out.end(), point);               // sort according 4D distance with respect to reference position
 
            for (double Q = p2->getQ();
                 p1 != out.end() && (p1->getQ() >= p2->getQ() - Qwatershed ||   // bottom
                                     p1->getQ() <= Q);                          // descend
                 ++p1) {
              Q = p1->getQ();
            }
 
            swap(*(__end++), *p2);                    // output head
          }
          
        } else {
 
          advance(__end = out.begin(), min(numberOfPostfits, out.size()));
 
          partial_sort(out.begin(), __end, out.end(), qualitySorter);
        }
 
        out.erase(__end, out.end());
 
      } else {
 
        sort(out.begin(), out.end(), qualitySorter);
      }
 
      copy(input.in.begin(), input.in.end(), back_inserter(out));
 
      return out;
    }

reset()

void JRECONSTRUCTION::JShowerPointSimplexParameters_t::reset ( )

inlineinherited

Reset fit parameters.

Definition at line 35 of file JShowerPointSimplexParameters_t.hh.

    {
      numberOfPrefits  =     0;  // 0 means process all inputs
      numberOfPostfits =     0;  // 0 means do not select
      Qwatershed       =     0;  // 0 means no watershed algorithm
      sigma_ns         =     1;
      TMaxLocal_ns     =    10;
      TMax_ns          =    60.0;  
      TMin_ns          =   -60.0;   
      DMax_m           =    50.0;
      ctMin            =     0.2;
      mestimator       =  JFIT::EM_LORENTZIAN;
      NMax             =  1000;
    }

equals()

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

inlineinherited

Equality.

Parameters

parameters

fit parameters

Returns

true if equals; else false

Definition at line 56 of file JShowerPointSimplexParameters_t.hh.

    {
      return (this->numberOfPrefits  == parameters.numberOfPrefits  &&          
              this->numberOfPostfits == parameters.numberOfPostfits &&
              this->Qwatershed       == parameters.Qwatershed       &&
              this->sigma_ns         == parameters.sigma_ns         &&
              this->TMaxLocal_ns     == parameters.TMaxLocal_ns     &&
              this->TMax_ns          == parameters.TMax_ns          &&
              this->TMin_ns          == parameters.TMin_ns          &&
              this->DMax_m           == parameters.DMax_m           &&
              this->ctMin            == parameters.ctMin            &&
              this->mestimator       == parameters.mestimator       &&
              this->NMax             == parameters.NMax
              );
    }

ClassDef()

JRECONSTRUCTION::JShowerPointSimplexParameters_t::ClassDef ( JShowerPointSimplexParameters_t , 3  )

inherited

Member Data Documentation

numberOfPrefits

size_t JRECONSTRUCTION::JShowerPointSimplexParameters_t::numberOfPrefits

inherited

number of prefits

Definition at line 74 of file JShowerPointSimplexParameters_t.hh.

numberOfPostfits

size_t JRECONSTRUCTION::JShowerPointSimplexParameters_t::numberOfPostfits

inherited

number of postfits

Definition at line 75 of file JShowerPointSimplexParameters_t.hh.

Qwatershed

double JRECONSTRUCTION::JShowerPointSimplexParameters_t::Qwatershed

inherited

quality threshold

Definition at line 76 of file JShowerPointSimplexParameters_t.hh.

sigma_ns

double JRECONSTRUCTION::JShowerPointSimplexParameters_t::sigma_ns

inherited

time resolution [ns]

Definition at line 77 of file JShowerPointSimplexParameters_t.hh.

TMaxLocal_ns

double JRECONSTRUCTION::JShowerPointSimplexParameters_t::TMaxLocal_ns

inherited

time window for local coincidences [ns]

Definition at line 78 of file JShowerPointSimplexParameters_t.hh.

TMax_ns

double JRECONSTRUCTION::JShowerPointSimplexParameters_t::TMax_ns

inherited

maximum time for local coincidences [ns]

Definition at line 79 of file JShowerPointSimplexParameters_t.hh.

TMin_ns

double JRECONSTRUCTION::JShowerPointSimplexParameters_t::TMin_ns

inherited

minimum time for local coincidences [ns]

Definition at line 80 of file JShowerPointSimplexParameters_t.hh.

DMax_m

double JRECONSTRUCTION::JShowerPointSimplexParameters_t::DMax_m

inherited

maximal distance to optical module [m]

Definition at line 81 of file JShowerPointSimplexParameters_t.hh.

ctMin

double JRECONSTRUCTION::JShowerPointSimplexParameters_t::ctMin

inherited

minimal cosine space angle between PMT axes

Definition at line 82 of file JShowerPointSimplexParameters_t.hh.

mestimator

int JRECONSTRUCTION::JShowerPointSimplexParameters_t::mestimator

inherited

m-estimator (see JFIT::JMEstimator_t)

Definition at line 83 of file JShowerPointSimplexParameters_t.hh.

NMax

int JRECONSTRUCTION::JShowerPointSimplexParameters_t::NMax

inherited

maximum number of iterations

Definition at line 84 of file JShowerPointSimplexParameters_t.hh.

---

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

• JShowerPointSimplex.hh