JRECONSTRUCTION::JMuonStart Class Reference

Auxiliary class to determine start and end position of muon trajectory. More...

#include <JMuonStart.hh>

Inheritance diagram for JRECONSTRUCTION::JMuonStart:

Classes
struct	input_type
	Input data type. More...

Public Types
typedef JRegressor< JEnergy >	JRegressor_t
typedef JModuleL0	module_type
typedef std::vector< module_type >	detector_type

Public Member Functions
	JMuonStart (const JMuonStartParameters_t &parameters, const storage_type &storage, const JK40Rates &rates_Hz, 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 JMuonStartParameters_t &parameters) const
	Equality.
	ClassDef (JMuonStartParameters_t, 3)

Public Attributes
const JK40Rates	rates_Hz
double	roadWidth_m
	road width [m]
size_t	numberOfPrefits
	number of prefits
size_t	numberOfPostfits
	number of postfits
double	TMin_ns
	minimal time w.r.t. Cherenkov hypothesis [ns]
double	TMax_ns
	maximal time w.r.t. Cherenkov hypothesis [ns]
double	Pmin
	minimal probability
double	Pmin1
	minimal probability single observation
double	Pmin2
	minimal probability for twofold observations
int	Nmax2
	maximal number for twofold observations

Detailed Description

Auxiliary class to determine start and end position of muon trajectory.

Definition at line 71 of file JMuonStart.hh.

Member Typedef Documentation

JRegressor_t

JRegressor<JEnergy> JRECONSTRUCTION::JMuonStart::JRegressor_t

Definition at line 76 of file JMuonStart.hh.

module_type

JModuleL0 JRECONSTRUCTION::JMuonStart::module_type

Definition at line 77 of file JMuonStart.hh.

detector_type

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

Definition at line 78 of file JMuonStart.hh.

Constructor & Destructor Documentation

JMuonStart()

JRECONSTRUCTION::JMuonStart::JMuonStart ( const JMuonStartParameters_t & parameters, const storage_type & storage, const JK40Rates & rates_Hz, const int debug = 0 )

inline

Constructor.

Parameters

parameters	parameters
storage	storage
rates_Hz	K40 rates [Hz]
debug	debug

Definition at line 125 of file JMuonStart.hh.

                                                       :
      JMuonStartParameters_t(parameters),
      JRegressor_t(storage),
      rates_Hz(rates_Hz)
    {
      if (this->getRmax() < parameters.roadWidth_m) {
        roadWidth_m = this->getRmax();
      }
 
      JRegressor_t::debug = debug;
    }

Member Function Documentation

getInput()

input_type JRECONSTRUCTION::JMuonStart::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 JMuonStart.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(), rates_Hz.getSinglesRate()), data[module.getID()]));
        }
      }
 
      return input;
    }

operator()()

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

inline

Fit function.

Parameters

input

input data

Returns

fit results

Definition at line 195 of file JMuonStart.hh.

    {
      using namespace std;
      using namespace JPP;
 
      JEvent event(JMUONSTART);
 
      JEvt out;
 
      // select start values
 
      JEvt in  = input.in;
 
      sort(in.begin(), in.end(), qualitySorter);
 
      if (numberOfPrefits != 0 && numberOfPrefits < in.size()) {
 
        JEvt::iterator __end = next(in.begin(), numberOfPrefits);
 
        if (numberOfPostfits > 0) {
 
          __end = gridify(__end, in.end(), numberOfPostfits);
        }
 
        in.erase(__end, in.end());
      }
 
      if (!in.empty()) {
        in.select(JHistory::is_event(in.begin()->getHistory()));
      }
 
      const JStart start(Pmin1, Pmin2, Nmax2);  
      
      struct JHit_t {
          
        double getZ()  const { return z; }
        double getP()  const { return p; }
          
        double z;
        double p;
      };
        
      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());
        
        vector<JHit_t> data;
 
        for (const auto& module : input.data) {
 
          JPosition3D pos(module->getPosition());
            
          pos.transform(R, tz.getPosition());
            
          if (pos.getX() <= roadWidth_m) {
 
            const double z = pos.getZ()  -  pos.getX() / getTanThetaC();
            const double t = tz .getT()  + (pos.getZ() + pos.getX() * getKappaC()) * getInverseSpeedOfLight();
 
            const double p = module.getProbability(rates_Hz.getMultiplesRates(), T_ns + t);
 
            data.push_back({ z, p });
          }
        }
 
        double Zmin = 0.0;     // start position
        double Zmax = 0.0;     // end   position
        
        if (!data.empty()) {
          
          sort(data.begin(), data.end(), make_comparator(&JHit_t::getZ));
 
          vector<JHit_t>::const_iterator         q1 = start.find(data. begin(), data. end());
          vector<JHit_t>::const_reverse_iterator q2 = start.find(data.rbegin(), data.rend());
 
          if (q1 != data.end() && q2 != data.rend()) {
 
            vector<JHit_t>::const_iterator         p1 = q1;  if (p1 != data. begin()) { --p1; }
            vector<JHit_t>::const_reverse_iterator p2 = q2;  if (p2 != data.rbegin()) { --p2; }
 
            Zmin = 0.5 * (p1->getZ() + q1->getZ());
            Zmax = 0.5 * (p2->getZ() + q2->getZ());
          }
        }
 
        JFit fit = *track;
          
        fit.push_back(event());
 
        out.push_back(fit);
 
        // set additional values
 
        out.rbegin()->setW(JSTART_ZMIN_M,            Zmin);
        out.rbegin()->setW(JSTART_ZMAX_M,            Zmax);
        out.rbegin()->setW(JSTART_LENGTH_METRES,     Zmax - Zmin);
        out.rbegin()->setW(JPP_COVERAGE_ORIENTATION, input.coverage.orientation);
        out.rbegin()->setW(JPP_COVERAGE_POSITION,    input.coverage.position);
      }
 
      // apply default sorter
 
      sort(out.begin(), out.end(), qualitySorter);
 
      copy(input.in.begin(), input.in.end(), back_inserter(out));
      
      return out;
    }

