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. More...

JEvt	operator() (const KM3NETDAQ::JDAQEvent &event, const JFIT::JEvt &in)
	Declaration of the member function that actually performs the reconstruction. More...

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

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

	ClassDef (JShowerPointSimplexParameters_t, 2)

Public Attributes
const JModuleRouter &	router

size_t	numberOfPrefits
	number of prefits More...

double	sigma_ns
	time resolution [ns] More...

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

double	TMax_ns
	maximum time for local coincidences [ns] More...

double	TMin_ns
	minimum time for local coincidences [ns] More...

double	DMax_m
	maximal distance to optical module [m] More...

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

int	mestimator
	m-estimator More...

int	NMax
	maximum number of iterations More...

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

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

private

Definition at line 75 of file JShowerPointSimplex.hh.

typedef JTRIGGER::JHitR1 JRECONSTRUCTION::JShowerPointSimplex::hit_type

private

Definition at line 77 of file JShowerPointSimplex.hh.

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

private

Definition at line 78 of file JShowerPointSimplex.hh.

Constructor & Destructor Documentation

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

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

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

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