[Halld-cpp] FMWPCHits in REST

Alexander Austregesilo aaustreg at jlab.org
Tue Jan 9 10:31:03 EST 2024


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



More information about the Halld-cpp mailing list