JRECONSTRUCTION::JMuonFeatures Struct Reference

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

#include <JMuonFeatures.hh>

Inheritance diagram for JRECONSTRUCTION::JMuonFeatures:

Classes
struct	input_type
	Input data type. More...

Public Types
typedef JRegressor< JLine3Z, JGandalf >	JRegressor_t
typedef JHitW0	hit_type
typedef std::vector< hit_type >	buffer_type

Public Member Functions
	JMuonFeatures (const JMuonGandalfParameters_t &parameters, const storage_type &storage, 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 JMuonGandalfParameters_t &parameters) const
	Equality.
	ClassDef (JMuonGandalfParameters_t, 3)

Public Attributes
double	roadWidth_m
	road width [m]
double	R_Hz
	default rate [Hz]
size_t	numberOfPrefits
	number of prefits
double	TTS_ns
	transition-time spread [ns]
double	E_GeV
	energy [GeV]
double	TMin_ns
	minimal time w.r.t. Cherenkov hypothesis [ns]
double	TMax_ns
	maximal time w.r.t. Cherenkov hypothesis [ns]
double	ZMin_m
	minimal z-positon [m]
double	ZMax_m
	maximal z-positon [m]
double	VMax_npe
	maximum number of of photo-electrons
double	cosLR
	maximal cosine space angle likelihood ratio test
int	NMax
	maximum number of iterations

Detailed Description

Wrapper class to add features after the final fit of muon trajectory.

JMuonFeatures computes value of number of hists, number of doms and number of lines from the hits of an event. Note that no values computed at previous steps are modified.

Definition at line 75 of file JMuonFeatures.hh.

Member Typedef Documentation

JRegressor_t

JRegressor<JLine3Z, JGandalf> JRECONSTRUCTION::JMuonFeatures::JRegressor_t

Definition at line 79 of file JMuonFeatures.hh.

hit_type

JHitW0 JRECONSTRUCTION::JMuonFeatures::hit_type

Definition at line 80 of file JMuonFeatures.hh.

buffer_type

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

Definition at line 81 of file JMuonFeatures.hh.

Constructor & Destructor Documentation

JMuonFeatures()

JRECONSTRUCTION::JMuonFeatures::JMuonFeatures ( const JMuonGandalfParameters_t & parameters, const storage_type & storage, const int debug = 0 )

inline

Constructor.

Parameters

parameters	parameters
storage	storage
debug	debug

Definition at line 126 of file JMuonFeatures.hh.

                                                             :
      JMuonGandalfParameters_t(parameters),
      JRegressor_t(storage)
    {
      using namespace JFIT;
      
      if (this->getRmax() < roadWidth_m) {
        roadWidth_m = this->getRmax();
      }
 
      JRegressor_t::debug              = debug;
      JRegressor_t::T_ns.setRange(TMin_ns, TMax_ns);
      JRegressor_t::Vmax_npe           = VMax_npe;
      JRegressor_t::MAXIMUM_ITERATIONS = NMax;
      
      this->parameters.resize(5);
      
      this->parameters[0] = JLine3Z::pX();
      this->parameters[1] = JLine3Z::pY();
      this->parameters[2] = JLine3Z::pT();
      this->parameters[3] = JLine3Z::pDX();
      this->parameters[4] = JLine3Z::pDY();
        
    }

Member Function Documentation

getInput()

input_type JRECONSTRUCTION::JMuonFeatures::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 163 of file JMuonFeatures.hh.

    {
      using namespace std;
      using namespace JTRIGGER;
 
      const JBuildL0<JHitL0>  buildL0;
 
      input_type     input(event.getDAQEventHeader(), in, coverage);
 
      vector<JHitL0> data;
      
      buildL0(event, router, true, back_inserter(data));
 
      for (const auto& hit : data) {
        input.data.push_back(hit_type(hit, summary.getRate(hit.getPMTIdentifier(), this->R_Hz)));
        input.mapString[hit.getModuleID()] = router.getModule(hit.getModuleID()).getString();
      }
      
      return input;
    }

operator()()

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

inline

Fit function.

Parameters

input

input data

Returns

fit results

Definition at line 195 of file JMuonFeatures.hh.

    {
      using namespace std;
      using namespace JFIT;
      using namespace JGEOMETRY3D;
 
      JEvent event(JMUONFEATURES);
      
      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 track = in.begin(); track != in.end(); ++track) {
 
        const JRotation3D  R (getDirection(*track));
        const JLine1Z      tz(getPosition (*track).rotate(R), track->getT());
        JRange<double>     Z_m;
      
        if (track->hasW(JSTART_LENGTH_METRES)       && 
            track->getW(JSTART_LENGTH_METRES) > 0.0) {
          Z_m = JZRange(ZMin_m, ZMax_m + track->getW(JSTART_LENGTH_METRES));
        }
 
        const JModel<JLine1Z> match(tz, roadWidth_m, JRegressor_t::T_ns, Z_m);
 
        // hit selection based on start value
 
        buffer_type buffer;
        
        for (buffer_type::const_iterator i = data.begin(); i != data.end(); ++i) {
 
          hit_type hit(*i);
            
          hit.rotate(R);
 
          if (match(hit)) {
            buffer.push_back(hit);
          }
        }
        
        // select first hit
 
        std::map<int, int> mapString = input.mapString;
        
        std::set<int> stringSet;
        for (buffer_type::const_iterator i = buffer.begin(); i != buffer.end(); ++i) {
          auto it = mapString.find((*i).getModuleID());
          stringSet.insert(it->second);
        }
        
        const int number_of_hits = JLANG::getCount(JLANG::make_array(buffer.begin(), buffer.end(), &JHitW0::getPMTIdentifier));
        const int number_of_doms = JLANG::getCount(JLANG::make_array(buffer.begin(), buffer.end(), &JHitW0::getModuleIdentifier));
        const int number_of_lines = stringSet.size();
 
        const int NDF = number_of_hits - this->parameters.size();
        
        if (NDF > 0) {
 
          out.push_back(JFit(*track).add(JMUONFEATURES));
 
          // set additional values
 
          out.rbegin()->setW(JMUONFEATURES_NUMBER_OF_HITS      , number_of_hits);
          out.rbegin()->setW(JMUONFEATURES_NUMBER_OF_DOMS      , number_of_doms);
          out.rbegin()->setW(JMUONFEATURES_NUMBER_OF_LINES     , number_of_lines);
          
        }
      }
 
      // apply default sorter
      
      sort(out.begin(), out.end(), qualitySorter);
 
      copy(input.in.begin(), input.in.end(), back_inserter(out));
      
      return out;
    }

