Nid = 1;
RA= linspace(0,24,37);
DEC = linspace(-90,90,19);
for indRA=1:length(RA)
for indDEC = 1:length(DEC)
idlist(Nid,1:3) = [length(DEC)*(indRA-1)+jj,RA(indRA),DEC(indDEC)];
Nid = Nid + 1;
end
end
ID = 234;
RAs = idlist(ID,2);
DECs = idlist(ID,3);
jbname='20131202_SN_Rp1burst_5kpc_4096_injallsky';
fname= fullfile([jbname,'_results'],['hphctseratreq_',jbname,num2str(ID)])
load(fname);
gpstime=injection.injSiginfo.timestamps(1);
[thetaS, phiS] = convertcelest2earth(RAs, DECs, gpstime)
detectors = {'H1','L1','V1','LCGT'};
phiV=[-180:4:180];
thetaV=[-90:4:90];
[Fp,Fc,Ts] = genfpfcts(detectors,phiV,thetaV,'');
AntennaFp = zeros(length(thetaV),length(phiV));
AntennaFc = zeros(length(thetaV),length(phiV));
for idet = 1:length(detectors)
AntennaFp = Fp(:,:,idet).^2 + AntennaFp;
AntennaFc = Fc(:,:,idet).^2 + AntennaFc;
end
AntennaFp = sqrt(AntennaFp);
AntennaFc = sqrt(AntennaFc);
[val,idxphi] = min(abs(phiV-phiS));
[val,idxtheta] = min(abs(thetaV-thetaS));
AntennaFp(idxtheta, idxphi)
AntennaFc(idxtheta, idxphi)
More than 5 years have passed since last update.
Register as a new user and use Qiita more conveniently
- You get articles that match your needs
- You can efficiently read back useful information
- You can use dark theme