determines energy deposit profile and fits a four-parameter Gaisser-Hillas to it (telescope-based). More...

#include <FdEnergyDepositFinder.h>

Inheritance diagram for FdEnergyDepositFinderKG::FdEnergyDepositFinder:

Public Types
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)

Protected Types
enum	InfoLevel { eInfoNone = 0, eInfoFinal = 1, eInfoIntermediate = 2, eInfoDebug = 3 }

Protected Attributes
int	fInfoLevel = 0

Private Types
enum	EErrorPropStage { eDefaultReco, eRp, eChi0, eT0, eSDPTheta, eSDPPhi, eMie, eNErrorPropagation, eFirstGeometryProp = eRp, eFirstErrProp = eRp, eLastGeometryProp = eSDPPhi }

enum	EErrorPropType { eMinus = -1, ePlus = 1 }

typedef std::map < EErrorPropStage, evt::GaisserHillas4Parameter >	RecResults

Private Member Functions
void	DumpCurrentParameters (const evt::ShowerFRecData &, bool, const EFitFunctionType) const

bool	FitProfile (fevt::Eye &eye, const EFitFunctionType type)

bool	GuessGHParameters (evt::ShowerFRecData &shower) const

bool	IsEventCandidate (const fevt::Eye &eye) const

bool	PrepareEyeCopy (fevt::Eye &eye, const EErrorPropType type) const

void	PropagateUncertainties (fevt::Eye &eye) const

	REGISTER_MODULE ("FdEnergyDepositFinderKG", FdEnergyDepositFinder)

Private Attributes
int	fAafCorrection = 0

CFMatrixCalculator	fCFMatrixCalculator

CFMatrixCalculator	fCFMatrixCalculatorDense

utl::InvisibleEnergy::ECompositionModel	fComposition = utl::InvisibleEnergy::eData

unsigned int	fDenseMatrixDim = 0

double	fElectronEnergyThreshold = 1*utl::MeV

EErrorPropStage	fErrPropStage

int	fErrPropVerbosity = 0

GHShapeParameters	fGHShapePars

utl::InvisibleEnergy::EInteractionModel	fInvisibleEnergyModel = utl::InvisibleEnergy::eDATA

int	fMaxErrPropPoints = 0

double	fMaxZenith = 0

bool	fOnlyNonMono = false

ProfileCalculator	fProfileCalculator

ProfileFitter	fProfileFitter

RecResults	fRecResultsMinus

RecResults	fRecResultsPlus

bool	fRescaleErrors = false

bool	fUseNoiseBins = false

Detailed Description

determines energy deposit profile and fits a four-parameter Gaisser-Hillas to it (telescope-based).

Module to reconstruct the energy deposit profiles based on the light at aperture of telescopes. Further details can be found in

Cherenkov-Fluorescence matrix approach: M. Unger et al., NIMA 588 (2008) 433
lateral distribution of fluorescence light: D. Gora et al., Astropart. Phys. 24 (2006) 484
lateral distribution of Cherenkov light: B.R. Dawson, M. Giller, and G. Wieczorek, Proc. 30th ICRC (2007)
multiple scattering correction: M.D. Roberts, J. Phys. G 31 (2005) 1291
invisible energy correction: H.M.J. Barbosa et al., Astropart. Phys. 22 (2004) 159
implementation of lateral distributions: M. Unger et al., GAP-2008-052
log-likelihood and telescope-based fit: S. Falk et al., GAP-2010-123

Author: Michael Unger

Date: May 2010 first version; Feb 2011 numerical error propagation

Definition at line 56 of file FdEnergyDepositFinder.h.

Member Typedef Documentation

typedef std::map<EErrorPropStage, evt::GaisserHillas4Parameter> FdEnergyDepositFinderKG::FdEnergyDepositFinder::RecResults

private

Definition at line 101 of file FdEnergyDepositFinder.h.

Member Enumeration Documentation

enum FdEnergyDepositFinderKG::FdEnergyDepositFinder::EErrorPropStage

private

Enumerator
eDefaultReco
eRp
eChi0
eT0
eSDPTheta
eSDPPhi
eMie
eNErrorPropagation
eFirstGeometryProp
eFirstErrProp
eLastGeometryProp

