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
	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	DMax_m
	maximal distance to optical module [m] More...

double	R_Hz
	default rate [Hz] More...

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

double	scanAngle_deg
	scanning angle step in [deg] More...

double	Emin_GeV
	minimum energy to scan More...

double	Emax_GeV
	maximum energy to scan More...

int	En
	number of points to scan in energy range More...

Private Types
typedef JRegressor< JShower3EZ, JAbstractMinimiser >	JRegressor_t

Private Attributes
vector< double >	Ev

JPP::JOmega3D	omega

Detailed Description

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

Definition at line 69 of file JShowerDirectionPrefit.hh.

Member Typedef Documentation

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

private

Definition at line 74 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),
       omega(parameters.scanAngle_deg * JMATH::PI / 180.0),
       router(router),
       summary(summary)
     {
       using namespace JPP;
       
       JRegressor_t::debug  = debug;
       JRegressor_t::T_ns.setRange(parameters.TMin_ns, parameters.TMax_ns);
       JRegressor_t::Vmax_npe = VMax_npe;
           
       if (Emin_GeV > Emax_GeV || En <= 1) {
         THROW(JException, "Invalid energy input " << Emin_GeV << ' ' << Emax_GeV << ' ' << En);
       }
 
       const double base = std::pow((Emax_GeV / Emin_GeV), 1.0 / (En - 1));
       
       for (int i = 0; i != En; ++i) {
               Ev.push_back(Emin_GeV * std::pow(base, i));
       }
     }

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	DAQ event
in	input fits

Returns: result fits

Definition at line 125 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) {
 
         const JVertex3D Vertex(getPosition(*shower),shower->getT());
         multiset<JDAQPMTIdentifier> top;
         
         const JFIT::JModel<JPoint4D> match(Vertex, DMax_m, JRegressor_t::T_ns);
         
         for (JDataL0_t::const_iterator i = dataL0.begin(); i != dataL0.end(); ++i) {
           if (match(*i)) {
             top.insert(i->getPMTIdentifier());
           }
         }
 
         const JDetectorSubset_t subdetector(detector, Vertex.getPosition(), DMax_m);
         
         vector<JPMTW0> data;
         for (JDetectorSubset_t::const_iterator module = subdetector.begin(); module != subdetector.end(); ++module) {
           const JDAQSummaryFrame& frame = summary.getSummaryFrame(module->getID());
 
           JModule dom(*module);
 
           for (size_t 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);
 
               data.push_back(JPMTW0(dom.getPMT(i), rate_Hz, count));
             }
           }
         }
 
         const JShower3EZ sh(JVertex3D(JVector3D(0,0,0), sh.getT()), JVersor3Z(), 1);
 
         for (JOmega3D_t::const_iterator dir = omega.begin(); dir != omega.end(); ++dir) {
           const JRotation3D R(*dir);
           vector<double> chi2v(En, 0.0);
           for (vector<JPMTW0>::const_iterator i = data.begin(); i != data.end(); ++i) {  
             JPMTW0 pmt = *i;
             pmt.rotate(R);                              
             JNPE_t::result_type H1 = (*this).getH1(sh, pmt);
             JNPE_t::result_type H0 = (*this).getH0(pmt.getR());  // getH0 = Get background hypothesis value for time integrated PDF. 
             const bool hit = pmt.getN() != 0;
             for(size_t j=0;j!=chi2v.size();++j){
               chi2v[j]+= getChi2(get_value(H0+Ev[j]*H1), hit); // H1 is linear with E
             }
           }
           //Store only the lowest chi2 for a given direction
           auto minChi2 = std::min_element(chi2v.begin(), chi2v.end());
           JShower3E sh_fit(Vertex, JDirection3D(*dir), Ev[minChi2-chi2v.begin()]);
           out.push_back(getFit(JHistory(shower->getHistory()).add(JSHOWERDIRECTIONPREFIT), sh_fit, getQuality(*minChi2), data.size(), sh_fit.getE()));
         }
       }
         
       return out;
     
     }

void JRECONSTRUCTION::JShowerDirectionPrefitParameters_t::reset ( )

inlineinherited

Reset fit parameters.

Definition at line 35 of file JShowerDirectionPrefitParameters_t.hh.

     {
       numberOfPrefits             =    1;
       TMax_ns                     =   30;
       TMin_ns                     =  -30;
       DMax_m                      =   80;
       scanAngle_deg               =   10;
       R_Hz                        = 10.0e3;
       VMax_npe                    =   20.0;
       Emin_GeV                    = 1;
       Emax_GeV                    = 1000;
       En                          = 50;
     }

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 55 of file JShowerDirectionPrefitParameters_t.hh.

     {
       return (this->TMax_ns            == parameters.TMax_ns                  &&         
               this->TMin_ns                  == parameters.TMin_ns                  &&         
               this->numberOfPrefits          == parameters.numberOfPrefits          &&                 
               this->DMax_m                   == parameters.DMax_m                   &&
               this->R_Hz                     == parameters.R_Hz                     && 
               this->scanAngle_deg            == parameters.scanAngle_deg            &&
               this->VMax_npe                 == parameters.VMax_npe                 &&
         this->Emin_GeV                 == parameters.Emin_GeV                 &&
         this->Emax_GeV                 == parameters.Emax_GeV                 &&
         this->En                       == parameters.En                         ); 
     }