Action #3549
Implement TOFCOR when reading events from "old" EVP datasets
| Status: | Closed | Start date: | 02/19/2021 | |
|---|---|---|---|---|
| Priority: | Urgent | Due date: | ||
| Assigned To: | % Done: | 100% | ||
| Category: | - | |||
| Target version: | 2.0.0 | |||
| Duration: |
Description
According to com-rp-rol-drg-057.pdf a TOF corrections needs to be applied to EVP datasets of version 2 and older (see DSD_REP keyword in the EVP file). An example of an old EVP dataset is the one in viewing period 1 of the Crab.
The TOF correction is implemented in TOFCOR in the DAL library DALAAA18 (see dalaaa19.tofcor.f). The correction is automatically implemented when reading an EVP file (see dalaaa19.pevpsr.f):
C VERNO integer data set representation version PEVPRI *
C-- De-scale parameter and transform from INTEGER*2 to
C argument spec
C Copy data from COMMON PEVPD1 to return variables
TIMTAG(1) = TIME(1)
TIMTAG(2) = TIME(2)
MPAR( 1) = EN(1)
MPAR( 2) = EN(2)
MPAR( 3) = EN(1) + EN(2)
MPAR( 4) = 0
MPAR( 5) = PHI
MPAR( 6) = 0
MPAR( 7) = LOC(1)/32.
MPAR( 8) = LOC(2)/32.
MPAR( 9) = LOC(3)/256.
MPAR(10) = LOC(4)/32.
MPAR(11) = LOC(5)/32.
MPAR(12) = LOC(6)/256.
MPAR(13) = LOC(7)/128.
MPAR(14) = LOC(8)/128.
CLASS = CL
MODCOM = MC
REFLAG = RF
VETO = VF
EVSCT(1) = CO(1)
EVSCT(2) = CO(2)
EVSCT(3) = CO(3)
EVSCT(4) = CO(4)
EVSCT(5) = EH
CTYP = CLASS(1:3)
C
C Now do the TOF correction for dataset versions < 3
C
IF(VERNO.LT.3)THEN
C
C Do not correct for SIM events
C
IF (CTYP .NE. 'SIM') THEN
IF (REFLAG .GE. 4) THEN
CALL TOFCOR(MPAR(1),MPAR(2),MPAR(14),IDUM)
ENDIF
ENDIF
ENDIF
ELSE
An equivalent TOF correction should be implemented in
GCOMEventList::read_events.
Recurrence
No recurrence.
History
#1
Updated by Knödlseder Jürgen over 5 years ago
- Status changed from New to In Progress
- % Done changed from 0 to 10
I recognised that VP 0001.0 includes a newer EVP with TOFCOR applied. However a valid dataset for testing is VP 0002.5 which has an old data representation.
#2
Updated by Knödlseder Jürgen over 5 years ago
- File vp0002_5.png added
- File vp0002_5-tofcor.png added
- % Done changed from 10 to 50
I implemented the TOFCOR correction, the results are displayed in the panels below. The top panel is the original panel without ToF correction which shows some structure in the ToF distribution due to oversampling. The bottom panel is after application of the TOFCOR correction. Note that the correction considerably modifies the ToF distribution, lacking events with small ToF. Note that the nominal ToF window is 115-130.
Before TOFCOR correction
After TOFCOR correction
#3
Updated by Knödlseder Jürgen over 5 years ago
- % Done changed from 50 to 60
vp0002_0vp0002_5vp0003_0vp0004_0vp0006_0vp0007_0vp0008_0vp0009_5vp0011_0vp0013_5vp0015_0vp0020_0vp0021_0vp0023_0vp0207_0vp0208_0vp0209_0vp0211_0vp0212_0vp0213_0vp0216_0vp0218_0vp0219_4vp0220_0vp0222_0vp0223_0vp0224_0vp0226_0vp0229_0
#4
Updated by Knödlseder Jürgen over 5 years ago
- Status changed from In Progress to Pull request
- % Done changed from 60 to 90
#5
Updated by Knödlseder Jürgen over 5 years ago
- Status changed from Pull request to Closed
- % Done changed from 90 to 100


