%%%%%%%%%%%%%% S4 performance given SPT %%%%%%%%%%%%%%%

% BK achieved survey weight per rx-year, det count per rx-year, and per-detector NET
wbk=[24000, 30000, 2000];
ndetbk=[288, 512, 512];
netbk=[287.6, 313.1, 837.1];
% BK survey weight per det-year
wbkdet = wbk./ndetbk;

% S4 HR frequencies and per-detector NET
nu=[20, 30, 40, 95, 145, 220, 270];
nets4=[214, 177, 224, 238, 309, 747, 1281];

% SPT performance
n145=1.3; % calculated for 25k*4 det-yrs
% Number of detectors for S4 HR
ndets4=[130,250,500,25000,25000,8700,8700].*4;

% S4 performance for cross check
%n145=2.19;
%ndets4=[533,900,1600,73200,63075,118130,186300];

% Calculate S4 performance at the other bands using the performance and
% detector count in the band specified above.
%
% 20, 30, 40, 95 are scaled from BK95
% 145 is scaled from BK150
% 220, 270 are scaled from BK220
for i = 1:length(nu)
  if ismember(nu(i),nu(1:4))
    n(i) = n145 .* sqrt(wbkdet(2)/wbkdet(1)) .* (netbk(2)/netbk(1)) .* (nets4(i)/nets4(5)) .* sqrt(ndets4(5)/ndets4(i));
  elseif ismember(nu(i),nu(5))
    n(i) = n145 .* sqrt(wbkdet(2)/wbkdet(2)) .* (netbk(2)/netbk(2)) .* (nets4(i)/nets4(5)) .* sqrt(ndets4(5)/ndets4(i));
  else
    n(i) = n145 .* sqrt(wbkdet(2)/wbkdet(3)) .* (netbk(2)/netbk(3)) .* (nets4(i)/nets4(5)) .* sqrt(ndets4(5)/ndets4(i));
  end
end