From b432d1258860af217e12fdcc77b27d84da9d5ce6 Mon Sep 17 00:00:00 2001 From: Andrey Pozdnyakov Date: Thu, 23 Aug 2018 16:30:53 +0200 Subject: [PATCH 1/2] Adding missing weight files --- Weights/MVAID/CSEVroot.py | 121 ++++++++++++ Weights/MVAID/MVAIDroot.py | 121 ++++++++++++ .../MVAID/ScalingFactors_80X_Summer16.root | Bin 0 -> 5612 bytes Weights/MVAID/csevsf.root | Bin 0 -> 3899 bytes Weights/MVAID/egammaEffi.txt_EGM2D.root | Bin 0 -> 9622 bytes Weights/MVAID/mvaIDsf.root | Bin 0 -> 4114 bytes Weights/TriggerSF/MakeTSFHistograms.py | 183 ++++++++++++++++++ Weights/TriggerSF/TriggerSFs.root | Bin 0 -> 6587 bytes 8 files changed, 425 insertions(+) create mode 100644 Weights/MVAID/CSEVroot.py create mode 100644 Weights/MVAID/MVAIDroot.py create mode 100644 Weights/MVAID/ScalingFactors_80X_Summer16.root create mode 100644 Weights/MVAID/csevsf.root create mode 100644 Weights/MVAID/egammaEffi.txt_EGM2D.root create mode 100644 Weights/MVAID/mvaIDsf.root create mode 100644 Weights/TriggerSF/MakeTSFHistograms.py create mode 100644 Weights/TriggerSF/TriggerSFs.root diff --git a/Weights/MVAID/CSEVroot.py b/Weights/MVAID/CSEVroot.py new file mode 100644 index 0000000..d83babf --- /dev/null +++ b/Weights/MVAID/CSEVroot.py @@ -0,0 +1,121 @@ +from ROOT import * +from array import array + +etabins = [-2.5, -1.566, -1.4442, -1.0, 0.0, 1.0, 1.4442, 1.566, 2.5] +ptbins = [20, 30, 40, 50, 200] + +nbins_eta = len(etabins)-1 +nbins_pt = len(ptbins)-1 + +sfs = [[ [0,0] for x in range(nbins_pt)] for y in range(nbins_eta)] + +##20-30, -2.5, -1.56 +sfs[0][0][0] = 1#0.989453 +sfs[0][0][1] = 1#0.0400067 +##20-30, -1.56, -1.44 +sfs[1][0][0] = 0 +sfs[1][0][1] = 0 +##20-30, -1.44, -1 +sfs[2][0][0] = 1#0.997501 +sfs[2][0][1] = 1#0.0744126 +##20-30, -1, 0 +sfs[3][0][0] = 1#0.999754 +sfs[3][0][1] = 1#0.0258297 +##20-30, 0, 1 +sfs[4][0][0] = 1#1.01528 +sfs[4][0][1] = 1#0.0218661 +##20-30, 1, 1.44 +sfs[5][0][0] = 1#0.955474 +sfs[5][0][1] = 1#0.1803399 +##20-30, 1.44, 1.566 +sfs[6][0][0] = 0 +sfs[6][0][1] = 0 +##20-30, 1.566, +sfs[7][0][0] = 1#0.977028 +sfs[7][0][1] = 1#0.0965378 + +##30-40, -2.5, -1.56 +sfs[0][1][0] = 1#0.99054 +sfs[0][1][1] = 1#0.0131795 +##30-40, -1.56, -1.44 +sfs[1][1][0] = 0 +sfs[1][1][1] = 0 +##30-40, -1.44, -1 +sfs[2][1][0] = 1#0.994976 +sfs[2][1][1] = 1#0.0245861 +##30-40, -1, 0 +sfs[3][1][0] = 1#0.997789 +sfs[3][1][1] = 1#0.0103248 +##30-40, 0, 1 +sfs[4][1][0] = 1#0.997776 +sfs[4][1][1] = 1#0.0271587 +##30-40, 1, 1.44 +sfs[5][1][0] = 1#0.996409 +sfs[5][1][1] = 1#0.0219291 +##30-40, 1.44, 1.566 +sfs[6][1][0] = 0 +sfs[6][1][1] = 0 +##30-40, 1.566, +sfs[7][1][0] = 1#0.994255 +sfs[7][1][1] = 1#0.0162298 + +##40-50, -2.5, -1.56 +sfs[0][2][0] = 1#0.992323 +sfs[0][2][1] = 1#0.0106739 +##40-50, -1.56, -1.44 +sfs[1][2][0] = 0 +sfs[1][2][1] = 0 +##40-50, -1.44, -1 +sfs[2][2][0] = 1#0.996538 +sfs[2][2][1] = 1#0.0076656 +##40-50, -1, 0 +sfs[3][2][0] = 1#0.997529 +sfs[3][2][1] = 1#0.0058141 +##40-50, 0, 1 +sfs[4][2][0] = 1#1.00076 +sfs[4][2][1] = 1#0.0063664 +##40-50, 1, 1.44 +sfs[5][2][0] = 1#1.00054 +sfs[5][2][1] = 1#0.0076592 +##40-50, 1.44, 1.566 +sfs[6][2][0] = 0 +sfs[6][2][1] = 0 +##40-50, 1.566, +sfs[7][2][0] = 1#0.991665 +sfs[7][2][1] = 1#0.0109673 + +##50-200, -2.5, -1.56 +sfs[0][3][0] = 1#0.984755 +sfs[0][3][1] = 1#0.0112026 +##50-200, -1.56, -1.44 +sfs[1][3][0] = 0 +sfs[1][3][1] = 0 +##50-200, -1.44, -1 +sfs[2][3][0] = 1#0.989592 +sfs[2][3][1] = 1#0.0121205 +##50-200, -1, 0 +sfs[3][3][0] = 1#0.994825 +sfs[3][3][1] = 1#0.0125226 +##50-200, 0, 1 +sfs[4][3][0] = 1#0.998808 +sfs[4][3][1] = 1#0.0135373 +##50-200, 1, 1.44 +sfs[5][3][0] = 1#0.990307 +sfs[5][3][1] = 1#0.0264489 +##50-200, 1.44, 1.566 +sfs[6][3][0] = 0 +sfs[6][3][1] = 0 +##50-200, 1.566, +sfs[7][3][0] = 1#0.994873 +sfs[7][3][1] = 1#0.0096489 + +fout = TFile("csevsf.root", "RECREATE") +sf2d = TH2F("csev_sfs", "csev_sfs", len(etabins)-1, array('d',etabins), len(ptbins)-1, array('d',ptbins)) +for i,pt in enumerate(sfs): + for j,eta in enumerate(pt): + sf2d.SetBinContent(i+1, j+1, eta[0]) + sf2d.SetBinError(i+1, j+1, eta[1]) + +fout.cd() +sf2d.Write() +fout.Close() diff --git a/Weights/MVAID/MVAIDroot.py b/Weights/MVAID/MVAIDroot.py new file mode 100644 index 0000000..23deb03 --- /dev/null +++ b/Weights/MVAID/MVAIDroot.py @@ -0,0 +1,121 @@ +from ROOT import * +from array import array + +etabins = [-2.5, -1.566, -1.4442, -1.0, 0.0, 1.0, 1.4442, 1.566, 2.5] +ptbins = [20, 30, 40, 50, 200] + +nbins_eta = len(etabins)-1 +nbins_pt = len(ptbins)-1 + +sfs = [[ [0,0] for x in range(nbins_pt)] for y in range(nbins_eta)] + +##20-30, -2.5, -1.56 +sfs[0][0][0] = 0.989453 +sfs[0][0][1] = 0.0400067 +##20-30, -1.56, -1.44 +sfs[1][0][0] = 0 +sfs[1][0][1] = 0 +##20-30, -1.44, -1 +sfs[2][0][0] = 0.997501 +sfs[2][0][1] = 0.0744126 +##20-30, -1, 0 +sfs[3][0][0] = 0.999754 +sfs[3][0][1] = 0.0258297 +##20-30, 0, 1 +sfs[4][0][0] = 1.01528 +sfs[4][0][1] = 0.0218661 +##20-30, 1, 1.44 +sfs[5][0][0] = 0.955474 +sfs[5][0][1] = 0.1803399 +##20-30, 1.44, 1.566 +sfs[6][0][0] = 0 +sfs[6][0][1] = 0 +##20-30, 1.566, +sfs[7][0][0] = 0.977028 +sfs[7][0][1] = 0.0965378 + +##30-40, -2.5, -1.56 +sfs[0][1][0] = 0.99054 +sfs[0][1][1] = 0.0131795 +##30-40, -1.56, -1.44 +sfs[1][1][0] = 0 +sfs[1][1][1] = 0 +##30-40, -1.44, -1 +sfs[2][1][0] = 0.994976 +sfs[2][1][1] = 0.0245861 +##30-40, -1, 0 +sfs[3][1][0] = 0.997789 +sfs[3][1][1] = 0.0103248 +##30-40, 0, 1 +sfs[4][1][0] = 0.997776 +sfs[4][1][1] = 0.0271587 +##30-40, 1, 1.44 +sfs[5][1][0] = 0.996409 +sfs[5][1][1] = 0.0219291 +##30-40, 1.44, 1.566 +sfs[6][1][0] = 0 +sfs[6][1][1] = 0 +##30-40, 1.566, +sfs[7][1][0] = 0.994255 +sfs[7][1][1] = 0.0162298 + +##40-50, -2.5, -1.56 +sfs[0][2][0] = 0.992323 +sfs[0][2][1] = 0.0106739 +##40-50, -1.56, -1.44 +sfs[1][2][0] = 0 +sfs[1][2][1] = 0 +##40-50, -1.44, -1 +sfs[2][2][0] = 0.996538 +sfs[2][2][1] = 0.0076656 +##40-50, -1, 0 +sfs[3][2][0] = 0.997529 +sfs[3][2][1] = 0.0058141 +##40-50, 0, 1 +sfs[4][2][0] = 1.00076 +sfs[4][2][1] = 0.0063664 +##40-50, 1, 1.44 +sfs[5][2][0] = 1.00054 +sfs[5][2][1] = 0.0076592 +##40-50, 1.44, 1.566 +sfs[6][2][0] = 0 +sfs[6][2][1] = 0 +##40-50, 1.566, +sfs[7][2][0] = 0.991665 +sfs[7][2][1] = 0.0109673 + +##50-200, -2.5, -1.56 +sfs[0][3][0] = 0.984755 +sfs[0][3][1] = 0.0112026 +##50-200, -1.56, -1.44 +sfs[1][3][0] = 0 +sfs[1][3][1] = 0 +##50-200, -1.44, -1 +sfs[2][3][0] = 0.989592 +sfs[2][3][1] = 0.0121205 +##50-200, -1, 0 +sfs[3][3][0] = 0.994825 +sfs[3][3][1] = 0.0125226 +##50-200, 0, 1 +sfs[4][3][0] = 0.998808 +sfs[4][3][1] = 0.0135373 +##50-200, 1, 1.44 +sfs[5][3][0] = 0.990307 +sfs[5][3][1] = 0.0264489 +##50-200, 1.44, 1.566 +sfs[6][3][0] = 0 +sfs[6][3][1] = 0 +##50-200, 1.566, +sfs[7][3][0] = 0.994873 +sfs[7][3][1] = 0.0096489 + +fout = TFile("mvaIDsf.root", "RECREATE") +sf2d = TH2F("mva_id_sfs", "mva_id_sfs", len(etabins)-1, array('d',etabins), len(ptbins)-1, array('d',ptbins)) +for i,pt in enumerate(sfs): + for j,eta in enumerate(pt): + sf2d.SetBinContent(i+1, j+1, eta[0]) + sf2d.SetBinError(i+1, j+1, eta[1]) + +fout.cd() +sf2d.Write() +fout.Close() diff --git a/Weights/MVAID/ScalingFactors_80X_Summer16.root b/Weights/MVAID/ScalingFactors_80X_Summer16.root new file mode 100644 index 0000000000000000000000000000000000000000..20b03afec35a9978344fcbc6016e135dab32491e GIT binary patch literal 5612 zcmd6rcTiJ%l*dB}p#(yaDkVmGP3VH4RKd`jp(qKx3DN|CfQU3fDbhkwLU|}im)=5A zq$(f^2m;cTB2_@j242~jH@o}XKiMse}Czh4LDio2I{__ zHbC|Fl~0}nnrJ$^*+Y$OF>cQ8j+z)-9}h2YYk4@z+SvEzO?xkdthm=ng8v)=0G_Ok zx{sY8CO$}9n)o)S2l2$z|4)gNJUPja?A~d9n3McSPV%)QH_?(*|8sNeKl9VVcpEtT zTN_<~DnsGo7o`6&z(R*ongno<4B$`VjYgN)qHm+o`|Si$-cm%FA{@SpfRmo1hfBe4 z36cv4%-hZ9cQ)|xr_Qo6cLwsZfb7gVK<2_DLTroX>GUQ5>4dP*1R-iW;rKoZ$&5X| z_*bt8{JwO+gtbpTKo!OP%q_z@Y>#ayDf-@b<@ z*Lic>;L;la{aL|W;4mGvhWAjlA| z8a5gZ2r>gycxZTN=sFL|s3&0ESKRW5O*f-7a$kkw?#N|!v1j7YSA;2m2?Tg6<vxzx(7}{_pAcjfzIZ^$DC;Dm0VciHBJ=O45L1o zF1;OEBag!*YqtqvhO8gDMEUmmUAu>eAlwk0R*(%2V7Y3UO=MGa$|j1a93;T;L^5z>6!GyKle;n#3!;`cYMMZZkAAfwoG2EFb2%n?-KHtk0t+RPJC<(#6re~QZX1Ql@UFI4CLj!Ml~!|cT3{sQ&S5kyR~ zB!K&WVnUzVTqn`yJm4obM-KnP<~#){{$_KEpoWN7M3eKM=UYf(aZJf)G$WCI{JYPI zpZHwXfA4eupOln;#CKNUHeq-ayMCvufO3KH?JYF=Z$g@UDs!E|&>2+Ob#$mhp9Y*^ zLXrX6N--{{6_!|b_{L7l*yR&KCv_KK2Q<-f_VxiA`*_)7h*eSB-NA!V-NV<$?GBV! zH(l-Rh$n9!;^*I|0gUkxtIV*4z8mK-`|<7NZT=fG9nnnCB*{6CFT5SHNe{@|Mm(*I zRKfmaw@J@cHb{Mo_H)d{jgCZAyA7))EaG1#3!CU|jUb1S2^OX|Ef(=wD^$_Qz(o>M z)$a00*udnErqT7@3fHytTZ@&S1EnfGKCfXzCVw=StNg%fnRnAacum_a9wxaK9AT6O z!S=i`KVtv5L#aQB&!QL*jOd}k;wGdGf* z&yn*FPEvLU=LAqYjDbai_a(RSO8SRRYjLeng4YNy%!D{wwF>es$;n-Nr*!8YPl2JF z*x``zW&wiKH4hF1D{$4cy^~0c*trwhEJu}9h3jVs|4q@phnDM8c@QjJwTZUMwhM8d z1WkY4-C&-4A2|3$eDe>Q|mO&rap;8ApNaT)(nQDc^HB{eg*|rS65uwEb*>#&%)J- z!Scm_IDeS9mvH!utB0^B<-Y6XL~qKGCx#EfNxzJ_=MRrgt18HXk6iOn<+}GgWg9Bp z9*|uc5RA(nel~F#(YhvkKXgk3GgXYZ%CjhZ(U}72m<7#)GzZSdr6pPKlH35$ELEAY zP+wXPdr0zez&c0yCD8m3yWXgtB=$Ws(4jG)yuRcbm*1#xE}2fom&P9jll(k#GKxK} zcUtf)$ARVtpV`P`u0hWVLsOfbQhUn zQXegR4sEEPOA)O^Rc#p^SUpWw4j*_b^3+Vf@oCTqug1uJO$2OCLbdyv%}Db^aoPLK z&RTZqIhwQNSl<~VxfqqjjC=!$HshuILU!`5D!EeY-H+Cj?l-(A^PWS#bR9I}V>bm< zhgP)H6U)x%lKOUuPBY8C>nT>frxmWx46bz_m+sulRXn+4ZBU|K`_jizaolmE(AP*F0nvovzbS6UdsT{{u@call zmUW7lg;BjPb&Lsyv(qCs`vct30{l7WEc>WUHxHD)rl+p$W#_sTW1?70eO8dNsEwf8 zNpqB9b;%QxbT^f)IyIT_z@RW(L#L3YJgb&ey0}=C{rvetl=>b&xFj@DW5R^Pyx}WO z6r%9SKpCw$TV9S+$?n4@3o?Cv#4Z!I*dO8L^wT0-#R%LzNuXSQv=?XJmXIcuR#y6b zmR!kmX~G)qMyGZZi;s6Z@+*Ff;E-6A0oOhXP97=#s9qs3!ycVJJ%A!{02EA<*VE$D z6*8_|cyN6yEXYqiwdsidhC+C{@bN)k4L4F*iz9Es0J&IOxr{DW`d0fv+}%e?UDW7B zq?zT%?V<4{sg-!?FsBQu!c^rhAZ{e1h*qUul|uJJcbNY;>1t)%>)>%JI=T}^|$p%(&6k(Y>Ay|QR_6Vqp{EU)Ko_= z{(SvG4o_z0g*+bR5VzIzIkk{$xb}muzoqvFcS5O+5FawBUDgj^&Q7MHnbd3)%Ruiw zDoqxj)vXHVq##GWEqNH0b9?O>7qS{2o(j$}3fthxo+NiLM6aB3v4{+|K^6E!xu)Pqj`;BH4}D!CP|hfMoB*N!vmZyvz&3$oS7ZE};< zDB=1Y`|Mu)b9*0#=ov86d)J~|poDbMeo)hE#m&MwlN414Cy0TH`g=)dK=u7xh(L6= z(6g=vBUJ0KLM9}yRu@Hnp{>U*{_%zn!}3`3nHMlmH@qB25d+HHsgYNjTX(QlZJ<{6 zj;rO%8|!00(;4T+0!&O*z{=(K-#+%bOfWK;KM%mXah!Y7F+!U+Bdj$Xf_WoWmtOPs zPVh~ZfyRZnjeK+88;Cv1+)o>2a*WmLR*b@S4DySy>rl zn#ftCWy2DXT~BN9thPYX110aBaq3NfJE|*|rC%}#$siqDi7R0>pT7TOhXe}Qc+($p z4ZbU(UaHg{1g(tLiwF!)wdM-s%B5`veQ@9&>3g~zCc8mZzm&2afZpRwQwr_=Tx820 zZ+wsz+J4n;I3{n}{RT_oSj*;)&qTG@)7%;Sq6Ngn|KjBdAgxD!8Z;V54!ObXh(n%p$ks0J5k0)h{AoqJ==m7cJEjh~YS z*&HorFmBe&s>`3ri3`SlVd!mP3Y$`p|N1t?GWX~l7(KV&61|VBxIXU25Ye)TeQm~_ z>VO{ifZ8$(ZEWUEud*}DOlZe&l)d$YUTkDA{VwKvu1IjJSe-J%ZimI?Q3Hgb$~{>m z(X5P%XPz25RhH;p`0Om6Mp7u}ajLddlsF|qY&0f3y{@{Yx#C`sDiiq5lNC%ftDa3h zV@lK%&Nz3-;C95lBza8#a~isH0qYHOQ6;mu7NBd&1zCn8llTwJkE!yrYcS~5j zf2!|g*Ea@1g|J0t60Gp^1+_Q5A78Lnsp$K!A~?Cyf6BmMoRbObiCWQ>`YKL1`KeEt z85KbV9Tept)BeYcpHhx;7EEb#G+#a}Nr+Ii{*b5;J8v&m$ys$D(wm0&ZOJ}ga~E{Q z$91Jv`Dn2`y300Bxlmt+J7jTsW@foRa1^~lv9fn&v8(vt3GXE@dS8atF1SDmZm@8i zl)~etdwiHDBNiN1*-kr`1?o^-i&M!?PBfwiaqS7#<{gEgFLL@EM;6m5K+u|?=)9Rpo~;>9`Yr-xXk|7o#T;P?Fx^gZC824$EhBX)XOk0!zv=|P{(L5 z+A!H6KB1OF;g2A0ulh4xA$iXJ#xO%0T5Km{!QnEcgivV`%^Q(BR8u8va2AeSpGZ3m zEr|pN`X=ZYPQjT@IM2MJtiesXBHx4TSyqkSUg$Ys9^F>Mu%U8(^!0gm70nAU;UyL0 zCt|P%?Nu1Mq3Y-`-t=c0(bCz?yuHY*ospNw9!IenM)8Mc`L==`Dm?+Tb&ek%tagnV z(To;NGv(S$FH+E$ADd`a2sDam$Ugg?_AoHrZJhr@a#S+c(AVCp@Ax@dbcSA$KXT73)2E@3xeYF z&UFzu3hUZ(KfD!uywxR9iUofuw|3RX2#(}tx!6N~A{-!%XRcp~V`8;nPAPDUiAw_U zKkAp4!JhXl+q;5t5$coDrK3!akc^m(K0mh{xqOFHQ7oB~Bk>z|g{tDLNHgDY^}GJb zf}JwKRMWu%&wa!jeDfFE?uDnULrPmgC|WMN&IOGaXM)uPQOlDh+mtu5o1rD|99tM x#0~!4kq&W#|BJZH1wN;h~603`lXwed7k?2Hc5iLZo zi54ZS=t9)sY|NZ<=FGe^=Y5}yK2Jooqd-ZR(z&&S)t!v_G^5C#AM2mpXqo*)uJ zElm&+L2`c=gn}3e08mQ;h)O3_HIzU%R#Ta}MC+SzJA}Xg*5e2G+Xx14lHdjiUv(y= z0RW`N8qRL^RCXwPKa_)r_usSpp9O&QuY({S!e(27V4%G zyGot@>nAbR64#*k8)}VmK>bf)p-W*-1n?#WpowaXE44dp=q^UN()|oN)#*suy9V7g zP+_sLd0`A{Y8#cE5l}!9V#jFXK|TW#WOwA3e1a#C00N}$K>c4@u*zio=o3idXdJk%!E*; zkFp|Ef(&Qv;7sZz5?zlrVCLJBZ*#%H8=FFhQXNB6bat8YQn)$z+j&sllS4 zK*5*i2d=-XLRv49f)S$-wSstEsQBqC)0|{X=wcGqBVA-`#AI8?D;{+Id`7ZGG3<Ex|P_4r04f&mn1BR=83sxB@dPB@BY5i3yT?pR$J_of(%u{9?%+YM=HYAW7QpW2 z?CxrhAiPjMgb8g@0p_>~6PmD?^8lGLLgf#qqo-=TmhzF>KMhV3Wngwq;39tHRJ@&J z0uOe+_Wo#o)?AC7zM}kce}A{%HJ_6NixH;2NHz8E$K(dqMHtNU&|k{jXl1RZ&&U(d zD)*Y(bkvWtB>a0MX681_W*}arZ{#JO;4ik%vrt2d__$b~qqC$!P!4m@Jr>ai+tC|* zSS+7`i5_X8!VIrM7IscOSJ6pClzKtIEYQ|3aX$Dgl6-d+Fdf&M;rK;#fb>~Sc-G$A{MWwOMv|-V;_(AtRFfbEaVH)MWhDTF|2P zrCKY@nRw+6=kf+kE>%$-1{X~NlvtKcWaB2YSBNcuDnBhqu2%53;_@B|)qKjC;&8qH zplftFb^b&wP+uFh+o7I2V8SAsJUZD(u?}Hz+Rdc%*B9r9u-2psq!haq)1R~#Jop#`JL9e(^JQ&P2u$$vh34`2gdU_a($CQmr6 zgMO2e*^-tB*rbwdIZy%ew*v9E|o?#)$H%1Lfj1HHaCsw zt$O=s@=4iWDfBmCf@VhTY7c!?%1#X}@9vJpC3E0*>FHEJ}L~*dB88w5zdQ@JW*>%sWYj1y8+YXZ$8yhpi&4OOpo7WOo zSGkZHI@u;Rx^y{GvYskW;C^(#p*Ae;&c|O63+{S|LoL=|=?*B~IB%lK_J9U=fo^+D z-sx^lzl&~|%^>@3|MsBb050RjGh)LC&B6IaF}^r>&fev>V8iOnH_yQ{vq`J#1eUih z%9$k5rp%+w&p8;BexgYx*!(bZ?-@fzhRxfG7WeuLr1b2{#p>6Mh^BAB9ah=vH;1M& z%P0y{Jms=heNxZ4CT&X}$-i^mk6D>a7E1J;dI{;V{OJ+XqAd4_nq>;iQYBbSdKL#B zzttt3@k8f}Z<$5n@)P-ABad9|TJZ{dXRh}SV(89oQr2$bm`;8##DA@h52{|rY?Rg5 zn=FBCFlBtHHAiBgx`j2fkEu8VuEx=KBOW=2pW6c9pok0WS-(S_P>B7E7tUn?saV>^$H^vBPe0@6s2i5n@*|&hW#oe%pOgYL<-nwZh3*+2INw5 zMxA{m;Fz5<(>YmjPbN)s7U@!(p(&&#sm>r=@O(f~Cj1y`1KOjv3! zpEsRY)PCh%@($MRw(xOgk@RL+w>S1Ij5qX_>lYD`%c1tywj#!=RotY3pG_r!!+ll#Ov74oMWV{TIkJF+t zhkB1+#kEsQeNmT8`}P0}&wdQAPlR6faQLe@_md)R4dw^PB+AwuR%Hrsr zW0xukll+e4i(--^()IzgT0oH>!?vJ&BBgmf+;BC}$^BS`(PsH6kvPgzQa1TjlTiPb zzC>#0ll(B*f^CNB>T5$a52lPCf3}?M`PjUp+t_ui;dULYe_RHRs|RVB9MePZoIn8M zNS!6{QB)-ZvO%Dm;)Fgcc$0Yv)5E+Jb`;;f`F$knTWbUG40Q1-9@(yR4%{;Lx}5C8 z?M8r;9H3d+eKI>=3x*3TpzgfwzL$`-8(L}!#Op@a`$sm?BwW{|!}BcZ;h(1OlWynR zc>{%q7HK9&MC+=?#|;rAYMNg#s}-TmiML^LJooFn0sxE7F9#_pU}i<|nrNOD*9Y)C ztT3{>NO}ZQV)$Z`*C8)ML_hu#b?e-2h|^m zA*w$rrE_S0MO70WT+F;wz+5Cu4>mbb8~#YNQdecTL9}j&pQ za9{OkyLFOX10oX8;TQL{3?%PsQ_1mI1VhGKCUIxUh$^Y0si~>N6_@9W9iIj(xaeSZ z_J`_Rm1*yL1CITvN%{At=u1pq=S>ISd-);TG_ax8Mp)eLwcl&$!(YD5q(YOsDxll3 zx8O}H>0(p?XLmnO&mi0!`BHzofy84zk4E>|dzY#%8TN0}GT!SF7-x6c<40qwDmTo! z>3Iv@v;aq0+r=R_u1|^A*S!*ThrQug~T^VJgv@<6P4! zUv9r)gk+O01S?5NJlR}~9C-6>`3UM6l+KwZ3u@_3@(b%moWC~YiZOC|tMl?PbFK)d z5_p;KrQ>>^F#GUDH4;MT4+!$ zGcn9w5>m2e9*#iv&~5XD^qy>UM3@Jq?64G15v`R!?40m-BlbU zHs=mdM4ZG;`ObEKF2mSsEXpl>I`aGBiK%|}TOWU!Z_W$3pvlhfl!~oMLQa^@XtI>y?dxox0qsJglwKBgfPhI~LzXgtJeE zNE9U|Pg|lVcSz#iG4I^%s(k-DN3Gzy|Vtljk*sK4&!F2YpW38xOF-U`G z^w8>JWu>9@&ta(_iZvAN4UNGkon`m`oIJ=Jt;m-r+K?4;z4c(b%YKG|z`r(Sox_IOh@`+ilZ^o_f>zA?bUCG5R zTJlx2Y2o%C8c)CWul&>~K~X{8rgXl$C{b>Rx}@pgqo13PVj9Mod`n2WBo0!f?!O{G z0yF?5ub|=Xou}|4$-2KnVNLgcL%(=#T*YBD0dEYtxyuJMafQfei0Kf?V z;L8C3c#rYNMfeW^{!LE&=K_xTCY5mzblgcmN!K0wxY~^bT-x_Hz^u3JMa} z)P{O`LdD^ZP!DGj!ZOH5 zUP9d2<+7b74CW^0AMEd-sjV-mA?9~^%m0dN0!S|J|91=kgy75Q{=YZ?0N39Bcfloc zxQrp`^6&Pye?Ogy%l5>kcOiI$`Wv$Od1-}lY^J6_zdKLcj1N6ttDSSAGo z6Tl&m_7NloiTs2@67o8-%yED~A}k;TEqNXbNKQj4X(cg}0dBlRpQ*c-m|(ZW0EaX5 z8y{p(sR2pq3vaV(&{yO~Bd*ls(+$Ch8d&li2(4@(_vo;sq_pI|5J+nu@(Qn6TJ8PX zJ2gWh%s?^GJEDGLH$=sBM2(FMM0JR(9ib4&8JT&l^SN~Q@a#Ius^E3(>DJWbq=zoe z0g-c4BtM)3P(?&QCP?Mv*V{)O$r}m58ForA#}AibyL3m<7J>IiE^?*RAj}#$h;GjJ zEwS$x)Mpe`nYi<*rU=Usd%DQZ4E3L8d_xwCBtK_(-cWqlu7v)f8XyF0J!_0a5Dcr- zQ1e4dS=fId-f~A$iLK26n0EpuaGO#k{C@Q5IhVR11N8hWW70JaTt)* zx7A9@IEa`cAJW6paMNFfL@WaWXvLT$e6gw8)~Z_bT{p1h00xZ-9&&z9tHaxG9yn4h z?aT>hP2;2_y&j*8+vg0wzreLQ9%^mv)S3^C$4p81y2>FB)JD~UYCB$dW^7I87ajJL zdzvgI!g@!qs&M?gwbl5v1I2rzm=z)A!^%%B#Z^w@PW|LW2%hEWLPoJ5>(&!WcGu_>W4z1Jf;a7HI_?{nzo9y~x zR6_YrwvrRNT##8I>n|uFx({7l3N-ivnd%yb!jSsF3ec z`IAt5!to0N-e}k2FkgWR|2N{`A$Ny&_7k1Qxeeh~)~uEy?drW?$6PS#8n%h6@bLR< zN2DKL+vs5#lM7~~taCmWEJ>F5&92a~IH;AV4m_2;06y5c8+Rfid)SXH#OSUnZRt{H zNA>$@^lwoLp?ha?Ux#}Yw}}#qSOlR5J5xqco%RP=;f~C$r1j#pA#P}iWiI7uB7hNX z+$jJ+H1qG|rF|(c60*O@i}wG@3qO)*K>eZrmBt?bug3oQ8cX{UsY|bUhU^1LaA}p^ zA)2EISl3`a6BWJg zzQkH7(dP5^Vp{aJQRw%U&8P0=`}6yvr|*0xm*#!KdotTtbHXwjbMNFIz@6P&J7SXG zb>?h~g)W~%({+_z_0IhCbi=MRx;^Vi@{XH6JE#M#31epSPoUSjI#4-!`{IqZUrwdJ z)!hv5fOv@y>>E$?TWMQ61);5d zpH;_u;OEBrq0#4oR-8AseFDa&8Z4yAwp)U1FspELQ--Cmr$6zN0WpguqK$$)^05GG zq&$s2C|}o{z9s09%J^sJLrJ6W?py=c+_)M9M!;Wpu43baamf3h$EGldZ+xh0p*Za9 z;e|0J((jJ~|E)|hr8kB}pInOp(=?LTo6<+Y*LJCF!_&lyJ6!tB>TUJRl`ZXVH>QOz zjxRuNc1bW}1@mQm4sNKn$i00{G*Y(iAN}4HJ3YvAXx01}ax*PdGC>wO41>3v%fg>LSFjWbE+`0p{z*6Qp|LB#yTWoA03byCJ5g(1ikf`!FQTUYH&N@W z|5uWhBK+TJ^PfI#{xwN^*`8gf%SBrTkaN&7|8W{;cFP_7kbEM32Z`h!LQqAk2=q{-^QgYew&UD;c4 zr!T0u<+0KwEZ;YgGLhlCNBf80<(GU8+AfZ^91qI6f2eug`3P$mKYlm|j_sutVLn&- z2*!2F=qcZj28i^JCgB0py#$a%|F>zAy#$cb?h-(tskR@~$5qqM4{v43UHnY6ZfZD4 zko-Q0H)=pafIi8eM(w}}5(t6F4fYxT(6)zSXu7&Vksf1sWz5 zfZ4b!uI@JUL_kXv+p)yWg6@z|8fM`7Q zu0jPEC@Z?*$}t6b1N}$%seTWK;O-M>Y>gopVsL^yKDn~OTf)81Q%_H$ev8&cn5dM9 zC{uz+(MZOaTmni%K$_<6LoNZif{-q;Wg-$tGeyMmF%XeG(7dk+`vcH zau8kap{mI=`BF?b10fmzKqLoIDa+Ptl>MpPF;*T+NvqwH_;hM8r>nSKhiojiGCE#; zM(R_Gfk!3oIll+GV%_3EcXBG?#pAe*9zNl%1+a5yewH=zq*fIDsPV>=o_;pbsA#J< z>;!%HF(i1D9{fj?sDAtZzlze$zn~=YH=%@*LNw^9%CZqL5MAHb#(?a|v6jR^?juVJ%&_Rk!F>44cg+%wA^Z5AV^0S>?u)(r`FQ>VSENG8K_77t}=mCe-w@RWHY_WJEk7{e&;ApM;s| z4{r#7fO_DZOT`;*7yKiVh`}-viW|-6mKeD(Y?oz7T^9=laaDP#n4SSll}pdtV-BlU z#ilUmp9}U7H@cQb77;y-)kCNYy+^dfk|YQdFl7eNK*?JN z5b%re0A#OJu(A^Z>uCWAJ-4~31TFe6v)#xSY4uL`zd-)BlvfSvZi64o1WWso$TYvm zZ$R?GfdBC@#umRaJeWthh1S+@6L8pP{r@_Q9mSHngE$nO(DIIlPpzHRu{l{57VIMO zTZ%vCq(zltFg8_3UW{;C7qqv}$f}ifYq<8sthM;Z_`_$~7?zNnqAe37DiaJbXC9pc&?98W|PefEf({O7QU7!^;8%QKC3Q296%t)wbc^lwCZzLl}41(|g7#}}YWi_fxKrq~?0%-C;hx?G- z*6Eh%Rw5u<>on`Hg(lXv8P+dL3avB9%?Eps$YYus>ycBpsLa~&{6<*^nuF+`>E>TX zDn?p~_^QS?zz7Si+$Nw5aUX<4AfAIf+C>%vAvN{pQ6B4lIwCrz6BfGx<^UlRDEFAP z7WUQ+&e6^{+FL5p@r!1|im6ZmBY>+t|C}|H_W_G#wd!mXZ+tTl)_P6`qGTIABPF&PgKUlWdR*O1UD_5F!l-BE~wYsM}qoAQh12_uocFWa1 zOF21=u=G`q%=NeK=vIdX(cMiq?6!9q&)DgQ!)K%{lx@tU1OU1D?Xlk#OuQQtEB@`q zR4?7wWnB~CL}IGv2KNV=`1?6SJ)Qk@ykOq+8r}i-JwmuW+`Jw*JK=xe{&-*ZgbGkB zhWBM#)^;B37EEBppPy2;2IPDvXFb2uo-!!vIJTEf{PNzC<3~u}I4BRu;2*_oCn*nO zQlzBV7QUl1z(yOol@)?bCt}+DB*2>x6GKZLlMo}6y;B$(q|`W{ z9HeetN&oR_o8qeCk>XJyd|Y{bv7ofw{*{O`dH)x-_iOCV)M0TM#fgLqLJQs939AgB z4JHTJFH+a`x=9z?_sIBNzK`B=VC&9hn}B78g}n5E36a;Bg=un4kS=lRw~`yBY}}@B z`b6RBWxwLJAs4XR$-pHfczY+3HGffgHD^@wR?^vHq)-7#m@~7oFz~}l!`|e!_C;3M zyo|C;3zvASce|GVl>?pbg3O>6zMl{GBxB~)D)rDeImTOFl&cP6(KxcxWL*xyXgU^$ z*L2TR&&{A{qa;N7&%8S(S6(EUxSk6A7-S2{`Bn$l;AEAUu~?bx8`99WeRw==qoT3n zTgiaVd~of0OP+%ry_&KnDJwg}z!DuzLL6yf+lz1WDP~D$5Y>o+9|Sne%+;#3vY>(2 z%u3iwL2dtLEv+j7>A*!pifK#ws=?vdnZ3Oqv+ua9isUUs8w+zhgP=EA6}32e%%lsi z`@ACY3Zk#G=+(DJGwzhC+gf=h$f0?zNxMQ`@_&-6Z>y9wNSh@k@51f4Qg~u;@KV{& zi0hV7`K zWTnuoxWnbU*(zsWb!m7VrkSKOuluDrR{nrscnvau<6AR|MB$wK!8abiha$n-RF(Se027mwTF4uRP~XnGM0k| zX$b6x5t4;1l#Q^LoV2%}a}K5+)T%2Xuknk!4!!Db)o^knEGycr^4^j6JqLodJS01| z0-8dJ`iZ-;W+gfx!~)6`*9(=nF#vY?M3EwA;HurSlLzAB;@bZPHRXqOQ=nE}FljuG zcl$9<2_Bi4xO~JeULR<(7>j3#W1PPJU|zo9C1A^s6j(HF3NGkzSVMoLh`85V%pCkp zV=+HwC*aTm0>Joe_|yzGKUC+aW3C%$Al{7GeW*=eZBxtKt3OKtd!Q{0HH*$B@lZo< z>qLFe_Su_dv<|MLkRo=Qz1p9GSg?z-cq-Pt|170(i7&Y9gkDKy_U4NGsFj7B4dYmf zfc(OyQdGWGU5>WDXHm<$Dp3a5@P6q~t%$z@Do}MPfza$g24I)lay|xC5gboI%c4|7 z$2%aJbq;NJaa;ZTvf)+N?x`84tK`Ctz+LnwGhfWV1!;F$`StRRTgxQ=;F8c=_e_;J z;jz8-#oeljMkVvQfmJYu zGW#vgusMm|T~c-d?237C@j0T?7M2x)XF;p(<8OkM`$>tPdi|6XtBCIs7=Fx2CFgfd ztzRcIH_@z-+aD3*Yiq&2J@8Of>Pe3BFdffn!4O?Xg($6%k{zw@%Xb_4V@)>si8Zc| zXjP?}g`JIrn{95YpOFCr==}2oX^+)hjoFGm(oP(0FQwLfs2%4}7m&RAfS%Jy>`t%* z`#w=hN#L4J`|56|R7jInOq%>Oh^HPsDAXKDnroz@$<8AvE>^o4QkJ8Bjd}l%>)Lw? zV%#s>Y%BNn^m79R9C)uIjFaoI_fj@=Bz+#)&gD^TG*XjYS6BGqqf7g8vo(+6nSRw2 ztBbu4Cj-0G8iPXb+wte9NYaefat0uCud9CivedOFwZ2ER>ZxlWH<&8*FD~g=_Ko6X zKM8UYoGT1#^1kzfytVJR*;cgFE!VfzdmiP-{u~e6R|+@5)<^3+A^8i#(rWLO-`}n- z^?A^ptaqndIG)o&9X?NUH5j=BxvGcKz*rWtEM6xjKF+J_2CMWg49Jq$(eUxT(QLV)7i+F?A z@=uOeDbyE#9^TBuw#+gfz4q=7PkEvf`&KLE&MT&khNh@pDk_iB^US`{;xv|=y_G+} z&kdI^N*D$|JxeUn=cZF5tWcMiNVTuJ!P&vB^p?k*54EnrZm6I!ST8lWOR~9PVf|f@ zp7=#6$!)l0huC6JKL^A5vZ=M|ag)_6@)(_Rn;Y}Ou9Kbo+0xdNZ_9xp9Qt`OB_9;N zGHMx3$PEr7GvBYjtO@SoWNIfq4*-3mb8_R@%*qCeeN8`)(8qr@4>{_hWQWwtNH>;l z?Q~r|=UIJ5H351tROH6k4`3~6F*M*z%H)g=vk@@+^%BPu+Ui^Ry0OM_^HZ(Yy2Ps| zvbCCXXK|f#Zv*?oWl!h?qHSwQL>)J4q#Oc;ZPJ}3wvN4 zS<*)HX0zO!wBxE?==gHgMB0lqxv$bJZKq&ucRLOH#1k1Ms*fq%+`^ylb{ERXw{G$7 z`AD;6rfU(+ioiYRD)&>Kg>iThREHz7rHys>(*mdwP@5mN?bd|>Q*=u!W_n(fo*IdK z73$MpEOrZ>8$UMu;A1wwI%inA}+73NTfgejgu3k!`7kc9s zD#r0HLRLqehS*LuxHCRliL+WDH4cA(k&A8AUnP>|{6)?hiR3Qh6&x(1={kY2-*`I@S*PMhl*O5WWX`H%cKhfv z{)|_S**a$xdOp(&7p6a@nfRWsWEuHdWF&xNwfU;c0K6 z6{QX}z6rXIYlN}bsMWAM1J+tu=1pogd&fC5U61_;8n&(o#3*RkTKI=kr02J_jiDYS zWwVOXr7HfJwd-COuJA@96+l%@}Cvdk~^qH>$=u!u=xZzeGUV|$mu-a zkI%1jJJA#*NNnF)u5`CZXnveY%KmoqXl?ozxT&)y?8_1^aL+khu;&eWX4)u1RJ~39 z35f+8r1G=jwslgaN|$Vj;p=nh+S?$Kq2*9oNh~%VTy|vtVB8;Pki|R7xH1_QliDi_Xp@xs_rw^ z1uDOQg~|sFva4%b#U`_0!yc|K9%3XjQZaPp(0nsb;yX$R>Y1$B(YP2^`k05Q2uZr) zX$e1#mg#B6PoD%iG6>r`qL~GZaK;Mk@niri>aSfOOWni_erh*N{*m0RLCHr{=z5fc z7v_x^ocuO_dH-ktzW_&MMO`cjvM;@^OtZDW*3jP>cP>5_*KwISf(PGif{JjDW(A7qzl;f-)Qo2wN?siv z;#ub3caoH8eQS(<$4pcy?+D&48D3cSm_R=Evh|!3F0q}4n>@CEMI!Yrj>Vw5vN8=Z zJ#S#sSUFG1R0~^^c}Ylm`?jm|tE!$vQ_$17h^2#yPj#L_`%E>Ren}n8CkY-=Y3V3_ zJ%O2$ae08Y-q?H!-;m`()`4+lfB92R$SitfcgwW zUTQ#dshj57sc;%?q9X51b=n>qVYevQ{$hs7>y6mKfe5h@8s`T!oausJF3)mVJYu|hC=X$Ptt$jToex7b_-T=V51ONa)0|2aS z007G*xtAjs1-X={{<{IH001Bv0EKOXPmc=Mh2<2^4#}E^rfu@qf9t{m{%)jWTrcJi zxO%=b*$n`oHqvr*wgtIhtdJUBc2_YLx})`%XD(8j(jgK3Rj)MEV-iF7l%*2;~902+4DBf045BujjX_c{03zs zlE=9%9`^w{ zW;mw}cumj9k$>dCzQ;&Qdz~LV4ZLX6B7KcV?==(8+s6-&MRf5Boxubvn#q)6IK0fk~*=g1jVu~I~2#$E#@nJEoMflmG zQ=WTB$6D29b?IPAMd=`+x|Gw5{(@6!L}4KjOLwU-C8@GBo>=XcfOY-OS3@su5^%s; zu!dUO9Rg%(3}v1l3Ub1PX3i6)`Zr;Jj|u=yZFJkw%Numh+tb#{#nuz)YUjqL;pSuQ z?8oQq=;~yPCf~fg$zj`I1k4DL!!~9<;|9LZflwq(MomQTx{U6C+4fZ;>fzvhLsA3?lb2D3AZ>i(rVXmiOp=LGK84w zzUNM1{N{oA4MF?e&AAZffkBDPcHq_nnaGl?CK0Ue&yAAnv+l^;GP$csKc;dwf18Y6 zS`D1HJkw}KI#MoO6CkcLWi#eizVfN-LGBCimQtx|L@LiI^*5<$t=LeO(4+#edAo zNE&lk1OK9?v8FB(L#0q{IPf^9(4;rkN=(ASs%ytCxz$Q4Ioz4|SN~YHvt4B_r*6*3 zxBm(9h27i~TWZYcMx9p_R{D;@X`W`+yqMlB9^cI)micr(|4xz#hom9wx(ni(#0W7j z*D?rc4Bg}VepyNY<~ha>Rd_0J@AS`zWDx(!yMokOXyQy{ggY|a;06D5Chw<`;fL&* z7$unj*vRwj7SRcTG0aI&@Zn;?9c?F`0;u*1m1s`*vyZDaTY_(IGinhSr+*#fHBHfP z2{GNjRAc*8@ztBJbov{}E2j5aH6Hq?DS98T2fIbsu}3v!kEuYL_NbUlz`-Iwx$sRm z$5P$ARq3NEg&~}#-*!KR!WkptQh!GKxnt%N{0?8+*q7!5ggryKm+#Yhc#8w$I60y| z#XM=PqUb;|)vMAB3n$^Lmf)DW+x z-hA^3Fx89FYwKmWyyUeN1Y2LPc_cq~>b(XqM32#<{Y~{VLe=FDt6QP)k&%&mP1E2P zwq{k7go?-1h7Q(=b&osj>3M#ujx{Z|86!~Ku8xPj=s8zCG)bE$NVW~JVw5vhZ@tgd zbc$`UqHTAzWZS{kPN&m+w|%=`wr?!wAv}cQKx8#xA~|CBpYe5I=1@G2BD4I$QyVO| z9UzL9`B~#sI`la#@3Zy?WCL!|Txob1vHO%gJstJ7tkJb5-9>hK>2&$)dU*Y}pf-y) zYZnJ6GD_(3lHCZuTP12A~8Vm=g=c3o8}dz-D9Vl`_EZVP|2$on>det&Bc8s#03)8 zaO)+Nwk8V@6s|-hrSd=uTs^mH`Z1-+K(KE3-H=-r>2r&paR@{Q$nwqof~?*4ZLmK! zzQJ^|B)usCF3#s{ZQRN=tR+f>vP{`eyxxo2)b327Ra%h;TP%as0rSNm&xT-?B^x(m z3YNJ-^jrhKk$7$aS2^eUqp%3pS>JuBlW=-u4P9~S#rh+H-xo|7-1N~{QG^)k)}{y8 zxnChUYuM2nZj9S5xj!Q>?aq12j9a!CH8|n&fLeFJSmr}%KdY?9bzBZ7Y@r`SPbZkE z^E@=ATf~!9Dsr#<2G%apQZ+X;p~pyVj6m>|;+W;M(sT?ZQ-Dnk{E;Z>4y2iP5TaYDo(|1nb<Y$)eC`(eB!OYD^6rF@)#nG(T^}2^Z`XAkYB@7Q^bboQF*RUqlS4Md zwIe8C)J1mzau8Y0?oun(Nq5AS8MMK*fa~H~2t9~v-S|Ef`K`GYcnm&$5$DpXasu2m z^EeytY1#>gys|@N>h#EMf6e`pSB7wfwfdYvS1)L5Bv7v1W9}PTPlY?J%7*1UU^D(S zd53!Iy^SYOVsM^md`PmoVsz9HO{HCBfM!W9|Dj>GPys+kcpU_f6FN|?l>@W_}Z0i=_N=21`HxJRH?y>shS6VA!|C!f~@hY?SZut|%5? zQNC{4$p*`P-R!Rynk}|a_$nkwREmAofnCK?PoiCFEoS?&I;jZB4^>8&If=zht4Cgo zc@+k(c7fV_pv?QkXwAYbC6Ygt+*+|+ubCZALdbTvSn%1j&X&2vGDs(DGEe6Oa29FL z2&`(C5?e3abK#ZEgQ(nq2W-qo^uPW_JV3YyrU|6VgBv?vVM066Coc^JpWl1@RyW>? zEBmT|3WO*UZ@<7azsbb!7ZC`H}9M*3aoXz;ks97firHLb-@dopy(hmP#Szh*eB$6rLL4WOhwY3?O|qb zJa{r(x-yh%;l+FSPbN)M6phl34^NByfK#9fKQ6LUM|Vvcx3N&j5Blv=Q+TB`Tl^1@ zf_Ck9t-Ne7Wii4?#U0`&dr~IPVYTTSixq25EC^J1N?j;XHoB2_=+ReTF!|<-eU+qM zAy|uIm}GIfv{c*thg6L5;t9ceBA&0jI)2mlV|+hjxa_?m#k#z>)1|vx9VYy{u8}Hd zDM(du4J2jm8UOo^h!w%Ppwni`@{F%{t1AI*nQcRboQqCgnWvAugFkU8^h`PTeC-|i z-jP)JSX;4zB{j_UL*4JMeM>*Iio6)1ZGWnM9>Wnf_ z8PrF{1N(CvmOqDva~vk4>^U_2ZyYASCqu=50J3DE+_??Izu5d`dor8+A2wtYRyJRm HSAhQl^Uh+7 literal 0 HcmV?d00001 diff --git a/Weights/TriggerSF/MakeTSFHistograms.py b/Weights/TriggerSF/MakeTSFHistograms.py new file mode 100644 index 0000000..b88f141 --- /dev/null +++ b/Weights/TriggerSF/MakeTSFHistograms.py @@ -0,0 +1,183 @@ +from ROOT import * +from array import array + +leadPhotonSF = [ +[ [0.0,0.0,0.0], [0.54,1.5,33.3333], 0.6810978055, [0.00519262468286,0.00519262468286] ], +[ [0.0,0.0,33.3333], [0.54,1.5,35.0], 0.762424051762, [0.00606810658022,0.00606810658022] ], +[ [0.0,0.0,35.0], [0.54,1.5,40.0], 0.778700113297, [0.00311642304626,0.00311642304626] ], +[ [0.0,0.0,40.0], [0.54,1.5,45.0], 0.795413970947, [0.00303555330709,0.00303555330709] ], +[ [0.0,0.0,45.0], [0.54,1.5,50.0], 0.808078825474, [0.00422688722282,0.00422688722282] ], +[ [0.0,0.0,50.0], [0.54,1.5,60.0], 0.819607198238, [0.0074540959624,0.0074540959624] ], +[ [0.0,0.0,60.0], [0.54,1.5,70.0], 0.851548194885, [0.0158679720146,0.0158679720146] ], +[ [0.0,0.0,70.0], [0.54,1.5,90.0], 0.868812680244, [0.0302805485333,0.0302805485333] ], +[ [0.0,0.0,90.0], [0.54,1.5,99999999], 0.84904062748, [0.0736341456908,0.0736341456908] ], +[ [0.54,0.0,0.0], [0.85,1.5,33.3333], 0.898902833462, [0.00100959242125,0.00100959242125] ], +[ [0.54,0.0,33.3333], [0.85,1.5,35.0], 0.958884656429, [0.000808069907818,0.000808069907818] ], +[ [0.54,0.0,35.0], [0.85,1.5,40.0], 0.965837776661, [0.000305692172361,0.000305692172361] ], +[ [0.54,0.0,40.0], [0.85,1.5,45.0], 0.970546543598, [0.00026519245806,0.00026519245806] ], +[ [0.54,0.0,45.0], [0.85,1.5,50.0], 0.973720788956, [0.000322496232078,0.000322496232078] ], +[ [0.54,0.0,50.0], [0.85,1.5,60.0], 0.976857423782, [0.00049804973141,0.00049804973141] ], +[ [0.54,0.0,60.0], [0.85,1.5,70.0], 0.980784952641, [0.000875908397042,0.000875908397042] ], +[ [0.54,0.0,70.0], [0.85,1.5,90.0], 0.985170722008, [0.00134766564404,0.00134766564404] ], +[ [0.54,0.0,90.0], [0.85,1.5,99999999], 0.989234626293, [0.00161217085916,0.00161217085916] ], +[ [0.85,0.0,0.0], [999,1.5,33.3333], 0.924440562725, [0.000377342171653,0.000377342171653] ], +[ [0.85,0.0,33.3333], [999,1.5,35.0], 0.971459209919, [0.000281199348832,0.000281199348832] ], +[ [0.85,0.0,35.0], [999,1.5,40.0], 0.979482769966, [0.000115534010648,0.000115534010648] ], +[ [0.85,0.0,40.0], [999,1.5,45.0], 0.982116281986, [9.05590478798e-05,9.05590478798e-05] ], +[ [0.85,0.0,45.0], [999,1.5,50.0], 0.98511326313, [0.000105097340197,0.000105097340197] ], +[ [0.85,0.0,50.0], [999,1.5,60.0], 0.984966993332, [0.00111781327877,0.00111781327877] ], +[ [0.85,0.0,60.0], [999,1.5,70.0], 0.986653208733, [0.000295146360226,0.000295146360226] ], +[ [0.85,0.0,70.0], [999,1.5,90.0], 0.98967808485, [0.000409470102893,0.000409470102893] ], +[ [0.85,0.0,90.0], [999,1.5,99999999], 0.991826117039, [0.000323844893609,0.000323844893609] ], +[ [0.0,1.5,0.0], [0.84,3.0,33.3333], 0.632706999779, [0.00549277181579,0.00549277181579] ], +[ [0.0,1.5,33.3333], [0.84,3.0,35.0], 0.705228030682, [0.00671890259184,0.00671890259184] ], +[ [0.0,1.5,35.0], [0.84,3.0,40.0], 0.729334890842, [0.00316854098441,0.00316854098441] ], +[ [0.0,1.5,40.0], [0.84,3.0,45.0], 0.752555906773, [0.00315251887441,0.00315251887441] ], +[ [0.0,1.5,45.0], [0.84,3.0,50.0], 0.770330071449, [0.00465503941806,0.00465503941806] ], +[ [0.0,1.5,50.0], [0.84,3.0,60.0], 0.786793231964, [0.00682801354333,0.00682801354333] ], +[ [0.0,1.5,60.0], [0.84,3.0,70.0], 0.816314458847, [0.0133762782917,0.0133762782917] ], +[ [0.0,1.5,70.0], [0.84,3.0,90.0], 0.857941806316, [0.0154562893515,0.0154562893515] ], +[ [0.,1.5,90.0], [0.84,3.0,99999999], 0.886136233807, [0.0226681113534,0.0226681113534] ], +[ [0.84,1.5,0.0], [0.9,3.0,33.3333], 0.872963190079, [0.00255836970587,0.00255836970587] ], +[ [0.84,1.5,33.3333], [0.9,3.0,35.0], 0.971260070801, [0.00169839848133,0.00169839848133] ], +[ [0.84,1.5,35.0], [0.9,3.0,40.0], 0.98152244091, [0.000588619026864,0.000588619026864] ], +[ [0.84,1.5,40.0], [0.9,3.0,45.0], 0.986895382404, [0.000503406678232,0.000503406678232] ], +[ [0.84,1.5,45.0], [0.9,3.0,50.0], 0.98783403635, [0.000744728019948,0.000744728019948] ], +[ [0.84,1.5,50.0], [0.9,3.0,60.0], 0.988175034523, [0.00106780215648,0.00106780215648] ], +[ [0.84,1.5,60.0], [0.9,3.0,70.0], 0.990262210369, [0.0019675187127,0.0019675187127] ], +[ [0.84,1.5,70.0], [0.9,3.0,90.0], 0.990519165993, [0.00309011400271,0.00309011400271] ], +[ [0.84,1.5,90.0], [0.9,3.0,99999999], 0.997375607491, [0.00160821560062,0.00160821560062] ], +[ [0.9,1.5,0.0], [999,3.0,33.3333], 0.827793955803, [0.000860672672767,0.000860672672767] ], +[ [0.9,1.5,33.3333], [999,3.0,35.0], 0.966853737831, [0.00050528634985,0.00050528634985] ], +[ [0.9,1.5,35.0], [999,3.0,40.0], 0.980059444904, [0.000211483705081,0.000211483705081] ], +[ [0.9,1.5,40.0], [999,3.0,45.0], 0.98436075449, [0.000193397430742,0.000193397430742] ], +[ [0.9,1.5,45.0], [999,3.0,50.0], 0.985809862614, [0.000283829305427,0.000283829305427] ], +[ [0.9,1.5,50.0], [999,3.0,60.0], 0.985918521881, [0.0011583820968,0.0011583820968] ], +[ [0.9,1.5,60.0], [999,3.0,70.0], 0.990703165531, [0.000662821684507,0.000662821684507] ], +[ [0.9,1.5,70.0], [999,3.0,90.0], 0.994281888008, [0.000630687802869,0.000630687802869] ], +[ [0.9,1.5,90.0], [999,3.0,99999999], 0.997075855732, [0.000553451337113,0.000553451337113] ], +[ [0.0,3.,0], [999,999,999999], 1.0, [1.,1.] ] +] + +subleadingPhotonSF = [ +[ [0.0,0.0,0.0], [0.54,1.5,22.5], 0.735150516033, [0.0124989687013,0.0124989687013] ], +[ [0.0,0.0,22.5], [0.54,1.5,25.0], 0.741944074631, [0.0104878349021,0.0104878349021] ], +[ [0.0,0.0,25.0], [0.54,1.5,27.5], 0.759739041328, [0.0085182165305,0.0085182165305] ], +[ [0.0,0.0,27.5], [0.54,1.5,30.0], 0.768184542656, [0.00742963937731,0.00742963937731] ], +[ [0.0,0.0,30.0], [0.54,1.5,32.5], 0.781085252762, [0.00607655931058,0.00607655931058] ], +[ [0.0,0.0,32.5], [0.54,1.5,35.0], 0.787775456905, [0.00543688578024,0.00543688578024] ], +[ [0.0,0.0,35.0], [0.54,1.5,40.0], 0.804114878178, [0.00298680600149,0.00298680600149] ], +[ [0.0,0.0,40.0], [0.54,1.5,45.0], 0.819701194763, [0.00290229435272,0.00290229435272] ], +[ [0.0,0.0,45.0], [0.54,1.5,50.0], 0.832450091839, [0.00402395578062,0.00402395578062] ], +[ [0.0,0.0,50.0], [0.54,1.5,60.0], 0.840809226036, [0.00712945852573,0.00712945852573] ], +[ [0.0,0.0,60.0], [0.54,1.5,70.0], 0.863569915295, [0.0158928774505,0.0158928774505] ], +[ [0.0,0.0,70.0], [0.54,1.5,90.0], 0.873844146729, [0.0302822319315,0.0302822319315] ], +[ [0.0,0.0,90.0], [0.54,1.5,99999999], 0.880566418171, [0.0771217355214,0.0771217355214] ], +[ [0.54,0.0,0.0], [0.85,1.5,22.5], 0.986665070057, [0.00290137137345,0.00290137137345] ], +[ [0.54,0.0,22.5], [0.85,1.5,25.0], 0.991021811962, [0.000881010903373,0.000881010903373] ], +[ [0.54,0.0,25.0], [0.85,1.5,27.5], 0.993968784809, [0.000365351913748,0.000365351913748] ], +[ [0.54,0.0,27.5], [0.85,1.5,30.0], 0.992568671703, [0.00080863027601,0.00080863027601] ], +[ [0.54,0.0,30.0], [0.85,1.5,32.5], 0.992429018021, [0.000927701357281,0.000927701357281] ], +[ [0.54,0.0,32.5], [0.85,1.5,35.0], 0.993385374546, [0.000273290379079,0.000273290379079] ], +[ [0.54,0.0,35.0], [0.85,1.5,40.0], 0.993640601635, [0.000134285419182,0.000134285419182] ], +[ [0.54,0.0,40.0], [0.85,1.5,45.0], 0.994162678719, [0.000117730435885,0.000117730435885] ], +[ [0.54,0.0,45.0], [0.85,1.5,50.0], 0.993752419949, [0.000159345143113,0.000159345143113] ], +[ [0.54,0.0,50.0], [0.85,1.5,60.0], 0.993596553802, [0.000408447927548,0.000408447927548] ], +[ [0.54,0.0,60.0], [0.85,1.5,70.0], 0.993239223957, [0.000583930534644,0.000583930534644] ], +[ [0.54,0.0,70.0], [0.85,1.5,90.0], 0.992464900017, [0.00420925695041,0.00420925695041] ], +[ [0.54,0.0,90.0], [0.85,1.5,99999999], 0.992560684681, [0.001409051152,0.001409051152] ], +[ [0.85,0.0,0.0], [999,1.5,22.5], 0.99289894104, [0.000459394846531,0.000459394846531] ], +[ [0.85,0.0,22.5], [999,1.5,25.0], 0.994199693203, [0.000459394846531,0.000459394846531] ], +[ [0.85,0.0,25.0], [999,1.5,27.5], 0.996753513813, [0.00019635266302,0.00019635266302] ], +[ [0.85,0.0,27.5], [999,1.5,30.0], 0.997437953949, [0.000105276812567,0.000105276812567] ], +[ [0.85,0.0,30.0], [999,1.5,32.5], 0.998195171356, [0.000120250429564,0.000120250429564] ], +[ [0.85,0.0,32.5], [999,1.5,35.0], 0.998485267162, [5.53143156827e-05,5.53143156827e-05] ], +[ [0.85,0.0,35.0], [999,1.5,40.0], 0.999016582966, [2.49504923609e-05,2.49504923609e-05] ], +[ [0.85,0.0,40.0], [999,1.5,45.0], 0.999399483204, [1.67820877013e-05,1.67820877013e-05] ], +[ [0.85,0.0,45.0], [999,1.5,50.0], 0.999610304832, [1.71324623735e-05,1.71324623735e-05] ], +[ [0.85,0.0,50.0], [999,1.5,60.0], 0.999711513519, [2.13249403322e-05,2.13249403322e-05] ], +[ [0.85,0.0,60.0], [999,1.5,70.0], 0.99971395731, [4.15814070615e-05,4.15814070615e-05] ], +[ [0.85,0.0,70.0], [999,1.5,90.0], 0.999810993671, [4.80188315009e-05,4.80188315009e-05] ], +[ [0.85,0.0,90.0], [999,1.5,99999999], 0.999862849712, [4.13267813275e-05,4.13267813275e-05] ], +[ [0.0,1.5,0.0], [0.84,3.0,22.5], 0.660493433475, [0.0189502448237,0.0189502448237] ], +[ [0.0,1.5,22.5], [0.84,3.0,25.0], 0.675408184528, [0.0141567263708,0.0141567263708] ], +[ [0.0,1.5,25.0], [0.84,3.0,27.5], 0.694644391537, [0.0109544486394,0.0109544486394] ], +[ [0.0,1.5,27.5], [0.84,3.0,30.0], 0.705103397369, [0.00740787194485,0.00740787194485] ], +[ [0.0,1.5,30.0], [0.84,3.0,32.5], 0.712502062321, [0.00673718942098,0.00673718942098] ], +[ [0.0,1.5,32.5], [0.84,3.0,35.0], 0.719795048237, [0.00535038292038,0.00535038292038] ], +[ [0.0,1.5,35.0], [0.84,3.0,40.0], 0.736499130726, [0.00315590127157,0.00315590127157] ], +[ [0.0,1.5,40.0], [0.84,3.0,45.0], 0.757362127304, [0.00315107032637,0.00315107032637] ], +[ [0.0,1.5,45.0], [0.84,3.0,50.0], 0.77392333746, [0.00466395309879,0.00466395309879] ], +[ [0.0,1.5,50.0], [0.84,3.0,60.0], 0.78893083334, [0.00687048443179,0.00687048443179] ], +[ [0.0,1.5,60.0], [0.84,3.0,70.0], 0.818382263184, [0.0135113622786,0.0135113622786] ], +[ [0.0,1.5,70.0], [0.84,3.0,90.0], 0.841334581375, [0.0180361620001,0.0180361620001] ], +[ [0.0,1.5,90.0], [0.84,3.0,99999999], 0.868076682091, [0.0243523219982,0.0243523219982] ], +[ [0.84,1.5,0.0], [0.9,3.0,22.5], 0.976537823677, [0.00489157641841,0.00489157641841] ], +[ [0.84,1.5,22.5], [0.9,3.0,25.0], 0.987386524677, [0.00432063989045,0.00432063989045] ], +[ [0.84,1.5,25.0], [0.9,3.0,27.5], 0.986037313938, [0.00208121288909,0.00208121288909] ], +[ [0.84,1.5,27.5], [0.9,3.0,30.0], 0.989064216614, [0.00131255284358,0.00131255284358] ], +[ [0.84,1.5,30.0], [0.9,3.0,32.5], 0.990115821362, [0.000987952486947,0.000987952486947] ], +[ [0.84,1.5,32.5], [0.9,3.0,35.0], 0.991321563721, [0.00064998825176,0.00064998825176] ], +[ [0.84,1.5,35.0], [0.9,3.0,40.0], 0.992424309254, [0.000392538029826,0.000392538029826] ], +[ [0.84,1.5,40.0], [0.9,3.0,45.0], 0.993506371975, [0.000357896360151,0.000357896360151] ], +[ [0.84,1.5,45.0], [0.9,3.0,50.0], 0.992881655693, [0.000667513509479,0.000667513509479] ], +[ [0.84,1.5,50.0], [0.9,3.0,60.0], 0.99330842495, [0.000823890634755,0.000823890634755] ], +[ [0.84,1.5,60.0], [0.9,3.0,70.0], 0.991472125053, [0.00186896480055,0.00186896480055] ], +[ [0.84,1.5,70.0], [0.9,3.0,90.0], 0.999999046326, [0.000123650174015,0.000123650174015] ], +[ [0.84,1.5,90.0], [0.9,3.0,99999999], 0.992282807827, [0.00316821950182,0.00316821950182] ], +[ [0.9,1.5,0.0], [999,3.0,22.5], 0.99313390255, [0.000574329725221,0.000574329725221] ], +[ [0.9,1.5,22.5], [999,3.0,25.0], 0.998826920986, [0.000196940017181,0.000196940017181] ], +[ [0.9,1.5,25.0], [999,3.0,27.5], 0.998663425446, [0.000196940017181,0.000196940017181] ], +[ [0.9,1.5,27.5], [999,3.0,30.0], 0.998753978656, [0.000196940017181,0.000196940017181] ], +[ [0.9,1.5,30.0], [999,3.0,32.5], 0.998868584633, [9.19086602694e-05,9.19086602694e-05] ], +[ [0.9,1.5,32.5], [999,3.0,35.0], 0.999282956123, [6.61051993856e-05,6.61051993856e-05] ], +[ [0.9,1.5,35.0], [999,3.0,40.0], 0.999356508255, [4.28383970185e-05,4.28383970185e-05] ], +[ [0.9,1.5,40.0], [999,3.0,45.0], 0.999553442001, [3.74810811141e-05,3.74810811141e-05] ], +[ [0.9,1.5,45.0], [999,3.0,50.0], 0.999633073807, [4.60404995528e-05,4.60404995528e-05] ], +[ [0.9,1.5,50.0], [999,3.0,60.0], 0.999295592308, [8.84144904368e-05,8.84144904368e-05] ], +[ [0.9,1.5,60.0], [999,3.0,70.0], 0.999351680279, [0.00015872505313,0.00015872505313] ], +[ [0.9,1.5,70.0], [999,3.0,90.0], 0.999379456043, [0.000199756405112,0.000199756405112] ], +[ [0.9,1.5,90.0], [999,3.0,99999999], 0.99986666441, [0.000145852311423,0.000145852311423] ], +[ [0.0,3.,0], [999,999,999999] , 1.0, [1.0,1.0] ] +] + +def GetSF(pt, r9, eta, arr, db=0): + for aa in arr: + if r9 >= aa[0][0] and r9 <= aa[1][0]: + if db: print 'Found r9:', r9, aa[0][0], aa[1][0] + if abs(eta) >= aa[0][1] and abs(eta) <= aa[1][1]: + if db: print 'Found eta:', eta, aa[0][1], aa[1][1], aa[0][2], aa[1][2] + if pt >= aa[0][2] and pt <= aa[1][2]: + if db: print 'Found pt:', pt, aa[0][2], aa[1][2] + return [aa[2], aa[3][0]] + return [1, 1] + +tfile = TFile("TriggerSFs.root", "RECREATE") + +r9Bins = [0, 0.54, 0.84, 0.9, 999] +etaBins = [0, 1.5, 3.0, 999] +ptBins = [0, 22.5, 25, 27.5, 30.0, 32.5, 35, 40, 45, 50, 60, 70, 90, 99999999] +lSFs = TH3F("leadingPhotonTSF", "leadingPhotonTSF", len(r9Bins)-1, array('d', r9Bins), len(etaBins)-1, array('d', etaBins), len(ptBins)-1, array('d', ptBins)) +sSFs = TH3F("subleadingPhotonTSF", "subleadingPhotonTSF", len(r9Bins)-1, array('d', r9Bins), len(etaBins)-1, array('d', etaBins), len(ptBins)-1, array('d', ptBins)) + +xbins = lSFs.GetNbinsX() +ybins = lSFs.GetNbinsY() +zbins = lSFs.GetNbinsZ() + +for xx in range(1, xbins+1): + mR9 = lSFs.GetXaxis().GetBinCenter(xx) + print mR9 + for yy in range(1, ybins+1): + mEta = lSFs.GetYaxis().GetBinCenter(yy) + for zz in range(1, zbins+1): + mPt = lSFs.GetZaxis().GetBinCenter(zz) + lSf = GetSF(mPt, mR9, mEta, leadPhotonSF) + sSf = GetSF(mPt, mR9, mEta, subleadingPhotonSF) + lSFs.SetBinContent(xx, yy, zz, lSf[0]) + lSFs.SetBinError(xx, yy, zz, lSf[1]) + sSFs.SetBinContent(xx, yy, zz, sSf[0]) + sSFs.SetBinError(xx, yy, zz, sSf[1]) + +tfile.cd() +lSFs.Write() +sSFs.Write() diff --git a/Weights/TriggerSF/TriggerSFs.root b/Weights/TriggerSF/TriggerSFs.root new file mode 100644 index 0000000000000000000000000000000000000000..734e1fe39f90797abaf14c788a466e6674326415 GIT binary patch literal 6587 zcmb_>1x#Gg_vTRC9frXNDDGaoSdqa7r?^XTx50`RDHJd6E-g^py|}w;DNx*(cC)+5 z{(s4Cvip*g_i|6ZbHDfUa?bhgadvQU1ptm%0RVs*06^0GVgWDDLN5mQVs8=t`vZnF z0D$2SfWsTG11bvQ?j;kC=r^{soxFJeuP%SU-#`$ZxKU&P=fB7=y8!@*>N3`L=Gf}a z)|Qs$&QKW_Hs`;5{>uUY@h|cp8USd00ls~)e-QuxVE6xt79!E#-GLnc@Dl%a-~5Y5 zP=CiQgJWlIY-Vk5sbb~e>R_)9{jUqUiYOxp01hC)6Yf~sMLuaxS=%F2m-T9_ln@X(4GhEI>Y5;btUx?=(S^oj*Y)J(I)ZuQab3oVcuJkRREUSH9u`Z}ntZ<-~p6 zUx(}WmMR99_-KRJ*0$>r-QIs;Ft%oSz9li+){vM*{e=ArhaQQ5UQW8vY1&CEd`a~^o+!|t)=)xq z{*?fq&q|&{TQax}cy~24Ms#DeB01qtN-%CP8CwPS_3r(;Ijf7JjO^26ZqYD78teA* zlFp&C-WiAwy5oDjpdPn9y2!3O&)j2YGEe05d6g3nxAGR^WU?ok-2E-_oXp?DAhb&t z^58P%$DhirLT5sszQ7FDh@Pqw#_l6sJNQca^V^zNGi)A-ffSrdy>o|sy<>#eXkY403 zto}-tzN{Bl^=m!Nhstg3`v{`UFk5Fo)~+vUdXPQlX#&OQPH{5&x(dQJhE@l$Tlm- z;K|AVaG$N}AAQr2Rzq+7daO0SMvl*o;zfn}tS&LmFQkwu^>AdI1rwmA-J%=t$tJ5O zKc|yvKo(7i`gGN@33*4d^;AeBMUj_6Ix|MJca)~|3+6Lw`?f15znozE4rFrSFRrao zHL*=;ldbNPw#UoYq~%14p#s`Lj&xb*xp}~)zZ6#_*+Po6?G$z>Yf`2(7B~2V$k3pK zw{^}f+W^n0VamZ!B;vWM>4pyr52{U3QmGht-M1Hus(@@#iYrpg_0{Xem(YF>BFHSs z{&Vuzx5fl-QYuHRvctCLr8fgfcSkPJ5-Pr;dx=R8QG(;VyfRRswE64%>09$k-S{bc zRT~mL)?U+zA70|Y-$|-b_BZI%q;fEbS;CJd6__PEL%ZhyevZCzYtcZB?KeHmfi<`;+o*4hr+DmSHLdBq%F7>G?p1SKGh3*_toY#8A5GN7{bl2JyW3K;33f=sgmQu!i7 za*nvSKVf4~GDKFc)cK5QoHW74S-1tJ%-bGEPlVQPMcmni9Ib~;xae{i>#iB$(?G6Y zK;qCuk%In5$K;S)jGuX?$ucb1Sg29&5AZ>R_=G=<+V+lkerHe@%(3HS?aV$pAq`^M zn@e4GgJ_}41zK{@X^TVF-!o)o86jmVNKTzj#uEheu{jrl>1yAtw;o!%24;AyG88P$ zUn>R7EyrMb1i6&3zj1k^1ByeO_8APLogNEGKSW#>4X!FpZs{pCNz0uyoKp8yPVdSq zu1rQl{q&Z&-LiN?75#x3qMB?lrTA7kos(K;rJhkOt28)JDzAvHEa_XGZqC7wsMC2J{&Bm)?aEozHv_q zZI`SFAR@Z=5enVd^AP${-WH0sdUZ~^ah5nDr~RJNF5kR;evlj`?|x1DK89^Xl6ud* zA+2?(cw>%>KGpEo$uV#QQ}Y7`xgcu6GR9Y@kq0z1LG3*oEhG6Cbh!z{+!zTGG9j4G z-2H+5*vE92_jj>?->n(#&~b5PHKlGuODR$k^WrT2Y|7k?l$9_13Zih{0$SC91rB z>AbLE&EFKSv6wK)QH4`uxRfD9{PG!PE)hrVb%wGgGPdR-X`hlRui>xx3>^*L{sn(! z{*f@DL**E`QjBhz3NkbPn&x*E**>qlJK1LfbLcPLlGB6^FN)KiYhO9qZK=&6Z|zJD zRbk;1&Bs25-PRxN|OI*bk7yEbO`l5GUp!XB+!2-b`9 zM>+X5oXn&gXgdHE+zdy8!0qyZD+V7pfklDaOGQ(eNm+nHkZD10QIlM~Tu@d3YC0&mjTV zQSWl?avycSa6e(vVbt=&--d0!<#%ppLqL4J3-1a?$U__BU_iV&_vd?GD5AdBIU8GH zAc-*U{C)F_@2Rslpi}e&pjo3b?xw*nkQ;fw8jg4sk03el9Fah8#PJw7C4!(Y>n7K> zx({A)APHVF(r;mcKgdxJO6q9wgA2?*Mr01(evuEtizVcWU|a3>q4x_3w@)H6oJh z@v!|$)HO;MiKj0-MwiLQ?A+^`1Jw_|Tl^lq8GT%NX}L3fMpFl;B(v0#xI%KzPjMJgv_~Dh>oMV32d&@x)osH58Z}WUHA=sQ-wJ5P)aU&zm zzRB3S4U`BVn)_l%iSm@%vLGZ>=Lq7j^#a$)xDsAUck1D^+$VTGBBpjdgsV4M9wcp% z7(=snmfN2K{kPPRhm2%e>20KlRg}^F@;H$!`dq-g@))Jp?mUlhhk-Jq?oM-FU}Q|s zX?_!H=Vvo577?$5V9K29#%biWEPAF{aErO|fsXOPswa|Xeia9P=#)zXWb8l&c;{XN zN`o+hOYQq{Zd7FlwAH@zRF^&?xux=Cu%pgT(9Q?EN^i5!*)>tFwMZMATUciuVYo+iO#T5zMx)#t)vNJ})-ieUs1$mNRVj3kOSR80} zv!IXs*$eW5pmI`Sxz;mc8D>DGn9R%>5z`>wBV_w@8M(i4Z7=@4NWvR8oof&1F| zF+VnI5P03PZ47EVOZVz98VDZDCmp}P<`RkzS?fJzU)*uLd_wUHhQT=3tKn-I_3=)v zeh4QbP*y87pE8UXUA!;vv~Ue;T8!5um>dowSO^T3%Wt#tp|q|2!JVV;i5D28i_ZWB z`q92L6G)ssTdzb+WWDcS= zz(|C1T3c#_A(wYyq?Jwxs5dDg(RCUD2pTRDfW@@8p_MOG= zW8J=mu(DaPMuQZnGWSMcr*Nlyq0~I}DX8@WR>_)o0Qy_?@H~9_%~k@`cX0Cs2*_3O zN&Vs>!KQLpcZADN=nTQh&O9|B-AHqGsa~*&$RaA{%~UZQ+~5yu~KS7M|2*dK?dHoR2JjwOQ8qF$)!5k@=9cq%wo~ zi;x)i{dS8MK_z1&JG%W4&CuKL-9`=*SZ)Jdj8R-#@P6Ir!vN9)w7AsZvO&4(5p!dE z!_PkiR%r-FINdDAGc`O5fF0=OWTRGi$4!H=!#y7-kd&JJLUmGF> zG+9SHZeei}b+PtGbG&+rP-=I?aO$GLIO`PZH@d=}249_6b&7_y=bPwNE5l@cG?&}? zYS~*hC<{%V!Pbr5sss!>EvD&jv$~$b36bYkqOv8aa3x+9D+qGPo=k&d+E~d}NXqbL z-NBCZgbAjLsdU1fZ|*+HBptqqj1JhhiXSH#6n<7<;(>pj@q7vu8&&`03t7aiRh z@w;w~3nD6W+xv2t3_Z{)>h$gBH;x|aeAd#`-!LZXW=(qCufcK!={tl_GX15nesrq! z>MRdQLkIeW_R?$vYF0zoM56o@XJXBds;0J>avgGvsuX)Hg;(jcnl@R|YSF2n^pEWW zDsg!2{N0Chg7Uv#Co-0K(Ivy+XN9@>EM~m6iY0m-I!rK6dM9cbRC)ED7`YR^j|P>C zxb@j%?bN1yg6LA{h;e=`S^kja$pF27g`Z-`+Z3S=ZjOqjl>{}vYb`_VKna1<_OqBHAPtK8G>XH!|x z0z1PR_?Ud`6JX^{Cao}i--`8Pb)mUvy4I<1Q5od_%#aierVc!&TxGHks)tMJ|aYiFNC0c?bjTJOX1dQr&9y@ z>0Bw03j5HcQmg%LfxfPgZ4>jxVG0eQB=wi-LalXn=loOO#_c|8RmG~!1@-$?|o=*QU~v-}EmMQy68+@=+4__mKf18@np zf_1k6&)SNckM182;e|oK-)Prtb zgHAif}6;DAAFEGfyqH zrMPe^rHspvukx;|ma~$aw%XU1dYr_I&Uh}L7w5?@mGwNF4<^c`B}}0F#&dI;p}yR| z#u&cSlhl-l2>gu^qtr^JoO-l6-X`KCvdbXtgdqMn?ja;Mc*!0-c%dw_qEa z1%FM)zvZC?&Df$1g*@nlB{`5jJ`p2=Eh}GiM(e-vLi9iN@o&8FG8FiyKK|o);l=gp q<+SVnZ7T5p=>YA8i_kyvX#b1+zX{BXbo$TaU$hlTvj=Y?;J*P$XCvkS literal 0 HcmV?d00001 From 84825898ea6c8ce63361d50d239bb3a58e366302 Mon Sep 17 00:00:00 2001 From: Andrey Pozdnyakov Date: Thu, 23 Aug 2018 16:56:31 +0200 Subject: [PATCH 2/2] Add missing template --- templates/models_2D_higgs_mjj70_vars.rs | 4 ++++ 1 file changed, 4 insertions(+) create mode 100644 templates/models_2D_higgs_mjj70_vars.rs diff --git a/templates/models_2D_higgs_mjj70_vars.rs b/templates/models_2D_higgs_mjj70_vars.rs new file mode 100644 index 0000000..c395398 --- /dev/null +++ b/templates/models_2D_higgs_mjj70_vars.rs @@ -0,0 +1,4 @@ +mgg[100,180]; +mjj[70,190]; +mtot[200,1600]; +ttHTagger[-1,1];