ecd.cpp

Go to the documentation of this file.

//=============================================================================================================
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "ecd.h"
 
//=============================================================================================================
// USED NAMESPACES
//=============================================================================================================
 
using namespace Eigen;
using namespace INVERSELIB;
 
//=============================================================================================================
// DEFINE MEMBER METHODS
//=============================================================================================================
 
ECD::ECD()
: valid(false)
, time(-1)
, rd(Vector3f::Zero(3))
, Q(Vector3f::Zero(3))
, good(0)
, khi2(0)
, nfree(0)
, neval(-1)
{
}
 
//=============================================================================================================
 
ECD::ECD(const ECD& p_ECD)
: valid(p_ECD.valid)
, time(p_ECD.time)
, rd(p_ECD.rd)
, Q(p_ECD.Q)
, good(p_ECD.good)
, khi2(p_ECD.khi2)
, nfree(p_ECD.nfree)
, neval(p_ECD.neval)
{
}
 
//=============================================================================================================
 
ECD::~ECD()
{
}
 
//=============================================================================================================
 
void ECD::print(FILE *f) const
{
    if (!f || !this->valid)
        return;
 
    fprintf(f,"%6.1f %7.2f %7.2f %7.2f %6.2f %6.2f %6.2f %6.2f %6.2f %6.2f %d\n",
            1000*this->time,                              /* Time */
            1000*this->rd[0],                             /* Dipole location */
            1000*this->rd[1],
            1000*this->rd[2],
            1e9*this->Q.norm(),                           /* Dipole moment */
            1e9*this->Q[0],1e9*this->Q[1],1e9*this->Q[2],
            this->khi2/this->nfree,                       /* This is the reduced khi^2 value */
            100*this->good,                               /* Goodness of fit */
            this->neval);                                 /* Number of function evaluations required */
}