Changeset 3920 in Sophya for trunk/AddOn/TAcq/racquproc.cc
- Timestamp:
- Dec 3, 2010, 4:51:54 AM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/AddOn/TAcq/racquproc.cc
r3918 r3920 181 181 Histo* phexp[2*MAXNBFIB]; 182 182 TMatrix<r_4>* psatur[2*MAXNBFIB]; 183 TVector<uint_4>* pvsatfreq[2*MAXNBFIB]; 183 184 for(int ifib=0; ifib<2*MAXNBFIB; ifib++) { 184 185 phist[ifib]=NULL; ph2[ifib]=NULL; phexp[ifib]=NULL; psatur[ifib]=NULL; … … 190 191 phexp[lc] = new Histo(-128.5,+128.5,257); 191 192 psatur[lc] = new TMatrix<r_4>(3,3); 193 pvsatfreq[lc] = new TVector<uint_4>(spectre_.NCols()); 192 194 *psatur[lc] = 0.; 193 195 } … … 254 256 else if (tbcp[j].imagI()==127) sii=1; 255 257 (*psatur[lc])(sir+1,sii+1)+=1.; 256 258 if ((sir!=0)||(sii!=0)) (*pvsatfreq[lc])(j)+=1; 257 259 spectre_(lc,j) += Zmod2(tbcp[j]); 258 260 } … … 271 273 else if (tbcp[j].imagI()==127) sii=1; 272 274 (*psatur[lc])(sir+1,sii+1)+=1.; 275 if ((sir!=0)||(sii!=0)) (*pvsatfreq[lc])(j)+=1; 273 276 274 277 spectre_(lc,j) += Zmod2(tbcp[j]); … … 333 336 sprintf(buftag,"saturV%d",(int)lc); 334 337 po << PPFNameTag(buftag) << (*psatur[lc]); 335 cout << " Matrice Saturation Voie[" << lc << "] : " << (*psatur[lc]);336 if (nzm_[lc]>0) {337 (*psatur[lc]) /= ((r_4)nzm_[lc]*0.01);338 cout << " Matrice Saturation Pourcentage Voie[" << lc << "] : " << (*psatur[lc]);339 }340 338 (*psatur[lc])=0.; 339 sprintf(buftag,"satfreqV%d",(int)lc); 340 po << PPFNameTag(buftag) << (*pvsatfreq[lc]); 341 (*pvsatfreq[lc])=0; 341 342 } 342 343 nzm_[lc]=0; … … 389 390 sprintf(buftag,"saturV%d",(int)lc); 390 391 po << PPFNameTag(buftag) << (*psatur[lc]); 391 cout << " Matrice Saturation Voie[" << lc << "] : " << (*psatur[lc]); 392 if (nzm_[lc]>0) { 393 (*psatur[lc]) /= ((r_4)nzm_[lc]*0.01); 394 cout << " Matrice Saturation Pourcentage Voie[" << lc << "] : " << (*psatur[lc]); 395 } 396 (*psatur[lc])=0.; 392 sprintf(buftag,"satfreqV%d",(int)lc); 393 po << PPFNameTag(buftag) << (*pvsatfreq[lc]); 394 delete psatur[lc]; 395 delete pvsatfreq[lc]; 397 396 } 398 397 nzm_[lc]=0; … … 413 412 << " -> " << fracsame << " % )" << endl; 414 413 } 415 416 for(sa_size_t lc=0; lc<2*memgr.NbFibres(); lc++) {417 if (phist[lc]) delete phist[lc];418 if (ph2[lc]) delete ph2[lc];419 if (phexp[lc]) delete phexp[lc];420 if (psatur[lc]) delete psatur[lc];421 }422 423 414 return 0; 424 415 }
Note:
See TracChangeset
for help on using the changeset viewer.