energy reconstruction in esaf by firenze group
DESCRIPTION
Data loaded from: Patern recognition (cluster of RecoPixels) TrackDirectionModules(1,2): Theta and Phi (and hopefully the track line on FS – not yet ). Energy Reconstruction in ESAF by Firenze group. The basic idea is to relate number of p.e. at a given time to Energy of EAS. - PowerPoint PPT PresentationTRANSCRIPT
![Page 1: Energy Reconstruction in ESAF by Firenze group](https://reader036.vdocument.in/reader036/viewer/2022070414/56814e44550346895dbbbc40/html5/thumbnails/1.jpg)
Energy Reconstruction in ESAFby Firenze group
Data loaded from: Patern recognition (cluster of RecoPixels) TrackDirectionModules(1,2): Theta and Phi (and
hopefully the track line on FS – not yet)
6-7 October 2005 Grenoble
The basic idea is to relate number of p.e. at a given time to Energy of EAS
Np.e.(T) = E/E1 x Q x Efficacy(x,y)IsInsidePixel() x Flux(T),
where Flux(T) = Ne(T) x Fluo Yield x Track Length x Transmission x Solid Angle
![Page 2: Energy Reconstruction in ESAF by Firenze group](https://reader036.vdocument.in/reader036/viewer/2022070414/56814e44550346895dbbbc40/html5/thumbnails/2.jpg)
● An example of event to deal with
● … from EnergyViewer
Energy Reconstruction in ESAF
![Page 3: Energy Reconstruction in ESAF by Firenze group](https://reader036.vdocument.in/reader036/viewer/2022070414/56814e44550346895dbbbc40/html5/thumbnails/3.jpg)
Energy Reconstruction in ESAF Reconstruction algorithm:
Sort pixels in GTU and prepare maps <GTU, vector<RecoPixelData*>>
Fit the track on the FS as a straight line and get ThetaFoV and PhiFoV as a function of the track line position
Find ThetaFoV at the shower maximum Given Theta, Phi of the ShowerTrack and ThetaFoV
and estimate for the shower maximum altitude (from HmaxXXXModules) estimate fluorescence yield at the shower maximum and mean attenuation. Make also a spectrum of fluorescent photons at the entrance pupil taking into account the spectrum at the production point and attenuation.
Make the background estimation: various methods (examing pixels around the track, looking track pixels but at different moments of time) – gives agreement results with simulated background
6-7 October 2005 Grenoble
OK
OK but should be taken from TDM
OK, need to iteratedOK but I want to improve it
OK
![Page 4: Energy Reconstruction in ESAF by Firenze group](https://reader036.vdocument.in/reader036/viewer/2022070414/56814e44550346895dbbbc40/html5/thumbnails/4.jpg)
Energy Reconstruction in ESAF Reconstruction algorithm:
Compute efficacy for the track (take into account dead spaces of PMTs and Pixels, holes in FS)
Simulate ShowerTrack object with E=10^21 eV and reconstructed Theta and Phi and with the maximum of around the found one.
Search for the Cherenkov peak Fit the time distrubution (taking into account photon to
p.e. conversion, optics and FS, atmosphere attenuation, shower production in space.
All steps are extensively and optionaly debugged (created the gzipped PS file with all relevant information).
6-7 October 2005 Grenoble
OK
OK in general but to be extensively testedOK but not usedOK, should be improved to take all steps of production OK
![Page 5: Energy Reconstruction in ESAF by Firenze group](https://reader036.vdocument.in/reader036/viewer/2022070414/56814e44550346895dbbbc40/html5/thumbnails/5.jpg)
● Shower Track on FS
Energy Reconstruction in ESAF
![Page 6: Energy Reconstruction in ESAF by Firenze group](https://reader036.vdocument.in/reader036/viewer/2022070414/56814e44550346895dbbbc40/html5/thumbnails/6.jpg)
● ThetaFoV and PhiFoV vs GTU
Energy Reconstruction in ESAF
![Page 7: Energy Reconstruction in ESAF by Firenze group](https://reader036.vdocument.in/reader036/viewer/2022070414/56814e44550346895dbbbc40/html5/thumbnails/7.jpg)
● X,Y vs GTU
Energy Reconstruction in ESAF
![Page 8: Energy Reconstruction in ESAF by Firenze group](https://reader036.vdocument.in/reader036/viewer/2022070414/56814e44550346895dbbbc40/html5/thumbnails/8.jpg)
● Recovering efficacy
Energy Reconstruction in ESAF
![Page 9: Energy Reconstruction in ESAF by Firenze group](https://reader036.vdocument.in/reader036/viewer/2022070414/56814e44550346895dbbbc40/html5/thumbnails/9.jpg)
● Recovering efficacy
Energy Reconstruction in ESAF
![Page 10: Energy Reconstruction in ESAF by Firenze group](https://reader036.vdocument.in/reader036/viewer/2022070414/56814e44550346895dbbbc40/html5/thumbnails/10.jpg)
● Recovering efficacy
Energy Reconstruction in ESAF
![Page 11: Energy Reconstruction in ESAF by Firenze group](https://reader036.vdocument.in/reader036/viewer/2022070414/56814e44550346895dbbbc40/html5/thumbnails/11.jpg)
● Recovering efficacy
Energy Reconstruction in ESAF
![Page 12: Energy Reconstruction in ESAF by Firenze group](https://reader036.vdocument.in/reader036/viewer/2022070414/56814e44550346895dbbbc40/html5/thumbnails/12.jpg)
● Another example
Energy Reconstruction in ESAF
![Page 13: Energy Reconstruction in ESAF by Firenze group](https://reader036.vdocument.in/reader036/viewer/2022070414/56814e44550346895dbbbc40/html5/thumbnails/13.jpg)
Energy Reconstruction in ESAF
![Page 14: Energy Reconstruction in ESAF by Firenze group](https://reader036.vdocument.in/reader036/viewer/2022070414/56814e44550346895dbbbc40/html5/thumbnails/14.jpg)
Curent and Future work
• Debuging all steps• Make better fluorescence shower (for all GTU own, and not the mean as it is now• Make backscattered Cherenkov light• Improve Fitting procedure (make it more complex).
Energy Module is now splitted:1. EnergyModule class2. EnergyViewer class3. EnergyFitUtils – collection of C routines