DONTWARN = dangling-else parentheses stringop-truncation
CFLAGS = -O3 -Wall -g -fno-exceptions $(patsubst %,-Wno-%,$(DONTWARN))
LDLIBS = -lm
default: irenafile

UNIT = CHAOS

irenafile:
	make -C ../adam
	ln -sf ../adam/irenafile

%.nE: %.dat
	./irenafile <$< | grep -v ^E >$@

%.EI: %.dat
	./irenafile <$< | ./IRENA.awk 'DROPX();1' >$@

%.Itime: %.EI
	./Itime.awk $< > $@

%.hist: %.EI
	@[ -n "$(HIST)" ] || ( echo please provide HIST= && false )
	./IRENA.awk -v UNIT=$(HIST) '/^EI/ $(CUT){doHIST()}' ${A} $< >$@

%.histb: %.EI IRENA.awk
	@rm -f $@
	@[ -n "$(HIST)" ] || ( echo please provide HIST= && false )
	./IRENA.awk -v UNIT=$(HIST) 'isPHA()$(CUT){doHISTB()}' ${A} $< >$@

%_coinc.hist: %.EI
	@[ -n "$(HIST)" ] || ( echo please provide HIST= && false )
	./IRENA.awk -v UNIT=$(HIST) '/^EI/ $(CUT){doHISTMuon()}' $< >$@

%_xrays.hist: %.EI
	@[ -n "$(HIST)" ] || ( echo please provide HIST= && false )
	./IRENA.awk -v UNIT=$(HIST) 'isXray()$(CUT){doHIST1(C)}' $< >$@

%_muons.hist: %.EI
	@[ -n "$(HIST)" ] || ( echo please provide HIST= && false )
	./IRENA.awk -v UNIT=$(HIST) 'isMuon()$(CUT){doHIST()}' $< >$@

%_musmuons.hist: %.EI
	./IRENA.awk -v UNIT=muMM 'isMuon()$(CUT){doHIST()}' resV=4 $< >$@

%_$(CH).ABPB: %.EI
	@[ -n "$(HIST)" ] || ( echo please provide HIST= && false )
	./IRENA.awk -v UNIT=$(HIST) \
           'isTime(){};Time && isPhase($(CH)) && A>10 {print Time, C, A, B, P, A/Banana(C)}' \
           $< > $@

%_$(CH).SS: %.EI
	./IRENA.awk -v UNIT=$(HIST) 'isSsync($(CH)){printSsync()}' $< >$@

%.PT: %.EI
	./IRENA.awk -i ../pressure.awk -v UNIT=$(HIST) \
           'isTime(){}; \
            Time && isBATE(){print "Pressure", Time, TEMP, P}' \
            $< >$@

%.HK: %.EI IRENA.awk
	./IRENA.awk  -v UNIT=CHAOS 'isHK(){doHK()}' $< >$@

%.HCHAOS: %.EI IRENA.awk
	./IRENA.awk  -v UNIT=CHAOS 'CHAOSHK(){doCHAOSHK()}' $< >$@

%.EH: %.EI
	grep -a '^[EH]' $< > $@

ifneq ($(HIST),)
%.PH: %.EI
	./IRENA.awk '$(BEGIN_HIST)isTime(){$$2=Time;print;next};/^[^E]/' $< >$@
else
%.PH: %.EI
	grep -v ^E $< > $@
endif

data/2019-08-1x-TANOS-tv_%vT.2dhist: data/2019-08-14-TANOS-tv-ambient-1sd.EI data/2019-08-15-TANOS-tv-ambient-2sd.EI
	./IRENA.awk -v UNIT=TANOS_TV 'isTime(){};isPHA($*)==2{doHIST2D(A)}' minV=-10 maxV=100 $^ > $@

.PRECIOUS: %.EI

%.L2: %.EI
	./IRENA.awk '{doL2()}' $< >$@

%.dat: %.b64
	tr -d '\0' < $< | awk '/^=B64 [A-Za-z0-9+/]{43}=$$/{print $$2}' | base64 -d > $@

# µM
%_A1A2.2dhist: %.EI
	./IRENA.awk -v UNIT=muM 'is2D(A1,A2)&&A>-50&&B>-50&&A<500&&B<500{print A,B}' $< |hist.py -s 1/0.838214 -S 1/0.838214 > $@
%_B1B2.2dhist: %.EI
	./IRENA.awk -v UNIT=muM 'is2D(B1,B2)&&A>-50&&B>-50&&A<500&&B<500{print A,B}' $< |hist.py -s 1/0.838214 -S 1/0.838214 > $@
