Changeset 2014 in Sophya for trunk/ArchTOIPipe/TestPipes/simofftst.cc

Timestamp:

May 28, 2002, 7:30:45 PM (23 years ago)

Author:

ansari

Message:

Amelioration de calcul d'offset (possibilite de coupure sur coordonnees
galactique)
Nouveau processeur nettoyeur (SimpleCleaner)
Ajout programme de traitement aksj02.cc , donnees Archeops KS1/KS3

Reza 28/5/2002

File:

• trunk/ArchTOIPipe/TestPipes/simofftst.cc (modified) (11 diffs)

trunk/ArchTOIPipe/TestPipes/simofftst.cc

@@ -3 +3 @@
 //                               Christophe Magneville
 //                               Reza Ansari
-// $Id: simofftst.cc,v 1.2 2002-05-16 13:13:00 ansari Exp $
+// $Id: simofftst.cc,v 1.3 2002-05-28 17:30:45 ansari Exp $
 
 /*   Test de processeurs ds simtoipr.cc   - Reza Avril 2001
@@ -11 +11 @@
             -intoi boloMuV_26 -wtoi 8192 -wdegli 1024 \
             inputbolo.fits filt.fits xx.ppf
-# Avec Filtre de Fourier
+# Avec Filtre 
 csh> simofftst -start 104385384 -end 104399964 -range -500,500 \
             -intoi boloMuV_26 -wtoi 8192 -wdegli 1024 \
-            -wfft 4096 -keepfft 5 -killfreq 15,4,2.5 \ 
             inputbolo.fits filtfft.fits spectre.ppf
-# Autre exemple b545k02
-cool> ./simofftst -start 104389122 -end 104649122 -range -1000.,1000. -intoi boloMu
-V_23 -wtoi 8192 -wdegli 1024 -wfft 4096 -keepfft 5 -killfreq 15,4,2.5 Cmv/b545k02
-2_16.00_16.49.fits b545.fits b545.ppf
+csh> $exedir/simofftst -start $sn0 -end $sn1  -intoi $cleanname \
+            -wtoi 8192 -wdegli 1024 -degli 3.,2.,3,2,0 \
+            -soff 128,9,3 -gfilt 16,2.5 -prstat \
+            -soffnt $outppf -bgalcut -20.,20. -bgal $glat \
+            -bgalfile $point $cleanfile $outfile
 
 */
@@ -38 +38 @@
 #include "timing.h"
 #include "histinit.h"
+#include "psighand.h"
 #include <stdexcept>
 
@@ -52 +53 @@
          << "         [-degli ns,ns2,mxpt,mnpt,wrec] [-soff wsz,nptfit,degpol] \n"
          << "         [-soffnt PPFName] [-gfilt wsz,sigma] \n"
+         << "         [-bgalcut bmin,bmax] [-bgal toiname] [-bgalfile pointFitsName] \n"
          << "         [-nooutflg] [-nowrtms] [-nowrticd] [-prstat] [-useseqbuff] \n"
          << "         inFitsName outFitsName \n" 
@@ -77 +79 @@
          << endl;
   }
+  if (fgerr) exit(1);
 }
 
@@ -124 +127 @@
   double gfilt_sigma = 2.;
 
+  // Coordinate file for galactic cut
+  string pointfile = "pointgal.fits";
+  double bmin = -20;
+  double bmax = 20.;
+  string bgaltoi = "bgal";
+  bool bgalcut = false;
+
   // File names 
   string infile;
@@ -183 +193 @@
       ia++;
     }    
+    else if (strcmp(arg[ia],"-bgalcut") == 0) {
+      if (ia == narg-1) Usage(true);
+      sscanf(arg[ia+1],"%lf,%lf", &bmin, &bmax);
+      bgalcut = true;
+      ia++;
+    }
+    else if (strcmp(arg[ia],"-bgalfile") == 0) {
+      if (ia == narg-1) Usage(true);
+      pointfile = arg[ia+1];
+      ia++;
+    }
+    else if (strcmp(arg[ia],"-bgal") == 0) {
+      if (ia == narg-1) Usage(true);
+      bgaltoi = arg[ia+1];
+      ia++;
+    }
     else if (strcmp(arg[ia],"-wfft") == 0) {
       if (ia == narg-1) Usage(true);  
@@ -226 +252 @@
   cout << " Initializing SOPHYA ... " << endl;
   SophyaInit();
+  SophyaConfigureSignalhandling(true);
+
   InitTim();
 
@@ -263 +291 @@
     double G_a = 1./(G_sigma*sqrt(M_PI*2.));
     SimpleFilter filt(gfilt_wsz, SimpleFilter::GaussFilter, G_a, G_sigma);
-    
+
+         
+    FITSTOIReader* rbgal = NULL;
+    if (bgalcut) {  // if Galactic cut 
+      cout << "> Creating bgal FITSTOIReader object - InFile=" << pointfile 
+           << " bgaltoiname= " << bgaltoi << endl;
+      cout << "  offes.SetBGalCut( " << bmin << "," << bmax << ")" << endl;
+      rbgal = new FITSTOIReader(pointfile); 
+      offes.SetBGalCut(bmin, bmax);      
+    }
+
     cout << "> Creating FITSTOIWriter OutFitsName= " << outfile << endl;
     FITSTOIWriter w(outfile);
@@ -281 +319 @@
       plombier.Connect(offes, "out", w, "degoff", "", 0, fg_wrtflag);
       
+    if (bgalcut) {
+      inname = "bgal";
+      plombier.Connect(*rbgal, bgaltoi, offes, inname);      
+    }
 
     if (fg_wrticd) {
@@ -315 +357 @@
     mgr->joinAll();
     PrtTim("End threads");
+    if (bgalcut)  delete rbgal;
 
   }