Fit the LDF and the universe. More...

#include <LDFFinderKG/LDFFinder.h>

Inheritance diagram for LDFFinderKG::LDFFinder:

Classes
struct	CurvFitResult

struct	LDFFitResult

struct	Silents

struct	Stage0

struct	Stage1

struct	Stage2

struct	Stage3

Public Types
enum	CurvatureParameter { eU = 0, eV = 1, eRCore = 2, eCTCore = 3 }

enum	InfoLevel { eNone = 0, eFinal = 1, eIntermediate = 2, eObscure = 3, eMinuit = 4 }

enum	LDFParameter { eShowerSize = 0, eCoreX = 1, eCoreY = 2, eShapeParameters = 3 }

enum	ResultFlag { eSuccess, eFailure, eBreakLoop, eContinueLoop }
	Flag returned by module methods to the RunController. More...

enum	VersionInfoType { eFilename = 1, eRevisionNumber = 2, eDate = 3, eTime = 4, eLastEditor = 5 }
	Different types of version info that can be retrieved from GetVersionInfo. More...

Public Member Functions
fwk::VModule::ResultFlag	Finish () override
	Finish: invoked at end of the run (NOT end of the event) More...

utl::Stopwatch &	GetStopwatch ()

const utl::Stopwatch &	GetStopwatch () const

std::string	GetVersionInfo (const VersionInfoType v) const
	Retrieve different sorts of module version info. More...

fwk::VModule::ResultFlag	Init () override
	Initialize: invoked at beginning of run (NOT beginning of event) More...

void	InitTiming ()

fwk::VModule::ResultFlag	Run (evt::Event &event) override
	Run: invoked once per event. More...

ResultFlag	RunWithTiming (evt::Event &event)

Static Public Member Functions
static std::string	GetResultFlagByName (const ResultFlag flag)

Private Member Functions
std::vector< int >	CleanSilentStation (const LDFFitResult &ldfres, const LDFFitConfig &config, const utl::Vector &showerAxis, std::vector< StationFitData > &stationFitData) const

double	FitCurvature (CurvFitResult &result, const utl::Point &core, const std::vector< StationFitData > &data) const

bool	FitGlobal (CurvFitResult &cres, LDFFitResult &lres, const LDFFitConfig &config, const std::vector< StationFitData > &data) const

bool	FitLDF (LDFFitResult &result, LDFFitConfig &config, const utl::Vector &showerAxis, const std::vector< StationFitData > &data) const

void	FitLDFSimplified (LDFFitResult &result, LDFFitConfig &config, const utl::Vector &showerAxis, const std::vector< StationFitData > &data) const
	Early estimate of shower size. More...

bool	FixBeta (const LDFFitResult &result, const utl::Vector &showerAxis, const std::vector< StationFitData > &data) const

bool	FixCore (utl::Point &core, utl::TimeStamp &coreTime, utl::Vector &showerAxis, const evt::Event &event) const

bool	FixGamma (const LDFFitResult &result, const utl::Vector &showerAxis, const std::vector< StationFitData > &data) const
	Criterion if gamma can be fitted analogue to beta must be checked/updated. More...

std::vector< StationFitData >	MakeStationFitData (const sevt::SEvent &sEvent) const
	Init station data used by LDF fit. More...

void	OutputResults (const evt::Event &event) const

double	ParameterizedRc (const double cosTheta, const double showerSize) const

	REGISTER_MODULE ("LDFFinderKG", LDFFinder)

void	SetRecData (evt::Event &event, const LDFFitResult &lresult, const CurvFitResult &cresult) const

Private Attributes
sdet::Array::Type	fArrayType = sdet::Array::e1500

std::vector< int >	fBadSilents

utl::Point	fBarycenter

utl::TimeStamp	fBaryTime

CoreType	fCoreType

EnergyConversion	fEnergyConversion

bool	fFixedAxis = false

bool	fFixedCurvature = false

bool	fGlobalFit = false

int	fInfoLevel = 0

LDFFitConfig	fLDFFitConfig

double	fMaxBaryToCoreDistance = 0

double	fMaxChi2 = 0

double	fMaxTheta = 0

int	fRoptN = 0

struct LDFFinderKG::LDFFinder::Silents	fSilentsVars

struct LDFFinderKG::LDFFinder::Stage0	fStage0

