Fork me on GitHub

Changeset 4ac0049 in git for modules/PileUpSubtractor4D.cc


Ignore:
Timestamp:
Jan 24, 2020, 3:51:00 PM (5 years ago)
Author:
Kaan Yüksel Oyulmaz <kaanyukseloyulmaz@…>
Branches:
Timing
Children:
79a7b3e
Parents:
6049672
Message:

VertexFinderDA4D.cc module is cleaned from graphs, FCChh_PileUpVtx.tcl card, PileUpSubtractor4D.cc and TimeSmearing.cc were updated in meeting

File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/PileUpSubtractor4D.cc

    r6049672 r4ac0049  
    170170      t_err = particle->PositionError.T();
    171171
    172       distanceCharged = pow((zvtx - z),2)/pow((zvtx_err - z_err),2) + pow((tvtx - t),2)/pow((tvtx_err - t_err),2);
    173       distanceNeutral = pow((tvtx - t),2)/pow((tvtx_err - t_err),2);
     172      distanceCharged = TMath::Sqrt(pow((zvtx - z),2)/pow((zvtx_err),2) + pow((tvtx - t),2)/pow((tvtx_err),2));
     173      distanceNeutral = TMath::Sqrt(pow((tvtx - t),2)/pow((tvtx_err),2));
    174174
    175175      if(candidate->Charge != 0 && distanceCharged < fChargedMinSignificance)
    176176      {
    177         candidate->IsRecoPU = 1;
     177        candidate->IsRecoPU = 0;
    178178      }
    179179      else if(candidate->Charge == 0 && distanceNeutral < fNeutralMinSignificance)
    180180      {
    181         candidate->IsRecoPU = 1;
     181        candidate->IsRecoPU = 0;
    182182      } 
    183183      else
    184184      {
    185         candidate->IsRecoPU = 0;
     185        candidate->IsRecoPU = 1;
    186186        if(candidate->Momentum.Pt() > fPTMin) array->Add(candidate);
    187187      }
Note: See TracChangeset for help on using the changeset viewer.