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
	number of prefits More...

size_t	numberOfOutfits
	number of fits to be saved in output 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...

double	E_GeV
	energy [GeV] More...

double	R_Hz
	default rate [Hz] More...

double	TTS_ns
	transition-time spread [ns] More...

double	VMax_npe
	maximum number of of photo-electrons More...

int	NMax
	maximum number of iterations 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 = VMax_npe;
       JRegressor_t::MAXIMUM_ITERATIONS = NMax;
 
       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 JFIT::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(), JHitL0::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;
       VMax_npe         =    20.0;
       NMax             =  1000;
     }

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 52 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           &&
               this->VMax_npe         == parameters.VMax_npe         &&
               this->NMax             == parameters.NMax);
     }