Definition at line 64 of file FdEnergyDepositFinder.h.

enum FdEnergyDepositFinderKG::FdEnergyDepositFinder::EErrorPropType

private

Enumerator
eMinus
ePlus

Definition at line 78 of file FdEnergyDepositFinder.h.

enum fwk::VModule::InfoLevel

protectedinherited

Enumerator
eInfoNone
eInfoFinal
eInfoIntermediate
eInfoDebug

Definition at line 125 of file VModule.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 FdEnergyDepositFinderKG::FdEnergyDepositFinder::DumpCurrentParameters	(	const evt::ShowerFRecData &	shower,
		bool	cfmUpdate,
		const EFitFunctionType	funcType
	)		const

private

Definition at line 802 of file FdEnergyDepositFinder.cc.

References utl::cm2, FdEnergyDepositFinderKG::eGH2dEdXChi2, FdEnergyDepositFinderKG::eGH4dEdXChi2, FdEnergyDepositFinderKG::eGH4LightLogLike, FdEnergyDepositFinderKG::eOne, FdEnergyDepositFinderKG::eTwo, fwk::CoordinateSystemRegistry::eUndefined, utl::g, evt::VGaisserHillasParameter::GetChiSquare(), evt::ShowerFRecData::GetGHParameters(), evt::VGaisserHillasParameter::GetNdof(), evt::VGaisserHillasParameter::GetNMax(), evt::GaisserHillas4Parameter::GetShapeParameter(), evt::VGaisserHillasParameter::GetXMax(), FdEnergyDepositFinderKG::gHRule, evt::ShowerFRecData::HasGHParameters(), and PeV.

fwk::VModule::ResultFlag FdEnergyDepositFinderKG::FdEnergyDepositFinder::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 61 of file FdEnergyDepositFinder.h.

References fwk::VModule::eSuccess.

bool FdEnergyDepositFinderKG::FdEnergyDepositFinder::FitProfile	(	fevt::Eye &	eye,
		const EFitFunctionType	type
	)

private

Definition at line 251 of file FdEnergyDepositFinder.cc.

References utl::cm2, evt::ShowerFRecData::eAtmospheric, FdEnergyDepositFinderKG::eFinalFitStep, evt::ShowerFRecData::eStatistical, fwk::CoordinateSystemRegistry::eUndefined, utl::g, fevt::EyeRecData::GetFRecShower(), evt::GaisserHillas4Parameter::GetIntegral(), evt::GaisserHillas4Parameter::GetIntegralError(), fevt::Eye::GetRecData(), evt::VGaisserHillasParameter::GetXMax(), and evt::ShowerFRecData::SetEnergyCutoff().

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

bool FdEnergyDepositFinderKG::FdEnergyDepositFinder::GuessGHParameters ( evt::ShowerFRecData & shower ) const

private

Definition at line 735 of file FdEnergyDepositFinder.cc.

References utl::cm2, FdEnergyDepositFinderKG::eOne, FdEnergyDepositFinderKG::eTwo, galactic::eV, utl::g, evt::ShowerFRecData::GetEnergyDeposit(), evt::ShowerFRecData::GetGHParameters(), evt::ShowerFRecData::HasGHParameters(), evt::ShowerFRecData::MakeGHParameters(), evt::VGaisserHillasParameter::SetNMax(), evt::GaisserHillas4Parameter::SetShapeParameter(), evt::VGaisserHillasParameter::SetXMax(), utl::Sqr(), and WARNING.

VModule::ResultFlag FdEnergyDepositFinderKG::FdEnergyDepositFinder::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 50 of file FdEnergyDepositFinder.cc.

References utl::deg, utl::InvisibleEnergy::eDATA, utl::InvisibleEnergy::eData, RdREASSimPreparator::eIron, utl::InvisibleEnergy::eMixedComposition, RdREASSimPreparator::eProton, utl::InvisibleEnergy::eQGSJET01, ERROR, utl::InvisibleEnergy::eSIBYLL21, io::eSuccess, and INFO.

void fwk::VModule::InitTiming ( )

inlineinherited

Definition at line 95 of file VModule.h.

References fwk::VModule::fStopwatch, and utl::Stopwatch::Reset().

