Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Riemann fit multiple scattering #174

Conversation

rovere
Copy link

@rovere rovere commented Sep 25, 2018

This PR implements the multiple scattering treatments in the Riemann Fit. In particular:

  • It modifies the previous implementation of the multiple scattering in the circle fit to correctly cover both the barrel and the forward case
  • It implements the multiple scattering in the line fit in the S-Z plane both for the barrel and the forward case.

The effective radiation length is still an approximate value since the phi angle is not taken into account (it is not known on a layer-by-layer case). Ad ad-hoc correction based on the inverse of the pt has been added, with a cut-off of 1 GeV.

The pulls are ok-ish, the material could be further tuned.
The Chi2 is flat on all eta range.

The usual MTV collection of plots are available at this link for single muon of pt_1. I verified that the behaviour of the fit is ok-ish also using ttbar and singleMu pt_10.

@cmsbot
Copy link

cmsbot commented Sep 25, 2018

A new Pull Request was created by @rovere (Marco Rovere) for CMSSW_10_2_X_Patatrack.

It involves the following packages:

RecoPixelVertexing/PixelTrackFitting

@cmsbot, @fwyzard can you please review it and eventually sign? Thanks.

cms-bot commands are listed here

@fwyzard
Copy link

fwyzard commented Sep 25, 2018

Validation summary

Reference release CMSSW_10_2_4 at 6e8f631
Development branch CMSSW_10_2_X_Patatrack at 68a8252
Testing PRs:

makeTrackValidationPlots.py plots

/RelValTTbar_13/CMSSW_10_2_2-PU25ns_102X_upgrade2018_realistic_v11-v2/GEN-SIM-DIGI-RAW

/RelValZMM_13/CMSSW_10_2_2-102X_upgrade2018_realistic_v11-v1/GEN-SIM-DIGI-RAW

DQM GUI plots

/RelValTTbar_13/CMSSW_10_2_2-PU25ns_102X_upgrade2018_realistic_v11-v2/GEN-SIM-DIGI-RAW

/RelValZMM_13/CMSSW_10_2_2-102X_upgrade2018_realistic_v11-v1/GEN-SIM-DIGI-RAW

logs and nvprof/nvvp profiles

/RelValTTbar_13/CMSSW_10_2_2-PU25ns_102X_upgrade2018_realistic_v11-v2/GEN-SIM-DIGI-RAW

/RelValZMM_13/CMSSW_10_2_2-102X_upgrade2018_realistic_v11-v1/GEN-SIM-DIGI-RAW

Logs

The full log is available at https://fwyzard.web.cern.ch/fwyzard/patatrack/pulls/fb1828c2c01ab1c354daa40f9a056b39ec072dc1/log .

@cmsbot
Copy link

cmsbot commented Sep 25, 2018

Pull request #174 was updated. @cmsbot, @fwyzard can you please check and sign again.

@fwyzard
Copy link

fwyzard commented Sep 25, 2018

@rovere can you check what broke with the 10824.9 workflow (Riemann fit on GPU) ?

@rovere
Copy link
Author

rovere commented Sep 25, 2018

Sure, I'll have a look.

@fwyzard
Copy link

fwyzard commented Sep 25, 2018

Looking at TTbar Riemann fit workflow on the CPU (10824.7)

  • resolutions are unchanged or improved (d_z)
  • pulls are more in line with the standard workflow (10824.5), especially in the core of the distribution
  • χ²/ndf is much flatter: image

@fwyzard
Copy link

fwyzard commented Sep 25, 2018

The 10824.9 workflow fails with

RecoPixelVertexing/PixelTrackFitting/plugins/PixelTrackReconstructionGPU.cc, line 118: cudaErrorLaunchFailure: unspecified launch failure

 114      // LAUNCH THE KERNEL FIT
 115      launchKernelFit(hits_and_covariancesGPU, 12*4*total_seeds, 4,
 116          bField, helix_fit_resultsGPU);
 117      // CUDA MEMCOPY DEVICE2HOST OF HELIX_FIT
 118      cudaCheck(cudaDeviceSynchronize());

As it is unlikely that the problem is in the cudaDeviceSynchronize(), I would guess at PixelTrackReconstructionGPU::launchKernelFit()

@fwyzard
Copy link

fwyzard commented Sep 25, 2018

cuda-gdb points to KernelLineFitAllHits():

