From ebc0f100371942e99c433cf967b55524a9765c87 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Mon, 29 Jul 2019 12:18:40 -0500 Subject: switch to using the multiphoton PCA time This commit updates the likelihood to use the multiphoton PCA time instead of the usual pt time. When looking at the reconstruction of muons in run 10,000 I noticed that the PMT hit times for the PMTs which had really high charge was all over the place. There were PMTs that were very close to each other with hit times differing by ~ 20 ns. I'm still not entirely sure what causes this (is it some hardware issue with the discriminator or is it a problem with the charge walk correction which always assumes a single PE?), but the multiphoton PCA times looked a lot more reasonable. Eventually I need to look into the ptms variable which is the multiphoton PCA transit time RMS. --- src/zdab_utils.c | 7 ++++++- 1 file changed, 6 insertions(+), 1 deletion(-) (limited to 'src/zdab_utils.c') diff --git a/src/zdab_utils.c b/src/zdab_utils.c index 48ad9a1..a7669fa 100644 --- a/src/zdab_utils.c +++ b/src/zdab_utils.c @@ -112,7 +112,6 @@ int get_event(zebraFile *f, event *ev, zebraBank *bev) ev->pmt_hits[id].hit = 1; ev->pmt_hits[id].ept = bpmt.ept; - ev->pmt_hits[id].t = bpmt.pt; ev->pmt_hits[id].qihl = bpmt.pihl; ev->pmt_hits[id].qihs = bpmt.pihs; ev->pmt_hits[id].qilx = bpmt.pilx; @@ -123,7 +122,13 @@ int get_event(zebraFile *f, event *ev, zebraBank *bev) ev->pmt_hits[id].qhs = bpmt.phs; ev->pmt_hits[id].qlx = bpmt.plx; ev->pmt_hits[id].pf = bpmt.pf; + ev->pmt_hits[id].pt = bpmt.pt; ev->pmt_hits[id].pt1 = bpmt.pt1; + ev->pmt_hits[id].ptm = bpmt.ptm; + + /* Use the multiphoton PCA time since we are looking at high energy + * events. */ + ev->pmt_hits[id].t = bpmt.ptm; /* Clear the PMT_FLAG_DIS bit. */ ev->pmt_hits[id].flags &= ~PMT_FLAG_DIS; -- cgit