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:

Public Member Functions
	JShowerPointSimplex (const JShowerPointSimplexParameters_t &parameters, const JModuleRouter &router, const int debug=0)
	Parameterized constructor.
JEvt	operator() (const KM3NETDAQ::JDAQEvent &event, const JFIT::JEvt &in)
	Declaration of the member function that actually performs the reconstruction.
void	reset ()
	Reset fit parameters.
bool	equals (const JShowerPointSimplexParameters_t &parameters) const
	Equality.
	ClassDef (JShowerPointSimplexParameters_t, 2)

Public Attributes
const JModuleRouter &	router
size_t	numberOfPrefits
	number of prefits
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
int	NMax
	maximum number of iterations

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

Detailed Description

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

Definition at line 68 of file JShowerPointSimplex.hh.

Member Typedef Documentation

JRegressor_t

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

private

Definition at line 75 of file JShowerPointSimplex.hh.

hit_type

JTRIGGER::JHitR1 JRECONSTRUCTION::JShowerPointSimplex::hit_type

private

Definition at line 77 of file JShowerPointSimplex.hh.

buffer_type

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

private

Definition at line 78 of file JShowerPointSimplex.hh.

Constructor & Destructor Documentation

JShowerPointSimplex()

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

inline

Parameterized constructor.

Parameters

parameters	struct that holds default-optimized parameters for the reconstruction.
router	module router, this is built via detector file.
debug	debug

Definition at line 89 of file JShowerPointSimplex.hh.

                                                                          :
      JShowerPointSimplexParameters_t(parameters),
      JRegressor_t(parameters.sigma_ns),
      router(router)
    {
      using namespace JPP;
      
      JRegressor_t::debug  = debug;
      JRegressor_t::MAXIMUM_ITERATIONS = NMax;
      JRegressor_t::EPSILON = 1e-4;
      //this->estimator.reset(new JMEstimatorTukey(tukey_std_dev));
      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

operator()()

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

inline

Declaration of the member function that actually performs the reconstruction.

Parameters

event	which is a JDAQEvent
in	input fits

Definition at line 116 of file JShowerPointSimplex.hh.

    {
      using namespace std;
      using namespace JPP;
        
      const JBuildL0<hit_type> buildL0;
      const JBuildL2<hit_type> buildL2(JL2Parameters(2, TMaxLocal_ns, ctMin));
      
      buffer_type dataL0;
      buffer_type dataL1;
      
      JEvt out;
        
      buildL0(JDAQTimeslice(event, true), router, back_inserter(dataL0));
      buildL2(JDAQTimeslice(event, false), router, back_inserter(dataL1));  // false = triggered  
 
      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());
          
          JShower3E sh_fit(this->value.getPosition(),
                           JDirection3D(),
                           this->value.getT(),
                           shower->getE());
          out.push_back(getFit(JHistory(shower->getHistory()).add(JSHOWERPOINTSIMPLEX), sh_fit, getQuality(chi2, NDF), NDF, sh_fit.getE())); 
        }
      }
        
    
    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
      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 54 of file JShowerPointSimplexParameters_t.hh.

    {
      return (this->numberOfPrefits  == parameters.numberOfPrefits  &&          
              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 , 2  )

inherited

Member Data Documentation

router

const JModuleRouter& JRECONSTRUCTION::JShowerPointSimplex::router

Definition at line 175 of file JShowerPointSimplex.hh.

numberOfPrefits

size_t JRECONSTRUCTION::JShowerPointSimplexParameters_t::numberOfPrefits

inherited

number of prefits

Definition at line 70 of file JShowerPointSimplexParameters_t.hh.

sigma_ns

double JRECONSTRUCTION::JShowerPointSimplexParameters_t::sigma_ns

inherited

time resolution [ns]

Definition at line 71 of file JShowerPointSimplexParameters_t.hh.

TMaxLocal_ns

double JRECONSTRUCTION::JShowerPointSimplexParameters_t::TMaxLocal_ns

inherited

time window for local coincidences [ns]

Definition at line 72 of file JShowerPointSimplexParameters_t.hh.

TMax_ns

double JRECONSTRUCTION::JShowerPointSimplexParameters_t::TMax_ns

inherited

maximum time for local coincidences [ns]

Definition at line 73 of file JShowerPointSimplexParameters_t.hh.

TMin_ns

double JRECONSTRUCTION::JShowerPointSimplexParameters_t::TMin_ns

inherited

minimum time for local coincidences [ns]

Definition at line 74 of file JShowerPointSimplexParameters_t.hh.

DMax_m

double JRECONSTRUCTION::JShowerPointSimplexParameters_t::DMax_m

inherited

maximal distance to optical module [m]

Definition at line 75 of file JShowerPointSimplexParameters_t.hh.

ctMin

double JRECONSTRUCTION::JShowerPointSimplexParameters_t::ctMin

inherited

minimal cosine space angle between PMT axes

Definition at line 76 of file JShowerPointSimplexParameters_t.hh.

mestimator

int JRECONSTRUCTION::JShowerPointSimplexParameters_t::mestimator

inherited

m-estimator

Definition at line 77 of file JShowerPointSimplexParameters_t.hh.

NMax

int JRECONSTRUCTION::JShowerPointSimplexParameters_t::NMax

inherited

maximum number of iterations

Definition at line 78 of file JShowerPointSimplexParameters_t.hh.

---

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

• JShowerPointSimplex.hh