bool FdEnergyDepositFinderKG::FdEnergyDepositFinder::IsEventCandidate ( const fevt::Eye & eye ) const

private

Definition at line 866 of file FdEnergyDepositFinder.cc.

References fwk::LocalCoordinateSystemConstructor< Policy >::Create(), degree, fevt::EyeRecData::GetFRecShower(), fevt::Eye::GetId(), fevt::EyeRecData::GetPCGF(), fevt::Eye::GetRecData(), fevt::EyeRecData::GetRp(), fevt::EyeRecData::GetRpError(), fevt::EyeRecData::GetTimeFitChiSquare(), fevt::EyeRecData::GetTimeFitNDof(), fevt::EyeRecData::HasFRecShower(), fevt::Eye::HasRecData(), and INFO.

bool FdEnergyDepositFinderKG::FdEnergyDepositFinder::PrepareEyeCopy	(	fevt::Eye &	eye,
		const EErrorPropType	type
	)		const

private

void FdEnergyDepositFinderKG::FdEnergyDepositFinder::PropagateUncertainties ( fevt::Eye & eye ) const

private

Definition at line 331 of file FdEnergyDepositFinder.cc.

References utl::cm2, fwk::LocalCoordinateSystemConstructor< Policy >::Create(), evt::ShowerFRecData::eAtmospheric, galactic::EeV, FdEnergyDepositFinderKG::eOne, evt::ShowerFRecData::eStatistical, FdEnergyDepositFinderKG::eTwo, utl::InvisibleEnergy::Factor(), utl::InvisibleEnergy::FactorDerivative(), utl::InvisibleEnergy::FactorVariance(), utl::g, gAtmVariability, gEnergyEyeTel, gEnergyResolutionConexMissingFactor, gEnergyTimeDrif, fevt::EyeRecData::GetFRecShower(), evt::GaisserHillas4Parameter::GetIntegral(), evt::VGaisserHillasParameter::GetNMax(), evt::VGaisserHillasParameter::GetNMaxError(), fevt::Eye::GetRecData(), evt::GaisserHillas4Parameter::GetShapeParameter(), evt::GaisserHillas4Parameter::GetShapeParameterError(), evt::VGaisserHillasParameter::GetXMax(), evt::VGaisserHillasParameter::GetXMaxError(), gHorizontalUniformityError_par, gNightlyRelativeCalibration, max, evt::VGaisserHillasParameter::SetNMax(), evt::GaisserHillas4Parameter::SetShapeParameter(), evt::VGaisserHillasParameter::SetXMax(), utl::Sqr(), and sqrt().

FdEnergyDepositFinderKG::FdEnergyDepositFinder::REGISTER_MODULE	(	"FdEnergyDepositFinderKG"	,
		FdEnergyDepositFinder
	)

private

VModule::ResultFlag FdEnergyDepositFinderKG::FdEnergyDepositFinder::Run ( evt::Event & event )

overridevirtual

Run: invoked once per event.

This method is for things that should be done once per event {You must override this method in your concrete module}

Implements fwk::VModule.

Definition at line 120 of file FdEnergyDepositFinder.cc.

References utl::cm2, FdEnergyDepositFinderKG::eGH2dEdXChi2, FdEnergyDepositFinderKG::eGH4dEdXChi2, FdEnergyDepositFinderKG::eGH4LightLogLike, evt::ShowerFRecData::eStatistical, io::eSuccess, galactic::eV, utl::g, fevt::FEvent::GetEye(), evt::Event::GetFEvent(), fevt::EyeRecData::GetFRecShower(), fevt::Eye::GetRecData(), FdEnergyDepositFinderKG::gHRule, evt::Event::HasFEvent(), INFO, fevt::FEvent::MakeEye(), evt::Event::MakeFEvent(), evt::ShowerFRecData::SetEnergyCutoff(), and atm::Atmosphere::SetUncertaintyBound().

ResultFlag fwk::VModule::RunWithTiming ( evt::Event & event )

inlineinherited

Definition at line 98 of file VModule.h.

References fwk::VModule::fStopwatch, fwk::VModule::Run(), utl::Stopwatch::Start(), and utl::Stopwatch::Stop().

Referenced by fwk::RunController::DoNextInSequence().

Member Data Documentation