Thread 1 "cmsRun" received signal SIGTRAP, Trace/breakpoint trap.
[Switching focus to CUDA kernel 0, grid 8, block (0,0,0), thread (0,0,0), device 0, sm 14, warp 0, lane 0]
0x00007fff7cd24138 in KernelLineFitAllHits<<<(23,1,1),(32,1,1)>>> ()

@fwyzard
Copy link

fwyzard commented Sep 25, 2018

Also, is this orthogonal to #169 and #170 ?

@rovere
Copy link
Author

rovere commented Sep 25, 2018

Yes, I kept them separate

@fwyzard
Copy link

fwyzard commented Sep 25, 2018

Validation summary

Reference release CMSSW_10_2_4 at 6e8f631
Development branch CMSSW_10_2_X_Patatrack at a64cea6
Testing PRs:

🚧 Validation running at vinavx2.cern.ch:/data/user/fwyzard/patatrack/validation/run.LXXxfUzJi7 ...

@fwyzard
Copy link

fwyzard commented Sep 27, 2018

This should finally work on the GPU together with #178 and the latest patches to Eigen.

@fwyzard fwyzard merged commit 728bcab into cms-patatrack:CMSSW_10_2_X_Patatrack Sep 27, 2018
@fwyzard fwyzard added this to the CMSSW_10_2_5_Patatrack milestone Sep 27, 2018
fwyzard pushed a commit that referenced this pull request Oct 8, 2020
Implement the multiple scattering treatments in the Riemann Fit. In particular:
  - modify the previous implementation of the multiple scattering in the circle fit to correctly cover both the barrel and the forward case;
  - implement the multiple scattering in the line fit in the S-Z plane both for the barrel and the forward case.

The effective radiation length is still an approximate value since the phi angle is not taken into account (it is not known on a layer-by-layer case). Ad ad-hoc correction based on the inverse of the pt has been added, with a cut-off of 1 GeV.

The pulls are ok-ish, the material could be further tuned.
The Chi2 is flat on all eta range.
fwyzard pushed a commit that referenced this pull request Oct 20, 2020
Implement the multiple scattering treatments in the Riemann Fit. In particular:
  - modify the previous implementation of the multiple scattering in the circle fit to correctly cover both the barrel and the forward case;
  - implement the multiple scattering in the line fit in the S-Z plane both for the barrel and the forward case.

The effective radiation length is still an approximate value since the phi angle is not taken into account (it is not known on a layer-by-layer case). Ad ad-hoc correction based on the inverse of the pt has been added, with a cut-off of 1 GeV.

The pulls are ok-ish, the material could be further tuned.
The Chi2 is flat on all eta range.
fwyzard pushed a commit that referenced this pull request Oct 20, 2020
Implement the multiple scattering treatments in the Riemann Fit. In particular:
  - modify the previous implementation of the multiple scattering in the circle fit to correctly cover both the barrel and the forward case;
  - implement the multiple scattering in the line fit in the S-Z plane both for the barrel and the forward case.

The effective radiation length is still an approximate value since the phi angle is not taken into account (it is not known on a layer-by-layer case). Ad ad-hoc correction based on the inverse of the pt has been added, with a cut-off of 1 GeV.

The pulls are ok-ish, the material could be further tuned.
The Chi2 is flat on all eta range.
fwyzard pushed a commit that referenced this pull request Oct 23, 2020
Implement the multiple scattering treatments in the Riemann Fit. In particular:
  - modify the previous implementation of the multiple scattering in the circle fit to correctly cover both the barrel and the forward case;
  - implement the multiple scattering in the line fit in the S-Z plane both for the barrel and the forward case.

The effective radiation length is still an approximate value since the phi angle is not taken into account (it is not known on a layer-by-layer case). Ad ad-hoc correction based on the inverse of the pt has been added, with a cut-off of 1 GeV.

The pulls are ok-ish, the material could be further tuned.
The Chi2 is flat on all eta range.
fwyzard pushed a commit that referenced this pull request Nov 6, 2020
Implement the multiple scattering treatments in the Riemann Fit. In particular:
  - modify the previous implementation of the multiple scattering in the circle fit to correctly cover both the barrel and the forward case;
  - implement the multiple scattering in the line fit in the S-Z plane both for the barrel and the forward case.

The effective radiation length is still an approximate value since the phi angle is not taken into account (it is not known on a layer-by-layer case). Ad ad-hoc correction based on the inverse of the pt has been added, with a cut-off of 1 GeV.

