Ensure GEANT simulation of exotic long-lived particles that decay outside beampipe

[wulsin]

To accomplish this, for any particle with status > 3, set its status to 2. All such particles are then checked to see if their decay vertex is outside the beampipe; if so it will be passed to GEANT for simulation.

This fix addresses the problem of a long-lived exotic particle with status code > 3 (e.g., an R-hadron, or chargino), which decays outside the beampipe. Before this fix, its detector interactions would not be simulated by GEANT, because this module was looking specifically for status values of 1, 2, or 3.

This fix has been tested for the case of di-chargino pair production, with the chargino ctau set to 1000 cm (decaying to neutralino + pion). Before the fix, there are no tracks associated to the chargino because no energy deposits have been simulated by GEANT. After the fix, there are reconstructed tracks for each chargino.

[cmsbuild]

A new Pull Request was created by @wulsin (Wells Wulsin) for CMSSW_8_1_X.

It involves the following packages:

SimG4Core/Generators

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

cms-bot commands are listed here #13028

[civanch]

please test

[cmsbuild]

The tests are being triggered in jenkins.
https://cmssdt.cern.ch/jenkins/job/ib-any-integration/15994/console

[cmsbuild]

-1

Tested at: 3f66f8a

You can see the results of the tests here:
https://cmssdt.cern.ch/SDT/jenkins-artifacts/pull-request-integration/PR-16366/15994/summary.html

I found follow errors while testing this PR

Failed tests: RelVals AddOn

• RelVals:

When I ran the RelVals I found an error in the following worklfows:
25202.0 step1

runTheMatrix-results/25202.0_TTbar_13+TTbar_13+DIGIUP15_PU25+RECOUP15_PU25+HARVESTUP15_PU25/step1_TTbar_13+TTbar_13+DIGIUP15_PU25+RECOUP15_PU25+HARVESTUP15_PU25.log

50202.0 step1

runTheMatrix-results/50202.0_TTbar_13+TTbar_13+DIGIUP15_PU50+RECOUP15_PU50+HARVESTUP15_PU50/step1_TTbar_13+TTbar_13+DIGIUP15_PU50+RECOUP15_PU50+HARVESTUP15_PU50.log

10024.0 step1

runTheMatrix-results/10024.0_TTbar_13+TTbar_13TeV_TuneCUETP8M1_2017_GenSimFull+DigiFull_2017+RecoFull_2017+ALCAFull_2017+HARVESTFull_2017/step1_TTbar_13+TTbar_13TeV_TuneCUETP8M1_2017_GenSimFull+DigiFull_2017+RecoFull_2017+ALCAFull_2017+HARVESTFull_2017.log

20024.0 step1

runTheMatrix-results/20024.0_TTbar_13+TTbar_13TeV_TuneCUETP8M1_2023D1_GenSimHLBeamSpotFull+DigiFull_2023D1+RecoFullGlobal_2023D1+HARVESTFullGlobal_2023D1/step1_TTbar_13+TTbar_13TeV_TuneCUETP8M1_2023D1_GenSimHLBeamSpotFull+DigiFull_2023D1+RecoFullGlobal_2023D1+HARVESTFullGlobal_2023D1.log

• AddOn:

I found errors in the following addon tests:

cmsDriver.py TTbar_Tauola_8TeV_cfi -s GEN,SIM,DIGI,L1,DIGI2RAW --mc --scenario=pp -n 10 --conditions auto:run1_mc_Fake --relval 9000,50 --datatier "GEN-SIM-RAW" --eventcontent RAWSIM --customise=HLTrigger/Configuration/CustomConfigs.L1T --fileout file:RelVal_Raw_Fake_MC.root : FAILED - time: date Thu Oct 27 02:38:43 2016-date Thu Oct 27 02:34:50 2016 s - exit: 34304
cmsRun /cvmfs/cms-ib.cern.ch/week0/slc6_amd64_gcc530/cms/cmssw/CMSSW_8_1_X_2016-10-26-1100/src/HLTrigger/Configuration/test/OnLine_HLT_Fake.py realData=False globalTag=@ inputFiles=@ : FAILED - time: date Thu Oct 27 02:38:43 2016-date Thu Oct 27 02:34:50 2016 s - exit: 21504
cmsDriver.py RelVal -s HLT:Fake,RAW2DIGI,L1Reco,RECO --mc --scenario=pp -n 10 --conditions auto:run1_mc_Fake --relval 9000,50 --datatier "RAW-HLT-RECO" --eventcontent FEVTDEBUGHLT --customise=HLTrigger/Configuration/CustomConfigs.L1THLT --processName=HLTRECO --filein file:RelVal_Raw_Fake_MC.root --fileout file:RelVal_Raw_Fake_MC_HLT_RECO.root : FAILED - time: date Thu Oct 27 02:38:43 2016-date Thu Oct 27 02:34:50 2016 s - exit: 21504

[cmsbuild]

Comparison not run due to runTheMatrix errors (RelVals and Igprof tests were also skipped)

[civanch]

@wulsin , my interpretation of the problem: some particles produced in TTbar_Tauola cannot be tracked by Geant4, so cannot have status=2.

[wulsin]

@civanch , Do you have a suggestion for how to avoid this problem? For certain particles, we could leave their status unchanged. Or perhaps there is another solution?

[cmsbuild]

Pull request #16366 was updated. @cmsbuild, @civanch, @mdhildreth, @davidlange6 can you please check and sign again.

[wulsin]

I changed this so that only exotic particles have their status modified. This is to avoid propagating certain Standard Model particles (like K0), that GEANT does not know how to handle. Let's see if this change solves the errors that we saw before.

[civanch]

please test

[cmsbuild]

The tests are being triggered in jenkins.
https://cmssdt.cern.ch/jenkins/job/ib-any-integration/16057/console

[cmsbuild]

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

[civanch]

+1

[cmsbuild]

This pull request is fully signed and it will be integrated in one of the next CMSSW_8_1_X IBs (tests are also fine). This pull request requires discussion in the ORP meeting before it's merged. @slava77, @davidlange6, @smuzaffar

[davidlange6]

assign generators
(as I'd like to be sure the exit codes of pythia8 are correctly interpreted)

[cmsbuild]

New categories assigned: generators

@thuer,@perrozzi,@govoni you have been requested to review this Pull request/Issue and eventually sign? Thanks

[perrozzi]

The PR looks good to me from the GEN point of view. However let me tag also @bendavid to crosscheck as I never specifically used this module. Josh can you please crosscheck?

[civanch]

@thuer,@perrozzi,@govoni, @bendavid, please, sign or comment.

[civanch]

@davidlange6 , it seems that GEN has "soft" approval and this PR cannot wait more.

[perrozzi]

+1

[cmsbuild]

This pull request is fully signed and it will be integrated in one of the next CMSSW_8_1_X IBs (tests are also fine). This pull request requires discussion in the ORP meeting before it's merged. @slava77, @davidlange6, @smuzaffar

[davidlange6]

+1