Skip to content
Snippets Groups Projects

Clean-up, resurrection MC-Checker, all configurables in

7 files
+ 19
23
Compare changes
  • Side-by-side
  • Inline
Files
7
@@ -399,12 +399,10 @@ void PrPixelTracking::doPairSearch( std::vector<PrPixelModuleHits>& modulehits,
if( x1*( x0-x1)*invdz + y1*( y0-y1)*invdz <0. ) continue;
}
if ( configuration == Conf::Backward){
if( configuration == Conf::Backward){
// info()<<"Continuing backward -----"<<endmsg;
if(x1*( x0-x1)*invdz + y1*( y0-y1)*invdz >0.) continue;
}
//zCross Cut not applied
if(x1*( x0-x1)*invdz + y1*( y0-y1)*invdz >0.) continue;
}
// Make a seed track out of these two hits.
// Extend the seed track towards smaller Z.
hitbuffer[1] = sens0;
@@ -415,7 +413,6 @@ void PrPixelTracking::doPairSearch( std::vector<PrPixelModuleHits>& modulehits,
extendTrack<configuration>(modulehits, hitbuffer );
if ( hitbuffer[0] < 6 ) continue;
//if track has 3 hits && !m_hardFlagging
if ( hitbuffer[0] == 6 && !m_hardFlagging){
ThreeHitVec.insert(ThreeHitVec.end(), hitbuffer.begin()+1, hitbuffer.begin() + 1 + hitbuffer[0]);
@@ -617,6 +614,7 @@ size_t PrPixelTracking::bestHit( std::vector<PrPixelModuleHits>& modulehits, con
const float z1 = m1.getZ(hit1);
//info() << "BestHit start: " << x0 << " " << m0.getY(hit0) << " " << z0 << " " << x1 << " " << m1.getY(hit1) << " " << z1 << endmsg;
const float td = 1.0 / (z1 - z0);
const float txn = (x1 - x0);
const float tx = txn * td;
@@ -679,6 +677,7 @@ size_t PrPixelTracking::bestHit( std::vector<PrPixelModuleHits>& modulehits, con
if (fabs(dy) > m_extraTol) continue;
//info() << hit_x << endmsg;
if(m_hardFlagging && m2.IsUsed(i)) continue;
const float scatterDenom = 1.0 / (hit_z - z1);
const float dx = xPred - hit_x;
const float scatterNum = (dx * dx) + (dy * dy);
Loading