#!/usr/bin/awk -f

@include "RPi.awk"

BEGIN {
    pmin = 900
    pmax = 1100
    resp = 1

    minp = 10000
    maxp = 0

    muonthr = 4
    maxt = 60
}

isH(){}

isP() {
    t = Time - pTime;
    pTime = Time
    if (p && t>0) {
            if (t>maxt) t = maxt;
            Np[p] += t;
    }
    if (Pressure < pmax && Pressure > pmin) {
            p = int(Pressure/resp)
            if (p<minp) minp = p
            if (p>maxp) maxp = p
    } else {
        p = 0
    }
}

p && isE() {
    t = Time - pTime;
    if (t<0 || t>maxt) next;
    if (A[0] < muonthr && A[3] < muonthr) next;
    if (A[1] >= muonthr) R1[p] += 1
    if (A[2] >= muonthr) R2[p] += 1
}

END {
    for (p=minp; p<=maxp; p++) {
            pr = (p+0.5)*resp
            N = Np[p]+0
            if (N) printf "%.1f %g %g %g %g\n", pr, R1[p]/N, sqrt(R1[p]+1)/N, R2[p]/N, sqrt(R2[p]+1)/N
            else printf "%.1f 0 0 0 0\n", pr
    }
}