--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/makeSpecE2.sh	Fri Apr 25 19:33:20 2025 +0000
@@ -0,0 +1,19 @@
+#!/bin/bash
+#prints normalized spectrum dN(E)/dE * E^2
+inidir=$PWD
+
+if [ $# -lt 1 ]
+then
+    echo "usage ./makeSpecE2.sh file"
+    exit 1
+fi
+
+if [ -f $1 ]
+then
+	ofile=$1.spec
+	cat $1 | awk '!/^($|[[:space:]]*#)/{ printf("%.1f %g\n",log($1)/log(10)+0.05,$2)}' | sort -g -k 1 | awk '{ if($1!=PREV&&NR>1) printf("%g %g %d\n",PREV,SUM,SUMN); if($1!=PREV||NR==1) {SUM=0; SUMN=0;}  SUM=SUM+$2; SUMN=SUMN+1; PREV=$1; } END {printf("%g %g %d\n",PREV,SUM,SUMN)}' | awk 'NF==3 { E=10^($1-0.05); printf("%g %g %d\n", E, $2*E/(10^0.05-10^(-0.05)),$3)}'>$ofile
+else
+	echo $1 not found
+	exit 1
+fi
+