Changeset a0431dc in git for modules/ParticlePropagator.cc

Timestamp:

Apr 16, 2014, 4:08:33 PM (11 years ago)

Author:

mselvaggi <mselvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

9687203

Parents:

64a4950

Message:

track counting btagging

File:

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

modules/ParticlePropagator.cc

@@ -106 +106 @@
   Double_t tmp, discr, discr2;
   Double_t delta, gammam, omega, asinrho;
-
+  Double_t ang_mom, rcu, rc2, dxy, xd, yd, zd;
+  
   const Double_t c_light = 2.99792458E8;
 
@@ -214 +215 @@
       if(x_c < 0.0) phi += TMath::Pi();
 
+      rcu = TMath::Abs(r);
+      rc2 = r_c*r_c;
+     
+      // calculate coordinates of closest approach to track circle in transverse plane xd, yd, zd
+      xd = x_c*x_c*x_c - x_c*rcu*r_c + x_c*y_c*y_c; 
+      xd  = ( rc2 > 0.0 ) ? xd / rc2 : -999;
+      yd  = y_c*(-rcu*r_c + rc2);
+      yd  = ( rc2 > 0.0 ) ? yd / rc2 : -999;
+      zd  = z + (TMath::Sqrt(xd*xd+yd*yd) - TMath::Sqrt(x*x+y*y))*pz/pt;
+
+      // calculate impact paramater
+      ang_mom = (xd*py - yd*px);
+      dxy = ang_mom/pt;
+    
+         
       // 3. time evaluation t = TMath::Min(t_r, t_z)
       //    t_r : time to exit from the sides
@@ -267 +283 @@
 
         candidate->Momentum = candidateMomentum;
-        candidate->AddCandidate(mother);
+        candidate->Xd = xd*1.0E3; 
+        candidate->Yd = yd*1.0E3;
+        candidate->Zd = zd*1.0E3;
+        
+        candidate->AddCandidate(mother);
 
         fOutputArray->Add(candidate);
@@ -287 +307 @@
 
 //------------------------------------------------------------------------------
+