Read in data(2382 candidate tracks in 2 wedges)
dat=Table[Read["float_102831.txt",Real],{i,1,2382},{j,1,5}];Close["float_102831.txt"]
I am going to use wedge 3 only, because it has more data
As with the previous smaller sample, the smallest eigenvalue is too big to be just resolution
| 0.1877440906638008` | 0.31033913234003846` | 0.5019103800074148` | 0.7851928105886267` |
| 0.45820164928869167` | |||
| 0.5831392870921455` | 0.3918209464881586` | ||
| 0.7471787327501067` | 0.14142080710492277` |
I can, however, now plot the distribution of its projection since I have some more data
| 0.4754713729161473` | |||
| 0.6070203400012769` | 0.3766605208496238` | ||
| 0.7516288324672467` | 0.13457025381347892` |
This is (intuitively) the best candidate for the impact parameter wrt to the beam for this sample
Plot its distribution:
This is rather wide, but we used tracks of all Pts. Let's try to use high-Pt only
This is the best candidate for curvature:
Select tracks of high Pt (unclear how much, at this point). Find 123 tracks.
Now plot "impact parameter" for those:
Now plot the (single, because we lack XFT) constraint:
This looks vaguely similar to the one obtained "in the dark" , but much cleaner ! Note it comes out centered at zero without ad-hoc corrections.
now select the central peak:
| 0.47482059037024443` | |||
| 0.5804496434264329` | 0.3843728485385215` | ||
| 0.770647881938822` | 0.10569412662198875` |
Calculate curvature for this cleaned-up sample:
Select a sample of high-Pt tracks:
Plot their impact parameter:
This is kind of wide, but we need to correct for beam position ! Let's plot d vs. phi:
THIS LINEAR RELATIONHIP IS FURTHER CLEAR EVIDENCE WE ARE OBSERVING REAL TRACKS !!!
now I want to fit the beam "by eye", and evaluate d wrt the beam :
OK, so this is the expression of d wrt the beam:
And this is its distribution:
This is rather clean now.
This is VERY CLEAN now !!! The beam constraint cleans up the sample a lot !!!
...so there is no trick,the correlation between the two variables is given by physics!!!
Finally, select a sample of good tracks, by cutting both on beam and K:
From this, I get the following estimates: