Wrapper class to make final fit of muon trajectory. More...

#include <JMuonGandalf.hh>

Inheritance diagram for JRECONSTRUCTION::JMuonGandalf:

Public Types
typedef JRegressor< JLine3Z, JGandalf >	JRegressor_t

typedef JTRIGGER::JHitL0	hit_type

typedef std::vector< hit_type >	buffer_type

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

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

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

void	reset ()
	Reset fit parameters. More...

bool	equals (const JMuonGandalfParameters_t &parameters) const
	Equality. More...

	ClassDef (JMuonGandalfParameters_t, 2)

Public Attributes
const JModuleRouter &	router

const JSummaryRouter &	summary

int	debug

double	roadWidth_m
	road width [m] More...

double	R_Hz
	default rate [Hz] More...

size_t	numberOfPrefits
	number of prefits More...

double	TTS_ns
	transition-time spread [ns] More...

double	E_GeV
	energy [GeV] More...

double	TMin_ns
	minimal time w.r.t. Cherenkov hypothesis [ns] More...

double	TMax_ns
	maximal time w.r.t. Cherenkov hypothesis [ns] More...

double	ZMin_m
	minimal z-positon [m] More...

double	ZMax_m
	maximal z-positon [m] More...

double	VMax_npe
	maximum number of of photo-electrons More...

int	NMax
	maximum number of iterations More...

bool	reprocess
	reprocess More...

Detailed Description

Wrapper class to make 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 72 of file JMuonGandalf.hh.

Member Typedef Documentation

◆ JRegressor_t

typedef JRegressor<JLine3Z, JGandalf> JRECONSTRUCTION::JMuonGandalf::JRegressor_t

Definition at line 76 of file JMuonGandalf.hh.

◆ hit_type

typedef JTRIGGER::JHitL0 JRECONSTRUCTION::JMuonGandalf::hit_type

Definition at line 77 of file JMuonGandalf.hh.

◆ buffer_type

typedef std::vector<hit_type> JRECONSTRUCTION::JMuonGandalf::buffer_type

Definition at line 78 of file JMuonGandalf.hh.

Constructor & Destructor Documentation

◆ JMuonGandalf()

JRECONSTRUCTION::JMuonGandalf::JMuonGandalf	(	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 91 of file JMuonGandalf.hh.

                                                             :
       JMuonGandalfParameters_t(parameters),
       JRegressor_t(pdf_file, 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::JMuonGandalf::operator()	(	const KM3NETDAQ::JDAQEvent &	event,
		const JEvt &	in
	)

inline

Fit function.

Parameters

event	event
in	start values

Returns: fit results

Definition at line 129 of file JMuonGandalf.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::JMuonGandalf::operator()	(	const buffer_type &	data,
		const JEvt &	in
	)

inline

Fit function.

Parameters

data	hit data
in	start values

Returns: fit results

Definition at line 152 of file JMuonGandalf.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()));
             
           hit.rotate(R);
  
           if (match(hit)) {
             buffer.push_back(hit);
           }
         }
           
         // select first hit
           
         sort(buffer.begin(), buffer.end(), JHitL0::compare);
           
         vector<JHitW0>::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
  
           if (track->getE() > 0.1)
             JRegressor_t::E_GeV   = track->getE();
           else
             JRegressor_t::E_GeV   = this->JMuonGandalfParameters_t::E_GeV;
  
           const double chi2 = (*this)(JLine3Z(tz), buffer.begin(), __end);
  
           // check error matrix
  
           bool status = true;
  
           for (size_t i = 0; i != this->V.size(); ++i) {
             if (std::isnan(this->V(i,i)) || this->V(i,i) < 0.0) {
               status = false;
             }
           }
  
           if (status) {
  
             JTrack3D tb(this->value);
             
             tb.rotate_back(R);
  
             out.push_back(getFit(JHistory(track->getHistory()).add(JMUONGANDALF), tb, getQuality(chi2), NDF));
  
             // set additional values
  
             out.rbegin()->setV(this->V.size(), this->V);
  
             out.rbegin()->setW(track->getW());
             out.rbegin()->setW(JGANDALF_BETA0_RAD,            sqrt(this->error.getDX() * this->error.getDX()  +
                                                                    this->error.getDY() * this->error.getDY()));
             out.rbegin()->setW(JGANDALF_BETA1_RAD,            sqrt(this->error.getDX() * this->error.getDY()));
             out.rbegin()->setW(JGANDALF_CHI2,                 chi2);
             out.rbegin()->setW(JGANDALF_NUMBER_OF_HITS,       distance(buffer.begin(), __end));
             out.rbegin()->setW(JGANDALF_LAMBDA,               this->lambda);
             out.rbegin()->setW(JGANDALF_NUMBER_OF_ITERATIONS, this->numberOfIterations);
           }
         }
       }
  
       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;
       reprocess       = false;
     }

◆ 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->NMax             == parameters.NMax            &&
               this->reprocess        == parameters.reprocess);
     }

◆ ClassDef()

JRECONSTRUCTION::JMuonGandalfParameters_t::ClassDef	(	JMuonGandalfParameters_t	,
		2
	)

inherited

Member Data Documentation

◆ router

const JModuleRouter& JRECONSTRUCTION::JMuonGandalf::router

Definition at line 246 of file JMuonGandalf.hh.

◆ summary

const JSummaryRouter& JRECONSTRUCTION::JMuonGandalf::summary

Definition at line 247 of file JMuonGandalf.hh.

◆ debug

int JRECONSTRUCTION::JMuonGandalf::debug

Definition at line 248 of file JMuonGandalf.hh.

◆ 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.

◆ NMax

int JRECONSTRUCTION::JMuonGandalfParameters_t::NMax

inherited

maximum number of iterations

Definition at line 91 of file JMuonGandalfParameters_t.hh.

◆ reprocess

bool JRECONSTRUCTION::JMuonGandalfParameters_t::reprocess

inherited

reprocess

Definition at line 92 of file JMuonGandalfParameters_t.hh.

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

JMuonGandalf.hh

Public Types

Public Member Functions

Public Attributes

Detailed Description

Member Typedef Documentation

◆ JRegressor_t

◆ hit_type

◆ buffer_type

Constructor & Destructor Documentation

◆ JMuonGandalf()

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

◆ reprocess