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 JShowerEnergyCorrection &correct, const int debug=0)
	Parameterized constructor. More...

JEvt	operator() (const KM3NETDAQ::JDAQEvent &event, const JFIT::JEvt &in)
	Declaration of the member function that actually performs the reconstruction. More...

double	getFinalBjY (double E_em, double E_h)

void	reset ()
	Reset fit parameters. More...

bool	equals (const JShowerBjorkenYParameters_t &parameters) const
	Equality. More...

	ClassDef (JShowerBjorkenYParameters_t, 1)

Public Attributes
const JModuleRouter &	router

const JSummaryRouter &	summary

const JShowerEnergyCorrection &	correct

size_t	numberOfPrefits
	number of prefits More...

double	TMax_ns
	maximum time for local coincidences [ns] More...

double	TMin_ns
	minimum time for local coincidences [ns] More...

double	roadWidth_m
	road width [m] More...

double	R_Hz
	default rate [Hz] More...

int	mestimator
	M-estimator. More...

double	fit_step
	simplex initial step More...

double	VMax_npe
	maximum number of of photo-electrons More...

int	NMax
	maximum number of iterations More...

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 74 of file JShowerBjorkenY.hh.

Member Typedef Documentation

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

private

Definition at line 79 of file JShowerBjorkenY.hh.

Constructor & Destructor Documentation

JRECONSTRUCTION::JShowerBjorkenY::JShowerBjorkenY	(	const JShowerBjorkenYParameters_t &	parameters,
		const JModuleRouter &	router,
		const JSummaryRouter &	summary,
		const std::string	pdfFile,
		const JShowerEnergyCorrection &	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 94 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

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 122 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(getPosition(*shower), getDirection(*shower), shower->getT(),
                       shower->getE(), 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());
 
           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 rate_Hz = summary.getRate(id);
               const size_t count   = top.count(id);
 
               buffer.push_back(JPMTW0(dom.getPMT(i), 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 sh_fit(this->value.getPosition(), this->value.getDirection(),
                           this->value.getT(), correct(this->value.getEem() + this->value.getEh()), this->value.getBy());
 
         double y = getFinalBjY(this->value.getEem(), this->value.getEh());
         
         sh_fit.rotate_back(R);
         
         sh_fit.add(sh.getPosition());
           
         out.push_back(getFit(JHistory(shower->getHistory()).add(JSHOWER_BJORKEN_Y), sh_fit, getQuality(chi2), 
                              NDF, sh_fit.getE()));
 
         out.rbegin()->setW(5, y);
         out.rbegin()->setW(6, this->value.getEem());
         out.rbegin()->setW(7, this->value.getEh());
 
       }
 
       return out;
     }

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

inline

Definition at line 229 of file JShowerBjorkenY.hh.

                                                {
       return E_h / (E_em + E_h);
     }

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;
     }

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);
     }