TWiki> Main Web>TWikiUsers>MarkElliotSmith>B2PpbarMunu>B2PPbarTauNu>ImperialKstMuMu (2021-04-09, MatthewBirch) EditAttachPDF

Repository

Data location:

  • 2011 - /LHCb/Collision11/Beam3500GeV-VeloClosed-Mag[Up,Down]/Real Data/Reco14/Stripping21r1p1a/90000000/LEPTONIC.MDST
  • 2012 - /LHCb/Collision12/Beam4000GeV-VeloClosed-Mag[Up,Down]/Real Data/Reco14/Stripping21r0p1a/90000000/LEPTONIC.MDST
  • 2016 - /LHCb/Collision16/Beam6500GeV-VeloClosed-Mag[Up,Down]/Real Data/Reco16/Stripping28r1/90000000/LEPTONIC.MDST
  • 2017 - /LHCb/Collision17/Beam6500GeV-VeloClosed-Mag[Up,Down]/Real Data/Reco17/Stripping29r2/90000000/LEPTONIC.MDST
  • 2018 - /LHCb/Collision18/Beam6500GeV-VeloClosed-Mag[Up,Down]/Real Data/Reco18/Stripping34/90000000/LEPTONIC.MDST
The stripping lines of interest are B2XMuMu_Line and B2LLXBDT_Bd2mumuKstarLine .

MC location

  • 2011 $J/\psi K^*$ - /MC/2011/Beam3500GeV-2011-MagDown-Nu2-Pythia8/Sim09c/Trig0x40760037/Reco14c/Stripping21r1p1Filtered/11144001/BU2LLK_NOPID.STRIP.DST - 578333 + 598755
  • 2011 $K^* \mu\mu$ - /MC/2011/Beam3500GeV-2011-MagDown-Nu2-Pythia8/Sim08e/Digi13/Trig0x40760037/Reco14a/Stripping20r1NoPrescalingFlagged/11114001/ALLSTREAMS.DST - 525851 + 506402 - don't use this, it is too small!
  • 2012 $J/\psi K^*$ - /MC/2012/Beam4000GeV-2012-MagDown-Nu2.5-Pythia8/Sim09c/Trig0x409f0045/Reco14c/Stripping21r0p1Filtered/11144001/BU2LLK_NOPID.STRIP.DST - 1279726 + 1208524
  • 2012 $K^* \mu\mu$ PHSP - /MC/2012/Beam4000GeV-2012-MagDown-Nu2.5-Pythia8/Sim08i/Digi13/Trig0x409f0045/Reco14c/Stripping21r0p1NoPrescalingFlagged/11114003/ALLSTREAMS.DST - 1028324 + 1008308 - too small
  • 2012 $K^* \mu\mu$ PHSP - /MC/2012/Beam4000GeV-2012-MagDown-Nu2.5-Pythia6/Sim08b/Digi13/Trig0x409f0045/Reco14a/Stripping20Filtered/11114005/B2XMUMU.STRIP.DST - 800589 + 781112
  • 2012 $K^* \mu\mu$ PHSP - /MC/2012/Beam4000GeV-2012-MagDown-Nu2.5-Pythia6/Sim08e/Digi13/Trig0x409f0045/Reco14a/Stripping20Filtered/11114005/B2XMUMU.STRIP.DST - 617964 + 674851
  • 2016 $J/\psi K^*$ - /MC/2016/Beam6500GeV-2016-Mag[Up,Down]-Nu1.6-25ns-Pythia8/Sim09b/Trig0x6138160F/Reco16/Turbo03/Stripping26NoPrescalingFlagged/11144001/ALLSTREAMS.DST - 1252070 + 1197537
  • 2017 $J/\psi K^*$ - /MC/2017/Beam6500GeV-2017-MagDown-Nu1.6-25ns-Pythia8/Sim09f/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2NoPrescalingFlagged/11144001/ALLSTREAMS.DST - 2032752 + 2133675
  • 2018 $J/\psi K^*$ - /MC/2018/Beam6500GeV-2018-MagDown-Nu1.6-25ns-Pythia8/Sim09f/Trig0x617d18a4/Reco18/Turbo05-WithTurcal/Stripping34NoPrescalingFlagged/11144001/ALLSTREAMS.DST - 2008602 + 2001400
  • 2016 $K^{*}\mu\mu$ - /MC/2016/Beam6500GeV-2016-MagDown-Nu1.6-25ns-Pythia8/Sim09b/Trig0x6138160F/Reco16/Turbo03/Stripping26NoPrescalingFlagged/11114001/ALLSTREAMS.DST - 657410 + 707852
  • 2016 $K^{*}\mu\mu$ b2sllball - /MC/2016/Beam6500GeV-2016-MagDown-Nu1.6-25ns-Pythia8/Sim09c/Trig0x6138160F/Reco16/Turbo03/Stripping28r1p1Filtered/11114002/BU2LLK_NOPID.STRIP.DST - 1132264 + 1156791
  • 2016 $K^{*}\mu\mu$ flat $q^{2}$ - /MC/2016/Beam6500GeV-2016-MagDown-Nu1.6-25ns-Pythia8/Sim09b/Trig0x6138160F/Reco16/Turbo03/ReconstructibleFiltered/11114014/MCTRUE.DST - 21157238 + 21039842
  • 2016 $K^{*}\mu\mu$ flat $q^{2}$ - /MC/2016/Beam6500GeV-2016-MagDown-Nu1.6-25ns-Pythia8/Sim09c/Trig0x6138160F/Reco16/Turbo03/ReconstructibleFiltered/11114014/MCTRUE.LDST - 4263449 + 4264509
  • 2017 $K^{*}\mu\mu$ flat $q^{2}$ - /MC/2017/Beam6500GeV-2017-MagDown-Nu1.6-25ns-Pythia8/Sim09f/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2Filtered/Merge14/11114014/MCTRUE.DST - 5135615 + 5015779

