Fit the LDF and the universe. More...

#include <LDFFinderOG/LDFFinder.h>

Inheritance diagram for LDFFinderOG::LDFFinder:

Classes
struct	CurvatureFitInterface

struct	LDFFitInterface

struct	Stage0

struct	Stage1

struct	Stage2

struct	Stage3

Public Types
enum	CoreType { eMC, eHybrid, eFit }

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 ()
	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 ()
	Initialize: invoked at beginning of run (NOT beginning of event) More...

void	InitTiming ()

	REGISTER_MODULE ("LDFFinderOG", LDFFinder)

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

ResultFlag	RunWithTiming (evt::Event &event)

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

Static Public Attributes
static LDFFunctionType	fgLDFType

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

Private Member Functions
double	EstimateChi2 (const LDFFitInterface &fi) const
	Estimate Chi2 using Minuit. More...

bool	EstimateCurvature (CurvatureFitInterface &fi) const

void	EstimateEnergy (LDFFitInterface &fi) const

void	EstimateLDF (LDFFitInterface &fi, const utl::Point &core) const

void	EstimateLDFbySeedStations (LDFFitInterface &fi, const utl::Point &core) const

double	EstimateNStationsInFit (const LDFFitInterface &fi) const

void	EstimateS1000 (LDFFitInterface &fi) const
	First guess of S1000 to give a starting point. More...

double	FitCurvatureDriver (CurvatureFitInterface &fi) const

bool	FitLDF (LDFFitInterface &fi, const int nStations) const

double	FitLDFDriver (LDFFitInterface &fi) const
	Fit using Minuit. More...

bool	FixBeta (const LDFFitInterface &fi) const

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

bool	FixGamma (const LDFFitInterface &fi) const

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

void	OutputStations (const LDFFitInterface &fi) const

double	ParameterizedRc (const LDFFitInterface &fi) const

void	SetRecData (evt::Event &event, const LDFFitInterface &lfi, const CurvatureFitInterface &cfi) const

Static Private Member Functions
static void	CurvatureFitFnc (int &nPar, double const grad, double &value, double const par, const int iFlag)

static void	LDFFitChi2Fnc (int &nPar, double const grad, double &value, double const par, const int iFlag)

static void	LDFFitMaxLikeFnc (int &nPar, double const grad, double &value, double const par, const int iFlag)
	Likelihood function to minimize. More...

Private Attributes
int	fCoreType

EnergyConversion	fEnergyConversion

int	fInfoLevel

double	fMaxBaryToCoreDistance

double	fMaxChi2

double	fMaxTheta

bool	fMinuitOutput

utl::RandomEngine *	fRandomEngine

int	fRoptN

double	fSilentRadiusTransition

struct LDFFinderOG::LDFFinder::Stage0	fStage0

struct LDFFinderOG::LDFFinder::Stage1	fStage1

struct LDFFinderOG::LDFFinder::Stage2	fStage2

struct LDFFinderOG::LDFFinder::Stage3	fStage3

bool	fUseMaxLike

bool	fUseSilentStations

Static Private Attributes
static utl::Point	fgBarycenter

static utl::TimeStamp	fgBaryTime

static const sevt::SEvent *	fgCurrentSEvent

static bool	fgFixLDF

static utl::CoordinateSystemPtr	fgLocalCS

static bool	fgLowerLimit

static double	fgNKGFermiMu

static double	fgNKGFermiTau

static double	fgNoRecovery

static utl::Vector	fgShowerAxis

static double	fgSilentMaxRadius

static double	fgSilentRadiusTransition

static double	fgSilentSignalThreshold

static bool	fgUseSilentStations

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 62 of file ZombieGraveyard/LDFFinderOG/LDFFinder.h.

Member Enumeration Documentation

enum LDFFinderOG::LDFFinder::CoreType

Enumerator
eMC
eHybrid
eFit

Definition at line 208 of file ZombieGraveyard/LDFFinderOG/LDFFinder.h.

enum LDFFinderOG::LDFFinder::InfoLevel

private

Enumerator
eNone
eFinal
eIntermediate
eObscure
eMinuit

Definition at line 196 of file ZombieGraveyard/LDFFinderOG/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

void LDFFinder::CurvatureFitFnc	(	int &	nPar,
		double *const	grad,
		double &	value,
		double *const	par,
		const int	iFlag
	)

staticprivate

double LDFFinder::EstimateChi2 ( const LDFFitInterface & fi ) const

private

Estimate Chi2 using Minuit.

Definition at line 1387 of file ZombieGraveyard/LDFFinderOG/LDFFinder.cc.

References eObscure, LDFFinderOG::LDFFinder::LDFFitInterface::fBeta, LDFFinderOG::LDFFinder::LDFFitInterface::fCore, LDFFinderOG::LDFFinder::LDFFitInterface::fGamma, fgLocalCS, fgUseSilentStations, fMinuitOutput, LDFFinderOG::LDFFinder::LDFFitInterface::fS1000, utl::BasicVector< HepVector >::GetX(), utl::BasicVector< HepVector >::GetY(), LDFFitChi2Fnc(), and OUT.

Referenced by FitLDF().

bool LDFFinder::EstimateCurvature ( CurvatureFitInterface & fi ) const

private

Definition at line 1697 of file ZombieGraveyard/LDFFinderOG/LDFFinder.cc.

void LDFFinder::EstimateEnergy ( LDFFitInterface & fi ) const

private

Definition at line 1082 of file ZombieGraveyard/LDFFinderOG/LDFFinder.cc.

References LDFFinderOG::LDFFinder::LDFFitInterface::fEnergy, fEnergyConversion, LDFFinderOG::LDFFinder::LDFFitInterface::fEnergyErr, fgLocalCS, fgShowerAxis, LDFFinderOG::LDFFinder::LDFFitInterface::fS1000, LDFFinderOG::LDFFinder::LDFFitInterface::fS1000Err, utl::BasicVector< HepVector >::GetCosTheta(), and LDFFinderOG::EnergyConversion::GetEnergy().

Referenced by EstimateLDF(), and FitLDF().

void LDFFinder::EstimateLDF	(	LDFFitInterface &	fi,
		const utl::Point &	core
	)		const

inlineprivate

void LDFFinderOG::LDFFinder::EstimateLDFbySeedStations	(	LDFFitInterface &	fi,
		const utl::Point &	core
	)		const

private

double LDFFinder::EstimateNStationsInFit ( const LDFFitInterface & fi ) const

private

Definition at line 1056 of file ZombieGraveyard/LDFFinderOG/LDFFinder.cc.

References LDFFinderOG::LDFFinder::LDFFitInterface::fCore, utl::Fermi(), fgCurrentSEvent, fgShowerAxis, fgSilentMaxRadius, fgSilentRadiusTransition, fgUseSilentStations, sdet::Station::GetPosition(), sdet::SDetector::GetStation(), LDFFinderOG::RPerp(), sevt::SEvent::StationsBegin(), and sevt::SEvent::StationsEnd().

void LDFFinder::EstimateS1000 ( LDFFitInterface & fi ) const

private

First guess of S1000 to give a starting point.

Definition at line 1026 of file ZombieGraveyard/LDFFinderOG/LDFFinder.cc.

References utl::abs(), sevt::SEvent::CandidateStationsBegin(), sevt::SEvent::CandidateStationsEnd(), LDFFinderOG::LDFFinder::LDFFitInterface::fBeta, LDFFinderOG::LDFFinder::LDFFitInterface::fGamma, fgBarycenter, fgCurrentSEvent, fgLDFType, fgNKGFermiMu, fgNKGFermiTau, fgShowerAxis, LDFFinderOG::LDFFinder::LDFFitInterface::fS1000, sdet::Station::GetPosition(), sdet::SDetector::GetStation(), LDFFinderOG::LDFFunction(), galactic::meter, and LDFFinderOG::RPerp().

Referenced by EstimateLDF().

VModule::ResultFlag LDFFinder::Finish ( )

virtual

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 1282 of file ZombieGraveyard/LDFFinderOG/LDFFinder.cc.

References fwk::VModule::eSuccess.

double LDFFinder::FitCurvatureDriver ( CurvatureFitInterface & fi ) const

private

Definition at line 1837 of file ZombieGraveyard/LDFFinderOG/LDFFinder.cc.

Referenced by Run().

bool LDFFinder::FitLDF	(	LDFFitInterface &	fi,
		const int	nStations
	)		const

inlineprivate

Definition at line 381 of file ZombieGraveyard/LDFFinderOG/LDFFinder.cc.

Referenced by Run().

double LDFFinder::FitLDFDriver ( LDFFitInterface & fi ) const

private

Fit using Minuit.

Definition at line 1290 of file ZombieGraveyard/LDFFinderOG/LDFFinder.cc.

Referenced by FitLDF().

bool LDFFinder::FixBeta ( const LDFFitInterface & fi ) const

private

Definition at line 904 of file ZombieGraveyard/LDFFinderOG/LDFFinder.cc.

References sevt::SEvent::CandidateStationsBegin(), sevt::SEvent::CandidateStationsEnd(), eIntermediate, LDFFinderOG::LDFFinder::LDFFitInterface::fCore, fgCurrentSEvent, fgShowerAxis, sdet::Station::GetPosition(), sdet::SDetector::GetStation(), galactic::meter, OUT, and LDFFinderOG::RPerp().

Referenced by Run().

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

private

Definition at line 773 of file ZombieGraveyard/LDFFinderOG/LDFFinder.cc.

Referenced by Run().

bool LDFFinder::FixGamma ( const LDFFitInterface & fi ) const

private

Definition at line 965 of file ZombieGraveyard/LDFFinderOG/LDFFinder.cc.

References sevt::SEvent::CandidateStationsBegin(), sevt::SEvent::CandidateStationsEnd(), eIntermediate, LDFFinderOG::LDFFinder::LDFFitInterface::fCore, fgCurrentSEvent, fgShowerAxis, sdet::Station::GetPosition(), sdet::SDetector::GetStation(), galactic::meter, OUT, and LDFFinderOG::RPerp().

Referenced by Run().

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 ( )

virtual

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 235 of file ZombieGraveyard/LDFFinderOG/LDFFinder.cc.