My Links
miniAOD -> ntuple
1) Robie tak jak do ntupli synchronizacyjnych. Biore próbki:
#recoilcorrections
#'/WJetsToLNu_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISummer16MiniAODv2-PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6-v1/MINIAODSIM',
#'/W1JetsToLNu_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISummer16MiniAODv2-PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6-v1/MINIAODSIM',
#'/W2JetsToLNu_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISummer16MiniAODv2-PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6-v1/MINIAODSIM',
#'/W3JetsToLNu_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISummer16MiniAODv2-PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6-v1/MINIAODSIM',
#'/W4JetsToLNu_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISummer16MiniAODv2-PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6-v1/MINIAODSIM',
#'/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISummer16MiniAODv2-PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6_ext1-v2/MINIAODSIM',
'/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISummer16MiniAODv2-PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6_ext2-v1/MINIAODSIM',
#'/DY1JetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISummer16MiniAODv2-PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6-v1/MINIAODSIM',
#'/DYJetsToLL_M-10to50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/RunIISummer16MiniAODv2-PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6-v1/MINIAODSIM',
#no recoilcorrections
#'/TT_TuneCUETP8M2T4_13TeV-powheg-pythia8/RunIISummer16MiniAODv2-PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6-v1/MINIAODSIM',
#'/ZZTo4L_13TeV_powheg_pythia8/RunIISummer16MiniAODv2-PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6-v1/MINIAODSIM',
#'/VVTo2L2Nu_13TeV_amcatnloFXFX_madspin_pythia8/RunIISummer16MiniAODv2-PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6-v1/MINIAODSIM',
#Data
#"/SingleMuon/Run2016B-23Sep2016-v3/MINIAOD",
#"/SingleMuon/Run2016C-23Sep2016-v1/MINIAOD",
#"/SingleMuon/Run2016D-23Sep2016-v1/MINIAOD",
#"/SingleMuon/Run2016E-23Sep2016-v1/MINIAOD",
#"/SingleMuon/Run2016F-23Sep2016-v1/MINIAOD",
#"/SingleMuon/Run2016G-23Sep2016-v1/MINIAOD",
2)
GlobalTagi: dla MC - 80X_mcRun2_asymptotic_2016_TrancheIV_v8; dla
RunB - 80X_dataRun2_2016SeptRepro_v4; dla Run C,D,E,F,G - 80X_dataRun2_2016SeptRepro_v3
3) Robię cięcia:
muon.pt() > 20 && fabs(muon.eta()) < 2.1
fabs(muon.muonBestTrack()->dxy(vtx.position())) < 0.045 &&
fabs(muon.muonBestTrack()->dz(vtx.position())) < 0.2
ID:
POG medium muon ID
tau.pt() > 20 and fabs(tau.eta()) < 2.3
ID:
TauIDRecommendation13TeV
#Then apply:
tau.tauID('decayModeFinding') > 0.5 and
fabs(packedLeadTauCand->dz()) < 0.2 # The PackedCandidate::dz() method is wrt. the first PV by default
Tau charge is equal to plus or minus 1.
DR(muon, tau) > 0.5
4) Wybieram najlepszą parę i robie ntuple;
ntuple -> rysunki
Robie następujące kroki:
- licze pileup. Dla
RunB generuje histogram pileupCalc.py. Jeśli robie wszystkie dane to biore Data_Pileup_2016_271036-284044_80bins.root. Dla MC: MC_Moriond17_PU25ns_V1.root.
- Wagi dla MC (mnoże wszystkie):
1. Pileup
2. Jeśli tau.isGenHadronicTau(): 0.95
3. Jeśli tau..isGenPromptMuon() or tau..isGenMuonTau(): to sf = [1.4, 1.72, 1.26, 2.6, 2.3] dla eta [0.4, 0.8, 1.2, 1.7, 2.3] kolejno
3. Jeśli WJety lub DY M > 50 to waga z tabeliki dla stitchingu
4. Jeśli pozdostałe MC to : xsec*lumi/sum(weight)
- Cięcia:
pair.diq() == -1 and mu.mt() < 50 and pair.PostSynchSelection() and pair.veto() and (pair.trigger('HLT_IsoMu24') or pair.trigger('HLT_IsoTkMu24')) \
and mu.iso() < 0.1 and mu.pt() > 26 and tau.pt() > 30 \
-. Histogramy MC stakuje
-. Histogramy qcd robie rządając pair.diq == 1. Dodaje oddzielnie dla Data i MC i potem robie Data-M. Dodaje do histogramu z pk.8
-. Rysuje Data i MC.
Related Topics