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

#include <JShowerDirectionPrefit.hh>

Inheritance diagram for JRECONSTRUCTION::JShowerDirectionPrefit:

Public Member Functions
	JShowerDirectionPrefit (const JShowerDirectionPrefitParameters_t &parameters, const JModuleRouter &router, const JSummaryRouter &summary, const std::string pdfFile, 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...

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

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

	ClassDef (JShowerDirectionPrefitParameters_t, 1)

Public Attributes
const JModuleRouter &	router

const JSummaryRouter &	summary

size_t	numberOfPrefits

size_t	numberOfOutfits

double	TMax_ns

double	TMin_ns

double	roadWidth_m

double	DMin_m

double	E_GeV

double	R_Hz

int	mestimator

double	MaxAngle_deg_lowE

double	scanAngle_deg_lowE

double	MaxAngle_deg_highE

double	scanAngle_deg_highE

Private Types
typedef JRegressor< JShower3EZ, JGandalf >	JRegressor_t

Detailed Description

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

Definition at line 71 of file JShowerDirectionPrefit.hh.

Member Typedef Documentation

typedef JRegressor<JShower3EZ, JGandalf> JRECONSTRUCTION::JShowerDirectionPrefit::JRegressor_t

private

Definition at line 76 of file JShowerDirectionPrefit.hh.

Constructor & Destructor Documentation

JRECONSTRUCTION::JShowerDirectionPrefit::JShowerDirectionPrefit	(	const JShowerDirectionPrefitParameters_t &	parameters,
		const JModuleRouter &	router,
		const JSummaryRouter &	summary,
		const std::string	pdfFile,
		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
debug	debug

Definition at line 90 of file JShowerDirectionPrefit.hh.

                                                                                  :
       JShowerDirectionPrefitParameters_t(parameters),
       JRegressor_t(pdfFile),
       router(router),
       summary(summary)
     {
       using namespace JPP;
       
       JRegressor_t::debug  = 1;
       JRegressor_t::T_ns.setRange(parameters.TMin_ns, parameters.TMax_ns);
       JRegressor_t::Vmax_npe = 20.0;
       JRegressor_t::MAXIMUM_ITERATIONS = 1000;
 
       this->parameters.resize(2);
       
       this->parameters[0] = JShower3EZ::pDX();
       this->parameters[1] = JShower3EZ::pDY();
       
       this->estimator.reset(getMEstimator(parameters.mestimator));
     }

Member Function Documentation

JEvt JRECONSTRUCTION::JShowerDirectionPrefit::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 121 of file JShowerDirectionPrefit.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) {
 
         JShower3E sh(getPosition(*shower), getDirection(*shower), shower->getT(), shower->getE());
         
         multiset<JDAQPMTIdentifier> top;
 
         double distance   = roadWidth_m;
         double max_angle  = MaxAngle_deg_lowE;
         double scan_angle = scanAngle_deg_lowE;
         
         if(sh.getE() > E_GeV) {
           distance   = DMin_m;  // to speed up the code 
           max_angle  = MaxAngle_deg_highE;
           scan_angle = scanAngle_deg_highE;
         }
 
         const JModel<JPoint4D> match(JPoint4D(sh.getPosition(), sh.getT()), distance, JRegressor_t::T_ns);
 
         JVector3D start_dir(0,0,0);
         
         for (JDataL0_t::const_iterator i = dataL0.begin(); i != dataL0.end(); ++i) {
 
           if (match(*i)) {
 
             top.insert(i->getPMTIdentifier());
 
             start_dir.add(JDirection3D(i->getPosition() - sh.getPosition()));     
 
           }
         }
 
         sh.setDirection(JDirection3D(start_dir));
 
         JDetectorSubset_t subdetector(detector, sh.getPosition(), distance);
 
         JOmega3D scan_directions(sh.getDirection(), 
                                  JOmega3D::range_type(0, max_angle * PI/180.0), 
                                  scan_angle * PI/180.0);
 
         for (JOmega3D_t::const_iterator dir = scan_directions.begin(); dir != scan_directions.end(); ++dir) {
           
           const JDirection3D conversion(*dir);      
           const JRotation3D R(conversion);
         
           vector<JPMTW0> buffer;
           
           for (JDetectorSubset_t::iterator module = subdetector.begin();
                module != subdetector.end(); ++module) {
             
             JModule dom(*module); 
             dom.rotate(R);  
 
             for (unsigned int i = 0; i != dom.size(); ++i) {
 
               const JDAQPMTIdentifier id(dom.getID(), i);
 
               const double rate_Hz = summary.getRate(id);
 
               buffer.push_back(JPMTW0(dom.getPMT(i), rate_Hz, top.count(id)));
             }      
           }
 
           double chi2 = (*this)(JShower3EZ(JVertex3D(JVector3D(0,0,0), sh.getT()), JVersor3Z(), 
                                            JEnergy(log10(sh.getE()))), buffer.begin(), buffer.end());
           
           double NDF = getCount(buffer.begin(), buffer.end()) - this->parameters.size();
 
           JShower3E sh_fit(this->value.getPosition(), this->value.getDirection(),
                            this->value.getT(), this->value.getE());
           
           sh_fit.rotate_back(R);
 
           sh_fit.add(sh.getPosition());
           
           out.push_back(getFit(JHistory(shower->getHistory()).add(JSHOWERDIRECTIONPREFIT), sh_fit, getQuality(chi2), 
                                NDF, sh_fit.getE())); 
         
         }
       }
 
       return out;
 
     }

void JRECONSTRUCTION::JShowerDirectionPrefitParameters_t::reset ( )

inlineinherited

Reset fit parameters.

Definition at line 35 of file JShowerDirectionPrefitParameters_t.hh.

     {
       numberOfPrefits             =   0;
       numberOfOutfits             =   4;
       TMax_ns                     =  30;
       TMin_ns                     = -30;
       roadWidth_m                 =  80;
       DMin_m                      =  80;
       E_GeV                       =  13;
       MaxAngle_deg_lowE           =  40;
       scanAngle_deg_lowE          =  10;
       MaxAngle_deg_highE          =  20;
       scanAngle_deg_highE         =   5;
       R_Hz                        =  10.0e3;
       mestimator                  =  JFIT::EM_NULL;
     }

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

inlineinherited

Equality.

Parameters

parameters fit parameters

Returns: true if equals; else false

Definition at line 58 of file JShowerDirectionPrefitParameters_t.hh.

     {
       return (this->TMax_ns                  == parameters.TMax_ns                  &&         
               this->TMin_ns                  == parameters.TMin_ns                  &&         
               this->numberOfPrefits          == parameters.numberOfPrefits          &&         
               this->numberOfOutfits          == parameters.numberOfOutfits          &&         
               this->roadWidth_m              == parameters.roadWidth_m              &&
               this->DMin_m                   == parameters.DMin_m                   &&
               this->E_GeV                    == parameters.E_GeV                    &&
               this->R_Hz                     == parameters.R_Hz                     && 
               this->mestimator               == parameters.mestimator               &&
               this->MaxAngle_deg_lowE        == parameters.MaxAngle_deg_lowE        &&
               this->scanAngle_deg_lowE       == parameters.scanAngle_deg_lowE       &&
               this->MaxAngle_deg_highE       == parameters.MaxAngle_deg_highE       &&
               this->scanAngle_deg_highE      == parameters.scanAngle_deg_highE); 
     }