Total scattering was calculated using the formula:
Theta_0 = 13.6 MeV/(beta*c*p) * Z * sqrt(SUM) * {1 + 0.038 * ln(SUM)}
where SUM = sum(x_i / X0_i)
Taking into account electrons, than beta = Z = 1
Total scattering angle as a function of polar angle and momenta. |
---|
![]() |
Total scattering angle as a function of momenta for different polar angles. |
---|
![]() |
In the calculation, the polar angle dependant deflection between the perpendicular direction to the surface of the RICH components (shortest distance through the materials) and particle incident angle (real particle path lengt length through the material) was taken into account.
Dependance of the deflection angle on the polar angle. |
---|
![]() |
Small angle scattering in RICH radiator as a function of polar angle and momenta. |
---|
![]() |
Dependance of the particle trajectory length in the RICH radiator on the polar angle. |
---|
![]() |
Scattering angle in RICH radiator as a function of momenta for different polar angles. | Scattering angle in RICH radiator as a function of momenta for polar angle Theta = 50deg. |
---|---|
![]() | ![]() |
In the following, only the mean deflection to the perpendicular direction of 22.25deg was taken into account to make the pictures simple
Scattering angle in RICH mirror substrate material SIGRADUR as a function of momenta for mean deflection angle 22.25deg |
---|
![]() |
Scattering in Aluminium. | Scattering in MgF2. |
---|---|
![]() | ![]() |
Scattering in CFK |
---|
![]() |
Area difference between these two cuts is:
I have compared the 1D distribution for (theta_rich - theta_RK) (in phi it looks similar) for electrons comming coming from conversion and the electrons comming coming from all other sources. the conversion signal I have chosen according to the following condition: "pid.genInfo1==7001 || pid.genInfo1==17001 || pid.genInfo1==52001", the "real" signal is negation of the previous condition. the determination of the mean and sigma of the distribution was done with the "standard" 2-gaussian fit for momentum bins (0,150>; (150,350> and (350,1000>
in the picture "compareSim.png" you can see mean and sigma of the delta_theta distributions for the 3 momenta bins. the upper two pictures are done for electrons and positrons separately (red color "real" signal, blue conversion signal). in the lowest picture the comparision between electrons and positrons for the "real" signal is done.
as one can see, the differences in the simulation are only minor. the yields for signal and conversion are comparable (almost 1:1) for momenta (0,150>, it drops to 2:1 for momenta (350,1000> (see thetaNeg.jpg)
The sigma of the delta_theta distributions for negative tracks is:
The sigma of the delta_theta distributions for positive tracks is:
==================================================
in phi are the distributions much more narrower:
The sigma of the delta_phi distributions for negative tracks is:
The sigma of the delta_phi distributions for positive tracks is:
At this point i have to add, that i set following limits for the fit parameters: for signal is sigma in the range (0.4; 2.0) and for the background (2.0; infinity), which you can see in the fits not very well describing the distributions. in the most cases, the background had sigma of 2.0 in few cases just a touch higher than 2.0! (not a very good model :-(...), since for experimental data this limits worked very well, i decided not to introduce new limits for the fits.
I hope, the selection of the conversion was done properly (according to the documentation i hope i have not overseen something...)
-- MartinJurkovic - 29 May 2007
I | Attachment | Action | Size | Date | Who | Comment |
---|---|---|---|---|---|---|
![]() | scattering50.txt | manage | 3.3 K | 12 Jun 2007 - 12:32 | MartinJurkovic | Total scattering and scattering in radiator as a function of momenta at polar angle 50deg |
![]() | theta2length.txt | manage | 1.3 K | 12 Jun 2007 - 12:31 | MartinJurkovic | Path length and deflection angle as a function of polar angle |