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:

Classes
struct	input_type
	Input data type. More...

Public Member Functions
	JShowerPositionFit (const JShowerPositionFitParameters_t &parameters, const storage_type &storage, const JPMTParametersMap &pmtParameters, 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 JShowerPositionFitParameters_t &parameters) const
	Equality.
	ClassDef (JShowerPositionFitParameters_t, 2)

Public Attributes
const JPMTParametersMap &	pmtParameters
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 JHitW0	hit_type
typedef std::vector< hit_type >	buffer_type
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 82 of file JShowerPositionFit.hh.

Member Typedef Documentation

hit_type

JHitW0 JRECONSTRUCTION::JShowerPositionFit::hit_type

private

Definition at line 88 of file JShowerPositionFit.hh.

buffer_type

std::vector<hit_type> JRECONSTRUCTION::JShowerPositionFit::buffer_type

private

Definition at line 89 of file JShowerPositionFit.hh.

JRegressor_t

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

private

Definition at line 90 of file JShowerPositionFit.hh.

Constructor & Destructor Documentation

JShowerPositionFit()

JRECONSTRUCTION::JShowerPositionFit::JShowerPositionFit ( const JShowerPositionFitParameters_t & parameters, const storage_type & storage, const JPMTParametersMap & pmtParameters, const int debug = 0 )

inline

Parameterized constructor.

Parameters

parameters	struct that holds default-optimized parameters for the reconstruction
storage	storage
pmtParameters	PMT parameters
debug	debug

Definition at line 135 of file JShowerPositionFit.hh.

                                                                        :
      JShowerPositionFitParameters_t(parameters),
      JRegressor_t(storage),
      pmtParameters(pmtParameters)
    {
      using namespace JPP;
      
      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;
      JRegressor_t::debug              = debug;
      
      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

getInput()

input_type JRECONSTRUCTION::JShowerPositionFit::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 180 of file JShowerPositionFit.hh.

    { 
      using namespace std;
      using namespace JPP;
      using namespace JTRIGGER;
      
      const JBuildL0<JHitL0> buildL0;
 
      input_type     input(event.getDAQEventHeader(), in, coverage);
 
      vector<JHitL0> data;
 
      buildL0(JDAQTimeslice(event, true), router, back_inserter(data));
 
      for (const auto& hit : data) {
 
        const JPMTIdentifier id(hit.getModuleID(), hit.getPMTAddress());
 
        const JPMTParameters& wip = pmtParameters.getPMTParameters(id);
 
        const int    type = wip.getType();
        const double QE   = wip.QE;
        const double R_Hz = summary.getRate(hit.getPMTIdentifier(), this->R_Hz);
 
        input.data.push_back(hit_type(hit, type, QE, R_Hz));
      }
 
      return input;
    }

operator()()

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

inline

Fit function.

Parameters

input

input data

Returns

fit results

Definition at line 219 of file JShowerPositionFit.hh.

                                            {
 
      using namespace std;
      using namespace JFIT;
      using namespace JGEOMETRY3D;
 
      JEvent event(JSHOWERPOSITIONFIT);
      JEvt out;
 
      const buffer_type& data = input.data;
 
      // select start values
 
      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) {
 
        JPoint4D vx(getPosition(*shower), shower->getT());
 
        const JFIT::JModel<JPoint4D> match(vx, DMax_m, JRegressor_t::T_ns);
 
        buffer_type buffer;
        
        for (buffer_type::const_iterator i = data.begin(); i != data.end(); ++i) {
  
          if (match(*i)) {
            buffer.push_back(*i);
          }
        }
 
        // select first hit
        
        sort(buffer.begin(), buffer.end(), JHitL0::compare);
        
        vector<hit_type>::iterator __end = unique(buffer.begin(), buffer.end(), equal_to<JDAQPMTIdentifier>());
 
        const int NDF = distance(buffer.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, buffer.begin(), __end);
          
            JShower3D result(JVertex3D(this->value.getPosition(), this->value.getT()), JDirection3D());
  
         
            out.push_back(getFit(JHistory(shower->getHistory(), event()), result, getQuality(chi2), NDF, this->value.getE()));
 
            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(JSHOWERPOSITIONFIT_ISOTROPIC_ENERGY, this->value.getE());  
            out.rbegin()->setW(JSHOWERPOSITIONFIT_LAMBDA,           Y.GetNrows() != 0 ? Y[0] : 0.0);
 
            // set additional values
            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::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 93 of file JShowerPositionFit.hh.

pmtParameters

const JPMTParametersMap& JRECONSTRUCTION::JShowerPositionFit::pmtParameters

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