Linear Collider Forum



Home » Software Tools » org.lcsim » Track ErrorMatrix
Re: Track ErrorMatrix [message #442 is a reply to message #441] Tue, 21 March 2006 11:52 Go to previous messageGo to previous message
tonyj
Messages: 138
Registered: January 2004

Hi Jan, I just added a useful little utility to org.lcsim, which makes events available from the pnuts or jython console.

I opened a stdhep file, and loaded the FastMC example which comes with org.lcsim, stepped to the first event, opened a pnuts console and typed the following:

e = lcsim.getCurrentEvent()
tracks = e.get("Tracks")
mat = tracks.get(0).getErrorMatrix()


Here is the result:
> e = lcsim.getCurrentEvent()
Run 0 Event 1 (Tue Mar 21 11:26:42 PST 2006) Detector: sdjan03
> tracks = e.get("Tracks")
[ReconTrack Parameters:  d0          phi0          omega         z0        tan(lambda)

  0.004371123604 -3.094400588213 -0.000038320436 -0.009565853223 -3.634806948174

Error Matrix:

  0.000031970110 -0.000000367025  0.000000001058  0.000000005490 -0.000000000031
 -0.000000367025  0.000000007752 -0.000000000020 -0.000000000051  0.000000000001
  0.000000001058 -0.000000000020  0.000000000000  0.000000000000 -0.000000000000
  0.000000005490 -0.000000000051  0.000000000000  0.000056281122 -0.000001312831
 -0.000000000031  0.000000000001 -0.000000000000 -0.000001312831  0.000000051777

, ReconTrack Parameters:  d0          phi0          omega         z0        tan(lambda)

 -0.001509352603  0.357422921987  0.000033598599  0.000949223046 -1.127703899207

Error Matrix:

  0.000007383962 -0.000000058269  0.000000000045  0.000000001241 -0.000000000005
 -0.000000058269  0.000000001326 -0.000000000001 -0.000000000005  0.000000000000
  0.000000000045 -0.000000000001  0.000000000000  0.000000000000 -0.000000000000
  0.000000001241 -0.000000000005  0.000000000000  0.000008122098 -0.000000084347
 -0.000000000005  0.000000000000 -0.000000000000 -0.000000084347  0.000000002367

]
> mat = tracks.get(0).getErrorMatrix()
[[3.197011047506327E-5, -3.6702525661545337E-7, 1.057507040728786E-9, 5.4897903484976655E-9, -3.113593994858509E-11],
 [-3.6702525661545337E-7, 7.751717752114958E-9, -2.0461588842961317E-11, -5.145716975496867E-11, 7.302574334288315E-13],
 [1.057507040728786E-9, -2.0461588842961317E-11, 1.1597072584178297E-13, 1.6441955199527996E-13, -2.7366390939866518E-15],
 [5.4897903484976655E-9, -5.145716975496867E-11, 1.6441955199527996E-13, 5.6281122020588585E-5, -1.3128312841628168E-6], 
 [-3.113593994858509E-11, 7.302574334288315E-13, -2.7366390939866518E-15, -1.3128312841628168E-6, 5.1777470049936985E-8]]


so at least in this example the error matrix appears to be set fine.

Tony

[Updated on: Tue, 21 March 2006 17:32]

 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic:RefinedCluster
Next Topic:ReconTrack and associated MCParticle
Goto Forum:
  


Current Time: Sat Feb 22 06:09:52 Pacific Standard Time 2020
.:: Contact :: Home ::.

Powered by: FUDforum 3.0.1.
Copyright ©2001-2010 FUDforum Bulletin Board Software