Changeset 19560


Ignore:
Timestamp:
Jul 7, 2019, 11:53:52 AM (10 days ago)
Author:
dorner
Message:
corrected conversion factors from CU to fluxes
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/DataCheck/Tools/get_data.sh

    r19448 r19560  
    156156   
    157157   # these values use the Crab spectrum from HESS 2006
     158   # hess 2006 http://www.aanda.org/articles/aa/pdf/2006/39/aa5351-06.pdf
     159   # PL.SetParameter(1,-2.63);
     160   # PL.SetParameter(0, 3.45e-11);
    158161   case "$table" in
    159162   "AnalysisResultsRunISDC")
    160163      zdfactor="pow(cos(fZenithDistanceMean*PI()/180)*exp(1-cos(fZenithDistanceMean*PI()/180)),4.2)"
    161164      thfactor="(1.18-IF(ISNULL(fThresholdMinSet),fThresholdMedian,fThresholdMinSet)*0.00062)"
    162       # ETh: 764 GeV
    163       crabflux="3.3"
     165      # ETh: 1060 GeV
     166      crabflux="1.9"
    164167      ;;
    165168   "AnalysisResultsRunCutsLC")
    166169      zdfactor="pow(cos(fZenithDistanceMean*PI()/180)*exp(1-cos(fZenithDistanceMean*PI()/180)),4.5)"
    167170      thfactor="(1.37-IF(ISNULL(fThresholdMinSet),fThresholdMedian,fThresholdMinSet)*0.00118)"
    168       # ETh: 543 GeV
    169       crabflux="5.7"
     171      # ETh: 775 GeV
     172      crabflux="3.2"
    170173      ;;
    171174   "AnalysisResultsRunLP")
     
    173176      #zdfactor="(1/(pow(cos(fZenithDistanceMean*PI()/180),3)+14.8/21.9*pow(sin(2*fZenithDistanceMean*PI()/180),5)))"
    174177      thfactor="(1-0.00124/1.21*(fThresholdMinSet-500)*(fThresholdMinSet>=500))"                                 
    175       # ETh: 715 GeV
    176       crabflux="3.6"
     178      # ETh: 1090 GeV
     179      crabflux="1.8"
    177180      ;;
    178181   "*")
     
    569572# path to file with mysql setup (user, password, host, db-name)
    570573sqlpw=/home/$USER/.mysql.pw.local
     574sqlpw=/home/$USER/.mysql.pw.ethz.fact
    571575#
    572576# -------------
     
    639643# -------------------------------------------------------------------------------------- #
    640644
     645
     646# LC for INTEGRAL Proposal
     647table="AnalysisResultsRunCutsLC" # CutsLC
     648bin=-1
     649# Mrk 421
     650source=1
     651name="Mrk421_all_nightly"
     652get_results
     653# Mrk 501
     654source=2
     655name="Mrk501_all_nightly"
     656get_results
     657# 1959
     658source=7
     659name="1959_all_nightly"
     660get_results
     661# for other INTEGRAL proposal
     662# Mrk 421
     663bin=0
     664source=1
     665name="Mrk421_all_period"
     666get_results
     667
     668exit
     669
     670
     671# Crab
     672source=5
     673nightmin=20121214
     674nightmax=20180418
     675table="AnalysisResultsRunCutsLC" # CutsLC
     676bin=-1
     677zdmax=30
     678thmax=330
     679name="testCrabnightlycheck"
     680get_results
     681
     682exit
     683
     684# Crab
     685source=5
     686nightmin=20121214
     687nightmax=20180418
     688table="AnalysisResultsRunCutsLC" # CutsLC
     689bin=20
     690name="testCrab20min"
     691get_results
     692
     693exit
     694
     695# Crab
     696source=5
     697nightmin=20121214
     698nightmax=20180418
     699table="AnalysisResultsRunCutsLC" # CutsLC
     700bin=-1
     701name="testCrabnightly"
     702get_results
     703
     704exit
     705
     706# Crab
     707source=5
     708#nightmin=20121214
     709#nightmax=20180418
     710table="AnalysisResultsRunCutsLC" # CutsLC
     711bin=-1
     712name="testCraball"
     713get_results
     714
     715exit
     716
     717
     718# Crab
     719source=5
     720nightmin=20121214
     721nightmax=20180418
     722table="AnalysisResultsRunCutsLC" # CutsLC
     723bin=60
     724name="testCrab"
     725get_results
     726
     727exit
     728
     729
     730# Mrk 501
     731source=2
     732zdmax=90
     733thmax=1500
     734# new analysis
     735table="AnalysisResultsRunCutsLC" # CutsLC
     736nightmin=20140622
     737nightmax=20140624
     738name="Mrk501_5min_Flare_newAnalysis_forHESS"
     739bin=5
     740get_results
     741
     742exit
     743
    641744# for 2344 paper with MAGIC
    642745
Note: See TracChangeset for help on using the changeset viewer.