Changeset a07b54c in git for modules/ParticlePropagator.cc

Timestamp:

Apr 14, 2020, 7:45:19 PM (5 years ago)

Author:

Roberto Preghenella <preghenella@…>

Branches:

master

Children:

974f5bc

Parents:

09ffaee

git-author:

Roberto Preghenella <preghenella@…> (04/05/20 19:17:47)

git-committer:

Roberto Preghenella <preghenella@…> (04/14/20 19:45:19)

Message:

Proper calculation of the DCAz coordinate

File:

• modules/ParticlePropagator.cc (modified) (2 diffs)

modules/ParticlePropagator.cc

@@ -135 +135 @@
   Double_t l, d0, dz, p, ctgTheta, phip, etap, alpha;
   Double_t bsx, bsy, bsz;
-
+  Double_t s0, s1, sd;
+  
   const Double_t c_light = 2.99792458E8;
 
@@ -297 +298 @@
       zd = z + (TMath::Sqrt(xd * xd + yd * yd) - TMath::Sqrt(x * x + y * y)) * pz / pt;
 
+      // proper calculation of the DCAz coordinate
+      // s0: track circle parameter at the track origin
+      // s1: track circle parameter at the closest approach to beam pipe
+      // sd: s1-s0 signed angular difference
+      s0 = atan2(y - y_c, x - x_c);
+      s1 = atan2(yd - y_c, xd - x_c);
+      sd = atan2(sin(s1 - s0), cos(s1 - s0));
+      zd = z - r * pz / pt * sd;
+      
       // use perigee momentum rather than original particle
       // momentum, since the orignal particle momentum isn't known