int FdEnergyDepositFinderKG::FdEnergyDepositFinder::fAafCorrection = 0

private

Definition at line 125 of file FdEnergyDepositFinder.h.

CFMatrixCalculator FdEnergyDepositFinderKG::FdEnergyDepositFinder::fCFMatrixCalculator

private

Definition at line 94 of file FdEnergyDepositFinder.h.

CFMatrixCalculator FdEnergyDepositFinderKG::FdEnergyDepositFinder::fCFMatrixCalculatorDense

private

Definition at line 95 of file FdEnergyDepositFinder.h.

utl::InvisibleEnergy::ECompositionModel FdEnergyDepositFinderKG::FdEnergyDepositFinder::fComposition = utl::InvisibleEnergy::eData

private

Definition at line 112 of file FdEnergyDepositFinder.h.

unsigned int FdEnergyDepositFinderKG::FdEnergyDepositFinder::fDenseMatrixDim = 0

private

Definition at line 124 of file FdEnergyDepositFinder.h.

double FdEnergyDepositFinderKG::FdEnergyDepositFinder::fElectronEnergyThreshold = 1*utl::MeV

private

Definition at line 109 of file FdEnergyDepositFinder.h.

EErrorPropStage FdEnergyDepositFinderKG::FdEnergyDepositFinder::fErrPropStage

private

Definition at line 99 of file FdEnergyDepositFinder.h.

int FdEnergyDepositFinderKG::FdEnergyDepositFinder::fErrPropVerbosity = 0

private

Definition at line 116 of file FdEnergyDepositFinder.h.

GHShapeParameters FdEnergyDepositFinderKG::FdEnergyDepositFinder::fGHShapePars

private

Definition at line 106 of file FdEnergyDepositFinder.h.

int fwk::VModule::fInfoLevel = 0

protectedinherited

Definition at line 123 of file VModule.h.

Referenced by RdChannelASCIINoiseImporterRD::RdChannelASCIINoiseImporterRD::Init(), Rd2dLDFFitter::Rd2dLDFFitter::Init(), RdChannelNoiseImporter_AERA::RdChannelNoiseImporter_AERA::Init(), and Rd2dLDFFitter::Rd2dLDFFitter::Run().

utl::InvisibleEnergy::EInteractionModel FdEnergyDepositFinderKG::FdEnergyDepositFinder::fInvisibleEnergyModel = utl::InvisibleEnergy::eDATA

private

Definition at line 113 of file FdEnergyDepositFinder.h.

int FdEnergyDepositFinderKG::FdEnergyDepositFinder::fMaxErrPropPoints = 0

private

Definition at line 118 of file FdEnergyDepositFinder.h.

double FdEnergyDepositFinderKG::FdEnergyDepositFinder::fMaxZenith = 0

private

Definition at line 121 of file FdEnergyDepositFinder.h.

bool FdEnergyDepositFinderKG::FdEnergyDepositFinder::fOnlyNonMono = false

private

Definition at line 122 of file FdEnergyDepositFinder.h.

ProfileCalculator FdEnergyDepositFinderKG::FdEnergyDepositFinder::fProfileCalculator

private

Definition at line 96 of file FdEnergyDepositFinder.h.

ProfileFitter FdEnergyDepositFinderKG::FdEnergyDepositFinder::fProfileFitter

private

Definition at line 97 of file FdEnergyDepositFinder.h.

RecResults FdEnergyDepositFinderKG::FdEnergyDepositFinder::fRecResultsMinus

private

Definition at line 103 of file FdEnergyDepositFinder.h.

RecResults FdEnergyDepositFinderKG::FdEnergyDepositFinder::fRecResultsPlus

private

Definition at line 102 of file FdEnergyDepositFinder.h.

bool FdEnergyDepositFinderKG::FdEnergyDepositFinder::fRescaleErrors = false

private

Definition at line 117 of file FdEnergyDepositFinder.h.

bool FdEnergyDepositFinderKG::FdEnergyDepositFinder::fUseNoiseBins = false

private

Definition at line 126 of file FdEnergyDepositFinder.h.

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

Public Types

Public Member Functions

Static Public Member Functions

Protected Types

Protected Attributes

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Member Enumeration Documentation

Member Function Documentation

Member Data Documentation