The pulls are ok-ish, the material could be further tuned.
The Chi2 is flat on all eta range.
fwyzard pushed a commit that referenced this pull request Nov 6, 2020
Implement the multiple scattering treatments in the Riemann Fit. In particular:
  - modify the previous implementation of the multiple scattering in the circle fit to correctly cover both the barrel and the forward case;
  - implement the multiple scattering in the line fit in the S-Z plane both for the barrel and the forward case.

The effective radiation length is still an approximate value since the phi angle is not taken into account (it is not known on a layer-by-layer case). Ad ad-hoc correction based on the inverse of the pt has been added, with a cut-off of 1 GeV.

The pulls are ok-ish, the material could be further tuned.
The Chi2 is flat on all eta range.
fwyzard pushed a commit that referenced this pull request Nov 16, 2020
Implement the multiple scattering treatments in the Riemann Fit. In particular:
  - modify the previous implementation of the multiple scattering in the circle fit to correctly cover both the barrel and the forward case;
  - implement the multiple scattering in the line fit in the S-Z plane both for the barrel and the forward case.

The effective radiation length is still an approximate value since the phi angle is not taken into account (it is not known on a layer-by-layer case). Ad ad-hoc correction based on the inverse of the pt has been added, with a cut-off of 1 GeV.

The pulls are ok-ish, the material could be further tuned.
The Chi2 is flat on all eta range.
fwyzard added a commit that referenced this pull request Nov 27, 2020
Implement the multiple scattering treatments in the Riemann Fit. In particular:
  - modify the previous implementation of the multiple scattering in the circle fit to correctly cover both the barrel and the forward case;
  - implement the multiple scattering in the line fit in the S-Z plane both for the barrel and the forward case.

The effective radiation length is still an approximate value since the phi angle is not taken into account (it is not known on a layer-by-layer case). Ad ad-hoc correction based on the inverse of the pt has been added, with a cut-off of 1 GeV.

The pulls are ok-ish, the material could be further tuned.
The Chi2 is flat on all eta range.
fwyzard pushed a commit that referenced this pull request Dec 26, 2020
Implement the multiple scattering treatments in the Riemann Fit. In particular:
  - modify the previous implementation of the multiple scattering in the circle fit to correctly cover both the barrel and the forward case;
  - implement the multiple scattering in the line fit in the S-Z plane both for the barrel and the forward case.

The effective radiation length is still an approximate value since the phi angle is not taken into account (it is not known on a layer-by-layer case). Ad ad-hoc correction based on the inverse of the pt has been added, with a cut-off of 1 GeV.

The pulls are ok-ish, the material could be further tuned.
The Chi2 is flat on all eta range.
fwyzard pushed a commit that referenced this pull request Jan 15, 2021
Implement the multiple scattering treatments in the Riemann Fit. In particular:
  - modify the previous implementation of the multiple scattering in the circle fit to correctly cover both the barrel and the forward case;
  - implement the multiple scattering in the line fit in the S-Z plane both for the barrel and the forward case.

The effective radiation length is still an approximate value since the phi angle is not taken into account (it is not known on a layer-by-layer case). Ad ad-hoc correction based on the inverse of the pt has been added, with a cut-off of 1 GeV.

The pulls are ok-ish, the material could be further tuned.
The Chi2 is flat on all eta range.
fwyzard pushed a commit that referenced this pull request Mar 23, 2021
Implement the multiple scattering treatments in the Riemann Fit. In particular:
  - modify the previous implementation of the multiple scattering in the circle fit to correctly cover both the barrel and the forward case;
  - implement the multiple scattering in the line fit in the S-Z plane both for the barrel and the forward case.

The effective radiation length is still an approximate value since the phi angle is not taken into account (it is not known on a layer-by-layer case). Ad ad-hoc correction based on the inverse of the pt has been added, with a cut-off of 1 GeV.

The pulls are ok-ish, the material could be further tuned.
The Chi2 is flat on all eta range.
fwyzard pushed a commit that referenced this pull request Apr 1, 2021
Implement the multiple scattering treatments in the Riemann Fit. In particular:
  - modify the previous implementation of the multiple scattering in the circle fit to correctly cover both the barrel and the forward case;
  - implement the multiple scattering in the line fit in the S-Z plane both for the barrel and the forward case.

The effective radiation length is still an approximate value since the phi angle is not taken into account (it is not known on a layer-by-layer case). Ad ad-hoc correction based on the inverse of the pt has been added, with a cut-off of 1 GeV.

The pulls are ok-ish, the material could be further tuned.
The Chi2 is flat on all eta range.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants