Fork me on GitHub

Changeset 641cb3d in git


Ignore:
Timestamp:
May 18, 2016, 11:25:44 AM (8 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
5658083
Parents:
786028a
Message:

removed couts

File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/VertexFinder.cc

    r786028a r641cb3d  
    121121      if (!clusterIDToInt.at (cluster->first).at ("ndf"))
    122122        continue;
     123       
    123124      // Grow the cluster if GrowSeeds is true
    124125      if (fGrowSeeds)
     
    127128      // If the cluster still has fewer than MinNDF tracks, release the tracks;
    128129      // otherwise, mark the seed track as claimed
     130     
    129131      if ((int) clusterIDToInt.at (cluster->first).at ("ndf") < fMinNDF)
    130132        {
     
    162164  for (map<unsigned, map<string, int> >::const_iterator cluster = clusterIDToInt.begin (); cluster != clusterIDToInt.end (); cluster++)
    163165  {
     166   
    164167    if (cluster->second.at ("ndf") < fMinNDF)
    165168      continue;
     
    167170  }
    168171  sort (clusterSumPT2.begin (), clusterSumPT2.end (), secondDescending);
     172 
    169173  for (vector<pair<unsigned, double> >::const_iterator cluster = clusterSumPT2.begin (); cluster != clusterSumPT2.end (); cluster++)
    170174  {
     
    272276      if (!nearTracks.size ())
    273277        {
    274           for (map<unsigned, map<string, double> >::const_iterator track = trackIDToDouble.begin (); track != trackIDToDouble.end (); track++)
     278       
     279            for (map<unsigned, map<string, double> >::const_iterator track = trackIDToDouble.begin (); track != trackIDToDouble.end (); track++)
    275280            {
    276281              if (trackIDToBool.at (track->first).at ("claimed") || trackIDToInt.at (track->first).at ("clusterIndex") == (int) clusterIndex)
    277282                continue;
     283               
     284              Double_t sz_tr =  track->second.at ("ez") *  track->second.at ("z");
     285              Double_t sz_vt =  clusterIDToDouble.at (clusterIndex).at ("ez") *  clusterIDToDouble.at (clusterIndex).at ("z");
     286               
    278287              double distance = fabs (clusterIDToDouble.at (clusterIndex).at ("z") - track->second.at ("z")) / hypot (clusterIDToDouble.at (clusterIndex).at ("ez"), track->second.at ("ez"));
     288              //double distance = fabs (clusterIDToDouble.at (clusterIndex).at ("z") - track->second.at ("z")) / hypot (sz_vt, sz_tr);
    279289              if (nearestDistance < 0.0 || distance < nearestDistance)
    280290                {
     
    286296            }
    287297        }
     298     
    288299      else
    289300        {
     
    300311            }
    301312        }
    302 
     313     
    303314      // If no tracks within Sigma of the cluster were found, stop growing.
    304315      done = nearestDistance > fSigma || nearestDistance < 0.0;
Note: See TracChangeset for help on using the changeset viewer.