hj<vdZddlmZddlZddlZddlmZddlZdZ dZ dZ d dZ e d k(r ee y) aScore AI load cycles against the furnace log by PHASE, not by event matching. The furnace is close to a metronome. In the pilot hour its charge interval was 212-229 s, and load cycles arrive at a similar rate, so a +-90 s matching window gives ~0.77 recall to *any* detector firing every ~210 s -- including one firing at random. Measured on the pilot: recall 0.938 at the true offset vs 0.773 mean at |offset| >= 180 s, and several wrong offsets scored 1.000. Recall/precision therefore say almost nothing here, and the old 0.84 headline should not be quoted. What does carry information is PHASE. If the detector really sees loading, each AI cycle should land at a repeatable point inside the furnace's charge-to-charge interval. If it is noise, phases spread uniformly. Because the interval is not constant across the day (breaks stretch it), phase is normalised per interval: phase = (t_ai - charge_prev) / (charge_next - charge_prev) in [0, 1) then tested for concentration with the Rayleigh test. `--control` re-runs the same test on random timestamps to confirm the test is calibrated -- it should reject at about the nominal 5%. Pilot result for reference: R=0.754, p<0.001, mean phase +41 s after charge, circular SD 26 s; the random control passed 5.5% of the time. ) annotationsN)Pathctjt|t}gg}}|D]s}tj||dz }|dks|dzt |k\r3||||dz}}||krF|j ||z ||z z |j ||z utj|tj|fS)zGNormalised position of each event inside its enclosing charge interval.r)npasarraysortedfloat searchsortedlenappend) eventschargeschoutspanstiabs W/home/brad/proj/sonice-ai-sentinel-edge/cursor-glass-factory/scripts/a1_phase_metric.pyphase_ofr$s F7OU +BRC  OOB "Q & q5AESW$ !ubQi1 6  AEa!e$% QU ::c?BJJu- --cdt|dkrddtdfSdtjz|z}tj|j tj |j }}ttj||}t|}||z|z}ttj| dd|z||zz d|zz zz}tt|dd}ttjtj||dzd z }|||fS) z=Return (R, p, mean_phase). R=0 uniform, R=1 perfectly locked.gg?nanrihgv@) r r rpicosmeansinhypotexpminmaxdegreesarctan2) phasesangCSRnZpmus rrayleighr24s 6{QCu%% bee)f C 66#;   rvvc{//1qA bhhq!nA F A A A bffaRjAQQ1q5 99:;A C3KA  2::a+,s2e; r)r rsqrtlogr&r)r-s r circ_sd_fracr7Cs7 bffSD\223q255yA BBrc V tj}|jddd|jddd|jdd|jd td d |jd dgd|jdtdd |jdtd|jdtd|jdd|j }t j t|jjdd}t j t|jjd}|d|j}|d}|Dcgc]}||j|jz }}|s td|j |jn t!|}|j" |j"n t%|} |D cgc]} || cxkr| ksnn| }} |s td |D cgc]} || cxkr| ksnn| } } t'||\} }t)| \}}}t+|rt-t/j0|n t-d!}t.j2j5d }d }t7|j8D]D}|j;|| t+|}t'||\}}t)|\}}}||d"kz }F|t%|j8d#z }id$|jd%|jd&t+|d't+| d(| r!t=t+|t+| z d)ndd*| r0t=d+t+|t+| z zt+| z d#ndd,tt+| d-t=|d#d.t=|d/d0t=|d1d2||k(rt=||zd#ndd3||k(rt=t?||zd#ndd4t=|d)d5tA|d"kd6| Dcgc]}t=t-|d/c}d7tC| |Dcgc]\}}t=t-||zd#c}}d8|d"krd9nd:}|jDr?t|jDjGt jH|d;d<=d>zdtKt jH|d;d<=y cc}wcc} wcc} wcc}wcc}}w)?Nz--cyclesTzload_cycles_v2.json)requiredhelpz --furnacezfurnace_20260617.jsonz--clip)r9z --video-t0rz(clip-local offset of the analysed window)typedefaultr:z--keyt_start)r=t_peakt_end)r<choicesz --win-startz'clip-local start of the analysed windowz --win-end)r;r<z --controlz--out)r<zutf-8)encodingcyclesper_clipcharge_video_tz no AI cyclesz(no AI cycles inside the requested windowrg?rclipanchor n_ai_cyclesn_furnace_charges_in_window count_ratiorcount_error_pctdn_phasedmedian_interval_secphase_Rr rayleigh_pmean_offset_sec_after_chargecircular_sd_seccontrol_false_positive_ratelockedr) offsets_secinterpretationzAI cycles are phase-locked to the machine-logged furnace charges; the random-timestamp control rejects at the nominal rate, so the test is calibrated.z`Phase is not distinguishable from uniform -- no evidence the detector tracks the furnace rhythm.Fr) ensure_asciiindent )&argparseArgumentParser add_argumentint parse_argsjsonloadsrrD read_textfurnacerGkeyvideo_t0 SystemExit win_startr%win_endr&rr2r r rmedianrandom default_rngrangecontroluniformroundr7boolzipr write_textdumpsprint)apargscycfurrGrcevlohirgin_winphrr-r0r1periodrngnpass_fakefppcctrlvsreps rmainrGs  "BOOJ4IOJOOK$5LOMOOHtO,OOLsA>53r7B1s2wB (2V!LaP`d   CBL  uVQ/  5A;  eAqk  '6AQb6k1(=W[  6AQ5a6!91=W[  &uT1~  $q4x. " r2!5q1%2# $ 3r5>J41aeAElA.J% & 4x "21 C8 xx TXX!!$**SuQ"ORV"Vah!i $**SuQ 78 u 4 )2@3Js*9#T5TT T1TT 3"T% __main__)returnr_)__doc__ __future__rr\rapathlibrnumpyrrr2r7r__name__rgrrrsN2#  . CL ^ z TV r