Project

General

Profile

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

  1. Download ExtraVirtualPlanes.dat
  2. Place it in the current working directory
  3. Modify /path/to/geometry/ParentGeometryFile.dat. Instead of
                            Module 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 like
                            Module 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
                            {
                            }
    
  4. Open ExtraVirtualPlanes.dat file
  5. 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