reset()

void JRECONSTRUCTION::JMuonStartParameters_t::reset ( )

inlineinherited

Reset fit parameters.

Definition at line 41 of file JMuonStartParameters_t.hh.

    {
      roadWidth_m      = std::numeric_limits<double>::max();
      numberOfPrefits  = 0;
      numberOfPostfits = 0;
      TMin_ns          = -25.0;
      TMax_ns          = +75.0;
      Pmin             = 0.0;
      Pmin1            = 0.0;
      Pmin2            = 0.0;
      Nmax2            = 2;
    }

equals()

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

inlineinherited

Equality.

Parameters

parameters

fit parameters

Returns

true if equals; else false

Definition at line 60 of file JMuonStartParameters_t.hh.

    {
      return (this->roadWidth_m      == parameters.roadWidth_m      &&
              this->numberOfPrefits  == parameters.numberOfPrefits  &&
              this->numberOfPostfits == parameters.numberOfPostfits &&
              this->TMin_ns          == parameters.TMin_ns          &&
              this->TMax_ns          == parameters.TMax_ns          &&
              this->Pmin             == parameters.Pmin             &&
              this->Pmin1            == parameters.Pmin1            &&
              this->Pmin2            == parameters.Pmin2            &&
              this->Nmax2            == parameters.Nmax2);
    }

ClassDef()

JRECONSTRUCTION::JMuonStartParameters_t::ClassDef ( JMuonStartParameters_t , 3  )

inherited

Member Data Documentation

rates_Hz

const JK40Rates JRECONSTRUCTION::JMuonStart::rates_Hz

Definition at line 306 of file JMuonStart.hh.

roadWidth_m

double JRECONSTRUCTION::JMuonStartParameters_t::roadWidth_m

inherited

road width [m]

Definition at line 75 of file JMuonStartParameters_t.hh.

numberOfPrefits

size_t JRECONSTRUCTION::JMuonStartParameters_t::numberOfPrefits

inherited

number of prefits

Definition at line 76 of file JMuonStartParameters_t.hh.

numberOfPostfits

size_t JRECONSTRUCTION::JMuonStartParameters_t::numberOfPostfits

inherited

number of postfits

Definition at line 77 of file JMuonStartParameters_t.hh.

TMin_ns

double JRECONSTRUCTION::JMuonStartParameters_t::TMin_ns

inherited

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

Definition at line 78 of file JMuonStartParameters_t.hh.

TMax_ns

double JRECONSTRUCTION::JMuonStartParameters_t::TMax_ns

inherited

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

Definition at line 79 of file JMuonStartParameters_t.hh.

Pmin

double JRECONSTRUCTION::JMuonStartParameters_t::Pmin

inherited

minimal probability

Definition at line 80 of file JMuonStartParameters_t.hh.

Pmin1

double JRECONSTRUCTION::JMuonStartParameters_t::Pmin1

inherited

minimal probability single observation

Definition at line 81 of file JMuonStartParameters_t.hh.

Pmin2

double JRECONSTRUCTION::JMuonStartParameters_t::Pmin2

inherited

minimal probability for twofold observations

Definition at line 82 of file JMuonStartParameters_t.hh.

Nmax2

int JRECONSTRUCTION::JMuonStartParameters_t::Nmax2

inherited

maximal number for twofold observations

Definition at line 83 of file JMuonStartParameters_t.hh.

---

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

• JMuonStart.hh