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
double	roadWidth_m
	road width [m]
double	R_Hz
	default rate [Hz]
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

Private Attributes
const JK40Rates	rates_Hz

Detailed Description

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

Definition at line 73 of file JMuonStart.hh.

Member Typedef Documentation

JRegressor_t

JRegressor<JEnergy> JRECONSTRUCTION::JMuonStart::JRegressor_t

Definition at line 78 of file JMuonStart.hh.

module_type

JModuleL0 JRECONSTRUCTION::JMuonStart::module_type

Definition at line 79 of file JMuonStart.hh.

detector_type

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

Definition at line 80 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 127 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 153 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(), R_Hz), 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 197 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;
        double y1;
        double yx;
      };
        
      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);
 
            double y1 = 0.0;
            double yx = 0.0;
 
            for (size_t i = 0; i != module->size(); ++i) {
 
              if (module.getStatus(i)) {
 
                JPMT pmt = module->getPMT(i);
                    
                pmt.transform(R, tz.getPosition());
                    
                const double npe = this->getY1(pmt);
 
                y1 += npe;
 
                if (count_if(module.begin(), module.end(), make_predicate(&JHitR0::getPMT, (JDAQHit::JPMT_t) i)) == 0) {
                  yx += npe;
                }
              }
            }
 
            data.push_back({ z, p, y1, yx });
          }
        }
 
        double Zmin = 0.0;     // start position
        double Zmax = 0.0;     // end   position
        double y1   = 0.0;     // npe due to MIP between start and end position
        double yx   = 0.0;     // npe due to MIP between start and end position that are not detected
        double pb   = 0.0;     // background probability
        
        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());
 
            for (vector<JHit_t>::const_iterator i = q1; i != q2.base(); ++i) {
              y1 += i->y1;
              yx += i->yx;
              pb += log(i->p > Pmin ? i->p : Pmin);
            }
          }
        }
 
        JFit fit = *track;
          
        fit.push_back(event());
 
        // move track
 
        fit.move(Zmin, getSpeedOfLight());
          
        out.push_back(fit);
 
        // set additional values
 
        out.rbegin()->setW(track->getW());
        out.rbegin()->setW(JSTART_NPE_MIP_TOTAL,     y1);
        out.rbegin()->setW(JSTART_NPE_MIP_MISSED,    yx);
        out.rbegin()->setW(JSTART_LENGTH_METRES,     Zmax - Zmin);
        out.rbegin()->setW(JSTART_BACKGROUND_LOGP,   pb);
        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();
      R_Hz             = 6.0e3;
      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 61 of file JMuonStartParameters_t.hh.

    {
      return (this->roadWidth_m      == parameters.roadWidth_m      &&
              this->R_Hz             == parameters.R_Hz             &&
              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

private

Definition at line 346 of file JMuonStart.hh.

roadWidth_m

double JRECONSTRUCTION::JMuonStartParameters_t::roadWidth_m

inherited

road width [m]

Definition at line 77 of file JMuonStartParameters_t.hh.

R_Hz

double JRECONSTRUCTION::JMuonStartParameters_t::R_Hz

inherited

default rate [Hz]

Definition at line 78 of file JMuonStartParameters_t.hh.

numberOfPrefits

size_t JRECONSTRUCTION::JMuonStartParameters_t::numberOfPrefits

inherited

number of prefits

Definition at line 79 of file JMuonStartParameters_t.hh.

numberOfPostfits

size_t JRECONSTRUCTION::JMuonStartParameters_t::numberOfPostfits

inherited

number of postfits

Definition at line 80 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 81 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 82 of file JMuonStartParameters_t.hh.

Pmin

double JRECONSTRUCTION::JMuonStartParameters_t::Pmin

inherited

minimal probability

Definition at line 83 of file JMuonStartParameters_t.hh.

Pmin1

double JRECONSTRUCTION::JMuonStartParameters_t::Pmin1

inherited

minimal probability single observation

Definition at line 84 of file JMuonStartParameters_t.hh.

Pmin2

double JRECONSTRUCTION::JMuonStartParameters_t::Pmin2

inherited

minimal probability for twofold observations

Definition at line 85 of file JMuonStartParameters_t.hh.

Nmax2

int JRECONSTRUCTION::JMuonStartParameters_t::Nmax2

inherited

maximal number for twofold observations

Definition at line 86 of file JMuonStartParameters_t.hh.

---

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

• JMuonStart.hh