Changeset 3792 in Sophya for trunk/Cosmo/RadioBeam/calcpk2.cc


Ignore:
Timestamp:
Jun 28, 2010, 6:06:00 PM (15 years ago)
Author:
ansari
Message:

Ajout du programme de calcul et sauvegarde reponse interfero, Reza 28/06/2010

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Cosmo/RadioBeam/calcpk2.cc

    r3789 r3792  
    6565   
    6666    bool fgcorrbeam=true;
    67     double tbeamarcmin=15.;
     67    double tbeamDoL=135;
    6868    if (narg>7) {
    69       tbeamarcmin=atof(arg[7]);
    70       if (tbeamarcmin<1.e-6)  fgcorrbeam=false;
     69      tbeamDoL=atof(arg[7]);
     70      if (tbeamDoL<1.)  fgcorrbeam=false;
    7171    }
    7272    bool fgclnsrc=true;
     
    120120    double lambda = conv.getLambda();
    121121    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.;
    123124    Four2DResponse tbeam(2, DoL, DoL );
    124125    ForegroundCleaner  cleaner(arep, tbeam, skycube);
    125126    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;
    128129      cleaner.BeamCorrections();
    129130    }
Note: See TracChangeset for help on using the changeset viewer.