| 20 |
c common noise computation |
c common noise computation |
| 21 |
c |
c |
| 22 |
c------------------------------------------------------------------------ |
c------------------------------------------------------------------------ |
| 23 |
parameter (scut=3.) |
parameter (scut=4.)!(scut=3.) |
| 24 |
|
parameter (NSTRIPMIN=64)!(NSTRIPMIN=10)!TEMPORANEO |
| 25 |
|
parameter (NSTRIPWARNING=64) |
| 26 |
|
|
| 27 |
|
real cn(nviews,nva1_view) ! common noise values |
| 28 |
|
integer cnflag(nviews,nva1_view) ! common noise flags: |
| 29 |
|
c 0 = ok |
| 30 |
|
c -1 = failed (less than NSTRIPMIN strips used) |
| 31 |
|
c 1 = warning (less than NSTRIPWARNING strips used) |
| 32 |
|
integer mask_vk_ev(nviews,nva1_view) ! mask of VA1s, local to the event (set to 0 if CN computation fails) |
| 33 |
|
|
|
real cn(nviews,nva1_view) !common noise value |
|
| 34 |
integer clstr(nviews,nva1_view,nstrips_va1) !flag matrix to mark |
integer clstr(nviews,nva1_view,nstrips_va1) !flag matrix to mark |
| 35 |
! strips with signal and exclude them from |
! strips with signal and exclude them from |
| 36 |
! common noise computation |
! common noise computation |
| 38 |
! strips with unusually high or low signal and |
! strips with unusually high or low signal and |
| 39 |
! exclude them from common noise computation |
! exclude them from common noise computation |
| 40 |
|
|
| 41 |
common/cn/cn,clstr,strange |
common/cn/cn,cnflag,mask_vk_ev,clstr,strange |
| 42 |
|
|
| 43 |
|
|
| 44 |
c$$$c------------------------------------------------------------------------ |
c$$$c------------------------------------------------------------------------ |