Choose timezone
Your profile timezone:
Powered by Indico
v3.3.3
Decription of the first part of the pulse shape code which Tapas and I have been looking at :
#preprocessing/getStepTimes.cc
Loop over each PMT (482 !TBC)
For each PMT, check all the entries (all the photons)
A function calculates one Volume ID for one PMT number
This calculated VolumeID is matched to the one given in the ROOT file (A mapping of PMT numbers to Volume-ID)
As it seems from the code, the VolumeID in the ROOT file should be the ID of the first volume of a PMT
e.g.
PMT # 1 (Top array)
Maps to 2250 (VolumeID)
All other volumes under PMT#1 have ids (2250 + i), i=1..8
TOP PMT (247 TBC) corresponds to Volume-ID (2250 TBC)
BOTTOM PMT (241 TBC) corresponds to Volume-ID (74 - 2234 TBC)
The above Volume-ID numbers are based on 9 volumes per PMT
In current Simulation, it seems to be 6 (TBC)
For each PMT,
Loop over all the photons in the event
once there is a match between PMT# and Volume-ID
Get the StepTime corresponding to that PMT for that photon
e.g. StepTime[0][i] = <val>
- 0 corresponds to the first PMT in the array
- i corresponds to the ith photon hitting PMT-0
## preprocessing/light_distribution_s2.cc
# Calculation of light collection efficiency
Loop over each photon in the event
Look for the Top/Bottom PMT
Match the VolumeID (Cathode part of the PMT corresponds to the first volume of PMT) with the PMT number
Store an incremetal histogram for each PMT that reflects the number of photons hitting the PMT
Total bin content of the histogram reflects the total number of photons collected
Each bin number is equal to PMT#
Each bin content is the number of photons collected at that PMT#
The number of photons are divided by the number of initial photons \
in order to calculate light collection efficiency scaled by 0.3 (? TBC)
Also, the number of initial photons is 1e6 ( ? TBC)