Workflow for MC

  • Make your tuples.
  • Preselection in the Selection folder. This applies minimal things but it will collate your tuples. It also calculates the angles and $q^{2}$. You should be able to start from the tuples produced here subsequently.
  • Apply PID resampling in the PID PID_sampling folder.
  • Apply the L0 emulation, TISTOS weights and tracking efficiency correction weights with the apply_weights.py script in the track_eff folder.
  • Apply final selection in the Final_Selection folder. This adds the BDT and does most of the selection, including vetoes and (corrected) PID.
  • If you are doing the $J/\psi K^{*}$ MC, then you also need to add the S-wave weights with the script in the S-wave folder.
  • Apply kinematic reweighting with the python script in kinematic_reweighter . This will apply the L0 efficiency weights and the tracking correction.

Acceptance

  • All the code is in the Acceptance folder. Should be fairly self explanatory.
Optimum orders:
Sample $q^2$ range O($\cos{\theta_l}$) O($\cos{\theta_K}$) O($\phi$) O($q^2$) O($m_{K\pi}$)
2012 0.098 - 1.2 4 4 0 2 1
2012 0.9 - 19.1 4 6 6 3 3
2016 0.098 - 1.2 4 4 0 2 2
2016 0.9 - 19.1 4 4 6 4 2
2017/2018 0.098 - 1.2 4 4 0 2 1
2017/2018 0.9 - 19.1 4 6 6 4 2

Workflow for data

  • Make your tuples.
  • Preselection in the Selection folder. This adds the BDT and applies most of the cuts. NOTE: It does not apply a cut on the BDT or L0!
  • Use perform_mass_cuts macro to apply the necessary mass cuts to make the rare-mode / $J/\psi$ samples.

MC tags and DaVinci versions:

Sample Event No DDDB CondDB
2011 $K^* \mu\mu$ for acceptance 11114001 dddb-20130929 sim-20130522-vc-md100
2011 $J/\psi K^*$ filtered, noPID 11144001 dddb-20170721-1 sim-20160614-1-vc-md100
2012 $K^* \mu\mu$ PHSP - sim08b for acceptance 11114005 Sim08-20130503-1 Sim08-20130503-1-vc-md100
2012 $K^* \mu\mu$ PHSP - sim08e for acceptance 11114005 dddb-20130929-1 sim-20130522-1-vc-md100
2012 $J/\psi K^*$ filtered, noPID 11144001 dddb-20170721-2 sim-20160321-2-vc-md100
2016 $J/\psi K^*$ 11144001 dddb-20150724 sim-20161124-2-vc-mu100
2017 $J/\psi K^*$ 11144001 dddb-20170721-3 sim-20180411-vc-md100
2018 $J/\psi K^*$ 11144001 dddb-20170721-3 sim-20190128-vc-md100
2016 $K^{*}\mu\mu$ 11114001 dddb-20150724 sim-20161124-2-vc-md100
2016 $K^{*}\mu\mu$ b2sllball 11114002 dddb-20170721-3 sim-20170721-2-vc-md100
2016 $K^{*}\mu\mu$ flat $q^{2}$ sim09b 11114014 dddb-20150724 sim-20161124-2-vc-md100
2016 $K^{*}\mu\mu$ flat $q^{2}$ sim09c 11114014 dddb-20170721-3 sim-20170721-2-vc-md100
2017 $K^{*}\mu\mu$ flat $q^{2}$ sim09f 11114014 dddb-20170721-3 sim-20180411-vc-md100
       

Year Stripping DaVinci
2011 v21r1p1a v39r1p1
2012 v21r0p1a v39r1p1
2016 28r1 v41r4p4
2017 29r2 v42r7p2
2018 34 v44r4

MC productions etc:

Sample Prod No Gauss L0 TCK L0 Muon L0 DiMuon
2017 $K^{*}\mu\mu$ flat $q^{2}$ sim09f 86034, 86035 Gauss v49r12 0x1709 'SumEtPrev<1000', 'Spd_Had(Mult)<450', 'Muon1(Pt)>28' 'Spd_DiMu(Mult)<900', 'Muon12(Pt)>676'
2016 $K^{*}\mu\mu$ flat $q^{2}$ sim09c 66612   0x160F 'Spd_Had(Mult)<450', 'Muon1(Pt)>36' 'Spd_DiMu(Mult)<900', 'Muon12(Pt)>900'
2016 $K^{*}\mu\mu$ flat $q^{2}$ sim09b 57822   0x160F 'Spd_Had(Mult)<450', 'Muon1(Pt)>36' 'Spd_DiMu(Mult)<900', 'Muon12(Pt)>900'

L0 Thresholds

Year TCK % of data L0 Muon L0 DiMuon
2016   2.05 23 400
    1.5 27 576
    4.6 31 676
    47.2 27 676
    2.3 31 676
    34 37 900
    2.7 31 784
    5.1 33 900
2017   0.7 15 324
    1.3 23 400
    1.5 27 576
    7.3 31 676
    1.6 39 1296
    38.7 35 1296
    14.6 23 400
    34.2 29 676
2018 0x1803 100 [SumEtPrev<1000] && [Spd Had(Mult)<450] && [Muon1(Pt)>35] [Spd DiMu (Mult)<900] && [Muon12(Pt)>1296]
I only list one 2018 TCK as the thresholds didn't change throughout the year.

nTracks scaling

2012  
2016 1.19
2017  
2018  

Amazing plots

Useful things

Notes

  • The best fit value to the published LHCb analysis, using all variables, is: $C_{9}=-1.027$, $C_{10}=0.498$
  • The best fit value to the published LHCb analysis using only the P variables is: $C_{9}=-0.970$, $C_{10}=0.324$
-- MarkElliotSmith - 2019-05-09
Attachments Attachments
Topic attachments
I Attachment History Action Size Date Who Comment
PDFpdf 2015_2017_data_17_09_2018.pdf r1 manage 3595.5 K 2019-07-19 - 15:06 MarkElliotSmith List of TCKs for the various years
PDFpdf Kstmumu_note_v0.pdf r1 manage 3647.9 K 2020-11-25 - 15:30 MarkElliotSmith note version 0
PDFpdf observables.pdf r1 manage 970.9 K 2020-06-29 - 23:21 MichaelMcCann  
Edit | Attach | Watch | Print version | History: r33 < r32 < r31 < r30 < r29 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r33 - 2021-04-09 - MatthewBirch
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright &© 2008-2021 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback