Changeset 929 in Sophya

Timestamp:

Apr 14, 2000, 3:35:55 PM (25 years ago)

Author:

ansari

Message:

Sophie: adapt to new class MapOperations of SkyMap

File:

• trunk/SophyaProg/PMixer/extractRS.cc (modified) (3 diffs)

trunk/SophyaProg/PMixer/extractRS.cc

@@ -1 @@
-#include "machdefs.h"
-#include <stdlib.h>
-#include <stdio.h>
-#include <fstream.h>
-#include <iostream.h>
-#include <math.h>
-
-#include <string>
-#include <vector>
-
-#include "tmatrix.h"
-#include "tvector.h"
-#include "fioarr.h"
-#include "vector3d.h"
-#include "srandgen.h"
-#include "timing.h"
-#include "sambainit.h"
-#include "pexceptions.h"
-#include "datacards.h"
-#include "fitsioserver.h"
-#include "fitsspherehealpix.h"
-#include "radspecvector.h"
-#include "blackbody.h"
-#include "nupower.h"
-
-#include "squarefilt.h"
-#include "trianglefilt.h"
-#include "specrespvector.h"
-#include "gaussfilt.h"
-#include "fitstarray.h"
-
+#include "pmixer.h"
 
 /*! \ingroup PMixer
-    \file extractRS.cc
- 
-  \brief  \b PROGRAM \b extractRS  <BR>
-  Program to create a map (SphereHEALPix) with point source
-  from small map (matrix) and convolving their radiation spectrum
-  RadSpectra
-  with a given filter response SpectralResponse
-*/
+ *  \file extractRS.cc 
+ * \brief  \b PROGRAM \b extractRS  <BR>
+ * Program to create a map (SphereHEALPix) with point source
+ * from small map (matrix) and convolving their radiation spectrum
+ * RadSpectra
+ * with a given filter response SpectralResponse
+ */
 
 // -----------------------------------------------------------------
@@ -46 +15 @@
 char * BuildFITSFileName(string const & fname);
 SpectralResponse * getSpectralResponse(DataCards & dc);
-template <class T>
-void  MeanSig(NDataBlock<T> const & dbl, double& gmoy, double& gsig);
 float ComputeFrequency(DataCards &,string&);
 
@@ -408 +375 @@
 
 
-/* Nouvelle-Fonction */
-template <class T>
-void  MeanSig(NDataBlock<T> const & dbl, double& gmoy, double& gsig)
-
-{
-  gmoy=0.;
-  gsig = 0.; 
-  double valok;
-  for(int k=0; k<(int)dbl.Size(); k++) {
-    valok = dbl(k);
-    gmoy += valok;  gsig += valok*valok; 
-  }
-  gmoy /= (double)dbl.Size();
-  gsig = gsig/(double)dbl.Size() - gmoy*gmoy;
-  if (gsig >= 0.) gsig = sqrt(gsig);
-}
 
 /* Nouvelle-Fonction */