| 1 | #!/bin/bash
|
|---|
| 2 | #
|
|---|
| 3 |
|
|---|
| 4 | source `dirname $0`/../Sourcefile.sh
|
|---|
| 5 | printprocesslog "INFO starting $0"
|
|---|
| 6 |
|
|---|
| 7 | echo "Sent data"
|
|---|
| 8 |
|
|---|
| 9 | start=$1
|
|---|
| 10 | stop=$2
|
|---|
| 11 | source=$3
|
|---|
| 12 | bin=$4
|
|---|
| 13 | email=$5
|
|---|
| 14 | table=$6
|
|---|
| 15 | time=$7
|
|---|
| 16 | expert=$8
|
|---|
| 17 |
|
|---|
| 18 | query="SELECT fSourceName FROM Source WHERE fSourceKey="$source
|
|---|
| 19 | sourcename=`sendquery`
|
|---|
| 20 |
|
|---|
| 21 | function printinfo()
|
|---|
| 22 | {
|
|---|
| 23 | echo "# "
|
|---|
| 24 | echo "# Please cite the FACT design paper and the QLA webpage when using these data."
|
|---|
| 25 | echo "# FACT design paper: H. Anderhub et al. JINST 8 (2013) P6008 "
|
|---|
| 26 | echo "# http://iopscience.iop.org/1748-0221/8/06/P06008 "
|
|---|
| 27 | echo "# QLA webpage: http://www.fact-project.org/monitoring"
|
|---|
| 28 | echo "# If you intent to use the data, please let us know for reference. "
|
|---|
| 29 | echo "# "
|
|---|
| 30 | echo "# Remarks:"
|
|---|
| 31 | echo "# * These are the results of a fast quick look analyis "
|
|---|
| 32 | echo "# on site, i.e. they are preliminary. "
|
|---|
| 33 | echo "# * The quick look analysis includes all data, "
|
|---|
| 34 | echo "# i.e. no data selection done."
|
|---|
| 35 | echo "# * The given values are not fluxes but excess rates "
|
|---|
| 36 | echo "# (number of excess events per effective ontime), "
|
|---|
| 37 | echo "# i.e. there is a dependence on trigger threshold and "
|
|---|
| 38 | echo "# zenith distance of the observation (with the current "
|
|---|
| 39 | echo "# analysis for zenith distance > 40 degree and trigger "
|
|---|
| 40 | echo "# threshold > 500 DAC counts)."
|
|---|
| 41 | echo "# * Nights with less than 20 minutes of data are neglected "
|
|---|
| 42 | echo "# for nightly binning. "
|
|---|
| 43 | echo "# * The QLA results are not reprocessed when a new software "
|
|---|
| 44 | echo "# version is introduced. "
|
|---|
| 45 | echo "# * In case, you need further details about the data or a"
|
|---|
| 46 | echo "# different binning, please do not hesitate to contact us."
|
|---|
| 47 | echo "# * The QLA contains all data since 12.12.2012. "
|
|---|
| 48 | echo "# For older data, please contact us. "
|
|---|
| 49 | echo "# "
|
|---|
| 50 | echo "# Contact: Daniela Dorner dorner@astro.uni-wuerzburg.de "
|
|---|
| 51 | echo "# "
|
|---|
| 52 | echo "# This file was created at "`date`
|
|---|
| 53 | echo "# Source: "$sourcename
|
|---|
| 54 | echo "# Timerange: "$start" - "$stop
|
|---|
| 55 | if [ $bin -lt 0 ]
|
|---|
| 56 | then
|
|---|
| 57 | echo "# Binning: nightly "
|
|---|
| 58 | else
|
|---|
| 59 | echo "# Binning: "$bin" min "
|
|---|
| 60 | fi
|
|---|
| 61 | echo "# "
|
|---|
| 62 | #echo "# start("$time") stop("$time") excess-rate(evts/h) error-excess-rate(evts/h) ontime(min) significance[if Nsig>10 and Nbg>10 (or Nexc>10)] "
|
|---|
| 63 | echo "# start("$time") stop("$time") excess-rate(evts/h) error-excess-rate(evts/h) "
|
|---|
| 64 | getdbsetup
|
|---|
| 65 | printprocesslog "DEBUG sendquery QUERY: "$query
|
|---|
| 66 | mysql -s -u $us --password=$pw --host=$ho $db -e " $query "
|
|---|
| 67 | }
|
|---|
| 68 |
|
|---|
| 69 | ontime1=" TIME_TO_SEC(TIMEDIFF(fRunStop,fRunStart))*fEffectiveOn "
|
|---|
| 70 | ontime2=" fOnTimeAfterCuts "
|
|---|
| 71 | ontimeif=" IF(ISNULL(fEffectiveOn), fOnTimeAfterCuts, TIME_TO_SEC(TIMEDIFF(fRunStop,fRunStart))*fEffectiveOn) "
|
|---|
| 72 | from=" FROM RunInfo LEFT JOIN "$table" USING (fNight, fRunID) "
|
|---|
| 73 | # time range and source
|
|---|
| 74 | where=" WHERE fSourceKey="$source" AND fNight BETWEEN "$start" AND "$stop
|
|---|
| 75 | # some sanity checks
|
|---|
| 76 | where=$where" AND fRunTypeKey=1 AND NOT ISNULL(fNumSigEvts) AND NOT ISNULL(fNumBgEvts) AND NOT fRunStart='0000-00-00 00:00:00' AND NOT fRunStop='0000-00-00 00:00:00' "
|
|---|
| 77 |
|
|---|
| 78 | case $time in
|
|---|
| 79 | mjd) timeselect="Mjd(Min(fRunStart)) AS start, Mjd(MAX(fRunStop)) AS stop, "
|
|---|
| 80 | timeselect2=" Mjd(MIN(o.start)) AS start, Mjd(MAX(o.stop)) AS stop, "
|
|---|
| 81 | break;;
|
|---|
| 82 | unix) timeselect="Unix_timestamp(CONVERT_TZ(Min(fRunStart), '+00:00', 'SYSTEM')) AS start, Unix_timestamp(CONVERT_TZ(Max(fRunStop), '+00:00', 'SYSTEM')) AS stop, "
|
|---|
| 83 | timeselect2=" Unix_timestamp(CONVERT_TZ(MIN(o.start), '+00:00', 'SYSTEM')) AS start, Unix_timestamp(CONVERT_TZ(MAX(o.stop), '+00:00', 'SYSTEM')) AS stop, "
|
|---|
| 84 | break;;
|
|---|
| 85 | *) timeselect="MIN(fRunStart) AS start, MAX(fRunStop) AS stop, "
|
|---|
| 86 | timeselect2=" MIN(o.start) AS start, MAX(o.stop) as stop, "
|
|---|
| 87 | break;;
|
|---|
| 88 | esac
|
|---|
| 89 |
|
|---|
| 90 | # todo:
|
|---|
| 91 | # implement login
|
|---|
| 92 | # implement for internal use
|
|---|
| 93 | # table
|
|---|
| 94 | # more binnings (insert box for minutes? period, 10d?, weekly?)
|
|---|
| 95 | # zd limit
|
|---|
| 96 | # threshold limit
|
|---|
| 97 | # datacheck
|
|---|
| 98 |
|
|---|
| 99 | if [ $bin -eq -1 ]
|
|---|
| 100 | then
|
|---|
| 101 | #query="SELECT Min(fRunStart), Max(fRunStop), Sum(fNumExcEvts)/Sum(fOnTimeAfterCuts)*3600, if ((Sum(fNumSigEvts) + 0.2*0.2*Sum(fNumBgEvts))<0, 0, Sqrt(Sum(fNumSigEvts) + 0.2*0.2*Sum(fNumBgEvts))) FROM RunInfo LEFT JOIN "$table" USING (fNight, fRunID) WHERE fSourceKey="$source" AND fNight BETWEEN "$start" AND "$stop" GROUP BY fNight HAVING SUM(fOnTimeAfterCuts)*60>20"
|
|---|
| 102 | query="SELECT "$timeselect" SUM(fNumExcEvts)/SUM("$ontimeif")*3600 AS excrate, "
|
|---|
| 103 | query=$query" ExcErr(Sum(fNumSigEvts), SUM(fNumBgEvts))/SUM("$ontimeif")*3600 AS excerr "
|
|---|
| 104 | #query=$query" IF (SUM(fNumSigEvts)>10 AND SUM(fNumBgEvts)>10, LiMa(SUM(fNumSigEvts), SUM(fNumBgEvts)), IF (SUM(fNumSigEvts)-SUM(fNumBgEvts)>10, concat('(', LiMa(SUM(fNumSigEvts), SUM(fNumBgEvts)), ')'), '-')) AS significance "
|
|---|
| 105 | #query=$query" Avg(fThresholdMedian) as th "
|
|---|
| 106 | query=$query$from$where" GROUP BY fNight HAVING SUM("$ontimeif")/60>20 ORDER BY fNight "
|
|---|
| 107 | else
|
|---|
| 108 | # query to get excess rate and significance from DB
|
|---|
| 109 | query="SELECT "$timeselect2" ROUND((SUM(o.sigevts)-SUM(o.bgevts))/SUM(o.ot)*3600, 2) AS excrate, "
|
|---|
| 110 | query=$query" ROUND(ExcErr(SUM(o.sigevts),SUM(o.bgevts))/SUM(o.ot)*3600,2) AS err "
|
|---|
| 111 | #query=$query" ROUND(SUM(o.ot)/60.,1) AS 'ontime [min]' "
|
|---|
| 112 | #query=$query" IF (SUM(o.sigevts)>10 AND SUM(o.bgevts)>10, LiMa(SUM(o.sigevts), SUM(o.bgevts)), IF (SUM(o.sigevts)-SUM(o.bgevts)>10, CONCAT('(', LiMa(SUM(o.sigevts), SUM(o.bgevts)), ')'), '-')) AS significance "
|
|---|
| 113 | query=$query" FROM (SELECT fNight, @ot:="$ontimeif" AS ot, fRunStart AS start, fRunStop AS stop, fNumSigEvts AS sigevts, fNumBgEvts AS bgevts, "
|
|---|
| 114 | query=$query" IF (@night=fNight AND FLOOR((@os+@ot)/20./60.)<1, @bl, @bl := @bl + 1) AS block, "
|
|---|
| 115 | query=$query" IF (@night=fNight AND FLOOR((@os+@ot)/20./60.)<1, @os:=@os + @ot, @os := @ot) AS os, @night :=fNight AS night "
|
|---|
| 116 | query=$query$from" CROSS JOIN (SELECT @night :=0, @ot :=0, @os :=0, @bl:=0) PARAMS "
|
|---|
| 117 | query=$query$where" ORDER BY fRunStart) o GROUP BY block ORDER BY start"
|
|---|
| 118 | fi
|
|---|
| 119 | #echo $query
|
|---|
| 120 |
|
|---|
| 121 | #printinfo | mail -a "From: dorner@astro.uni-wuerzburg.de" -s 'test downloading data ' -b dorner@astro.uni-wuerzburg.de $email
|
|---|
| 122 | printinfo | mail -s 'FACT data (test download)' -c dorner@astro.uni-wuerzburg.de -r dorner@astro.uni-wuerzburg.de $email
|
|---|
| 123 | #printinfo
|
|---|
| 124 |
|
|---|
| 125 | finish
|
|---|
| 126 |
|
|---|
| 127 | # new query:
|
|---|
| 128 | # select o.night, Min(o.start), Max(o.stop), Sum(o.sigevts), Sum(o.bgevts), Round(Sum(o.ot)/60.,1) FROM (Select fNight, fRunID, fOnTimeAfterCuts as ot, fRunStart as start, fRunStop as stop, fNumSigEvts as sigevts, fNumBgEvts as bgevts, if (@night=fNight AND floor((@os+fOnTimeAfterCuts)/20./60.)<1, @bl, @bl := @bl + 1) as block, if (@night=fNight AND floor((@os+fOnTimeAfterCuts)/20./60.)<1, @os:=@os + fOnTimeAfterCuts, @os := fOnTimeAfterCuts) as os, @os/20./60. as a, floor(@os/20./60.) as b, @night :=fNight as night from RunInfo left join AnalysisResultsRunLP using(fNight,fRunID) CROSS JOIN (SELECT @night :=0, @os :=0, @bl:=0) PARAMS where fNight between 20130414 and 20130415 and fRunTypeKey=1 and fSourceKey=1) o group by block;
|
|---|
| 129 |
|
|---|
| 130 | # input
|
|---|
| 131 | # - binning
|
|---|
| 132 | # - source
|
|---|
| 133 | # - email address
|
|---|
| 134 | # - night range
|
|---|
| 135 |
|
|---|
| 136 | # binning
|
|---|
| 137 | bin=20
|
|---|
| 138 |
|
|---|
| 139 | # get date
|
|---|
| 140 | # add funtion
|
|---|
| 141 | night=20160102
|
|---|
| 142 |
|
|---|
| 143 | # get sources for current night from DB (RunInfo)
|
|---|
| 144 | sourcekey=1
|
|---|
| 145 |
|
|---|
| 146 | # get limits from DB
|
|---|
| 147 | siglimit=3.0
|
|---|
| 148 | exclimit=20
|
|---|
| 149 |
|
|---|
| 150 | # query to get excess rate and significance from DB
|
|---|
| 151 | query="SELECT o.* FROM"
|
|---|
| 152 | query=$query" (SELECT fRunID, IF (@os+fOnTimeAfterCuts/60.>"$bin", @sig := fNumSigEvts, @sig := @sig + fNumSigEvts) AS sigevts, "
|
|---|
| 153 | query=$query" IF (@os+fOnTimeAfterCuts/60.>"$bin", @bg := fNumBgEvts, @bg := @bg + fNumBgEvts) AS bgevts, "
|
|---|
| 154 | query=$query" IF (@os+fOnTimeAfterCuts/60.>"$bin", @exc := fNumExcEvts, @exc := @exc + fNumExcEvts) AS excevts, "
|
|---|
| 155 | query=$query" IF (@os+fOnTimeAfterCuts/60.>"$bin", @os := fOnTimeAfterCuts/60., @os := @os + fOnTimeAfterCuts/60.) AS ontime, "
|
|---|
| 156 | query=$query" @exc/@os*60 AS excrate, "
|
|---|
| 157 | query=$query" IF (@sig>0 AND @bg*5>0 AND (@sig*Log(@sig/(@sig+@bg*5)*1.2/0.2)+@bg*5*Log(@bg*5/(@sig+@bg*5)*1.2))> 0, Sqrt((@sig*Log(@sig/(@sig+@bg*5)*1.2/0.2)+@bg*5*Log(@bg*5/(@sig+@bg*5)*1.2))*2), 0) AS sig "
|
|---|
| 158 | query=$query" FROM AnalysisResultsRunLP O LEFT JOIN RunInfo USING (fNight, fRunID) "
|
|---|
| 159 | query=$query" CROSS JOIN (SELECT @os :=0, @sig :=0, @bg := 0, @exc :=0, @rate :=0) PARAMS "
|
|---|
| 160 | query=$query" WHERE fNight="$night" AND fSourceKey="$sourcekey") o HAVING ontime > "$bin"*0.9"
|
|---|
| 161 | results=( `sendquery` )
|
|---|
| 162 | #echo ${results[@]}
|
|---|
| 163 |
|
|---|
| 164 | i=0
|
|---|
| 165 | while [ 0 -lt 1 ]
|
|---|
| 166 | do
|
|---|
| 167 | num=`echo "$i * 7" | bc -l`
|
|---|
| 168 | if [ "${results[$num]}" = "" ]
|
|---|
| 169 | then
|
|---|
| 170 | break
|
|---|
| 171 | fi
|
|---|
| 172 | runid=${results[$num]}
|
|---|
| 173 | sig=${results[$num+6]}
|
|---|
| 174 | exc=${results[$num+5]}
|
|---|
| 175 |
|
|---|
| 176 | trigger=`echo " $exc > $exclimit && $sig > $siglimit " | bc -l`
|
|---|
| 177 |
|
|---|
| 178 | # echo $i" "$num" "$exc" "$sig
|
|---|
| 179 | if [ $trigger -eq 1 ]
|
|---|
| 180 | then
|
|---|
| 181 | # do whatever to be done to trigger
|
|---|
| 182 | # - send email/sms / call
|
|---|
| 183 | # - create amon file
|
|---|
| 184 | # - prepare email for alert
|
|---|
| 185 | echo $runid" "$exc" "$sig" -> "$trigger
|
|---|
| 186 | triggerfile="/home/fact/flare.alerts/"$night"_"$runid"_"$sourcekey".trigger"
|
|---|
| 187 | touch $triggerfile
|
|---|
| 188 | fi
|
|---|
| 189 | # counter
|
|---|
| 190 | i=`echo $i +1 | bc -l`
|
|---|
| 191 | done
|
|---|
| 192 |
|
|---|
| 193 | finish
|
|---|
| 194 |
|
|---|
| 195 |
|
|---|
| 196 | # old query:
|
|---|
| 197 | # #query="SELECT "$timeselect2", Round(o.excrate,2), Round(ExcErr(o.sigevts, o.bgevts),2), IF (o.sigevts>10 AND o.bgevts>10, LiMa(o.sigevts, o.bgevts), '-'), Round(o.ontime,1) as ontime FROM"
|
|---|
| 198 | # query="SELECT "$timeselect2" AS time, ROUND((o.sigevts-o.bgevts)/o.ontime*60,2) AS excrate, ROUND(ExcErr(o.sigevts, o.bgevts)/o.ontime*60,2) AS excerr, ROUND(o.ontime,1) AS ontime, IF (o.sigevts>10 AND o.bgevts>10, LiMa(o.sigevts, o.bgevts), IF (o.sigevts-o.bgevts>10, CONCAT('(', LiMa(o.sigevts, o.bgevts), ')'), '-')) AS significance FROM"
|
|---|
| 199 | # query=$query" (SELECT IF (@os+fOnTimeAfterCuts/60.>20, @start:=fRunStart, IF (@os=0, @start:=fRunStart, @start)) as start, fRunStop as stop, IF (@os+fOnTimeAfterCuts/60.>"$bin", @sig := fNumSigEvts, @sig := @sig + fNumSigEvts) AS sigevts, "
|
|---|
| 200 | # #query=$query" (SELECT fRunID, IF (@os+fOnTimeAfterCuts/60.>"$bin", @sig := fNumSigEvts, @sig := @sig + fNumSigEvts) AS sigevts, "
|
|---|
| 201 | # query=$query" IF (@os+fOnTimeAfterCuts/60.>"$bin", @bg := fNumBgEvts, @bg := @bg + fNumBgEvts) AS bgevts, "
|
|---|
| 202 | # query=$query" IF (@os+fOnTimeAfterCuts/60.>"$bin", @exc := fNumExcEvts, @exc := @exc + fNumExcEvts) AS excevts, "
|
|---|
| 203 | # query=$query" IF (@os+fOnTimeAfterCuts/60.>"$bin", @os := fOnTimeAfterCuts/60., @os := @os + fOnTimeAfterCuts/60.) AS ontime "
|
|---|
| 204 | # #query=$query" @exc/@os*60 AS excrate, "
|
|---|
| 205 | # #query=$query" ExcErr(@sig, @bg) AS excerr, "
|
|---|
| 206 | # #query=$query" IF (@sig>0 AND @bg*5>0 AND (@sig*Log(@sig/(@sig+@bg*5)*1.2/0.2)+@bg*5*Log(@bg*5/(@sig+@bg*5)*1.2))> 0, Sqrt((@sig*Log(@sig/(@sig+@bg*5)*1.2/0.2)+@bg*5*Log(@bg*5/(@sig+@bg*5)*1.2))*2), 0) AS sig "
|
|---|
| 207 | # #query=$query" LiMa(@sig,@bg) AS signif "
|
|---|
| 208 | # query=$query" FROM AnalysisResultsRunLP O LEFT JOIN RunInfo USING (fNight, fRunID) "
|
|---|
| 209 | # query=$query" CROSS JOIN (SELECT @os :=0, @sig :=0, @bg := 0, @exc :=0, @rate :=0, @start :=0) PARAMS "
|
|---|
| 210 | # query=$query" WHERE fNight between "$start" AND "$stop" AND fSourceKey="$source" ORDER BY fNight, fRunID) o HAVING ontime > "$bin"*0.75"
|
|---|
| 211 | #
|
|---|
| 212 |
|
|---|
