JRECONSTRUCTION::JMuonGandalf Struct Reference

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

#include <JMuonGandalf.hh>

Inheritance diagram for JRECONSTRUCTION::JMuonGandalf:

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
	JMuonGandalf (const JMuonGandalfParameters_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 JMuonGandalfParameters_t &parameters) const
	Equality.
	ClassDef (JMuonGandalfParameters_t, 3)

Public Attributes
const JPMTParametersMap &	pmtParameters
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 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 82 of file JMuonGandalf.hh.

Member Typedef Documentation

JRegressor_t

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

Definition at line 86 of file JMuonGandalf.hh.

hit_type

JHitW0 JRECONSTRUCTION::JMuonGandalf::hit_type

Definition at line 87 of file JMuonGandalf.hh.

buffer_type

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

Definition at line 88 of file JMuonGandalf.hh.

Constructor & Destructor Documentation

JMuonGandalf()

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

                                                            :
      JMuonGandalfParameters_t(parameters),
      JRegressor_t(storage),
      pmtParameters(pmtParameters)
    {
      if (this->getRmax() < roadWidth_m) {
        roadWidth_m = this->getRmax();
      }
      
      JRegressor_t::debug              = debug;
      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::JMuonGandalf::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 170 of file JMuonGandalf.hh.

    {
      using namespace std;
      using namespace JPP;
 
      const JBuildL0<JHitL0>  buildL0;
 
      input_type     input(event.getDAQEventHeader(), in, coverage);
 
      vector<JHitL0> data;
 
      buildL0(event, router, true, back_inserter(data));
 
      for (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::JMuonGandalf::operator() ( const input_type & input )

inline

Fit function.

Parameters

input

input data

Returns

fit results

Definition at line 210 of file JMuonGandalf.hh.

    {
      using namespace std;
      using namespace JFIT;
      using namespace JGEOMETRY3D;
 
      JEvent event(JMUONGANDALF);
 
      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->getW(JSTART_LENGTH_METRES, 0.0) > 0.0) {
          Z_m = JZRange(track->getW(JSTART_ZMIN_M) + ZMin_m,
                        track->getW(JSTART_ZMAX_M) + ZMax_m);
        }
 
        const JModel<JLine1Z> match(tz, roadWidth_m, 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
          
        sort(buffer.begin(), buffer.end(), JHitL0::compare);
          
        buffer_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
 
          if (this->JMuonGandalfParameters_t::E_GeV > 0.1)
            JRegressor_t::E_GeV  =  this->JMuonGandalfParameters_t::E_GeV;
          else
            JRegressor_t::E_GeV  =  track->getE();
 
          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(), event()), tb, getQuality(chi2), NDF));
 
            // set additional values
 
            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(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_NUMBER_OF_HITS,       distance(buffer.begin(), __end));
            out.rbegin()->setW(JGANDALF_LAMBDA,               Y.GetNrows() != 0 ? Y[0] : 0.0);
            out.rbegin()->setW(JGANDALF_NUMBER_OF_ITERATIONS, this->numberOfIterations);
            out.rbegin()->setW(JPP_COVERAGE_ORIENTATION,      input.coverage.orientation);
            out.rbegin()->setW(JPP_COVERAGE_POSITION,         input.coverage.position);
          }
        }
      }
 
      {
        for (JEvt::iterator i = out.begin(); i != out.end(); ++i) {
 
          double Q = numeric_limits<double>::max();
 
          for (JEvt::iterator p = out.begin(); p != out.end(); ++p) {
            if (p != i) {
              if (getDot(getDirection(*i), getDirection(*p)) <= cosLR) {
                if (i->getQ() - p->getQ() < Q) {
                  Q = i->getQ() - p->getQ();
                }
              }
            }
          }
 
          i->setW(JGANDALF_LIKELIHOOD_RATIO, Q);
        }
      }
 
      // 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

pmtParameters

const JPMTParametersMap& JRECONSTRUCTION::JMuonGandalf::pmtParameters

Definition at line 356 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.

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:

• JMuonGandalf.hh