Changeset e55f5b0 in git for modules/ParticlePropagator.cc

Timestamp:

Aug 23, 2016, 3:20:20 PM (8 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

91acc76

Parents:

443f093

Message:

remove empty spaces

File:

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

modules/ParticlePropagator.cc

@@ -100 +100 @@
   {
     fBeamSpotInputArray = 0;
-  }  
+  }
   // create output arrays
 
@@ -133 +133 @@
   Double_t l, d0, dz, p, ctgTheta, phip, etap, alpha;
   Double_t bsx, bsy, bsz;
-         
+
   const Double_t c_light = 2.99792458E8;
- 
-  if (!fBeamSpotInputArray || fBeamSpotInputArray->GetSize () == 0) 
+
+  if (!fBeamSpotInputArray || fBeamSpotInputArray->GetSize () == 0)
     beamSpotPosition.SetXYZT(0.0, 0.0, 0.0, 0.0);
   else
@@ -143 +143 @@
     beamSpotPosition = beamSpotCandidate.Position;
   }
- 
+
   fItInputArray->Reset();
   while((candidate = static_cast<Candidate*>(fItInputArray->Next())))
@@ -152 +152 @@
     y = candidatePosition.Y()*1.0E-3;
     z = candidatePosition.Z()*1.0E-3;
-    
+
     bsx = beamSpotPosition.X()*1.0E-3;
     bsy = beamSpotPosition.Y()*1.0E-3;
     bsz = beamSpotPosition.Z()*1.0E-3;
-   
+
     q = candidate->Charge;
 
@@ -206 +206 @@
       y_t = y + py*t;
       z_t = z + pz*t;
-      
+
       l = TMath::Sqrt( (x_t - x)*(x_t - x) + (y_t - y)*(y_t - y) + (z_t - z)*(z_t - z));
 
@@ -215 +215 @@
       candidate->Position.SetXYZT(x_t*1.0E3, y_t*1.0E3, z_t*1.0E3, candidatePosition.T() + t*e*1.0E3);
       candidate->L = l*1.0E3;
-   
+
       candidate->Momentum = candidateMomentum;
       candidate->AddCandidate(mother);
@@ -267 +267 @@
       yd = (rc2 > 0.0) ? yd / rc2 : -999;
       zd = z + (TMath::Sqrt(xd*xd + yd*yd) - TMath::Sqrt(x*x + y*y))*pz/pt;
-      
+
       // use perigee momentum rather than original particle
-      // momentum, since the orignal particle momentum isn't known 
-      
+      // 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);
       etap = candidateMomentum.Eta();
       phip = TMath::ATan2(py, px);
-       
+
       candidateMomentum.SetPtEtaPhiE(pt, etap, phip, candidateMomentum.E());
-      
+
       // calculate additional track parameters (correct for beamspot position)
-       
+
       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 ());
-          
+
+
       // 3. time evaluation t = TMath::Min(t_r, t_z)
       //    t_r : time to exit from the sides
@@ -330 +331 @@
       r_t = TMath::Hypot(x_t, y_t);
 
-      
+
       // compute path length for an helix
-      
+
       alpha = pz*1.0E9 / c_light / gammam;
       l = t * TMath::Sqrt(alpha*alpha + r*r*omega*omega);
-          
+
       if(r_t > 0.0)
       {
-        
+
         // store these variables before cloning
-           
         candidate->D0 = d0*1.0E3;
         candidate->DZ = dz*1.0E3;
@@ -347 +347 @@
         candidate->CtgTheta = ctgTheta;
         candidate->Phi = phip;
-      
+
         mother = candidate;
         candidate = static_cast<Candidate*>(candidate->Clone());
@@ -355 +355 @@
 
         candidate->Momentum = candidateMomentum;
-        
+
             candidate->L  =  l*1.0E3;
-        
+
             candidate->Xd = xd*1.0E3;
         candidate->Yd = yd*1.0E3;