JRECONSTRUCTION::JShowerDirectionPrefit Class Reference

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.
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 JShowerDirectionPrefitParameters_t &parameters) const
	Equality.
	ClassDef (JShowerDirectionPrefitParameters_t, 2)

Public Attributes
const JModuleRouter &	router
const JSummaryRouter &	summary
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	DMax_m
	maximal distance to optical module [m]
double	R_Hz
	default rate [Hz]
double	VMax_npe
	maximum number of of photo-electrons
double	scanAngle_deg
	scanning angle step in [deg]
double	Emin_GeV
	minimum energy to scan
double	Emax_GeV
	maximum energy to scan
int	En
	number of points to scan in energy range

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 70 of file JShowerDirectionPrefit.hh.

Member Typedef Documentation

JRegressor_t

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

private

Definition at line 75 of file JShowerDirectionPrefit.hh.

Constructor & Destructor Documentation

JShowerDirectionPrefit()

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 91 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

operator()()

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

reset()

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

equals()

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

ClassDef()

JRECONSTRUCTION::JShowerDirectionPrefitParameters_t::ClassDef ( JShowerDirectionPrefitParameters_t , 2  )

inherited

Member Data Documentation

Ev

vector<double> JRECONSTRUCTION::JShowerDirectionPrefit::Ev

private

Definition at line 77 of file JShowerDirectionPrefit.hh.

omega

JPP::JOmega3D JRECONSTRUCTION::JShowerDirectionPrefit::omega

private

Definition at line 78 of file JShowerDirectionPrefit.hh.

router

const JModuleRouter& JRECONSTRUCTION::JShowerDirectionPrefit::router

Definition at line 205 of file JShowerDirectionPrefit.hh.

summary

const JSummaryRouter& JRECONSTRUCTION::JShowerDirectionPrefit::summary

Definition at line 206 of file JShowerDirectionPrefit.hh.

numberOfPrefits

size_t JRECONSTRUCTION::JShowerDirectionPrefitParameters_t::numberOfPrefits

inherited

number of prefits

Definition at line 71 of file JShowerDirectionPrefitParameters_t.hh.

TMax_ns

double JRECONSTRUCTION::JShowerDirectionPrefitParameters_t::TMax_ns

inherited

maximum time for local coincidences [ns]

Definition at line 72 of file JShowerDirectionPrefitParameters_t.hh.

TMin_ns

double JRECONSTRUCTION::JShowerDirectionPrefitParameters_t::TMin_ns

inherited

minimum time for local coincidences [ns]

Definition at line 73 of file JShowerDirectionPrefitParameters_t.hh.

DMax_m

double JRECONSTRUCTION::JShowerDirectionPrefitParameters_t::DMax_m

inherited

maximal distance to optical module [m]

Definition at line 74 of file JShowerDirectionPrefitParameters_t.hh.

R_Hz

double JRECONSTRUCTION::JShowerDirectionPrefitParameters_t::R_Hz

inherited

default rate [Hz]

Definition at line 75 of file JShowerDirectionPrefitParameters_t.hh.

VMax_npe

double JRECONSTRUCTION::JShowerDirectionPrefitParameters_t::VMax_npe

inherited

maximum number of of photo-electrons

Definition at line 76 of file JShowerDirectionPrefitParameters_t.hh.

scanAngle_deg

double JRECONSTRUCTION::JShowerDirectionPrefitParameters_t::scanAngle_deg

inherited

scanning angle step in [deg]

Definition at line 77 of file JShowerDirectionPrefitParameters_t.hh.

Emin_GeV

double JRECONSTRUCTION::JShowerDirectionPrefitParameters_t::Emin_GeV

inherited

minimum energy to scan

Definition at line 78 of file JShowerDirectionPrefitParameters_t.hh.

Emax_GeV

double JRECONSTRUCTION::JShowerDirectionPrefitParameters_t::Emax_GeV

inherited

maximum energy to scan

Definition at line 79 of file JShowerDirectionPrefitParameters_t.hh.

En

int JRECONSTRUCTION::JShowerDirectionPrefitParameters_t::En

inherited

number of points to scan in energy range

Definition at line 80 of file JShowerDirectionPrefitParameters_t.hh.

---

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

• JShowerDirectionPrefit.hh