MC for displaced low mass dimuons (scouting data) UCSB : N. Amin, C. - - PowerPoint PPT Presentation

mc for displaced low mass dimuons scouting data
SMART_READER_LITE
LIVE PREVIEW

MC for displaced low mass dimuons (scouting data) UCSB : N. Amin, C. - - PowerPoint PPT Presentation

MC for displaced low mass dimuons (scouting data) UCSB : N. Amin, C. Campagnari, U. Sarica UCSD : V. Krutelyov, M. Masciovecchio, A. Yagil Boston : I. Suarez, A. Tsatsos Rutgers : A. Gandrakota, A. Lath, H. Routray Nebraska : F . Golf July 6, 2020


slide-1
SLIDE 1

MC for displaced low mass dimuons (scouting data)

UCSB: N. Amin, C. Campagnari, U. Sarica UCSD: V. Krutelyov, M. Masciovecchio, A. Yagil Boston: I. Suarez, A. Tsatsos Rutgers: A. Gandrakota, A. Lath, H. Routray Nebraska: F . Golf July 6, 2020

slide-2
SLIDE 2

⚫ Look for events with at least 2 displaced muons and an associated

displaced vertex (DV)

  • Datasets /ScoutingCaloMuon/Run201{7,8}*/RAW
  • HLT: DST_DoubleMu3_noVtx_CaloScouting (muon pT > 3GeV, |𝜃| < 2.4)
  • Consider full mass range (down to 2mmuon) and displacements

between 0 and 11cm

⚫ In addition to model-independent results, provide interpretations for 3

main signal models

Introduction to analysis

2

  • Example

t with EXO MC):

erpretation.

H

Zd Zd

gg → φ

H → ZdZd

B → φ

slide-3
SLIDE 3

⚫ Select events with 2 good ScoutingCaloMuons and

associated displaced ScoutingVertex with selections below

Analysis selections

3

Baseline selections

At least 2 OS muons and 1 DV in collections, then

  • DV
  • (x, y, z) errors < (0.05, 0.05, 0.1) cm
  • chi2/ndof < 5
  • 𝜍 < 11 cm
  • Muon (pT > 3 GeV, |𝜃| < 2.4)
  • ID
  • Num. tracker layers with meas. > 5
  • chi2/ndof < 3
  • Isolation
  • Track isolation < 0.1
  • 𝛦R with closest jet > 0.3
  • Dimuon kinematics
  • |𝛦𝜚(muon 1, muon 2)| < 2.8
  • |𝛦𝜚(dimuon, DV vector)| < 0.02

2017 L1 trigger seeds

OR of

  • DoubleMu4_SQ_OS_dR_Max1p2
  • DoubleMu0er1p4_SQ_OS_dR_Max1p4
  • DoubleMu_15_7

2018 L1 trigger seeds

OR of

  • DoubleMu4p5_SQ_OS_dR_Max1p2
  • DoubleMu0er1p4_SQ_OS_dR_Max1p4
  • DoubleMu_15_7

Additional selections

⚫ distance(DV, nearest pixel module plane) > 0.05cm ⚫ if DV 𝜍 > 3.5cm, expected pixel hits ≤ observed pixel hits, for both muons ⚫ log10(|𝛦𝜃𝜈𝜈|/|𝛦𝜚𝜈𝜈|) < 1.25 ⚫ |dxy/𝜏dxy| > 2 for both muons ⚫ |dxy|/(𝜍m/pT,𝜈𝜈) > 0.1 for both muons
slide-4
SLIDE 4

⚫ After making baseline selections, bin events by

  • DV 𝜍: [0, 0.2, 1, 2.4, 3.1, 7, 11] cm
  • dimuon pT: [0, 25, ∞] GeV
  • both muons isolated/one muon non-isolated

⚫ Fit the background distribution to evaluate exclusion limits ⚫ More analysis details in a previous LL EXO presentation

Signal extraction

4

𝜍>1 cm 𝜍<1 cm isolated category

slide-5
SLIDE 5

⚫ gg → 𝜚

  • MadGraph gridpack generating 𝜚 with one extra parton
  • Subsequent Pythia fragment gives 𝜚 lifetime and forces 𝜈𝜈 decay

⚫ H → ZdZd

  • One Zd decays to 𝜈𝜈 and the other according to SM Z BRs
  • Using latest version of Higgs Powheg gridpacks (/cvmfs/cms.cern.ch/phys_generator/

gridpacks/2017/13TeV/powheg/V2/gg_H_quark-mass-effects_ZZ_NNPDF31_13TeV/gg_H_quark-mass- effects_NNPDF31_13TeV_M125/v4/gg_H_quark-mass-effects_NNPDF31_13TeV_M125.tgz) with modified

