Lorentzian M-estimator.
More...
#include <JMEstimator.hh>
|
virtual double | getRho (const double z) const |
| Get maximum likelihood estimate. More...
|
|
virtual double | getPsi (const double z) const |
| Get derivative of maximum likelihood estimate. More...
|
|
Lorentzian M-estimator.
This estimator prouces a logarithmic dependence for large deviations.
Definition at line 79 of file JMEstimator.hh.
virtual double JFIT::JMEstimatorLorentzian::getRho |
( |
const double |
z | ) |
const |
|
inlinevirtual |
Get maximum likelihood estimate.
- Parameters
-
- Returns
- likelihood
Implements JFIT::JMEstimator.
Definition at line 82 of file JMEstimator.hh.
82 {
return log (1.0 + 0.5*z*z); }
then cat $TRIPOD_INITIAL<< EOF1 256877.5 4743716.7-2438.42 256815.5 4743395.0-2435.53 257096.2 4743636.0-2439.5EOFfiJEditDetector-a $DETECTOR_INITIAL-s"-1 addz -6.9"-o $DETECTOReval`JPrintDetector-a $DETECTOR-O SUMMARY`for STRING in ${STRINGS[*]};do set_variable MODULE`getModule-a $DETECTOR-L"$STRING 0"`JEditDetector-a $DETECTOR-M"$MODULE setz -2.9"-o $DETECTORdonecp-p $TRIPOD_INITIAL $TRIPODJAcoustics.sh $DETECTOR_IDcat > acoustics_trigger_parameters txt<< EOFQ=0.0;TMax_s=0.020;numberOfHits=90;EOFJAcousticsEventBuilder.sh $DETECTOR $RUNS[*]INPUT_FILES=(`ls KM3NeT_ ${(l:8::0::0:) DETECTOR_ID}_0 *${^RUNS}_event.root`) cd $WORKDIRif[!$HOMEDIR-ef $WORKDIR];then cp-p $HOMEDIR/$DETECTOR $WORKDIR cp-p $HOMEDIR/$TRIPOD $WORKDIR cp-p $HOMEDIR/${^INPUT_FILES}$WORKDIR cp-p $HOMEDIR/{acoustics_fit_parameters, acoustics_trigger_parameters, disable, hydrophone, mechanics, sound_velocity, tripod, waveform}.txt $WORKDIRfisource $JPP_DIR/examples/JAcoustics/acoustics-fit-toolkit.shtimer_startinitialise stage_1B > &stage log
virtual double JFIT::JMEstimatorLorentzian::getPsi |
( |
const double |
z | ) |
const |
|
inlinevirtual |
Get derivative of maximum likelihood estimate.
- Parameters
-
- Returns
- derivative
Implements JFIT::JMEstimator.
Definition at line 83 of file JMEstimator.hh.
83 {
return z / (1.0 + 0.5*z*z); }
The documentation for this struct was generated from the following file: