Changeset 9330b7b in git

Timestamp:

Apr 3, 2017, 4:52:11 PM (7 years ago)

Author:

Pavel Demin <pavel-demin@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

accffb9

Parents:

781edf86

Message:

add RadiusMax and HalfLengthMax to ParticlePropagator

Location:

modules

Files:

• ParticlePropagator.cc (modified) (7 diffs)
• ParticlePropagator.h (modified) (1 diff)

modules/ParticlePropagator.cc

@@ -86 +86 @@
     return;
   }
+
+  fRadiusMax = GetDouble("RadiusMax", fRadius);
+  fHalfLengthMax = GetDouble("HalfLengthMax", fHalfLength);
 
   // import array with output from filter/classifier module
@@ -160 +163 @@
 
     // check that particle position is inside the cylinder
-    if(TMath::Hypot(x, y) > fRadius || TMath::Abs(z) > fHalfLength)
+    if(TMath::Hypot(x, y) > fRadiusMax || TMath::Abs(z) > fHalfLengthMax)
     {
       continue;
@@ -177 +180 @@
     }
 
-    if(TMath::Abs(q) < 1.0E-9 || TMath::Abs(fBz) < 1.0E-9)
+    if(TMath::Hypot(x, y) > fRadius || TMath::Abs(z) > fHalfLength)
+    {
+      mother = candidate;
+      candidate = static_cast<Candidate*>(candidate->Clone());
+
+      candidate->InitialPosition = candidatePosition;
+      candidate->Position = candidatePosition;
+      candidate->L = 0.0;
+
+      candidate->Momentum = candidateMomentum;
+      candidate->AddCandidate(mother);
+
+      fOutputArray->Add(candidate);
+    }
+    else if(TMath::Abs(q) < 1.0E-9 || TMath::Abs(fBz) < 1.0E-9)
     {
       // solve pt2*t^2 + 2*(px*x + py*y)*t - (fRadius2 - x*x - y*y) = 0
@@ -271 +288 @@
       // momentum, since the orignal particle momentum isn't known
 
-      px = TMath::Sign(1.0,r) * pt * (-y_c / r_c);
-      py = TMath::Sign(1.0,r) * pt * (x_c / r_c);
+      px = TMath::Sign(1.0, r) * pt * (-y_c / r_c);
+      py = TMath::Sign(1.0, r) * pt * (x_c / r_c);
       etap = candidateMomentum.Eta();
       phip = TMath::ATan2(py, px);
@@ -280 +297 @@
       // calculate additional track parameters (correct for beamspot position)
 
-      d0        = (  (x - bsx) * py - (y - bsy) * px) / pt;
+      d0        = ((x - bsx) * py - (y - bsy) * px) / pt;
       dz        = z - ((x - bsx) * px + (y - bsy) * py) / pt * (pz / pt);
       p         = candidateMomentum.P();
-      ctgTheta  = 1.0 / TMath::Tan (candidateMomentum.Theta ());
+      ctgTheta  = 1.0 / TMath::Tan (candidateMomentum.Theta());
 
 
@@ -301 +318 @@
       else
       {
-        asinrho = TMath::ASin( (fRadius*fRadius - r_c*r_c - r*r) / (2*TMath::Abs(r)*r_c)  );
+        asinrho = TMath::ASin((fRadius*fRadius - r_c*r_c - r*r) / (2*TMath::Abs(r)*r_c));
         delta = phi_0 - phi;
         if(delta <-TMath::Pi()) delta += 2*TMath::Pi();
@@ -356 +373 @@
         candidate->Momentum = candidateMomentum;
 
-            candidate->L  =  l*1.0E3;
-
-            candidate->Xd = xd*1.0E3;
+        candidate->L = l*1.0E3;
+
+        candidate->Xd = xd*1.0E3;
         candidate->Yd = yd*1.0E3;
         candidate->Zd = zd*1.0E3;

modules/ParticlePropagator.h

@@ -50 +50 @@
 private:
 
-  Double_t fRadius, fRadius2, fHalfLength;
+  Double_t fRadius, fRadius2, fRadiusMax, fHalfLength, fHalfLengthMax;
   Double_t fBz;