Uploading a new reprocessing from MAUS 2.1.

Also adding a script and pdf plot of the emittance across the tracker in 10MeV slices -- particles are randomly assigned to 'station 1' or 'station 5', then emittance is calculated and plotted as a function of pz. Particles in this sample are equivalent to the "cut14" passed samples.

(No corrections or mean-position subtraction is done, the emittance here is just the 4th root of the determinant of the covariance matrix. Vertical error bars are statistical (emittance/sqrt(N)) and horizontal error bars are the bin width)

Attaching the extracted reconstructed MC, going through the same program (linked from GitHub), and the pass/failed cuts plots

There's a new variable in this root file (cut_tof1_tku_momentum) that attempts to define a region on a plot of momentum at TKU vs momentum at TOF1. However, this cut does not make its way into the overall 'cut_allPassed' variable. The variable is there to look at, but it's not yet used (nor does it appear in the many pass/fail plots).

The distributions in the MC still look odd in the tracker (not as odd as before, but still odd)

Sure. (There's no 'truth' MC here, just the reconstructed MC, for the record)

The spill number comes from the MAUS spill objects GetSpillNumber(), but the reconstructed event number comes from counting the loops over the reconstructed events object (Adam told me a while ago not to trust the event number function, and instead to just use the loop counter -- I don't know if that's still true, but it's what I'm still doing)

The spill number is in the branch 'SpillNumber', and the reconstructed event number is in the branch 'ReconstructedEventNumber'

(The git repository should be up to date -- so if you see something silly you can track it down)

Hum, I tried loading event and spill number as in attached file and get output

Loading cuts ROOT file blackmore_data/run7469_extracted_data__MAUS2pt1b.root
0 0 -2.18275327939e-60 0.0
1 1 -0.000114954425953 4.94065645841e-324
2 1 -0.000114954425953 9.88131291682e-324
3 1 -2.18275327939e-60 1.48219693752e-323
4 1 -3.29014704019e+110 1.97626258336e-323
5 1 -1.30120806415e+197 2.47032822921e-323
6 1 -5.38012113181e+68 2.96439387505e-323
7 1 -0.000114954425953 3.45845952089e-323
8 1 -0.000114954425953 3.95252516673e-323
9 1 -0.000114954425953 4.44659081257e-323
10 1 -2.18275327939e-60 4.94065645841e-323
11 0 2.95821624332e+18 5.43472210425e-323

The first column is entry number, second column is cut_allData which looks okay, third and fourth column are spill/event number which look messed up. Probably I am doing something silly...

Nope, I just assigned ints as doubles when declaring the spill and event number branches (who'd have thought it'd go THAT crazy... silly mistake nonetheless.)

Here's the fixed one, which works fine with your script. MC one will be along shortly.

Here's the MC with the fixed spill & event number branches. I've also removed the old root files, so we'll hopefully not get confused.

Adds new branches for:

- Pulls at upstream tracker stations: TKU_s1_pull, TKU_s2_pull, TKU_s3_pull, TKU_s4_pull, TKU_s5_pull
- Errors on values reported by scifi track point: TKU_s1_x_error, TKU_s1_y_error, TKU_s1_px_error, TKU_s1_py_error, TKU_s1_kappa_error (and similar for stations 2--5 ). The branches for TKU_s*_kappa_error may be replaced with an error on Pz in the future, but right now just returning what the track point tells me.

No changes to the cuts, so no new pass/fail cut plots.

Updating data file with correct electron time-of-flight peak information

(When I added the tof-tracker momentum cut, Scott's suggested electron TOF was 25.95ns, whereas it's more like 25.48 -- I used Scott's TOF when starting to implement his cuts, but should have stuck with the 'real' value -- this restores the real value)

Adding latest MC after running it through this program along with its pass/fail cut plots.

There's something odd about the behaviour at TOF0 in the horizontal which is throwing off the Rayner reconstruction (compare x pixel positions at TOF0 to 'x' according to the PMT times). The distribution becomes very offset, despite using the same calibration files. TOF0 y and TOF1 x and y look believable. (But without believing TOF0 xPixel --> x, the rest of the Rayner reconstruction isn't believable)

--> Any thoughts, Durga?

Updated github repository: https://github.com/victoriablackmore/Extract-MAUS-data

-- Now includes Scott's muon mass calculation cut in 'cut_allPassed'. (Though I've corrected the equation to (p_track + p_correction)/(beta*gamma))
-- Simplified cuts for whether tracks hit 1--5 tracker stations, using pattern recognitions get_num_points().
-- Removed some now redundant cut variables

