JRECONSTRUCTION::JShowerBjorkenY Class Reference

class to handle the direction fit of the shower reconstruction, mainly dedicated for ORCA More...

#include <JShowerBjorkenY.hh>

Inheritance diagram for JRECONSTRUCTION::JShowerBjorkenY:

Public Member Functions
	JShowerBjorkenY (const JShowerBjorkenYParameters_t &parameters, const JModuleRouter &router, const JSummaryRouter &summary, const std::string pdfFile, const JEnergyCorrection &correct, const int debug=0)
	Parameterized constructor.
JEvt	operator() (const KM3NETDAQ::JDAQEvent &event, const JFIT::JEvt &in)
	Declaration of the member function that actually performs the reconstruction.
void	reset ()
	Reset fit parameters.
bool	equals (const JShowerBjorkenYParameters_t &parameters) const
	Equality.
	ClassDef (JShowerBjorkenYParameters_t, 1)

Static Public Member Functions
static double	getFinalBjY (double E_em, double E_h)

Public Attributes
const JModuleRouter &	router
const JSummaryRouter &	summary
const JEnergyCorrection &	correct
size_t	numberOfPrefits
	number of prefits
double	TMax_ns
	maximum time for local coincidences [ns]
double	TMin_ns
	minimum time for local coincidences [ns]
double	roadWidth_m
	road width [m]
double	R_Hz
	default rate [Hz]
int	mestimator
	M-estimator (see JFIT::JMEstimator_t)
double	fit_step
	simplex initial step
double	VMax_npe
	maximum number of of photo-electrons
int	NMax
	maximum number of iterations

Private Types
typedef JRegressor< JShowerEH, JSimplex >	JRegressor_t

Detailed Description

class to handle the direction fit of the shower reconstruction, mainly dedicated for ORCA

Definition at line 71 of file JShowerBjorkenY.hh.

Member Typedef Documentation

JRegressor_t

JRegressor<JShowerEH, JSimplex> JRECONSTRUCTION::JShowerBjorkenY::JRegressor_t

private

Definition at line 76 of file JShowerBjorkenY.hh.

Constructor & Destructor Documentation

JShowerBjorkenY()

JRECONSTRUCTION::JShowerBjorkenY::JShowerBjorkenY ( const JShowerBjorkenYParameters_t & parameters, const JModuleRouter & router, const JSummaryRouter & summary, const std::string pdfFile, const JEnergyCorrection & correct, const int debug = 0 )

inline

Parameterized constructor.

Parameters

parameters	struct that holds default-optimized parameters for the reconstruction, available in $JPP_DATA.
router	module router, this is built via detector file.
summary	summary router
pdfFile	PDF file
correct	energy correction
debug	debug

Definition at line 91 of file JShowerBjorkenY.hh.

                                                                   :
      JShowerBjorkenYParameters_t(parameters),
      JRegressor_t(pdfFile),
      router(router),
      summary(summary),
      correct(correct)
    {
      using namespace JPP;
      
      JRegressor_t::debug  = debug;
      JRegressor_t::T_ns.setRange(parameters.TMin_ns, parameters.TMax_ns);
      JRegressor_t::Vmax_npe = VMax_npe;
      JRegressor_t::MAXIMUM_ITERATIONS = NMax;
      
      this->estimator.reset(getMEstimator(parameters.mestimator));
    }

Member Function Documentation

operator()()

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

inline

Declaration of the member function that actually performs the reconstruction.

Parameters

event	= JDAQEvent
in	= input fits

