[Halld-cpp] FMWPCHits in REST

[Alexander Austregesilo]

I just updated the branch in order to save the raw integral as "charge" 
q. In addition, I set the threshold to 0 and removed the requirement on 
the quality factor so that really all hits are saved to the output.

I noticed that you create a DFMWPCMatchedTrack In the DCPPEpEm_factory 
to determine the closest wire and the number of FMWPC hits for each 
track. At the moment, it is not clear t me how this would work with a 
REST file. Simon, do you have a recommendation how we can handle this? I 
assume, we have to save more than just the raw hits to the REST file.


On 1/9/2024 10:32 AM, Alexander Austregesilo wrote:
> The branch can be found here: 
> https://github.com/JeffersonLab/halld_recon/tree/aaust_fmwpc_reco2
>
> On 1/9/24 10:31, Alexander Austregesilo wrote:
>> Hi Ilya at el,
>>
>> I have a halld_recon branch where the uncalibrated FMWPCDigiHits are 
>> basically copied to FMWPCHit objects and written to the REST output 
>> by default. There is currently no threshold on the amplitude and the 
>> default timing cut is in the range [-60,900] where everything should 
>> be accepted. The only difference between digihits and hits is the 
>> removal of objects with quality factor != 0 and the scaling of the 
>> time from samples to nanoseconds:
>>
>>>   t_scale = 8.0/10.0;    // 8 ns/count and integer time is in 1/10th 
>>> of sample
>>
>> The member "q" currently has the same value as the raw amplitude 
>> "amp", as I was not sure if the calculation of the integral and the 
>> pedestal subtraction is implemented correctly.
>>
>> Reading in a REST files, the FMWPCHits can be accessed by this piece 
>> of code:
>>
>>> vector<const DFMWPCHit*> myfmwpchits;
>>> loop->Get(myfmwpchits);
>>
>> Please have a look and let me know if anything else is needed in the 
>> REST files.
>>
>> Cheers,
>>
>> Alex
>>
-- 
Alexander Austregesilo

Staff Scientist - Experimental Nuclear Physics
Thomas Jefferson National Accelerator Facility
Newport News, VA
aaustreg at jlab.org
(757) 269-6982