sourceCorrespondance

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  
  function sourceNum = sourceCorrespondance(sourceName, mjd)

   this function is just a simple correspondance between a source name and
   anumber, for posterity, just to make saving things as a text file
   easier.

  sjcm
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

0001 function [sourceNum sourceFlux fluxErr polang polfrac isoff ] = sourceCorrespondance(sourceName, mjd)
0002 
0003 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0004 %
0005 %  function sourceNum = sourceCorrespondance(sourceName, mjd)
0006 %
0007 %   this function is just a simple correspondance between a source name and
0008 %   anumber, for posterity, just to make saving things as a text file
0009 %   easier.
0010 %
0011 %  sjcm
0012 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0013 
0014 if(nargin<2)
0015   year = 2012;
0016 else
0017   year = mjd2fracyear(mjd);
0018 end
0019 %[casS casSerr cygS cygSerr tauS tauSerr] = CalibFluxDen(5.0, year);
0020 
0021 s1 = 'CasA';
0022 s2 = 'TauA';
0023 s3 = 'CygA';
0024 
0025 [casS, casSerr] = calcFlux(s1, year);
0026 [tauS, tauSerr] = calcFlux(s2, year);
0027 [cygS, cygSerr] = calcFlux(s3, year);
0028 
0029 
0030 switch(sourceName{1})
0031   case 'CygA'
0032     sourceNum = 1;
0033     sourceFlux = cygS;
0034     fluxErr    = cygSerr; %
0035 %    sourceFlux = 365;
0036 %    fluxErr    = 0.02; %
0037     
0038   case 'cygaOff'
0039     sourceNum = 1;
0040     sourceFlux = cygS;
0041     fluxErr    = cygSerr; %
0042     
0043   case 'TauA'
0044     sourceNum = 2;
0045     sourceFlux = tauS;
0046     fluxErr    = tauSerr;
0047     polang = 38;% Need to get an accurate number for this
0048     polfrac = 0.05;% Need to get an accurate number for this
0049 %    sourceFlux = 615;
0050 %    fluxErr    = 0.02; %
0051     
0052   case 'TauAoff'
0053     sourceNum = 2;
0054     sourceFlux = tauS;
0055     fluxErr    = tauSerr;
0056     
0057   case 'J0437+296' %3c123
0058     sourceNum = 3;
0059     sourceFlux = 16.56;
0060     fluxErr    = 0.03; %
0061     
0062   case 'J0915-119' %3C218
0063     sourceNum = 4;
0064     sourceFlux = 13.65;
0065     fluxErr    = 0.03; %
0066     
0067   case 'J1230+123' %3c274
0068     sourceNum = 5;
0069     sourceFlux = 68;
0070     fluxErr    = 2; %
0071     
0072   case '3C274' %3c274; virgoA
0073     sourceNum = 5;
0074     sourceFlux = 68;
0075     fluxErr    = 2; %
0076     
0077   case 'J1331+305' %3C286'
0078     sourceNum = 6;
0079     sourceFlux = 7.6;
0080     fluxErr    = 0.03; %
0081     
0082   case 'J1720-009' %3C353
0083     sourceNum = 7;
0084     sourceFlux = 23.49;
0085     fluxErr    = 0.03; %
0086     
0087   case '3C353' %3C353
0088     sourceNum = 7;
0089     sourceFlux = 23.49;
0090     fluxErr    = 0.03; %
0091     
0092   case 'CasA'
0093     sourceNum = 8;
0094     sourceFlux = casS;
0095     fluxErr    = casSerr;
0096 %    sourceFlux = 670;
0097 %    fluxErr    = 0.02; %
0098     
0099   case 'CasAoff'
0100     sourceNum = 8;
0101     sourceFlux = casS;
0102     fluxErr    = casSerr;
0103       
0104   case 'J0225+621'
0105     sourceNum = 9;
0106     sourceFlux = 135.5;  % my best estimate
0107     fluxErr    = 0.05; %
0108     
0109   case 'J0927+390'
0110     sourceNum = 10;
0111     sourceFlux = 8.5;
0112     fluxErr    = 0.50; %
0113     
0114   case 'M42'
0115     sourceNum = 11;
0116     sourceFlux = 480; % Need to get an accurate number for this
0117     fluxErr    = 0.01; %  Need to get an accurate number for this
0118     polang = 0;% Need to get an accurate number for this
0119     polfrac = 0;% Need to get an accurate number for this
0120     
0121   otherwise
0122     sourceNum = 99;
0123     sourceFlux = 0;
0124     fluxErr    = 0;
0125 end
0126 
0127 return;