struct LDFFinderKG::LDFFinder::Stage1	fStage1

struct LDFFinderKG::LDFFinder::Stage2	fStage2

struct LDFFinderKG::LDFFinder::Stage3	fStage3

bool	fUseMaxLikelihood = false

bool	fUseSilentStationsFromOtherGrids = false

Detailed Description

Fit the LDF and the universe.

Fit is performed in several stages:

Stage 0: (no actual minimization)
Estimate LDF slope $\beta$ from the shower zenith angle $\vartheta$ .
Estimate $S_{1000}$ from the station closest to 1000m distance to core (or barycenter).
Stage 1: Fit $S_{1000}$ .
Fit $\beta$ if enough stations available.
Fit ...

Author: Markus Roth; Darko Veberic

Date: 06 Apr 2004; 25 Nov 2004 rewrite

Definition at line 71 of file SdReconstruction/LDFFinderKG/LDFFinder.h.

Member Enumeration Documentation

enum LDFFinderKG::LDFFinder::CurvatureParameter

Enumerator
eU
eV
eRCore
eCTCore

Definition at line 117 of file SdReconstruction/LDFFinderKG/LDFFinder.h.

enum LDFFinderKG::LDFFinder::InfoLevel

Enumerator
eNone
eFinal
eIntermediate
eObscure
eMinuit

Definition at line 77 of file SdReconstruction/LDFFinderKG/LDFFinder.h.

enum LDFFinderKG::LDFFinder::LDFParameter

Enumerator
eShowerSize
eCoreX
eCoreY
eShapeParameters

Definition at line 85 of file SdReconstruction/LDFFinderKG/LDFFinder.h.

enum fwk::VModule::ResultFlag

inherited

Flag returned by module methods to the RunController.

Enumerator
eSuccess	Report success to RunController.
eFailure	Report failure to RunController, causing RunController to terminate execution.
eBreakLoop	Break current loop. It works for nested loops too!
eContinueLoop	Skip remaining modules in the current loop and continue with next iteration of the loop.

Definition at line 60 of file VModule.h.

enum fwk::VModule::VersionInfoType

inherited

Different types of version info that can be retrieved from GetVersionInfo.

Enumerator
eFilename
eRevisionNumber
eDate
eTime
eLastEditor

Definition at line 110 of file VModule.h.

Member Function Documentation

vector< int > LDFFinder::CleanSilentStation	(	const LDFFitResult &	ldfres,
		const LDFFitConfig &	config,
		const utl::Vector &	showerAxis,
		std::vector< StationFitData > &	stationFitData
	)		const

private

Definition at line 1891 of file SdReconstruction/LDFFinderKG/LDFFinder.cc.

References LDFFinderKG::LDFFitConfig::fLDF, LDFFinderKG::LDFFinder::LDFFitResult::fPar, LDFFinderKG::gBaryCS, and LDFFinderKG::RPerp().

fwk::VModule::ResultFlag LDFFinderKG::LDFFinder::Finish ( )

inlineoverridevirtual

Finish: invoked at end of the run (NOT end of the event)

This method is for things that should be done at the end of the run (for example, closing files or writing out histograms) {You must override this method in your concrete module}

Implements fwk::VModule.

Definition at line 75 of file SdReconstruction/LDFFinderKG/LDFFinder.h.

References fwk::VModule::eSuccess.

double LDFFinder::FitCurvature	(	CurvFitResult &	result,
		const utl::Point &	core,
		const std::vector< StationFitData > &	data
	)		const

private

Definition at line 1755 of file SdReconstruction/LDFFinderKG/LDFFinder.cc.

References galactic::eV, LDFFinderKG::LDFFinder::CurvFitResult::fChi2, LDFFinderKG::LDFFinder::CurvFitResult::fCov, LDFFinderKG::LDFFinder::CurvFitResult::fNdof, LDFFinderKG::LDFFinder::CurvFitResult::fPar, LDFFinderKG::HasNaN(), INFO, jacobian(), utl::kilometer, utl::m, galactic::meter, OUT, and sqrt().

bool LDFFinder::FitGlobal	(	CurvFitResult &	cres,
		LDFFitResult &	lres,
		const LDFFitConfig &	config,
		const std::vector< StationFitData > &	data
	)		const

private

