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:

Public Member Functions
	JShowerFit (const JShowerFitParameters_t &parameters, const JModuleRouter &router, const JSummaryRouter &summary, const std::string pdfFile, const JShowerEnergyCorrection &correct, 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 JShowerFitParameters_t &parameters) const
	Equality. More...
	ClassDef (JShowerFitParameters_t, 2)

Public Attributes
const JModuleRouter &	router
const JSummaryRouter &	summary
const JShowerEnergyCorrection &	correct
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	DStep_m
	step increase for the distance to optical module [m] More...
double	R_Hz
	default rate [Hz] More...
int	mestimator
	M-estimator More...
double	Vmax_npe
	maximum number of of photo-electrons More...
double	DMax_m
	maximal distance to optical module [m] More...

Private Types
typedef JRegressor< JShower3EZ, JGandalf >	JRegressor_t

Detailed Description

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

Definition at line 72 of file JShowerFit.hh.

Member Typedef Documentation

JRegressor_t

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

private

Definition at line 77 of file JShowerFit.hh.

Constructor & Destructor Documentation

JShowerFit()

JRECONSTRUCTION::JShowerFit::JShowerFit ( const JShowerFitParameters_t &  parameters, const JModuleRouter &  router, const JSummaryRouter &  summary, const std::string  pdfFile, const JShowerEnergyCorrection &  correct, const int  debug = 0  )

inline

Parameterized constructor.

Parameters

parameters	struct that holds default-optimized parameters for the reconstruction, available in $JPP_DATA.
router	module router, this is built via detector file.
summary	summary router
pdfFile	PDF file
correct	energy correction
debug	debug

Definition at line 92 of file JShowerFit.hh.

                                                         :
      JShowerFitParameters_t(parameters),
      JRegressor_t(pdfFile),
      router(router),
      summary(summary),
      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

operator()()

JEvt JRECONSTRUCTION::JShowerFit::operator() ( const KM3NETDAQ::JDAQEvent &  event, const JFIT::JEvt &  in  )

inline

Declaration of the member function that actually performs the reconstruction.

Parameters

event	= JDAQEvent
in	= input fits

Definition at line 128 of file JShowerFit.hh.

    {
      using namespace std;
      using namespace JPP;
 
      typedef vector<JHitL0>  JDataL0_t;
      JBuildL0<JHitL0>        buildL0;
 
      JEvt out;
 
      const JDetector &detector = router.getReference();
 
      JDataL0_t dataL0;
 
      buildL0(JDAQTimeslice(event, true), router, back_inserter(dataL0));
  
      for (JEvt::const_iterator shower = in.begin(); shower != in.end(); ++shower) {
 
        const JVertex3D Vertex(getPosition(*shower), shower->getT());
        
        multiset<JDAQPMTIdentifier> top;
 
        const double distance = DMax_m  +  DStep_m * log10(shower->getE());
 
        const JFIT::JModel<JPoint4D> match(Vertex, distance, JRegressor_t::T_ns);
 
        for (JDataL0_t::const_iterator i = dataL0.begin(); i != dataL0.end(); ++i) {
          if (match(*i)) {
            top.insert(i->getPMTIdentifier());
          }
        }
        
        const JDetectorSubset_t subdetector(detector, Vertex.getPosition(), distance);
            
        const JRotation3D R(getDirection(*shower));
        
        vector<JPMTW0> buffer;
          
        for (JDetectorSubset_t::const_iterator module = subdetector.begin();
             module != subdetector.end(); ++module) {
            
          const JDAQSummaryFrame& frame = summary.getSummaryFrame(module->getID());
 
          JModule dom(*module); 
 
          dom.rotate(R);  
 
          for (size_t i = 0; i != dom.size(); ++i) {
 
            if (getDAQStatus(frame, *module, i)            &&
                getPMTStatus(frame, *module, i)            &&
                frame[i].is_valid()                        &&
                !module->getPMT(i).has(PMT_DISABLE)) {
              const JDAQPMTIdentifier id(module->getID(), i);
                
              const double rate_Hz = summary.getRate(id);
              const size_t count   = top.count(id);
            
              buffer.push_back(JPMTW0(dom.getPMT(i), rate_Hz, count));
            }
          } 
        }
 
        double chi2 = (*this)(JShower3EZ(JVertex3D(JVector3D(0,0,0), shower->getT()), JVersor3Z(), 
                                         shower->getE()), buffer.begin(), buffer.end());
          
        double NDF = getCount(buffer.begin(), buffer.end()) - this->parameters.size(); 
 
        JShower3E sh_fit(this->value.getPosition(), this->value.getDirection(),
                         this->value.getT(), correct(this->value.getE()));
        
        sh_fit.rotate_back(R);
 
        sh_fit.add(Vertex.getPosition());
          
        out.push_back(getFit(JHistory(shower->getHistory()).add(JSHOWERCOMPLETEFIT), sh_fit, getQuality(chi2), 
                             NDF, sh_fit.getE())); 
        out.rbegin()->setV(this->V.size(), this->V);  
        out.rbegin()->setW(JSHOWERFIT_ENERGY,  this->value.getE());  // Uncorrected Energy
      }
 
      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

router

const JModuleRouter& JRECONSTRUCTION::JShowerFit::router

Definition at line 212 of file JShowerFit.hh.

summary

const JSummaryRouter& JRECONSTRUCTION::JShowerFit::summary

Definition at line 213 of file JShowerFit.hh.

correct

const JShowerEnergyCorrection& JRECONSTRUCTION::JShowerFit::correct

Definition at line 214 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

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