Changeset 1735 in Sophya for trunk/ArchTOIPipe/ProcWSophya

Timestamp:

Nov 6, 2001, 7:01:46 PM (24 years ago)

Author:

cmv

Message:

add FlagAround + enum FlagToiAll cmv 6/11/2001

Location:

trunk/ArchTOIPipe/ProcWSophya

Files:

• ktoibad.cc (modified) (2 diffs)
• ktoibad.h (modified) (2 diffs)

trunk/ArchTOIPipe/ProcWSophya/ktoibad.cc

@@ -104 +104 @@
 
   getData(0,sn,bolog,fgbolog);
-  if(fgbolog && flgNotLookAt) {
+  if(fgbolog & flgNotLookAt) {
     putData(0,sn,bolog,fgbolog);
     continue;
@@ -218 +218 @@
 return;                                                                            
 }
+
+
+
+
+
+
+
+////////////////////////////////////////////////////////////////////////
+// Flag des sample autour d'un sample deja flaggue :
+// INPUT:
+// lm : longueur a flagguer de part et d'autre du mauvais sample
+//      <---lm---><bad_sample><---lm--->
+////////////////////////////////////////////////////////////////////////
+
+////////////////////////////////////////////////////////////////////////
+FlagAround::FlagAround(uint_4 lm,uint_8 flgs)
+{
+ // Mise en forme des arguments
+ LSide = lm;
+
+ SetFlag();
+ SetFlagAroundFlag(flgs);
+}
+
+FlagAround::~FlagAround()
+{
+}
+
+////////////////////////////////////////////////////////////////////////
+void FlagAround::Print(ostream & os)
+{
+ os<<"FlagAround::Print"
+   <<" LSide="<<LSide
+   <<" flgSample="<<flgSample<<" flgAround="<<flgAround
+   <<endl;
+}
+
+////////////////////////////////////////////////////////////////////////
+void FlagAround::init() {
+  cout << "FlagAround::init" << endl;
+  declareInput("BoloIn");       // input index 0
+  declareOutput("BoloOut");     // output index 0
+
+  int_4 lbuff = 2*LSide +10;
+  if(neededHistory<lbuff) setNeededHistory(lbuff);
+}
+
+////////////////////////////////////////////////////////////////////////
+void FlagAround::run()
+{
+long snb = getMinIn();
+long sne = getMaxIn();
+
+if(snb>sne) {
+  cout<<"FlagAround::run() - Bad sample interval"<<snb<<" , "<<sne<<endl;
+  throw ParmError("FlagAround::run() - Bad sample interval");
+}
+if(!checkInputTOIIndex(0)) {
+  cout<<"FlagAround::run() - Input TOI (BoloIn) not connected! "<<endl;
+  throw ParmError("FlagAround::run() Input TOI (BoloIn) not connected!");
+}
+if(!checkOutputTOIIndex(0)) {
+  cout<<"FlagAround::run() - Output TOI (BoloOut) not connected! "<<endl;
+  throw ParmError("FlagAround::run() Output TOI (BoloOut) not connected!");
+}
+
+//---------------------------------------------------------
+try {
+
+cout<<"FlagAround::run(): SampleNum de "<<snb<<" a "<<sne<<endl;
+uint_4 nflag=0, nflagged=0, nflagalready=0;
+
+// Un peu boeuf! On pourrait optimiser!
+// Si on ne tue pas (LSide=0), alors juste un transfert de donnees
+for(int_4 sn=snb;sn<=sne;sn++) {
+  r_8 bolo; uint_8 fgbolo=0;
+  getData(0,sn,bolo,fgbolo);
+  if(fgbolo&flgSample) { // Pixel flaggue flgSample
+    nflag++;
+  } else if(fgbolo&flgAround) { // Pixel deja flagge around?
+    nflagalready++;
+  } else if(LSide>0) { // Pixel a flagguer?
+    int_4 ideb = sn-LSide, ifin = sn+LSide;
+    for(int_4 i=ideb;i<=ifin;i++) { // Y a t-il un "mauvais" pixel autour?
+      if(i<snb || i>sne || i==sn) continue;
+      r_8 b; uint_8 fg=0;
+      getData(0,i,b,fg);
+      if(fg&flgSample) {fgbolo |= flgAround; nflagged++; break;}
+    }
+  }
+  putData(0,sn,bolo,fgbolo);
+}
+cout<<"FlagAround::run(): Fin de boucle: nflag="<<nflag
+    <<" nflagged="<<nflagged<<" (already="<<nflagalready
+    <<") / tot="<<sne-snb+1<<endl;
+
+//---------------------------------------------------------
+} catch (PException & exc) {
+  cout<<"FlagAround: Catched Exception "<<(string)typeid(exc).name()
+      <<"\n .... Msg= "<<exc.Msg()<<endl;
+}
+
+return;                                                                            
+}

trunk/ArchTOIPipe/ProcWSophya/ktoibad.h

@@ -9 +9 @@
 //-- Pour rechercher les glitches
 //   toi BoloIn ---> |FlagGlitch | ---> toi BoloOut
-
 class FlagGlitch : public TOIProcessor {
 public:
@@ -43 +42 @@
 };
 
+//-- Pour flagguer les samples autour de samples deja flaggues
+//   toi BoloIn ---> |FlagAround | ---> toi BoloOut
+class FlagAround : public TOIProcessor {
+public:
+  FlagAround(uint_4 lm,uint_8 flgs=FlgToiAll);
+  virtual ~FlagAround();
+
+  virtual void  init(void);
+  virtual void  run(void);
+
+  //! Flagguer les samples avec la valeur flag
+  inline void SetFlag(uint_8 flag=FlgToiAround) {flgAround = flag;}
+  //! Flagguer autour des samples ayant la valeur flag
+  inline void SetFlagAroundFlag(uint_8 flag=FlgToiAll) {flgSample = flag;}
+  //! Print
+  void Print(ostream & os);
+  inline void Print(void) {Print(cout);}
+
+protected:
+  uint_8 flgSample,flgAround;
+  uint_4 LSide;
+};
+
 #endif