JFIT::JMuonGandalf Class Reference

class to handle Muon Gandalf reconstruction first, two other angular reconstructionS should be ran : JMuonPrefit and JMuonSimplex More...

#include <JMuonGandalf.hh>

Public Member Functions
	JMuonGandalf ()
	Default constructor. More...
	JMuonGandalf (const JLANG::JSharedPointer< const JDETECTOR::JModuleRouter > &router, const JFIT::JMuonGandalfParameters_t &parameters, std::string pdfFile)
	Parameterized constructor. More...
	~JMuonGandalf ()
void	getJEvt (const KM3NETDAQ::JDAQTimeslice &timeSlice, JFIT::JEvt &InPreFits, JFIT::JEvt &OutFits) const
	Declaration of Member function that actually performs the reconstruction. More...

Static Public Member Functions
static bool	compare (const JTRIGGER::JHitL0 &first, const JTRIGGER::JHitL0 &second)
	Compare hits by PMT identifier and time. More...

Private Types
typedef JFIT::JRegressor< JFIT::JLine3Z, JFIT::JGandalf >	JRegressor_t

Private Attributes
JLANG::JSharedPointer< const JDETECTOR::JModuleRouter >	router_
JFIT::JMuonGandalfParameters_t	parameters_
JRegressor_t	fitRegresor_

Detailed Description

class to handle Muon Gandalf reconstruction first, two other angular reconstructionS should be ran : JMuonPrefit and JMuonSimplex

Definition at line 47 of file JMuonGandalf.hh.

Member Typedef Documentation

JRegressor_t

typedef JFIT::JRegressor<JFIT::JLine3Z, JFIT::JGandalf> JFIT::JMuonGandalf::JRegressor_t

private

Definition at line 50 of file JMuonGandalf.hh.

Constructor & Destructor Documentation

JMuonGandalf() [1/2]

JFIT::JMuonGandalf::JMuonGandalf ( )

inline

Default constructor.

Definition at line 62 of file JMuonGandalf.hh.

    {}

JMuonGandalf() [2/2]

JFIT::JMuonGandalf::JMuonGandalf ( const JLANG::JSharedPointer< const JDETECTOR::JModuleRouter > &  router, const JFIT::JMuonGandalfParameters_t &  parameters, std::string  pdfFile  )

inline

Parameterized constructor.

Parameters

router	JSharedPointer of JModuleRouter, this is built via detector file.
parameters	struct that holds default-optimized parameters for the reconstruction, for ARCA and ORCA configuration. These default parameters could be found in $JPP_DATA.
pdfFile	PDF file descriptor

Definition at line 72 of file JMuonGandalf.hh.

                                   :
      router_(router), 
      parameters_(parameters)
    {
      JRegressor_t::debug              = 1;//see what to do with this
      JRegressor_t::T_ns.setRange(-50.0, +450.0);// ns
      JRegressor_t::Vmax_npe           =   10.0;   
      JRegressor_t::MAXIMUM_ITERATIONS = 10000;
      
      fitRegresor_ = {pdfFile, parameters.TTS_ns}; 
      
      fitRegresor_.parameters.resize(5);
      
      fitRegresor_.parameters[0] = JFIT::JLine3Z::pX();
      fitRegresor_.parameters[1] = JFIT::JLine3Z::pY();
      fitRegresor_.parameters[2] = JFIT::JLine3Z::pT();
      fitRegresor_.parameters[3] = JFIT::JLine3Z::pDX();
      fitRegresor_.parameters[4] = JFIT::JLine3Z::pDY();
    }

~JMuonGandalf()

JFIT::JMuonGandalf::~JMuonGandalf ( )

inline

Definition at line 94 of file JMuonGandalf.hh.

{}

Member Function Documentation

getJEvt()

void JFIT::JMuonGandalf::getJEvt	(	const KM3NETDAQ::JDAQTimeslice &	timeSlice,
		JFIT::JEvt &	InPreFits,
		JFIT::JEvt &	OutFits
	)		const

Declaration of Member function that actually performs the reconstruction.

Member function definition.

Parameters

timeSlice	which is contructed via a JDAQEvent
InPreFits	input JEvt which is a container of prefits, normally JMuonSimplex
OutFits	non const (output) JEvt.

Definition at line 133 of file JMuonGandalf.hh.

