R3BRootGroup / R3BRoot

Framework for Simulations and Data Analysis of R3B Experiment
https://github.com/R3BRootGroup/R3BRoot/wiki
GNU General Public License v3.0
18 stars 102 forks source link

Updated IncomingBeta reconstruction for s091 experiment. #1063

Closed ryotani closed 5 months ago

ryotani commented 6 months ago

Clean-up and minor improvements in los classes.

Include the timing of the S2 hit after subtraction of the reference time in FrsSciPosCalData.

Bug fix in R3BFrsSciTcal2Cal when obtaining good pairs of hits of left and right PMTs.

Update IncomingData to reconstruct S2(FrsSci) to LosCal for the s091 experiment. Los CFD out was used for the Tref signal for S2, so LosCal is used for reconstructing ToF specific to the experiment. When multiple hits in Los are seen, only good hits wrt the trigger timing are chosen.

Describe your proposal.

Mention any issue this PR is resolved or is related to.


Checklist:

nmozumdar commented 6 months ago

Hi Ryo, just looking through it fast, your updates look ok from my side. However I should give you a heads up that in the next step you will have to match the Musli hit for the charge to the corresponding los hit. Which means if you take your current approach to selecting the hit within the task you will have to copy the same stuff multiple times. In fact it will be even more complicated within the Musli Task. This is why I took the approach of selecting the hit within the TProvideTStart and then matching it to the musli hit hence selecting for the correct Frs Z. Please have a look at he MusliMapped2Cal task.

ryotani commented 6 months ago

Thanks @nmozumdar for the comment. Actually, your point is one of the difficult parts for me to write down in the code.

As you may remember, the trigger timing distributed to the DAQ systems was one of the big confusions during the experiment last month. The tref signal sent to S2 was for some reason having a jittering. Thus we should live with another trigger reference signal; that is a CFD out from a channel of the LOS PMT. While TStart is deduced by the difference from the average timing of 8 PMTs in LOS to the trigger timing, the definition of reference timing for S2 DAQ is slightly different. As I didn't want to change TStartProvider too much, I decided to code in the R3BIncomingBeta instead. In the actual analysis for this experiment, I think one can set the same condition in TStartProvider with the SetWindow function to select the same LOS hit for all the analysis: https://github.com/R3BRootGroup/R3BRoot/blob/4f4b89370d2bbf0c0deb25aec87f6394542381be/los/R3BLosProvideTStart.h#L34-L39

ryotani commented 5 months ago

@jose-luis-rs Hi, this is also fine to be merged.

jose-luis-rs commented 5 months ago

Hello @ryotani , did you check that the code still works fine for the S509 PID?

jose-luis-rs commented 5 months ago

for example for run0190

ryotani commented 5 months ago

@jose-luis-rs I fixed a bug. Now it should be alight with s509.

jose-luis-rs commented 5 months ago

Thanks @ryotani Please, could you squash your commits?

ryotani commented 5 months ago

@jose-luis-rs Is this fine?

jose-luis-rs commented 5 months ago

@jose-luis-rs Is this fine?

Yes, thanks