Fork me on GitHub

Changeset e55f5b0 in git for modules


Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/ParticlePropagator.cc

    r443f093 re55f5b0  
    100100  {
    101101    fBeamSpotInputArray = 0;
    102   } 
     102  }
    103103  // create output arrays
    104104
     
    133133  Double_t l, d0, dz, p, ctgTheta, phip, etap, alpha;
    134134  Double_t bsx, bsy, bsz;
    135          
     135
    136136  const Double_t c_light = 2.99792458E8;
    137  
    138   if (!fBeamSpotInputArray || fBeamSpotInputArray->GetSize () == 0) 
     137
     138  if (!fBeamSpotInputArray || fBeamSpotInputArray->GetSize () == 0)
    139139    beamSpotPosition.SetXYZT(0.0, 0.0, 0.0, 0.0);
    140140  else
     
    143143    beamSpotPosition = beamSpotCandidate.Position;
    144144  }
    145  
     145
    146146  fItInputArray->Reset();
    147147  while((candidate = static_cast<Candidate*>(fItInputArray->Next())))
     
    152152    y = candidatePosition.Y()*1.0E-3;
    153153    z = candidatePosition.Z()*1.0E-3;
    154    
     154
    155155    bsx = beamSpotPosition.X()*1.0E-3;
    156156    bsy = beamSpotPosition.Y()*1.0E-3;
    157157    bsz = beamSpotPosition.Z()*1.0E-3;
    158    
     158
    159159    q = candidate->Charge;
    160160
     
    206206      y_t = y + py*t;
    207207      z_t = z + pz*t;
    208      
     208
    209209      l = TMath::Sqrt( (x_t - x)*(x_t - x) + (y_t - y)*(y_t - y) + (z_t - z)*(z_t - z));
    210210
     
    215215      candidate->Position.SetXYZT(x_t*1.0E3, y_t*1.0E3, z_t*1.0E3, candidatePosition.T() + t*e*1.0E3);
    216216      candidate->L = l*1.0E3;
    217    
     217
    218218      candidate->Momentum = candidateMomentum;
    219219      candidate->AddCandidate(mother);
     
    267267      yd = (rc2 > 0.0) ? yd / rc2 : -999;
    268268      zd = z + (TMath::Sqrt(xd*xd + yd*yd) - TMath::Sqrt(x*x + y*y))*pz/pt;
    269      
     269
    270270      // use perigee momentum rather than original particle
    271       // momentum, since the orignal particle momentum isn't known 
    272      
     271      // momentum, since the orignal particle momentum isn't known
     272
    273273      px = TMath::Sign(1.0,r) * pt * (-y_c / r_c);
    274274      py = TMath::Sign(1.0,r) * pt * (x_c / r_c);
    275275      etap = candidateMomentum.Eta();
    276276      phip = TMath::ATan2(py, px);
    277        
     277
    278278      candidateMomentum.SetPtEtaPhiE(pt, etap, phip, candidateMomentum.E());
    279      
     279
    280280      // calculate additional track parameters (correct for beamspot position)
    281        
     281
    282282      d0        = (  (x - bsx) * py - (y - bsy) * px) / pt;
    283283      dz        = z - ((x - bsx) * px + (y - bsy) * py) / pt * (pz / pt);
    284284      p         = candidateMomentum.P();
    285285      ctgTheta  = 1.0 / TMath::Tan (candidateMomentum.Theta ());
    286          
     286
     287
    287288      // 3. time evaluation t = TMath::Min(t_r, t_z)
    288289      //    t_r : time to exit from the sides
     
    330331      r_t = TMath::Hypot(x_t, y_t);
    331332
    332      
     333
    333334      // compute path length for an helix
    334      
     335
    335336      alpha = pz*1.0E9 / c_light / gammam;
    336337      l = t * TMath::Sqrt(alpha*alpha + r*r*omega*omega);
    337          
     338
    338339      if(r_t > 0.0)
    339340      {
    340        
     341
    341342        // store these variables before cloning
    342            
    343343        candidate->D0 = d0*1.0E3;
    344344        candidate->DZ = dz*1.0E3;
     
    347347        candidate->CtgTheta = ctgTheta;
    348348        candidate->Phi = phip;
    349      
     349
    350350        mother = candidate;
    351351        candidate = static_cast<Candidate*>(candidate->Clone());
     
    355355
    356356        candidate->Momentum = candidateMomentum;
    357        
     357
    358358            candidate->L  =  l*1.0E3;
    359        
     359
    360360            candidate->Xd = xd*1.0E3;
    361361        candidate->Yd = yd*1.0E3;
Note: See TracChangeset for help on using the changeset viewer.