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

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:

• trunk/Cosmo/RadioBeam/mdish.cc (modified) (3 diffs)

trunk/Cosmo/RadioBeam/mdish.cc

@@ -59 +59 @@
 // -- Four2DRespTable : Reponse tabulee instrumentale ds le plan k_x,k_y (angles theta,phi) 
 //---------------------------------------------------------------
-Four2DRespTable::Four2DRespTable(Histo2D const & hrep, double d)
-  : Four2DResponse(0,d,d) , hrep_(hrep)
+Four2DRespTable::Four2DRespTable()
+  : Four2DResponse(0,1.,1.)
+{
+}
+
+Four2DRespTable::Four2DRespTable(Histo2D const & hrep, double d, double lambda)
+  : Four2DResponse(0,d,d,lambda) , hrep_(hrep)
 {
 }
@@ -75 +80 @@
 }
 
+void Four2DRespTable::writeToPPF(string flnm)
+{
+  DVList dvinfo;
+  dvinfo["DoL"] = dx_;
+  dvinfo["LambdaRef"] = lambdaref_;
+  dvinfo["Lambda"] = lambda_;
+  POutPersist po(flnm);
+  po << hrep_;
+  po << dvinfo;
+}
+
+void Four2DRespTable::readFromPPF(string flnm)
+{
+  PInPersist pin(flnm);
+  DVList dvinfo;
+  pin >> hrep_;
+  pin >> dvinfo;
+  dx_ = dy_ = dvinfo["DoL"];
+  setLambdaRef((double)dvinfo["LambdaRef"]);
+  setLambda((double)dvinfo["Lambda"]); 
+}
+
+
+
 //---------------------------------------------------------------
 // -- Four2DRespRatio : rapport de la reponse entre deux objets Four2DResponse
@@ -87 +116 @@
   double ra = a_.Value(kx,ky);
   double rb = b_.Value(kx,ky);
-  if (rb<divzthr_) {
-    if (ra<rb)  return 0.;
-    else rb=divzthr_;
-  }
+  if (ra<rb) {
+    if (rb>1.e-39)  return(ra/rb);  
+    else return 0.;
+  }
+  if (rb<divzthr_)  rb=divzthr_;
   return (ra/rb);
 }