JRECONSTRUCTION::JShowerFeatures Struct Reference

Wrapper class to add features after the final fit of shower. More...

#include <JShowerFeatures.hh>

Inheritance diagram for JRECONSTRUCTION::JShowerFeatures:

Classes
struct	input_type
	Input data type. More...

Public Types
typedef JRegressor< JShower3EZ, JAbstractMinimiser >	JRegressor_t
typedef JModuleL0	module_type
typedef std::vector< module_type >	detector_type
typedef JRegressorStorage< JShower3EZ, JFIT::JAbstractMinimiser >	storage_type

Public Member Functions
	JShowerFeatures (const JShowerFitParameters_t &parameters, const storage_type &storage, const JPMTParametersMap &pmtParameters, const int debug=0)
	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
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]

Detailed Description

Wrapper class to add features after the final fit of shower.

Note that no values computed at previous steps are modified.

Definition at line 77 of file JShowerFeatures.hh.

Member Typedef Documentation

JRegressor_t

JRegressor<JShower3EZ, JAbstractMinimiser> JRECONSTRUCTION::JShowerFeatures::JRegressor_t

Definition at line 81 of file JShowerFeatures.hh.

module_type

JModuleL0 JRECONSTRUCTION::JShowerFeatures::module_type

Definition at line 83 of file JShowerFeatures.hh.

detector_type

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

Definition at line 84 of file JShowerFeatures.hh.

storage_type

JRegressorStorage<JShower3EZ, JFIT::JAbstractMinimiser> JRECONSTRUCTION::JShowerFeatures::storage_type

Definition at line 85 of file JShowerFeatures.hh.

Constructor & Destructor Documentation

JShowerFeatures()

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

inline

Constructor.

Parameters

parameters	parameters
storage	storage
pmtParameters	PMT parameters
debug	debug

Definition at line 127 of file JShowerFeatures.hh.

                                                             :
      JShowerFitParameters_t(parameters),
      JRegressor_t(storage),
      pmtParameters(pmtParameters)
    {
      using namespace JFIT;
 
      JRegressor_t::debug              = debug;
      JRegressor_t::T_ns.setRange(TMin_ns, TMax_ns);
    }

Member Function Documentation

getInput()

input_type JRECONSTRUCTION::JShowerFeatures::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 151 of file JShowerFeatures.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::JShowerFeatures::operator() ( const input_type & input )

inline

Fit function.

Parameters

input

input data

Returns

fit results

Definition at line 195 of file JShowerFeatures.hh.

    {
      using namespace std;
      using namespace JFIT;
      using namespace JGEOMETRY3D;
 
      JEvent event(JSHOWERFEATURES);
 
      JEvt out;
      
      JEvt in = input.in;
      
      in.select(numberOfPrefits, qualitySorter);
 
      if (!in.empty()) {
        in.select(JHistory::is_event(in.begin()->getHistory()));
      }
 
      // Computation of max pmt getting a hit in a DOM for a given event, snapshot                                                   
 
      struct {
        size_t nhits =  0;
        int    id    = -1;
      } snapshot;
 
      for (const auto& module : input.data) {
        if (module.size() != 0) {
          
          size_t snapshot_hits = 0;
 
          for (size_t i = 0; i != module->size(); ++i) {
            
            if (module.getStatus(i)) {
 
              const size_t nhits = std::count_if(
                                                 module.begin(),
                                                 module.end(),
                                                 [i](const JHitR0& hit) {
                                                   return hit.getPMT() == i;
                                                 }
                                                 );
 
              if (nhits > 0) snapshot_hits++;
            }
          }
          
          if (snapshot_hits > snapshot.nhits){
            snapshot.nhits  = snapshot_hits;
            snapshot.id     = module->getID();
          }
        }
      }
      
      for (JEvt::const_iterator shower = in.begin(); shower != in.end(); ++shower) {
        
        //Computation of the minimal distance from jshower vertex to a module
        struct {
          double Dmin = numeric_limits<double>::max();
          int    id   = -1;
        } vertex;
 
        for (const auto& module : input.data) {
          if (module->getFloor() != 0) {
 
            const double D = getDistance(module->getPosition(), getPosition(*shower));
    
            if (D < vertex.Dmin) {
              
              vertex.Dmin = D;
              vertex.id   = module->getID();
              
            }
          }
        }
      
        out.push_back(JFit(*shower).add(JSHOWERFEATURES));
        out.rbegin()->setW(JSHOWERFEATURES_DMIN,               vertex.Dmin);
        out.rbegin()->setW(JSHOWERFEATURES_DMIN_MODULEID,      vertex.id);
        out.rbegin()->setW(JSHOWERFEATURES_NSNAPSHOT,          snapshot.nhits);
        out.rbegin()->setW(JSHOWERFEATURES_NSNAPSHOT_MODULEID, snapshot.id);
 
      }
 
      // 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::JShowerFeatures::pmtParameters

Definition at line 287 of file JShowerFeatures.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 struct was generated from the following file:

• JShowerFeatures.hh