Definition at line 119 of file JShowerBjorkenY.hh.

    {
      using namespace std;
      using namespace JPP;
 
      typedef vector<JHitL0>  JDataL0_t;
      JBuildL0<JHitL0>        buildL0;
 
      JEvt out;
 
      const JDetector &detector = router.getReference();
 
      JDataL0_t dataL0;
 
      buildL0(JDAQTimeslice(event, true), router, back_inserter(dataL0));
  
      for (JEvt::const_iterator shower = in.begin(); shower != in.end(); ++shower) {
        
        JShower3EY sh(getShower(*shower), 0.0);
 
        multiset<JDAQPMTIdentifier> top;
 
        const JModel<JPoint4D> match(JPoint4D(sh.getPosition(), sh.getT()), roadWidth_m, JRegressor_t::T_ns);
 
        for (JDataL0_t::const_iterator i = dataL0.begin(); i != dataL0.end(); ++i) {
          
          if (match(*i)) {
            top.insert(i->getPMTIdentifier());
          }
        }
        
        JDetectorSubset_t subdetector(detector, sh.getPosition(), roadWidth_m);
 
        const JDirection3D conversion(sh.getDirection());           
        const JRotation3D R(conversion);
        
        vector<JPMTW0> buffer;
          
        for (JDetectorSubset_t::iterator module = subdetector.begin();
             module != subdetector.end(); ++module) {
            
          const JDAQSummaryFrame& frame = summary.getSummaryFrame(module->getID(), R_Hz);
 
          JModule dom(*module); 
 
          dom.rotate(R);  
 
          for (unsigned int i = 0; i != dom.size(); ++i) {
 
            if (getDAQStatus(frame, *module, i)            &&
                getPMTStatus(frame, *module, i)            &&
                frame[i].is_valid()                        &&
                !module->getPMT(i).has(PMT_DISABLE)) {
 
              const JDAQPMTIdentifier id(module->getID(), i);
                
              const double QE      = 1.0;
              const double rate_Hz = summary.getRate(id, this->R_Hz);
              const size_t count   = top.count(id);
 
              buffer.push_back(JPMTW0(dom.getPMT(i), QE, rate_Hz, count));
            }
          }   
        }
 
        this->step.resize(2);
        this->step[0] = JShowerEH(JPoint4D(JVector3D(), 0.0), JVersor3Z(), fit_step, 0.0, 0.0);
        this->step[1] = JShowerEH(JPoint4D(JVector3D(), 0.0), JVersor3Z(), 0.0, fit_step, 0.0);
 
        double f_h = 1 - 0.681 * (std::pow(shower->getE()/0.863, -0.207));
        
        double chi2 = (*this)(JShowerEH(JVertex3D(JVector3D(0,0,0), sh.getT()),
                                        JVersor3Z(), 
                                        log10(sh.getE()),
                                        log10(f_h*sh.getE()),
                                        sh.getBjY()),
                              buffer.begin(), buffer.end());
 
        double NDF = getCount(buffer.begin(), buffer.end()) - this->step.size(); 
        
        JShower3EY result(JShower3E(JShower3D(this->value.getVertex(),
                                              this->value.getDirection()),
                                    correct(this->value.getEem() + this->value.getEh())),
                          this->value.getBy());
 
        double y = getFinalBjY(this->value.getEem(), this->value.getEh());
        
        result.rotate_back(R);
        
        result.add(sh.getPosition());
          
        out.push_back(getFit(JHistory(shower->getHistory()).add(JSHOWER_BJORKEN_Y), result, getQuality(chi2), NDF, result.getE()));
 
        out.rbegin()->setW(5, y);
        out.rbegin()->setW(6, this->value.getEem());
        out.rbegin()->setW(7, this->value.getEh());
 
      }
 
      return out;
    }

getFinalBjY()

static double JRECONSTRUCTION::JShowerBjorkenY::getFinalBjY ( double E_em, double E_h )

inlinestatic

Definition at line 231 of file JShowerBjorkenY.hh.

                                                             {
      return E_h / (E_em + E_h);
    }

reset()

void JRECONSTRUCTION::JShowerBjorkenYParameters_t::reset ( )

inlineinherited

Reset fit parameters.

