Actions
2018-01-05 emittance¶
11:00 GMT on Friday 5th January
http://mice.iit.edu/phonebridge.html
Agenda:- Status of paper (Blackmore)
- Status of tracker work (Hunt)
Instructions to add global track matching to recon¶
Prerequisites¶
- You will need a reasonably recent version of MAUS
- You will need a geometry downloaded in the usual way
Geometry Modifications¶
The aim is to add extra virtual planes that tells track matching where it should read out
- Download ExtraVirtualPlanes.dat
- Place it in the current working directory
- Modify
/path/to/geometry/ParentGeometryFile.dat
. Instead ofModule Virtuals { Volume None Position 0.0 0.0 0.0+0.5*@RepeatNumber m PropertyString SensitiveDetector Virtual PropertyBool RepeatModule2 True PropertyInt NumberOfRepeats 48 }
you should add a few lines likeModule Virtuals { Volume None Position 0.0 0.0 0.0+0.5*@RepeatNumber m PropertyString SensitiveDetector Virtual PropertyBool RepeatModule2 True PropertyInt NumberOfRepeats 48 } Module ExtraVirtualPlanes.dat { }
- Open ExtraVirtualPlanes.dat file
- Feel free to modify ExtraVirtualPlanes.dat to account for the actual diffuser position; Victoria may be able to advise you on the correct position (it should be close, but will need modifying for survey etc)
Data cards modifications¶
Add the following data cards
track_matching_tolerances = { "TOF0t":1e6, # ns between actual and expected TOF0-1 Delta t "TOF1x":1000.0, "TOF1y":1000.0, "TOF2x":1000.0, "TOF2y":1000.0, "KLy":1000.0, "EMRx":1000.0, # Multiplier for the standard tolerance which is the reconstructed error*sqrt(12) "EMRy":1000.0, "TOF12maxSpeed":1.0, # fraction of c to calculate travel time between TOFs for through matching "TOF12minSpeed":0.5, "TKDpos":1e6, # position "TKDp":1e6, # momentum } # kEPlus, kEMinus, kMuPlus, kMuMinus, kPiPlus, kPiMinus track_matching_pid_hypothesis = "kMuPlus" track_matching_through_matching_logic = "propagate" track_matching_geometry_algorithm = "axial"
Run script modifications¶
Add some lines to your analyze_data_offline.py
script after all other reconstruction mappers. Instead of
my_map.append(MAUS.MapCppEMRPlaneHits()) my_map.append(MAUS.MapCppEMRSpacePoints()) my_map.append(MAUS.MapCppEMRRecon())
do
my_map.append(MAUS.MapCppEMRPlaneHits()) my_map.append(MAUS.MapCppEMRSpacePoints()) my_map.append(MAUS.MapCppEMRRecon()) my_map.append(MAUS.MapCppGlobalReconImport()) my_map.append(MAUS.MapCppGlobalTrackMatching())
Updated by Rogers, Chris over 5 years ago · 3 revisions