reset()

void JRECONSTRUCTION::JMuonGandalfParameters_t::reset ( )

inlineinherited

Reset fit parameters.

Definition at line 41 of file JMuonGandalfParameters_t.hh.

    {
      roadWidth_m     = std::numeric_limits<double>::max();
      R_Hz            = 6.0e3;
      numberOfPrefits = 0;
      TTS_ns          = 2;
      E_GeV           = 1.0e3;
      TMin_ns         =  -50.0;     
      TMax_ns         = +450.0;     
      ZMin_m          =    0.0;     
      ZMax_m          =    0.0;
      VMax_npe        =   10.0;
      cosLR           =    0.7;
      NMax            = 1000;
    }

equals()

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

inlineinherited

Equality.

Parameters

parameters

fit parameters

Returns

true if equals; else false

Definition at line 63 of file JMuonGandalfParameters_t.hh.

    {
      return (this->roadWidth_m      == parameters.roadWidth_m     &&
              this->R_Hz             == parameters.R_Hz            &&
              this->numberOfPrefits  == parameters.numberOfPrefits &&
              this->TTS_ns           == parameters.TTS_ns          &&
              this->E_GeV            == parameters.E_GeV           &&
              this->TMin_ns          == parameters.TMin_ns         &&
              this->TMax_ns          == parameters.TMax_ns         &&
              this->ZMin_m           == parameters.ZMin_m          &&
              this->ZMax_m           == parameters.ZMax_m          &&
              this->VMax_npe         == parameters.VMax_npe        &&
              this->cosLR            == parameters.cosLR           &&
              this->NMax             == parameters.NMax);
    }

ClassDef()

JRECONSTRUCTION::JMuonGandalfParameters_t::ClassDef ( JMuonGandalfParameters_t , 3  )

inherited

Member Data Documentation

roadWidth_m

double JRECONSTRUCTION::JMuonGandalfParameters_t::roadWidth_m

inherited

road width [m]

Definition at line 81 of file JMuonGandalfParameters_t.hh.

R_Hz

double JRECONSTRUCTION::JMuonGandalfParameters_t::R_Hz

inherited

default rate [Hz]

Definition at line 82 of file JMuonGandalfParameters_t.hh.

numberOfPrefits

size_t JRECONSTRUCTION::JMuonGandalfParameters_t::numberOfPrefits

inherited

number of prefits

Definition at line 83 of file JMuonGandalfParameters_t.hh.

TTS_ns

double JRECONSTRUCTION::JMuonGandalfParameters_t::TTS_ns

inherited

transition-time spread [ns]

Definition at line 84 of file JMuonGandalfParameters_t.hh.

E_GeV

double JRECONSTRUCTION::JMuonGandalfParameters_t::E_GeV

inherited

energy [GeV]

Definition at line 85 of file JMuonGandalfParameters_t.hh.

TMin_ns

double JRECONSTRUCTION::JMuonGandalfParameters_t::TMin_ns

inherited

minimal time w.r.t. Cherenkov hypothesis [ns]

Definition at line 86 of file JMuonGandalfParameters_t.hh.

TMax_ns

double JRECONSTRUCTION::JMuonGandalfParameters_t::TMax_ns

inherited

maximal time w.r.t. Cherenkov hypothesis [ns]

Definition at line 87 of file JMuonGandalfParameters_t.hh.

ZMin_m

double JRECONSTRUCTION::JMuonGandalfParameters_t::ZMin_m

inherited

minimal z-positon [m]

Definition at line 88 of file JMuonGandalfParameters_t.hh.

ZMax_m

double JRECONSTRUCTION::JMuonGandalfParameters_t::ZMax_m

inherited

maximal z-positon [m]

Definition at line 89 of file JMuonGandalfParameters_t.hh.

VMax_npe

double JRECONSTRUCTION::JMuonGandalfParameters_t::VMax_npe

inherited

maximum number of of photo-electrons

Definition at line 90 of file JMuonGandalfParameters_t.hh.

cosLR

double JRECONSTRUCTION::JMuonGandalfParameters_t::cosLR

inherited

maximal cosine space angle likelihood ratio test

Definition at line 91 of file JMuonGandalfParameters_t.hh.

NMax

int JRECONSTRUCTION::JMuonGandalfParameters_t::NMax

inherited

maximum number of iterations

Definition at line 92 of file JMuonGandalfParameters_t.hh.

---

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

• JMuonFeatures.hh