| 1 | #!/bin/bash
|
|---|
| 2 |
|
|---|
| 3 | # add:
|
|---|
| 4 | # - fkt f dch value?
|
|---|
| 5 | # - function for correction factor
|
|---|
| 6 | # - x-day binning
|
|---|
| 7 | # < 20121212 data
|
|---|
| 8 | # zd, th for internal
|
|---|
| 9 |
|
|---|
| 10 | # todo: check cu-factor for <20121212)
|
|---|
| 11 |
|
|---|
| 12 | # comparison:
|
|---|
| 13 | # function with select/table: 18.4 sec
|
|---|
| 14 | # function with if: 5.7 sec
|
|---|
| 15 | # function with if:
|
|---|
| 16 | # vorteil: schneller
|
|---|
| 17 | # nachteil: keine history
|
|---|
| 18 | # nachteil: es braucht 2 fkt (qla, isdc)
|
|---|
| 19 |
|
|---|
| 20 | # create function CU(night int, ana tinyint) returns double(4,1) deterministic begin declare cu double; set cu=0; select fCU into cu from CU where fValid=1 and fAnalysis=ana and fNight<night order by fNight desc limit 0,1; return cu; end $$
|
|---|
| 21 |
|
|---|
| 22 | # floor((mjd(timestamp)-mjd(startnight)+0.5)/numdays)
|
|---|
| 23 | # floor((mjd(timestamp)-mjd(startnight)-0.5)/numdays)
|
|---|
| 24 |
|
|---|
| 25 |
|
|---|
| 26 | function get_results()
|
|---|
| 27 | {
|
|---|
| 28 | # some basic query parts
|
|---|
| 29 |
|
|---|
| 30 | # DataCheck (old)
|
|---|
| 31 | ##data with old feedback and/or different bias voltage
|
|---|
| 32 | #query=$query" AND fNight>20120420 AND NOT fNight IN (20120406,20120410,20120503) AND"
|
|---|
| 33 | ## broken bias channel
|
|---|
| 34 | #query=$query" NOT fNight BETWEEN 20121206 AND 20130110"
|
|---|
| 35 | # bg-rate cut
|
|---|
| 36 | zdparam=" pow(0.753833*cos(Radians(fZenithDistanceMean)), 7.647435)*exp(-5.753686*pow(Radians(fZenithDistanceMean),2.089609))"
|
|---|
| 37 | thparam=" pow((if(isnull(fThresholdMinSet),fThresholdMedian,fThresholdMinSet)-329.4203),2)*(-0.0000002044803) "
|
|---|
| 38 | param=" (fNumEvtsAfterBgCuts/5-fNumSigEvts)/fOnTimeAfterCuts - "$zdparam" - "$thparam" "
|
|---|
| 39 | dchold=" -0.085 < ("$param") "
|
|---|
| 40 | dchold=$dchold" AND 0.25 > ("$param") "
|
|---|
| 41 | # Datacheck (new) -> combine
|
|---|
| 42 | dchval=" fNumEvtsAfterBgCuts/(1.41*POW(fZenithDistanceMean*PI()/180,2)+0.975)/(-7.53e-12*POW(10, LOG10(fThresholdMinSet)*3.69)+1.035)/TIME_TO_SEC(TIMEDIFF(fRunStop,fRunStart))/fEffectiveOn "
|
|---|
| 43 | # some semi-automatic datacheck
|
|---|
| 44 | dch=" AND (("$dchval" BETWEEN 0.8 AND 1.7 AND fNight BETWEEN 20140520 AND 20150131) " #A
|
|---|
| 45 | dch=$dch" OR ("$dchval" BETWEEN 0.4 AND 1.6 AND fNight BETWEEN 20150201 AND 20150715) " #B
|
|---|
| 46 | dch=$dch" OR ("$dchval" BETWEEN 0.7 AND 1.4 AND fNight BETWEEN 20150716 AND 20160218) " #C
|
|---|
| 47 | dch=$dch" OR ("$dchval" BETWEEN 0.5 AND 1.0 AND fNight > 20160220) " #D
|
|---|
| 48 | dch=$dch" OR ("$dchold" AND fNight<20140520)) " #old
|
|---|
| 49 |
|
|---|
| 50 | ontime1=" TIME_TO_SEC(TIMEDIFF(fRunStop,fRunStart))*fEffectiveOn "
|
|---|
| 51 | ontime2=" fOnTimeAfterCuts "
|
|---|
| 52 | ontimeif=" IF(ISNULL(fEffectiveOn), "$ontime2", "$ontime1") "
|
|---|
| 53 | from=" FROM RunInfo LEFT JOIN "$table" USING (fNight, fRunID) "
|
|---|
| 54 | # time range and source
|
|---|
| 55 | where=" WHERE fSourceKey="$source" AND fNight BETWEEN "$nightmin" AND "$nightmax
|
|---|
| 56 | # some sanity checks
|
|---|
| 57 | where=$where" AND fRunTypeKey=1 "
|
|---|
| 58 | # where=$where" AND NOT ISNULL(fNumSigEvts) AND NOT ISNULL(fNumBgEvts) "
|
|---|
| 59 | # where=$where" AND NOT fRunStart='0000-00-00 00:00:00' AND NOT fRunStop='0000-00-00 00:00:00' "
|
|---|
| 60 | # zd cut
|
|---|
| 61 | where=$where" AND fZenithDistanceMax < "$zdmax
|
|---|
| 62 | # th cut
|
|---|
| 63 | where=$where" AND fThresholdMedian < "$thmax
|
|---|
| 64 | where=$where" "$dch
|
|---|
| 65 |
|
|---|
| 66 | cufactor=" Avg(25.2) "
|
|---|
| 67 | crabflux="3.37e-11"
|
|---|
| 68 | fluxprec=13
|
|---|
| 69 | crabflux="3.37"
|
|---|
| 70 | fluxprec=2
|
|---|
| 71 | # crabflux="2.8e-11"
|
|---|
| 72 | # crabflux="3.9e-11"
|
|---|
| 73 | # range of crab fluxes:
|
|---|
| 74 | # Dortmund: 2.8e-11
|
|---|
| 75 | # HESS: 3.37e-11
|
|---|
| 76 | # HAWC: 3.01e-11
|
|---|
| 77 | # Wue: 3.39e-11 - 3.9e-11 (ISDC analysis)
|
|---|
| 78 | # 15-20% difference
|
|---|
| 79 |
|
|---|
| 80 | case $timeunit in
|
|---|
| 81 | mjd) start=" Mjd(Min(fRunStart)) "
|
|---|
| 82 | stop=" Mjd(MAX(fRunStop)) "
|
|---|
| 83 | deltat=" (Mjd(MAX(fRunStop))-Mjd(Min(fRunStart)))/2 "
|
|---|
| 84 | time=" Mjd(Min(fRunStart))+"$deltat
|
|---|
| 85 | start2=" Mjd(MIN(o.start)) "
|
|---|
| 86 | stop2=" Mjd(MAX(o.stop)) "
|
|---|
| 87 | deltat2=" (Mjd(MAX(o.stop))-Mjd(MIN(o.start)))/2 "
|
|---|
| 88 | time2=" Mjd(MIN(o.start))+"$deltat2
|
|---|
| 89 | ;;
|
|---|
| 90 | unix) start="Unix_timestamp(CONVERT_TZ(Min(fRunStart), '+00:00', 'SYSTEM')) "
|
|---|
| 91 | stop="Unix_timestamp(CONVERT_TZ(Max(fRunStop), '+00:00', 'SYSTEM')) "
|
|---|
| 92 | deltat=" (Unix_timestamp(CONVERT_TZ(Max(fRunStop), '+00:00', 'SYSTEM')) - Unix_timestamp(CONVERT_TZ(Min(fRunStart), '+00:00', 'SYSTEM')))/2 "
|
|---|
| 93 | time=" Unix_timestamp(CONVERT_TZ(Min(fRunStart), '+00:00', 'SYSTEM'))+"$deltat
|
|---|
| 94 | startstop2=" Unix_timestamp(CONVERT_TZ(MIN(o.start), '+00:00', 'SYSTEM')) AS start, "
|
|---|
| 95 | startstop2=$starstop2" Unix_timestamp(CONVERT_TZ(MAX(o.stop), '+00:00', 'SYSTEM')) AS stop, "
|
|---|
| 96 | time2=" (Unix_timestamp(CONVERT_TZ(Max(o.stop), '+00:00', 'SYSTEM')) - Unix_timestamp(CONVERT_TZ(Min(o.start), '+00:00', 'SYSTEM')))/2 "
|
|---|
| 97 | time2=" Unix_timestamp(CONVERT_TZ(Min(o.start), '+00:00', 'SYSTEM'))+"$deltat2
|
|---|
| 98 | ;;
|
|---|
| 99 | *) start=" MIN(fRunStart) "
|
|---|
| 100 | stop=" MAX(fRunStop) "
|
|---|
| 101 | deltat=" sec_to_time(time_to_sec(timediff(MAX(fRunStop), Min(fRunStart)))/2) "
|
|---|
| 102 | time=" addtime(Min(fRunStart), "$deltat") "
|
|---|
| 103 | start2=" MIN(o.start) "
|
|---|
| 104 | stop2=" MAX(o.stop) "
|
|---|
| 105 | deltat2=" sec_to_time(time_to_sec(timediff(MAX(o.stop), Min(o.start)))/2) "
|
|---|
| 106 | time2=" addtime(Min(o.start), "$deltat2") "
|
|---|
| 107 | ;;
|
|---|
| 108 | esac
|
|---|
| 109 | ontime=" SUM("$ontimeif")/60."
|
|---|
| 110 | ontime2=" SUM(o.ot)/60. "
|
|---|
| 111 |
|
|---|
| 112 | excrate=" SUM(fNumExcEvts)/SUM("$ontimeif")*3600 "
|
|---|
| 113 | excerr="ExcErr(Sum(fNumSigEvts), SUM(fNumBgEvts))"
|
|---|
| 114 | significance="LiMa(Sum(fNumSigEvts), SUM(fNumBgEvts))"
|
|---|
| 115 | numexc="Sum(fNumExcEvts)"
|
|---|
| 116 | numsig="Sum(fNumSigEvts)"
|
|---|
| 117 | numbg="Sum(fNumBgEvts)"
|
|---|
| 118 | excrateerr=" "$excerr"/SUM("$ontimeif")*3600 "
|
|---|
| 119 | # thomas correction factor
|
|---|
| 120 | correvts=" fNumExcEvts*(pow(cos(fZenithDistanceMean*PI()/180),3)+14.8/21.9*pow(sin(2*fZenithDistanceMean*PI()/180),5))/((1-0.00124/1.21*(if(isnull(fThresholdMinSet),fThresholdMedian,fThresholdMinSet)-500)*(if(isnull(fThresholdMinSet),fThresholdMedian,fThresholdMinSet)>=500))) "
|
|---|
| 121 | correxcrate=" SUM("$correvts")/SUM("$ontimeif")*3600 "
|
|---|
| 122 | # corerr = MMath::ErrorExc(excevtssum+bgevtssum, bgevtssum*5, 0.2)/ontimesum*3600.*corrate/excrate;
|
|---|
| 123 | correxcrateerr=" "$excerr"/SUM("$ontimeif")*3600*SUM("$correvts")/SUM(fNumExcEvts) "
|
|---|
| 124 | # correction on run basis (not relevant for hess)
|
|---|
| 125 | #cu=$correxcrate"/"$cufactor
|
|---|
| 126 | cu=" SUM("$correvts"/CUQLA(fNight))/SUM("$ontimeif")*3600 "
|
|---|
| 127 | #cuerr=$correxcrateerr"/"$cufactor
|
|---|
| 128 | cuerr=" "$excerr"/SUM("$ontimeif")*3600*SUM("$correvts"/CUQLA(fNight))/SUM(fNumExcEvts) "
|
|---|
| 129 | flux=$cu" * "$crabflux
|
|---|
| 130 | fluxerr=$cuerr" * "$crabflux
|
|---|
| 131 |
|
|---|
| 132 |
|
|---|
| 133 | excrate2=" (SUM(o.sigevts)-SUM(o.bgevts))/SUM(o.ot)*3600 "
|
|---|
| 134 | excerr2="ExcErr(SUM(o.sigevts),SUM(o.bgevts))"
|
|---|
| 135 | significance2="LiMa(SUM(o.sigevts),SUM(o.bgevts))"
|
|---|
| 136 | numexc2="Sum(o.sigevts-o.bgevts)"
|
|---|
| 137 | numsig2="Sum(o.sigevts)"
|
|---|
| 138 | numbg2="Sum(o.bgevts)"
|
|---|
| 139 | excrateerr2=" "$excerr2"/SUM(o.ot)*3600 "
|
|---|
| 140 | correxcrate2=" SUM(o.corevts)/SUM(o.ot)*3600 "
|
|---|
| 141 | correxcrateerr2=" "$excerr2"/SUM(o.ot)*3600*SUM(o.corevts)/(SUM(o.sigevts)-SUM(o.bgevts)) "
|
|---|
| 142 | #cu2=$correxcrate2"/"$cufactor
|
|---|
| 143 | cu2=" SUM(o.corevts/o.cu)/SUM(o.ot)*3600 "
|
|---|
| 144 | #cuerr2=$correxcrateerr2"/"$cufactor
|
|---|
| 145 | cuerr2=" "$excerr2"/SUM(o.ot)*3600*SUM(o.corevts/o.cu)/(SUM(o.sigevts)-SUM(o.bgevts)) "
|
|---|
| 146 | flux2="$cu2*"$crabflux
|
|---|
| 147 | fluxerr2="$cuerr2*"$crabflux
|
|---|
| 148 |
|
|---|
| 149 | # order information such that it is easily readable by tgraph
|
|---|
| 150 | # tgraph: X, Y
|
|---|
| 151 | # tgraph errors: X, Y, EX, EY
|
|---|
| 152 |
|
|---|
| 153 | # internal
|
|---|
| 154 | # --------
|
|---|
| 155 | # timeselect:
|
|---|
| 156 | # mjdstar, mjdstop, mjdmean, ontime
|
|---|
| 157 | # excselect:
|
|---|
| 158 | # excrate, excerr
|
|---|
| 159 | # corrected: excrate, excerr
|
|---|
| 160 | # CU CUerr
|
|---|
| 161 | # flux, fluxerr
|
|---|
| 162 | # addselect:
|
|---|
| 163 | # signif
|
|---|
| 164 | # num exc, num sig, num bg
|
|---|
| 165 | # other info: zd? th?
|
|---|
| 166 | #
|
|---|
| 167 | #
|
|---|
| 168 | # exter nal
|
|---|
| 169 | # --------
|
|---|
| 170 | # time, delta time, start, stop
|
|---|
| 171 | # corr-excrate, corr-excerr
|
|---|
| 172 | # flux, flux-err
|
|---|
| 173 |
|
|---|
| 174 | if [ $bin -le 0 ]
|
|---|
| 175 | then
|
|---|
| 176 | queryint="SELECT "
|
|---|
| 177 | if [ $bin -eq 0 ]
|
|---|
| 178 | then
|
|---|
| 179 | queryint=$queryint" fPeriod as num, "
|
|---|
| 180 | else
|
|---|
| 181 | queryint=$queryint" FLOOR((Mjd(fRunStart)-Mjd("$nightmin")-0.5)/"`echo $bin | sed -e 's/-//'`".) as num, "
|
|---|
| 182 | fi
|
|---|
| 183 | queryint=$queryint" "$time" as time, "$start" as start, "$stop" as stop, "
|
|---|
| 184 | queryint=$queryint" round("$excrate", 1) as excrate, round("$correxcrate", 1) as correxcrate, "
|
|---|
| 185 | queryint=$queryint" round($cu, 1) as cu, $flux as flux, "
|
|---|
| 186 | queryint=$queryint" "$deltat" as deltat, round("$ontime", 1) as ontime, "
|
|---|
| 187 | queryint=$queryint" round("$excrateerr", 1) as excrateerr, round("$correxcrateerr", 1) as correxcrateerr, "
|
|---|
| 188 | queryint=$queryint" round($cuerr, 1) as cuerr, $fluxerr as fluxerr, "
|
|---|
| 189 | queryint=$queryint" round("$significance", 1) as significance, "
|
|---|
| 190 | queryint=$queryint" Min(fNight) as nightmin, "
|
|---|
| 191 | queryint=$queryint" Max(fNight) as nightmax, "
|
|---|
| 192 | queryint=$queryint" "$numexc" as numexc, "
|
|---|
| 193 | queryint=$queryint" "$numsig" as numsig, "
|
|---|
| 194 | queryint=$queryint" "$numbg" as numbg "
|
|---|
| 195 | #queryint=$queryint" now() as now "
|
|---|
| 196 |
|
|---|
| 197 | queryext="SELECT "
|
|---|
| 198 | if [ $bin -eq 0 ]
|
|---|
| 199 | then
|
|---|
| 200 | queryext=$queryext" fPeriod as num, "
|
|---|
| 201 | else
|
|---|
| 202 | queryext=$queryext" FLOOR((Mjd(fRunStart)-Mjd("$nightmin")-0.5)/"`echo $bin | sed -e 's/-//'`".) as num, "
|
|---|
| 203 | fi
|
|---|
| 204 | queryext=$queryext" "$time" as time, "$start" as start, "$stop" as stop, "
|
|---|
| 205 | queryext=$queryext" round("$correxcrate", 1) as correxcrate, round($flux, "$fluxprec") as flux, "
|
|---|
| 206 | queryext=$queryext" "$deltat" as deltat, round("$ontime", 1) as ontime, "
|
|---|
| 207 | queryext=$queryext" round("$correxcrateerr", 1) as correxcrateerr, round($fluxerr, "$fluxprec") as fluxerr, "
|
|---|
| 208 | queryext=$queryext" round("$significance", 1) as significance "
|
|---|
| 209 | #queryext=$queryext" now() as now "
|
|---|
| 210 |
|
|---|
| 211 | querybase=$from$where
|
|---|
| 212 | #if [ $bin -eq -1 ]
|
|---|
| 213 | #then
|
|---|
| 214 | # querybase=$querybase" GROUP BY fNight "
|
|---|
| 215 | #fi
|
|---|
| 216 | #querybase=$querybase" GROUP BY FLOOR((Mjd(fRunStart)-Mjd("$nightmin")-0.5)/"`echo $bin | sed -e 's/-//'`".)"
|
|---|
| 217 | querybase=$querybase" GROUP BY num "
|
|---|
| 218 | #echo "-----------------"$querybase
|
|---|
| 219 | if [ "$ontimelimit" = "" ]
|
|---|
| 220 | then
|
|---|
| 221 | querybase=$querybase" HAVING SUM("$ontimeif")>1200 AND NOT SUM(fNumBgEvts)=SUM(fNumSigEvts) ORDER BY num " # 20 min
|
|---|
| 222 | else
|
|---|
| 223 | querybase=$querybase" HAVING SUM("$ontimeif")>"$ontimelimit" AND NOT SUM(fNumBgEvts)=SUM(fNumSigEvts) ORDER BY num "
|
|---|
| 224 | fi
|
|---|
| 225 |
|
|---|
| 226 | queryint=$queryint" "$querybase
|
|---|
| 227 | queryext=$queryext" "$querybase
|
|---|
| 228 | else
|
|---|
| 229 | queryint="SELECT "
|
|---|
| 230 | queryint=$queryint" "$time2" as time, "$start2" as start, "$stop2" as stop, "
|
|---|
| 231 | queryint=$queryint" round("$excrate2", 1) as excrate, round("$correxcrate2", 1) as correxcrate, "
|
|---|
| 232 | queryint=$queryint" round("$cu2", 1) as cu, round("$flux2", "$fluxprec") as flux, "
|
|---|
| 233 | queryint=$queryint" round("$excrateerr2", 1) as excrateerr, round("$correxcrateerr2", 1) as correxcrateerr, "
|
|---|
| 234 | queryint=$queryint" "$deltat2" as deltat, round("$ontime2", 1) as ontime, "
|
|---|
| 235 | queryint=$queryint" round("$cuerr2", 1) as cuerr, round("$fluxerr2", "$fluxprec") as fluxerr, "
|
|---|
| 236 | queryint=$queryint" round("$significance2", 1) as significance, "
|
|---|
| 237 | queryint=$queryint" avg(o.night) as night, "
|
|---|
| 238 | queryint=$queryint" "$numexc2" as numexc, "
|
|---|
| 239 | queryint=$queryint" "$numsig2" as numsig, "
|
|---|
| 240 | queryint=$queryint" "$numbg2" as numbg "
|
|---|
| 241 | #queryint=$queryint" now() as now "
|
|---|
| 242 |
|
|---|
| 243 | queryext="SELECT "
|
|---|
| 244 | queryext=$queryext" "$time2" as time, "$start2" as start, "$stop2" as stop, "
|
|---|
| 245 | queryext=$queryext" round("$correxcrate2", 1) as correxcrate, round("$flux2", "$fluxprec") as flux, "
|
|---|
| 246 | queryext=$queryext" "$deltat2" as deltat, round("$ontime2", 1) as ontime, "
|
|---|
| 247 | queryext=$queryext" round("$correxcrateerr2", 1) as correxcrateerr, round("$fluxerr2", "$fluxprec") as fluxerr, "
|
|---|
| 248 | queryext=$queryext" round("$significance2", 1) as significance "
|
|---|
| 249 | #queryext=$queryext" now() as now "
|
|---|
| 250 |
|
|---|
| 251 | querybase=" FROM (SELECT fNight, @ot:="$ontimeif" AS ot, fRunStart AS start, fRunStop AS stop, fNumSigEvts AS sigevts, fNumBgEvts AS bgevts, "
|
|---|
| 252 | querybase=$querybase" "$correvts" AS corevts, CUQLA(fNight) AS cu, "
|
|---|
| 253 | querybase=$querybase" IF (@night=fNight AND FLOOR((@os+@ot)/"$bin"./60.)<1, @bl, @bl := @bl + 1) AS block, "
|
|---|
| 254 | querybase=$querybase" IF (@night=fNight AND FLOOR((@os+@ot)/"$bin"./60.)<1, @os:=@os + @ot, @os := @ot) AS os, @night :=fNight AS night "
|
|---|
| 255 | querybase=$querybase$from" CROSS JOIN (SELECT @night :=0, @ot :=0, @os :=0, @bl:=0) PARAMS "
|
|---|
| 256 | querybase=$querybase$where" ORDER BY fRunStart) o GROUP BY block HAVING NOT sum(o.bgevts)=sum(o.sigevts) AND ontime>0.75*"$bin
|
|---|
| 257 | # AND NOT SUM(fNumBgEvts)=SUM(fNumSigEvts)
|
|---|
| 258 |
|
|---|
| 259 | queryint=$queryint" "$querybase
|
|---|
| 260 | queryext=$queryext" "$querybase
|
|---|
| 261 | fi
|
|---|
| 262 |
|
|---|
| 263 |
|
|---|
| 264 | fileint=$path"/data/FACT_preliminary_"$name"_internal.dat"
|
|---|
| 265 | if [ "$overwrite" = "yes" ]
|
|---|
| 266 | then
|
|---|
| 267 | echo "internal: "$fileint
|
|---|
| 268 | echo "# this file was created at "`date` > $fileint
|
|---|
| 269 | fi
|
|---|
| 270 | if [ $bin -le 0 ]
|
|---|
| 271 | then
|
|---|
| 272 | echo "# numbin time[mjd] start[mjd] stop[mjd] excrate[evts/h] corr.excrate[evts/h] flux[CU] flux[e-11/cm2/s] delta_time[mjd] ontime[min] excrate_err[evts/h] corr.excrate_err[evts/h] flux_err[CU] flux_err[e-11/cm2/s] significance nightmin, nightmax num_exc num_sig num_bg " >> $fileint
|
|---|
| 273 | else
|
|---|
| 274 | echo "# time[mjd] start[mjd] stop[mjd] excrate[evts/h] corr.excrate[evts/h] flux[CU] flux[e-11/cm2/s] delta_time[mjd] ontime[min] excrate_err[evts/h] corr.excrate_err[evts/h] flux_err[CU] flux_err[e-11/cm2/s] significance nightmin, nightmax num_exc num_sig num_bg " >> $fileint
|
|---|
| 275 | fi
|
|---|
| 276 | #echo "$queryint"
|
|---|
| 277 | mysql --defaults-file=$sqlpw -u root $dbname -s -e "$queryint" >> $fileint
|
|---|
| 278 | #mysql --defaults-file=$sqlpw -u root $dbname -e "$queryint"
|
|---|
| 279 |
|
|---|
| 280 |
|
|---|
| 281 | fileext=$path"/data/FACT_preliminary_"$name".dat"
|
|---|
| 282 | if [ "$overwrite" = "yes" ]
|
|---|
| 283 | then
|
|---|
| 284 | echo "external: "$fileext
|
|---|
| 285 | echo "# this file was created at "`date` > $fileext
|
|---|
| 286 | fi
|
|---|
| 287 | if [ $bin -lt 0 ]
|
|---|
| 288 | then
|
|---|
| 289 | echo "# numbin time[mjd] start[mjd] stop[mjd] corr.excrate[evts/h] flux[e-11/cm2/s] delta_time[mjd] ontime[min] corr.excrate_err[evts/h] flux_err[e-11/cm2/s] significance " >> $fileext
|
|---|
| 290 | else
|
|---|
| 291 | echo "# time[mjd] start[mjd] stop[mjd] corr.excrate[evts/h] flux[e-11/cm2/s] delta_time[mjd] ontime[min] corr.excrate_err[evts/h] flux_err[e-11/cm2/s] significance " >> $fileext
|
|---|
| 292 | fi
|
|---|
| 293 | #echo "$queryext"
|
|---|
| 294 | mysql --defaults-file=$sqlpw -u root $dbname -s -e "$queryext" >> $fileext
|
|---|
| 295 | #mysql --defaults-file=$sqlpw -u root $dbname -e "$queryext"
|
|---|
| 296 | }
|
|---|
| 297 |
|
|---|
| 298 | # setup
|
|---|
| 299 | # db
|
|---|
| 300 | sqlpw=/home/$USER/.mysql.pw
|
|---|
| 301 | dbname=factdata20170804
|
|---|
| 302 | # selection
|
|---|
| 303 | timeunit=mjd
|
|---|
| 304 | #bin=20 # min
|
|---|
| 305 | #bin=0 # period
|
|---|
| 306 | #bin=-365 # yearly
|
|---|
| 307 | bin=-1 # nightly
|
|---|
| 308 | zdmax=90
|
|---|
| 309 | thmax=1500
|
|---|
| 310 | path=`dirname $0`
|
|---|
| 311 | table="AnalysisResultsRunLP"
|
|---|
| 312 | table="AnalysisResultsAllQLA"
|
|---|
| 313 | overwrite="yes"
|
|---|
| 314 | bin=-1 # nightly
|
|---|
| 315 | nightmin=20111115
|
|---|
| 316 | nightmax=20171231
|
|---|
| 317 |
|
|---|
| 318 |
|
|---|
| 319 | # Mrk 421
|
|---|
| 320 | source=1
|
|---|
| 321 | name="Mrk421_nightly"
|
|---|
| 322 | bin=-1
|
|---|
| 323 | get_results
|
|---|
| 324 | name="Mrk421_20min"
|
|---|
| 325 | bin=20
|
|---|
| 326 | get_results
|
|---|
| 327 | name="Mrk421_3d"
|
|---|
| 328 | bin=-3
|
|---|
| 329 | get_results
|
|---|
| 330 | name="Mrk421_10d"
|
|---|
| 331 | bin=-10
|
|---|
| 332 | get_results
|
|---|
| 333 | name="Mrk421_period"
|
|---|
| 334 | bin=0
|
|---|
| 335 | get_results
|
|---|
| 336 |
|
|---|
| 337 |
|
|---|
| 338 |
|
|---|
| 339 | # Mrk 501
|
|---|
| 340 | source=2
|
|---|
| 341 | name="Mrk501_nightly"
|
|---|
| 342 | bin=-1
|
|---|
| 343 | get_results
|
|---|
| 344 | name="Mrk501_20min"
|
|---|
| 345 | bin=20
|
|---|
| 346 | get_results
|
|---|
| 347 | name="Mrk501_3d"
|
|---|
| 348 | bin=-3
|
|---|
| 349 | get_results
|
|---|
| 350 | name="Mrk501_10d"
|
|---|
| 351 | bin=-10
|
|---|
| 352 | get_results
|
|---|
| 353 | name="Mrk501_period"
|
|---|
| 354 | bin=0
|
|---|
| 355 | get_results
|
|---|
| 356 |
|
|---|
| 357 |
|
|---|
| 358 |
|
|---|
| 359 | # 2344
|
|---|
| 360 | source=3
|
|---|
| 361 | name="2344_nightly"
|
|---|
| 362 | bin=-1
|
|---|
| 363 | get_results
|
|---|
| 364 | name="2344_20min"
|
|---|
| 365 | bin=20
|
|---|
| 366 | get_results
|
|---|
| 367 | name="2344_period"
|
|---|
| 368 | bin=0
|
|---|
| 369 | get_results
|
|---|
| 370 |
|
|---|
| 371 |
|
|---|
| 372 |
|
|---|
| 373 | # 1959
|
|---|
| 374 | source=7
|
|---|
| 375 | name="1959_nightly"
|
|---|
| 376 | bin=-1
|
|---|
| 377 | get_results
|
|---|
| 378 | name="1959_20min"
|
|---|
| 379 | bin=20
|
|---|
| 380 | get_results
|
|---|
| 381 | name="1959_period"
|
|---|
| 382 | bin=0
|
|---|
| 383 | get_results
|
|---|
| 384 |
|
|---|
| 385 |
|
|---|
| 386 |
|
|---|
| 387 | # 0323
|
|---|
| 388 | source=12
|
|---|
| 389 | name="0323_nightly"
|
|---|
| 390 | bin=-1
|
|---|
| 391 | get_results
|
|---|
| 392 | name="0323_20min"
|
|---|
| 393 | bin=20
|
|---|
| 394 | get_results
|
|---|
| 395 | name="0323_period"
|
|---|
| 396 | bin=0
|
|---|
| 397 | get_results
|
|---|
| 398 |
|
|---|
| 399 |
|
|---|
| 400 |
|
|---|
| 401 | # crab
|
|---|
| 402 | source=5
|
|---|
| 403 | name="Crab_nightly"
|
|---|
| 404 | bin=-1
|
|---|
| 405 | get_results
|
|---|
| 406 | name="Crab_20min"
|
|---|
| 407 | bin=20
|
|---|
| 408 | get_results
|
|---|
| 409 | name="Crab_period"
|
|---|
| 410 | bin=0
|
|---|
| 411 | get_results
|
|---|
| 412 | name="Crab_season"
|
|---|
| 413 | bin=-365
|
|---|
| 414 | nightmin=20110716
|
|---|
| 415 | nightmax=20180716
|
|---|
| 416 | get_results
|
|---|
| 417 |
|
|---|
| 418 |
|
|---|
| 419 |
|
|---|
| 420 | exit
|
|---|
| 421 |
|
|---|
| 422 | # more examples
|
|---|
| 423 |
|
|---|
| 424 | name="1959_2016"
|
|---|
| 425 | source=7
|
|---|
| 426 | bin=-1
|
|---|
| 427 | nightmin=20160201
|
|---|
| 428 | nightmax=20161105
|
|---|
| 429 | get_results
|
|---|
| 430 |
|
|---|
| 431 | name="1959_all_variable"
|
|---|
| 432 | overwrite="no"
|
|---|
| 433 | source=7
|
|---|
| 434 | bin=-365
|
|---|
| 435 | nightmin=20120201
|
|---|
| 436 | nightmax=20130131
|
|---|
| 437 | get_results
|
|---|
| 438 | nightmin=20130201
|
|---|
| 439 | nightmax=20140131
|
|---|
| 440 | get_results
|
|---|
| 441 | nightmin=20140201
|
|---|
| 442 | nightmax=20150131
|
|---|
| 443 | get_results
|
|---|
| 444 | bin=0
|
|---|
| 445 | nightmin=20150201
|
|---|
| 446 | nightmax=20160131
|
|---|
| 447 | get_results
|
|---|
| 448 | bin=-1
|
|---|
| 449 | nightmin=20160201
|
|---|
| 450 | nightmax=20170131
|
|---|
| 451 | get_results
|
|---|
| 452 | bin=0
|
|---|
| 453 | nightmin=20170201
|
|---|
| 454 | nightmax=20180131
|
|---|
| 455 | get_results
|
|---|
| 456 |
|
|---|
| 457 |
|
|---|
| 458 |
|
|---|
| 459 | overwrite="yes"
|
|---|
| 460 | name="1959_all_variable2"
|
|---|
| 461 | overwrite="no"
|
|---|
| 462 | source=7
|
|---|
| 463 | bin=-365
|
|---|
| 464 | nightmin=20120201
|
|---|
| 465 | nightmax=20130131
|
|---|
| 466 | get_results
|
|---|
| 467 | nightmin=20130201
|
|---|
| 468 | nightmax=20140131
|
|---|
| 469 | get_results
|
|---|
| 470 | nightmin=20140201
|
|---|
| 471 | nightmax=20150131
|
|---|
| 472 | get_results
|
|---|
| 473 | bin=0
|
|---|
| 474 | nightmin=20150201
|
|---|
| 475 | nightmax=20160131
|
|---|
| 476 | get_results
|
|---|
| 477 | bin=-1
|
|---|
| 478 | nightmin=20160201
|
|---|
| 479 | nightmax=20160817
|
|---|
| 480 | get_results
|
|---|
| 481 | bin=0
|
|---|
| 482 | nightmin=20160818
|
|---|
| 483 | nightmax=20180131
|
|---|
| 484 | get_results
|
|---|
| 485 |
|
|---|
| 486 |
|
|---|
| 487 |
|
|---|
| 488 | exit
|
|---|
| 489 |
|
|---|
| 490 | overwrite="yes"
|
|---|
| 491 |
|
|---|
| 492 | bin=0
|
|---|
| 493 | source=3
|
|---|
| 494 | name="2344period"
|
|---|
| 495 | get_results
|
|---|
| 496 |
|
|---|
| 497 | exit
|
|---|
| 498 |
|
|---|
| 499 | # flare night (HESS)
|
|---|
| 500 | name="Mrk501_10min_flarenight"
|
|---|
| 501 | source=2
|
|---|
| 502 | bin=10
|
|---|
| 503 | nightmin=20140623
|
|---|
| 504 | nightmax=20140623
|
|---|
| 505 | get_results
|
|---|
| 506 |
|
|---|
| 507 |
|
|---|
| 508 | exit
|
|---|
| 509 |
|
|---|
| 510 | # flare night (HESS)
|
|---|
| 511 | name="Mrk501_5min_flarenight"
|
|---|
| 512 | source=2
|
|---|
| 513 | bin=5
|
|---|
| 514 | nightmin=20140623
|
|---|
| 515 | nightmax=20140623
|
|---|
| 516 | get_results
|
|---|
| 517 |
|
|---|
| 518 |
|
|---|
| 519 | exit
|
|---|
| 520 |
|
|---|
| 521 |
|
|---|
| 522 | # full sample
|
|---|
| 523 | name="Mrk421_all_nightly"
|
|---|
| 524 | source=1
|
|---|
| 525 | get_results
|
|---|
| 526 |
|
|---|
| 527 | name="Mrk501_all_nightly"
|
|---|
| 528 | source=2
|
|---|
| 529 | get_results
|
|---|
| 530 |
|
|---|
| 531 | name="1959_all_nightly"
|
|---|
| 532 | source=7
|
|---|
| 533 | get_results
|
|---|
| 534 |
|
|---|
| 535 | name="2344_all_nightly"
|
|---|
| 536 | source=3
|
|---|
| 537 | get_results
|
|---|
| 538 |
|
|---|
| 539 | exit
|
|---|
| 540 |
|
|---|
| 541 |
|
|---|
| 542 | name="HESE20160427"
|
|---|
| 543 | source=19
|
|---|
| 544 | nightmin=20160425
|
|---|
| 545 | bin=-10
|
|---|
| 546 | get_results
|
|---|
| 547 |
|
|---|
| 548 | name="AMON20160731"
|
|---|
| 549 | source=21
|
|---|
| 550 | nightmin=20160730
|
|---|
| 551 | bin=-10
|
|---|
| 552 | get_results
|
|---|
| 553 |
|
|---|
| 554 |
|
|---|
| 555 | # full sample
|
|---|
| 556 | name="Mrk501_all_3d"
|
|---|
| 557 | source=2
|
|---|
| 558 | bin=-3
|
|---|
| 559 | get_results
|
|---|
| 560 |
|
|---|
| 561 | exit
|
|---|
| 562 |
|
|---|
| 563 | # full sample
|
|---|
| 564 | name="Mrk421_all_3d"
|
|---|
| 565 | source=1
|
|---|
| 566 | bin=-3
|
|---|
| 567 | get_results
|
|---|
| 568 |
|
|---|
| 569 | name="1959_2016_nightly"
|
|---|
| 570 | source=7
|
|---|
| 571 | nightmin=20160301
|
|---|
| 572 | nightmax=20160831
|
|---|
| 573 | get_results
|
|---|
| 574 |
|
|---|
| 575 | # full sample
|
|---|
| 576 | name="Mrk421_all_nightly"
|
|---|
| 577 | source=1
|
|---|
| 578 | get_results
|
|---|
| 579 |
|
|---|
| 580 | name="Mrk501_all_nightly"
|
|---|
| 581 | source=2
|
|---|
| 582 | get_results
|
|---|
| 583 |
|
|---|
| 584 | name="1959_all_nightly"
|
|---|
| 585 | source=7
|
|---|
| 586 | get_results
|
|---|
| 587 |
|
|---|
| 588 | bin=20 # min
|
|---|
| 589 |
|
|---|
| 590 | # full sample
|
|---|
| 591 | name="Mrk421_all_20min"
|
|---|
| 592 | source=1
|
|---|
| 593 | get_results
|
|---|
| 594 |
|
|---|
| 595 | name="Mrk501_all_20min"
|
|---|
| 596 | source=2
|
|---|
| 597 | get_results
|
|---|
| 598 |
|
|---|
| 599 | name="1959_all_20min"
|
|---|
| 600 | source=7
|
|---|
| 601 | get_results
|
|---|
| 602 |
|
|---|
| 603 | exit
|
|---|
| 604 |
|
|---|
| 605 | name="Mrk501_all_weekly"
|
|---|
| 606 | nightmin=20140524
|
|---|
| 607 | nightmax=20140930
|
|---|
| 608 | source=2
|
|---|
| 609 | bin=-7
|
|---|
| 610 | get_results
|
|---|
| 611 |
|
|---|
| 612 | exit
|
|---|
| 613 |
|
|---|
| 614 |
|
|---|
