| 1 | #include "machdefs.h"
|
|---|
| 2 | #include "toimanager.h"
|
|---|
| 3 | #include "pexceptions.h"
|
|---|
| 4 | #include "ntuple.h"
|
|---|
| 5 | #include "ktoibad.h"
|
|---|
| 6 |
|
|---|
| 7 |
|
|---|
| 8 | ////////////////////////////////////////////////////////////////////////
|
|---|
| 9 | // Flag des sample hors dynamique
|
|---|
| 10 | ////////////////////////////////////////////////////////////////////////
|
|---|
| 11 |
|
|---|
| 12 | FlagOutOfRange::FlagOutOfRange(r_8 vmin,r_8 vmax,uint_8 flag)
|
|---|
| 13 | : VMin(vmin), VMax(vmax), flgBad(flag)
|
|---|
| 14 | {
|
|---|
| 15 | if(VMax<VMin) {
|
|---|
| 16 | cout<<"FlagOutOfRange::run() - Bad range interval "<<VMin<<" , "<<VMax<<endl;
|
|---|
| 17 | throw ParmError("FlagOutOfRange::FlagOutOfRange() - Bad range interval");
|
|---|
| 18 | }
|
|---|
| 19 | }
|
|---|
| 20 |
|
|---|
| 21 | FlagOutOfRange::~FlagOutOfRange()
|
|---|
| 22 | {
|
|---|
| 23 | }
|
|---|
| 24 |
|
|---|
| 25 | void FlagOutOfRange::Print(::ostream & os)
|
|---|
| 26 | {
|
|---|
| 27 | os<<"FlagOutOfRange::Print Range: ["<<VMin<<","<<VMax<<"] flgBad="<<flgBad<<endl;
|
|---|
| 28 | }
|
|---|
| 29 |
|
|---|
| 30 | void FlagOutOfRange::init()
|
|---|
| 31 | {
|
|---|
| 32 | cout << "FlagOutOfRange::init" << endl;
|
|---|
| 33 | declareInput("DataIn"); // input index 0
|
|---|
| 34 | declareOutput("DataOut"); // output index 0
|
|---|
| 35 | }
|
|---|
| 36 |
|
|---|
| 37 | void FlagOutOfRange::run()
|
|---|
| 38 | {
|
|---|
| 39 | long snb = getMinIn();
|
|---|
| 40 | long sne = getMaxIn();
|
|---|
| 41 |
|
|---|
| 42 | if(snb>sne) {
|
|---|
| 43 | cout<<"FlagOutOfRange::run() - Bad sample interval"<<snb<<" , "<<sne<<endl;
|
|---|
| 44 | throw ParmError("FlagOutOfRange::run() - Bad sample interval");
|
|---|
| 45 | }
|
|---|
| 46 | if(!checkInputTOIIndex(0)) {
|
|---|
| 47 | cout<<"FlagOutOfRange::run() - Input TOI (DataIn) not connected! "<<endl;
|
|---|
| 48 | throw ParmError("FlagOutOfRange::run() Input TOI (DataIn) not connected!");
|
|---|
| 49 | }
|
|---|
| 50 | if(!checkOutputTOIIndex(0)) {
|
|---|
| 51 | cout<<"FlagOutOfRange::run() - Output TOI (DataOut) not connected! "<<endl;
|
|---|
| 52 | throw ParmError("FlagOutOfRange::run() Output TOI (DataOut) not connected!");
|
|---|
| 53 | }
|
|---|
| 54 |
|
|---|
| 55 | //---------------------------------------------------------
|
|---|
| 56 | try {
|
|---|
| 57 |
|
|---|
| 58 | cout<<"FlagOutOfRange::run(): SampleNum de "<<snb<<" a "<<sne<<endl;
|
|---|
| 59 | uint_4 nbad=0;
|
|---|
| 60 |
|
|---|
| 61 | for(int_4 sn=snb;sn<=sne;sn++) {
|
|---|
| 62 | r_8 bolo; uint_8 fgbolo=0;
|
|---|
| 63 | getData(0,sn,bolo,fgbolo);
|
|---|
| 64 | if(bolo<VMin || bolo>VMax) { // Pixel out of range
|
|---|
| 65 | nbad++;
|
|---|
| 66 | fgbolo |= flgBad;
|
|---|
| 67 | }
|
|---|
| 68 | putData(0,sn,bolo,fgbolo);
|
|---|
| 69 | }
|
|---|
| 70 | cout<<"FlagOutOfRange::run(): Fin de boucle: nbad="<<nbad
|
|---|
| 71 | <<" / tot="<<sne-snb+1<<endl;
|
|---|
| 72 |
|
|---|
| 73 | //---------------------------------------------------------
|
|---|
| 74 | } catch (PException & exc) {
|
|---|
| 75 | cout<<"FlagOutOfRange: Catched Exception "<<(string)typeid(exc).name()
|
|---|
| 76 | <<"\n .... Msg= "<<exc.Msg()<<endl;
|
|---|
| 77 | }
|
|---|
| 78 |
|
|---|
| 79 | return;
|
|---|
| 80 | }
|
|---|
| 81 |
|
|---|
| 82 |
|
|---|
| 83 | ////////////////////////////////////////////////////////////////////////
|
|---|
| 84 | // Recherche des glitches :
|
|---|
| 85 | // INPUT:
|
|---|
| 86 | // lg : longeur maxi d'un glitch (impair, si pair lg++)
|
|---|
| 87 | // doit etre suffisamment long pour que les glitches soient <lg
|
|---|
| 88 | // suffisamment court pour ne pas tuer la galaxie
|
|---|
| 89 | // lm : longueur de la fenetre glissante (de chaque cote)
|
|---|
| 90 | // <---lm---><--lg--><---lm--->
|
|---|
| 91 | // doit etre suffisamment long pour avoir une bonne stabilite
|
|---|
| 92 | // des calculs (a0,a1) et suffisamment court pour ne pas
|
|---|
| 93 | // tuer la galaxie.
|
|---|
| 94 | // nsg : nombre de sigmas pour valider un glitch
|
|---|
| 95 | ////////////////////////////////////////////////////////////////////////
|
|---|
| 96 |
|
|---|
| 97 | FlagGlitch::FlagGlitch(uint_4 lg,uint_4 lm,r_8 nsg)
|
|---|
| 98 | : mLP(0), mLPsnb(-1), mLPsne(-1)
|
|---|
| 99 | {
|
|---|
| 100 | // Mise en forme des arguments
|
|---|
| 101 | if(nsg<=0) nsg = 3.;
|
|---|
| 102 | if(lg<=0) lg=3; if(lg%2==0) lg++;
|
|---|
| 103 | if(lm<3) lm = 5*lg;
|
|---|
| 104 |
|
|---|
| 105 | LGlitch = lg;
|
|---|
| 106 | LSide = lm;
|
|---|
| 107 | NSCut = nsg;
|
|---|
| 108 | DoNotLookAt();
|
|---|
| 109 | SetFlag();
|
|---|
| 110 | SetBuffUpd();
|
|---|
| 111 | }
|
|---|
| 112 |
|
|---|
| 113 | FlagGlitch::~FlagGlitch()
|
|---|
| 114 | {
|
|---|
| 115 | }
|
|---|
| 116 |
|
|---|
| 117 | void FlagGlitch::Print(::ostream & os)
|
|---|
| 118 | {
|
|---|
| 119 | os<<"FlagGlitch::Print"
|
|---|
| 120 | <<" LGlitch="<<LGlitch<<" LSide="<<LSide<<" NSCut="<<NSCut
|
|---|
| 121 | <<" flgNotLookAt="<<flgNotLookAt<<" flgGlitch="<<flgGlitch
|
|---|
| 122 | <<" BuffUpdate="<<BuffUpdate<<endl;
|
|---|
| 123 | }
|
|---|
| 124 |
|
|---|
| 125 | ////////////////////////////////////////////////////////////////////////
|
|---|
| 126 | void FlagGlitch::init()
|
|---|
| 127 | {
|
|---|
| 128 | if(mLP) cout << "FlagGlitch::init" << endl;
|
|---|
| 129 | declareInput("DataIn"); // input index 0
|
|---|
| 130 | declareOutput("DataOut"); // output index 0
|
|---|
| 131 |
|
|---|
| 132 | int_4 lbuff = 2*LSide + LGlitch +10;
|
|---|
| 133 | if(neededHistory<lbuff) setNeededHistory(lbuff);
|
|---|
| 134 | }
|
|---|
| 135 |
|
|---|
| 136 | void FlagGlitch::run()
|
|---|
| 137 | {
|
|---|
| 138 | long snb = getMinIn();
|
|---|
| 139 | long sne = getMaxIn();
|
|---|
| 140 |
|
|---|
| 141 | if(snb>sne) {
|
|---|
| 142 | cout<<"FlagGlitch::run() - Bad sample interval"<<snb<<" , "<<sne<<endl;
|
|---|
| 143 | throw ParmError("FlagGlitch::run() - Bad sample interval");
|
|---|
| 144 | }
|
|---|
| 145 | if(!checkInputTOIIndex(0)) {
|
|---|
| 146 | cout<<"FlagGlitch::run() - Input TOI (DataIn) not connected! "<<endl;
|
|---|
| 147 | throw ParmError("FlagGlitch::run() Input TOI (DataIn) not connected!");
|
|---|
| 148 | }
|
|---|
| 149 | if(!checkOutputTOIIndex(0)) {
|
|---|
| 150 | cout<<"FlagGlitch::run() - Output TOI (DataOut) not connected! "<<endl;
|
|---|
| 151 | throw ParmError("FlagGlitch::run() Output TOI (DataOut) not connected!");
|
|---|
| 152 | }
|
|---|
| 153 |
|
|---|
| 154 | //---------------------------------------------------------
|
|---|
| 155 | try {
|
|---|
| 156 |
|
|---|
| 157 | SLinParBuff Slpb0(LSide,BuffUpdate,0.,0.,true);
|
|---|
| 158 | SLinParBuff Slpb1(LSide,BuffUpdate,0.,0.,true);
|
|---|
| 159 | SLinParBuff* Slpb[2] = {&Slpb0,&Slpb1};
|
|---|
| 160 |
|
|---|
| 161 | char* nament_[8] = {"s","bol","sg","sd","slg","sld","ypg","ypd"};
|
|---|
| 162 | r_4 xnt_[8];
|
|---|
| 163 | NTuple nt_(8,nament_);
|
|---|
| 164 | int intuple = mLP/10;
|
|---|
| 165 |
|
|---|
| 166 | if(mLP) cout<<"FlagGlitch::run(): SampleNum de "<<snb<<" a "<<sne<<endl;
|
|---|
| 167 | uint_4 nflag=0;
|
|---|
| 168 | r_8 nsg = NSCut;
|
|---|
| 169 | int_4 lg = LGlitch;
|
|---|
| 170 | int_4 lgs2 = (lg-1)/2;
|
|---|
| 171 | int_4 lm = LSide;
|
|---|
| 172 | int_4 lmmin = lm/2; if(lmmin<2) lmmin=2;
|
|---|
| 173 |
|
|---|
| 174 | // Initialisation
|
|---|
| 175 | if(LSide>0) for(int_4 sn=snb;sn<snb+LSide && sn<=sne;sn++) {
|
|---|
| 176 | r_8 bolog; uint_8 fgbolog=0;
|
|---|
| 177 | getData(0,sn,bolog,fgbolog);
|
|---|
| 178 | if(!(fgbolog & flgNotLookAt)) Slpb[1]->Push((r_8)sn,bolog);
|
|---|
| 179 | }
|
|---|
| 180 |
|
|---|
| 181 | for(int_4 sn=snb;sn<=sne;sn++) {
|
|---|
| 182 | int_4 lp = (mLPsnb<=sn && sn<=mLPsne) ? mLP%10: 0;
|
|---|
| 183 | r_8 bolog; uint_8 fgbolog=0;
|
|---|
| 184 |
|
|---|
| 185 | getData(0,sn,bolog,fgbolog);
|
|---|
| 186 | if(fgbolog & flgNotLookAt) {
|
|---|
| 187 | putData(0,sn,bolog,fgbolog);
|
|---|
| 188 | continue;
|
|---|
| 189 | }
|
|---|
| 190 |
|
|---|
| 191 | // i
|
|---|
| 192 | // <---lm---><--lg--><---lm--->
|
|---|
| 193 | // | c=0 | | c=1 |
|
|---|
| 194 | // ideb | | ifin
|
|---|
| 195 | // ifin ideb
|
|---|
| 196 | r_8 bolo; uint_8 fgbolo=0;
|
|---|
| 197 | int_4 ifill[2] = {sn-lgs2-1,sn+lgs2+lm};
|
|---|
| 198 | if(lp>1) printf(">>>> sn=%d v=(%f,%ld) ifill=[%d,%d]\n"
|
|---|
| 199 | ,sn,bolog,fgbolog,ifill[0],ifill[1]);
|
|---|
| 200 |
|
|---|
| 201 | if(ifill[0]<snb) {
|
|---|
| 202 | putData(0,sn,bolog,fgbolog);
|
|---|
| 203 | continue;
|
|---|
| 204 | }
|
|---|
| 205 | getData(0,ifill[0],bolo,fgbolo);
|
|---|
| 206 | if(fgbolo&flgNotLookAt) Slpb[0]->Pop();
|
|---|
| 207 | else Slpb[0]->Push((r_8)ifill[0],bolo);
|
|---|
| 208 | if((int_4)Slpb[0]->NPoints()<lmmin) {
|
|---|
| 209 | putData(0,sn,bolog,fgbolog);
|
|---|
| 210 | continue;
|
|---|
| 211 | }
|
|---|
| 212 |
|
|---|
| 213 | if(ifill[1]<=sne) {
|
|---|
| 214 | getData(0,ifill[1],bolo,fgbolo);
|
|---|
| 215 | if(fgbolo&flgNotLookAt) Slpb[1]->Pop();
|
|---|
| 216 | else Slpb[1]->Push((r_8)ifill[1],bolo);
|
|---|
| 217 | } else Slpb[1]->Pop();
|
|---|
| 218 | if((int_4)Slpb[1]->NPoints()<lmmin) {
|
|---|
| 219 | putData(0,sn,bolog,fgbolog);
|
|---|
| 220 | continue;
|
|---|
| 221 | }
|
|---|
| 222 |
|
|---|
| 223 | // Calcul des valeurs pour decider
|
|---|
| 224 | double s[2],slin[2],ypred[2];
|
|---|
| 225 | for(int c=0;c<2;c++) {
|
|---|
| 226 | r_8 mean,a0,a1;
|
|---|
| 227 | s[c] = Slpb[c]->Compute(mean);
|
|---|
| 228 | slin[c] = Slpb[c]->Compute(a0,a1);
|
|---|
| 229 | ypred[c] = a0+a1*sn;
|
|---|
| 230 | if(lp>1) printf("... c=%d m=%f s=%f ypred=%f slin=%f (%d)\n"
|
|---|
| 231 | ,c,mean,s[c],ypred[c],slin[c],Slpb[c]->NPoints());
|
|---|
| 232 | if(slin[c]<0. || s[c]<0.) {
|
|---|
| 233 | putData(0,sn,bolog,fgbolog);
|
|---|
| 234 | continue;
|
|---|
| 235 | }
|
|---|
| 236 | }
|
|---|
| 237 |
|
|---|
| 238 | if(intuple>0) {
|
|---|
| 239 | xnt_[0] = sn; xnt_[1] = bolog;
|
|---|
| 240 | xnt_[2] = s[0]; xnt_[3] = s[1];
|
|---|
| 241 | xnt_[4] = slin[0]; xnt_[5] = slin[1];
|
|---|
| 242 | xnt_[6] = ypred[0]; xnt_[7] = ypred[1];
|
|---|
| 243 | nt_.Fill(xnt_);
|
|---|
| 244 | }
|
|---|
| 245 |
|
|---|
| 246 | // Si les predictions ne sont pas compatibles (sig), pas de deglitch
|
|---|
| 247 | r_8 yy = ypred[1]-ypred[0];
|
|---|
| 248 | r_8 sy = nsg*sqrt(s[0]*s[0]+s[1]*s[1]);
|
|---|
| 249 | if(lp>1) printf("... yy=|%f| > sy=%f ???\n",yy,sy);
|
|---|
| 250 | if(fabs(yy) > sy) {
|
|---|
| 251 | putData(0,sn,bolog,fgbolog);
|
|---|
| 252 | continue;
|
|---|
| 253 | }
|
|---|
| 254 |
|
|---|
| 255 | // Si la valeur est entre les 2 predictions, pas de deglitch
|
|---|
| 256 | r_8 dp0 = ypred[0];
|
|---|
| 257 | r_8 dp1 = ypred[1];
|
|---|
| 258 | if(ypred[1]<ypred[0]) {dp0 = ypred[1]; dp1 = ypred[0];}
|
|---|
| 259 | if(lp>1) printf("... %f dans [%f,%f] ???\n",bolog,dp0,dp1);
|
|---|
| 260 | if(dp0<bolog && bolog<dp1) {
|
|---|
| 261 | putData(0,sn,bolog,fgbolog);
|
|---|
| 262 | continue;
|
|---|
| 263 | }
|
|---|
| 264 |
|
|---|
| 265 | // A + de N sigmas (lin) des 2 predictions
|
|---|
| 266 | dp0 = fabs(bolog-ypred[0]); slin[0] *= nsg;
|
|---|
| 267 | dp1 = fabs(bolog-ypred[1]); slin[1] *= nsg;
|
|---|
| 268 | if(lp>1) printf("... dp0=%f > %f et dp1=%f > %f ???\n"
|
|---|
| 269 | ,dp0,slin[0],dp1,slin[1]);
|
|---|
| 270 | if(dp0 < slin[0] || dp1 < slin[1]) {
|
|---|
| 271 | putData(0,sn,bolog,fgbolog);
|
|---|
| 272 | continue;
|
|---|
| 273 | }
|
|---|
| 274 |
|
|---|
| 275 | // C'est un glitch
|
|---|
| 276 | nflag++;
|
|---|
| 277 | fgbolog |= flgGlitch;
|
|---|
| 278 | if(lp>1)
|
|---|
| 279 | printf("......... Glitch en sn=%d ... %d glitches found\n",sn,nflag);
|
|---|
| 280 | putData(0,sn,bolog,fgbolog);
|
|---|
| 281 |
|
|---|
| 282 | }
|
|---|
| 283 | if(mLP) cout<<"FlagGlitch::run(): Fin de boucle: nflag="<<nflag
|
|---|
| 284 | <<" / "<<sne-snb+1<<endl;
|
|---|
| 285 |
|
|---|
| 286 | if(intuple>0) {
|
|---|
| 287 | char *fileout_ = "flagglitch.ppf";
|
|---|
| 288 | string tag_ = "nfg";
|
|---|
| 289 | POutPersist pos_(fileout_); pos_.PutObject(nt_,tag_);
|
|---|
| 290 | }
|
|---|
| 291 |
|
|---|
| 292 | //---------------------------------------------------------
|
|---|
| 293 | } catch (PException & exc) {
|
|---|
| 294 | cout<<"FlagGlitch: Catched Exception "<<(string)typeid(exc).name()
|
|---|
| 295 | <<"\n .... Msg= "<<exc.Msg()<<endl;
|
|---|
| 296 | }
|
|---|
| 297 |
|
|---|
| 298 | return;
|
|---|
| 299 | }
|
|---|
| 300 |
|
|---|
| 301 |
|
|---|
| 302 | ////////////////////////////////////////////////////////////////////////
|
|---|
| 303 | // Flag des sample autour d'un sample deja flaggue :
|
|---|
| 304 | // INPUT:
|
|---|
| 305 | // lm : longueur a flagguer de part et d'autre du mauvais sample
|
|---|
| 306 | // <---lm---><bad_sample><---lm--->
|
|---|
| 307 | ////////////////////////////////////////////////////////////////////////
|
|---|
| 308 |
|
|---|
| 309 | FlagAroundFlag::FlagAroundFlag(uint_4 lm,uint_8 flgs)
|
|---|
| 310 | {
|
|---|
| 311 | // Mise en forme des arguments
|
|---|
| 312 | LSide = lm;
|
|---|
| 313 |
|
|---|
| 314 | SetFlag();
|
|---|
| 315 | SetFlagAroundFlag(flgs);
|
|---|
| 316 | SetLimits();
|
|---|
| 317 | }
|
|---|
| 318 |
|
|---|
| 319 | FlagAroundFlag::~FlagAroundFlag()
|
|---|
| 320 | {
|
|---|
| 321 | }
|
|---|
| 322 |
|
|---|
| 323 | void FlagAroundFlag::Print(::ostream & os)
|
|---|
| 324 | {
|
|---|
| 325 | os<<"FlagAroundFlag::Print"
|
|---|
| 326 | <<" LSide="<<LSide
|
|---|
| 327 | <<" flgSample="<<flgSample<<" flgAround="<<flgAround<<endl
|
|---|
| 328 | <<" ... Range: ["<<VMin<<","<<VMax<<"]"<<" flgBad="<<flgBad<<endl;
|
|---|
| 329 | }
|
|---|
| 330 |
|
|---|
| 331 | void FlagAroundFlag::init()
|
|---|
| 332 | {
|
|---|
| 333 | cout << "FlagAroundFlag::init" << endl;
|
|---|
| 334 | declareInput("DataIn"); // input index 0
|
|---|
| 335 | declareOutput("DataOut"); // output index 0
|
|---|
| 336 |
|
|---|
| 337 | int_4 lbuff = 2*LSide +10;
|
|---|
| 338 | if(neededHistory<lbuff) setNeededHistory(lbuff);
|
|---|
| 339 | }
|
|---|
| 340 |
|
|---|
| 341 | void FlagAroundFlag::run()
|
|---|
| 342 | {
|
|---|
| 343 | long snb = getMinIn();
|
|---|
| 344 | long sne = getMaxIn();
|
|---|
| 345 |
|
|---|
| 346 | if(snb>sne) {
|
|---|
| 347 | cout<<"FlagAroundFlag::run() - Bad sample interval"<<snb<<" , "<<sne<<endl;
|
|---|
| 348 | throw ParmError("FlagAroundFlag::run() - Bad sample interval");
|
|---|
| 349 | }
|
|---|
| 350 | if(!checkInputTOIIndex(0)) {
|
|---|
| 351 | cout<<"FlagAroundFlag::run() - Input TOI (DataIn) not connected! "<<endl;
|
|---|
| 352 | throw ParmError("FlagAroundFlag::run() Input TOI (DataIn) not connected!");
|
|---|
| 353 | }
|
|---|
| 354 | if(!checkOutputTOIIndex(0)) {
|
|---|
| 355 | cout<<"FlagAroundFlag::run() - Output TOI (DataOut) not connected! "<<endl;
|
|---|
| 356 | throw ParmError("FlagAroundFlag::run() Output TOI (DataOut) not connected!");
|
|---|
| 357 | }
|
|---|
| 358 |
|
|---|
| 359 | //---------------------------------------------------------
|
|---|
| 360 | try {
|
|---|
| 361 |
|
|---|
| 362 | cout<<"FlagAroundFlag::run(): SampleNum de "<<snb<<" a "<<sne<<endl;
|
|---|
| 363 | uint_4 nflag=0, nflagged=0, nflagalready=0, nbad=0;
|
|---|
| 364 |
|
|---|
| 365 | // Un peu lourdingue! On pourrait optimiser!
|
|---|
| 366 | // Si on ne tue pas (LSide=0), alors juste un transfert de donnees
|
|---|
| 367 | for(int_4 sn=snb;sn<=sne;sn++) {
|
|---|
| 368 | r_8 bolo; uint_8 fgbolo=0;
|
|---|
| 369 | getData(0,sn,bolo,fgbolo);
|
|---|
| 370 | if(VMin<=VMax && (bolo<VMin || bolo>VMax)) { // Pixel out of range
|
|---|
| 371 | nbad++;
|
|---|
| 372 | fgbolo |= flgBad;
|
|---|
| 373 | } else if(fgbolo&flgSample) { // Pixel flaggue flgSample
|
|---|
| 374 | nflag++;
|
|---|
| 375 | } else if(fgbolo&flgAround) { // Pixel deja flagge around?
|
|---|
| 376 | nflagalready++;
|
|---|
| 377 | } else if(LSide>0) { // Pixel a flagguer?
|
|---|
| 378 | int_4 ideb = sn-LSide, ifin = sn+LSide;
|
|---|
| 379 | for(int_4 i=ideb;i<=ifin;i++) { // Y a t-il un "mauvais" pixel autour?
|
|---|
| 380 | if(i<snb || i>sne || i==sn) continue;
|
|---|
| 381 | r_8 b; uint_8 fg=0;
|
|---|
| 382 | getData(0,i,b,fg);
|
|---|
| 383 | if(fg&flgSample) {fgbolo |= flgAround; nflagged++; break;}
|
|---|
| 384 | }
|
|---|
| 385 | }
|
|---|
| 386 | putData(0,sn,bolo,fgbolo);
|
|---|
| 387 | }
|
|---|
| 388 | cout<<"FlagAroundFlag::run(): Fin de boucle: nflag="<<nflag
|
|---|
| 389 | <<" nflagged="<<nflagged<<" (already="<<nflagalready<<")"<<endl
|
|---|
| 390 | <<" ...... bad="<<nbad<<" / tot="<<sne-snb+1<<endl;
|
|---|
| 391 |
|
|---|
| 392 | //---------------------------------------------------------
|
|---|
| 393 | } catch (PException & exc) {
|
|---|
| 394 | cout<<"FlagAroundFlag: Catched Exception "<<(string)typeid(exc).name()
|
|---|
| 395 | <<"\n .... Msg= "<<exc.Msg()<<endl;
|
|---|
| 396 | }
|
|---|
| 397 |
|
|---|
| 398 | return;
|
|---|
| 399 | }
|
|---|