#!/bin/bash

# The LMA file is for a single sector, so filter the twiss file to cover the same interval
sddsprocess lma.twi lma.twi1 -filter=col,s,0,27.61

for slice in run-*.sa00 ; do 
    root=${slice/.sa00/}
    sddssplit $slice -digits=3 -rootname=$slice- -exten=
    for sample in ${slice}-??? ; do 
        index=${sample/$slice-/}
        echo processing $sample $index
        touschekLifetime -twiss=lma.twi1 -aperture=lma.mmap -sliceAnalysis=$sample ${root}.tl$index -ignoreMismatch
        \rm $sample
    done
    echo combining $root
    sddscombine ${root}.tl??? -collapse  -pipe=out \
        | sddsprocess -pipe=in ${root}.tl -process=tLifetime,ave,%sAve -process=tLifetime,stand,%sStDev
    \rm ${root}.tl???
done

for result in *.tl ; do 
    sddsprintout -param=tLifetime* $result
done