{
  if ( InPreFits.empty() ) return;
  
  const bool reprocess          = parameters_.reprocess;
  double roadWidth_m            = parameters_.roadWidth_m;
  const double R_Hz             = parameters_.R_Hz; 
  const size_t numberOfPrefits  = parameters_.numberOfPrefits; 
  const double E_GeV            = parameters_.E_GeV; 
 
  typedef std::vector<JTRIGGER::JHitL0>  JDataL0_t;
  typedef std::vector<JHitW0>  JDataW0_t;
 
  const JTRIGGER::JBuildL0<JTRIGGER::JHitL0>  buildL0;
  
  if (fitRegresor_.getRmax() < roadWidth_m) {
    roadWidth_m = fitRegresor_.getRmax();
  }
 
  JEvt::iterator __end = InPreFits.end();
 
  if (reprocess) {
    __end = std::partition( InPreFits.begin(), __end, 
                            JFIT::JHistory::is_not_event(JFIT::JFitApplication_t::JMUONGANDALF)
                            );
  }
 
  if (InPreFits.begin() != __end) {
 
    std::copy(InPreFits.begin(), __end, std::back_inserter(OutFits));
 
    if (numberOfPrefits > 0) {
      std::advance(__end = InPreFits.begin(), std::min(numberOfPrefits, OutFits.size()));
    }
 
    std::partial_sort(InPreFits.begin(), __end, InPreFits.end(), qualitySorter);
 
    __end = std::partition(InPreFits.begin(), __end, 
                           JFIT::JHistory::is_event(InPreFits.begin()->getHistory())
                           );
 
    JDataL0_t dataL0;
 
    buildL0(timeSlice, *router_, std::back_inserter(dataL0));
 
    if (dataL0.size() >= fitRegresor_.parameters.size()) {
 
      for (JEvt::const_iterator track = InPreFits.begin(); track != __end; ++track) {
 
        const JGEOMETRY3D::JRotation3D     R (JFIT::getDirection(*track));
        const JFIT::JLine1Z         tz(getPosition (*track).rotate(R), track->getT());
        const JFIT::JModel<JFIT::JLine1Z> match(tz, roadWidth_m, 
                                                JFIT::JRegressor<JFIT::JLine3Z, JFIT::JGandalf>::T_ns);
 
        //hit selection based on start value
 
        JDataW0_t data;
 
        for (JDataL0_t::const_iterator i = dataL0.begin(); i != dataL0.end(); ++i) {
 
          JFIT::JHitW0 hit(*i, R_Hz);
 
          hit.rotate(R);
 
          if (match(hit)) {
            data.push_back(hit);
          }
        }
 
        //select first hit
 
        std::sort(data.begin(), data.end(), compare);
 
        JDataW0_t::iterator __end = std::unique( data.begin(), data.end(), 
                                                 std::equal_to<KM3NETDAQ::JDAQPMTIdentifier>()
                                                 );
 
        const int NDF = std::distance(data.begin(), __end) - fitRegresor_.parameters.size();
 
        if (NDF >= 0) {
          
          //set fit parameters  
 
          if (track->getE() > 0.1)
            fitRegresor_.E_GeV   = track->getE();
          else
            fitRegresor_.E_GeV   = E_GeV;
 
          const double chi2 = fitRegresor_(JLine3Z(tz), data.begin(), __end);
 
          JGEOMETRY3D::JTrack3D tb(fitRegresor_.value);
 
          tb.rotate_back(R);
 
          const double energy(0);
          OutFits.push_back(JFIT::getFit(JFIT::JHistory(track->getHistory()).add(JFIT::JFitApplication_t::JMUONGANDALF), tb, 
                                     JFIT::getQuality(chi2), NDF,energy,JFIT::JFitStatus_t::OKAY));
 
          //set additional weights
 
          OutFits.rbegin()->setW(JGANDALF_BETA0_RAD,sqrt(fitRegresor_.error.getDX() * fitRegresor_.error.getDX() +
                                                     fitRegresor_.error.getDY() * fitRegresor_.error.getDY()));
          OutFits.rbegin()->setW(JGANDALF_BETA1_RAD,sqrt(fitRegresor_.error.getDX() * fitRegresor_.error.getDY()));
          OutFits.rbegin()->setW(JGANDALF_CHI2,                 chi2);
          OutFits.rbegin()->setW(JGANDALF_NUMBER_OF_HITS,       distance(data.begin(), __end));
          OutFits.rbegin()->setW(JGANDALF_LAMBDA,               fitRegresor_.lambda);
          OutFits.rbegin()->setW(JGANDALF_NUMBER_OF_ITERATIONS, fitRegresor_.numberOfIterations);
        }
      }
    }
    //apply default sorter
 
    std::sort(OutFits.begin(), OutFits.end(), qualitySorter);
  }
    
}

compare()

static bool JFIT::JMuonGandalf::compare ( const JTRIGGER::JHitL0 &  first, const JTRIGGER::JHitL0 &  second  )

inlinestatic

Compare hits by PMT identifier and time.

Parameters

first	first hit
second	second hit

Returns

true if first before second; else false

Definition at line 116 of file JMuonGandalf.hh.

    {
      if (std::equal_to<KM3NETDAQ::JDAQPMTIdentifier>()(first, second))
        return std::less<JTRIGGER::JHit>()(first, second);
      else
        return std::less<KM3NETDAQ::JDAQPMTIdentifier>()(first, second);
    }

Member Data Documentation

router_

JLANG::JSharedPointer<const JDETECTOR::JModuleRouter> JFIT::JMuonGandalf::router_

private

Definition at line 52 of file JMuonGandalf.hh.

parameters_

JFIT::JMuonGandalfParameters_t JFIT::JMuonGandalf::parameters_

private

Definition at line 53 of file JMuonGandalf.hh.

fitRegresor_

JRegressor_t JFIT::JMuonGandalf::fitRegresor_

mutableprivate

Definition at line 55 of file JMuonGandalf.hh.

---

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

• JMuonGandalf.hh