Changeset 2047 in Sophya

Timestamp:

Jun 10, 2002, 10:03:24 AM (23 years ago)

Author:

ansari

Message:

Ajout programme toistat.cc - TOI 2 StatNTuple(Mean-sigma) - Reza 10/6/2002

Location:

trunk/ArchTOIPipe/TestPipes

Files:

• SMakefile (modified) (2 diffs)
• aksj02.cc (modified) (6 diffs)
• toistat.cc (added)

trunk/ArchTOIPipe/TestPipes/SMakefile

@@ -32 +32 @@
 endif
 ifeq ($(MACHEROS),Linux)
-LIBS := $(LIBS) -ldl -lf2c
+LIBS := $(LIBS) -ldl 
 endif
 
@@ -67 +67 @@
         ../cleantstpipes.csh
 
-all :  mesovh mesovh2 tsttoi2map simtst simofftst aksj02 tstmap2toi
+all :  mesovh mesovh2 tsttoi2map simtst simofftst aksj02 tstmap2toi toistat

trunk/ArchTOIPipe/TestPipes/aksj02.cc

@@ -42 +42 @@
          << "         [-dbg] [-wtoi sz] [-wclean wsz,nbw] [-range min,max] \n"
          << "         [-soff wsz,nptfit,degpol] [-soffnt PPFName] \n"
-         << "         [-bgalcut bmin,bmax] [-bgal toiname] [-bgalfile pointFitsName] \n"
          << "         [-soff2 wsz,nptfit,degpol] [-soff2nt PPFName] \n"
-         << "         [-bgal2cut bmin,bmax] \n"
+         << "         [-bgal toiname] [-bgalfile pointFitsName] [-bgalnf N]\n"
+         << "         [-bgal toiname] [-bgalfile pointFitsName] \n"
+         << "         [-bgalcut bmin,bmax]  [-bgal2cut bmin,bmax] \n"
          << "         [-gfilt wsz,sigma] \n" 
          << "         [-nooutflg] [-wrtms] [-wrtifo] [-prstat] [-useseqbuff] \n"
@@ -61 +62 @@
          << "   -bgal toiname: Galactic latitude toiname \n"
          << "   -bgalfile pointFitsName: Pointing FITS file name (for gal. latitude) \n"
+         << "   -bgalnf n: Number of bgalfiles pointFitsName%d.fits \n"
          << "   -soff2 wsz,nptfit,degpol : set SimpleOffsetEstimator_2 parameters\n"
          << "   -soff2nt PPFName : Writes out SimpleOffsetEstimator_2 NTuple \n"
@@ -106 +108 @@
 
   // cleaner parameters
-  double clean_wsz = 256;
-  double clean_nbw = 5;
+  int clean_wsz = 256;
+  int clean_nbw = 5;
 
   // SimpleOffsetEstimator parameters
@@ -139 +141 @@
   double bmax2 = 10.;
   bool bgal2cut = false;
+  int bgalNbFiles = 0;
 
   // File names 
@@ -222 +225 @@
       ia++;
     }
+    else if (strcmp(arg[ia],"-bgalnf") == 0) {
+      if (ia == narg-1) Usage(true);
+      bgalNbFiles = atoi(arg[ia+1]);
+      ia++;
+    }
     else if (strcmp(arg[ia],"-bgal") == 0) {
       if (ia == narg-1) Usage(true);
@@ -307 +315 @@
     FITSTOIReader* rbgal = NULL;
     if (bgalfile) {  // if Galactic cut 
-      cout << "> Creating bgal FITSTOIReader object - InFile=" << pointfile 
+      string pointfileI = pointfile;
+      MuTyV numf=0;
+      if(bgalNbFiles > 0) {
+        pointfileI += (string)numf + ".fits";
+      }
+      cout << "> Creating bgal FITSTOIReader object - InFile=" << pointfileI
            << " bgaltoiname= " << bgaltoi << endl;
-      rbgal = new FITSTOIReader(pointfile); 
+      rbgal = new FITSTOIReader(pointfileI);
+      if(bgalNbFiles > 0) 
+        for(int kff=1; kff<bgalNbFiles; kff++) {
+          numf=kff;
+          pointfileI = pointfile;
+          pointfileI += (string)numf + ".fits";
+          cout << " bgal FITSTOIReader.addFile( " << pointfileI << " )" << endl;
+          rbgal->addFile( pointfileI );
+        }
+      
     }