JRECONSTRUCTION::JShowerPositionFit Class Reference

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.
JEvt	operator() (const KM3NETDAQ::JDAQEvent &event, const JFIT::JEvt &in)
	Declaration of the member function that actually performs the reconstruction.
void	reset ()
	Reset fit parameters.
bool	equals (const JShowerPositionFitParameters_t &parameters) const
	Equality.
	ClassDef (JShowerPositionFitParameters_t, 2)

Public Attributes
const JModuleRouter &	router
const JSummaryRouter &	summary
size_t	numberOfPrefits
	number of prefits
double	TMax_ns
	maximum time for local coincidences [ns]
double	TMin_ns
	minimum time for local coincidences [ns]
double	DMax_m
	maximal distance to optical module [m]
double	Emin_GeV
	minimum energy to scan
double	Emax_GeV
	maximum energy to scan
int	En
	number of points to scan in energy range
double	R_Hz
	default rate [Hz]
double	TTS_ns
	transition-time spread [ns]
double	VMax_npe
	maximum number of of photo-electrons
int	NMax
	maximum number of iterations

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

JRegressor_t

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

private

Definition at line 80 of file JShowerPositionFit.hh.

Constructor & Destructor Documentation

JShowerPositionFit()

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(5);
       
      this->parameters[0] = JPoint4E::pX();
      this->parameters[1] = JPoint4E::pY();
      this->parameters[2] = JPoint4E::pZ();
      this->parameters[3] = JPoint4E::pT();
      this->parameters[4] = JPoint4E::pE();
    }

Member Function Documentation

operator()()

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 137 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 e = Ev.begin(); e != Ev.end(); ++e) {
            JPoint4E sh(vx, *e);
            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;
    }

reset()

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

equals()

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

ClassDef()

JRECONSTRUCTION::JShowerPositionFitParameters_t::ClassDef ( JShowerPositionFitParameters_t , 2  )

inherited

Member Data Documentation

Ev

std::vector<double> JRECONSTRUCTION::JShowerPositionFit::Ev

private

Definition at line 82 of file JShowerPositionFit.hh.

router

const JModuleRouter& JRECONSTRUCTION::JShowerPositionFit::router

Definition at line 196 of file JShowerPositionFit.hh.

summary

const JSummaryRouter& JRECONSTRUCTION::JShowerPositionFit::summary

Definition at line 197 of file JShowerPositionFit.hh.

numberOfPrefits

size_t JRECONSTRUCTION::JShowerPositionFitParameters_t::numberOfPrefits

inherited

number of prefits

Definition at line 70 of file JShowerPositionFitParameters_t.hh.

TMax_ns

double JRECONSTRUCTION::JShowerPositionFitParameters_t::TMax_ns

inherited

maximum time for local coincidences [ns]

Definition at line 71 of file JShowerPositionFitParameters_t.hh.

TMin_ns

double JRECONSTRUCTION::JShowerPositionFitParameters_t::TMin_ns

inherited

minimum time for local coincidences [ns]

Definition at line 72 of file JShowerPositionFitParameters_t.hh.

DMax_m

double JRECONSTRUCTION::JShowerPositionFitParameters_t::DMax_m

inherited

maximal distance to optical module [m]

Definition at line 73 of file JShowerPositionFitParameters_t.hh.

Emin_GeV

double JRECONSTRUCTION::JShowerPositionFitParameters_t::Emin_GeV

inherited

minimum energy to scan

Definition at line 74 of file JShowerPositionFitParameters_t.hh.

Emax_GeV

double JRECONSTRUCTION::JShowerPositionFitParameters_t::Emax_GeV

inherited

maximum energy to scan

Definition at line 75 of file JShowerPositionFitParameters_t.hh.

En

int JRECONSTRUCTION::JShowerPositionFitParameters_t::En

inherited

number of points to scan in energy range

Definition at line 76 of file JShowerPositionFitParameters_t.hh.

R_Hz

double JRECONSTRUCTION::JShowerPositionFitParameters_t::R_Hz

inherited

default rate [Hz]

Definition at line 77 of file JShowerPositionFitParameters_t.hh.

TTS_ns

double JRECONSTRUCTION::JShowerPositionFitParameters_t::TTS_ns

inherited

transition-time spread [ns]

Definition at line 78 of file JShowerPositionFitParameters_t.hh.

VMax_npe

double JRECONSTRUCTION::JShowerPositionFitParameters_t::VMax_npe

inherited

maximum number of of photo-electrons

Definition at line 79 of file JShowerPositionFitParameters_t.hh.

NMax

int JRECONSTRUCTION::JShowerPositionFitParameters_t::NMax

inherited

maximum number of iterations

Definition at line 80 of file JShowerPositionFitParameters_t.hh.

---

The documentation for this class was generated from the following file:

• JShowerPositionFit.hh