#!/usr/bin/awk -f BEGIN {i=0; sum=0.; integral=0.; binwidth=0.; underflow=0.; overflow=0. } { if (i==0) { underflow=$5; binwidthh=(0.5/($2-2.))*($4-$3); } if (i>0) { sum+=$2*($1+binwidthh)*binwidthh*2.; integral+=$2*binwidthh*2.; } ++i; } END { printf "mean=" sum " integral=" integral+(underflow*binwidthh*2.) " 1/i=" 1./(integral+(underflow*binwidthh*2.))" binwidth=" binwidthh*2. "\n" > "/dev/stderr" ; print sum/(integral+(underflow*binwidthh*2.)); } # FILES=`find . -name jet_1_1_pt_1.dat`; for I in $FILES; do echo $I; ../calc_mean < $I; done