Fit function.
193 const double STANDARD_DEVIATIONS = 3.0;
198 data.reserve(dataL0.size() +
203 for (JOmega3D_t::const_iterator dir =
omega.begin(); dir !=
omega.end(); ++dir) {
208 buffer_type::iterator __end =
copy(dataL1.begin(), dataL1.end(),
data.begin());
210 for (buffer_type::iterator
i =
data.begin();
i != __end; ++
i) {
221 partial_sort(
data.begin(), __end,
data.end(),
cmz);
231 buffer_type::iterator p = __end;
232 buffer_type::iterator q =
copy(dataL0.begin(), dataL0.end(), p);
234 for (buffer_type::iterator
i = p;
i != q; ++
i) {
240 if (match1D.count(*
i,
data.begin(), __end) != 0) {
251 if (
distance(
data.begin(), __end) <= NUMBER_OF_PARAMETERS) {
259 double chi2 = numeric_limits<double>::max();
260 int NDF =
distance(
data.begin(), __end) - NUMBER_OF_PARAMETERS;
268 if (number_of_outliers > NDF - 1) {
269 number_of_outliers = NDF - 1;
272 double ymin = numeric_limits<double>::max();
274 buffer_type::iterator __end1 = __end;
276 for (
int n = 0;
n <= number_of_outliers; ++
n, --__end1) {
288 (*this)(
data.begin(), __end1);
297 if (y <= -(STANDARD_DEVIATIONS * STANDARD_DEVIATIONS)) {
299 WARNING(endl <<
"chi2(1) " << y << endl);
311 NDF =
distance(
data.begin(), __end1) - NUMBER_OF_PARAMETERS;
320 ymin -= STANDARD_DEVIATIONS * STANDARD_DEVIATIONS;
325 const int number_of_outliers = NDF - 1;
329 (*this)(
data.begin(), __end);
336 for (
int n = 0;
n <= number_of_outliers; ++
n) {
341 for (
size_t i = 0;
i !=
Y.size(); ++
i) {
351 if (ymax < STANDARD_DEVIATIONS * STANDARD_DEVIATIONS) {
357 this->update(
data.begin(), __end,
V);
371 if (chi2 != numeric_limits<double>::max()) {
385 JEvt::iterator __end = out.end();
387 if (
parameters.numberOfPrefits < out.size()) {
404 JEvt::iterator __p = __end;
419 out.erase(__end, out.end());
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
JPredicate< JResult_t T::*, JComparison::eq > make_predicate(JResult_t T::*member, const JResult_t value)
Helper method to create predicate for data member.
static struct JTRIGGER::clusterizeWeight clusterizeWeight
int getModuleID() const
Get module identifier.
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.
Data structure for fit parameters.
JFit getFit(const JHistory &history, const JTrack3D &track, const double Q, const int NDF, const double energy=0.0, const int status=SINGLE_STAGE)
Get fit.
*fatal Wrong number of arguments esac JCookie sh typeset Z DETECTOR typeset Z SOURCE_RUN typeset Z TARGET_RUN set_variable PARAMETERS_FILE $WORKDIR parameters
void update(const size_t k, const double value)
Update inverted matrix at given diagonal element.
static struct JRECONSTRUCTION::JMuonPrefit::cmz cmz
double getQuality(const double chi2, const int N, const int NDF)
Get quality of fit.
bool qualitySorter(const JFit &first, const JFit &second)
Comparison of fit results.
int numberOfOutliers
maximum number of outliers
static struct JTRIGGER::JHitR1::compare compare
double roadWidth_m
road width [m]
static const int JMUONPREFIT
double gridAngle_deg
grid angle for directions [deg]
void set(const JLine1Z &track, T __begin, T __end)
Set time residual vector.
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.
double getDZ() const
Get Z-slope.
int factoryLimit
factory limit for combinatorics
then JCookie sh JDataQuality D $DETECTOR_ID R
then usage $script< input file >[option[primary[working directory]]] nWhere option can be N
then if[[!-f $DETECTOR]] then JDetector sh $DETECTOR fi cat $WORKDIR trigger_parameters txt<< EOFtrigger3DMuon.enabled=1;trigger3DMuon.numberOfHits=5;trigger3DMuon.gridAngle_deg=1;ctMin=0.0;TMaxLocal_ns=15.0;EOF set_variable TRIGGEREFFICIENCY_TRIGGERED_EVENTS_ONLY INPUT_FILES=() for((i=1;$i<=$NUMBER_OF_RUNS;++i));do JSirene.sh $DETECTOR $JPP_DATA/genhen.km3net_wpd_V2_0.evt.gz $WORKDIR/sirene_ ${i}.root JTriggerEfficiency.sh $DETECTOR $DETECTOR $WORKDIR/sirene_ ${i}.root $WORKDIR/trigger_efficiency_ ${i}.root $WORKDIR/trigger_parameters.txt $JPP_DATA/PMT_parameters.txt INPUT_FILES+=($WORKDIR/trigger_efficiency_ ${i}.root) done for ANGLE_DEG in $ANGLES_DEG[*];do set_variable SIGMA_NS 3.0 set_variable OUTLIERS 3 set_variable OUTPUT_FILE $WORKDIR/matrix\[${ANGLE_DEG}\deg\].root $JPP_DIR/examples/JReconstruction-f"$INPUT_FILES[*]"-o $OUTPUT_FILE-S ${SIGMA_NS}-A ${ANGLE_DEG}-O ${OUTLIERS}-d ${DEBUG}--!fiif[[$OPTION=="plot"]];then if((0));then for H1 in h0 h1;do JPlot1D-f"$WORKDIR/matrix["${^ANGLES_DEG}" deg].root:${H1}"-y"1 2e3"-Y-L TR-T""-\^"number of events [a.u.]"-> o chi2
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.
double getChi2(const double P)
Get chi2 corresponding to given probability.
static struct JTRIGGER::clusterize clusterize
double sigma_ns
time resolution [ns]
int NMaxHits
maximal number of hits
1.8.5