#include <RdLDFFitter.h>
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 () |
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 () |
RdLDFFitter () | |
fwk::VModule::ResultFlag | Run (evt::Event &event) |
Run: invoked once per event. More... | |
ResultFlag | RunWithTiming (evt::Event &event) |
virtual | ~RdLDFFitter () |
Static Public Member Functions | |
static std::string | GetResultFlagByName (const ResultFlag flag) |
Protected Types | |
enum | InfoLevel { eInfoNone = 0, eInfoFinal = 1, eInfoIntermediate = 2, eInfoDebug = 3 } |
Private Member Functions | |
std::vector< std::pair< double, double > > | GetAnglesToEFieldExpectation (const utl::Point core, const double a) |
std::vector< std::pair< double, double > > | GetAnglesToLorentzVector () |
double | GetEnergy (const double energy_estimator) |
fwk::VModule::ResultFlag | PlotGoodnessOfFit (const std::string folder, const std::string eventIdentifier, const utl::Point core) |
fwk::VModule::ResultFlag | PlotScan (const std::vector< std::vector< double > > &scanResult, const std::string filename, const utl::Point scanCenter, const utl::CoordinateSystemPtr cs, const utl::Point coreSd, const utl::Point coreRd, const unsigned int nStepsX, const unsigned int nStepsY, const double widthX, const double widthY, const std::vector< StationFitData > &vStationFitData, const EventFitData &eventFitData, const double a) |
ParameterLDFModel | PreLDFFit (int LDFModel, const std::vector< StationFitData > &vStationFitData, utl::Vector ShowerAxis, utl::Point CorePosition) |
ParameterScintillatorLDFModel | PreScintillatorLDFFit (const std::vector< ScintillatorFitData > &vScintillatorData, utl::Vector ShowerAxis, utl::Point CorePosition) |
REGISTER_MODULE ("RdLDFFitter", RdLDFFitter) | |
fwk::VModule::ResultFlag | SaveContours (unsigned int runNumber, const int eventId, const std::vector< std::pair< double, double > > contour, const std::vector< std::pair< double, double > > contour2, ROOT::Minuit2::MnUserCovariance cov, const double coreX, const double coreY) |
fwk::VModule::ResultFlag | SaveScan (const std::vector< std::vector< double > > &scanResult, const std::string filename, const double coreX, const double coreY, const unsigned int nStepsX, const unsigned int nStepsY, const double widthX, const double widthY) |
utl::Point | Scan (std::vector< std::vector< double > > &scanResult, const LDFLikelihoodFunction L, const utl::Point core, const unsigned int nStepsX, const unsigned int nStepsY, const double widthX, const double widthY, const double a) |
nSteps is considered to be odd. If it is even, nSteps will be increased by 1 to be odd. More... | |
utl::Point | Scan (std::vector< std::vector< double > > &scanResult, const LDFLikelihoodFunction L, const utl::Point core, const unsigned int nSteps, const double width, const double a) |
Private Attributes | |
const evt::Event * | fConstEvent |
evt::Event * | fEvent |
FitConfig | fFitConfig |
int | fInfoLevel |
utl::CoordinateSystemPtr | fLocalCS |
utl::Vector | fMagneticFieldAxis |
double | fMCUncertaintyOnEFieldVector |
unsigned int | fNStepsX |
unsigned int | fNStepsY |
std::string | fOutputFolder |
bool | fPerformScan |
bool | fPlotSigmaContours |
const short | fPlotStyle |
const revt::REvent * | fREvent |
bool | fSaveContour |
ScanCenter | fScanCenter |
utl::Vector | fShowerAxis |
double | fStepWidthX |
double | fStepWidthY |
Definition at line 67 of file RdLDFFitter.h.
|
protectedinherited |
|
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. |
|
inherited |
RdLDFFitter::RdLDFFitter::RdLDFFitter | ( | ) |
Definition at line 82 of file RdLDFFitter.cc.
|
virtual |
Definition at line 103 of file RdLDFFitter.cc.
|
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 652 of file RdLDFFitter.cc.
References fwk::VModule::eSuccess, and INFO.
|
private |
Definition at line 1131 of file RdLDFFitter.cc.
References fLocalCS, fMagneticFieldAxis, fREvent, fShowerAxis, utl::RadioGeometryUtilities::GetAngleToEFieldExpectation2D(), revt::StationRecData::GetParameter(), revt::StationRecData::GetParameterError(), rdet::Station::GetPosition(), det::Detector::GetRDetector(), rdet::RDetector::GetStation(), utl::Normalized(), revt::REvent::SignalStationsBegin(), and revt::REvent::SignalStationsEnd().
Referenced by PlotScan(), and Run().
|
private |
Definition at line 1155 of file RdLDFFitter.cc.
References utl::Angle(), fLocalCS, fMagneticFieldAxis, fREvent, fShowerAxis, utl::RadioGeometryUtilities::GetLorentzVector(), revt::StationRecData::GetParameter(), revt::StationRecData::GetParameterError(), rdet::Station::GetPosition(), det::Detector::GetRDetector(), rdet::RDetector::GetStation(), utl::Normalized(), revt::REvent::SignalStationsBegin(), and revt::REvent::SignalStationsEnd().
Referenced by PlotScan().
|
private |
Definition at line 1257 of file RdLDFFitter.cc.
References utl::eV, utl::meter, utl::micro, std::pow(), and utl::volt.
Referenced by Run().
|
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().
|
inlineinherited |
Definition at line 106 of file VModule.h.
References fwk::VModule::fStopwatch.
|
inlineinherited |
Definition at line 107 of file VModule.h.
References fwk::VModule::fStopwatch.
|
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().
|
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 107 of file RdLDFFitter.cc.
References RdLDFFitter::FitConfig::chargeExcessStrength, RdLDFFitter::eFD, RdLDFFitter::eMC, RdLDFFitter::eRadio, RdLDFFitter::eSD, fwk::VModule::eSuccess, fFitConfig, fInfoLevel, RdLDFFitter::FitConfig::fitChargeExcess, fMCUncertaintyOnEFieldVector, fNStepsX, fNStepsY, fOutputFolder, fPerformScan, fSaveContour, fScanCenter, fStepWidthX, fStepWidthY, utl::Branch::GetChild(), utl::Branch::GetData(), utl::Branch::GetTopBranch(), INFO, RdLDFFitter::FitConfig::LDFModel, RdLDFFitter::FitConfig::useChargeExcessCorrectionInLDFFit, RdLDFFitter::FitConfig::usePolarisation, RdLDFFitter::FitConfig::useScintillator, RdLDFFitter::FitConfig::useSDCoreToImproveRadioCore, and RdLDFFitter::FitConfig::useStationsWithoutSignal.
|
inlineinherited |
Definition at line 95 of file VModule.h.
References fwk::VModule::fStopwatch, and utl::Stopwatch::Reset().
|
private |
Definition at line 1183 of file RdLDFFitter.cc.
References utl::Angle(), RdLDFFitter::FitConfig::chargeExcessStrength, fwk::VModule::eSuccess, fFitConfig, fLocalCS, fMagneticFieldAxis, fREvent, fShowerAxis, utl::RadioGeometryUtilities::GetExpectedEFieldVector(), revt::StationRecData::GetParameter(), revt::StationRecData::GetParameterError(), rdet::Station::GetPosition(), det::Detector::GetRDetector(), rdet::RDetector::GetStation(), utl::Normalized(), revt::REvent::SignalStationsBegin(), and revt::REvent::SignalStationsEnd().
Referenced by Run().
|
private |
Definition at line 734 of file RdLDFFitter.cc.
References RdGeoCeLDFFitter::a, utl::deg, RdLDFFitter::eDebug, utl::EeV, fwk::VModule::eSuccess, fConstEvent, fEvent, fLocalCS, fPlotSigmaContours, fPlotStyle, GetAnglesToEFieldExpectation(), GetAnglesToLorentzVector(), evt::ShowerSRecData::GetAxis(), evt::ShowerRRecData::GetAxis(), evt::ShowerSimData::GetDirection(), evt::ShowerSimData::GetEnergy(), evt::ShowerSRecData::GetEnergy(), utl::Probability::GetFisherCDF(), utl::Probability::GetFisherPDF(), RdLDFFitter::GetFormattedNumber(), rdet::RDetector::GetFullStationList(), evt::ShowerSimData::GetLocalCoordinateSystem(), utl::BasicVector< HepVector >::GetPhi(), rdet::Station::GetPosition(), det::Detector::GetRDetector(), evt::Event::GetRecShower(), evt::ShowerRecData::GetRRecShower(), evt::Event::GetSimShower(), evt::ShowerRecData::GetSRecShower(), rdet::RDetector::GetStation(), utl::BasicVector< HepVector >::GetTheta(), utl::BasicVector< HepVector >::GetX(), utl::BasicVector< HepVector >::GetY(), utl::BasicVector< HepVector >::GetZ(), evt::Event::HasSimShower(), evt::ShowerRecData::HasSRecShower(), INFO2, utl::m, RdLDFFitter::EventFitData::SDCoreXError, RdLDFFitter::EventFitData::SDCoreXYCorrelation, RdLDFFitter::EventFitData::SDCoreYError, galactic::second, evt::ShowerRRecData::SetParameter(), and sqrt().
Referenced by Run().
|
private |
Definition at line 1298 of file RdLDFFitter.cc.
References RdLDFFitter::ParameterLDFModel::A, RdLDFFitter::ParameterLDFModel::AError, RdLDFFitter::StationFitData::antennaPosition, data, utl::RadioGeometryUtilities::GetDistanceToAxis(), RdLDFFitter::ParameterLDFModel::R0, RdLDFFitter::ParameterLDFModel::R0Error, RdLDFFitter::StationFitData::signal, RdLDFFitter::StationFitData::signalError, and WARNING.
Referenced by Run().
|
private |
Definition at line 1347 of file RdLDFFitter.cc.
References utl::g, utl::RadioGeometryUtilities::GetDistanceToAxis(), RdLDFFitter::ParameterScintillatorLDFModel::N_e, RdLDFFitter::ParameterScintillatorLDFModel::N_eError, sqrt(), and RdLDFFitter::ParameterScintillatorLDFModel::status.
Referenced by Run().
|
private |
|
virtual |
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 160 of file RdLDFFitter.cc.
References RdLDFFitter::ParameterLDFModel::A, RdLDFFitter::ParameterLDFModel::AError, RdLDFFitter::StationFitData::antennaPosition, RdLDFFitter::FitConfig::chargeExcessStrength, utl::deg, utl::degree, SdHorizontalReconstructionNS::eCoreX, SdHorizontalReconstructionNS::eCoreY, RdLDFFitter::eDebug, fwk::VModule::eFailure, RdLDFFitter::eFD, RdLDFFitter::StationFitData::efield, RdLDFFitter::eFinal, RdLDFFitter::eIntermediate, RdLDFFitter::eMC, RdLDFFitter::eRadio, RdLDFFitter::eSD, fwk::VModule::eSuccess, fConstEvent, fEvent, fFitConfig, filename, fInfoLevel, RdLDFFitter::FitConfig::fitChargeExcess, RdLDFFitter::EventFitData::fLocalCS, fLocalCS, fMagneticFieldAxis, fMCUncertaintyOnEFieldVector, fNStepsX, fNStepsY, fOutputFolder, fPerformScan, fREvent, fSaveContour, fScanCenter, RdLDFFitter::EventFitData::fShowerAxis, fShowerAxis, fStepWidthX, fStepWidthY, GetAnglesToEFieldExpectation(), evt::ShowerRRecData::GetAxis(), evt::ShowerRRecData::GetCoordinateOrigin(), evt::ShowerSRecData::GetCoreError(), evt::ShowerSRecData::GetCorePosition(), evt::ShowerRRecData::GetCorePosition(), evt::ShowerFRecData::GetCorePosition(), evt::ShowerSRecData::GetCorrelationXY(), GetEnergy(), utl::Probability::GetFisherCDF(), revt::REvent::GetHeader(), revt::Header::GetId(), revt::StationRecData::GetParameter(), evt::ShowerRRecData::GetParameter(), revt::StationRecData::GetParameterError(), utl::BasicVector< HepVector >::GetPhi(), rdet::Station::GetPosition(), det::Detector::GetRDetector(), revt::Station::GetRecData(), evt::Event::GetRecShower(), det::Detector::GetReferenceCoordinateSystem(), revt::Header::GetRunNumber(), rdet::RDetector::GetStation(), utl::BasicVector< HepVector >::GetTheta(), utl::BasicVector< HepVector >::GetX(), utl::BasicVector< HepVector >::GetY(), utl::BasicVector< HepVector >::GetZ(), evt::ShowerRecData::HasFRecShower(), revt::StationRecData::HasParameter(), evt::ShowerRRecData::HasParameter(), revt::StationRecData::HasParameterError(), evt::Event::HasREvent(), evt::Event::HasSimShower(), evt::ShowerRecData::HasSRecShower(), INFO, INFO2, utl::BasicVector< HepVector >::kSpherical, RdLDFFitter::FitConfig::LDFModel, RdLDFFitter::StationFitData::lorentzAngleError, utl::m, RdLDFFitter::ParameterScintillatorLDFModel::N_e, RdLDFFitter::ParameterScintillatorLDFModel::N_eError, un2::NKG(), OUT, PlotGoodnessOfFit(), PlotScan(), std::pow(), PreLDFFit(), PreScintillatorLDFFit(), RdLDFFitter::ParameterLDFModel::R0, RdLDFFitter::ParameterLDFModel::R0Error, result, SaveContours(), SaveScan(), Scan(), RdLDFFitter::ScintillatorFitData::scintillatorPosition, RdLDFFitter::EventFitData::SDCore, RdLDFFitter::EventFitData::SDCoreXError, RdLDFFitter::EventFitData::SDCoreXYCorrelation, RdLDFFitter::EventFitData::SDCoreYError, galactic::second, revt::StationRecData::SetParameter(), evt::ShowerRRecData::SetParameter(), evt::ShowerRRecData::SetParameterCovariance(), RdLDFFitter::StationFitData::signal, RdLDFFitter::ScintillatorFitData::signal, RdLDFFitter::StationFitData::signalError, revt::REvent::SignalStationsBegin(), revt::REvent::SignalStationsEnd(), RdLDFFitter::ScintillatorFitData::silent, sqrt(), revt::REvent::StationsBegin(), revt::REvent::StationsEnd(), RdLDFFitter::ParameterScintillatorLDFModel::status, RdLDFFitter::FitConfig::useChargeExcessCorrectionInLDFFit, RdLDFFitter::FitConfig::usePolarisation, RdLDFFitter::FitConfig::useScintillator, RdLDFFitter::FitConfig::useSDCoreToImproveRadioCore, and WARNING.
|
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().
|
private |
Definition at line 1262 of file RdLDFFitter.cc.
References fwk::VModule::eFailure, fwk::VModule::eSuccess, exists, fOutputFolder, out, and WARNING.
Referenced by Run().
|
private |
Definition at line 1229 of file RdLDFFitter.cc.
References fwk::VModule::eFailure, fwk::VModule::eSuccess, and WARNING.
Referenced by Run().
|
private |
nSteps is considered to be odd. If it is even, nSteps will be increased by 1 to be odd.
Definition at line 665 of file RdLDFFitter.cc.
References RdLDFFitter::eDebug, RdLDFFitter::eIntermediate, fLocalCS, utl::BasicVector< HepVector >::GetX(), utl::BasicVector< HepVector >::GetY(), utl::BasicVector< HepVector >::GetZ(), and INFO2.
|
inlineprivate |
Definition at line 101 of file RdLDFFitter.h.
References Scan().
|
private |
Definition at line 156 of file RdLDFFitter.h.
Referenced by PlotScan(), and Run().
|
private |
Definition at line 157 of file RdLDFFitter.h.
Referenced by PlotScan(), and Run().
|
private |
Definition at line 139 of file RdLDFFitter.h.
Referenced by Init(), PlotGoodnessOfFit(), and Run().
|
private |
Definition at line 141 of file RdLDFFitter.h.
|
private |
Definition at line 143 of file RdLDFFitter.h.
Referenced by GetAnglesToEFieldExpectation(), GetAnglesToLorentzVector(), PlotGoodnessOfFit(), PlotScan(), Run(), and Scan().
|
private |
Definition at line 159 of file RdLDFFitter.h.
Referenced by GetAnglesToEFieldExpectation(), GetAnglesToLorentzVector(), PlotGoodnessOfFit(), and Run().
|
private |
Definition at line 151 of file RdLDFFitter.h.
|
private |
Definition at line 148 of file RdLDFFitter.h.
|
private |
Definition at line 148 of file RdLDFFitter.h.
|
private |
Definition at line 161 of file RdLDFFitter.h.
Referenced by Init(), Run(), and SaveContours().
|
private |
Definition at line 146 of file RdLDFFitter.h.
|
private |
Definition at line 147 of file RdLDFFitter.h.
Referenced by PlotScan().
|
private |
Definition at line 153 of file RdLDFFitter.h.
Referenced by PlotScan().
|
private |
Definition at line 155 of file RdLDFFitter.h.
Referenced by GetAnglesToEFieldExpectation(), GetAnglesToLorentzVector(), PlotGoodnessOfFit(), and Run().
|
private |
Definition at line 145 of file RdLDFFitter.h.
|
private |
Definition at line 150 of file RdLDFFitter.h.
|
private |
Definition at line 158 of file RdLDFFitter.h.
Referenced by GetAnglesToEFieldExpectation(), GetAnglesToLorentzVector(), PlotGoodnessOfFit(), and Run().
|
private |
Definition at line 149 of file RdLDFFitter.h.
|
private |
Definition at line 149 of file RdLDFFitter.h.