-- cut_allPassed now actually includes passing all cuts. This is true if and only if we have a good convergence on the Mark Rayner TOF reconstruction, have a muon-compatible time of flight, hit TOF0, TOF1 & 5 stations of TKU, have only one hit in TOF0, TOF1 and TKU, have a Pvalue > 0.01 and a mass between 90 & 120MeV.

-- run7469_extracted_data__MAUS2pt3pt1_b.root is the latest extracted root file from MAUS v2.3.1 (2.4 to come soon once I've done updating..)
-- cut_results.zip is an updated plotting script and the pass/fail plots for each individual cut and the cumulative 'cut_allPassed'. Note that since I've removed some redundant cuts, there's now fewer of them. Cuts are:

_noCut: Fill into histogram, regardless of whether it passes any cuts
        _cut1: Fill histogram if PMT-timed position at TOF0 is OK
        _cut2: Fill histogram if PMT-timed position at TOF1 is OK
        _cut3: Fill histogram if Rayner Reconstruction at TOFs is OK
        _cut4: Fill histogram if time-of-flight is OK
        _cut5: Fill histogram if all 5 tracker planes were hit
        _cut6: Fill histogram if all detectors hit
        _cut7: Fill histogram if there was only one spacepoint at TOF0
        _cut8: Fill histogram if there was only one spacepoint at TOF1
        _cut9: Fill histogram if there was only one track in TKU
        _cut10: Fill histogram if particle mass is OK
        _cut11: Fill histogram if Pvalue is OK
        _cut12: Fill if all cuts are passed

(a mean calculated muon mass for cut_allPassed of 105.4 isn't too shabby...)

I've tightened up the mass cut to accepting between 99 and 112MeV, which cleans up the tof-tracker momentum plot quite nicely (also attached).

A new extracted data file for MAUS 2.5. Plots for pass/fail cuts coming shortly(ish).

Attached is the tracker tracks extrapolated to TOF1. The file name holds version information i.e. geometry #160, MAUS v2.5.0 for recon, and version 1 (in case we find a bug and I need to do an update)

Attaching MC root file for MAUS v2.5

Updating extracted data file:

• cut_allPassed has been modified to remove the mass cut (added in the last update) and replace it with a cut on the momentum lost between TOF1 and TKU.
  • removed: cut_muon_mass branch
  • added: cut_momentum_loss branch
  • modified: cut_allPassed now depends on cut_momentum_loss and not cut_muon_mass

MC file will follow shortly.

Cut plots shown at CM45 for run 7469 data file: run7469_extracted_data__MAUS2pt5_b.root (zip file includes the .cpp used to create them)

Updating with MC, using a MC calibration for the TOF positioning by PMT

• run7417_calibration_file_MC.txt -- speed of light in scintillator + cable length corrections for MC (derived from run 7469, contrary to file name)

Attached new version for tracker at tof1; now includes tracks extrapolated to the diffuser. The z position is in the listing (so you can see which track is where). Let me know if you need anything else.

Attached extracted MC file; my routines may have gotten confused by the fact that all spill numbers are 1. I think this arises due to the way MC is stitched together from many reconstructions each with 1 spill...

Victoria asked me to run the track extrapolation script for run 8590 real data; extrapolated data files attached. Victoria, I have made some quite significant changes to the bureaucracy in my script - so can you check that you are happy (and let me know if you are not). Thanks!

Attached the extrapolation for 8590 mc. Note I only include the first 10 MC files in here, i.e. a glob like 000025/0000?/*_sim.root - I can run over all 1000 files or so, but grateful to hear back about whether the data is useful first!

Attached extrapolated tracks for the MC data set. I think the 7469 geometry is the right one to use here. I had a slightly modified set of z positions for diffuser and TOF, based on input from Victoria. All of my spill numbers are 1, as this is what is recorded in the originating file. I can force to increment the spill number or something if required.

Attached extrapolated tracks mc_3mm200_07469_MAUS-v284_2_extrapolated_tracks.csv.gz from following data set:

http://reco.mice.rl.ac.uk/MAUS-v2.8.4/MC/mc_3mm200_07469_MAUS-v284_2.tar

extrapolated using geometry 160 downloaded for run 07469

For some reason it failed to upload. Let me try again...

Last file got corrupted during upload. Let's try again.