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

#include <JMuonFeatures.hh>

Inheritance diagram for JRECONSTRUCTION::JMuonFeatures:

Public Types
typedef JRegressor< JLine3Z, JGandalf >	JRegressor_t

typedef JTRIGGER::JHitL0	hit_type

typedef std::vector< hit_type >	buffer_type

Public Member Functions
	JMuonFeatures (const JMuonGandalfParameters_t &parameters, const JModuleRouter &router, const JSummaryRouter &summary, const std::string &pdf_file, const int debug=0)
	Constructor.

JEvt	operator() (const KM3NETDAQ::JDAQEvent &event, const JEvt &in)
	Fit function.

JEvt	operator() (const buffer_type &data, const JEvt &in)
	Fit function.

void	reset ()
	Reset fit parameters.

bool	equals (const JMuonGandalfParameters_t &parameters) const
	Equality.

	ClassDef (JMuonGandalfParameters_t, 3)

Public Attributes
const JModuleRouter &	router

const JSummaryRouter &	summary

int	debug

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

int	NMax
	maximum number of iterations

Detailed Description

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

The JMuonGandalf fit uses one or more start values (usually taken from the output of JMuonSimplex).
All hits of which the PMT position lies within a set road width (JMuonGandalfParameters_t::roadWidth_m) and time is within a set window (JMuonGandalfParameters_t::TMin_ns, JMuonGandalfParameters_t::TMax_ns) around the Cherenkov hypothesis are taken.
In case there are multiple hits from the same PMT is the specified window, the first hit is taken and the other hits are discarded.
The PDF is accordingly evaluated, i.e. the normalised probability for a first hit at the given time of the hit is taken. The normalisation is consistently based on the specified time window.
Note that this hit selection is unbiased with respect to the PDF of a single PMT.

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

JTRIGGER::JHitL0 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 JModuleRouter &	router,
		const JSummaryRouter &	summary,
		const std::string &	pdf_file,
		const int	debug = 0 )

inline

Constructor.

Parameters

parameters	parameters
router	module router
summary	summary file router
pdf_file	PDF file
debug	debug

Definition at line 94 of file JMuonFeatures.hh.

                                                             :
      JMuonGandalfParameters_t(parameters),
      JRegressor_t(pdf_file, JTimeRange(parameters.TMin_ns, parameters.TMax_ns), parameters.TTS_ns),
      router (router),
      summary(summary)
    {
      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

◆ operator()() [1/2]

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

inline

Fit function.

Parameters

event	event
in	start values

Returns: fit results

Definition at line 132 of file JMuonFeatures.hh.

    {
      using namespace std;
      using namespace JFIT;
      using namespace JTRIGGER;
 
      const JBuildL0<hit_type>  buildL0;
 
      buffer_type dataL0;
 
      buildL0(event, router, true, back_inserter(dataL0));
 
      return (*this)(dataL0, in);
    }

◆ operator()() [2/2]

JEvt JRECONSTRUCTION::JMuonFeatures::operator()	(	const buffer_type &	data,
		const JEvt &	in )

inline

Fit function.

Parameters

data	hit data
in	start values

Returns: fit results

Definition at line 155 of file JMuonFeatures.hh.

    {
      using namespace std;
      using namespace JFIT;
      using namespace JGEOMETRY3D;
 
      JEvt out;
      
      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
 
        vector<JHitW0> buffer;
          
        for (buffer_type::const_iterator i = data.begin(); i != data.end(); ++i) {
 
          JHitW0 hit(*i, summary.getRate(i->getPMTIdentifier(), this->R_Hz));
            
          hit.rotate(R);
 
          if (match(hit)) {
            buffer.push_back(hit);
          }
        }
          
        // select first hit
 
        std::set<int> string_ids;
        for (vector<JHitW0>::const_iterator i = buffer.begin(); i != buffer.end(); ++i) {
          string_ids.insert(router.getModule(((*i).getModuleID())).getString());
        }
        
        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 = string_ids.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);
          
        }
      }
 
      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;
      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 62 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->NMax             == parameters.NMax);
    }

◆ ClassDef()

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

inherited

Member Data Documentation

◆ router

const JModuleRouter& JRECONSTRUCTION::JMuonFeatures::router

Definition at line 221 of file JMuonFeatures.hh.

◆ summary

const JSummaryRouter& JRECONSTRUCTION::JMuonFeatures::summary

Definition at line 222 of file JMuonFeatures.hh.

◆ debug

int JRECONSTRUCTION::JMuonFeatures::debug

Definition at line 223 of file JMuonFeatures.hh.

◆ roadWidth_m

double JRECONSTRUCTION::JMuonGandalfParameters_t::roadWidth_m

inherited

road width [m]

Definition at line 79 of file JMuonGandalfParameters_t.hh.

◆ R_Hz

double JRECONSTRUCTION::JMuonGandalfParameters_t::R_Hz

inherited

default rate [Hz]

Definition at line 80 of file JMuonGandalfParameters_t.hh.

◆ numberOfPrefits

size_t JRECONSTRUCTION::JMuonGandalfParameters_t::numberOfPrefits

inherited

number of prefits

Definition at line 81 of file JMuonGandalfParameters_t.hh.

◆ TTS_ns

double JRECONSTRUCTION::JMuonGandalfParameters_t::TTS_ns

inherited

transition-time spread [ns]

Definition at line 82 of file JMuonGandalfParameters_t.hh.

◆ E_GeV

double JRECONSTRUCTION::JMuonGandalfParameters_t::E_GeV

inherited

energy [GeV]

Definition at line 83 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 84 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 85 of file JMuonGandalfParameters_t.hh.

◆ ZMin_m

double JRECONSTRUCTION::JMuonGandalfParameters_t::ZMin_m

inherited

minimal z-positon [m]

Definition at line 86 of file JMuonGandalfParameters_t.hh.

◆ ZMax_m

double JRECONSTRUCTION::JMuonGandalfParameters_t::ZMax_m

inherited

maximal z-positon [m]

Definition at line 87 of file JMuonGandalfParameters_t.hh.

◆ VMax_npe

double JRECONSTRUCTION::JMuonGandalfParameters_t::VMax_npe

inherited

maximum number of of photo-electrons

Definition at line 88 of file JMuonGandalfParameters_t.hh.

◆ NMax

int JRECONSTRUCTION::JMuonGandalfParameters_t::NMax

inherited

maximum number of iterations

Definition at line 89 of file JMuonGandalfParameters_t.hh.

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

JMuonFeatures.hh

Public Types

Public Member Functions

Public Attributes

Detailed Description

Member Typedef Documentation

◆ JRegressor_t

◆ hit_type

◆ buffer_type

Constructor & Destructor Documentation

◆ JMuonFeatures()

Member Function Documentation

◆ operator()() [1/2]

◆ operator()() [2/2]

◆ reset()

◆ equals()

◆ ClassDef()

Member Data Documentation

◆ router

◆ summary

◆ debug

◆ roadWidth_m

◆ R_Hz

◆ numberOfPrefits

◆ TTS_ns

◆ E_GeV

◆ TMin_ns

◆ TMax_ns

◆ ZMin_m

◆ ZMax_m

◆ VMax_npe

◆ NMax