Fit function.
191 const double STANDARD_DEVIATIONS = 3.0;
196 data.reserve(dataL0.size() +
201 for (JOmega3D_t::const_iterator dir =
omega.begin(); dir !=
omega.end(); ++dir) {
206 buffer_type::iterator __end =
copy(dataL1.begin(), dataL1.end(),
data.begin());
208 for (buffer_type::iterator i =
data.begin(); i != __end; ++i) {
219 partial_sort(
data.begin(), __end,
data.end(),
cmz);
229 buffer_type::iterator p = __end;
230 buffer_type::iterator q =
copy(dataL0.begin(), dataL0.end(), p);
232 for (buffer_type::iterator i = p; i != q; ++i) {
234 if (find_if(
data.begin(), __end, bind2nd(equal_to<JDAQModuleIdentifier>(), i->getModuleID())) == __end) {
238 if (match1D.count(*i,
data.begin(), __end) != 0) {
249 if (
distance(
data.begin(), __end) <= NUMBER_OF_PARAMETERS) {
257 double chi2 = numeric_limits<double>::max();
258 int NDF =
distance(
data.begin(), __end) - NUMBER_OF_PARAMETERS;
266 if (number_of_outliers > NDF - 1) {
267 number_of_outliers = NDF - 1;
270 double ymin = numeric_limits<double>::max();
272 buffer_type::iterator __end1 = __end;
274 for (
int n = 0;
n <= number_of_outliers; ++
n, --__end1) {
286 (*this)(
data.begin(), __end1);
295 if (y <= -(STANDARD_DEVIATIONS * STANDARD_DEVIATIONS)) {
297 WARNING(endl <<
"chi2(1) " << y << endl);
309 NDF =
distance(
data.begin(), __end1) - NUMBER_OF_PARAMETERS;
318 ymin -= STANDARD_DEVIATIONS * STANDARD_DEVIATIONS;
323 const int number_of_outliers = NDF - 1;
327 (*this)(
data.begin(), __end);
334 for (
int n = 0;
n <= number_of_outliers; ++
n) {
339 for (
size_t i = 0; i !=
Y.size(); ++i) {
349 if (ymax < STANDARD_DEVIATIONS * STANDARD_DEVIATIONS) {
355 this->update(
data.begin(), __end,
V);
369 if (chi2 != numeric_limits<double>::max()) {
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
void set(const JVector3D &pos, T __begin, T __end, const double alpha, const double sigma)
Set co-variance matrix.
JPosition3D & rotate_back(const JRotation3D &R)
Rotate back.
bool useL0
option for L0 hit use
std::vector< T >::difference_type distance(typename std::vector< T >::const_iterator first, typename PhysicsEvent::const_iterator< T > second)
Specialisation of STL distance.
then JShowerPostfit f $INPUT_FILE o $OUTPUT_FILE N
static struct JTRIGGER::@82 clusterizeWeight
Anonymous struct for weighed clustering of hits.
void update(const size_t k, const double value)
Update inverted matrix at given diagonal element.
double getQuality(const double chi2, const int N, const int NDF)
Get quality of fit.
static struct JRECONSTRUCTION::JMuonPrefit::@64 cmz
Auxiliary data structure for sorting of hits.
int numberOfOutliers
maximum number of outliers
double roadWidth_m
road width [m]
static const int JMUONPREFIT
static struct JTRIGGER::@80 clusterize
Anonymous structure for clustering of hits.
double gridAngle_deg
grid angle for directions [deg]
void set(const JLine1Z &track, T __begin, T __end)
Set time residual vector.
then usage $script[distance] fi case set_variable R
counter_type advance(counter_type &counter, const counter_type value, const counter_type limit=std::numeric_limits< counter_type >::max())
Advance counter.
double TMaxLocal_ns
time window for local coincidences [ns]
void invert()
Invert matrix according LDU decomposition.
JFit getFit(const JHistory &history, const JTrack3D &track, const double Q, const int NDF, const double energy=0.0, const int status=0)
Get fit.
int factoryLimit
factory limit for combinatorics
Data structure for fit of straight line paralel to z-axis.
int getCount(const T &hit)
Get hit count.
void copy(const Head &from, JHead &to)
Copy header from from to to.
static struct JTRIGGER::JHitR1::@86 compare
Auxiliary data structure for sorting of hits.
double getChi2(const double P)
Get chi2 corresponding to given probability.
double sigma_ns
time resolution [ns]
int NMaxHits
maximal number of hits
1.8.5