86 double _diff = _sv1 - _sv2;
87 if(_diff < 0)
return -1;
88 if(_diff == 0)
return 0;
89 if(_diff > 0)
return 1;
107 double _diff = _sv1 - _sv2;
108 if(_diff < 0)
return -1;
109 if(_diff == 0)
return 0;
110 if(_diff > 0)
return 1;
138 double _diff = _sv1 - _sv2;
139 if(_diff < 0)
return -1;
140 if(_diff == 0)
return 0;
141 if(_diff > 0)
return 1;
159 double _diff = _sv1 - _sv2;
160 if(_diff < 0)
return -1;
161 if(_diff == 0)
return 0;
162 if(_diff > 0)
return 1;
195 double _diff = _sv1 - _sv2;
196 if(_diff < 0)
return -1;
197 if(_diff == 0)
return 0;
198 if(_diff > 0)
return 1;
214 double _diff = _sv1 - _sv2;
215 if(_diff < 0)
return -1;
216 if(_diff == 0)
return 0;
217 if(_diff > 0)
return 1;
244 double _diff = _sv1 - _sv2;
245 if(_diff < 0)
return -1;
246 if(_diff == 0)
return 0;
247 if(_diff > 0)
return 1;
265 double _diff = _sv1 - _sv2;
266 if(_diff < 0)
return -1;
267 if(_diff == 0)
return 0;
268 if(_diff > 0)
return 1;
295 double _diff = _sv1 - _sv2;
296 if(_diff < 0)
return -1;
297 if(_diff == 0)
return 0;
298 if(_diff > 0)
return 1;
316 double _diff = _sv1 - _sv2;
317 if(_diff < 0)
return -1;
318 if(_diff == 0)
return 0;
319 if(_diff > 0)
return 1;
346 double _diff = _sv1 - _sv2;
347 if(_diff < 0)
return -1;
348 if(_diff == 0)
return 0;
349 if(_diff > 0)
return 1;
367 double _diff = _sv1 - _sv2;
368 if(_diff < 0)
return -1;
369 if(_diff == 0)
return 0;
370 if(_diff > 0)
return 1;
398 double _diff = _sv1 - _sv2;
399 if(_diff < 0)
return -1;
400 if(_diff == 0)
return 0;
401 if(_diff > 0)
return 1;
421 double _diff = _sv1 - _sv2;
422 if(_diff < 0)
return -1;
423 if(_diff == 0)
return 0;
424 if(_diff > 0)
return 1;
450 double _diff = _sv1 - _sv2;
451 if(_diff < 0)
return -1;
452 if(_diff == 0)
return 0;
453 if(_diff > 0)
return 1;
473 double _diff = _sv1 - _sv2;
474 if(_diff < 0)
return -1;
475 if(_diff == 0)
return 0;
476 if(_diff > 0)
return 1;
504 double _diff = _sv1 - _sv2;
505 if(_diff < 0)
return -1;
506 if(_diff == 0)
return 0;
507 if(_diff > 0)
return 1;
529 double _diff = _sv1 - _sv2;
530 if(_diff < 0)
return -1;
531 if(_diff == 0)
return 0;
532 if(_diff > 0)
return 1;
571 double _diff = _sv1 - _sv2;
572 if(_diff < 0)
return -1;
573 if(_diff == 0)
return 0;
574 if(_diff > 0)
return 1;
596 double _diff = _sv1 - _sv2;
597 if(_diff < 0)
return -1;
598 if(_diff == 0)
return 0;
599 if(_diff > 0)
return 1;
630 double _diff = _sv1 - _sv2;
631 if(_diff < 0)
return -1;
632 if(_diff == 0)
return 0;
633 if(_diff > 0)
return 1;
655 double _diff = _sv1 - _sv2;
656 if(_diff < 0)
return -1;
657 if(_diff == 0)
return 0;
658 if(_diff > 0)
return 1;
693 double _diff = _sv1 - _sv2;
694 if(_diff < 0)
return -1;
695 if(_diff == 0)
return 0;
696 if(_diff > 0)
return 1;
714 double _diff = _sv1 - _sv2;
715 if(_diff < 0)
return -1;
716 if(_diff == 0)
return 0;
717 if(_diff > 0)
return 1;
745 double _diff = _sv1 - _sv2;
746 if(_diff < 0)
return -1;
747 if(_diff == 0)
return 0;
748 if(_diff > 0)
return 1;
766 double _diff = _sv1 - _sv2;
767 if(_diff < 0)
return -1;
768 if(_diff == 0)
return 0;
769 if(_diff > 0)
return 1;
802 double _diff = _sv1 - _sv2;
803 if(_diff < 0)
return -1;
804 if(_diff == 0)
return 0;
805 if(_diff > 0)
return 1;
821 double _diff = _sv1 - _sv2;
822 if(_diff < 0)
return -1;
823 if(_diff == 0)
return 0;
824 if(_diff > 0)
return 1;
851 double _diff = _sv1 - _sv2;
852 if(_diff < 0)
return -1;
853 if(_diff == 0)
return 0;
854 if(_diff > 0)
return 1;
872 double _diff = _sv1 - _sv2;
873 if(_diff < 0)
return -1;
874 if(_diff == 0)
return 0;
875 if(_diff > 0)
return 1;
902 double _diff = _sv1 - _sv2;
903 if(_diff < 0)
return -1;
904 if(_diff == 0)
return 0;
905 if(_diff > 0)
return 1;
923 double _diff = _sv1 - _sv2;
924 if(_diff < 0)
return -1;
925 if(_diff == 0)
return 0;
926 if(_diff > 0)
return 1;
953 double _diff = _sv1 - _sv2;
954 if(_diff < 0)
return -1;
955 if(_diff == 0)
return 0;
956 if(_diff > 0)
return 1;
974 double _diff = _sv1 - _sv2;
975 if(_diff < 0)
return -1;
976 if(_diff == 0)
return 0;
977 if(_diff > 0)
return 1;
1005 double _diff = _sv1 - _sv2;
1006 if(_diff < 0)
return -1;
1007 if(_diff == 0)
return 0;
1008 if(_diff > 0)
return 1;
1028 double _diff = _sv1 - _sv2;
1029 if(_diff < 0)
return -1;
1030 if(_diff == 0)
return 0;
1031 if(_diff > 0)
return 1;
1057 double _diff = _sv1 - _sv2;
1058 if(_diff < 0)
return -1;
1059 if(_diff == 0)
return 0;
1060 if(_diff > 0)
return 1;
1080 double _diff = _sv1 - _sv2;
1081 if(_diff < 0)
return -1;
1082 if(_diff == 0)
return 0;
1083 if(_diff > 0)
return 1;
1111 double _diff = _sv1 - _sv2;
1112 if(_diff < 0)
return -1;
1113 if(_diff == 0)
return 0;
1114 if(_diff > 0)
return 1;
1136 double _diff = _sv1 - _sv2;
1137 if(_diff < 0)
return -1;
1138 if(_diff == 0)
return 0;
1139 if(_diff > 0)
return 1;
1178 double _diff = _sv1 - _sv2;
1179 if(_diff < 0)
return -1;
1180 if(_diff == 0)
return 0;
1181 if(_diff > 0)
return 1;
1203 double _diff = _sv1 - _sv2;
1204 if(_diff < 0)
return -1;
1205 if(_diff == 0)
return 0;
1206 if(_diff > 0)
return 1;
1237 double _diff = _sv1 - _sv2;
1238 if(_diff < 0)
return -1;
1239 if(_diff == 0)
return 0;
1240 if(_diff > 0)
return 1;
1262 double _diff = _sv1 - _sv2;
1263 if(_diff < 0)
return -1;
1264 if(_diff == 0)
return 0;
1265 if(_diff > 0)
return 1;
1297 double _diff = _sv1 - _sv2;
1298 if(_diff < 0)
return -1;
1299 if(_diff == 0)
return 0;
1300 if(_diff > 0)
return 1;
1322 double _diff = _sv1 - _sv2;
1323 if(_diff < 0)
return -1;
1324 if(_diff == 0)
return 0;
1325 if(_diff > 0)
return 1;
1347 double _diff = _sv1 - _sv2;
1348 if(_diff < 0)
return -1;
1349 if(_diff == 0)
return 0;
1350 if(_diff > 0)
return 1;
1387 double _diff = _sv1 - _sv2;
1388 if(_diff < 0)
return -1;
1389 if(_diff == 0)
return 0;
1390 if(_diff > 0)
return 1;
1410 double _diff = _sv1 - _sv2;
1411 if(_diff < 0)
return -1;
1412 if(_diff == 0)
return 0;
1413 if(_diff > 0)
return 1;
1445 double _diff = _sv1 - _sv2;
1446 if(_diff < 0)
return -1;
1447 if(_diff == 0)
return 0;
1448 if(_diff > 0)
return 1;
1468 double _diff = _sv1 - _sv2;
1469 if(_diff < 0)
return -1;
1470 if(_diff == 0)
return 0;
1471 if(_diff > 0)
return 1;
1495 double _diff = _sv1 - _sv2;
1496 if(_diff < 0)
return -1;
1497 if(_diff == 0)
return 0;
1498 if(_diff > 0)
return 1;
1521 double _diff = _sv1 - _sv2;
1522 if(_diff < 0)
return -1;
1523 if(_diff == 0)
return 0;
1524 if(_diff > 0)
return 1;
1541 101, 3999.5, 4100.5);
1555 11001, -5002.5, 50002.5);
1561 11001, -5002.5, 50002.5,
1581 11001, -5002.5, 50002.5);
1587 11001, -5002.5, 50002.5,
1595 11001, -5002.5, 50002.5);
1601 11001, -5002.5, 50002.5,
1633 11001, -5002.5, 50002.5);
1653 11001, -5002.5, 50002.5);
1659 11001, -5002.5, 50002.5,
1679 11001, -5002.5, 50002.5);
1685 11001, -5002.5, 50002.5,
1693 11001, -5002.5, 50002.5);
1699 11001, -5002.5, 50002.5,
1731 11001, -5002.5, 50002.5);
1743 961, -600.625, 600.625);
1757 11001, -5002.5, 50002.5);
1761 501, -10020, 10020);
1794 _t_muPC1X_auto.
reset();
1795 _t_muPC1_sparks_union.
reset();
1796 _t_muPC1X_cluster_contiguous.
reset();
1797 _t_muPC1Y_temporal_auto.
reset();
1798 _t_muPC1X_raw_auto.
reset();
1799 _t_muPC1X_AD.
reset();
1800 _t_muPC1Y_spark_search.
reset();
1801 _t_muPC1X_cluster_contiguous_sep.
reset();
1802 _t_muPC1X_spatial_auto.
reset();
1803 _t_muPC1X_export.
reset();
1804 _t_muPC1X_SC.
reset();
1805 _t_muPC1Y_export.
reset();
1806 _t_muPC1Y_raw_auto.
reset();
1810 _t_muPC1Y_cluster_contiguous_sep.
reset();
1811 _t_muPC1_auto.
reset();
1812 _t_muPC1Y_SC_auto.
reset();
1814 _t_muPC1X_spark_search.
reset();
1815 _t_muPC1X_AD_auto.
reset();
1816 _t_muPC1X_sparks.
reset();
1817 _t_muPC1Y_sparks.
reset();
1819 _t_muPC1Y_AD.
reset();
1820 _t_muPC1_sparks_clustered.
reset();
1821 _t_muPC1_hits.
reset();
1822 _t_muPC1Y_raw.
reset();
1823 _t_muPC1Y_SC.
reset();
1824 _t_muSC_muPC1.
reset();
1825 _t_muPC1_sparks_chamberwide.
reset();
1826 _t_muPC1Y_AD_auto.
reset();
1827 _t_muPC1Y_cluster_contiguous.
reset();
1829 _t_muPC1Y_spatial_auto.
reset();
1830 _t_muPC1X_raw.
reset();
1831 _t_muPC1X_SC_auto.
reset();
1832 int _u_MMuPC1AnalysisMQL_92 = 0;
1833 double _u_MMuPC1AnalysisMQL_93 = 0;
1834 bool _u_MMuPC1AnalysisMQL_94 =
false;
1835 bool _u_MMuPC1AnalysisMQL_95 =
false;
1836 int _u_MMuPC1AnalysisMQL_96 = 0;
1837 int _u_MMuPC1AnalysisMQL_97 = 0;
1838 double _u_MMuPC1AnalysisMQL_98 = 0;
1839 bool _u_MMuPC1AnalysisMQL_99 =
false;
1840 bool _u_MMuPC1AnalysisMQL_100 =
false;
1841 int _u_MMuPC1AnalysisMQL_101 = 0;
1842 int _u_MMuPC1AnalysisMQL_102 = 0;
1843 double _u_MMuPC1AnalysisMQL_103 = 0;
1844 bool _u_MMuPC1AnalysisMQL_104 =
false;
1845 bool _u_MMuPC1AnalysisMQL_105 =
false;
1846 int _u_MMuPC1AnalysisMQL_106 = 0;
1847 int _u_MMuPC1AnalysisMQL_107 = 0;
1848 double _u_MMuPC1AnalysisMQL_108 = 0;
1849 bool _u_MMuPC1AnalysisMQL_109 =
false;
1850 bool _u_MMuPC1AnalysisMQL_110 =
false;
1851 int _u_MMuPC1AnalysisMQL_111 = 0;
1852 int _u_MMuPC1AnalysisMQL_112 = 0;
1853 double _u_MMuPC1AnalysisMQL_113 = 0;
1854 bool _u_MMuPC1AnalysisMQL_114 =
false;
1855 bool _u_MMuPC1AnalysisMQL_115 =
false;
1856 int _u_MMuPC1AnalysisMQL_116 = 0;
1857 int _u_MMuPC1AnalysisMQL_117 = 0;
1858 double _u_MMuPC1AnalysisMQL_118 = 0;
1859 bool _u_MMuPC1AnalysisMQL_119 =
false;
1860 bool _u_MMuPC1AnalysisMQL_120 =
false;
1861 int _u_MMuPC1AnalysisMQL_121 = 0;
1862 double _u_MMuPC1AnalysisMQL_122 = 0;
1863 bool _u_MMuPC1AnalysisMQL_123 =
false;
1864 bool _u_MMuPC1AnalysisMQL_124 =
false;
1865 int _u_MMuPC1AnalysisMQL_125 = 0;
1866 double _u_MMuPC1AnalysisMQL_126 = 0;
1867 bool _u_MMuPC1AnalysisMQL_127 =
false;
1868 bool _u_MMuPC1AnalysisMQL_128 =
false;
1869 int _u_MMuPC1AnalysisMQL_129 = 0;
1870 int _u_MMuPC1AnalysisMQL_130 = 0;
1871 double _u_MMuPC1AnalysisMQL_131 = 0;
1872 bool _u_MMuPC1AnalysisMQL_132 =
false;
1873 bool _u_MMuPC1AnalysisMQL_133 =
false;
1874 int _u_MMuPC1AnalysisMQL_134 = 0;
1875 int _u_MMuPC1AnalysisMQL_135 = 0;
1876 double _u_MMuPC1AnalysisMQL_136 = 0;
1877 bool _u_MMuPC1AnalysisMQL_137 =
false;
1878 bool _u_MMuPC1AnalysisMQL_138 =
false;
1879 int _u_MMuPC1AnalysisMQL_139 = 0;
1880 int _u_MMuPC1AnalysisMQL_140 = 0;
1881 double _u_MMuPC1AnalysisMQL_141 = 0;
1882 bool _u_MMuPC1AnalysisMQL_142 =
false;
1883 bool _u_MMuPC1AnalysisMQL_143 =
false;
1884 int _u_MMuPC1AnalysisMQL_144 = 0;
1885 int _u_MMuPC1AnalysisMQL_145 = 0;
1886 double _u_MMuPC1AnalysisMQL_146 = 0;
1887 bool _u_MMuPC1AnalysisMQL_147 =
false;
1888 bool _u_MMuPC1AnalysisMQL_148 =
false;
1889 int _u_MMuPC1AnalysisMQL_149 = 0;
1890 int _u_MMuPC1AnalysisMQL_150 = 0;
1891 double _u_MMuPC1AnalysisMQL_151 = 0;
1892 bool _u_MMuPC1AnalysisMQL_152 =
false;
1893 bool _u_MMuPC1AnalysisMQL_153 =
false;
1894 int _u_MMuPC1AnalysisMQL_154 = 0;
1895 double _u_MMuPC1AnalysisMQL_155 = 0;
1896 bool _u_MMuPC1AnalysisMQL_156 =
false;
1897 bool _u_MMuPC1AnalysisMQL_157 =
false;
1898 int _u_MMuPC1AnalysisMQL_158 = 0;
1899 double _u_MMuPC1AnalysisMQL_159 = 0;
1900 bool _u_MMuPC1AnalysisMQL_160 =
false;
1901 bool _u_MMuPC1AnalysisMQL_161 =
false;
1902 int _u_MMuPC1AnalysisMQL_162 = 0;
1903 double _u_MMuPC1AnalysisMQL_163 = 0;
1904 bool _u_MMuPC1AnalysisMQL_164 =
false;
1905 bool _u_MMuPC1AnalysisMQL_165 =
false;
1906 int _u_MMuPC1AnalysisMQL_166 = 0;
1907 double _u_MMuPC1AnalysisMQL_167 = 0;
1908 bool _u_MMuPC1AnalysisMQL_168 =
false;
1909 bool _u_MMuPC1AnalysisMQL_169 =
false;
1910 int _u_MMuPC1AnalysisMQL_170 = 0;
1913 for(
int _u_MMuPC1AnalysisMQL_171 = 0; _u_MMuPC1AnalysisMQL_171 < _t_hits.
size(); _u_MMuPC1AnalysisMQL_171++)
1915 double time = _t_hits[_u_MMuPC1AnalysisMQL_171].time;
1916 int parameter = _t_hits[_u_MMuPC1AnalysisMQL_171].parameter;
1917 if((parameter >= 4001) && (parameter <= 4074))
1920 double time_0 =
time;
1922 double time = time_0;
1923 int parameter = parameter_0;
1927 double time_0 =
time;
1928 int Xwire = parameter - 4000;
1929 double time = time_0;
1930 int _u_MMuPC1AnalysisMQL_172 = _t_muPC1X_raw.
next();
1931 _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_172].time = time_0;
1932 _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_172].Xwire = Xwire;
1938 double time_0 =
time;
1939 int Ywire = parameter - 4050;
1940 double time = time_0;
1941 int _u_MMuPC1AnalysisMQL_173 = _t_muPC1Y_raw.
next();
1942 _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_173].time = time_0;
1943 _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_173].Ywire = Ywire;
1952 if((parameter == 6011))
1955 double time_0 =
time;
1956 double time = time_0;
1957 int _u_MMuPC1AnalysisMQL_174 = _t_muSC.
next();
1958 _t_muSC[_u_MMuPC1AnalysisMQL_174].time = time_0;
1964 for(
int _u_MMuPC1AnalysisMQL_175 = 0; _u_MMuPC1AnalysisMQL_175 < _t_muPC1X_raw.
size(); _u_MMuPC1AnalysisMQL_175++)
1966 double time = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_175].time;
1967 int Xwire = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_175].Xwire;
1970 double time_1 = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_175].time;
1971 int Xwire_1 = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_175].Xwire;
1973 if(!_u_MMuPC1AnalysisMQL_94)
1977 printf(
"Warning: table muPC1X_raw is not sorted on time_2 as required for join with muPC1X_raw (_u_MMuPC1AnalysisMQL_40)\n");
1979 _u_MMuPC1AnalysisMQL_94 =
true;
1981 else if(!_u_MMuPC1AnalysisMQL_95)
1983 if(!(time_1 >= _u_MMuPC1AnalysisMQL_93))
1985 printf(
"Warning: table muPC1X_raw is not sorted on time_1 as required for join with muPC1X_raw\n");
1986 _u_MMuPC1AnalysisMQL_95 =
true;
1989 _u_MMuPC1AnalysisMQL_93 = time_1;
1990 for( ; _u_MMuPC1AnalysisMQL_92 < _t_muPC1X_raw.
size(); _u_MMuPC1AnalysisMQL_92++)
1992 double time_2 = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_92].time;
1993 int Xwire_2 = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_92].Xwire;
1995 double _u_MMuPC1AnalysisMQL_176 = time_1 - time_2;
1996 if(_u_MMuPC1AnalysisMQL_176 < 1)
2001 for(
int _u_MMuPC1AnalysisMQL_178 = _u_MMuPC1AnalysisMQL_92; _u_MMuPC1AnalysisMQL_178 < _t_muPC1X_raw.
size(); _u_MMuPC1AnalysisMQL_178++)
2003 double time_2 = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_178].time;
2004 int Xwire_2 = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_178].Xwire;
2006 double _u_MMuPC1AnalysisMQL_176 = time_1 - time_2;
2007 if(_u_MMuPC1AnalysisMQL_176 < -50005)
2012 double time_1_0 = time_1;
2013 int Xwire_1_0 = Xwire_1;
2014 double time_2_0 = time_2;
2015 int Xwire_2_0 = Xwire_2;
2016 double time_1 = time_1_0;
2017 int Xwire_1 = Xwire_1_0;
2018 double time_2 = time_2_0;
2019 int Xwire_2 = Xwire_2_0;
2026 if((time_2 - time_1) >= 0)
2033 if(((time_2 - time_1) >= 0) && (Xwire_1==Xwire_2))
2043 bool _u_MMuPC1AnalysisMQL_179[_t_muPC1X_raw.
size()];
2044 memset(_u_MMuPC1AnalysisMQL_179, 0,
sizeof(_u_MMuPC1AnalysisMQL_179));
2047 printf(
"Warning: table muPC1X_raw is not sorted on time_1 as required for cluster.\n");
2049 for(
int _u_MMuPC1AnalysisMQL_180 = 0; _u_MMuPC1AnalysisMQL_180 < _t_muPC1X_raw.
size(); _u_MMuPC1AnalysisMQL_180++)
2051 double time = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_180].time;
2052 int Xwire = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_180].Xwire;
2053 if(_u_MMuPC1AnalysisMQL_179[_u_MMuPC1AnalysisMQL_180])
2057 _u_MMuPC1AnalysisMQL_179[_u_MMuPC1AnalysisMQL_180] =
true;
2058 int cluster_size = 1;
2059 int cluster_members[_t_muPC1X_raw.
size()];
2060 cluster_members[0] = _u_MMuPC1AnalysisMQL_180;
2062 for(
int _u_MMuPC1AnalysisMQL_181 = _u_MMuPC1AnalysisMQL_180 + 1; _u_MMuPC1AnalysisMQL_181 < _t_muPC1X_raw.
size(); _u_MMuPC1AnalysisMQL_181++)
2064 if(_u_MMuPC1AnalysisMQL_179[_u_MMuPC1AnalysisMQL_181])
2068 bool coincidence_seen =
false;
2069 for(
int _u_MMuPC1AnalysisMQL_182 = 0; _u_MMuPC1AnalysisMQL_182 < cluster_size; _u_MMuPC1AnalysisMQL_182++)
2071 int _u_MMuPC1AnalysisMQL_183 = cluster_members[_u_MMuPC1AnalysisMQL_182];
2072 double time_1 = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_183].time;
2073 double time_2 = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_181].time;
2074 int Xwire_1 = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_183].Xwire;
2075 int Xwire_2 = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_181].Xwire;
2077 double _u_MMuPC1AnalysisMQL_184 = time_1 - time_2;
2078 if(_u_MMuPC1AnalysisMQL_184 < -1000)
2082 coincidence_seen =
true;
2083 cluster_members[cluster_size++] = _u_MMuPC1AnalysisMQL_181;
2084 _u_MMuPC1AnalysisMQL_179[_u_MMuPC1AnalysisMQL_181] =
true;
2085 _u_MMuPC1AnalysisMQL_181 = _u_MMuPC1AnalysisMQL_180;
2088 if(!coincidence_seen)
2093 double time_min = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_180].time;
2094 double time_max = time_min;
2095 double time_sum = 0;
2096 int Xwire_min = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_180].Xwire;
2097 int Xwire_max = Xwire_min;
2099 for(
int _u_MMuPC1AnalysisMQL_182 = 0; _u_MMuPC1AnalysisMQL_182 < cluster_size; _u_MMuPC1AnalysisMQL_182++)
2101 int _u_MMuPC1AnalysisMQL_183 = cluster_members[_u_MMuPC1AnalysisMQL_182];
2102 time = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_183].time;
2103 time_min =
MIN(time_min, time);
2104 time_max =
MAX(time_max, time);
2105 time_sum = time_sum +
time;
2106 Xwire = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_183].Xwire;
2107 Xwire_min =
MIN(Xwire_min, Xwire);
2108 Xwire_max =
MAX(Xwire_max, Xwire);
2109 Xwire_sum = Xwire_sum + Xwire;
2111 double time_1 = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_180].time;
2112 double time_avg = time_sum/((double)cluster_size);
2113 double time_width = time_max - time_min + 1;
2114 int Xwire_1 = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_180].Xwire;
2115 double Xwire_avg = Xwire_sum/((double)cluster_size);
2116 int Xwire_width = Xwire_max - Xwire_min + 1;
2118 double time = time_min;
2119 int cluster_size_0 = cluster_size;
2120 int cluster_size = cluster_size_0;
2124 double time_0 =
time;
2125 double time = time_0;
2130 int _u_MMuPC1AnalysisMQL_185 = _t_muPC1X_sparks.
next();
2131 _t_muPC1X_sparks[_u_MMuPC1AnalysisMQL_185].time = time_0;
2141 for(
int _u_MMuPC1AnalysisMQL_186 = 0; _u_MMuPC1AnalysisMQL_186 < _t_muPC1X_raw.
size(); _u_MMuPC1AnalysisMQL_186++)
2143 double time = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_186].time;
2144 int Xwire = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_186].Xwire;
2147 double muPC1X_raw_time = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_186].time;
2148 int Xwire = _t_muPC1X_raw[_u_MMuPC1AnalysisMQL_186].Xwire;
2150 if(!_u_MMuPC1AnalysisMQL_99)
2154 printf(
"Warning: table muPC1X_sparks is not sorted on muPC1X_sparks_time as required for join with muPC1X_raw (_u_MMuPC1AnalysisMQL_44)\n");
2156 _u_MMuPC1AnalysisMQL_99 =
true;
2158 else if(!_u_MMuPC1AnalysisMQL_100)
2160 if(!(muPC1X_raw_time >= _u_MMuPC1AnalysisMQL_98))
2162 printf(
"Warning: table muPC1X_raw is not sorted on muPC1X_raw_time as required for join with muPC1X_sparks\n");
2163 _u_MMuPC1AnalysisMQL_100 =
true;
2166 _u_MMuPC1AnalysisMQL_98 = muPC1X_raw_time;
2167 for( ; _u_MMuPC1AnalysisMQL_97 < _t_muPC1X_sparks.
size(); _u_MMuPC1AnalysisMQL_97++)
2169 double muPC1X_sparks_time = _t_muPC1X_sparks[_u_MMuPC1AnalysisMQL_97].time;
2171 double _u_MMuPC1AnalysisMQL_187 = muPC1X_raw_time - muPC1X_sparks_time;
2177 _u_MMuPC1AnalysisMQL_101 = 0;
2178 for(
int _u_MMuPC1AnalysisMQL_189 = _u_MMuPC1AnalysisMQL_97; _u_MMuPC1AnalysisMQL_189 < _t_muPC1X_sparks.
size(); _u_MMuPC1AnalysisMQL_189++)
2180 double muPC1X_sparks_time = _t_muPC1X_sparks[_u_MMuPC1AnalysisMQL_189].time;
2182 double _u_MMuPC1AnalysisMQL_187 = muPC1X_raw_time - muPC1X_sparks_time;
2183 if(_u_MMuPC1AnalysisMQL_187 < -100)
2187 _u_MMuPC1AnalysisMQL_101++;
2190 if(_u_MMuPC1AnalysisMQL_101 == 0)
2193 double time = muPC1X_raw_time;
2194 int Xwire_0 = Xwire;
2195 int Xwire = Xwire_0;
2196 int _u_MMuPC1AnalysisMQL_190 = _t_muPC1X_SC.
next();
2197 _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_190].time =
time;
2198 _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_190].Xwire = Xwire_0;
2205 for(
int _u_MMuPC1AnalysisMQL_191 = 0; _u_MMuPC1AnalysisMQL_191 < _t_muPC1X_SC.
size(); _u_MMuPC1AnalysisMQL_191++)
2207 double time = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_191].time;
2208 int Xwire = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_191].Xwire;
2211 double time_1 = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_191].time;
2212 int Xwire_1 = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_191].Xwire;
2214 if(!_u_MMuPC1AnalysisMQL_104)
2218 printf(
"Warning: table muPC1X_SC is not sorted on time_2 as required for join with muPC1X_SC (_u_MMuPC1AnalysisMQL_46)\n");
2220 _u_MMuPC1AnalysisMQL_104 =
true;
2222 else if(!_u_MMuPC1AnalysisMQL_105)
2224 if(!(time_1 >= _u_MMuPC1AnalysisMQL_103))
2226 printf(
"Warning: table muPC1X_SC is not sorted on time_1 as required for join with muPC1X_SC\n");
2227 _u_MMuPC1AnalysisMQL_105 =
true;
2230 _u_MMuPC1AnalysisMQL_103 = time_1;
2231 for( ; _u_MMuPC1AnalysisMQL_102 < _t_muPC1X_SC.
size(); _u_MMuPC1AnalysisMQL_102++)
2233 double time_2 = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_102].time;
2234 int Xwire_2 = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_102].Xwire;
2236 double _u_MMuPC1AnalysisMQL_192 = time_1 - time_2;
2237 if(_u_MMuPC1AnalysisMQL_192 < 1)
2242 for(
int _u_MMuPC1AnalysisMQL_194 = _u_MMuPC1AnalysisMQL_102; _u_MMuPC1AnalysisMQL_194 < _t_muPC1X_SC.
size(); _u_MMuPC1AnalysisMQL_194++)
2244 double time_2 = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_194].time;
2245 int Xwire_2 = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_194].Xwire;
2247 double _u_MMuPC1AnalysisMQL_192 = time_1 - time_2;
2248 if(_u_MMuPC1AnalysisMQL_192 < -50005)
2253 double time_1_0 = time_1;
2254 int Xwire_1_0 = Xwire_1;
2255 double time_2_0 = time_2;
2256 int Xwire_2_0 = Xwire_2;
2257 double time_1 = time_1_0;
2258 int Xwire_1 = Xwire_1_0;
2259 double time_2 = time_2_0;
2260 int Xwire_2 = Xwire_2_0;
2263 if((time_2 - time_1) >= 0)
2270 if(((time_2 - time_1) >= 0) && (Xwire_1==Xwire_2))
2280 bool _u_MMuPC1AnalysisMQL_195[_t_muPC1X_SC.
size()];
2281 memset(_u_MMuPC1AnalysisMQL_195, 0,
sizeof(_u_MMuPC1AnalysisMQL_195));
2284 printf(
"Warning: table muPC1X_SC is not sorted on time_1 as required for cluster.\n");
2286 for(
int _u_MMuPC1AnalysisMQL_196 = 0; _u_MMuPC1AnalysisMQL_196 < _t_muPC1X_SC.
size(); _u_MMuPC1AnalysisMQL_196++)
2288 double time = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_196].time;
2289 int Xwire = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_196].Xwire;
2290 if(_u_MMuPC1AnalysisMQL_195[_u_MMuPC1AnalysisMQL_196])
2294 _u_MMuPC1AnalysisMQL_195[_u_MMuPC1AnalysisMQL_196] =
true;
2295 int cluster_size = 1;
2296 int cluster_members[_t_muPC1X_SC.
size()];
2297 cluster_members[0] = _u_MMuPC1AnalysisMQL_196;
2299 for(
int _u_MMuPC1AnalysisMQL_197 = _u_MMuPC1AnalysisMQL_196 + 1; _u_MMuPC1AnalysisMQL_197 < _t_muPC1X_SC.
size(); _u_MMuPC1AnalysisMQL_197++)
2301 if(_u_MMuPC1AnalysisMQL_195[_u_MMuPC1AnalysisMQL_197])
2305 bool coincidence_seen =
false;
2306 for(
int _u_MMuPC1AnalysisMQL_198 = 0; _u_MMuPC1AnalysisMQL_198 < cluster_size; _u_MMuPC1AnalysisMQL_198++)
2308 if(_u_MMuPC1AnalysisMQL_198 > 0)
break;
2309 int _u_MMuPC1AnalysisMQL_199 = cluster_members[_u_MMuPC1AnalysisMQL_198];
2310 double time_1 = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_199].time;
2311 double time_2 = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_197].time;
2312 int Xwire_1 = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_199].Xwire;
2313 int Xwire_2 = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_197].Xwire;
2315 double _u_MMuPC1AnalysisMQL_200 = time_1 - time_2;
2320 coincidence_seen =
true;
2321 if(!(Xwire_1 == Xwire_2))
2325 cluster_members[cluster_size++] = _u_MMuPC1AnalysisMQL_197;
2326 _u_MMuPC1AnalysisMQL_195[_u_MMuPC1AnalysisMQL_197] =
true;
2329 if(!coincidence_seen)
2334 double time_min = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_196].time;
2335 double time_max = time_min;
2336 double time_sum = 0;
2337 int Xwire_min = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_196].Xwire;
2338 int Xwire_max = Xwire_min;
2340 for(
int _u_MMuPC1AnalysisMQL_198 = 0; _u_MMuPC1AnalysisMQL_198 < cluster_size; _u_MMuPC1AnalysisMQL_198++)
2342 int _u_MMuPC1AnalysisMQL_199 = cluster_members[_u_MMuPC1AnalysisMQL_198];
2343 time = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_199].time;
2344 time_min =
MIN(time_min, time);
2345 time_max =
MAX(time_max, time);
2346 time_sum = time_sum +
time;
2347 Xwire = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_199].Xwire;
2348 Xwire_min =
MIN(Xwire_min, Xwire);
2349 Xwire_max =
MAX(Xwire_max, Xwire);
2350 Xwire_sum = Xwire_sum + Xwire;
2352 double time_1 = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_196].time;
2353 double time_avg = time_sum/((double)cluster_size);
2354 double time_width = time_max - time_min + 1;
2355 int Xwire_1 = _t_muPC1X_SC[_u_MMuPC1AnalysisMQL_196].Xwire;
2356 double Xwire_avg = Xwire_sum/((double)cluster_size);
2357 int Xwire_width = Xwire_max - Xwire_min + 1;
2359 double time = time_min;
2360 int Xwire = Xwire_1;
2361 int _u_MMuPC1AnalysisMQL_201 = _t_muPC1X_AD.
next();
2362 _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_201].time =
time;
2363 _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_201].Xwire = Xwire;
2369 for(
int _u_MMuPC1AnalysisMQL_202 = 0; _u_MMuPC1AnalysisMQL_202 < _t_muPC1X_AD.
size(); _u_MMuPC1AnalysisMQL_202++)
2371 double time = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_202].time;
2372 int Xwire = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_202].Xwire;
2375 double time_1 = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_202].time;
2376 int Xwire_1 = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_202].Xwire;
2378 if(!_u_MMuPC1AnalysisMQL_109)
2382 printf(
"Warning: table muPC1X_AD is not sorted on time_2 as required for join with muPC1X_AD (_u_MMuPC1AnalysisMQL_50)\n");
2384 _u_MMuPC1AnalysisMQL_109 =
true;
2386 else if(!_u_MMuPC1AnalysisMQL_110)
2388 if(!(time_1 >= _u_MMuPC1AnalysisMQL_108))
2390 printf(
"Warning: table muPC1X_AD is not sorted on time_1 as required for join with muPC1X_AD\n");
2391 _u_MMuPC1AnalysisMQL_110 =
true;
2394 _u_MMuPC1AnalysisMQL_108 = time_1;
2395 for( ; _u_MMuPC1AnalysisMQL_107 < _t_muPC1X_AD.
size(); _u_MMuPC1AnalysisMQL_107++)
2397 double time_2 = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_107].time;
2398 int Xwire_2 = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_107].Xwire;
2400 double _u_MMuPC1AnalysisMQL_203 = time_1 - time_2;
2401 if(_u_MMuPC1AnalysisMQL_203 < 1)
2406 for(
int _u_MMuPC1AnalysisMQL_205 = _u_MMuPC1AnalysisMQL_107; _u_MMuPC1AnalysisMQL_205 < _t_muPC1X_AD.
size(); _u_MMuPC1AnalysisMQL_205++)
2408 double time_2 = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_205].time;
2409 int Xwire_2 = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_205].Xwire;
2411 double _u_MMuPC1AnalysisMQL_203 = time_1 - time_2;
2412 if(_u_MMuPC1AnalysisMQL_203 < -50005)
2417 double time_1_0 = time_1;
2418 int Xwire_1_0 = Xwire_1;
2419 double time_2_0 = time_2;
2420 int Xwire_2_0 = Xwire_2;
2421 double time_1 = time_1_0;
2422 int Xwire_1 = Xwire_1_0;
2423 double time_2 = time_2_0;
2424 int Xwire_2 = Xwire_2_0;
2427 if((time_2 - time_1) >= 0)
2434 if(((time_2 - time_1) >= 0) && (Xwire_1==Xwire_2))
2444 bool _u_MMuPC1AnalysisMQL_206[_t_muPC1X_AD.
size()];
2445 memset(_u_MMuPC1AnalysisMQL_206, 0,
sizeof(_u_MMuPC1AnalysisMQL_206));
2448 printf(
"Warning: table muPC1X_AD is not sorted on time_1 as required for cluster.\n");
2450 for(
int _u_MMuPC1AnalysisMQL_207 = 0; _u_MMuPC1AnalysisMQL_207 < _t_muPC1X_AD.
size(); _u_MMuPC1AnalysisMQL_207++)
2452 double time = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_207].time;
2453 int Xwire = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_207].Xwire;
2454 if(_u_MMuPC1AnalysisMQL_206[_u_MMuPC1AnalysisMQL_207])
2458 _u_MMuPC1AnalysisMQL_206[_u_MMuPC1AnalysisMQL_207] =
true;
2459 int cluster_size = 1;
2460 int cluster_members[_t_muPC1X_AD.
size()];
2461 cluster_members[0] = _u_MMuPC1AnalysisMQL_207;
2463 for(
int _u_MMuPC1AnalysisMQL_208 = _u_MMuPC1AnalysisMQL_207 + 1; _u_MMuPC1AnalysisMQL_208 < _t_muPC1X_AD.
size(); _u_MMuPC1AnalysisMQL_208++)
2465 if(_u_MMuPC1AnalysisMQL_206[_u_MMuPC1AnalysisMQL_208])
2469 bool coincidence_seen =
false;
2470 for(
int _u_MMuPC1AnalysisMQL_209 = 0; _u_MMuPC1AnalysisMQL_209 < cluster_size; _u_MMuPC1AnalysisMQL_209++)
2472 int _u_MMuPC1AnalysisMQL_210 = cluster_members[_u_MMuPC1AnalysisMQL_209];
2473 double time_1 = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_210].time;
2474 double time_2 = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_208].time;
2475 int Xwire_1 = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_210].Xwire;
2476 int Xwire_2 = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_208].Xwire;
2478 double _u_MMuPC1AnalysisMQL_211 = time_1 - time_2;
2483 coincidence_seen =
true;
2484 if(!(((abs(Xwire_1-Xwire_2)-1) == 0)))
2488 cluster_members[cluster_size++] = _u_MMuPC1AnalysisMQL_208;
2489 _u_MMuPC1AnalysisMQL_206[_u_MMuPC1AnalysisMQL_208] =
true;
2490 _u_MMuPC1AnalysisMQL_208 = _u_MMuPC1AnalysisMQL_207;
2493 if(!coincidence_seen)
2498 double time_min = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_207].time;
2499 double time_max = time_min;
2500 double time_sum = 0;
2501 int Xwire_min = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_207].Xwire;
2502 int Xwire_max = Xwire_min;
2504 for(
int _u_MMuPC1AnalysisMQL_209 = 0; _u_MMuPC1AnalysisMQL_209 < cluster_size; _u_MMuPC1AnalysisMQL_209++)
2506 int _u_MMuPC1AnalysisMQL_210 = cluster_members[_u_MMuPC1AnalysisMQL_209];
2507 time = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_210].time;
2508 time_min =
MIN(time_min, time);
2509 time_max =
MAX(time_max, time);
2510 time_sum = time_sum +
time;
2511 Xwire = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_210].Xwire;
2512 Xwire_min =
MIN(Xwire_min, Xwire);
2513 Xwire_max =
MAX(Xwire_max, Xwire);
2514 Xwire_sum = Xwire_sum + Xwire;
2516 double time_1 = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_207].time;
2517 double time_avg = time_sum/((double)cluster_size);
2518 double time_width = time_max - time_min + 1;
2519 int Xwire_1 = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_207].Xwire;
2520 double Xwire_avg = Xwire_sum/((double)cluster_size);
2521 int Xwire_width = Xwire_max - Xwire_min + 1;
2523 double time = time_min;
2524 int Xwire_min_0 = Xwire_min;
2525 int Xwire_max_0 = Xwire_max;
2526 int Xwire_min = Xwire_min_0;
2527 int Xwire_max = Xwire_max_0;
2528 int _u_MMuPC1AnalysisMQL_212 = _t_muPC1X_cluster_contiguous.
next();
2529 _t_muPC1X_cluster_contiguous[_u_MMuPC1AnalysisMQL_212].time =
time;
2530 _t_muPC1X_cluster_contiguous[_u_MMuPC1AnalysisMQL_212].Xwire_min = Xwire_min_0;
2531 _t_muPC1X_cluster_contiguous[_u_MMuPC1AnalysisMQL_212].Xwire_max = Xwire_max_0;
2537 for(
int _u_MMuPC1AnalysisMQL_213 = 0; _u_MMuPC1AnalysisMQL_213 < _t_muPC1X_cluster_contiguous.
size(); _u_MMuPC1AnalysisMQL_213++)
2539 double time = _t_muPC1X_cluster_contiguous[_u_MMuPC1AnalysisMQL_213].time;
2540 int Xwire_min = _t_muPC1X_cluster_contiguous[_u_MMuPC1AnalysisMQL_213].Xwire_min;
2541 int Xwire_max = _t_muPC1X_cluster_contiguous[_u_MMuPC1AnalysisMQL_213].Xwire_max;
2544 double time_1 = _t_muPC1X_cluster_contiguous[_u_MMuPC1AnalysisMQL_213].time;
2545 int Xwire_min_1 = _t_muPC1X_cluster_contiguous[_u_MMuPC1AnalysisMQL_213].Xwire_min;
2546 int Xwire_max_1 = _t_muPC1X_cluster_contiguous[_u_MMuPC1AnalysisMQL_213].Xwire_max;
2548 if(!_u_MMuPC1AnalysisMQL_114)
2552 printf(
"Warning: table muPC1X_cluster_contiguous is not sorted on time_2 as required for join with muPC1X_cluster_contiguous (_u_MMuPC1AnalysisMQL_54)\n");
2554 _u_MMuPC1AnalysisMQL_114 =
true;
2556 else if(!_u_MMuPC1AnalysisMQL_115)
2558 if(!(time_1 >= _u_MMuPC1AnalysisMQL_113))
2560 printf(
"Warning: table muPC1X_cluster_contiguous is not sorted on time_1 as required for join with muPC1X_cluster_contiguous\n");
2561 _u_MMuPC1AnalysisMQL_115 =
true;
2564 _u_MMuPC1AnalysisMQL_113 = time_1;
2565 for( ; _u_MMuPC1AnalysisMQL_112 < _t_muPC1X_cluster_contiguous.
size(); _u_MMuPC1AnalysisMQL_112++)
2567 double time_2 = _t_muPC1X_cluster_contiguous[_u_MMuPC1AnalysisMQL_112].time;
2568 int Xwire_min_2 = _t_muPC1X_cluster_contiguous[_u_MMuPC1AnalysisMQL_112].Xwire_min;
2569 int Xwire_max_2 = _t_muPC1X_cluster_contiguous[_u_MMuPC1AnalysisMQL_112].Xwire_max;
2571 double _u_MMuPC1AnalysisMQL_214 = time_1 - time_2;
2572 if(_u_MMuPC1AnalysisMQL_214 < 5)
2577 for(
int _u_MMuPC1AnalysisMQL_216 = _u_MMuPC1AnalysisMQL_112; _u_MMuPC1AnalysisMQL_216 < _t_muPC1X_cluster_contiguous.
size(); _u_MMuPC1AnalysisMQL_216++)
2579 double time_2 = _t_muPC1X_cluster_contiguous[_u_MMuPC1AnalysisMQL_216].time;
2580 int Xwire_min_2 = _t_muPC1X_cluster_contiguous[_u_MMuPC1AnalysisMQL_216].Xwire_min;
2581 int Xwire_max_2 = _t_muPC1X_cluster_contiguous[_u_MMuPC1AnalysisMQL_216].Xwire_max;
2583 double _u_MMuPC1AnalysisMQL_214 = time_1 - time_2;
2588 if(!((time_1==time_2) && (Xwire_min_1==Xwire_min_2) && (Xwire_max_1==Xwire_max_2)) &&
2589 ((Xwire_min_1 > Xwire_max_2) || (Xwire_min_2 > Xwire_max_1)))
2592 int cluster_sep =
MIN(abs(Xwire_min_1-Xwire_max_2),abs(Xwire_min_2-Xwire_max_1));
2603 bool _u_MMuPC1AnalysisMQL_217[_t_muPC1X_AD.
size()];
2604 memset(_u_MMuPC1AnalysisMQL_217, 0,
sizeof(_u_MMuPC1AnalysisMQL_217));
2607 printf(
"Warning: table muPC1X_AD is not sorted on time_1 as required for cluster.\n");
2609 for(
int _u_MMuPC1AnalysisMQL_218 = 0; _u_MMuPC1AnalysisMQL_218 < _t_muPC1X_AD.
size(); _u_MMuPC1AnalysisMQL_218++)
2611 double time = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_218].time;
2612 int Xwire = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_218].Xwire;
2613 if(_u_MMuPC1AnalysisMQL_217[_u_MMuPC1AnalysisMQL_218])
2617 _u_MMuPC1AnalysisMQL_217[_u_MMuPC1AnalysisMQL_218] =
true;
2618 int cluster_size = 1;
2619 int cluster_members[_t_muPC1X_AD.
size()];
2620 cluster_members[0] = _u_MMuPC1AnalysisMQL_218;
2622 for(
int _u_MMuPC1AnalysisMQL_219 = _u_MMuPC1AnalysisMQL_218 + 1; _u_MMuPC1AnalysisMQL_219 < _t_muPC1X_AD.
size(); _u_MMuPC1AnalysisMQL_219++)
2624 if(_u_MMuPC1AnalysisMQL_217[_u_MMuPC1AnalysisMQL_219])
2628 bool coincidence_seen =
false;
2629 for(
int _u_MMuPC1AnalysisMQL_220 = 0; _u_MMuPC1AnalysisMQL_220 < cluster_size; _u_MMuPC1AnalysisMQL_220++)
2631 int _u_MMuPC1AnalysisMQL_221 = cluster_members[_u_MMuPC1AnalysisMQL_220];
2632 double time_1 = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_221].time;
2633 double time_2 = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_219].time;
2634 int Xwire_1 = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_221].Xwire;
2635 int Xwire_2 = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_219].Xwire;
2637 double _u_MMuPC1AnalysisMQL_222 = time_1 - time_2;
2642 coincidence_seen =
true;
2647 cluster_members[cluster_size++] = _u_MMuPC1AnalysisMQL_219;
2648 _u_MMuPC1AnalysisMQL_217[_u_MMuPC1AnalysisMQL_219] =
true;
2649 _u_MMuPC1AnalysisMQL_219 = _u_MMuPC1AnalysisMQL_218;
2652 if(!coincidence_seen)
2657 double time_min = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_218].time;
2658 double time_max = time_min;
2659 double time_sum = 0;
2660 int Xwire_min = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_218].Xwire;
2661 int Xwire_max = Xwire_min;
2663 for(
int _u_MMuPC1AnalysisMQL_220 = 0; _u_MMuPC1AnalysisMQL_220 < cluster_size; _u_MMuPC1AnalysisMQL_220++)
2665 int _u_MMuPC1AnalysisMQL_221 = cluster_members[_u_MMuPC1AnalysisMQL_220];
2666 time = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_221].time;
2667 time_min =
MIN(time_min, time);
2668 time_max =
MAX(time_max, time);
2669 time_sum = time_sum +
time;
2670 Xwire = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_221].Xwire;
2671 Xwire_min =
MIN(Xwire_min, Xwire);
2672 Xwire_max =
MAX(Xwire_max, Xwire);
2673 Xwire_sum = Xwire_sum + Xwire;
2675 double time_1 = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_218].time;
2676 double time_avg = time_sum/((double)cluster_size);
2677 double time_width = time_max - time_min + 1;
2678 int Xwire_1 = _t_muPC1X_AD[_u_MMuPC1AnalysisMQL_218].Xwire;
2679 double Xwire_avg = Xwire_sum/((double)cluster_size);
2680 int Xwire_width = Xwire_max - Xwire_min + 1;
2682 double time = time_min;
2683 double Xwire = Xwire_avg;
2684 int cluster_size_0 = cluster_size;
2685 int cluster_width = Xwire_width;
2686 int cluster_size = cluster_size_0;
2700 double Xwire_0 = Xwire;
2701 double muPC1Time =
time;
2702 double Xwire = Xwire_0;
2703 int _u_MMuPC1AnalysisMQL_223 = _t_muPC1X_export.
next();
2704 _t_muPC1X_export[_u_MMuPC1AnalysisMQL_223].Xwire = Xwire_0;
2705 _t_muPC1X_export[_u_MMuPC1AnalysisMQL_223].muPC1Time = muPC1Time;
2707 int _u_MMuPC1AnalysisMQL_224 = _t_muPC1X.
next();
2708 _t_muPC1X[_u_MMuPC1AnalysisMQL_224].time =
time;
2709 _t_muPC1X[_u_MMuPC1AnalysisMQL_224].Xwire = Xwire;
2710 _t_muPC1X[_u_MMuPC1AnalysisMQL_224].cluster_size = cluster_size_0;
2711 _t_muPC1X[_u_MMuPC1AnalysisMQL_224].cluster_width = cluster_width;
2717 for(
int _u_MMuPC1AnalysisMQL_225 = 0; _u_MMuPC1AnalysisMQL_225 < _t_muPC1X.
size(); _u_MMuPC1AnalysisMQL_225++)
2719 double time = _t_muPC1X[_u_MMuPC1AnalysisMQL_225].time;
2720 double Xwire = _t_muPC1X[_u_MMuPC1AnalysisMQL_225].Xwire;
2721 int cluster_size = _t_muPC1X[_u_MMuPC1AnalysisMQL_225].cluster_size;
2722 int cluster_width = _t_muPC1X[_u_MMuPC1AnalysisMQL_225].cluster_width;
2725 double time_1 = _t_muPC1X[_u_MMuPC1AnalysisMQL_225].time;
2726 double Xwire_1 = _t_muPC1X[_u_MMuPC1AnalysisMQL_225].Xwire;
2727 int cluster_size_1 = _t_muPC1X[_u_MMuPC1AnalysisMQL_225].cluster_size;
2728 int cluster_width_1 = _t_muPC1X[_u_MMuPC1AnalysisMQL_225].cluster_width;
2730 if(!_u_MMuPC1AnalysisMQL_119)
2734 printf(
"Warning: table muPC1X is not sorted on time_2 as required for join with muPC1X (_u_MMuPC1AnalysisMQL_58)\n");
2736 _u_MMuPC1AnalysisMQL_119 =
true;
2738 else if(!_u_MMuPC1AnalysisMQL_120)
2740 if(!(time_1 >= _u_MMuPC1AnalysisMQL_118))
2742 printf(
"Warning: table muPC1X is not sorted on time_1 as required for join with muPC1X\n");
2743 _u_MMuPC1AnalysisMQL_120 =
true;
2746 _u_MMuPC1AnalysisMQL_118 = time_1;
2747 for( ; _u_MMuPC1AnalysisMQL_117 < _t_muPC1X.
size(); _u_MMuPC1AnalysisMQL_117++)
2749 double time_2 = _t_muPC1X[_u_MMuPC1AnalysisMQL_117].time;
2750 double Xwire_2 = _t_muPC1X[_u_MMuPC1AnalysisMQL_117].Xwire;
2751 int cluster_size_2 = _t_muPC1X[_u_MMuPC1AnalysisMQL_117].cluster_size;
2752 int cluster_width_2 = _t_muPC1X[_u_MMuPC1AnalysisMQL_117].cluster_width;
2754 double _u_MMuPC1AnalysisMQL_226 = time_1 - time_2;
2755 if(_u_MMuPC1AnalysisMQL_226 < 1)
2760 for(
int _u_MMuPC1AnalysisMQL_228 = _u_MMuPC1AnalysisMQL_117; _u_MMuPC1AnalysisMQL_228 < _t_muPC1X.
size(); _u_MMuPC1AnalysisMQL_228++)
2762 double time_2 = _t_muPC1X[_u_MMuPC1AnalysisMQL_228].time;
2763 double Xwire_2 = _t_muPC1X[_u_MMuPC1AnalysisMQL_228].Xwire;
2764 int cluster_size_2 = _t_muPC1X[_u_MMuPC1AnalysisMQL_228].cluster_size;
2765 int cluster_width_2 = _t_muPC1X[_u_MMuPC1AnalysisMQL_228].cluster_width;
2767 double _u_MMuPC1AnalysisMQL_226 = time_1 - time_2;
2768 if(_u_MMuPC1AnalysisMQL_226 < -50005)
2773 double time_1_0 = time_1;
2774 double time_2_0 = time_2;
2775 double time_1 = time_1_0;
2776 double time_2 = time_2_0;
2779 if((time_2 - time_1) >= 0)
2791 for(
int _u_MMuPC1AnalysisMQL_229 = 0; _u_MMuPC1AnalysisMQL_229 < _t_muPC1X.
size(); _u_MMuPC1AnalysisMQL_229++)
2793 double time = _t_muPC1X[_u_MMuPC1AnalysisMQL_229].time;
2794 double Xwire = _t_muPC1X[_u_MMuPC1AnalysisMQL_229].Xwire;
2795 int cluster_size = _t_muPC1X[_u_MMuPC1AnalysisMQL_229].cluster_size;
2796 int cluster_width = _t_muPC1X[_u_MMuPC1AnalysisMQL_229].cluster_width;
2799 double time_1 = _t_muPC1X[_u_MMuPC1AnalysisMQL_229].time;
2800 double Xwire_1 = _t_muPC1X[_u_MMuPC1AnalysisMQL_229].Xwire;
2801 int cluster_size_1 = _t_muPC1X[_u_MMuPC1AnalysisMQL_229].cluster_size;
2802 int cluster_width_1 = _t_muPC1X[_u_MMuPC1AnalysisMQL_229].cluster_width;
2804 if(!_u_MMuPC1AnalysisMQL_123)
2808 printf(
"Warning: table muPC1X is not sorted on time_2 as required for join with muPC1X (_u_MMuPC1AnalysisMQL_60)\n");
2810 _u_MMuPC1AnalysisMQL_123 =
true;
2812 else if(!_u_MMuPC1AnalysisMQL_124)
2814 if(!(time_1 >= _u_MMuPC1AnalysisMQL_122))
2816 printf(
"Warning: table muPC1X is not sorted on time_1 as required for join with muPC1X\n");
2817 _u_MMuPC1AnalysisMQL_124 =
true;
2820 _u_MMuPC1AnalysisMQL_122 = time_1;
2821 for( ; _u_MMuPC1AnalysisMQL_121 < _t_muPC1X.
size(); _u_MMuPC1AnalysisMQL_121++)
2823 double time_2 = _t_muPC1X[_u_MMuPC1AnalysisMQL_121].time;
2824 double Xwire_2 = _t_muPC1X[_u_MMuPC1AnalysisMQL_121].Xwire;
2825 int cluster_size_2 = _t_muPC1X[_u_MMuPC1AnalysisMQL_121].cluster_size;
2826 int cluster_width_2 = _t_muPC1X[_u_MMuPC1AnalysisMQL_121].cluster_width;
2828 double _u_MMuPC1AnalysisMQL_230 = time_1 - time_2;
2829 if(_u_MMuPC1AnalysisMQL_230 < 1)
2834 for(
int _u_MMuPC1AnalysisMQL_232 = _u_MMuPC1AnalysisMQL_121; _u_MMuPC1AnalysisMQL_232 < _t_muPC1X.
size(); _u_MMuPC1AnalysisMQL_232++)
2836 double time_2 = _t_muPC1X[_u_MMuPC1AnalysisMQL_232].time;
2837 double Xwire_2 = _t_muPC1X[_u_MMuPC1AnalysisMQL_232].Xwire;
2838 int cluster_size_2 = _t_muPC1X[_u_MMuPC1AnalysisMQL_232].cluster_size;
2839 int cluster_width_2 = _t_muPC1X[_u_MMuPC1AnalysisMQL_232].cluster_width;
2841 double _u_MMuPC1AnalysisMQL_230 = time_1 - time_2;
2842 if(_u_MMuPC1AnalysisMQL_230 < -150)
2847 double time1 = time_1;
2848 double time2 = time_2;
2849 double wire_diff = fabs(Xwire_1-Xwire_2);
2852 if(fabs(time_2-time_1)>=0)
2864 for(
int _u_MMuPC1AnalysisMQL_233 = 0; _u_MMuPC1AnalysisMQL_233 < _t_muPC1Y_raw.
size(); _u_MMuPC1AnalysisMQL_233++)
2866 double time = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_233].time;
2867 int Ywire = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_233].Ywire;
2870 double time_1 = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_233].time;
2871 int Ywire_1 = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_233].Ywire;
2873 if(!_u_MMuPC1AnalysisMQL_127)
2877 printf(
"Warning: table muPC1Y_raw is not sorted on time_2 as required for join with muPC1Y_raw (_u_MMuPC1AnalysisMQL_62)\n");
2879 _u_MMuPC1AnalysisMQL_127 =
true;
2881 else if(!_u_MMuPC1AnalysisMQL_128)
2883 if(!(time_1 >= _u_MMuPC1AnalysisMQL_126))
2885 printf(
"Warning: table muPC1Y_raw is not sorted on time_1 as required for join with muPC1Y_raw\n");
2886 _u_MMuPC1AnalysisMQL_128 =
true;
2889 _u_MMuPC1AnalysisMQL_126 = time_1;
2890 for( ; _u_MMuPC1AnalysisMQL_125 < _t_muPC1Y_raw.
size(); _u_MMuPC1AnalysisMQL_125++)
2892 double time_2 = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_125].time;
2893 int Ywire_2 = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_125].Ywire;
2895 double _u_MMuPC1AnalysisMQL_234 = time_1 - time_2;
2896 if(_u_MMuPC1AnalysisMQL_234 < 1)
2901 for(
int _u_MMuPC1AnalysisMQL_236 = _u_MMuPC1AnalysisMQL_125; _u_MMuPC1AnalysisMQL_236 < _t_muPC1Y_raw.
size(); _u_MMuPC1AnalysisMQL_236++)
2903 double time_2 = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_236].time;
2904 int Ywire_2 = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_236].Ywire;
2906 double _u_MMuPC1AnalysisMQL_234 = time_1 - time_2;
2907 if(_u_MMuPC1AnalysisMQL_234 < -50005)
2912 double time_1_0 = time_1;
2913 int Ywire_1_0 = Ywire_1;
2914 double time_2_0 = time_2;
2915 int Ywire_2_0 = Ywire_2;
2916 double time_1 = time_1_0;
2917 int Ywire_1 = Ywire_1_0;
2918 double time_2 = time_2_0;
2919 int Ywire_2 = Ywire_2_0;
2926 if((time_2 - time_1) >= 0)
2933 if(((time_2 - time_1) >= 0) && (Ywire_1==Ywire_2))
2943 bool _u_MMuPC1AnalysisMQL_237[_t_muPC1Y_raw.
size()];
2944 memset(_u_MMuPC1AnalysisMQL_237, 0,
sizeof(_u_MMuPC1AnalysisMQL_237));
2947 printf(
"Warning: table muPC1Y_raw is not sorted on time_1 as required for cluster.\n");
2949 for(
int _u_MMuPC1AnalysisMQL_238 = 0; _u_MMuPC1AnalysisMQL_238 < _t_muPC1Y_raw.
size(); _u_MMuPC1AnalysisMQL_238++)
2951 double time = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_238].time;
2952 int Ywire = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_238].Ywire;
2953 if(_u_MMuPC1AnalysisMQL_237[_u_MMuPC1AnalysisMQL_238])
2957 _u_MMuPC1AnalysisMQL_237[_u_MMuPC1AnalysisMQL_238] =
true;
2958 int cluster_size = 1;
2959 int cluster_members[_t_muPC1Y_raw.
size()];
2960 cluster_members[0] = _u_MMuPC1AnalysisMQL_238;
2962 for(
int _u_MMuPC1AnalysisMQL_239 = _u_MMuPC1AnalysisMQL_238 + 1; _u_MMuPC1AnalysisMQL_239 < _t_muPC1Y_raw.
size(); _u_MMuPC1AnalysisMQL_239++)
2964 if(_u_MMuPC1AnalysisMQL_237[_u_MMuPC1AnalysisMQL_239])
2968 bool coincidence_seen =
false;
2969 for(
int _u_MMuPC1AnalysisMQL_240 = 0; _u_MMuPC1AnalysisMQL_240 < cluster_size; _u_MMuPC1AnalysisMQL_240++)
2971 int _u_MMuPC1AnalysisMQL_241 = cluster_members[_u_MMuPC1AnalysisMQL_240];
2972 double time_1 = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_241].time;
2973 double time_2 = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_239].time;
2974 int Ywire_1 = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_241].Ywire;
2975 int Ywire_2 = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_239].Ywire;
2977 double _u_MMuPC1AnalysisMQL_242 = time_1 - time_2;
2978 if(_u_MMuPC1AnalysisMQL_242 < -1000)
2982 coincidence_seen =
true;
2983 cluster_members[cluster_size++] = _u_MMuPC1AnalysisMQL_239;
2984 _u_MMuPC1AnalysisMQL_237[_u_MMuPC1AnalysisMQL_239] =
true;
2985 _u_MMuPC1AnalysisMQL_239 = _u_MMuPC1AnalysisMQL_238;
2988 if(!coincidence_seen)
2993 double time_min = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_238].time;
2994 double time_max = time_min;
2995 double time_sum = 0;
2996 int Ywire_min = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_238].Ywire;
2997 int Ywire_max = Ywire_min;
2999 for(
int _u_MMuPC1AnalysisMQL_240 = 0; _u_MMuPC1AnalysisMQL_240 < cluster_size; _u_MMuPC1AnalysisMQL_240++)
3001 int _u_MMuPC1AnalysisMQL_241 = cluster_members[_u_MMuPC1AnalysisMQL_240];
3002 time = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_241].time;
3003 time_min =
MIN(time_min, time);
3004 time_max =
MAX(time_max, time);
3005 time_sum = time_sum +
time;
3006 Ywire = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_241].Ywire;
3007 Ywire_min =
MIN(Ywire_min, Ywire);
3008 Ywire_max =
MAX(Ywire_max, Ywire);
3009 Ywire_sum = Ywire_sum + Ywire;
3011 double time_1 = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_238].time;
3012 double time_avg = time_sum/((double)cluster_size);
3013 double time_width = time_max - time_min + 1;
3014 int Ywire_1 = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_238].Ywire;
3015 double Ywire_avg = Ywire_sum/((double)cluster_size);
3016 int Ywire_width = Ywire_max - Ywire_min + 1;
3018 double time = time_min;
3019 int cluster_size_0 = cluster_size;
3020 int cluster_size = cluster_size_0;
3024 double time_0 =
time;
3025 double time = time_0;
3030 int _u_MMuPC1AnalysisMQL_243 = _t_muPC1Y_sparks.
next();
3031 _t_muPC1Y_sparks[_u_MMuPC1AnalysisMQL_243].time = time_0;
3041 for(
int _u_MMuPC1AnalysisMQL_244 = 0; _u_MMuPC1AnalysisMQL_244 < _t_muPC1Y_raw.
size(); _u_MMuPC1AnalysisMQL_244++)
3043 double time = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_244].time;
3044 int Ywire = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_244].Ywire;
3047 double muPC1Y_raw_time = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_244].time;
3048 int Ywire = _t_muPC1Y_raw[_u_MMuPC1AnalysisMQL_244].Ywire;
3050 if(!_u_MMuPC1AnalysisMQL_132)
3054 printf(
"Warning: table muPC1Y_sparks is not sorted on muPC1Y_sparks_time as required for join with muPC1Y_raw (_u_MMuPC1AnalysisMQL_66)\n");
3056 _u_MMuPC1AnalysisMQL_132 =
true;
3058 else if(!_u_MMuPC1AnalysisMQL_133)
3060 if(!(muPC1Y_raw_time >= _u_MMuPC1AnalysisMQL_131))
3062 printf(
"Warning: table muPC1Y_raw is not sorted on muPC1Y_raw_time as required for join with muPC1Y_sparks\n");
3063 _u_MMuPC1AnalysisMQL_133 =
true;
3066 _u_MMuPC1AnalysisMQL_131 = muPC1Y_raw_time;
3067 for( ; _u_MMuPC1AnalysisMQL_130 < _t_muPC1Y_sparks.
size(); _u_MMuPC1AnalysisMQL_130++)
3069 double muPC1Y_sparks_time = _t_muPC1Y_sparks[_u_MMuPC1AnalysisMQL_130].time;
3071 double _u_MMuPC1AnalysisMQL_245 = muPC1Y_raw_time - muPC1Y_sparks_time;
3077 _u_MMuPC1AnalysisMQL_134 = 0;
3078 for(
int _u_MMuPC1AnalysisMQL_247 = _u_MMuPC1AnalysisMQL_130; _u_MMuPC1AnalysisMQL_247 < _t_muPC1Y_sparks.
size(); _u_MMuPC1AnalysisMQL_247++)
3080 double muPC1Y_sparks_time = _t_muPC1Y_sparks[_u_MMuPC1AnalysisMQL_247].time;
3082 double _u_MMuPC1AnalysisMQL_245 = muPC1Y_raw_time - muPC1Y_sparks_time;
3083 if(_u_MMuPC1AnalysisMQL_245 < -100)
3087 _u_MMuPC1AnalysisMQL_134++;
3090 if(_u_MMuPC1AnalysisMQL_134 == 0)
3093 double time = muPC1Y_raw_time;
3094 int Ywire_0 = Ywire;
3095 int Ywire = Ywire_0;
3096 int _u_MMuPC1AnalysisMQL_248 = _t_muPC1Y_SC.
next();
3097 _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_248].time =
time;
3098 _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_248].Ywire = Ywire_0;
3105 for(
int _u_MMuPC1AnalysisMQL_249 = 0; _u_MMuPC1AnalysisMQL_249 < _t_muPC1Y_SC.
size(); _u_MMuPC1AnalysisMQL_249++)
3107 double time = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_249].time;
3108 int Ywire = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_249].Ywire;
3111 double time_1 = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_249].time;
3112 int Ywire_1 = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_249].Ywire;
3114 if(!_u_MMuPC1AnalysisMQL_137)
3118 printf(
"Warning: table muPC1Y_SC is not sorted on time_2 as required for join with muPC1Y_SC (_u_MMuPC1AnalysisMQL_68)\n");
3120 _u_MMuPC1AnalysisMQL_137 =
true;
3122 else if(!_u_MMuPC1AnalysisMQL_138)
3124 if(!(time_1 >= _u_MMuPC1AnalysisMQL_136))
3126 printf(
"Warning: table muPC1Y_SC is not sorted on time_1 as required for join with muPC1Y_SC\n");
3127 _u_MMuPC1AnalysisMQL_138 =
true;
3130 _u_MMuPC1AnalysisMQL_136 = time_1;
3131 for( ; _u_MMuPC1AnalysisMQL_135 < _t_muPC1Y_SC.
size(); _u_MMuPC1AnalysisMQL_135++)
3133 double time_2 = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_135].time;
3134 int Ywire_2 = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_135].Ywire;
3136 double _u_MMuPC1AnalysisMQL_250 = time_1 - time_2;
3137 if(_u_MMuPC1AnalysisMQL_250 < 1)
3142 for(
int _u_MMuPC1AnalysisMQL_252 = _u_MMuPC1AnalysisMQL_135; _u_MMuPC1AnalysisMQL_252 < _t_muPC1Y_SC.
size(); _u_MMuPC1AnalysisMQL_252++)
3144 double time_2 = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_252].time;
3145 int Ywire_2 = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_252].Ywire;
3147 double _u_MMuPC1AnalysisMQL_250 = time_1 - time_2;
3148 if(_u_MMuPC1AnalysisMQL_250 < -50005)
3153 double time_1_0 = time_1;
3154 int Ywire_1_0 = Ywire_1;
3155 double time_2_0 = time_2;
3156 int Ywire_2_0 = Ywire_2;
3157 double time_1 = time_1_0;
3158 int Ywire_1 = Ywire_1_0;
3159 double time_2 = time_2_0;
3160 int Ywire_2 = Ywire_2_0;
3163 if((time_2 - time_1) >= 0)
3170 if(((time_2 - time_1) >= 0) && (Ywire_1==Ywire_2))
3180 bool _u_MMuPC1AnalysisMQL_253[_t_muPC1Y_SC.
size()];
3181 memset(_u_MMuPC1AnalysisMQL_253, 0,
sizeof(_u_MMuPC1AnalysisMQL_253));
3184 printf(
"Warning: table muPC1Y_SC is not sorted on time_1 as required for cluster.\n");
3186 for(
int _u_MMuPC1AnalysisMQL_254 = 0; _u_MMuPC1AnalysisMQL_254 < _t_muPC1Y_SC.
size(); _u_MMuPC1AnalysisMQL_254++)
3188 double time = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_254].time;
3189 int Ywire = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_254].Ywire;
3190 if(_u_MMuPC1AnalysisMQL_253[_u_MMuPC1AnalysisMQL_254])
3194 _u_MMuPC1AnalysisMQL_253[_u_MMuPC1AnalysisMQL_254] =
true;
3195 int cluster_size = 1;
3196 int cluster_members[_t_muPC1Y_SC.
size()];
3197 cluster_members[0] = _u_MMuPC1AnalysisMQL_254;
3199 for(
int _u_MMuPC1AnalysisMQL_255 = _u_MMuPC1AnalysisMQL_254 + 1; _u_MMuPC1AnalysisMQL_255 < _t_muPC1Y_SC.
size(); _u_MMuPC1AnalysisMQL_255++)
3201 if(_u_MMuPC1AnalysisMQL_253[_u_MMuPC1AnalysisMQL_255])
3205 bool coincidence_seen =
false;
3206 for(
int _u_MMuPC1AnalysisMQL_256 = 0; _u_MMuPC1AnalysisMQL_256 < cluster_size; _u_MMuPC1AnalysisMQL_256++)
3208 if(_u_MMuPC1AnalysisMQL_256 > 0)
break;
3209 int _u_MMuPC1AnalysisMQL_257 = cluster_members[_u_MMuPC1AnalysisMQL_256];
3210 double time_1 = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_257].time;
3211 double time_2 = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_255].time;
3212 int Ywire_1 = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_257].Ywire;
3213 int Ywire_2 = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_255].Ywire;
3215 double _u_MMuPC1AnalysisMQL_258 = time_1 - time_2;
3220 coincidence_seen =
true;
3221 if(!(Ywire_1 == Ywire_2))
3225 cluster_members[cluster_size++] = _u_MMuPC1AnalysisMQL_255;
3226 _u_MMuPC1AnalysisMQL_253[_u_MMuPC1AnalysisMQL_255] =
true;
3229 if(!coincidence_seen)
3234 double time_min = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_254].time;
3235 double time_max = time_min;
3236 double time_sum = 0;
3237 int Ywire_min = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_254].Ywire;
3238 int Ywire_max = Ywire_min;
3240 for(
int _u_MMuPC1AnalysisMQL_256 = 0; _u_MMuPC1AnalysisMQL_256 < cluster_size; _u_MMuPC1AnalysisMQL_256++)
3242 int _u_MMuPC1AnalysisMQL_257 = cluster_members[_u_MMuPC1AnalysisMQL_256];
3243 time = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_257].time;
3244 time_min =
MIN(time_min, time);
3245 time_max =
MAX(time_max, time);
3246 time_sum = time_sum +
time;
3247 Ywire = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_257].Ywire;
3248 Ywire_min =
MIN(Ywire_min, Ywire);
3249 Ywire_max =
MAX(Ywire_max, Ywire);
3250 Ywire_sum = Ywire_sum + Ywire;
3252 double time_1 = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_254].time;
3253 double time_avg = time_sum/((double)cluster_size);
3254 double time_width = time_max - time_min + 1;
3255 int Ywire_1 = _t_muPC1Y_SC[_u_MMuPC1AnalysisMQL_254].Ywire;
3256 double Ywire_avg = Ywire_sum/((double)cluster_size);
3257 int Ywire_width = Ywire_max - Ywire_min + 1;
3259 double time = time_min;
3260 int Ywire = Ywire_1;
3261 int _u_MMuPC1AnalysisMQL_259 = _t_muPC1Y_AD.
next();
3262 _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_259].time =
time;
3263 _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_259].Ywire = Ywire;
3269 for(
int _u_MMuPC1AnalysisMQL_260 = 0; _u_MMuPC1AnalysisMQL_260 < _t_muPC1Y_AD.
size(); _u_MMuPC1AnalysisMQL_260++)
3271 double time = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_260].time;
3272 int Ywire = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_260].Ywire;
3275 double time_1 = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_260].time;
3276 int Ywire_1 = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_260].Ywire;
3278 if(!_u_MMuPC1AnalysisMQL_142)
3282 printf(
"Warning: table muPC1Y_AD is not sorted on time_2 as required for join with muPC1Y_AD (_u_MMuPC1AnalysisMQL_72)\n");
3284 _u_MMuPC1AnalysisMQL_142 =
true;
3286 else if(!_u_MMuPC1AnalysisMQL_143)
3288 if(!(time_1 >= _u_MMuPC1AnalysisMQL_141))
3290 printf(
"Warning: table muPC1Y_AD is not sorted on time_1 as required for join with muPC1Y_AD\n");
3291 _u_MMuPC1AnalysisMQL_143 =
true;
3294 _u_MMuPC1AnalysisMQL_141 = time_1;
3295 for( ; _u_MMuPC1AnalysisMQL_140 < _t_muPC1Y_AD.
size(); _u_MMuPC1AnalysisMQL_140++)
3297 double time_2 = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_140].time;
3298 int Ywire_2 = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_140].Ywire;
3300 double _u_MMuPC1AnalysisMQL_261 = time_1 - time_2;
3301 if(_u_MMuPC1AnalysisMQL_261 < 1)
3306 for(
int _u_MMuPC1AnalysisMQL_263 = _u_MMuPC1AnalysisMQL_140; _u_MMuPC1AnalysisMQL_263 < _t_muPC1Y_AD.
size(); _u_MMuPC1AnalysisMQL_263++)
3308 double time_2 = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_263].time;
3309 int Ywire_2 = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_263].Ywire;
3311 double _u_MMuPC1AnalysisMQL_261 = time_1 - time_2;
3312 if(_u_MMuPC1AnalysisMQL_261 < -50005)
3317 double time_1_0 = time_1;
3318 int Ywire_1_0 = Ywire_1;
3319 double time_2_0 = time_2;
3320 int Ywire_2_0 = Ywire_2;
3321 double time_1 = time_1_0;
3322 int Ywire_1 = Ywire_1_0;
3323 double time_2 = time_2_0;
3324 int Ywire_2 = Ywire_2_0;
3327 if((time_2 - time_1) >= 0)
3334 if(((time_2 - time_1) >= 0) && (Ywire_1==Ywire_2))
3344 bool _u_MMuPC1AnalysisMQL_264[_t_muPC1Y_AD.
size()];
3345 memset(_u_MMuPC1AnalysisMQL_264, 0,
sizeof(_u_MMuPC1AnalysisMQL_264));
3348 printf(
"Warning: table muPC1Y_AD is not sorted on time_1 as required for cluster.\n");
3350 for(
int _u_MMuPC1AnalysisMQL_265 = 0; _u_MMuPC1AnalysisMQL_265 < _t_muPC1Y_AD.
size(); _u_MMuPC1AnalysisMQL_265++)
3352 double time = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_265].time;
3353 int Ywire = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_265].Ywire;
3354 if(_u_MMuPC1AnalysisMQL_264[_u_MMuPC1AnalysisMQL_265])
3358 _u_MMuPC1AnalysisMQL_264[_u_MMuPC1AnalysisMQL_265] =
true;
3359 int cluster_size = 1;
3360 int cluster_members[_t_muPC1Y_AD.
size()];
3361 cluster_members[0] = _u_MMuPC1AnalysisMQL_265;
3363 for(
int _u_MMuPC1AnalysisMQL_266 = _u_MMuPC1AnalysisMQL_265 + 1; _u_MMuPC1AnalysisMQL_266 < _t_muPC1Y_AD.
size(); _u_MMuPC1AnalysisMQL_266++)
3365 if(_u_MMuPC1AnalysisMQL_264[_u_MMuPC1AnalysisMQL_266])
3369 bool coincidence_seen =
false;
3370 for(
int _u_MMuPC1AnalysisMQL_267 = 0; _u_MMuPC1AnalysisMQL_267 < cluster_size; _u_MMuPC1AnalysisMQL_267++)
3372 int _u_MMuPC1AnalysisMQL_268 = cluster_members[_u_MMuPC1AnalysisMQL_267];
3373 double time_1 = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_268].time;
3374 double time_2 = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_266].time;
3375 int Ywire_1 = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_268].Ywire;
3376 int Ywire_2 = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_266].Ywire;
3378 double _u_MMuPC1AnalysisMQL_269 = time_1 - time_2;
3383 coincidence_seen =
true;
3384 if(!(((abs(Ywire_1-Ywire_2)-1) == 0)))
3388 cluster_members[cluster_size++] = _u_MMuPC1AnalysisMQL_266;
3389 _u_MMuPC1AnalysisMQL_264[_u_MMuPC1AnalysisMQL_266] =
true;
3390 _u_MMuPC1AnalysisMQL_266 = _u_MMuPC1AnalysisMQL_265;
3393 if(!coincidence_seen)
3398 double time_min = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_265].time;
3399 double time_max = time_min;
3400 double time_sum = 0;
3401 int Ywire_min = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_265].Ywire;
3402 int Ywire_max = Ywire_min;
3404 for(
int _u_MMuPC1AnalysisMQL_267 = 0; _u_MMuPC1AnalysisMQL_267 < cluster_size; _u_MMuPC1AnalysisMQL_267++)
3406 int _u_MMuPC1AnalysisMQL_268 = cluster_members[_u_MMuPC1AnalysisMQL_267];
3407 time = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_268].time;
3408 time_min =
MIN(time_min, time);
3409 time_max =
MAX(time_max, time);
3410 time_sum = time_sum +
time;
3411 Ywire = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_268].Ywire;
3412 Ywire_min =
MIN(Ywire_min, Ywire);
3413 Ywire_max =
MAX(Ywire_max, Ywire);
3414 Ywire_sum = Ywire_sum + Ywire;
3416 double time_1 = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_265].time;
3417 double time_avg = time_sum/((double)cluster_size);
3418 double time_width = time_max - time_min + 1;
3419 int Ywire_1 = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_265].Ywire;
3420 double Ywire_avg = Ywire_sum/((double)cluster_size);
3421 int Ywire_width = Ywire_max - Ywire_min + 1;
3423 double time = time_min;
3424 int Ywire_min_0 = Ywire_min;
3425 int Ywire_max_0 = Ywire_max;
3426 int Ywire_min = Ywire_min_0;
3427 int Ywire_max = Ywire_max_0;
3428 int _u_MMuPC1AnalysisMQL_270 = _t_muPC1Y_cluster_contiguous.
next();
3429 _t_muPC1Y_cluster_contiguous[_u_MMuPC1AnalysisMQL_270].time =
time;
3430 _t_muPC1Y_cluster_contiguous[_u_MMuPC1AnalysisMQL_270].Ywire_min = Ywire_min_0;
3431 _t_muPC1Y_cluster_contiguous[_u_MMuPC1AnalysisMQL_270].Ywire_max = Ywire_max_0;
3437 for(
int _u_MMuPC1AnalysisMQL_271 = 0; _u_MMuPC1AnalysisMQL_271 < _t_muPC1Y_cluster_contiguous.
size(); _u_MMuPC1AnalysisMQL_271++)
3439 double time = _t_muPC1Y_cluster_contiguous[_u_MMuPC1AnalysisMQL_271].time;
3440 int Ywire_min = _t_muPC1Y_cluster_contiguous[_u_MMuPC1AnalysisMQL_271].Ywire_min;
3441 int Ywire_max = _t_muPC1Y_cluster_contiguous[_u_MMuPC1AnalysisMQL_271].Ywire_max;
3444 double time_1 = _t_muPC1Y_cluster_contiguous[_u_MMuPC1AnalysisMQL_271].time;
3445 int Ywire_min_1 = _t_muPC1Y_cluster_contiguous[_u_MMuPC1AnalysisMQL_271].Ywire_min;
3446 int Ywire_max_1 = _t_muPC1Y_cluster_contiguous[_u_MMuPC1AnalysisMQL_271].Ywire_max;
3448 if(!_u_MMuPC1AnalysisMQL_147)
3452 printf(
"Warning: table muPC1Y_cluster_contiguous is not sorted on time_2 as required for join with muPC1Y_cluster_contiguous (_u_MMuPC1AnalysisMQL_76)\n");
3454 _u_MMuPC1AnalysisMQL_147 =
true;
3456 else if(!_u_MMuPC1AnalysisMQL_148)
3458 if(!(time_1 >= _u_MMuPC1AnalysisMQL_146))
3460 printf(
"Warning: table muPC1Y_cluster_contiguous is not sorted on time_1 as required for join with muPC1Y_cluster_contiguous\n");
3461 _u_MMuPC1AnalysisMQL_148 =
true;
3464 _u_MMuPC1AnalysisMQL_146 = time_1;
3465 for( ; _u_MMuPC1AnalysisMQL_145 < _t_muPC1Y_cluster_contiguous.
size(); _u_MMuPC1AnalysisMQL_145++)
3467 double time_2 = _t_muPC1Y_cluster_contiguous[_u_MMuPC1AnalysisMQL_145].time;
3468 int Ywire_min_2 = _t_muPC1Y_cluster_contiguous[_u_MMuPC1AnalysisMQL_145].Ywire_min;
3469 int Ywire_max_2 = _t_muPC1Y_cluster_contiguous[_u_MMuPC1AnalysisMQL_145].Ywire_max;
3471 double _u_MMuPC1AnalysisMQL_272 = time_1 - time_2;
3472 if(_u_MMuPC1AnalysisMQL_272 < 5)
3477 for(
int _u_MMuPC1AnalysisMQL_274 = _u_MMuPC1AnalysisMQL_145; _u_MMuPC1AnalysisMQL_274 < _t_muPC1Y_cluster_contiguous.
size(); _u_MMuPC1AnalysisMQL_274++)
3479 double time_2 = _t_muPC1Y_cluster_contiguous[_u_MMuPC1AnalysisMQL_274].time;
3480 int Ywire_min_2 = _t_muPC1Y_cluster_contiguous[_u_MMuPC1AnalysisMQL_274].Ywire_min;
3481 int Ywire_max_2 = _t_muPC1Y_cluster_contiguous[_u_MMuPC1AnalysisMQL_274].Ywire_max;
3483 double _u_MMuPC1AnalysisMQL_272 = time_1 - time_2;
3488 if(!((time_1==time_2) && (Ywire_min_1==Ywire_min_2) && (Ywire_max_1==Ywire_max_2)) &&
3489 ((Ywire_min_1 > Ywire_max_2) || (Ywire_min_2 > Ywire_max_1)))
3492 int cluster_sep =
MIN(abs(Ywire_min_1-Ywire_max_2),abs(Ywire_min_2-Ywire_max_1));
3503 bool _u_MMuPC1AnalysisMQL_275[_t_muPC1Y_AD.
size()];
3504 memset(_u_MMuPC1AnalysisMQL_275, 0,
sizeof(_u_MMuPC1AnalysisMQL_275));
3507 printf(
"Warning: table muPC1Y_AD is not sorted on time_1 as required for cluster.\n");
3509 for(
int _u_MMuPC1AnalysisMQL_276 = 0; _u_MMuPC1AnalysisMQL_276 < _t_muPC1Y_AD.
size(); _u_MMuPC1AnalysisMQL_276++)
3511 double time = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_276].time;
3512 int Ywire = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_276].Ywire;
3513 if(_u_MMuPC1AnalysisMQL_275[_u_MMuPC1AnalysisMQL_276])
3517 _u_MMuPC1AnalysisMQL_275[_u_MMuPC1AnalysisMQL_276] =
true;
3518 int cluster_size = 1;
3519 int cluster_members[_t_muPC1Y_AD.
size()];
3520 cluster_members[0] = _u_MMuPC1AnalysisMQL_276;
3522 for(
int _u_MMuPC1AnalysisMQL_277 = _u_MMuPC1AnalysisMQL_276 + 1; _u_MMuPC1AnalysisMQL_277 < _t_muPC1Y_AD.
size(); _u_MMuPC1AnalysisMQL_277++)
3524 if(_u_MMuPC1AnalysisMQL_275[_u_MMuPC1AnalysisMQL_277])
3528 bool coincidence_seen =
false;
3529 for(
int _u_MMuPC1AnalysisMQL_278 = 0; _u_MMuPC1AnalysisMQL_278 < cluster_size; _u_MMuPC1AnalysisMQL_278++)
3531 int _u_MMuPC1AnalysisMQL_279 = cluster_members[_u_MMuPC1AnalysisMQL_278];
3532 double time_1 = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_279].time;
3533 double time_2 = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_277].time;
3534 int Ywire_1 = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_279].Ywire;
3535 int Ywire_2 = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_277].Ywire;
3537 double _u_MMuPC1AnalysisMQL_280 = time_1 - time_2;
3542 coincidence_seen =
true;
3547 cluster_members[cluster_size++] = _u_MMuPC1AnalysisMQL_277;
3548 _u_MMuPC1AnalysisMQL_275[_u_MMuPC1AnalysisMQL_277] =
true;
3549 _u_MMuPC1AnalysisMQL_277 = _u_MMuPC1AnalysisMQL_276;
3552 if(!coincidence_seen)
3557 double time_min = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_276].time;
3558 double time_max = time_min;
3559 double time_sum = 0;
3560 int Ywire_min = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_276].Ywire;
3561 int Ywire_max = Ywire_min;
3563 for(
int _u_MMuPC1AnalysisMQL_278 = 0; _u_MMuPC1AnalysisMQL_278 < cluster_size; _u_MMuPC1AnalysisMQL_278++)
3565 int _u_MMuPC1AnalysisMQL_279 = cluster_members[_u_MMuPC1AnalysisMQL_278];
3566 time = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_279].time;
3567 time_min =
MIN(time_min, time);
3568 time_max =
MAX(time_max, time);
3569 time_sum = time_sum +
time;
3570 Ywire = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_279].Ywire;
3571 Ywire_min =
MIN(Ywire_min, Ywire);
3572 Ywire_max =
MAX(Ywire_max, Ywire);
3573 Ywire_sum = Ywire_sum + Ywire;
3575 double time_1 = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_276].time;
3576 double time_avg = time_sum/((double)cluster_size);
3577 double time_width = time_max - time_min + 1;
3578 int Ywire_1 = _t_muPC1Y_AD[_u_MMuPC1AnalysisMQL_276].Ywire;
3579 double Ywire_avg = Ywire_sum/((double)cluster_size);
3580 int Ywire_width = Ywire_max - Ywire_min + 1;
3582 double time = time_min;
3583 double Ywire = Ywire_avg;
3584 int cluster_size_0 = cluster_size;
3585 int cluster_width = Ywire_width;
3586 int cluster_size = cluster_size_0;
3600 double Ywire_0 = Ywire;
3601 double muPC1Time =
time;
3602 double Ywire = Ywire_0;
3603 int _u_MMuPC1AnalysisMQL_281 = _t_muPC1Y_export.
next();
3604 _t_muPC1Y_export[_u_MMuPC1AnalysisMQL_281].Ywire = Ywire_0;
3605 _t_muPC1Y_export[_u_MMuPC1AnalysisMQL_281].muPC1Time = muPC1Time;
3607 int _u_MMuPC1AnalysisMQL_282 = _t_muPC1Y.
next();
3608 _t_muPC1Y[_u_MMuPC1AnalysisMQL_282].time =
time;
3609 _t_muPC1Y[_u_MMuPC1AnalysisMQL_282].Ywire = Ywire;
3610 _t_muPC1Y[_u_MMuPC1AnalysisMQL_282].cluster_size = cluster_size_0;
3611 _t_muPC1Y[_u_MMuPC1AnalysisMQL_282].cluster_width = cluster_width;
3617 for(
int _u_MMuPC1AnalysisMQL_283 = 0; _u_MMuPC1AnalysisMQL_283 < _t_muPC1Y.
size(); _u_MMuPC1AnalysisMQL_283++)
3619 double time = _t_muPC1Y[_u_MMuPC1AnalysisMQL_283].time;
3620 double Ywire = _t_muPC1Y[_u_MMuPC1AnalysisMQL_283].Ywire;
3621 int cluster_size = _t_muPC1Y[_u_MMuPC1AnalysisMQL_283].cluster_size;
3622 int cluster_width = _t_muPC1Y[_u_MMuPC1AnalysisMQL_283].cluster_width;
3625 double time_1 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_283].time;
3626 double Ywire_1 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_283].Ywire;
3627 int cluster_size_1 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_283].cluster_size;
3628 int cluster_width_1 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_283].cluster_width;
3630 if(!_u_MMuPC1AnalysisMQL_152)
3634 printf(
"Warning: table muPC1Y is not sorted on time_2 as required for join with muPC1Y (_u_MMuPC1AnalysisMQL_80)\n");
3636 _u_MMuPC1AnalysisMQL_152 =
true;
3638 else if(!_u_MMuPC1AnalysisMQL_153)
3640 if(!(time_1 >= _u_MMuPC1AnalysisMQL_151))
3642 printf(
"Warning: table muPC1Y is not sorted on time_1 as required for join with muPC1Y\n");
3643 _u_MMuPC1AnalysisMQL_153 =
true;
3646 _u_MMuPC1AnalysisMQL_151 = time_1;
3647 for( ; _u_MMuPC1AnalysisMQL_150 < _t_muPC1Y.
size(); _u_MMuPC1AnalysisMQL_150++)
3649 double time_2 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_150].time;
3650 double Ywire_2 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_150].Ywire;
3651 int cluster_size_2 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_150].cluster_size;
3652 int cluster_width_2 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_150].cluster_width;
3654 double _u_MMuPC1AnalysisMQL_284 = time_1 - time_2;
3655 if(_u_MMuPC1AnalysisMQL_284 < 1)
3660 for(
int _u_MMuPC1AnalysisMQL_286 = _u_MMuPC1AnalysisMQL_150; _u_MMuPC1AnalysisMQL_286 < _t_muPC1Y.
size(); _u_MMuPC1AnalysisMQL_286++)
3662 double time_2 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_286].time;
3663 double Ywire_2 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_286].Ywire;
3664 int cluster_size_2 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_286].cluster_size;
3665 int cluster_width_2 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_286].cluster_width;
3667 double _u_MMuPC1AnalysisMQL_284 = time_1 - time_2;
3668 if(_u_MMuPC1AnalysisMQL_284 < -50005)
3673 double time_1_0 = time_1;
3674 double time_2_0 = time_2;
3675 double time_1 = time_1_0;
3676 double time_2 = time_2_0;
3679 if((time_2 - time_1) >= 0)
3691 for(
int _u_MMuPC1AnalysisMQL_287 = 0; _u_MMuPC1AnalysisMQL_287 < _t_muPC1Y.
size(); _u_MMuPC1AnalysisMQL_287++)
3693 double time = _t_muPC1Y[_u_MMuPC1AnalysisMQL_287].time;
3694 double Ywire = _t_muPC1Y[_u_MMuPC1AnalysisMQL_287].Ywire;
3695 int cluster_size = _t_muPC1Y[_u_MMuPC1AnalysisMQL_287].cluster_size;
3696 int cluster_width = _t_muPC1Y[_u_MMuPC1AnalysisMQL_287].cluster_width;
3699 double time_1 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_287].time;
3700 double Ywire_1 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_287].Ywire;
3701 int cluster_size_1 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_287].cluster_size;
3702 int cluster_width_1 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_287].cluster_width;
3704 if(!_u_MMuPC1AnalysisMQL_156)
3708 printf(
"Warning: table muPC1Y is not sorted on time_2 as required for join with muPC1Y (_u_MMuPC1AnalysisMQL_82)\n");
3710 _u_MMuPC1AnalysisMQL_156 =
true;
3712 else if(!_u_MMuPC1AnalysisMQL_157)
3714 if(!(time_1 >= _u_MMuPC1AnalysisMQL_155))
3716 printf(
"Warning: table muPC1Y is not sorted on time_1 as required for join with muPC1Y\n");
3717 _u_MMuPC1AnalysisMQL_157 =
true;
3720 _u_MMuPC1AnalysisMQL_155 = time_1;
3721 for( ; _u_MMuPC1AnalysisMQL_154 < _t_muPC1Y.
size(); _u_MMuPC1AnalysisMQL_154++)
3723 double time_2 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_154].time;
3724 double Ywire_2 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_154].Ywire;
3725 int cluster_size_2 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_154].cluster_size;
3726 int cluster_width_2 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_154].cluster_width;
3728 double _u_MMuPC1AnalysisMQL_288 = time_1 - time_2;
3729 if(_u_MMuPC1AnalysisMQL_288 < 1)
3734 for(
int _u_MMuPC1AnalysisMQL_290 = _u_MMuPC1AnalysisMQL_154; _u_MMuPC1AnalysisMQL_290 < _t_muPC1Y.
size(); _u_MMuPC1AnalysisMQL_290++)
3736 double time_2 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_290].time;
3737 double Ywire_2 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_290].Ywire;
3738 int cluster_size_2 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_290].cluster_size;
3739 int cluster_width_2 = _t_muPC1Y[_u_MMuPC1AnalysisMQL_290].cluster_width;
3741 double _u_MMuPC1AnalysisMQL_288 = time_1 - time_2;
3742 if(_u_MMuPC1AnalysisMQL_288 < -150)
3747 double time1 = time_1;
3748 double time2 = time_2;
3749 double wire_diff = fabs(Ywire_1-Ywire_2);
3752 if(fabs(time_2-time_1)>=0)
3762 for(
int _u_MMuPC1AnalysisMQL_291 = 0; _u_MMuPC1AnalysisMQL_291 < _t_muPC1X.
size(); _u_MMuPC1AnalysisMQL_291++)
3764 double time = _t_muPC1X[_u_MMuPC1AnalysisMQL_291].time;
3765 double Xwire = _t_muPC1X[_u_MMuPC1AnalysisMQL_291].Xwire;
3766 int cluster_size = _t_muPC1X[_u_MMuPC1AnalysisMQL_291].cluster_size;
3767 int cluster_width = _t_muPC1X[_u_MMuPC1AnalysisMQL_291].cluster_width;
3770 double muPC1X_time = _t_muPC1X[_u_MMuPC1AnalysisMQL_291].time;
3771 double Xwire = _t_muPC1X[_u_MMuPC1AnalysisMQL_291].Xwire;
3772 int muPC1X_cluster_size = _t_muPC1X[_u_MMuPC1AnalysisMQL_291].cluster_size;
3773 int muPC1X_cluster_width = _t_muPC1X[_u_MMuPC1AnalysisMQL_291].cluster_width;
3775 if(!_u_MMuPC1AnalysisMQL_160)
3779 printf(
"Warning: table muPC1Y is not sorted on muPC1Y_time as required for join with muPC1X (_u_MMuPC1AnalysisMQL_84)\n");
3781 _u_MMuPC1AnalysisMQL_160 =
true;
3783 else if(!_u_MMuPC1AnalysisMQL_161)
3785 if(!(muPC1X_time >= _u_MMuPC1AnalysisMQL_159))
3787 printf(
"Warning: table muPC1X is not sorted on muPC1X_time as required for join with muPC1Y\n");
3788 _u_MMuPC1AnalysisMQL_161 =
true;
3791 _u_MMuPC1AnalysisMQL_159 = muPC1X_time;
3792 for( ; _u_MMuPC1AnalysisMQL_158 < _t_muPC1Y.
size(); _u_MMuPC1AnalysisMQL_158++)
3794 double muPC1Y_time = _t_muPC1Y[_u_MMuPC1AnalysisMQL_158].time;
3795 double Ywire = _t_muPC1Y[_u_MMuPC1AnalysisMQL_158].Ywire;
3796 int muPC1Y_cluster_size = _t_muPC1Y[_u_MMuPC1AnalysisMQL_158].cluster_size;
3797 int muPC1Y_cluster_width = _t_muPC1Y[_u_MMuPC1AnalysisMQL_158].cluster_width;
3799 double _u_MMuPC1AnalysisMQL_292 = muPC1X_time - muPC1Y_time;
3800 if(_u_MMuPC1AnalysisMQL_292 < 601)
3805 for(
int _u_MMuPC1AnalysisMQL_294 = _u_MMuPC1AnalysisMQL_158; _u_MMuPC1AnalysisMQL_294 < _t_muPC1Y.
size(); _u_MMuPC1AnalysisMQL_294++)
3807 double muPC1Y_time = _t_muPC1Y[_u_MMuPC1AnalysisMQL_294].time;
3808 double Ywire = _t_muPC1Y[_u_MMuPC1AnalysisMQL_294].Ywire;
3809 int muPC1Y_cluster_size = _t_muPC1Y[_u_MMuPC1AnalysisMQL_294].cluster_size;
3810 int muPC1Y_cluster_width = _t_muPC1Y[_u_MMuPC1AnalysisMQL_294].cluster_width;
3812 double _u_MMuPC1AnalysisMQL_292 = muPC1X_time - muPC1Y_time;
3813 if(_u_MMuPC1AnalysisMQL_292 < -601)
3818 double Xwire_0 = Xwire;
3819 double Ywire_0 = Ywire;
3820 double muPC1Time =
MIN(muPC1X_time, muPC1Y_time);
3821 double XYtdiff = muPC1Y_time - muPC1X_time;
3822 double Xwire = Xwire_0;
3823 double Ywire = Ywire_0;
3824 int _u_MMuPC1AnalysisMQL_295 = _t_muPC1XY.
next();
3825 _t_muPC1XY[_u_MMuPC1AnalysisMQL_295].Xwire = Xwire_0;
3826 _t_muPC1XY[_u_MMuPC1AnalysisMQL_295].Ywire = Ywire_0;
3827 _t_muPC1XY[_u_MMuPC1AnalysisMQL_295].muPC1Time = muPC1Time;
3828 _t_muPC1XY[_u_MMuPC1AnalysisMQL_295].XYtdiff = XYtdiff;
3834 for(
int _u_MMuPC1AnalysisMQL_296 = 0; _u_MMuPC1AnalysisMQL_296 < _t_muPC1XY.
size(); _u_MMuPC1AnalysisMQL_296++)
3836 double Xwire = _t_muPC1XY[_u_MMuPC1AnalysisMQL_296].Xwire;
3837 double Ywire = _t_muPC1XY[_u_MMuPC1AnalysisMQL_296].Ywire;
3838 double muPC1Time = _t_muPC1XY[_u_MMuPC1AnalysisMQL_296].muPC1Time;
3839 double XYtdiff = _t_muPC1XY[_u_MMuPC1AnalysisMQL_296].XYtdiff;
3844 if((XYtdiff >= -95) && (XYtdiff <= 95))
3847 double Xwire_0 = Xwire;
3848 double Ywire_0 = Ywire;
3849 double muPC1Time_0 = muPC1Time;
3850 double Xwire = Xwire_0;
3851 double Ywire = Ywire_0;
3852 double muPC1Time = muPC1Time_0;
3857 int _u_MMuPC1AnalysisMQL_297 = _t_muPC1.
next();
3858 _t_muPC1[_u_MMuPC1AnalysisMQL_297].Xwire = Xwire_0;
3859 _t_muPC1[_u_MMuPC1AnalysisMQL_297].Ywire = Ywire_0;
3860 _t_muPC1[_u_MMuPC1AnalysisMQL_297].muPC1Time = muPC1Time_0;
3866 for(
int _u_MMuPC1AnalysisMQL_298 = 0; _u_MMuPC1AnalysisMQL_298 < _t_muPC1.
size(); _u_MMuPC1AnalysisMQL_298++)
3868 double Xwire = _t_muPC1[_u_MMuPC1AnalysisMQL_298].Xwire;
3869 double Ywire = _t_muPC1[_u_MMuPC1AnalysisMQL_298].Ywire;
3870 double muPC1Time = _t_muPC1[_u_MMuPC1AnalysisMQL_298].muPC1Time;
3873 double Xwire_1 = _t_muPC1[_u_MMuPC1AnalysisMQL_298].Xwire;
3874 double Ywire_1 = _t_muPC1[_u_MMuPC1AnalysisMQL_298].Ywire;
3875 double muPC1Time_1 = _t_muPC1[_u_MMuPC1AnalysisMQL_298].muPC1Time;
3877 if(!_u_MMuPC1AnalysisMQL_164)
3881 printf(
"Warning: table muPC1 is not sorted on muPC1Time_2 as required for join with muPC1 (_u_MMuPC1AnalysisMQL_87)\n");
3883 _u_MMuPC1AnalysisMQL_164 =
true;
3885 else if(!_u_MMuPC1AnalysisMQL_165)
3887 if(!(muPC1Time_1 >= _u_MMuPC1AnalysisMQL_163))
3889 printf(
"Warning: table muPC1 is not sorted on muPC1Time_1 as required for join with muPC1\n");
3890 _u_MMuPC1AnalysisMQL_165 =
true;
3893 _u_MMuPC1AnalysisMQL_163 = muPC1Time_1;
3894 for( ; _u_MMuPC1AnalysisMQL_162 < _t_muPC1.
size(); _u_MMuPC1AnalysisMQL_162++)
3896 double Xwire_2 = _t_muPC1[_u_MMuPC1AnalysisMQL_162].Xwire;
3897 double Ywire_2 = _t_muPC1[_u_MMuPC1AnalysisMQL_162].Ywire;
3898 double muPC1Time_2 = _t_muPC1[_u_MMuPC1AnalysisMQL_162].muPC1Time;
3900 double _u_MMuPC1AnalysisMQL_299 = muPC1Time_1 - muPC1Time_2;
3901 if(_u_MMuPC1AnalysisMQL_299 < 1)
3906 for(
int _u_MMuPC1AnalysisMQL_301 = _u_MMuPC1AnalysisMQL_162; _u_MMuPC1AnalysisMQL_301 < _t_muPC1.
size(); _u_MMuPC1AnalysisMQL_301++)
3908 double Xwire_2 = _t_muPC1[_u_MMuPC1AnalysisMQL_301].Xwire;
3909 double Ywire_2 = _t_muPC1[_u_MMuPC1AnalysisMQL_301].Ywire;
3910 double muPC1Time_2 = _t_muPC1[_u_MMuPC1AnalysisMQL_301].muPC1Time;
3912 double _u_MMuPC1AnalysisMQL_299 = muPC1Time_1 - muPC1Time_2;
3913 if(_u_MMuPC1AnalysisMQL_299 < -50005)
3918 double muPC1Time_1_0 = muPC1Time_1;
3919 double muPC1Time_2_0 = muPC1Time_2;
3920 double muPC1Time_1 = muPC1Time_1_0;
3921 double muPC1Time_2 = muPC1Time_2_0;
3924 if((muPC1Time_2 - muPC1Time_1) >= 0)
3934 for(
int _u_MMuPC1AnalysisMQL_302 = 0; _u_MMuPC1AnalysisMQL_302 < _t_muSC.
size(); _u_MMuPC1AnalysisMQL_302++)
3936 double time = _t_muSC[_u_MMuPC1AnalysisMQL_302].time;
3939 double time = _t_muSC[_u_MMuPC1AnalysisMQL_302].time;
3941 if(!_u_MMuPC1AnalysisMQL_168)
3945 printf(
"Warning: table muPC1 is not sorted on muPC1Time as required for join with muSC (_u_MMuPC1AnalysisMQL_89)\n");
3947 _u_MMuPC1AnalysisMQL_168 =
true;
3949 else if(!_u_MMuPC1AnalysisMQL_169)
3951 if(!(time >= _u_MMuPC1AnalysisMQL_167))
3953 printf(
"Warning: table muSC is not sorted on time as required for join with muPC1\n");
3954 _u_MMuPC1AnalysisMQL_169 =
true;
3957 _u_MMuPC1AnalysisMQL_167 =
time;
3958 for( ; _u_MMuPC1AnalysisMQL_166 < _t_muPC1.
size(); _u_MMuPC1AnalysisMQL_166++)
3960 double Xwire = _t_muPC1[_u_MMuPC1AnalysisMQL_166].Xwire;
3961 double Ywire = _t_muPC1[_u_MMuPC1AnalysisMQL_166].Ywire;
3962 double muPC1Time = _t_muPC1[_u_MMuPC1AnalysisMQL_166].muPC1Time;
3964 double _u_MMuPC1AnalysisMQL_303 = time - muPC1Time;
3965 if(_u_MMuPC1AnalysisMQL_303 < 10000)
3970 for(
int _u_MMuPC1AnalysisMQL_305 = _u_MMuPC1AnalysisMQL_166; _u_MMuPC1AnalysisMQL_305 < _t_muPC1.
size(); _u_MMuPC1AnalysisMQL_305++)
3972 double Xwire = _t_muPC1[_u_MMuPC1AnalysisMQL_305].Xwire;
3973 double Ywire = _t_muPC1[_u_MMuPC1AnalysisMQL_305].Ywire;
3974 double muPC1Time = _t_muPC1[_u_MMuPC1AnalysisMQL_305].muPC1Time;
3976 double _u_MMuPC1AnalysisMQL_303 = time - muPC1Time;
3977 if(_u_MMuPC1AnalysisMQL_303 < -10000)
3982 double Xwire_0 = Xwire;
3983 double Ywire_0 = Ywire;
3984 double time_0 =
time;
3985 double muPC1Time_0 = muPC1Time;
3986 double Xwire = Xwire_0;
3987 double Ywire = Ywire_0;
3988 double time = time_0;
3989 double muPC1Time = muPC1Time_0;
3993 if(((time-muPC1Time) >= -160) && ((time-muPC1Time) <= 20))
4002 for(
int _u_MMuPC1AnalysisMQL_306 = 0; _u_MMuPC1AnalysisMQL_306 < _t_muPC1X_sparks.
size(); _u_MMuPC1AnalysisMQL_306++)
4004 double time = _t_muPC1X_sparks[_u_MMuPC1AnalysisMQL_306].time;
4006 double time_0 =
time;
4007 double time = time_0;
4008 int _u_MMuPC1AnalysisMQL_307 = _t_muPC1_sparks_union.
next();
4009 _t_muPC1_sparks_union[_u_MMuPC1AnalysisMQL_307].time = time_0;
4012 for(
int _u_MMuPC1AnalysisMQL_308 = 0; _u_MMuPC1AnalysisMQL_308 < _t_muPC1Y_sparks.
size(); _u_MMuPC1AnalysisMQL_308++)
4014 double time = _t_muPC1Y_sparks[_u_MMuPC1AnalysisMQL_308].time;
4016 double time_0 =
time;
4017 double time = time_0;
4018 int _u_MMuPC1AnalysisMQL_309 = _t_muPC1_sparks_union.
next();
4019 _t_muPC1_sparks_union[_u_MMuPC1AnalysisMQL_309].time = time_0;
4023 bool _u_MMuPC1AnalysisMQL_310[_t_muPC1_sparks_union.
size()];
4024 memset(_u_MMuPC1AnalysisMQL_310, 0,
sizeof(_u_MMuPC1AnalysisMQL_310));
4027 printf(
"Warning: table muPC1_sparks_union is not sorted on time_1 as required for cluster.\n");
4029 for(
int _u_MMuPC1AnalysisMQL_311 = 0; _u_MMuPC1AnalysisMQL_311 < _t_muPC1_sparks_union.
size(); _u_MMuPC1AnalysisMQL_311++)
4031 double time = _t_muPC1_sparks_union[_u_MMuPC1AnalysisMQL_311].time;
4032 if(_u_MMuPC1AnalysisMQL_310[_u_MMuPC1AnalysisMQL_311])
4036 _u_MMuPC1AnalysisMQL_310[_u_MMuPC1AnalysisMQL_311] =
true;
4037 int cluster_size = 1;
4038 int cluster_members[_t_muPC1_sparks_union.
size()];
4039 cluster_members[0] = _u_MMuPC1AnalysisMQL_311;
4041 for(
int _u_MMuPC1AnalysisMQL_312 = _u_MMuPC1AnalysisMQL_311 + 1; _u_MMuPC1AnalysisMQL_312 < _t_muPC1_sparks_union.
size(); _u_MMuPC1AnalysisMQL_312++)
4043 if(_u_MMuPC1AnalysisMQL_310[_u_MMuPC1AnalysisMQL_312])
4047 bool coincidence_seen =
false;
4048 for(
int _u_MMuPC1AnalysisMQL_313 = 0; _u_MMuPC1AnalysisMQL_313 < cluster_size; _u_MMuPC1AnalysisMQL_313++)
4050 int _u_MMuPC1AnalysisMQL_314 = cluster_members[_u_MMuPC1AnalysisMQL_313];
4051 double time_1 = _t_muPC1_sparks_union[_u_MMuPC1AnalysisMQL_314].time;
4052 double time_2 = _t_muPC1_sparks_union[_u_MMuPC1AnalysisMQL_312].time;
4054 double _u_MMuPC1AnalysisMQL_315 = time_1 - time_2;
4055 if(_u_MMuPC1AnalysisMQL_315 < -200)
4059 coincidence_seen =
true;
4060 cluster_members[cluster_size++] = _u_MMuPC1AnalysisMQL_312;
4061 _u_MMuPC1AnalysisMQL_310[_u_MMuPC1AnalysisMQL_312] =
true;
4062 _u_MMuPC1AnalysisMQL_312 = _u_MMuPC1AnalysisMQL_311;
4065 if(!coincidence_seen)
4070 double time_min = _t_muPC1_sparks_union[_u_MMuPC1AnalysisMQL_311].time;
4071 double time_max = time_min;
4072 double time_sum = 0;
4073 for(
int _u_MMuPC1AnalysisMQL_313 = 0; _u_MMuPC1AnalysisMQL_313 < cluster_size; _u_MMuPC1AnalysisMQL_313++)
4075 int _u_MMuPC1AnalysisMQL_314 = cluster_members[_u_MMuPC1AnalysisMQL_313];
4076 time = _t_muPC1_sparks_union[_u_MMuPC1AnalysisMQL_314].time;
4077 time_min =
MIN(time_min, time);
4078 time_max =
MAX(time_max, time);
4079 time_sum = time_sum +
time;
4081 double time_1 = _t_muPC1_sparks_union[_u_MMuPC1AnalysisMQL_311].time;
4082 double time_avg = time_sum/((double)cluster_size);
4083 double time_width = time_max - time_min + 1;
4085 double time = time_min;
4086 int cluster_size_0 = cluster_size;
4087 int cluster_size = cluster_size_0;
4092 if(cluster_size == 2)
4095 double time_0 =
time;
4096 double time = time_0;
4110 _t_muPC1_sparks_union.
toBank(pevent,
"M1SP");
4111 _t_muPC1_sparks_union.
finish();
4112 _t_muPC1X_cluster_contiguous.
finish();
4113 _t_muPC1Y_temporal_auto.
finish();
4114 _t_muPC1X_raw_auto.
finish();
4116 _t_muPC1Y_spark_search.
finish();
4117 _t_muPC1X_cluster_contiguous_sep.
finish();
4118 _t_muPC1X_spatial_auto.
finish();
4119 _t_muPC1X_export.
toBank(pevent,
"MP1X");
4120 _t_muPC1X_export.
finish();
4122 _t_muPC1Y_export.
toBank(pevent,
"MP1Y");
4123 _t_muPC1Y_export.
finish();
4124 _t_muPC1Y_raw_auto.
finish();
4128 _t_muPC1Y_cluster_contiguous_sep.
finish();
4130 _t_muPC1Y_SC_auto.
finish();
4132 _t_muPC1X_spark_search.
finish();
4133 _t_muPC1X_AD_auto.
finish();
4134 _t_muPC1X_sparks.
finish();
4135 _t_muPC1Y_sparks.
finish();
4136 _t_muPC1.
toBank(pevent,
"MPC1");
4139 _t_muPC1_sparks_clustered.
finish();
4144 _t_muPC1_sparks_chamberwide.
finish();
4145 _t_muPC1Y_AD_auto.
finish();
4146 _t_muPC1Y_cluster_contiguous.
finish();
4148 _t_muPC1Y_spatial_auto.
finish();
4150 _t_muPC1X_SC_auto.
finish();