-
Notifications
You must be signed in to change notification settings - Fork 10
/
Copy pathweights.txt
50 lines (37 loc) · 2.18 KB
/
weights.txt
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
#sf_ewkV
#
#process kz
sf_ewk = valueAt(,h1Corrs[cZEWK] = new THCorr1((TH1D*)fKFactor->Get("EWKcorr/Z"));,event.genParticles.bound(gen.pt(),genBosonPtMin,genBosonPtMax);)
h1Corrs[cZEWK]->GetHist()->Divide(h1Corrs[cZNLO]->GetHist());
#process kw
sf_ewkV = valueAt(corrDir,h1Corrs[cWEWK] = new THCorr1((TH1D*)fKFactor->Get("EWKcorr/W"));,events.genParticles.bound(gen.pt(),genBosonPtMin,genBosonPtMax);)
h1Corrs[cWEWK]->GetHist()->Divide(h1Corrs[cWNLO]->GetHist());
#process ka
sf_ewkV = valueAt(corrDict,h1Corrs[cAEWK] = new THCorr1((TH1D*)fKFactor->Get("EWKcorr/photon"));,event.genparticles.bound(gen.pt(),genBosonPtMin,genBosonPtMax);)
h1Corrs[cAEWK]->GetHist()->Divide(h1Corrs[cANLO]->GetHist());
#(when generated boson was not found)
#process kz same as before but gt->genBosonPt = bound(vpt.Pt(),genBosonPtMin,genBosonPtMax);
#process kw same as before but gt->genBosonPt = bound(vpt.Pt(),genBosonPtMin,genBosonPtMax);
#sf_qcdV
gt->genBosonPt = gt->genBosonPt = bound(gen.pt(),genBosonPtMin,genBosonPtMax);
#process kz
sf_qcdV = valueAt(corrDict,h1Corrs[cWNLO] = new THCorr1((TH1D*)fKFactor->Get("WJets_012j_NLO/nominal")); ,gt->genBosonPt)
#sf_metTrig
sf_metTrig = valueAt(corrDict,look below, module of the sum of muon momentum vectors)
OpenCorrection(cTrigMET,dirPath+"moriond17/metTriggerEfficiency_recoil_monojet_TH1F.root",
"hden_monojet_recoil_clone_passed",1);
#sf_pho
sf_pho = valueAT(corrDict,look below,event.photons eta and pt)
OpenCorrection(cPho,dirPath+"moriond17/scalefactors_80x_medium_photon_37ifb.root",
"EGamma_SF2D",2);
#sf_phoTrig
if simplephoton -> if number of loose photons = 1, pt = that photon pt
if complicated photon -> pt is the momentum of photon with greates momentum
sf_phoTrig = valueAt(corrDict, below, event.photon.pt)
OpenCorrection(cTrigPho,dirPath+"moriond17/photonTriggerEfficiency_photon_TH1F.root",
"hden_photonpt_clone_passed",1);
#sf_qcdV2j
gen = event.genParticles particle
gt->genBosonPt = bound(gen.pt(),genBosonPtMin,genBosonPtMax);
sf_qcdV2j = valueAt(corrDict, below, gt->genBosonPt)
OpenCorrection(cANLO2j,dirPath+"moriond17/histo_photons_2jet.root","Func",1);