diff --git a/Validation/inc/ValKalSeed.hh b/Validation/inc/ValKalSeed.hh index 0c9c753e97..4b32c52fd4 100644 --- a/Validation/inc/ValKalSeed.hh +++ b/Validation/inc/ValKalSeed.hh @@ -43,6 +43,7 @@ class ValKalSeed { TH1D* _hpT; TH1D* _hpce; TH1D* _hpcep; + TH1D* _hsignedp; TH1D* _hpe; TH1D* _hRho; TH1D* _hPhi; diff --git a/Validation/src/ValKalSeed.cc b/Validation/src/ValKalSeed.cc index b9f9e9444a..8704f73b54 100644 --- a/Validation/src/ValKalSeed.cc +++ b/Validation/src/ValKalSeed.cc @@ -37,6 +37,7 @@ namespace mu2e { _hpT = tfs.make("pT", "p TPR", 100, 0., 110.); _hpce = tfs.make("pce", "p CE", 100, 95.0, 110.); _hpcep = tfs.make("pcep", "p CE+", 100, 82.0, 97.); + _hsignedp = tfs.make("signedp", "signedp", 200, -110., 110.); _hpe = tfs.make("pe", "p error", 100, 0.0, 1.0); _hRho = tfs.make("rho", "Transverse radius", 100, 0.0, 800.); _hPhi = tfs.make("phi", "phi", 100, -M_PI, M_PI); @@ -115,6 +116,7 @@ namespace mu2e { double p_mc = mcTrkP(event,vdid,p_pri); SurfaceId sid = _vdmap[vdid]; auto ikinter = ks.intersection(sid); + double ksCharge = ks.intersections().front().pstate_.charge(); if(ikinter != ks.intersections().end()){ auto mom3 = ikinter->momentum3(); double p = mom3.R(); @@ -124,6 +126,7 @@ namespace mu2e { if (isTPR) _hpT->Fill(p); _hpce->Fill(p); _hpcep->Fill(p); + _hsignedp->Fill(p*ksCharge); _hpe->Fill(ikinter->momerr()); _hRho->Fill(ikinter->position3().Rho()); _hPhi->Fill(mom3.Phi());