Definition at line 1928 of file SdReconstruction/LDFFinderKG/LDFFinder.cc.

bool LDFFinder::FitLDF	(	LDFFitResult &	result,
		LDFFitConfig &	config,
		const utl::Vector &	showerAxis,
		const std::vector< StationFitData > &	data
	)		const

private

Definition at line 1508 of file SdReconstruction/LDFFinderKG/LDFFinder.cc.

void LDFFinder::FitLDFSimplified	(	LDFFitResult &	result,
		LDFFitConfig &	config,
		const utl::Vector &	showerAxis,
		const std::vector< StationFitData > &	data
	)		const

private

bool LDFFinder::FixBeta	(	const LDFFitResult &	result,
		const utl::Vector &	showerAxis,
		const std::vector< StationFitData > &	data
	)		const

private

Definition at line 1000 of file SdReconstruction/LDFFinderKG/LDFFinder.cc.

References LDFFinderKG::LDFFinder::LDFFitResult::GetCore(), galactic::meter, OUT, LDFFinderKG::RPerp(), and utl::s.

bool LDFFinder::FixCore	(	utl::Point &	core,
		utl::TimeStamp &	coreTime,
		utl::Vector &	showerAxis,
		const evt::Event &	event
	)		const

private

bool LDFFinder::FixGamma	(	const LDFFitResult &	result,
		const utl::Vector &	showerAxis,
		const std::vector< StationFitData > &	data
	)		const

private

Criterion if gamma can be fitted analogue to beta must be checked/updated.

Definition at line 1062 of file SdReconstruction/LDFFinderKG/LDFFinder.cc.

References LDFFinderKG::LDFFinder::LDFFitResult::GetCore(), galactic::meter, OUT, LDFFinderKG::RPerp(), and utl::s.

std::string fwk::VModule::GetResultFlagByName ( const ResultFlag flag )

staticinherited

Definition at line 8 of file VModule.cc.

References fwk::VModule::eBreakLoop, fwk::VModule::eContinueLoop, fwk::VModule::eFailure, and fwk::VModule::eSuccess.

Referenced by fwk::RunController::DoNextInSequence(), fwk::RunController::FinishBranch(), and fwk::RunController::InitBranch().

utl::Stopwatch& fwk::VModule::GetStopwatch ( )

inlineinherited

Definition at line 106 of file VModule.h.

References fwk::VModule::fStopwatch.

const utl::Stopwatch& fwk::VModule::GetStopwatch ( ) const

inlineinherited

Definition at line 107 of file VModule.h.

References fwk::VModule::fStopwatch.

std::string fwk::VModule::GetVersionInfo ( const VersionInfoType v ) const

inherited

Retrieve different sorts of module version info.

Definition at line 26 of file VModule.cc.

Referenced by fwk::CentralConfig::GetConfig(), ThresholdCalculatorKG::ThresholdCalculator::Init(), fdDoubleBumpFinder::FdDoubleBumpFinder::Init(), LaserGeneratorNA::LaserGenerator::Init(), LaserLightSimulatorNA::LaserLightSimulator::Init(), FdElectronicsSimulatorOG::FdElectronicsSimulator::Init(), TelescopeSimulatorKG::TelescopeSimulator::Init(), TelescopeSimulatorKG2::TelescopeSimulator::Init(), SdSimpleSimKG::SdSimpleSim::Init(), ShowerPhotonGeneratorOG::ShowerPhotonGenerator::Init(), and testRunController::testModuleVersionInfo().

VModule::ResultFlag LDFFinder::Init ( )

overridevirtual

Initialize: invoked at beginning of run (NOT beginning of event)

This method is for things that should be done once at the beginning of a run (for example, booking histograms, performing calculations that need to be done only once, initializing parameters) {You must override this method in your concrete module}

Implements fwk::VModule.

Definition at line 222 of file SdReconstruction/LDFFinderKG/LDFFinder.cc.

References sdet::Array::e1500, sdet::Array::e433, sdet::Array::e750, LDFFinderKG::eFit, utl::wcd::eGAP2012_012, utl::wcd::eGAP2014_035, LDFFinderKG::eHybrid, LDFFinderKG::eMC, LDFFinderKG::eModeled, ERROR, io::eSuccess, LDFFinderKG::GetShowerSizeLabel(), INFO, and G4StationSimulatorOG::p.