Changeset 3792 in Sophya for trunk/Cosmo/RadioBeam/calcpk2.cc
- Timestamp:
- Jun 28, 2010, 6:06:00 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/Cosmo/RadioBeam/calcpk2.cc
r3789 r3792 65 65 66 66 bool fgcorrbeam=true; 67 double tbeam arcmin=15.;67 double tbeamDoL=135; 68 68 if (narg>7) { 69 tbeam arcmin=atof(arg[7]);70 if (tbeam arcmin<1.e-6) fgcorrbeam=false;69 tbeamDoL=atof(arg[7]); 70 if (tbeamDoL<1.) fgcorrbeam=false; 71 71 } 72 72 bool fgclnsrc=true; … … 120 120 double lambda = conv.getLambda(); 121 121 Four2DResponse arep(2, InterfArrayDiametre/lambda, InterfArrayDiametre/lambda, lambda); 122 double DoL = 1.22/ArcminToRadian(tbeamarcmin); 122 double DoL = tbeamDoL; 123 double tbeamarcmin = RadianToDegree(1.22/DoL)*60.; 123 124 Four2DResponse tbeam(2, DoL, DoL ); 124 125 ForegroundCleaner cleaner(arep, tbeam, skycube); 125 126 if (fgcorrbeam) { 126 cout << "calcpk2[3.a] : calling cleaner.BeamCorrections() for target beam (" << tbeamarcmin << " arcmin)"127 << " Diam/Lambda=" << DoL<< endl;127 cout << "calcpk2[3.a] : calling cleaner.BeamCorrections() for target beam D/Lambda=" << DoL 128 << " -> arcmin " << tbeamarcmin << endl; 128 129 cleaner.BeamCorrections(); 129 130 }
Note:
See TracChangeset
for help on using the changeset viewer.