JRECONSTRUCTION::JShowerFit Class Reference

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

#include <JShowerFit.hh>

Inheritance diagram for JRECONSTRUCTION::JShowerFit:

Classes
struct	input_type
	Input data type. More...

Public Member Functions
	JShowerFit (const JShowerFitParameters_t &parameters, const storage_type &storage, const JPMTParametersMap &pmtParameters, const JEnergyCorrection &correct, const int debug=0)
	Parameterized constructor.
input_type	getInput (const JModuleRouter &router, const JSummaryRouter &summary, const JDAQEvent &event, const JEvt &in, const coverage_type &coverage) const
	Get input data.
JEvt	operator() (const input_type &input)
	Fit function.
void	reset ()
	Reset fit parameters.
bool	equals (const JShowerFitParameters_t &parameters) const
	Equality.
	ClassDef (JShowerFitParameters_t, 2)

Public Attributes
const JPMTParametersMap &	pmtParameters
const JEnergyCorrection &	correct
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	DStep_m
	step increase for the distance to optical module [m]
double	R_Hz
	default rate [Hz]
int	mestimator
	M-estimator (see JFIT::JMEstimator_t)
double	Vmax_npe
	maximum number of of photo-electrons
double	DMax_m
	maximal distance to optical module [m]

Private Types
typedef JRegressor< JShower3EZ, JGandalf >	JRegressor_t
typedef JModuleL0	module_type
typedef std::vector< module_type >	detector_type

Detailed Description

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

Definition at line 65 of file JShowerFit.hh.

Member Typedef Documentation

JRegressor_t

JRegressor<JShower3EZ, JGandalf> JRECONSTRUCTION::JShowerFit::JRegressor_t

private

Definition at line 70 of file JShowerFit.hh.

module_type

JModuleL0 JRECONSTRUCTION::JShowerFit::module_type

private

Definition at line 71 of file JShowerFit.hh.

detector_type

std::vector<module_type> JRECONSTRUCTION::JShowerFit::detector_type

private

Definition at line 72 of file JShowerFit.hh.

Constructor & Destructor Documentation

JShowerFit()

JRECONSTRUCTION::JShowerFit::JShowerFit ( const JShowerFitParameters_t & parameters, const storage_type & storage, const JPMTParametersMap & pmtParameters, const JEnergyCorrection & correct, const int debug = 0 )

inline

Parameterized constructor.

Parameters

parameters	parameters
storage	storage
pmtParameters	PMT parameters
correct	energy correction
debug	debug

Definition at line 120 of file JShowerFit.hh.

                                                       :
      JShowerFitParameters_t(parameters),
      JRegressor_t (storage),
      pmtParameters(pmtParameters),
      correct      (correct)
    {
      using namespace JPP;
      
      JRegressor_t::debug  = debug;
      JRegressor_t::T_ns.setRange(parameters.TMin_ns, parameters.TMax_ns);
      JRegressor_t::Vmax_npe = parameters.Vmax_npe;
      JRegressor_t::MAXIMUM_ITERATIONS = 1000;
      JRegressor_t::EPSILON = 1e-3;
      JRegressor_t::EPSILON_ABSOLUTE = true;
      
      this->parameters.resize(3);
 
      this->parameters[0] = JShower3EZ::pDX();
      this->parameters[1] = JShower3EZ::pDY();
      this->parameters[2] = JShower3EZ::pE();
 
      this->estimator.reset(getMEstimator(parameters.mestimator));
    }

Member Function Documentation

getInput()

input_type JRECONSTRUCTION::JShowerFit::getInput ( const JModuleRouter & router, const JSummaryRouter & summary, const JDAQEvent & event, const JEvt & in, const coverage_type & coverage ) const

inline

Get input data.

Parameters

router	module router
summary	summary data
event	event
in	start values
coverage	coverage

Returns

input data

