JRECONSTRUCTION::JMuonStart Class Reference

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

#include <JMuonStart.hh>

Inheritance diagram for JRECONSTRUCTION::JMuonStart:

Public Types
typedef JRegressor< JEnergy >	JRegressor_t
typedef JTRIGGER::JHitL0	hit_type
typedef std::vector< hit_type >	buffer_type

Public Member Functions
	JMuonStart (const JMuonStartParameters_t &parameters, const JModuleRouter &router, const JSummaryRouter &summary, const std::string &pdfFile, const JK40Rates &rates_Hz, 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 JMuonStartParameters_t &parameters) const
	Equality.
	ClassDef (JMuonStartParameters_t, 2)

Public Attributes
double	roadWidth_m
	road width [m]
double	R_Hz
	default rate [Hz]
size_t	numberOfPrefits
	number of prefits
double	TMin_ns
	minimal time w.r.t. Cherenkov hypothesis [ns]
double	TMax_ns
	maximal time w.r.t. Cherenkov hypothesis [ns]
double	Pmin1
	minimal probability single observation
double	Pmin2
	minimal probability for twofold observations
int	Nmax2
	maximal number for twofold observations
bool	reprocess
	reprocess

Private Attributes
const JModuleRouter &	router
const JSummaryRouter &	summary
const JK40Rates	rates_Hz
int	debug

Detailed Description

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

Definition at line 112 of file JMuonStart.hh.

Member Typedef Documentation

JRegressor_t

JRegressor<JEnergy> JRECONSTRUCTION::JMuonStart::JRegressor_t

Definition at line 117 of file JMuonStart.hh.

hit_type

JTRIGGER::JHitL0 JRECONSTRUCTION::JMuonStart::hit_type

Definition at line 118 of file JMuonStart.hh.

buffer_type

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

Definition at line 119 of file JMuonStart.hh.

Constructor & Destructor Documentation

JMuonStart()

JRECONSTRUCTION::JMuonStart::JMuonStart ( const JMuonStartParameters_t & parameters, const JModuleRouter & router, const JSummaryRouter & summary, const std::string & pdfFile, const JK40Rates & rates_Hz, const int debug = 0 )

inline

Constructor.

Parameters

parameters	parameters
router	module router
summary	summary router
pdfFile	PDF file
rates_Hz	K40 rates [Hz]
debug	debug

Definition at line 133 of file JMuonStart.hh.

                                                       :
      JMuonStartParameters_t(parameters),
      JRegressor_t(pdfFile),
      router(router),
      summary(summary),
      rates_Hz(rates_Hz),
      debug(debug)
    {
      JRegressor_t::debug = debug;
      JRegressor_t::T_ns.setRange(parameters.TMin_ns, parameters.TMax_ns);
 
      if (this->getRmax() < parameters.roadWidth_m) {
        roadWidth_m = this->getRmax();
      }
    }

Member Function Documentation

operator()() [1/2]

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

inline

Fit function.

Parameters

event	event
in	start values

Returns

fit results

Definition at line 162 of file JMuonStart.hh.

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

operator()() [2/2]

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

inline

Fit function.

Parameters

data	dataL0
in	start values

Returns

fit results

