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()) {
383 JEvt::iterator __end = out.end();
385 if (
parameters.numberOfPrefits < out.size()) {
402 JEvt::iterator __p = __end;
417 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.
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.
then JShowerPostfit f $INPUT_FILE o $OUTPUT_FILE N
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.
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
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.
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.
double getDZ() const
Get Z-slope.
static struct JTRIGGER::@81 clusterizeWeight
Anonymous struct for weighed clustering of hits.
int factoryLimit
factory limit for combinatorics
static struct JRECONSTRUCTION::JMuonPrefit::@63 cmz
Auxiliary data structure for sorting of hits.
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::@79 clusterize
Anonymous structure for clustering of hits.
static struct JTRIGGER::JHitR1::@85 compare
Auxiliary data structure for sorting of hits.
double sigma_ns
time resolution [ns]
int NMaxHits
maximal number of hits
1.8.5