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
Added skimming filter for disappearing muons #40666
Conversation
-code-checks Logs: https://cmssdt.cern.ch/SDT/code-checks/cms-sw-PR-40666/34006
Code check has found code style and quality issues which could be resolved by applying following patch(s)
|
-code-checks Logs: https://cmssdt.cern.ch/SDT/code-checks/cms-sw-PR-40666/34007
Code check has found code style and quality issues which could be resolved by applying following patch(s)
|
Applied code-checks Added code-checks and code-format
adb4853
to
d005f43
Compare
+code-checks Logs: https://cmssdt.cern.ch/SDT/code-checks/cms-sw-PR-40666/34008
|
A new Pull Request was created by @revering for master. It involves the following packages:
@cmsbuild, @bbilin, @sunilUIET, @kskovpen can you please review it and eventually sign? Thanks. cms-bot commands are listed here |
please test |
-1 Failed Tests: ClangBuild Clang BuildI found compilation warning while trying to compile with clang. Command used:
See details on the summary page. |
Added check of eta and phi when matching the track to a vertex.
+code-checks Logs: https://cmssdt.cern.ch/SDT/code-checks/cms-sw-PR-40666/34595 |
Pull request #40666 was updated. @cmsbuild, @bbilin, @sunilUIET, @kskovpen can you please check and sign again. |
please test with 10805.31 |
Hi @perrotta , could you please test again? The bot seems not to receive this command. |
please test workflow 10805.31 |
Thank you @youyingli : it was my fault, indeed, as I wrote a wrong syntax for the command |
+1 Summary: https://cmssdt.cern.ch/SDT/jenkins-artifacts/pull-request-integration/PR-5397a2/31289/summary.html Comparison SummarySummary:
|
+pdmv |
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. @perrotta, @dpiparo, @rappoccio (and backports should be raised in the release meeting by the corresponding L2) |
+1 |
PR description:
Created a skim filter used for a search for disappearing muons (https://icms.cern.ch/tools-api/restplus/relay/piggyback/notes/AN/2020/214/files/3/download). The filter uses a tag-and-probe method to select high-quality muons which pair with probe tracks to form an invariant mass near the Z-peak, then extrapolates the probe track to the muon chambers and selects events with missing or low energy standalone muons nearby. Options are included to also select for control regions with same-sign tag/probe pairs or invariant masses away from the Z mass.
In the full analysis we need hit level data from the HCAL to fully establish our signal region, so this skim selection is used to choose a subset of interesting RECO events to keep the complete HCAL hit information in. This selection is designed to apply to AOD level files so that in existing runs we can choose events to reconstruct from RAW data without needing to re-reco the entire dataset. With the default options, roughly 0.5% of events in the single muon 2018 UL dataset are selected.
For an additional check of the expected event rate and data size, we ran ran the filter on a run 3 relVal RECO MC dataset /RelValZMM_14/CMSSW_13_0_0_pre3-PU_130X_mcRun3_2022_realistic_v2-v1/GEN-SIM-RECO. This DY->mumu process should be the most significant source of events passing our filter.
After keeping the full hbhe and ho hit collections, the relVal file had 2574 events in 1.7 GB and the post-filter file had 77 events in 60 MB, corresponding to a pass rate of 3% and an expected dataset size of 10 TB at 300 fb^-1. As an additional cross-check, we ran the skim on run 2 AOD data and confirmed a 0.5% pass rate, likely due to the inclusion of non-DY processes in the overall singlemuon dataset with lower acceptance rates than Z->mumu.
PR validation:
I ran the skim on a set of DY MC files and plotted each cut to confirm that the selections are correctly applied. I also ran our signal analyzer on the skimmed and non-skimmed MC to verify that there were no events cut by the filter which would have otherwise passed our signal or control region selections.