Definition at line 184 of file JMuonStart.hh.

    {
      using namespace std;
      using namespace JPP;
 
      JEvt out;
 
      const JDetector& detector = router.getReference();
 
      const JStart start(Pmin1, Pmin2, Nmax2);  
      
      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());
        const JFIT::JModel<JLine1Z> match(tz, roadWidth_m, JRegressor_t::T_ns);
        
        map<int, multiset<int> > top;
        
        for (buffer_type::const_iterator i = data.begin(); i != data.end(); ++i) {
          
          hit_type hit(*i);
          
          hit.rotate(R);
          
          if (match(hit)) {
            top[hit.getModuleID()].insert(hit.getPMTAddress());
          }
        }
        
        struct JHit_t {
          
          double getZ() const { return z; }
          double getP() const { return p; }
          
          double z;
          double p;
        };
        
        vector<JHit_t> data;
        
        for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
          
          if (summary.hasSummaryFrame(module->getID())) {
 
            const JDAQSummaryFrame& frame = summary.getSummaryFrame(module->getID());
            
            JPosition3D pos(module->getPosition());
            
            pos.transform(R, tz.getPosition());
            
            if (pos.getX() <= roadWidth_m) {
              
              const double z = pos.getZ()  -  pos.getX() / getTanThetaC();
              
              const double p = getProbability(*module, frame, rates_Hz.getMultiplesRates(), JRegressor_t::T_ns.getLength(), top[module->getID()]);
              
              data.push_back({ z, p });
            }
          }
        }
        
        double Zmin       = 0.0;     // relative start position [m]
        double Zmax       = 0.0;     // relative end   position [m]
        double npe_total  = 0.0;     // npe due to mip
        double npe_missed = 0.0;     // npe due to mip not detected
        
        if (!data.empty()) {
          
          sort(data.begin(), data.end(), JLANG::make_comparator(&JHit_t::getZ));
          
          vector<JHit_t>::const_iterator         track_start = start.find(data. begin(), data. end());
          vector<JHit_t>::const_reverse_iterator track_end   = start.find(data.rbegin(), data.rend());
          
          if (track_start != data. end()) { Zmin = track_start->getZ(); }
          if (track_end   != data.rend()) { Zmax = track_end  ->getZ(); }
          
          // additional features
          
          for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
            
            JPosition3D pos(module->getPosition());
            
            pos.transform(R, tz.getPosition());
            
            if (pos.getX() <= roadWidth_m) {
              
              const double z = pos.getZ()  -  pos.getX() / getTanThetaC();
              
              if (z >= Zmin && z <= Zmax) {
                
                for (size_t i = 0; i != module->size(); ++i) {
                  
                  const JDAQSummaryFrame& frame = summary.getSummaryFrame(module->getID());
                  
                  if (getDAQStatus(frame, *module, i)      &&
                      getPMTStatus(frame, *module, i)      &&
                      frame[i].is_valid()                  &&
                      !module->getPMT(i).has(PMT_DISABLE)) {
                    
                    JPMT pmt = module->getPMT(i);
                    
                    pmt.transform(R, tz.getPosition());
                    
                    const double ya = this->getNPE(pmt, 0.0).getYA();
                    
                    npe_total += ya;
 
                    if (top[module->getID()].count(i) == 0) {
                      npe_missed += ya;
                    }
                  }
                }
              }
            }
          }
        }
 
        JFit fit = *track;
 
        // move track
          
        fit.move(Zmin, getSpeedOfLight());
          
        out.push_back(fit.add(JMUONSTART));
          
        out.rbegin()->setW(track->getW());
        out.rbegin()->setW(JSTART_NPE_MIP_TOTAL,  npe_total);
        out.rbegin()->setW(JSTART_NPE_MIP_MISSED, npe_missed);
        out.rbegin()->setW(JSTART_LENGTH_METRES,  Zmax - Zmin);
      }
      
      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();
      R_Hz            = 6.0e3;
      numberOfPrefits = 0;
      TMin_ns         = -25.0;     
      TMax_ns         = +75.0;     
      Pmin1           = 0.0;
      Pmin2           = 0.0;
      Nmax2           = 2;
      reprocess       = false;
    }

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->R_Hz             == parameters.R_Hz            &&
              this->numberOfPrefits  == parameters.numberOfPrefits &&
              this->TMin_ns          == parameters.TMin_ns         &&
              this->TMax_ns          == parameters.TMax_ns         &&
              this->Pmin1            == parameters.Pmin1           &&
              this->Pmin2            == parameters.Pmin2           &&
              this->Nmax2            == parameters.Nmax2           &&
              this->reprocess        == parameters.reprocess);
    }

ClassDef()

JRECONSTRUCTION::JMuonStartParameters_t::ClassDef ( JMuonStartParameters_t , 2  )

inherited

Member Data Documentation

router

const JModuleRouter& JRECONSTRUCTION::JMuonStart::router

private

Definition at line 320 of file JMuonStart.hh.

summary

const JSummaryRouter& JRECONSTRUCTION::JMuonStart::summary

private

Definition at line 321 of file JMuonStart.hh.

rates_Hz

const JK40Rates JRECONSTRUCTION::JMuonStart::rates_Hz

private

Definition at line 322 of file JMuonStart.hh.

debug

int JRECONSTRUCTION::JMuonStart::debug

private

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

R_Hz

double JRECONSTRUCTION::JMuonStartParameters_t::R_Hz

inherited

default rate [Hz]

Definition at line 76 of file JMuonStartParameters_t.hh.

numberOfPrefits

size_t JRECONSTRUCTION::JMuonStartParameters_t::numberOfPrefits

inherited

number of prefits

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.

Pmin1

double JRECONSTRUCTION::JMuonStartParameters_t::Pmin1

inherited

minimal probability single observation

Definition at line 80 of file JMuonStartParameters_t.hh.

Pmin2

double JRECONSTRUCTION::JMuonStartParameters_t::Pmin2

inherited

minimal probability for twofold observations

Definition at line 81 of file JMuonStartParameters_t.hh.

Nmax2

int JRECONSTRUCTION::JMuonStartParameters_t::Nmax2

inherited

maximal number for twofold observations

Definition at line 82 of file JMuonStartParameters_t.hh.

reprocess

bool JRECONSTRUCTION::JMuonStartParameters_t::reprocess

inherited

reprocess

Definition at line 83 of file JMuonStartParameters_t.hh.

---

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

• JMuonStart.hh