Definition at line 158 of file JShowerFit.hh.

    {
      using namespace std;
      using namespace JTRIGGER;
 
      input_type   input(event.getDAQEventHeader(), in, coverage);
 
      const JBuildL0 <JHitR0>   buildL0;
      map<int, vector<JHitR0> > data;
 
      const JDAQTimeslice timeslice(event, true);
 
      JSuperFrame2D<JHit> buffer;
 
      for (JDAQTimeslice::const_iterator i = timeslice.begin(); i != timeslice.end(); ++i) {
 
        if (router.hasModule(i->getModuleID())) {
 
          buffer(*i, router.getModule(i->getModuleID()));
 
          buildL0(buffer, back_inserter(data[i->getModuleID()]));
        }
      }
 
      for (const auto& module : router.getReference()) {
        if (!module.empty()) {
          input.data.push_back(module_type(module, summary.getSummaryFrame(module.getID(), R_Hz), data[module.getID()]));
        }
      }
 
      return input;
    }

operator()()

JEvt JRECONSTRUCTION::JShowerFit::operator() ( const input_type & input )

inline

Fit function.

Parameters

input

input data

Returns

fit results

Definition at line 201 of file JShowerFit.hh.

    {
      using namespace std;
      using namespace JPP;
 
      JEvent event(JSHOWERCOMPLETEFIT);
 
      JEvt out;
      
      JEvt in = input.in;
 
      in.select(numberOfPrefits, qualitySorter);
 
      if (!in.empty()) {
        in.select(JHistory::is_event(in.begin()->getHistory()));
      }
 
      for (JEvt::const_iterator shower = in.begin(); shower != in.end(); ++shower) {
 
        vector<JPMTW0> data;
 
        const JPosition3D vertex(getPosition(*shower));
        const double      time = shower->getT();
        const double      distance = DMax_m  +  DStep_m * log10(shower->getE());
        const JRotation3D R(getDirection(*shower));
        
        for (const auto& module : input.data) {
 
          JPosition3D pos(module->getPosition());
 
          pos.sub(vertex);
 
          if (pos.getLength() < distance) {
 
            for (size_t i = 0; i != module->size(); ++i) {
 
              if (module.getStatus(i)) {
 
                const double t1 = time + pos.getLength() * getInverseSpeedOfLight() * getIndexOfRefraction();
 
                struct {
 
                  bool operator()(const JHitR0& hit) const
                  {
                    return (hit.getPMT() == pmt && T_ns(hit.getT()));
                  }
 
                  const JTimeRange T_ns;
                  const size_t     pmt;
 
                } match = { JRegressor_t::T_ns + t1, i };
 
                const JPMTIdentifier id(module->getID(), i);
 
                const JPMTParameters& wip = pmtParameters.getPMTParameters(id);
 
                const size_t ns = count_if(module.begin(), module.end(), match);
                const double QE = wip.QE;
                
                JPMT pmt = module->getPMT(i);
 
                pmt.sub(vertex);
                pmt.rotate(R);
                
                data.push_back(JPMTW0(pmt, QE, module.frame.getRate(i), ns));
              }
            }
          }
        }
            
 
        double chi2 = (*this)(JShower3EZ(JVertex3D(JVector3D(0,0,0), shower->getT()), JVersor3Z(), 
                                         shower->getE()), data.begin(), data.end());
          
        double NDF = getCount(data.begin(), data.end()) - this->parameters.size(); 
 
        JShower3E result(JShower3D(this->value.getVertex(), this->value.getDirection()), correct(this->value.getE()));
 
        // check error matrix
        bool status = true;
 
        for (size_t i = 0; i != this->V.size(); ++i) {
          if (std::isnan(this->V(i,i)) || this->V(i,i) < 0.0) {
            status = false;
          }
        }
        
        if (status) {
 
          result.rotate_back(R);
 
          result.add(vertex.getPosition());
          
 
          out.push_back(getFit(JHistory(shower->getHistory(), event()), result, getQuality(chi2), NDF, result.getE())); 
 
          out.rbegin()->setW(shower->getW());
 
          const size_t N = this->V.size();
 
          TMatrixDSym  M(N);
 
          for (size_t row = 0; row != N; ++row) {
            for (size_t col = 0; col != N; ++col) {
              M(row,col) = this->V(row,col);
            }
          }
 
          const TMatrixDSymEigen E(M);
          const TVectorD&        Y = E.GetEigenValues();
 
          out.rbegin()->setV(this->V.size(), this->V);  
 
          out.rbegin()->setW(JSHOWERFIT_LAMBDA,  Y.GetNrows() != 0 ? Y[0] : 0.0);
          out.rbegin()->setW(JSHOWERFIT_ENERGY,  this->value.getE());  // Uncorrected Energy
          
          out.rbegin()->setW(JPP_COVERAGE_ORIENTATION, input.coverage.orientation);
          out.rbegin()->setW(JPP_COVERAGE_POSITION,    input.coverage.position);
 
        }
      }
      
      // apply default sorter
 
      sort(out.begin(), out.end(), qualitySorter);
 
      copy(input.in.begin(), input.in.end(), back_inserter(out));
      
      return out;
    }

