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

Geant4 Sensitive Detectors update 4: revised selection of particle type for shower libraries #22145

Merged
merged 5 commits into from Apr 22, 2018

Conversation

civanch
Copy link
Contributor

@civanch civanch commented Feb 7, 2018

This PR was discussed at the SIM meeting:
https://indico.cern.ch/event/702324/contributions/2885360/attachments/1596130/2528121/VI20180206.pdf

Several static methods are added to identify Geant4 particle type. These methods are used to select shower library branches with EM or hadronic hits. Hyperons and exotic particles will not any more treated as "hadrons", they should decay and only decay products will provide showers. Due to that, in many workflows one can expect different simulation histories and only statistical comparison of results is possible.

HF hits are sampled differently:

  1. for e+-, gamma and stable hadrons SL is called, if kinetic energy of a particle is below 50 MeV the particle is killed 2) Hyperons are ignored - no hits are created, allowing them to decay first
  2. If baryon number of a particle is above 0 then kinetic energy is used for SL interpolation and not total energy. DUe to these modifications low-energy ions does not produce any more extra hits

CastorSD class was changed more in order to use this opportunity when results are not identical:

  1. many commented lines removed; 2) fabs -> std::abs 3) early check of primary particle type and do not perform computations if a particle is not for shower library and cannot produce Cerenkov light.

Additionally some commented lines and unused variables are removed in other classes.

@cmsbuild
Copy link
Contributor

cmsbuild commented Feb 7, 2018

The code-checks are being triggered in jenkins.

@civanch
Copy link
Contributor Author

civanch commented Feb 7, 2018

please test

@cmsbuild
Copy link
Contributor

cmsbuild commented Feb 7, 2018

@cmsbuild
Copy link
Contributor

cmsbuild commented Feb 7, 2018

The tests are being triggered in jenkins.
https://cmssdt.cern.ch/jenkins/job/ib-any-integration/25948/console Started: 2018/02/07 19:42

@cmsbuild
Copy link
Contributor

cmsbuild commented Feb 7, 2018

A new Pull Request was created by @civanch (Vladimir Ivantchenko) for master.

It involves the following packages:

SimG4CMS/Calo
SimG4CMS/Forward
SimG4CMS/Tracker
SimG4Core/Notification

@cmsbuild, @civanch, @mdhildreth can you please review it and eventually sign? Thanks.
@makortel this is something you requested to watch as well.
@davidlange6, @slava77, @fabiocos you are the release manager for this.

cms-bot commands are listed here

@civanch
Copy link
Contributor Author

civanch commented Feb 7, 2018

@bsunanda , @kpedro88 - we discussed this at the meeting

@cmsbuild
Copy link
Contributor

cmsbuild commented Feb 7, 2018

@cmsbuild
Copy link
Contributor

cmsbuild commented Feb 7, 2018

Comparison job queued.

@cmsbuild
Copy link
Contributor

cmsbuild commented Feb 8, 2018

Comparison is ready
https://cmssdt.cern.ch/SDT/jenkins-artifacts/pull-request-integration/PR-22145/25948/summary.html

@slava77 comparisons for the following workflows were not done due to missing matrix map:

  • /build/cmsbld/jenkins/workspace/compare-root-files-short-matrix/results/JR-comparison/PR-22145/1325.7_TTbar_13_94XNanoAODINPUT+TTbar_13_94XNanoAODINPUT+NANOEDMMC2017+HARVESTNANOAODMC2017

Comparison Summary:

  • No significant changes to the logs found
  • Reco comparison results: 22451 differences found in the comparisons
  • DQMHistoTests: Total files compared: 27
  • DQMHistoTests: Total histograms compared: 2466491
  • DQMHistoTests: Total failures: 83949
  • DQMHistoTests: Total nulls: 4
  • DQMHistoTests: Total successes: 2382369
  • DQMHistoTests: Total skipped: 169
  • DQMHistoTests: Total Missing objects: 0
  • DQMHistoSizes: Histogram memory added: 0.909999999938 KiB( 22 files compared)
  • Checked 111 log files, 9 edm output root files, 27 DQM output files

