- Timestamp:
- 07/26/19 21:43:32 (5 years ago)
- File:
-
- 1 edited
-
trunk/DataCheck/Tools/get_data.sh (modified) (11 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/DataCheck/Tools/get_data.sh
r19568 r19569 49 49 # -------- 50 50 # time: time, delta time, start, stop, ontime 51 # flux: excrate, excerr, CU, CUerr, corrate, corerr, CUcor CUcorerr, fluxcor, fluxcorerr,51 # flux: excrate, excerr, CU, CUerr, flux, fluxerr corrate, corerr, CUcor CUcorerr, fluxcor, fluxcorerr, 52 52 # other info on flux: signif, cu-factor, num exc, num sig, num bg 53 53 # other info: zd th R750cor R750ref … … 231 231 cu=" ROUND("$CU", 2) AS cu" 232 232 cuerr=" ROUND("$CUerr", 2) AS cuerr" 233 cucor=" ROUND("$CUcor", 2) AS cu "234 cucorerr=" ROUND("$CUcorerr", 2) AS cu err"233 cucor=" ROUND("$CUcor", 2) AS cucor" 234 cucorerr=" ROUND("$CUcorerr", 2) AS cucorerr" 235 235 flux="ROUND("$CU" * "$crabflux", 2) AS flux" 236 236 fluxerr="ROUND("$CUerr" * "$crabflux", 2) AS fluxerr" 237 fluxcor="ROUND("$CUcor" * "$crabflux", 2) AS fluxcor" 238 fluxcorerr="ROUND("$CUcorerr" * "$crabflux", 2) AS fluxcorerr" 237 239 # for minute binning 238 240 ontime2=" ROUND(SUM(o.ot)/60., 1) AS ontime" … … 247 249 cu2=" ROUND("$CU2", 2) AS cu" 248 250 cuerr2=" ROUND("$CUerr2", 2) AS cuerr" 249 cucor2=" ROUND("$CUcor2", 2) AS cu "250 cucorerr2=" ROUND("$CUcorerr2", 2) AS cu err"251 cucor2=" ROUND("$CUcor2", 2) AS cucor" 252 cucorerr2=" ROUND("$CUcorerr2", 2) AS cucorerr" 251 253 flux2="ROUND("$CU2" * "$crabflux", "$fluxprec") AS flux" 252 254 fluxerr2="ROUND("$CUerr2" *"$crabflux", "$fluxprec") AS fluxerr" 255 fluxcor2="ROUND("$CUcor2" * "$crabflux", "$fluxprec") AS fluxcor" 256 fluxcorerr2="ROUND("$CUcorerr2" *"$crabflux", "$fluxprec") AS fluxcorerr" 253 257 254 258 case $timeunit in … … 373 377 # internal 374 378 queryint=$querystart 375 queryint=$queryint" "$excrate", "$cu", "$ correxcrate", "$cucor", "$flux", "379 queryint=$queryint" "$excrate", "$cu", "$flux", "$correxcrate", "$cucor", "$fluxcor", " 376 380 queryint=$queryint" "$deltat", "$ontime", " 377 queryint=$queryint" "$excrateerr", "$cuerr", "$ correxcrateerr", "$cucorerr", "$fluxerr", "381 queryint=$queryint" "$excrateerr", "$cuerr", "$fluxerr", "$correxcrateerr", "$cucorerr", "$fluxcorerr", " 378 382 queryint=$queryint" "$significance", " 379 383 queryint=$queryint" MIN(fNight) AS nightmin, MAX(fNight) AS nightmax, " … … 386 390 # for collaborators 387 391 querycol=$querystart 388 querycol=$querycol" "$excrate", "$correxcrate", "$cucor", "$flux ", "392 querycol=$querycol" "$excrate", "$correxcrate", "$cucor", "$fluxcor", " 389 393 querycol=$querycol" "$deltat", "$ontime", " 390 querycol=$querycol" "$excrateerr", "$correxcrateerr", "$cucorerr", "$flux err", "394 querycol=$querycol" "$excrateerr", "$correxcrateerr", "$cucorerr", "$fluxcorerr", " 391 395 querycol=$querycol" "$significance 392 396 querycol=$querycol" "$querybase" "$querydch" "$queryend … … 413 417 # internal 414 418 queryint=$querystart 415 queryint=$queryint" "$excrate2", "$cu2", "$ correxcrate2", "$cucor2", "$flux2", "419 queryint=$queryint" "$excrate2", "$cu2", "$flux2", "$correxcrate2", "$cucor2", "$fluxcor2", " 416 420 queryint=$queryint" "$deltat2", "$ontime2", " 417 queryint=$queryint" "$excrateerr2", "$cuerr2", "$ correxcrateerr2", "$cucorerr2", "$fluxerr2", "421 queryint=$queryint" "$excrateerr2", "$cuerr2", "$fluxerr2", "$correxcrateerr2", "$cucorerr2", "$fluxcorerr2", " 418 422 queryint=$queryint" "$significance2", " 419 423 queryint=$queryint" avg(o.night) AS night, " … … 425 429 # for collaborators 426 430 querycol=$querystart 427 querycol=$querycol" "$excrate2", "$correxcrate2", "$cu2", "$flux 2", "431 querycol=$querycol" "$excrate2", "$correxcrate2", "$cu2", "$fluxcor2", " 428 432 querycol=$querycol" "$deltat2", "$ontime2", " 429 querycol=$querycol" "$excrateerr2", "$correxcrateerr2", "$cuerr2", "$flux err2", "433 querycol=$querycol" "$excrateerr2", "$correxcrateerr2", "$cuerr2", "$fluxcorerr2", " 430 434 querycol=$querycol" "$significance2 431 435 querycol=$querycol" "$querybase … … 475 479 echo "# $queryint" >> $fileint 476 480 echo "#" >> $fileint 477 headerint="# time["$timeunit"] start["$timeunit"] stop["$timeunit"] excess-rate[evts/h] flux[CU] corrected_excess-rate[evts/h] corrected_flux[CU] corrected_flux[e-11/cm2/s] (stop-start)/2["$timeunit"] ontime[min]"478 headerint=$headerint" excess-rate_error[evts/h] flux_error[CU] corrected_excess-rate_error[evts/h] corrected_flux_error[CU] corrected_flux_error[e-11/cm2/s] significance night num_exc num_sig num_bg "481 headerint="# time["$timeunit"] start["$timeunit"] stop["$timeunit"] excess-rate[evts/h] flux[CU] flux[e-11TeV/cm2/s] corrected_excess-rate[evts/h] corrected_flux[CU] corrected_flux[e-11TeV/cm2/s] (stop-start)/2["$timeunit"] ontime[min]" 482 headerint=$headerint" excess-rate_error[evts/h] flux_error[CU] flux_error[e-11TeV/cm2/s] corrected_excess-rate_error[evts/h] corrected_flux_error[CU] corrected_flux_error[e-11/cm2/s] significance night num_exc num_sig num_bg " 479 483 headerint=$headerint" zdmin zdmax thmin thmax avg(cufactor) avg(R750cor) avg(R750ref) " 480 484 echo $headerint >> $fileint … … 498 502 #echo "# $querycol" >> $filecol 499 503 #echo "#" >> $filecol 500 headercol="# time["$timeunit"] start["$timeunit"] stop["$timeunit"] excess-rate[evts/h] corrected_excess-rate[evts/h] corrected_flux[CU] corrected_flux[e-11 /cm2/s] (stop-start)/2["$timeunit"] ontime[min]"501 headercol=$headercol" excess-rate_error[evts/h] corrected_excess-rate_error[evts/h] corrected_flux_error[CU] corrected_flux_error[e-11 /cm2/s] significance "504 headercol="# time["$timeunit"] start["$timeunit"] stop["$timeunit"] excess-rate[evts/h] corrected_excess-rate[evts/h] corrected_flux[CU] corrected_flux[e-11TeV/cm2/s] (stop-start)/2["$timeunit"] ontime[min]" 505 headercol=$headercol" excess-rate_error[evts/h] corrected_excess-rate_error[evts/h] corrected_flux_error[CU] corrected_flux_error[e-11TeV/cm2/s] significance " 502 506 echo $headercol >> $filecol 503 507 #echo "$querycol" … … 651 655 # -------------------------------------------------------------------------------------- # 652 656 657 # some test 658 table="AnalysisResultsRunCutsLC" # CutsLC 659 bin=0 653 660 # Mrk 501 661 source=2 662 name="Mrk501_all_P" 663 get_results 664 665 #exit 666 667 # LC for ICRC 668 table="AnalysisResultsRunCutsLC" # CutsLC 669 bin=-1 670 # Mrk 421 671 source=1 672 name="Mrk421_all_nightly" 673 get_results 674 # Mrk 501 675 source=2 676 name="Mrk501_all_nightly" 677 get_results 678 # 1959 679 source=7 680 name="1959_all_nightly" 681 get_results 682 683 # 2344 684 source=3 685 name="2344_all_nightly" 686 get_results 687 688 # Crab 689 source=5 690 name="Crab_all_nightly" 691 get_results 692 693 694 #exit 695 654 696 # Mrk 501 655 697 source=2 … … 676 718 bin=5 677 719 get_results 678 679 exit680 681 # LC for ICRC682 table="AnalysisResultsRunCutsLC" # CutsLC683 bin=-1684 # Mrk 421685 source=1686 name="Mrk421_all_nightly"687 get_results688 # Mrk 501689 source=2690 name="Mrk501_all_nightly"691 get_results692 # 1959693 source=7694 name="1959_all_nightly"695 get_results696 697 720 698 721 exit
Note:
See TracChangeset
for help on using the changeset viewer.
