Skip to content

Patch for ATLASRECTS-6442. Do not extrapolate to calo along momentum if we are already past it

Christos Anastopoulos requested to merge ATLAS-EGamma/athena:ATLASRECTS-6442_debug into master

This allows the event discussed in ATLASRECTS-6442 to be processed

it will print this warning

AWtoALL 18:50:03 LRTEMGSFCaloExtensionBui...EMLastCaloExtensionTool     0     0 WARNING Can extrapolate to calo along momentum if already past it. Probematic parameters : CurvilinearParametersT p      arameters:
14992 RAWtoALL 18:50:03  * TrackParameters on Surface
14993 RAWtoALL 18:50:03  * loc1  : 0.0000000
14994 RAWtoALL 18:50:03  * loc2  : 0.0000000
14995 RAWtoALL 18:50:03  * phi   : -1.9496552
14996 RAWtoALL 18:50:03  * Theta : 0.2369496
14997 RAWtoALL 18:50:03  * q/p   : 0.0067803
14998 RAWtoALL 18:50:03  * charge: 1.0000000
14999 RAWtoALL 18:50:03  * covariance matrix = 0x7fe5fbf3afa0
15000 RAWtoALL 18:50:03  * corresponding global parameters:
15001 RAWtoALL 18:50:03  *    position  (x,  y,  z ) = (-197.5742493, -7267.3413086, 2739.0637207)
15002 RAWtoALL 18:50:03  *    momentum  (px, py, pz) = (-12.8048811, -32.1657829, 143.3658295)
15003 RAWtoALL 18:50:03 associated surface:
15004 RAWtoALL 18:50:03 Trk::PlaneSurface
15005 RAWtoALL 18:50:03      Center position  (x, y, z) = (-197.5742, -7267.3413, 2739.0637)
15006 RAWtoALL 18:50:03      Rotation:             colX = (0.929087, -0.369860, 0.000000)
15007 RAWtoALL 18:50:03                            colY = (0.359526, 0.903127, 0.234739)
15008 RAWtoALL 18:50:03                            colZ = (-0.086821, -0.218093, 0.972059)
15009 RAWtoALL 18:50:03      Bounds  : Trk::NoBounds ... boundless surface

rather than go to an infinite loop.

Still will be good to not have these cases but this will need more time. So for now assume is rare , print the WARNING .

The fundamental issue is that you can not extrapolate to calo along momentum if you are past it. The even more fundamental issue is how we end up with what should be an ID track particle with last measurement in the MS ?

Merge request reports