Definition at line 35 of file JShowerBjorkenYParameters_t.hh.

    {
      numberOfPrefits             =    1;
      TMax_ns                     =   30;
      TMin_ns                     =  -30;
      roadWidth_m                 =   80;
      R_Hz                        = 10.0e3;
      mestimator                  = JFIT::EM_NULL;
      fit_step                    =    0.5;
      VMax_npe                    =   20.0;
      NMax                        = 1000;
    }

equals()

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

inlineinherited

Equality.

Parameters

parameters

fit parameters

Returns

true if equals; else false

Definition at line 54 of file JShowerBjorkenYParameters_t.hh.

    {
      return (this->TMax_ns                  == parameters.TMax_ns                  &&         
              this->TMin_ns                  == parameters.TMin_ns                  &&         
              this->numberOfPrefits          == parameters.numberOfPrefits          &&         
              this->roadWidth_m              == parameters.roadWidth_m              &&
              this->R_Hz                     == parameters.R_Hz                     && 
              this->mestimator               == parameters.mestimator               &&
              this->fit_step                 == parameters.fit_step                 &&
              this->VMax_npe                 == parameters.VMax_npe                 &&
              this->NMax                     == parameters.NMax);
    }

ClassDef()

JRECONSTRUCTION::JShowerBjorkenYParameters_t::ClassDef ( JShowerBjorkenYParameters_t , 1  )

inherited

Member Data Documentation

router

const JModuleRouter& JRECONSTRUCTION::JShowerBjorkenY::router

Definition at line 221 of file JShowerBjorkenY.hh.

summary

const JSummaryRouter& JRECONSTRUCTION::JShowerBjorkenY::summary

Definition at line 222 of file JShowerBjorkenY.hh.

correct

const JEnergyCorrection& JRECONSTRUCTION::JShowerBjorkenY::correct

Definition at line 223 of file JShowerBjorkenY.hh.

numberOfPrefits

size_t JRECONSTRUCTION::JShowerBjorkenYParameters_t::numberOfPrefits

inherited

number of prefits

Definition at line 69 of file JShowerBjorkenYParameters_t.hh.

TMax_ns

double JRECONSTRUCTION::JShowerBjorkenYParameters_t::TMax_ns

inherited

maximum time for local coincidences [ns]

Definition at line 70 of file JShowerBjorkenYParameters_t.hh.

TMin_ns

double JRECONSTRUCTION::JShowerBjorkenYParameters_t::TMin_ns

inherited

minimum time for local coincidences [ns]

Definition at line 71 of file JShowerBjorkenYParameters_t.hh.

roadWidth_m

double JRECONSTRUCTION::JShowerBjorkenYParameters_t::roadWidth_m

inherited

road width [m]

Definition at line 72 of file JShowerBjorkenYParameters_t.hh.

R_Hz

double JRECONSTRUCTION::JShowerBjorkenYParameters_t::R_Hz

inherited

default rate [Hz]

Definition at line 73 of file JShowerBjorkenYParameters_t.hh.

mestimator

int JRECONSTRUCTION::JShowerBjorkenYParameters_t::mestimator

inherited

M-estimator (see JFIT::JMEstimator_t)

Definition at line 74 of file JShowerBjorkenYParameters_t.hh.

fit_step

double JRECONSTRUCTION::JShowerBjorkenYParameters_t::fit_step

inherited

simplex initial step

Definition at line 75 of file JShowerBjorkenYParameters_t.hh.

VMax_npe

double JRECONSTRUCTION::JShowerBjorkenYParameters_t::VMax_npe

inherited

maximum number of of photo-electrons

Definition at line 76 of file JShowerBjorkenYParameters_t.hh.

NMax

int JRECONSTRUCTION::JShowerBjorkenYParameters_t::NMax

inherited

maximum number of iterations

Definition at line 77 of file JShowerBjorkenYParameters_t.hh.

---

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

• JShowerBjorkenY.hh