Changeset 542

Timestamp:

Oct 3, 2011, 6:53:43 PM (13 years ago)

Author:

campagne

Message:

forbid INTERACTIVE Mode

Location:

BAORadio/AmasNancay

Files:

• anarawonoffmaker.sh (modified) (1 diff)
• calibmaker.sh (modified) (2 diffs)
• gainmaker.sh (modified) (1 diff)
• mergeRawOnOff.cc (modified) (7 diffs)
• onoffmaker.sh (modified) (1 diff)
• rawonoffmaker.sh (modified) (1 diff)
• submit2bqs-anarawonoffmaker.sh (modified) (1 diff)

BAORadio/AmasNancay/anarawonoffmaker.sh

@@ -82 +82 @@
     iojobpath="/sps/baoradio/AmasNancay/JEC"
     $MKDIR -p $iojobpath
+    $ECHO "Usage in INTERACTIVE may be destructive... To be used with great care !!!!"
+    exit 0;
 elif [ ${ENVIRONMENT} = "BATCH" ] ; then
     iojobpath=${TMPBATCH}

BAORadio/AmasNancay/calibmaker.sh

@@ -84 +84 @@
     iojobpath="/sps/baoradio/AmasNancay/JEC"
     $MKDIR -p $iojobpath
+    $ECHO "Usage in INTERACTIVE may be destructive... To be used with great care !!!!"
+    exit 0;
 elif [ ${ENVIRONMENT} = "BATCH" ] ; then
     iojobpath=${TMPBATCH}
@@ -140 +142 @@
   cd ${topDir}
   $LS | $XARGS -i $RM -rf {}
+
   
   $ECHO "DEBUG: (${scriptName}) irods root dir $i"

BAORadio/AmasNancay/gainmaker.sh

@@ -82 +82 @@
     iojobpath="/sps/baoradio/AmasNancay/JEC"
     $MKDIR -p $iojobpath
+    $ECHO "Usage in INTERACTIVE may be destructive... To be used with great care !!!!"
+    exit 0;   
 elif [ ${ENVIRONMENT} = "BATCH" ] ; then
     iojobpath=${TMPBATCH}

BAORadio/AmasNancay/mergeRawOnOff.cc

@@ -37 +37 @@
 #include "ppersist.h"
 
+//-----------------------------------------------
+const sa_size_t NUMBER_OF_CHANNELS = 2;
+const sa_size_t  NUMBER_OF_FREQ = 8192;
+const r_4    LOWER_FREQUENCY = 1250.0; //MHz
+const r_4    TOTAL_BANDWIDTH = 250.0; //MHz
 //-----------------------------------------------
 char *regexp (const char *string, const char *patrn, int *begin, int *end) {    
@@ -61 +66 @@
 
 struct Param {
-  string debuglev_;
+  int debuglev_;
   string inPath_;
   string outPath_;
@@ -154 +159 @@
   directoryName = para.inPath_ + "/" + para.sourceName_;
 
+  //Make the listing of the directory
   listOfFiles = ListOfFileInDir(directoryName,para.ppfFile_);
   
@@ -159 +165 @@
   iFileEnd = listOfFiles.end();
   
-  StringMatch match("specONOFFRaw[0-9]+");
+  StringMatch match("specONOFFRaw[0-9]+"); //Tag of the PPF objects
+  TMatrix<r_4> sumOfSpectra(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ);
+  uint_4 nSpectra=0;
+  //Loop on files
   for (iFile = listOfFiles.begin(); iFile != iFileEnd; ++iFile) {
-    cout << "load file <" << *iFile << ">" << endl;
+    if (para.debuglev_>90){
+      cout << "load file <" << *iFile << ">" << endl;
+    }
     PInPersist fin(*iFile);
     vector<string> vec = fin.GetNameTags();
     list<string> listOfSpectra;
+    //Keep only required PPF objects
     std::remove_copy_if(
                         vec.begin(), vec.end(), back_inserter(listOfSpectra),
@@ -172 +184 @@
     iSpecEnd = listOfSpectra.end();
     listOfSpectra.sort(stringCompare);
-
+    //Loop of spectra matrix
     for (iSpec = listOfSpectra.begin(); iSpec !=iSpecEnd;  ++iSpec){
-      cout << " spactra <" << *iSpec << ">" << endl;
-    }
-
-  }
-  
+      if (para.debuglev_>90){
+        cout << " spactra <" << *iSpec << ">" << endl;
+      }
+      TMatrix<r_4> aSpec(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ);
+      fin.GetObject(aSpec,*iSpec);
+      //How to see if the GetObject is ok?? Ask Reza
+      nSpectra++;
+      sumOfSpectra+=aSpec;
+    }//eo loop on spectra in a file
+  }//eo loop on files
+  
+  //Normalisation
+  if(nSpectra>0)sumOfSpectra/=(r_4)(nSpectra);
+
+  {//Save the result
+    stringstream tmp;
+    tmp << nSpectra;
+    string fileName = para.outPath_+"/meanDiffOnOffRaw_"+StringToLower(para.sourceName_)+"-"+tmp.str()+"Cycles.ppf";
+    cout << "Save mean based on " <<  nSpectra << " cycles " << endl;
+    string tag = "mean";
+    POutPersist fos(fileName);
+    fos << PPFNameTag(tag) << sumOfSpectra;
+  }
 }
 //-------------------------------------------------------
@@ -226 +256 @@
   }//eo while
 
-  para.debuglev_ = debuglev;
+  para.debuglev_ = atoi(debuglev.c_str());
   para.inPath_   = inputPath;
   para.outPath_  = outputPath;
@@ -250 +280 @@
   try {
 
-    int b,e;
-    char *match=regexp("truc0machin","[a-z]+[0-9]*",&b,&e);
-    printf("->%s<-\n(b=%d e=%d)\n",match,b,e);
+//     int b,e;
+//     char *match=regexp("truc0machin","[a-z]+[0-9]*",&b,&e);
+//     printf("->%s<-\n(b=%d e=%d)\n",match,b,e);
 
     if ( action == "default" ) {

BAORadio/AmasNancay/onoffmaker.sh

@@ -84 +84 @@
     iojobpath="/sps/baoradio/AmasNancay/JEC"
     $MKDIR -p $iojobpath
+    $ECHO "Usage in INTERACTIVE may be destructive... To be used with great care !!!!"
+    exit 0;
 elif [ ${ENVIRONMENT} = "BATCH" ] ; then
     iojobpath=${TMPBATCH}

BAORadio/AmasNancay/rawonoffmaker.sh

@@ -82 +82 @@
     iojobpath="/sps/baoradio/AmasNancay/JEC"
     $MKDIR -p $iojobpath
+    $ECHO "Usage in INTERACTIVE may be destructive... To be used with great care !!!!"
+    exit 0;
 elif [ ${ENVIRONMENT} = "BATCH" ] ; then
     iojobpath=${TMPBATCH}

BAORadio/AmasNancay/submit2bqs-anarawonoffmaker.sh

@@ -34 +34 @@
 
 
-jobBatchName="rawonoff-${sourceRadio}${dateJob}"
+jobBatchName="anarawonoff-${sourceRadio}${dateJob}"
 
 jobLogName="${jobBatchName}.log.$$"