Changeset 2031 in Sophya for trunk/ArchTOIPipe
- Timestamp:
- May 30, 2002, 2:53:22 PM (23 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/ArchTOIPipe/TestPipes/aksj02.cc
r2028 r2031 289 289 offes.SavePolyNTuple(soff_nt, soff_nt_ppfname); 290 290 291 SimpleOffsetEstimator offes2;291 SimpleOffsetEstimator* offes2; 292 292 if (soff2_dofg) { 293 293 cout << "> Creating 2nd SimpleOffsetEstimator(" << soff2_mwsz << "," << soff2_nptfit 294 294 << "," << soff2_degpol << ")" << endl; 295 offes2.SetParams(soff2_mwsz, soff2_nptfit, soff2_degpol); 296 offes2.SavePolyNTuple(soff2_nt, soff2_nt_ppfname); 297 } 298 299 if (gfilt_fg) 295 offes2 = new SimpleOffsetEstimator(soff2_mwsz, soff2_nptfit, soff2_degpol); 296 offes2->SavePolyNTuple(soff2_nt, soff2_nt_ppfname); 297 } 298 299 SimpleFilter *filt = NULL; 300 if (gfilt_fg) { 300 301 cout << "> Creating a GaussianFilter SimpleFilter Object " << endl; 301 double G_sigma = gfilt_sigma;302 double G_a = 1./(G_sigma*sqrt(M_PI*2.));303 SimpleFilter filt(gfilt_wsz, SimpleFilter::GaussFilter, G_a, G_sigma);304 302 double G_sigma = gfilt_sigma; 303 double G_a = 1./(G_sigma*sqrt(M_PI*2.)); 304 filt = new SimpleFilter(gfilt_wsz, SimpleFilter::GaussFilter, G_a, G_sigma); 305 } 305 306 306 307 FITSTOIReader* rbgal = NULL; … … 317 318 if (soff2_dofg && bgalfile && bgal2cut) { 318 319 cout << " offes2.SetBGalCut( " << bmin2 << "," << bmax2 << ")" << endl; 319 offes2 .SetBGalCut(bmin, bmax);320 offes2->SetBGalCut(bmin, bmax); 320 321 } 321 322 … … 336 337 outname = "bgalcopie"; 337 338 plombier.Connect(*rbgal, bgaltoi, offes, inname); 338 plombier.Connect(offes, outname, offes2, inname);339 plombier.Connect(offes, outname, *offes2, inname); 339 340 } 340 341 else if (soff2_dofg && bgal2cut) { 341 342 inname = "bgal"; 342 plombier.Connect(*rbgal, bgaltoi, offes2, inname);343 plombier.Connect(*rbgal, bgaltoi, *offes2, inname); 343 344 } 344 345 else if (bgalcut) { … … 350 351 inname = "in"; 351 352 outname = "out"; 352 plombier.Connect(offes, outname, offes2, inname);353 plombier.Connect(offes, outname, *offes2, inname); 353 354 if (bgalcut) { 354 355 inname = "bgal"; 355 356 outname = "bgalcopie"; 356 plombier.Connect(offes, outname, offes2, inname);357 plombier.Connect(offes, outname, *offes2, inname); 357 358 } 358 359 } … … 360 361 if (gfilt_fg) { 361 362 if (soff2_dofg) 362 plombier.Connect( offes2, "out",filt, "in");363 plombier.Connect(*offes2, "out", *filt, "in"); 363 364 else 364 plombier.Connect(offes, "out", filt, "in");365 plombier.Connect( filt, "out", w, "outfcln", "", 0, fg_wrtflag);365 plombier.Connect(offes, "out", *filt, "in"); 366 plombier.Connect(*filt, "out", w, "outfcln", "", 0, fg_wrtflag); 366 367 } 367 368 else { 368 369 if (soff2_dofg) 369 plombier.Connect( offes2, "out", w, "outcln", "", 0, fg_wrtflag);370 plombier.Connect(*offes2, "out", w, "outcln", "", 0, fg_wrtflag); 370 371 else 371 372 plombier.Connect(offes, "out", w, "outcln", "", 0, fg_wrtflag); … … 376 377 plombier.Connect(offes, "offset", w, "offset"); 377 378 if (soff2_dofg) 378 plombier.Connect( offes2, "offset", w, "offset2");379 plombier.Connect(*offes2, "offset", w, "offset2"); 379 380 } 380 381 if (fg_wrtms) { … … 388 389 cout << cleaner; 389 390 cout << offes; 390 if (soff2_dofg) cout << offes2;391 if (soff2_dofg) cout << *offes2; 391 392 cout << w; 392 393 … … 405 406 cout << cleaner; 406 407 cout << offes; 407 if (soff2_dofg) cout << offes2;408 if (gfilt_fg) cout << filt;408 if (soff2_dofg) cout << *offes2; 409 if (gfilt_fg) cout << *filt; 409 410 cout << w; 410 411 … … 412 413 PrtTim("End threads"); 413 414 if (bgalcut) delete rbgal; 414 415 if (gfilt_fg) delete filt; 416 if (soff2_dofg) delete offes2; 415 417 } 416 418 catch (PThrowable & exc) {
Note:
See TracChangeset
for help on using the changeset viewer.