reset()

void JRECONSTRUCTION::JShowerFitParameters_t::reset ( )

inlineinherited

Reset fit parameters.

Definition at line 35 of file JShowerFitParameters_t.hh.

    {
      numberOfPrefits             =   4;
      TMax_ns                     =  30;
      TMin_ns                     = -30;
      DStep_m                     =  0;
      DMax_m                      =  60;
      R_Hz                        =  10.0e3;
      mestimator                  =  JFIT::EM_NULL;
      Vmax_npe                    =  20;
    }

equals()

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

inlineinherited

Equality.

Parameters

parameters

fit parameters

Returns

true if equals; else false

Definition at line 53 of file JShowerFitParameters_t.hh.

    {
      return (this->TMax_ns              == parameters.TMax_ns                  &&         
              this->TMin_ns                  == parameters.TMin_ns                  &&         
              this->numberOfPrefits          == parameters.numberOfPrefits          &&         
              this->DStep_m                  == parameters.DStep_m                  &&
              this->DMax_m                   == parameters.DMax_m                   &&
              this->R_Hz                     == parameters.R_Hz                     && 
              this->mestimator               == parameters.mestimator               &&
              this->Vmax_npe                 == parameters.Vmax_npe                 
              );
    }

ClassDef()

JRECONSTRUCTION::JShowerFitParameters_t::ClassDef ( JShowerFitParameters_t , 2  )

inherited

Member Data Documentation

pmtParameters

const JPMTParametersMap& JRECONSTRUCTION::JShowerFit::pmtParameters

Definition at line 333 of file JShowerFit.hh.

correct

const JEnergyCorrection& JRECONSTRUCTION::JShowerFit::correct

Definition at line 334 of file JShowerFit.hh.

numberOfPrefits

size_t JRECONSTRUCTION::JShowerFitParameters_t::numberOfPrefits

inherited

number of prefits

Definition at line 69 of file JShowerFitParameters_t.hh.

TMax_ns

double JRECONSTRUCTION::JShowerFitParameters_t::TMax_ns

inherited

maximum time for local coincidences [ns]

Definition at line 70 of file JShowerFitParameters_t.hh.

TMin_ns

double JRECONSTRUCTION::JShowerFitParameters_t::TMin_ns

inherited

minimum time for local coincidences [ns]

Definition at line 71 of file JShowerFitParameters_t.hh.

DStep_m

double JRECONSTRUCTION::JShowerFitParameters_t::DStep_m

inherited

step increase for the distance to optical module [m]

Definition at line 72 of file JShowerFitParameters_t.hh.

R_Hz

double JRECONSTRUCTION::JShowerFitParameters_t::R_Hz

inherited

default rate [Hz]

Definition at line 73 of file JShowerFitParameters_t.hh.

mestimator

int JRECONSTRUCTION::JShowerFitParameters_t::mestimator

inherited

M-estimator (see JFIT::JMEstimator_t)

Definition at line 74 of file JShowerFitParameters_t.hh.

Vmax_npe

double JRECONSTRUCTION::JShowerFitParameters_t::Vmax_npe

inherited

maximum number of of photo-electrons

Definition at line 75 of file JShowerFitParameters_t.hh.

DMax_m

double JRECONSTRUCTION::JShowerFitParameters_t::DMax_m

inherited

maximal distance to optical module [m]

Definition at line 76 of file JShowerFitParameters_t.hh.

---

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

• JShowerFit.hh