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

Public Attributes
const JModuleRouter &	router

const JSummaryRouter &	summary

size_t	numberOfPrefits
	number of prefits 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	Emin_GeV
	minimum energy to scan More...

double	Emax_GeV
	maximum energy to scan More...

int	En
	number of points to scan in energy range 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< JPoint4E, JGandalf >	JRegressor_t

Private Attributes
std::vector< double >	Ev

Detailed Description

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

Definition at line 73 of file JShowerPositionFit.hh.

Member Typedef Documentation

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

private

Definition at line 80 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 94 of file JShowerPositionFit.hh.

                                                                         :
       JShowerPositionFitParameters_t(parameters),
       JRegressor_t(pdfFile, parameters.TTS_ns),
       router(router),
       summary(summary)
     {
       using namespace JPP;
       
       JRegressor_t::debug  = debug;
       JRegressor_t::T_ns.setRange(parameters.TMin_ns, parameters.TMax_ns);
       JRegressor_t::Vmax_npe = VMax_npe;
       JRegressor_t::MAXIMUM_ITERATIONS = NMax;
       JRegressor_t::EPSILON = 1e-3;
       
       if (Emin_GeV > Emax_GeV || En <= 1) {
         THROW(JException, "Invalid energy input " << Emin_GeV << ' ' << Emax_GeV << ' ' << En);
       }
 
       const double base = std::pow((Emax_GeV / Emin_GeV), 1.0 / (En - 1));
 
       for (int i = 0; i != En; ++i) {
         Ev.push_back(Emin_GeV * std::pow(base, i));
       }
       
       this->parameters.resize(4);
        
       this->parameters[0] = JPoint4E::pX();
       this->parameters[1] = JPoint4E::pY();
       this->parameters[2] = JPoint4E::pZ();
       this->parameters[3] = JPoint4E::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 136 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, DMax_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
           for (vector<double>::iterator i = Ev.begin(); i != Ev.end(); ++i) {
             JPoint4E sh(vx, *i);
             double chi2 = (*this)(sh, data.begin(), __end);
           
             JShower3E sh_fit(this->value.getPosition(), JDirection3D(),
                              this->value.getT(), this->value.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 
       TMax_ns          =   100.0;  
       TMin_ns          =  -100.0;   
       DMax_m           =    80.0;
       R_Hz             = 10000;
       Emin_GeV         = 1;
       Emax_GeV         = 200;
       En               = 10;
       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 53 of file JShowerPositionFitParameters_t.hh.

     {
       return (this->numberOfPrefits  == parameters.numberOfPrefits  &&         
               this->TMax_ns          == parameters.TMax_ns          &&
               this->TMin_ns          == parameters.TMin_ns          &&
               this->DMax_m           == parameters.DMax_m           &&
               this->R_Hz             == parameters.R_Hz             &&
               this->TTS_ns           == parameters.TTS_ns           &&
               this->Emin_GeV         == parameters.Emin_GeV         &&
               this->Emax_GeV         == parameters.Emax_GeV         &&
               this->En               == parameters.En               &&
               this->VMax_npe         == parameters.VMax_npe         &&
               this->NMax             == parameters.NMax);
     }