CASP7 ?= /projects/compbio/experiments/protein-predict/casp7
TARGET_FILE ?= ${CASP7}/target_list.txt

# TO DO:
#
# Several of the targets here assume that all targets to evaluate are
# 	full-length sequences.
# 	Evaluation by domains will require rewriting some of this
#	Makefile (and probably some of the scripts).

%.all:	%.distribute-costfcn %.evaluate  %.summarize
	date
	echo $*.all done

ALPHA_POWERS ?= 0 2 3 4 5 10 20 50

%.summarize: 
	date
	make -k $(foreach c,gdt real_cost, $(foreach a,${ALPHA_POWERS}, $*-${c}-tau_${a}.rdb $*-${c}-tau_${a}.avg))
	make -k $(foreach c,gdt real_cost,  $*-${c}-btr.rdb $*-${c}-btr.avg)
	date
	echo $*.summarize done.

%.evaluate:   	
	date
	para-trickle-make -manyids -se2log -no2letter -modelsdir ${CASP7} \
		-makefile ./Makefile -target decoys/predictions.evaluate.$*.rdb \
		-wait -rmvdir \
		< ${TARGET_FILE}
	date
	echo $*.evaluate done

%.distribute-costfcn:	%.costfcn
	$(foreach t,$(shell cat ${TARGET_FILE}),\
		sed -e s/TARGET/${t}/g < $^ > ${CASP7}/${t}/$^ ;)

merged:
	make -k $(foreach c,gdt real_cost, $(foreach a,${ALPHA_POWERS}, ${c}-tau_${a}-merged.rdb ${c}-tau_${a}-merged.avg))
	make -k $(foreach c,gdt real_cost,  ${c}-btr-merged.rdb ${c}-btr-merged.avg)

define TAU_MAKE
%-gdt-tau_${1}.rdb:
	time ./cfneval.pl -s decoys/predictions.evaluate.$$*.rdb -f0 --alpha ${1} > $$@

%-real_cost-tau_$(1).rdb:
	time ./cfneval.pl -s decoys/predictions.evaluate.$$*.rdb -f20 --alpha ${1} -q real_cost > $$@
%-tau_${1}.avg: %-tau_${1}.rdb
	summ -m < $$^  \
	| sort -gr +7 \
	> $$@
endef

$(foreach a,${ALPHA_POWERS},$(eval $(call TAU_MAKE,${a})))

%-gdt-btr.rdb:
	time ./cfneval.pl -s decoys/predictions.evaluate.$*.rdb -f0 -r btr  > $@

%-real_cost-btr.rdb:
	time ./cfneval.pl -s decoys/predictions.evaluate.$*.rdb -f20 -q real_cost -r btr > $@

%-btr.avg: %-btr.rdb
	summ -m < $^  \
	| sort -g +7 \
	> $@


%-merged.rdb: builtins-%.rdb pred-%.rdb 	\
		ehl2+sheets-%.rdb undertaker-sheets-%.rdb \
		dssp-ehl2-%.rdb contact449a_45-%.rdb \
		constraint-optimized95-0.3-%.rdb constraint-fraction-local100-%.rdb
	cat  builtins-$*.rdb \
	| jointbl target pred-$*.rdb \
	> tmp1.rdb
	$(foreach x,ehl2+sheets undertaker-sheets dssp-ehl2 contact449a_45 constraint-optimized95-0.3 constraint-fraction-local100 , \
		column target -c constraints $(x) < $(x)-$*.rdb > tmp2.rdb; \
		jointbl target tmp2.rdb < tmp1.rdb > tmp3.rdb; \
		mv -f tmp3.rdb tmp1.rdb; )
	mv -f tmp1.rdb $@
	rm tmp2.rdb


%-tau-merged.avg: %-tau-merged.rdb
	summ -m < $^  \
	| sort -gr +7 \
	> $@
%-tau_0-merged.avg: %-tau_0-merged.rdb
	summ -m < $^  \
	| sort -gr +7 \
	> $@
%-tau_1-merged.avg: %-tau_1-merged.rdb
	summ -m < $^  \
	| sort -gr +7 \
	> $@
%-tau_2-merged.avg: %-tau_2-merged.rdb
	summ -m < $^  \
	| sort -gr +7 \
	> $@

%-btr-merged.avg: %-btr-merged.rdb
	summ -m < $^  \
	| sort -g +7 \
	> $@

# Sat Nov 17 15:03:06 PST 2007  John Archie
#  I've removed a couple of the TODO items as they have been done.