JHUGen decay files, making SM Z have LL properties, for a given mass and lifetime

  • Default pythia fragment used afterward

⚫ B → 𝜚

  • Generation starting with Pythia-only, for a given mass and lifetime
  • Gives rise to soft muons, so this is accompanied by a gen-level filtering module

requiring two pT>3GeV, |𝜃|<2.4 muons

  • Filter efficiency ~1% for m𝜚=2GeV, cτ=20mm point

Generation details

5

slide-6
SLIDE 6

Generated kinematics

6

slide-7
SLIDE 7

⚫ To get an idea of selection efficiencies wrt generation, consider H→ZdZd

(mZd=8GeV, cτ=10mm) events after requiring exactly 2 muons with pT>3GeV, |𝜃| <2.4 at gen. level, then plot efficiencies vs lxy or trailing muon pT

⚫ Due to HLT configuration requiring at least 2 pixel hits, 0 efficiency past lxy~11cm

Selection efficiencies

7

  • Gen. lxy (cm)

Trailing gen. muon pT (GeV)

slide-8
SLIDE 8

⚫ We would like to request a 2D grid (mass,cτ) of benchmark samples with 2017 and 2018 conditions

  • 100k events per point
  • For each model, this is 3x3x2=18 datasets, 1.8M events

⚫ Scouting collections only stored up to AODSIM, so in principle, we do not need any tiers beyond that ⚫ Links to gridpacks and fragments in the summary table below (subdirectories here) ⚫ If appropriate, will fill out request spreadsheet entries after this meeting ⚫ NOTE: We are still investigating some details with the B→φ model, but wish to proceed with the

generation of the other two

Production requests

8

model gridpack Pythia fragment mass points (GeV) cτ points (mm) conditions gg → φ MG gridpack (set mass) modified (set cτ) 1, 8, 15 1, 10, 100 2017, 2018 H → ZdZd powheg + JHUGen gridpack (set mass, cτ) default 1, 8, 15 1, 10, 100 2017, 2018 B → φ no gridpack modified + gen. filter (set mass, cτ) 1, 2, 4 1, 10, 100 2017, 2018

slide-9
SLIDE 9

Backup

9

slide-10
SLIDE 10

⚫ Wei Shi previously reported issues with LLP not being propagated past beam pipe

  • They use a patch to CMSSW to specially treat a particular pdgId

⚫ We observed the same issue

  • Worked around it by using pdgId>4000000 (e.g., 6000211) for the 𝜚 particle (no need to edit CMSSW)
  • For the Zd particle, which is a modified SM Z (pdgId 23), propagation works well out of the box

Displaced simulation issues

10 from Wei Shi

slide-11
SLIDE 11

⚫ Make matrix of reweighting

between c𝜐={1,10,50,100}mm samples

  • Plotting reco lxy

distributions after analysis selections

⚫ Lower diagonal (reweighting

high c𝜐 to low c𝜐) shows good agreement (compare orange and green) and doesn't suffer from extreme weights

Lifetime reweighting

11

w(t) = cτ0 cτ1 exp ( ct cτ0 − ct cτ1)

ct = c( ⃗ L xy ⋅ ⃗ pT) m p2

T

slide-12
SLIDE 12

Cutflow efficiencies

12

⚫ pass L1 OR: OR of 3 seeds ⚫ pass HLT: DST_DoubleMu3_noVtx_CaloScouting ⚫ pass object sel: event has at least 1 DV and 2 matching muons,

with some baseline cuts on the DV

⚫ pass iso. baseline: isolated baseline selection ⚫ pass extra cuts:

  • distance(DV, nearest pixel module plane) > 0.05cm
  • if DV 𝜍 > 3.5cm, expected pixel hits ≤ observed pixel hits
  • log10(|𝛦𝜃𝜈𝜈|/|𝛦𝜚𝜈𝜈|) < 1.25

⚫ pass dxy cuts: (|dxy/𝜏dxy| > 2) AND (|dxy|/(𝜍m/pT,𝜈𝜈) > 0.1) for both

muons

⚫ Cuts are cumulative starting with the denominator

(using gen. level quantities):

  • Exactly 1 Zd with 2 muons (pT > 5 GeV, |𝜃| < 2.4)
  • 𝜍 < 11cm
  • 𝛦R𝜈𝜈<1.2
  • distance(DV, nearest pixel module plane) > 0.05cm
  • log10(|𝛦𝜃𝜈𝜈|/|𝛦𝜚𝜈𝜈|) < 1.25

H → ZdZd