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, 1)

Public Attributes
const JModuleRouter &	router

size_t	numberOfPrefits
	number of prefits More...

size_t	numberOfOutfits
	number of fits to be saved in output More...

size_t	minPrefitsSize
	number of prefits under which an event is treated as a very low energy one 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	roadWidth_m
	road width [m] More...

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

JTOOLS::JRange< double >	TWindow_ns
	time window for local coincidences [ns] for events with few prefits More...

JTOOLS::JRange< double >	pos_grid_m
	edges in [m] of the position grid More...

double	pos_step_m
	step in [m] of position grid More...

JTOOLS::JRange< double >	time_grid_ns
	edges in [ns] of the time grid More...

double	time_step_ns
	step in [ns] of time grid More...

double	simplex_step_m
	step in [m] of JSimplex minimiser More...

double	simplex_step_ns
	step in [ns] of JSimplex minimiser More...

double	tukey_std_dev
	standard deviation of Tukey estimator More...

int	NMax
	maximum number of iterations More...

Private Types
typedef JRegressor< JPoint4D, JSimplex >	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 66 of file JShowerPointSimplex.hh.

Member Typedef Documentation

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

private

Definition at line 73 of file JShowerPointSimplex.hh.

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

private

Definition at line 75 of file JShowerPointSimplex.hh.

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

private

Definition at line 76 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 87 of file JShowerPointSimplex.hh.

                                                                           :
       JShowerPointSimplexParameters_t(parameters),
       JRegressor_t(parameters.sigma_ns),
       router(router)
     {
       using namespace JPP;
       
       JRegressor_t::debug  = 1;
       JRegressor_t::MAXIMUM_ITERATIONS = NMax;
 
       this->estimator.reset(new JMEstimatorTukey(tukey_std_dev));
     }

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 108 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());
 
         buffer_type data;
   
         data.reserve(dataL0.size() + dataL1.size());
 
         JFIT::JModel<JPoint4D> match(vx, roadWidth_m, JTimeRange(TMin_ns, TMax_ns));
 
         JRange<double>   posGrid_m(0, 0);
         JTimeRange       timeGrid_ns(0, 0);
 
         /*
          * part to help the reco of events with few hits:
          * in this case a bigger time window in the hit selection and
          * a grid in position and time are considered.
          */
         if(in.size() <= minPrefitsSize){
           
           match = JFIT::JModel<JPoint4D>(vx, roadWidth_m, TWindow_ns);
 
           posGrid_m   = pos_grid_m;
           timeGrid_ns = time_grid_ns;
         }
 
         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, bind2nd(equal_to<JDAQModuleIdentifier>(), i->getModuleID())) == __end) {
  
             if (match(*i)) {
               data.push_back(*i);
             }
           }
         }
        
         const int NDF = getCount(data.begin(), data.end()) - this->step.size();
 
         if (NDF > 0) {
 
           for(double x = posGrid_m.getLowerLimit(); x <= posGrid_m.getUpperLimit(); x += pos_step_m){
             for(double y = posGrid_m.getLowerLimit(); y <= posGrid_m.getUpperLimit(); y += pos_step_m){
               for(double z = posGrid_m.getLowerLimit(); z <= posGrid_m.getUpperLimit(); z += pos_step_m){
                 for(double t = timeGrid_ns.getLowerLimit(); t <= timeGrid_ns.getUpperLimit(); t += time_step_ns){
                   
                   JPoint4D vxs(JVector3D(vx.getX() + x, vx.getY() + y, vx.getZ() + z), vx.getT() + t);
               
                   this->step.resize(4);
                   this->step[0] = JPoint4D(JVector3D(simplex_step_m, 0.0, 0.0), 0.0);
                   this->step[1] = JPoint4D(JVector3D(0.0, simplex_step_m, 0.0), 0.0);
                   this->step[2] = JPoint4D(JVector3D(0.0, 0.0, simplex_step_m), 0.0);
                   this->step[3] = JPoint4D(JVector3D(), simplex_step_ns);
 
                   double chi2 = (*this)(vxs, 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 34 of file JShowerPointSimplexParameters_t.hh.

     {
       using namespace JPP;
       
       numberOfPrefits  =     0;  // 0 means process all inputs
       numberOfOutfits  =   100;
       minPrefitsSize   =     8;
       sigma_ns         =     1;
       TMaxLocal_ns     =    10;
       TMax_ns          =    60.0;  
       TMin_ns          =   -60.0;   
       roadWidth_m      =    50.0;
       ctMin            =     0.2;
       TWindow_ns       =  JRange<double>(-130, 120);
       pos_grid_m       =  JRange<double>( -20,  20);
       pos_step_m       =    20;
       time_grid_ns     =  JRange<double>( -90,  90);
       time_step_ns     =    20;
       simplex_step_m   =     5;
       simplex_step_ns  =    10;
       tukey_std_dev    =    20;
       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 64 of file JShowerPointSimplexParameters_t.hh.

     {
       return (this->numberOfPrefits  == parameters.numberOfPrefits  &&         
               this->numberOfOutfits  == parameters.numberOfOutfits  &&         
               this->minPrefitsSize   == parameters.minPrefitsSize   &&
               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->roadWidth_m      == parameters.roadWidth_m      &&
               this->ctMin            == parameters.ctMin            &&
               this->TWindow_ns       == parameters.TWindow_ns       &&
               this->pos_grid_m       == parameters.pos_grid_m       &&
               this->pos_step_m       == parameters.pos_step_m       &&
               this->time_grid_ns     == parameters.time_grid_ns     &&
               this->time_step_ns     == parameters.time_step_ns     &&
               this->simplex_step_m   == parameters.simplex_step_m   &&
               this->simplex_step_ns  == parameters.simplex_step_ns  &&
               this->tukey_std_dev    == parameters.tukey_std_dev    &&
               this->NMax             == parameters.NMax);
     }