@civanch
Copy link
Contributor Author

civanch commented Feb 8, 2018

+1
differences expected due to change of SIM histories

@cmsbuild
Copy link
Contributor

Comparison is ready
https://cmssdt.cern.ch/SDT/jenkins-artifacts/pull-request-integration/PR-22145/26234/summary.html

Comparison Summary:

  • No significant changes to the logs found
  • Reco comparison results: 26041 differences found in the comparisons
  • DQMHistoTests: Total files compared: 28
  • DQMHistoTests: Total histograms compared: 2498056
  • DQMHistoTests: Total failures: 105196
  • DQMHistoTests: Total nulls: 2
  • DQMHistoTests: Total successes: 2392682
  • DQMHistoTests: Total skipped: 176
  • DQMHistoTests: Total Missing objects: 0
  • DQMHistoSizes: Histogram memory added: 0.989999999969 KiB( 23 files compared)
  • Checked 115 log files, 9 edm output root files, 28 DQM output files

@civanch
Copy link
Contributor Author

civanch commented Feb 23, 2018

+1

@fabiocos
Copy link
Contributor

fabiocos commented Mar 1, 2018

@civanch could you please document (may be in the next simulation meeting) the differences that you have observed, so as we may quantify them? If they are really a few percent effect they do not look totally irrelevant to me.

@abdoulline
Copy link

Our (HCAL) 50 GeV single-pion Calo-scan test may not be very useful to spot fine details, as due to the change in SIM rndm sequence(s), the results at this level of stat (50k ev) fluctuate within several % easily...
But at least there is no major issues observed:
https://cms-cpt-software.web.cern.ch/cms-cpt-software/General/Validation/SVSuite/HCAL/calo_scan_single_pi/10_X/10_1_0_pre1_HFfix_vs_10_1_0_pre1_SinglePi/

@civanch
Copy link
Contributor Author

civanch commented Mar 2, 2018

10k MinBias events - a small suppression of the high energy tail:
http://test-civanch.web.cern.ch/test-civanch/run_matrix/res_mb/e_HF.png

@fabiocos
Copy link
Contributor

@civanch thanks, how were the histograms normalized? I would expect to see no difference, within statistics, for both EE and HE, right?

@civanch
Copy link
Contributor Author

civanch commented Mar 12, 2018

@fabiocos , statistically distributions are not in contrudiction I guess. There is an expresion of some systematics but I would not say it isstatistically significant. Also Hypérons and ions do not produce HF showers, this may slightly change hits in HE and EE. Normalisation of shown histograms is not done - the same number of events 50k is used for both cases.

@fabiocos
Copy link
Contributor

@civanch @franzoni @klute as discussed, this will enter into a test build on top of 10_2_0_pre2, and a few significant dedicated test samples with decent statistics should be agreed with JetMET to test its effect on jets in HF and MET

@fabiocos
Copy link
Contributor

unhold

@cmsbuild
Copy link
Contributor

This pull request is fully signed and it will be integrated in one of the next master IBs (tests are also fine). This pull request will now be reviewed by the release team before it's merged. @davidlange6, @slava77, @smuzaffar, @fabiocos (and backports should be raised in the release meeting by the corresponding L2)

@fabiocos
Copy link
Contributor

@civanch @franzoni @klute as agreed we will check this PR in 10_2_0_pre2 . As the expected effect of this PR is small, in order to make the validation work lighter we can "plan for success", as David likes, and integrate it into the pre-release, assuming that only in case unexpected effects are observed we make a dedicated build where this PR is removed for further investigations.

@fabiocos
Copy link
Contributor

+1

@cmsbuild cmsbuild merged commit caf9c33 into cms-sw:master Apr 22, 2018
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.

None yet

6 participants