1 #ifndef __JFIT__JMESTIMATOR__
2 #define __JFIT__JMESTIMATOR__
13 namespace JPP {
using namespace JFIT; }
34 virtual double getRho(
const double z)
const = 0;
43 virtual double getPsi(
const double z)
const = 0;
56 virtual double getRho(
const double z)
const {
return z; }
57 virtual double getPsi(
const double z)
const {
return 1.0; }
69 virtual double getRho(
const double z)
const {
return 0.5*z*z; }
70 virtual double getPsi(
const double z)
const {
return z; }
82 virtual double getRho(
const double z)
const {
return log (1.0 + 0.5*z*z); }
83 virtual double getPsi(
const double z)
const {
return z / (1.0 + 0.5*z*z); }
95 virtual double getRho(
const double z)
const {
return sqrt(1.0 + 0.5*z*z) - 1.0; }
96 virtual double getPsi(
const double z)
const {
return 0.5 * z / sqrt(1.0 + 0.5*z*z); }
117 virtual double getRho(
const double z)
const override
119 const double w = 0.5 *
k*
k / 3.0;
123 const double u = z/
k;
124 const double v = 1.0 - u*
u;
126 return w * (1.0 - v*v*
v);
132 virtual double getPsi(
const double z)
const override
136 const double u = z/
k;
137 const double v = 1.0 - u*
u;
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 getPsi(const double z) const
Get derivative of maximum likelihood estimate.
Interface for maximum likelihood estimator (M-estimator).
virtual double getPsi(const double z) const =0
Get derivative of maximum likelihood estimate.
JMEstimatorTukey(const double __k)
Constructor.
virtual double getPsi(const double z) const
Get derivative of maximum likelihood estimate.
virtual double getRho(const double z) const =0
Get maximum likelihood estimate.
virtual double getRho(const double z) const
Get maximum likelihood estimate.
virtual double getPsi(const double z) const override
Get derivative of maximum likelihood estimate.
virtual double getRho(const double z) const
Get maximum likelihood estimate.
virtual double getPsi(const double z) const
Get derivative of maximum likelihood estimate.
virtual double getPsi(const double z) const
Get derivative of maximum likelihood estimate.
virtual double getRho(const double z) const
Get maximum likelihood estimate.
virtual ~JMEstimator()
Virtual destructor.
virtual double getRho(const double z) const
Get maximum likelihood estimate.
JMEstimator * getMEstimator(const int type)
Get M-Estimator.
Tukey's biweight M-estimator.
JMEstimator_t
Definition of the various M-Estimators available to use.
virtual double getRho(const double z) const override
Get maximum likelihood estimate.