casa5
Public
Commits
| 991 991 | # eta_a = eta_p*eta_s*eta_b*eta_t |
| 992 992 | |
| 993 993 | # correlator quantization efficiency. |
| 994 994 | eta_q = cq[iobs] |
| 995 995 | |
| 996 996 | # Receiver radiation temperature in K. |
| 997 997 | if telescope=='ALMA' or telescope=='ACA' or telescope=='ALMASD': |
| 998 998 | # ALMA-40.00.00.00-001-A-SPE.pdf |
| 999 999 | # http://www.eso.org/sci/facilities/alma/system/frontend/ |
| 1000 1000 | |
| 1001 - | # limits instead of centers, go to higher band in gaps |
| 1002 - | f0=[31.3,45,84,116,163,211,275,373,500,720] |
| 1001 + | # lower limits |
| 1002 + | # f0=[ 31, 67, 84, 125, 162, 211, 275, 385, 602, 787, 950] |
| 1003 + | # go to higher band in gaps |
| 1004 + | f0=[ 31, 45, 84, 116, 162, 211, 275, 373, 500, 720, 950] |
| 1005 + | # 80% spec |
| 1006 + | # t0=[ 17, 30, 37, 51, 65, 83, 147, 196, 175, 230] |
| 1003 1007 | # cycle 1 OT values 7/12 |
| 1004 - | t0=[ 17, 30, 45, 51, 65, 55, 75, 196, 100, 230] |
| 1008 + | t0=[ 17, 30, 45, 51, 65, 55, 75, 196, 100, 230] |
| 1005 1009 | |
| 1006 1010 | flim=[31.3,950] |
| 1007 1011 | if self.verbose: self.msg("using ALMA/ACA Rx specs",origin="noisetemp") |
| 1008 1012 | else: |
| 1009 1013 | if telescope=='EVLA': |
| 1010 1014 | # 201009114 from rick perley: |
| 1011 1015 | # f0=[1.5,3,6,10,15,23,33,45] |
| 1012 1016 | t0=[10.,15,12,15,10,12,15,28] |
| 1013 1017 | # limits |
| 1014 1018 | f0=[1,2,4,8,12,18,26.5,40,50] |
| 1037 1041 | flim=[180.,720] |
| 1038 1042 | else: |
| 1039 1043 | self.msg("I don't know about the "+telescope+" receivers, using 200K",priority="warn",origin="noisetemp") |
| 1040 1044 | f0=[10,900] |
| 1041 1045 | t0=[200,200] |
| 1042 1046 | flim=[0,5000] |
| 1043 1047 | |
| 1044 1048 | obsfreq=freq_ghz.get("value") |
| 1045 1049 | # z=pl.where(abs(obsfreq-pl.array(f0)) == min(abs(obsfreq-pl.array(f0)))) |
| 1046 1050 | # t_rx=t0[z[0]] |
| 1047 - | z=0 |
| 1048 - | while(f0[z]<obsfreq and z<len(t0)): |
| 1049 - | z+=1 |
| 1050 - | t_rx=t0[z-1] |
| 1051 1051 | |
| 1052 1052 | if obsfreq<flim[0]: |
| 1053 + | t_rx=t0[0] |
| 1053 1054 | self.msg("observing freqency is lower than expected for "+telescope,priority="warn",origin="noise") |
| 1054 1055 | self.msg("proceeding with extrapolated receiver temp="+str(t_rx),priority="warn",origin="noise") |
| 1056 + | else: |
| 1057 + | z=0 |
| 1058 + | while(f0[z]<obsfreq and z<len(t0)): |
| 1059 + | z+=1 |
| 1060 + | t_rx=t0[z-1] |
| 1061 + | |
| 1055 1062 | if obsfreq>flim[1]: |
| 1056 1063 | self.msg("observing freqency is higher than expected for "+telescope,priority="warn",origin="noise") |
| 1057 1064 | self.msg("proceeding with extrapolated receiver temp="+str(t_rx),priority="warn",origin="noise") |
| 1058 1065 | if obsfreq<=flim[1] and obsfreq>=flim[0]: |
| 1059 1066 | self.msg("interpolated receiver temp="+str(t_rx),origin="noise") |
| 1060 1067 | |
| 1061 1068 | return eta_p, eta_s, eta_b, eta_t, eta_q, t_rx |
| 1062 1069 | |
| 1063 1070 | |
| 1064 1071 | |