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

#include <JShowerPositionFit.hh>

Inheritance diagram for JRECONSTRUCTION::JShowerPositionFit:

Public Member Functions
	JShowerPositionFit (const JShowerPositionFitParameters_t &parameters, const JModuleRouter &router, const JSummaryRouter &summary, const std::string pdfFile, 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 JShowerPositionFitParameters_t &parameters) const
	Equality. More...

	ClassDef (JShowerPositionFitParameters_t, 1)

Public Attributes
const JModuleRouter &	router

const JSummaryRouter &	summary

size_t	numberOfPrefits

size_t	numberOfOutfits

double	TMax_ns

double	TMin_ns

double	ctMin

double	roadWidth_m

double	E_GeV

double	R_Hz

double	TTS_ns

Static Public Attributes
struct {
}	compare
	Auxiliary data structure for sorting of hits. More...

Private Types
typedef JRegressor< JPoint4D, JGandalf >	JRegressor_t

Detailed Description

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

Definition at line 72 of file JShowerPositionFit.hh.

Member Typedef Documentation

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

private

Definition at line 79 of file JShowerPositionFit.hh.

Constructor & Destructor Documentation

JRECONSTRUCTION::JShowerPositionFit::JShowerPositionFit	(	const JShowerPositionFitParameters_t &	parameters,
		const JModuleRouter &	router,
		const JSummaryRouter &	summary,
		const std::string	pdfFile,
		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.
summary	summary router
pdfFile	file name with PDF
debug	debug

Definition at line 93 of file JShowerPositionFit.hh.

                                                                         :
       JShowerPositionFitParameters_t(parameters),
       JRegressor_t(pdfFile, parameters.TTS_ns),
       router(router),
       summary(summary)
     {
       using namespace JPP;
       
       JRegressor_t::debug  = 1;
       JRegressor_t::T_ns.setRange(parameters.TMin_ns, parameters.TMax_ns);
       JRegressor_t::Vmax_npe = 20.0;
       JRegressor_t::MAXIMUM_ITERATIONS = 1000;
 
       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::JShowerPositionFit::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 124 of file JShowerPositionFit.hh.

     {
       using namespace std;
       using namespace JPP;
         
       typedef vector<JHitL0> JDataL0_t;
       const JBuildL0<JHitL0> buildL0;
         
       JEvt out;
         
       JDataL0_t dataL0;
       buildL0(JDAQTimeslice(event, true), router, back_inserter(dataL0));
   
       for (JEvt::const_iterator shower = in.begin(); shower != in.end(); ++shower) {
 
         JPoint4D vx(getPosition(*shower), shower->getT());
 
         vector<JHitW0> data;
 
         const JModel<JPoint4D> match(vx, roadWidth_m, JRegressor_t::T_ns);
 
         for (JDataL0_t::const_iterator i = dataL0.begin(); i != dataL0.end(); ++i) {
 
           const double rate_Hz = summary.getRate(*i);
           
           if (match(*i)) {
             data.push_back(JHitW0(*i, rate_Hz));
           }
         }
 
         // select first hit
         
         sort(data.begin(), data.end(), compare);
         
         vector<JHitW0>::iterator __end = unique(data.begin(), data.end(), equal_to<JDAQPMTIdentifier>());
 
         const int NDF = distance(data.begin(), __end) - this->parameters.size();
 
         if (NDF > 0) {
         
           // set fit parameters
           
           if (shower->getE() > 0.1)
             JRegressor_t::E_GeV   = shower->getE();
           else
             JRegressor_t::E_GeV   = this->JShowerPositionFitParameters_t::E_GeV;
 
           double chi2 = (*this)(vx, data.begin(), __end);
           
           JShower3E sh_fit(this->value.getPosition(), JDirection3D(),
                            this->value.getT(), shower->getE());
   
           out.push_back(getFit(JHistory(shower->getHistory()).add(JSHOWERPOSITIONFIT), sh_fit, 
                                getQuality(chi2), NDF, sh_fit.getE())); 
           
         }    
       }
       
       return out;
     }

void JRECONSTRUCTION::JShowerPositionFitParameters_t::reset ( )

inlineinherited

Reset fit parameters.

Definition at line 32 of file JShowerPositionFitParameters_t.hh.

     {
       numberOfPrefits  =     0;  // process all inputs from PointSimplex, which should give a max of 100 
       numberOfOutfits  =     5;
       TMax_ns          =   100.0;  
       TMin_ns          =  -100.0;   
       roadWidth_m      =    80.0;
       E_GeV            =    10;
       R_Hz             = 10000;   
       TTS_ns           =     2;
     }

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

inlineinherited

Equality.

Parameters

parameters fit parameters

Returns: true if equals; else false

Definition at line 50 of file JShowerPositionFitParameters_t.hh.

     {
       return (this->numberOfPrefits  == parameters.numberOfPrefits  &&         
               this->numberOfOutfits  == parameters.numberOfOutfits  &&
               this->TMax_ns          == parameters.TMax_ns          &&
               this->TMin_ns          == parameters.TMin_ns          &&
               this->ctMin            == parameters.ctMin            &&
               this->roadWidth_m      == parameters.roadWidth_m      &&
               this->E_GeV            == parameters.E_GeV            &&
               this->R_Hz             == parameters.R_Hz             &&
               this->TTS_ns           == parameters.TTS_ns);
     }