%_AifB.2dhist: %.EI
	./IRENA.awk -v UNIT=muM 'isuMlight(B1,B2) && is2D(A1,A2) && A>-50 && B>-50 && A<500 && B<500 {print A,B}' $< |hist.py -s 1/0.838214 -S 1/0.838214 > $@
%_BifA.2dhist: %.EI
	./IRENA.awk -v UNIT=muM 'isuMlight(A1,A2) && is2D(B1,B2) && A>-50 && B>-50 && A<500 && B<500 {print A,B}' $< |hist.py -s 1/0.838214 -S 1/0.838214 > $@

%_BifA.1dhist: %.EI
	./IRENA.awk -v UNIT=muM 'isuMlight(A1,A2) && isuMlight(B1,B2) {print A+B}' $< |hist.py -s 1/0.838214 > $@
%_AifB.1dhist: %.EI
	./IRENA.awk -v UNIT=muM 'isuMlight(B1,B2) && isuMlight(A1,A2) {print A+B}' $< |hist.py -s 1/0.838214 > $@
%_Light.hist: %.EI
	./IRENA.awk -v UNIT=muM 'isuMlight(A1,A2) && isuMlight(B1,B2) {doHIST()}' $< > $@

CH=AH
PA_CUT=&& B/A<0 && 2000<A && A<2500
%_$(CH).pulse_avg: %.EI
	./IRENA.awk -v UNIT=$(HIST) 'isTime(){};isPhase($(CH)) $(PA_CUT) {print Time, A/Banana(C)}' $< \
      | derive.py -p3 -N1 -s300 -W60 \
      > $@
%_$(CH).dcorr: %.EI %_$(CH).pulse_avg
	./IRENA.awk -v UNIT=$(HIST) 'isTime(){};isPhase($(CH)) $(PA_CUT) {print Time, B/A, A, A*drift_corr(Time, \"$(word 2, $^)\",1)}' $< \
      > $@

%_$(CH)$(CUTN)_p.ihist: %.EH
	./IRENA.awk -v UNIT=$(HIST) 'isPhase($(CH)) && $(CUT) {print B/A}' $< | hist.py -I -s 10000 > $@

data/%.tres_samples: data/%.EI %_tres_results.awk tres.awk
	./tres.awk -i $(word 2, $^) 'T_samples(){}' $< > $@
data/%.tres_delays: data/%.EH %_tres_results.awk tres.awk
	./tres.awk -i $(word 2, $^) 'T_delays(){}' $< > $@

# CHAOS SUMO

SUMO_S1 = is2D(4,5) && 2*A>B && 2*B>A && (E1=A+B)
SUMO_S2 = is2D(6,7) && 2*A>B && 2*B>A && (E2=A+B)
SUMO_CUT= $(SUMO_S1) && $(SUMO_S2)  && Dtime()>10000 && isPHA(17)
SUMO2_CUT= is2D(9, 17) && Dtime()>10000 

%_S12g50_CH.1dhist: %.EI
	./IRENA.awk '$(SUMO_CUT) && E1+E2>50 {print A/16}' $< \
       | hist -s 0.838214 -n -10,250 \
       > $@

%_D3g50_CH.1dhist: %.EI
	./IRENA.awk '$(SUMO2_CUT) && A>50 {print B/16}' $< \
       | hist -s 0.838214 -n -10,250 \
       > $@

%_S12vCH.2dhist: %.EI
	./IRENA.awk '$(SUMO_CUT) {print E1+E2, A/16}' $< \
      | hist -s 0.838214 -n -10,250 -Q \
      > $@

%_D3vCH.2dhist: %.EI
	./IRENA.awk '$(SUMO2_CUT) {print A, B/16}' $< \
      | hist -s 0.838214 -n -10,250 -Q \
      > $@

%_D3vS12.2dhist: %.EI
	./IRENA.awk '$(SUMO_CUT) && $(SUMO2_CUT) {print E1+E2, A}' $< \
      | hist -s 0.838214 -n -10,250 -Q \
      > $@

%_S1vS2.2dhist: %.EI
	./IRENA.awk '$(SUMO_S1) && $(SUMO_S2) && Dtime()>10000 {print E1, E2}' $< \
      | hist -s 0.838214 -n -10,250 -Q \
      > $@

%_S1g50_D3g50_CH.1dhist: %.EI
	./IRENA.awk '$(SUMO_S1) && E1>50 && $(SUMO2_CUT) && A>50 {print B/16}' $< \
       | hist -s 0.838214 -n -10,250 \
       > $@

%_S2g50_D3g50_CH.1dhist: %.EI
	./IRENA.awk '$(SUMO_S2) && E2>50 && $(SUMO2_CUT) && A>50 {print B/16}' $< \
       | hist -s 0.838214 -n -10,250 \
       > $@