Motivations
Details
In dataMonitoring/Common Common-06-10-00 there's a first attempt at a piece of code to be run in the pipeline.
It takes the digi and merit trending root files for a given run and the basic usage is:
Usage: pBadTimeIntervalLogger.py [options] digitrend merittrend Options: -h, --help show this help message and exit -r R, --root-output=R path to the output root file -i, --interactive run in interactive mode (show the plots) -p P, --png-output=P path to the output png file -x X, --xml-output=X path to the output xml file -f F, --flare-threshold=F threshold on the normalized ACD hit rate for the flare intervals -c C, --count-threshold=C threshold on the normalized ACD tile count for the bad time intervals -t T, --min-time-padding=T minimum time padding (in s) for fitting the normalized transient rate outside any flare interval -M M, --min-int-loss=M minimum integral loss (in s) for defining a bad time interval -m M, --min-diff-loss=M minimum differential loss for defining a bad time interval -s S, --int-loss-start=S integrated loss (in s) defining the start of a bad time interval -e E, --frac-loss-end=E fraction of the total integrated loss defining the end of a bad time interval -n N, --num-plateau-points=N number of points for the evaluation of the total integrated loss Usage: pBadTimeIntervalLogger.py [options] digitrend merittrend
At first order you can forget about most of the parameters and do something:
./pBadTimeIntervalLogger.py /data/work/datamon/solartrend/runs/r0319424872_digitrend.root /data/work/datamon/solartrend/runs/r0319424872_merittrend.root -r test.root -p test.png -x test.xml
The png output looks like this:
test.png
And the xml output looks like:
<?xml version="1.0" encoding="ISO-8859-15"?> <!-- Generated by pXmlWriter on Thu Apr 7 16:25:26 2011. --> <solarFlareSummary> <flareIntervals> <interval duration="1995.0" source="NormAcdTileRate63" start_met="319427310.0" end_met="319429305.0"/> </flareIntervals> <badIntervals> <interval duration="1410.0" source="NormAcdTileCount" start_met="319427325.0" end_met="319428735.0"/> <interval duration="1575.0" source="NormTransientRate" integral_loss="1008.7" start_met="319427325.0" end_met="319428900.0"/> </badIntervals> </solarFlareSummary>