PARSER = ahepamfile
include ../irena/libirena.make
.PRECIOUS: %.dat %.AHA

%.AHA: %.dat ahepamfile
	./ahepamfile < $< > $@

%.doS: %.AHA
	./AHEPAM.awk 'doS(){}' $< > $@

%.AES: %.dat ahepamfile
	./ahepamfile < $< | grep '^[ES]' | grep -v 'X' | sort -nk2,2 -k1 > $@

%.CH_$(ch): %.dat AHEPAM.awk Makefile
	./ahepamfile < $< \
      | ./AHEPAM.awk 'isS($($(ch))){if(N==1)print Nix;print ELINE, S, SS, TT, T1}' \
      > $@

%.2D_$(ch1)$(ch2): %.dat AHEPAM.awk Makefile
	./ahepamfile < $< \
      | ./AHEPAM.awk 'isS2($($(ch1)),$($(ch2))){print A1, A2, B1, B2, P1, P2, Ba1, Ba2, E1, E2, S1, S2, SS1, SS2, TT1, TT2, T1, T2}' \
      > $@

%.HDORN: %.AHA HDORN.awk
	./HDORN.awk 'isHDORN(){doHDORN()}' $< > $@

SIGMA=600

%.HD4BIAS: %.HDORN
	derive.py -x 3 -y 4,6-18 -s $(SIGMA) -F %5.2f -N 1 $< > $@

%.HD4PEEK: %.HDORN
	derive.py -x 3 -y 4-16 -s $(SIGMA) -F %.2f -N 1 $< > $@

resV=1
EE=E
%.hist: %.AHA
	./AHEPAM.awk 'isEE(){doHIST(Ch,$(EE))}' resV=$(resV) $< > $@

%.SETH: %.AHA
	./AHEPAM.awk 'isSETH(){doSETHBGO()}' $< > $@

CUT=P
CUT_P=&&P>-0.398&&P<=-0.105
CUT_nP=&&(P<=-0.398||P>-0.105)

%_cut$(CUT).hist: %.AHA
	./AHEPAM.awk 'isE($$3,$$4)$(CUT_$(CUT)){doHIST(Ch,EE)}' resV=$(resV) $< > $@



A=0,2
B=0,1
C=0,0
D=1,2
E=1,1
F=1,0
G=2,2
H=2,1
I=2,0
J=3,2
K=3,1
L=3,0

ch=F
ch1=F
ch2=E

%_$(ch1)$(ch2).2dhist: %.dat ahepamfile
	./ahepamfile < $< \
      | ./AHEPAM.awk 'is2D($($(ch1)),$($(ch2)), 3) {do2DHIST(E1,E2)}' \
           minV=-10 maxV=250 resV=$(resV) \
      > $@

%_$(ch).AHAE: %.dat ahepamfile
	./ahepamfile < $< | ./AHEPAM.awk 'isE($($(ch))) && E>170 && E<250' > $@