#! /usr/bin/gawk -f # read the output of multiplicities.py and produce the good or # bad events # ./tanos-cut.awk BAD=1 data/tanos-jr-flight.EIDT BEGIN { BAD = 0 } /^EIDT/ { D2 = 65536 D3 = 65536 if (NF>=18) D2 = $18 if (NF>=19) D3 = $19 good = D2>100 && D3>1000 && (D2<1610 || D2>=1640) if (!BAD && good || BAD && !good) print }