From 5ec7e8318c98a69c38583846ece628aa9600e2c4 Mon Sep 17 00:00:00 2001 From: Michal Kreps Date: Tue, 19 May 2026 13:36:24 +0100 Subject: [PATCH] Add FONLL distributions for sqrt{s}=13.6 TeV together with information how these are obtained. --- README.md | 2 +- rootfiles/fonll/LHCb13.6.root | Bin 0 -> 13582 bytes rootfiles/fonll/LHCc13.6.root | Bin 0 -> 13715 bytes rootfiles/fonll/README.md | 3 +++ utils/makeFONLLHistograms.py | 39 ++++++++++++++++++++++++++++++++++ 5 files changed, 43 insertions(+), 1 deletion(-) create mode 100644 rootfiles/fonll/LHCb13.6.root create mode 100644 rootfiles/fonll/LHCc13.6.root create mode 100644 rootfiles/fonll/README.md create mode 100644 utils/makeFONLLHistograms.py diff --git a/README.md b/README.md index f003872..d07a5ab 100644 --- a/README.md +++ b/README.md @@ -101,7 +101,7 @@ Particle settings should be defined after the corresponding `@#` tag using the s * `energy`: * Sets the pp collision energy (in TeV) used to get the correct parent kinematics - * Supported: `7`, `8`, `13` and `14` + * Supported: `7`, `8`, `13`, `13.6` and `14` * more types may be added in $RAPIDSIM_ROOT/rootfiles/fonll or $RAPIDSIM_CONFIG/rootfiles/fonll * `parent`: diff --git a/rootfiles/fonll/LHCb13.6.root b/rootfiles/fonll/LHCb13.6.root new file mode 100644 index 0000000000000000000000000000000000000000..b88c0bc54cadda642f4b68ef837bc9db22b4590a GIT binary patch literal 13582 zcma)jbx<5n)Gbc%U_lo`&;$Yjf?I$92^xaCTX0!maVNM#a1Rok5MUuJi!Cm}b#YzX zSr%UMdsXlK`F&H>(_Jlh&h$)8-E;35H)m%LG_+e`G&D3zG&BKXG_RY2=dpHna&X{L{vd14%g-&y{ciy6f1ZJc{m=G)_R-Md{+j-K zH~*Et6^GLQ1e_%8A3pZ{|JJko!-xJ4uO*h|2j2IXRvu<(`pQpiNzn4y(Y(>)vhX}q zuAlWK^dxdb)7tN9X(*ExeWi9&{x}l%eEVzI*RWVJ7FYMB_N7RPIu9_i>GmSI8F#Q0pV|TZ{$Dfejva#P|h>o9j|43aUH>cK$ z{X8Qz-|p9LG1*6xJ_$P^YtFcAtGZk~>99q;a`rY@`Lc|NXi8hL1(C&_%{;21K=|Dp zDI(nNpZB{{-tsBZ$yc;BOis7132n8NWDFnpHCh@sPjz;3KYm#Ja_FkZzVOL`-;zMR z>+bkbmLU+j7YJEJJp!feAx_IRr?VS0XEy|+-U#f_jK?iJ0s<=3ZP?=@b~TkhljXw< z3_ilVvOX&x73pW$e0oBz$ve2)D=~JT(Oq?{@3UdzaBC-M=(@H>(x{$0y&59j)u&|> zdcNlTd6~^KMFtdOJ?V~aj--TM5%huXJ0w{itX=6krBFUTB76tfLM=Sn? z{1-tSrjYYf$=Q7h-LQk7&05de{0sBzWQn$9Zp3DKTmc#D$Dy=~TP+1~f@Wnky{kLw zU8_k{o%St|(Z*9RBifse?d!5d3(-=8!dxSF&Ftjq_}k&U9TUncT3w=h>e)S1)$~X;>i5p-&B~)B zi>Iw2nzy(H^-LM{(pB5ZtsnqeBy5)NQyWB}{RF8!kJj-51x4)qn7W>XQ_UVfCL6`m zxlqAKa=y36p)fP5(@~+78v!_RF5SjH6z!(u9pj=5!@4BTOrAwZO+1PHij&b!Y`pS< zp56mL2lX^Uk?djuB7;Agt=80Z6{Blg$L_1N0gAyBg22 zw##M%ZM{(~c3unUdml&o_|0E9z`<91N7o;&zBM@TP|ls*`8J3B%d+8|%_#~t4dA?L zVOm)=Y(I9*Px}fm+Q%wC6cjCffruLm(OGx`Yu?4`3&5-C@nr}kn~ekET+{Jye?ug^ zRK`d!BI}!pMvp|HamQjMtoxDg-^56mVJ%FC$+bL$F}n3g)vo~F(4rj8CvfJZ1gQE1 z6xw34&zKP;;E-l2Ja%k2#^Ktksvpc7er67fEdJeN|M*oZ>Rj0=>|9xOIC+B&lCYEQ zM*^29WItk;MyhnJak?3r=cbz~i{EhPajHX(Q=j@pzT}Y~dUy#J4%2Y#T`fVbZr^pW zU|wfWwMZcXilI~Q>k2Shin(vdYw`^r-0M>o{czPpK#mfWxjnDqDD+m z#WLb`2+}NxlZ)RrylFz1zd@T-14JB@Z5J$muUUvd;a-IbO&tCCxmDQ5TCG@mL8qEQjaeD?2zC8*Km{doiFPN<1tf?H6czm03DYE-KHmGXafHIP!{>k zlUeqjW?_ZT=e9HN?1gU7RdPM%#wTm2_ow1e1`^zD_}UObuSW-42_FaH(9rqSYLPC4k8?<&ZSB0qy%}VJaw&ReR5EgB zw8G`{y-2VQl=NmrgxK#mM_{Zouw1~qE!#yK6>M5tZ`!5lG9-S&Qe%JJ2VbC^p;6adV2sv*_B$#J)Ho!_soaSo&VAb(wqapm(d}ro%qE3*1 z7D9wbG@M(DG=(jgsvV$>dUiv0IS@^dphdCnK$kb(J|g#!HD8gzn`pRtH;&*T8qU^` zYrLd$FtZc}c?VkKSdjZ85qC94hu5?d>J1=29UD@Go7Z>kX&5-T;P_3o)};5iY8z8f zoFs6Sdrm^TlK~MV5kQZaa;nRX)a>6M`6e3)k-We#Rx?8D> zS9Q)b=Z8gY8B{2@Iv!v;>};x-YUcQzvzbd+Xqx*fy9w%b@?eCr}h{ z#d92V{e@=6di!=WDp$vm;HHZ4{-GGU>VFa`iN?GTte(;t*4!D}8_|>r8kJrBu?q*q zSebp+x=pV(c{_GNktuOITVm^g;13whQv1>8ODkoW%zVHE@0HHzGx}= z`4@bfEp-bw%>#*NNPbROAk1j4?FV$mb(>$mPU4Rx6xh51r7z83m_Dmpo2WYS`Svb*_kGn42XR^ z)3~&g6alqPwV9e=R6;WHPJRZ%S2FS7Z!C+%|#z%HUt=r)Y$)vw2UJER?l$RiA&Dc$j?Xv8qqc_@&#`NIv+G7z!`K zn0roGwTHnuELb(-;d@B2z5Jkj57KFe@^7t5X)|0sDAp-ecCKgwuE6?r89C@CQY4TX zld}e?5UF|MzANm=Su|uWj_geEA>TX5k?C^uY6|S#@s`YN_0`{}TBB6jNp2DpNfjcN zJzvfE-Db(44fD#tR;uRKma(RtKv}=o&O)03&v$NxGP;bHXt~<0nghPYhaA>XnKQU265MsG{ZTtR_GZY-m|X*c7~6SP z*pIWvpgZgCXRWskH^w)4H!C*`x5l@5w=1`Ri#Fj)F_GKl8^Ep6ZSL*zJnMV94mYv?|=?8om)5F9s4N-f}fu7w&b9*`S6A#Jkg93lqB;ZVC^$rfhJ) zwBj4-y+nxJ4Y!3yT(dV?;t%L|mg%oWNL-9p#i5|k^AF`JcQM)Lggbyy*_;evHP@AO z@jg)8IsJ~rXlM?Vh??8VruZNz?VM}JXtX;gLqshXR)hAy&%n;%1?e&J(%&<#1+?8O1 zx=$e1o%;F>vAJQ7aFi?g21lP(Y$LT7FR_JTpRme?U`C%ScONkJocfxVc(Vm2oaD;9 z@u4s6<*p8QUv8`{O$MK;)zYA__J)30Uu#}VhU=RR^S*)DY#Igr8ta^qmTXty4Y$6% z*hLx`f6c6LzAGeD6b*A zpYFey`dFU*NujWY%)U|8A~zCD0HogMo3hhvuKy8lQC1`Sw~GY$d|zfNs+>+;p{VA4 zZWNdt$T3}NU7;!v1EvIW?dw;%Z!}G%o#N4SmE|~5dZ5U@N9Ceb zMSO;wQGzHVP-;JD3R{akrvkjDy#IBIvV#A!mnOh1Vn3N7HyO+dRNl|3bXV0`o-7Xi zOjW8%OUeodf%mH>wRB2Hvns$pY#{p$nU#}19tP8{6$jRzc1n+w0J=c4{T_zpZ=wzH zjg#g48eZA}Q6r#Tf@<|ahR6QgR1@dxv>io_WwAI>VW8{&=F~_&4Ku8)vL;V4Hy7-y zYJ(`_R?9U8A}OX%HOdb(cFK;%#Dli6lLgcghuhK5pA2)p8S88h-DXOD%ac@sdKCpC z*yg>9OGjR8io{-;t4)jdPkp>LTZ|4L&h|dz$kH*(V@(24ZJWL2=@yS;cP<9eZv(|a z4BP&1NqX87+0}a5leUxC4tr3l=6N5_Rb8_pqb-c1L9a#(RHnq!*=u{nGkyteL&ZV7 z+XHW%pZx&+)ZL!`W;5At5WBr54ieek`?>rkK`0L-v3>jYn5r7oX}gV`$!@vw>HO<3 zXS_vT5=dc0R^HC|PfD5PN-{`!#GY?#pq>4^U^s?PO*9LnF;X1AK+1X;0-8w@lIdX> z5U<*{Qn1aNNok%I=ROAypYt6b3}D(Ddw<)GisW^!0&Ny%=`ds&n$!%npFbUuB{4JB zKc6DaYK}K$FgRZrb|-l@KQQQvo11(dZf`7}19ICQ5lDEDIkgfw!rq=2IF{TEQMnKO zmSK};e10{&L1Jlaa(+LIY$Of>qNCoz2f)D{`U;)k5XmeBt7Yi9oU5Dv3~QKOu8eNJ z`DrA>7*yZ^vb%xlr!~pHQ_*Rz|W&_ z;G3!#v)o~KeW7mAEKhx%Zqe+{V!>910~g(l$#kHoe)FXKfwwqNLVt(604P<3`K1~t zT}5ub{E4%g+1!ZqAXeXYne-r0H)BfqAXT@yN0d#?C}l>Z{wNCgNu_cFL_22kdaM^* zA+N7c3e+OGmT*Fn9DV=Yj)q3n@E3YopAJOGx-o4YrdjIK@q8!E8!&7C=r{n}Zo$u5UScDD;Vgxm28)ch2SB_h(I;^Oq z*LM^&8h9M|C}6*g9MPAAuasTQpCZnAOQsFVXcE8I+Od_`rDribss&uAT#&bYF|_Hz z4cE>O6naKZBzfX&orOL-UuPK?NJ!$Y5@Ox=Sk1PB&-zHM5}Pah_M^l4cbYF*#!>V+ z(_NRDdM7PPd#^i;6m32}JV1wFd0WhlH}aWofXB!_OVq)E)C<&kaLU7VFy-THu%Ds7 z?kJ=K#s5c&QTUdT#G&I5?$qLW2!3o`dwhI%w0r#%fsisdXI!W%^)J~rxnb;H%FVJa zrK#wCM2*G6LtAGugNCYD_%iLaFJ|=4bfKyw)>S+&7LIkZ^1*-}|nX9x$Z> zm{Q$mXXj48i#lLtjSpPZT5pL*ZF>Ld=Ksu)@{N6KBJZ12Wv9VDIrn-%w;Ds@*?1Wf z*~SmGBw_DAzEPM;vr4fBxlFF_$|A`u>_P13JYecjDl%kxU|w0^B=TcED>{S;zqgWpTBP>8W?9_{~)*`NWI%u(D^{P-y&~rsu$Mt z-G+cgP!GdFoLgxSz56*bfvo8Jxsq%0ZZ`Tw>u>2n$-4wr!5kf~iN(@{9vm;GSrZq{ z!Mq7fRv)bG4^=yoT19O!zk9%!q*OxqsSWlH?e8<>PLy&v9ZL;iPp>_X!jRL%P!WNc zVAY*^wYV*6FM1LW>R4B#aI6la>*phz40E0dK}=isloG>gLo?H`)ao)3rJrERCsjoA zZup2sVYSaDzhj1wN7-hw))(nD^c+KJ$cDs(mAL4h zzbpPjVLFgU{gTRte4J0oT-~Nsin3-n@vK>gx3e#i;7@qD`a=e_sb*1+cxE2$7Dy`5 z?Px~1KqIy%2dLMXMMDNr6dixEXnO(jX_yLJh(AqepclBaB}Bo&8E*I8%CQhgF6-7c zaiHe4`|^8sn(DRuw8(-NIl;%AD0^qWq$?|ctP$Pb zd4@Ly&kDY(!fpo(==l}?xD3{b$5SuntPaZ1qURy$=$7)@`+RE>dK<}IZ8#)#u$xS> z8FL+fb53GJpZWTds=DVVG!>Y$8KCFiy>-Wx=1jF-_izu|9h>a|?O3x6-P*v<_U_vo z0drTNB)`F1k9} z;L;FZ@$ExF=uHg~KB@o?qC58_{@ftB37`=Uj3hp9hD5X*{b`k7q?3UhR&-fbEzY1J z9sFzF5N#M_o1@Z6#+~QU#G2o~lDNs?<9;{#?dvIF>r%ep57D_L@&h3WXh$|!X61%r z$2k_q$vf(jF2a1}4rACLzsOQ4wstZ4_JDR{sj0cZL=5eGHdPv8O?hE|&r}KT@x$hC<&kZ5{3739On)j$;4wDdV1^Zz;>Q%8&KIL zD(t+%nj_^r?$EJHKEmI7--^pLb{OwI8Sr%@UUUE)S_MdI$Q+B z>m1LD_~svFpmo49!Vr0GJ$wC>z;1?7SID07Dt{W!bNDkzSfe^ygy?FRunRo(qQ3%G z@4PQu?@y)Zhyf{~*a&!87T#L?(n@ruG>oH_Y?~b(>Oa8QD#lZ);bszkG(Ps}YJ@t* zceFg#SB+}%1gT0A0hnx9+p!yrw-cb@H|meX~xJtLbF z5_xi^=7Sli{Q2%9w9a>(9l`G^f%rxp4-_6$?-SOI-6L^Q9v${uY$AeLxO7O&F;G#p#`*S6+bi)SLPX3;%*dKwpl+1d<3DoM@wKsO5LKs$)U zU6bs7&~9M}pq?SUruhy{kWGN9$worqvblu`fy9s=9_YjAAkQhzTwp)w#XgLf>0SwTHyi65{U-G364hd2 zs*vMCFotZwT9T@HNhKoHyT=5_Sqcy6q($(XK0OlaB=IC1U!9#p>Qc4gfyl(ju4A

Ej2c?BBDB?P(?OfCKGOpru$m4px>JQ zm!Dntr=J8I-S&>>J(9-w9NB0{@x+zo;o<*2Pl$AERgvh3*B$?2g?}mhvkbBP6?@8} z>S!HyAl+_I%TqinzL)HoJ4R^kM4tB&FU zE(UncT_l^78+4f+S1uD((3vzH9I)+wkGwQG{%WSP556=0Wu%X}PYY9XY>Eco4gX|l z^&yq>S1?vsj*u#v>0(qrO@Qn$I;*G$e^p4us~b9R<^_XPDq`MBrGk~$X4+8lMwPN=W!?1v)S-M>vJFl!4vkH{ULp7>NQRNB!Z+Zm&@J`+i-sMsxxm%&)XN_L2 ziy{){rgVklU-YdRlL9jO5(cVdNr8-8XT=rg*Yfe%V`>GXxI)o*Q=m>vGiu`%zl5@1#Ds=7}v!Zb3#;n z;&CI<9{L8t7_LM1^G~{{;hGcPq9WjPI1S#EGH6#uj?m5^&X;C1S6eFdz@LrM^o!OpKcAHTq8V=nFVlR37T&KiIMOFQR>SYv=-*dsF*-|xm-Y5_Pl2mhNj@<&{epVWaV7~8Q+(eJ zS2-Rf_ofdAh6Q%#BqkK*Cn^-$K&1*}Tmoz%$tH#j?WLa>H|5HMKm2&|PMY#fHpZSaqlVkupjXDFaiip-}}j2teP18 z`fx%22OyJou)VPY7T(OP!A4BYL#0wf-9K)L9EO>EJjJ+_3sC3Jk@a+1lBVfdS!Yo0 z4ZYlMj>BVOtDVl%AM{WzNZ$)Q39I<{S}X5Y6=5JQe0^)@EHPwA_At&UcJEE=x*wMG%fAY0og!SfNhz<%pPp;4bV23m`9vqw?(Mr-Tqq zW1P(OCF_JwHX}3LwIevtu@S#Sq2SWh#lOiWuhq{d?+Tuf)ONpPc(0DD5|bpQQrpF& zlF>c-z49UJ&fG|&S4xGU=poBNPU0;MrdW1Z2qbGRn1>F6`zRnuBu2}F<17>}=SQc+ zzg|)Qa-IC7?L$T$g-e*u!Nq6%PblsJ#yit@X!~#9^+F#hN$7TPA16-|AmIGRPd6|r zEYP-z?07ZaEoA`%O7rn#kqCMJS>*TOOLmB4&N9q0*+ioUR!e&qc1$JRc(=)EYqm*y zwme|B=8JwiibZnM7v$w1sxH+bkfGDMRRoc zdU)F(>PvjqDMo+0=NQ)l`W^JxZh05`f=8M_9P+ya^e|OGFFvvQ>g-2wWem{py0RQMl6xwuQ#+Js4N<;7hGwVmZ7 zxi0XojTPY^7(N`yw|Fg^K)rwoC%*y@UQFjE_?2@mFpuSfrB{E02Wy>PNjeqL$x)^( zvQ6bKVr^y^u zlzJPz8+WhdDCESVQ8MX#KEK$etor$mic&z8C!Nyh0NtY8_aj!DUiMKl^`mK?Ja!wB zNf_rLU1H<<$O7alE?dl&b$Vrn>tJK&b7l2e-7k~uy3X3T3hlifQtnf0%>%wW_3H3A z>$B|7)lRVF2T3xurY|w}GbLG!pB@DgB$vA+J_@8Hu?vVH=!~{bURCO$)>{GEbOU>F z-<*_GBP;2RZ`UYARE`@j%*vi^P8AMc>Hbl0b9-+zg+fy_H+l+=&rNc^(^g8& z?lK6M=^&joy8!`ST!s^(FSkTQjE|-E)ersL4EbLQ`^wXcypnWt2xvY%%1i05|H3K> z&7?QxJ`J!l4yBJL9N^#uNE3fCiQN4V>6>)`H@93&w=pv0o(DtWzn1o~k@6`c zwLCD&052xGJU{WI0c%R{F19`BlpZXAIr{^`g4s40a_lf{e-dg~z5p{j&Liwln|2Cl z9buN6_B^hSz;2xvce|7w}4ZRN!e z`kTpUo!{O(`cW5j4dGRqx%w)07Zu19T+8WNelG-i!yiYw`IB+VhRVeS2b&fWJ=bG< z61BuF_MV|Lv!FHjAyx7g?@=%JPQRN`iVkD0eo|?H*#mK7?9jBx_2`QsPXR~cd{ibE zo$7&PQGuXd^y01~Gr=i;*m+X;o5n3Q%lg%=&jPyPQB8V`_nQ{i8C`)sKS-DoAKT9A|;=n0~2ohZTJ|ErUva{ zmV7!O2kyb`FEuLzrgY?4LIi-RsIUaSi#Wgb7Lq}Ulv{rJitxxZH zHAV2Qbezs#rKz%LSe4-&L@3a|;0HM>@U=y>Wa80o9S1f}bO8)`Y?)SRGN(87h633$ z6}|3~DjF738+E78c?LA!Pj-bWYe{BC*JV{PYe@9;zQ2@8>88{nze`OTkXbwtKyfDF zc`4eQ>sfX!jJEvd`l?Lt=e;UXk?dh$Q|Y0@hNiVZ6j07mGX`o-UNX*qg$;^+U!?ep z@nJY8L%OdOLbA;079MM*Cz5h%ta`AWKOjuxXU4VI#6xp#g1PQhHX) zmY2_mSIqI&*}ng)P1Ba68lQfeMf>h14IkMw?7Xkset0D~BYNP2Au{36(h{Yt$xX8z zEJegP{fwBsTee-r8^Y3U)zxWa$q~?T0=JI{Ngq3WE}I$u%hJ}Khn9qfLU0vZD2y=g zFw1>FBAih5bVmu6KVa(r#+d`v|L(g5fNz)NL129AGzHM#8J96MB4YW;?$-GS0bSw) z+*@at6S6nt;V5&*>8-6W#)OLM;&T6U%#7ire1{bjgqnMV89hZ~h#FU1uAk-32p{Nb-36nZykI5~AV(Dp{8Gx3A1=BI=iQ zWGf?&A^XE}7?wC=b@aQZfbWl{Uz4a5mPQH=axOurXjO5JWPsOzSzpImZnfLHGFWcy!m9k{1&cYne6_jvW^>Yp8vfX9wJ)AN z3=CR=d9^)zkF^tw?5YSa`vyB7tCat=#HXFIbxCw9%4Prckkbb1ZW=y+6oOW1j(fLF zR%w2>c+N?bjNAD6_P*JSHTf5;AvrsFOgHN0TzTD78>7R5b)w<8fo$S+Xro?Kn(_zj zdp%+96v`jK-iF^jjpk5VbY?x*KnvAf3q9VOC`vKHFDK_)y4qnK!_ioA{)6Un$=vk1 zq=Bfa)`u+0HHm>}EILnyvFOgwst+Em^7*<<>^7rMs!p~MIqI^qdWGKpx4jKAnK=Bq zg=aQQfm;YY<8aEshT36#D+56bxNxTXiY5Ggv*KNU$013p-}w;W%JKAmB8PNNC{KAZ z;|s!|#1GcsTyAou(py@f+tI8GjnRPn06?AiqfNHvl%_rSeEYyP)-n9VW5P$I^&8H_ zF{E8pr*DY!otC89Essf6o@n-#%^0=IedpWsF%xMz<9M4}@9RcC#-+1^;i~dZGuERBTt2nSx{K2OIUPmCnYmS%rHu%9?9hu>c?JZ_%T$FE5ioxxYg^{ zHW*ZAn0DQ|ikCiQB1QnGX`rXPra$Fhv11$}ohHo~a!bt>X&gk~&+xdF&gYa6o95m3>1T=eaI~eT0+}df@=7^3IwaNo|9BU!$ zdtzn-javIUz$EaF`#$Td(st0!98=91>7DcAJ=G8|>EhqmeB6IBKc)X7_W#NJ{z~uu z%IN;LGQYobG0@QdXEEOYA~*1#T=B1lj|KQ2n0Nou2meHj|G+5QxqA?5c(_@aIa;}W zbh36Pd+!W1cK|(iuyeAvviw`Qd;HBU6cV5TS^nz%ruB^+=yfP1#jX~j91)_f0a?#P zQy*AV|IDc>rlqDgZRiRBH_BIjs{Vn@FK1r9ZREcic5o>`qWF4c*5ZS{v$C#{(s+?@EOiYS9uMVZKHE0m zo8Ei#Q}FZ|D#oV)4JULBIzreVOF8e?+%-9_wcmRyY^?byY#tqGdA3<)z^}8!L0nm= z#joT&#O13^VbOvU($qcs8 z4)GaiTzpZTi&`aQQ`IMbiJKk~F#a4I@zC}HGJ%Cv#j8;=fTbntY`)t(4w-oMlUW=X zyh}uA+AtfA$Vx<6@-g>aQ1Zywq@AuukNk0Ep2U-xUz0$EHP_Gi8>uTIgR>NC$e@XC&2GkpK&NIM; zzn%4niwgJ-iwE3JtrLlPw43IRG~p^ZyT^O_Cmg`Ai#TM-p7A<_#~snq2%9|f;j=Fp z)B7d3Wa}}%Ia)QRF5;oH@tmcI-Q^w~_UpAtPM;EXkz1%wzFZ*gH%<)qEeMgjmT64@ z%ePtDsdEh%@}%HK{X(^pG~v}@-L;~P_n!Bw69Z9EQK#>?;|CO9>lNx5*LMl?)jLVF z(bnM*`4`y3&n6U~HV_L5+|$yufVo9|#fcypYx!QBr-PXI{<(pH;;ZpEfgd}k@qL7e zI$*7<-E?ykyXn}A`@~S5156|2OqroeaCv2tF&n1Npa8uJ^QA~{xv>EwCrA!#@BttKQS3(ONrSv|d$P|$1C1*V-Lo3qErb`7BQ+<1qh zT`%0R1W3rk`SD~<3B!|bWP_o1qi?G4m+kU8Hf`-I@AOp_%hNole$}+Sf`{`;T|4-- zlI1kP`rPG@<@!0ik^{~U(}9-#fKt(2L&C0w9Ic2U`i=bbot2;FLMJK;PM)xJ1IPR8 zpFHZiDp*F8(o$KJ7q^`Wp+4UE(eCS=RV@>hHjY_VmuRZH5r?NmATpR4{o*g{Va?~@ ziGpDNU4K3oW75^J@4-lVuEEQ{e(4GRP#TreB5;z@O&Qk9#7^j_ z^0F(!W#G3R$sXJ~j6lXTsyOBO`f}{kr!#F^sf5PN#b`2EcbfV*DB7d`CJG{+tacHH zEY?DNEalY}=0?<)Na`I|J160BoRe&Q4Cln)A1X}DjUpss{522?gVXRgVOoCrdJzkm zM)24ZtS*nCISs@cLK(p%K=%*zpnNeebbgDOpdIWR->g)x{?2Kg@KhzcC_*u11u|G} zYY%#UG#oHi+2G%aGYmGV@79p!U*ezC^nAY@AZLU_xSfUkEc+*YNre_I*aFY)t2{7S zT2!Uq%!9Nzsde+2K0j4Du|Xgurc+^HrtqB>)Y6mY*>p4wWj1fVks;TYfGArWh422T zJOlPgp+7W=rP8611n+qggU{VxEZ~yeMR6RvB)_dJ-P1O;s@X74^K4!*nGhbbr1+Xg zNzWe_`aDddiiwEjR$GX=Dnfxyk%(Q62se4JFe2bJ@1HJHi3>B~i+7Js)1m05!wU7k z#1aS%MPxjzS&x9+$RqvbynKnOy_%t>%%n+nTU@U4wQg$)i$4vDq+h!L>c@ z2>})AyYz9%>SJVzoiseCo96holb{QIz<=nToJ9Ag2<^hrfDGdGwVuM+=U3O)8+JC- z952hqPt?iG9`VpWklMnMKwG*!_dZoofDrzBMpJRRBe;Gkk6SmUfzqQKpUY)C^(VI0 zPQhpqxjSZq%b+8zai+xb=h#!F47&0rP8>A9x%otMzLCv%R?fTw$Bu^zrJvv|?oL*J z{<2^zr}Xj=$>-$RuU<-!G?~b4;lWc9@A_U@C5PKe^Z6^gfSTlBPOY+`n=|PV5*Sx# z768z&X@|?0rV!q!)?9k~X3q zHTPOD_`Avg_C7Pz4Q`^{oG;M<@?=ENA5(yl1p-X=f!vnVYE z=GK~ivdD%6*|cuo;XVriPS$8T;$<#Bt>+I4f|70kiu(ElCjAm7YTa8v$pd7-3EMYE`+*g;jP=hJ5(tjL`#DUzq*Ll#~I64M{!P zXcW!w$z@kR3O>7VTy!Dpc*|+;Y&VWWgYEbG>3A8rLh#3kr-#!By)+3m?z^&2v-Z(X zMK}fz1q`MsUX;Iyb|kh7A|1|r zr%3fpR1L@Vb}a&db1Uh2?cuxG&yey{EsOwx;_tJN-9`ByM7((-pJ}kG^aA#?!+Nqo zhP8Zs#f7;By`&O3An|fqH5OSmvcdMHZ{vX?;3mZoZP&w2CR#>TKRThrj@soLPZR<1 zU?%X{hFCItzUw!eCR@Ub#a}GRz_-Q=l6>4YpXMXq#jylMCB0F!`S42~ENjV;E*u#$ z^CX`oACCR0ZJBAqD}q9)5*(PZ&zBB2`|!Cl^3#asf`yIv%VZ*BHTrA9bGr=Tn|?NH zbvDd$o@ma2)PVW4)LNdo00Cnfja+IRrz``}In{vb#NQll)~s#yuSxakOW*SE zj-GmI&Iff_;2qgVyJqO#1aTB`MRAxsRnPlS`{o&|H@N{kU)K+(+Wm)@fc_n(^ma_^ zb}Ugqp!sDRIx}@&0_{54p6@yt39F%5+d{Um{oM%FP9To^D%DC4!AX?IVwpob6DB=j zOE~&sK=7E-nj4%e&=YRS?vaZUxgyIeXljMr-zH36XOVu%vjag2zdlAig_A9xr%T^m zisL|HpvpHw`njN#lfT|omrez)jg~Zsi}VR85mg34pnSIbT0{84*F~!t)!?`W$Ni+` zydqR)Giq_+c$mMkaa`sss7G3X`c>YtCci zDg2e^jZW;&$gEswwCDDg?C@U@1w^Zw6JDvIHqg5*@ zf%;(FUie%xzf0@s^T$a!cctY9j6|`x;>q$K%QjV$Q9W6QPZ9r*vBcKr7b)Bm1!{N)A$#9jiaX#WSUgZs4r literal 0 HcmV?d00001 diff --git a/rootfiles/fonll/LHCc13.6.root b/rootfiles/fonll/LHCc13.6.root new file mode 100644 index 0000000000000000000000000000000000000000..a15fe1a7ba3b08c3c787a2a6146a4eebd14d404a GIT binary patch literal 13715 zcma)jbx<5n)Gbc%;BG+^5;VAbfB+%5yUSt=EV7Fe+}$mZ1a}LI1$Xx)KyY^!7GCmu zRqy@zeN%O(ds=RtnW^dOd+r%`0KgLo>0Ss43CRixi9h+TDgLd0{A>EeF1e|O98>ACS{dH3ii-E+Ei!zj0dfkG^Mag#O0 zbu9Zm0fgUSfkpyc&Y0R5OZ90XA-_5~KQ_yBL>(BD_4#2yXdnSH)%uyQzntkAa5nNX zT*GmCQmMHAK4jHgj#^<>{ynR=y<5R}2(gSuyH@y4f2`xeN%#Ghij_(o4LkcsR7l#_ zB=7?icucMKbd3jrLj+Tf;ar?O>I+~0wYBcxiQ1Z+LJ6V6zMlfCaT~j6hZ6VzdYHoX6Pzp);>Fvp z?*%h1D#Y4Y`fitIPJtPEve**V|pthVDB5!Wl$p0Y)KQzUJs zb;lJRUoIz($`!}C2<;{MobUJfyn)Rhw3Zt|J1Jooq;jM50nteiHgOGE{1fQLLBeLM zmvBj@GvmM1{_t<}>=+EM)Hl27D&bv@uBd<%zXa1=|XLu*EiA z=LR=fjHCS0w2!NXExs{!yL>rNjBO7=;MThbT3f$#`iTw^)=-)->tb;q z3h_JSuU9u+bLP~ummOyDuvv6)I&4Y{oJcIp|Q!#B!$O>Baf)&R_H%gUPcGDCj^~~C*@F^(l zLoA8joi#GyX8;E=Mlf48Y>_cL?yAo2)q0OP7n`CnZ^i_>r&@{X#Q$*O810Lg*)t*2 z;C2^em3TepE08!*nZN~ojk@Ib4O$ugf5_T z({)2kb2cSM{1KN!L`KRa$ub0_N_xA>EH6Ah6dC1{qu_7fYm%UH{bannCBVW~4%8Ak z7dCvesj_>dJNkwwK_x9~dAxT|ah>Zc)fC$oW9|1EQCZnx;-SkgH#0qXL7j zMQjKi2Sa)A=@b)vZw}ZsJe4+430zN2cyIWkCS)rxUhE_o3nU;U(T`?|=t1FAcH%Cp z#c~m~48m+)+1SS&DB~6h>lS$@RwE4{cbeyUv7VWL)-SOcNYG&xv^LGMef@Lpk{lK4 z&cN2ExA8^o?}OR5M^;-V!x5c`zwuSt zI4PJuh05;}kY`QiI~84(MTaCz;A4fP|9B-a;#l=COie2-HR_ z(`X|^(v%xJIvsqK zLLzXA3sgMF+^FqCy8FcfT-$#V*dA2E5nEuv>LAfUfa$`)g?kn&3lUc!qPBaM*UGKt z@$YMf7>&r`my)i=GzG(F=!+@$13jM)5#(~qR z_jpIZWwZyAhVCa(rj0Jt+Kqxn6#qIY8I$yvbzOz@A{#ZmC^YA>p-%Jx>nH;Xq6HEJ zvRiq$GI>IuVU*zNGnaowvBz;#b?Lh8R_yWZR|2$kjGj2qhp z0eosV9_8kqRU5fN$<_d}(3X!f;E9++rj2P~`;803Ss(}RUjwKIga0TY+eiRHen1iu z7Hml|x*CmNZHV-uER^|d82+^Ac50*#tLrr%q$bNSldiqxAnUAv15S3N)XC;PFWBo` zb;N*437WkSVA8d$C9?&IJ>ByKh(Us=4Q6*!#K-)?qTDPf9arBuNHPhwdvKZwVZ__A zZ)Z)pl|geqVI?0GCPfPAb$y;(;cy=G=RABfd+qPy$l7kV>tVd6VBG8wfP3Xn;Lqi+ z?ho*f-{(WwgYG4@fq(k|5yA7USIcxbrFd!~V!K-}?%;2+?*Fc;HQftQNnI(m_!<5Y z2a;Qw~JBo<0Zm`hwLdL_ys+EMCYSF zIy^o9-HQjrFws_<&3YF zb`Wehj#aZkNwD2H;bAuB27J~=20oXPYe9eYT6=@&^Vps_JsBt3j03BWlu^s#o4bd` zEV|v`IaHzoxF9Via{sTxX!*lYfgyCH8-X`ui@P6i7sXRB`Xy~ca#86P3yYIh-fP$0 zm<|J;WrmthbdK>nQj7rkRHxzZ9K-^^(exwlWS=-SKw?;8mnLYlPl-Z0f=KtDs4-U; zBww6R?gojH^#@wp&H?3PKMralBt-ICJq~iKNB4OaFfQ<=mcUShsT>OgT8dY1Ks~|4O%3 z5VbRXY?+3x8@csI_nZx6w~S%r1K1h=6fzmPe#*hT3Ft+B7Elh%ZC4ruJ7?>e9S*TB zc0lo_U7)Rb(+e7hb#LK30{||>J~Gx2j)Qm_n7-1=`8geY?@fL@FpPe;Bu?iu=-C~D z=ac~MRdg^;1G_3Wo@yM{doI|1hdrpjxO)D%4WxXl8_{aHI`H%ab{Qr-{zuJGo|1<|PhdhTZM=(b*M>oed$2%v&*F@(*JEu8k zAZH)N%}CHl(MZ?G*2vo^!YJK07yGBnj|H9u{srL$@mu5*oK?0KfGhp) z&Ypnq2mD*c6P4An7JwW5R%d@e+ymXM^hxw8TPwhw{!iy{K-vT6t@+8wYFR75gMQ(` zw=uszAb%HW0zHpdOwcFp%;>jpuSdl0%L$e|P1I$>Ctr?T;eLh9YKW=WZNA9lJM8Gm~1CAZrq)P#P#lCFXrIc2sJ%kuAr z5&NP-d4GjikAn@UgS1NmDlYl~cS!g&IA3+C~zKY&s zn*tQk%R{5U^K(z+2wu27Gol6d={|sTnV_P;^n+v~m`!rT+-cmiGASPa*Y|`Q5$HXy$#4ijk)AGNnd}nAvxh78RsU`k4LF!W#pD~1EcV^ zqDgX7>VHkyY(db3mQ}dE0wpe)yLsqrH4w z6Jp|*H=|pJ9TJz1-M3;f^j5S;hbVe%lYA5Br>4e=1&7${^FWL_AZCT^BoM1Yv%JdC zj5egoana{J{g#W>8tNf~r*R*KFGA&cCP4%l$(+Bk)Uxr%w9z#LGpqIU`?oY|(vmJAfF*qQ>0|L|hA zT-C04m;y?10F5V==do&CJ}eNnwp9(z?OFiW67$F_R1B&PQ}mkqL|EtfYsX*9E{s_1 ziLlSh*Q!MUH7kmzf!Y^S#+RoMyvrYJNJ!+-{}Bjz{}BjG{|E#^`Ts{Cptx%N;}1;9 zk@EiI5A-m#QiqbBhBR1|w3wAX$zer(xg%#U|MH^f3zpIwgUt&9Mp4>!et0zzJ3C91D9^DU6}QK`g}?4UsxDm|$vi3+WCe{D5E^~S zF8aZ4Kx*OU?zqt2&&V_}>WsJ^H3Ei0A%i_RTTa%z#0`o$mJs&#J5$qV!u@>GI0Ova z@rf(U$rwC44GDO#ErI9=ozRnF0PX|{(z#5}CF!1?N$aGbh3Es%OIgnwWD}=MLeD+V zdN#kcNIc}#-#i2IpPi{kew=hdrV4pl&maBvHuzM~;@mra{PxVN{8_-p=}PX|t^6)% zhUV(ayU_8P2lGp+G7=9F$s(D_tb(=R z@sUw$Cl=2~62gT^w}^QMA#oig;>Rg_g$crludIFxCgFWwqw^vey= zqj};CIzFYC;zj{&BFaoG7d8#|oPrrdktoCWor2l$W-m%|z3z*T%CyYf=R2Z70dE+q zz4tJsZ&&^0&4MvU$m>Nr@3Oj1mMYP`aN3x-t4@rAjDdvJdO!Jt<7!u71Cr9M>_+ZnQy> zHA=(6m1t5n`*Ed@5+Y{qZ6+@+-OS#~P;>uB90irr7W)b+s89f7t9V$FfdZJ1oCVMi zRnukj)kk72ptpv7KhEPUqy@9MO$v02vZW@ zb@*U=f*2=3pIW8LUJU|y2ZRS~^MXww({`*NOxYSuImDRnR?BH)+M0^ z02Ol2HAAX`Ct#3{xgwf@8i;9fGmUbR*8kE!htKbX`OYTygIVs;=b%Qv31Yn;(aBOX z<6$VOn|$b0WtRCGLCFPp&DmJ*TcUxeniYqNL=cISQAnyBAJq znCvQ<^cF6LNgC>!UfMx&nGkk)%r#$Q70iHdF~h!xlVCqeuql~Y6Jk}obT0Z7M!k2@ zP!??3gnCz5lvx0kC&hVeTf!>Gd`g$zfGFhbd*~(UJ$rr$k;(T@LIbZ#O93WJ{czBY zm%nRmz%GT@ucR+v>C9O0pFEghPp^pXT>NXhB-QN$!_`n={UP9baivLsCpLPg8u-oG zoG^+2a;MRg)+Cqe-_KV&8p$+6+_*lcEW$aqqch8fIzwp9iBbV(HQ!Oz>D zc|Jvi$i&MA#j`rzk^A~+7onbYyP6GwBCASj2-fm(J8s$7NZZJD-=3Y=M?O(WKX-N|vkM2cH=g?RKK5v5MHD#L#vhRaehcC~%aFecf|?{y z3@u^jLNVa}VtXFyEcmF;6%;V++|+Ae6XQHYIPThX^OVbW0!QQ$g!DX_iQXqTNeP)O zgK&{u?VqE0{(9Q{akTX3RQ;EwYPC_Tdur60EtXj0R%(LBcdwKJ^&KR2F$?Au?^nU} zmt|WK49_j(H}o860i$_8H1_t0Y{=>l%NujMimdgxzE=|-Dhf!CCw^V@)`NSr5O34^ zHkj@Sr?A+|=18{S@8UQ8B~b_-@MD5sc5jc97s(GFVhU_e z30=};{4ix;CVwB*t*hJUsraUT(PW-f#VHs?xV`kPt%H9u;!!Kxt^}%oC_Ev7_uA0Omce*OG-CP`B+&H)5APnmmZB=NGRytHH$j<}= zdOB5Ax?r1u@K@!*vS;-}b=6y*8?!v3(IbO1Fs(=r`5v>^$_3!0rMZ zXwQ02K+<@ZZy~{6BZ1k(a%UI+1as8k$HKW%A*cwsuP*iG>UWln>Wi+YR7FPb(ks7T zDrdT1Y@=m3oN7BC1LZQI*7UW2C6qzKd77(fqCNfZiF?Qhm8z}BB)3IqA8pH>z+E(( zIm`n2sbLR&p8Q?ds(dc9gE>q_AA-ZzWx5xMz{#@==%yMpo>aMT4!S(G2zdG&iB;p# zdDIIcs0El5_B|y?R`c9Uq+U9lyihunnaC+U_STQ}Q}CdIa=403r`qIW4OvOr+j((q zR7@P*s_rOc@h_ZzIBG}9V!0jb=wlyXdhbyqm(1t!-~%qEY`-HYdYeT!5a-;YGBaX6 zaV4`T`V)2K=9)a~#(tIWYv$L=odeOyh7Wt$2IGS_?RE!bd)dn&*~BuMu?&D=p1sVZ zAP>#e7Ja<(g3#!Z(qQxi&6ilB*q)4l3#G~A0Z<*RD8SAgXOVv&XOr8W-qvrb#A#7* zh#t-gMixJRxvieE%hnVGBL#o0Lg}KR*6(h( zh8x3I-JpSV>0v(c9>KW ztg-oXg6JSD;oGDjGXSIYZ90)(wh=4?mwn|j*|R&{BT&v_ZqBl&K&)F zguRpc;JryXKxz&4YGRFGK5n7W_`}-e^%@;Nu;3(j!PtVWyj^@8=+`?7UR7Jef0w)V zM~NHA@c}!RjzWH|t#b8BtSYCi6>t29&xW#_Y1}0jiMQH|q866X=>BDr^tijeM))8u zm4aIxU7JKmXBDPv)#3X)SR>b!_jSlqRrsLdk)_MdhbZwtp@DHRlFpN9`zY199pr*V zoke;_WN{vU+?rMz4>IsT1(p_USzlZ zX(QTaFjzz$B1NkJ_`Prcaq6||B8F&G)kIO!e6$FTx9Yb`i#u5BAzfKLhoGh{-yQvR zHexxsnDg^2%^e}>J2Xu#=f1`n0Gpz*&`Y`89d%o*Rps24#tM89eJHI*I~+$ot!9fH z|I;>AA3v`)d(j6jIvdS{yTV{mIa-2OVA=sDHQJSKXk49y=xVob1pgw85#e>cD@!=r zozKxwc*9>)jgapvLs-n(`7Z^qW5wz)xY(<&r)^)jvEFw!MLuQuX{fDtA!IpE=5blsZv{VWwaL6ZqUmODq$a+j zgOB~_)gcuSC9G1Bu2(L@{q@Iw@d_;X+tI+r8JHvc73Bqy?`ncDZ<#`)4f)BMkmTi;fg6&Wtotor=Pd%A0%&PVbR#ljQGzc&bdIAIx zd9+?zU3naSxYBU+$Q$$^@h7+@G#cz5oL(Zy8LzWSlj-iv;HrF+_^czw;Y8Pv6#{pf zBfc=K{nm3yq7@}fh(fHQJ^6O0XtUD+C ztM#jc;!*-QcMXPT0Xvs7_S-%`Ms6H0%qDqp%z1gqCZCX&2BH2Hi5P=XX;@}Gyuek>Rzue8f-0!*{%BUYJ%*Yyf!t;sj#1`saqg6$Yez_ZHo^0oF~_g*ulDLSKtbhHa-0&ydK1 zLw)Y%6louZ#qNjw;dK2!w?;ZF*oFi#5j)iBy%_E+_i%IYy}O%n-lux?@=F^*mY)QN zV|s~~h%t9=+bm#2hs}lf)F|2YOiH}>jo(*$!X0HzlRPIAe&>M|D{g+$9g1Xv5&X*> ze>vajK{mqc{jE5^h9Ova1ZxD3Zq^GR{;b@*GqDDp&h}WD`_J6%cRed{9f}rlei0U( zJ_|A?ROx2m_v^v;u#EQJliGgsy_J5OnzAY4(R79OLZ}AaLAA#F9c!jcT`r!R&7tUK z!b9@v`*PcAt=VQ>kG}4amjT=NNFHMz=tcX}ox4sBwDx8PK_p@*$$ z5X{HUIHQGV6b9X~n+Co=?}qr#cN!csySJ%=JKcdT%l?i%dhXxe!k0asf$pvGzV5H1 z9zXkX|JKd zcVc#zNEE4jE?V1udn#``7TYU;1xfH2whRkk+I3AdbH8D);HGj<$%)Mb1z`^DiysW` z!|#SL|5Ev72*?S&o}RRt3O{o1I6_UWRebr#jwdb$UP`mf&zP69TxH7d`|G9qMFZzt zzOkIz*Rq?Nlfus-yVAc-*L6RKCsE!^?6t^3#WB~ z2gW-8gkAa%CYm-M<^C?`<#F~(DBgF>U5kKyMi6z79{|=qKhH(2XhGgz_E$;H5pa!B z%(^hlRXeWOf$EdBPJhkMtZ}Lp`eJ_7A#HUj@Lk23i#hZAQ&hjq%O-7wivLipsk+s0-OA5>J-*aM*s8tcPfA{*>5x~;#CknAeNezc?r9A( z#-*!7ZaYDh@?b$>oILwDvHn2e9aFXvf4~mx3v0Gq-cFv?Oj6higmg*IX)@ZQndR{I zYy#e+o=)E2AU97P2RqQ_j;(+kbWvKZQC}{JBQhDPOWepQSMu*D=|ChQ5+e((3^p7 zo*J=@XqS78V7+w?s7|o9NY_A|9KTXY1^vF}&=-mKW{+%SQ^{TZuuw2#b?!cO!r920 z{Hq%Bv80!(ZsU0HV9i1R02Q!G#Cm?a$?m(Q(**tXitOCa$~+$AU=r}mrgLuf6ti}6 z^wS!Q547dna*@D1Mf1q}$~7}67E~r0JSuCGLK1K+zr)#K9G)mlZ+o<6ck^WbE}}xi z$}74I7C#y#K>z;RmXeL*DWb#DSuAttVZv$etD9wbTY|n37g^wd)%^E&>fD*aTk8UR za&(7c-K`?+Cg%!8BA^aj1I=7J)HDXR?^}CjEi!)0b5DV-Mh_$x#(s3`JdMe>AGxmf zlp7KH&D40M%M}JrYWXDw-pwC{xAN2*tvVc3ztF==gZSp6y8X-l*z|~If3z0QGEqQe zlihGO0nT$iE_w2<)5w2)i}UMLBEdf7q+@0i-ZA{mj#_uI1Ejx00A##+lg z%~uNqhHhZ-mx)Dll3K>RvngShoY!pCd>IzF=)Th8vW^DMzizCYpb`XvP^{d z^t_vvkz>S60?ZCSZ(kJ?HaR6MCurDi!v;fIzv-(gELVu2geU)&m1tT{_?-wRRs&8Z z6#A1>*s$*t*qkuSf+VRgY}p^gE!j`mLfEM1(#a2#gIc_dhs=a`8oJqMlY{)@sqgBl zPHIzeeOkHNrQoNnOX5nV{pyf*?v$&t`tog1k_6Oc@$@C)#tEJ=%NX!Y*A#M-SPacd zUgzjwyBfDt+dDU=8~6}!+D;pvv@S(fIjLTT|)L2xjY06bk8Au;={`#~`R&)SP8 zalF;s=`1~AT<7W2#j9NX2Qp|tQ~DXKP`_m6*aG>GVVZEk3VVO{ltR>PMdt4X&|@0> ztD><*WjTFmB6aTJgCBuP_O8>kNl(8T^r0j-djONge#FnT{pXT(bsvAp@mL+%e0IY- z)jWw#;oOEdbt@6gkJ$}Y5WAr}I|R+m$S(rj4~~iG-U$uvON&n)lJ-Nl^<;h}XE05J z&JHJo1TrpiuTzKf+Tjz8KX^ZaD-23jH!L}lQ(kq5lAN!OF<2|iZ%Iowv&he9PCwwZ zC#-(JS{R%7gpt@YV1ZaDS*%&GbESH1ut{6?v@tPqc6~qqjwMa`E-9D&?SP=l7!n6Z zMHI|;aW!Jxe|Eg{X%;S?uQ5+_@l>pg#NQ#nX~aN$)6tD;E!5dE(f#p@|DaN57j!}E z0Q81h3GmE00C)hi4lXuu*~SBN)0+gn3Z2NhNSz-xkq>&67V$U&>T}-(%zsLhh3aZVXjl|<@21NAK@eY_`oNH^#P1kJ`f{9?V- zR#!~>M(RI|iL1+TcYU>AYS3MM+o?N?&4WSx{g1`!_P0Nel2wILixmv|)Fk)P1n1 zGGM*b=^uDHrfU3{vf_HZQbz8f7g-K-@oP4@)UAc9ifO)`1sq+(XJ6_+F5WjwKBK3) zuBZMb-DCgL8mjz9#s8-@{Hy=_tL6LOY7PI6ML|OPpA~`s3)25T&E#JpAv5qlFjD{W z0{@hh|G<2-_wdBk^mMm2ceZv{aac%vh)+vssoB5vAO5iPjMqiY6@uLW(30m6AE@Lgo4ppnx*$J^oYlZw_Y5xS7GN zGP9k}V;8)8n&&r8Nxm1vKdc$|*||lk^`2lVO2y62kY6L zTAxCfkbzbg+ZBd`st9E_cab=P1*vf{jlZKLiWhu z>1fsLJ_KLS;n6)kN5<2cL6$nfzC-mZuPSp98@Q}$2E?y1)55`%3~k|0?XRFyXlNBY znk7SM+OhzPz2-^i)SCt-F<{Uh9X-yssDyZ~IT+<`=U$o+3}iS*9IGtE6IcfY6^w~h_?&5Qhx9C5PE zTJb9?C%I2Jm&XaX48{O!+EO)FCQxV`1=M@ltIDlTn zE=xql;}jNqLPaTLO6JS!STdpCB(Q4dxwJiAv8W;Jsk_C%T*T)3h&<5r&NQcA1--~U z#5Z3af|#@1!m4zpqjnZbR|vK z6+0jo0ukGY!+faa~+}@Kg^dIxdxSgPc&gg)g9)i zQDwRo{#9;bNY8;$`zKT7nu;q8xKdfw_jZ!A;1wRwxWFlKFY~+P+7mF$=Uc6&lbicOxXA7I`);7_ps^9pe zZKztE=XD!Y%sDE00@iCf+TN*H&Eae=UjJCDUBoIm+e!Okq($G^yGbWXm%p$$I?}`uc^~sO&*zBrknfh+)oMnBD zq_!7+d|m`18Zf6>X|fsBVrZKx2nyH>;B_@2+?eD!IxEaE8=f z0sDo;l9rV0@@w&wuO%>jBkRm~iPXN-pjETHGPk3y`e39UKT4}chzWL0^D-mxnFys& zbu|0VG#aju2M{{SEEk{S{~2m7LrnL=B2#E=U`J!%*SOg|X+7N%o z*ZVE>S1Q$Q99A0+N-u_vc0egq7AG@ouEo53Ip6*#1=obB z7)b^Vu4RCMpgI|BB*Ej#>JWCyVkyK%`?%4i7Ev8E_arP%Y~3a^;3H4NH-wO)x)ee(h2YxVR$i23 zb5WF}**y8iMsFR!k#-o0Ujq<%hHNu})0E1k(jgHza$E^PmmaT{F^S-j?5D0t@4uJs z>zG+rZCRvxHLsga35{5he90rF;foDn2-U1$#ACkK5u~UHSEN?PW0S|jOxiCD2fyQ) z?lBX;G8ej%dUlx$K{gvztZfpF$2Ag`@w8z%0dhW^7_8;xi&yMdjWlK^&al~GdUI>t zGBqDC?Q0Fmt8Og^Kh$m5ylN6lGohhnZ>_~L_94}Gx*jwl=)vw%XhX3YVkFt=5N&2n z63s|6N-9bkm&A+79`F(Py{|LHuSx+=o0Ox~HMlg(VO>>V1!XH4$1tY&m)}G3pbT zK`EDeFRG!+vplcsbqB>K*48e8C}DXxli_vXiS{I8LOBEaY$=_df~gBT;o zd9UyDQf)DlA%baV7%O`-m7g7!?c`NnA3yNAcnxZl;v~)_a9VnDmBhKpeN@fiwANyH zV-K!M3gXZ%8@aoX9wQid8djY_H3mGAG9UQDk@OSW-p&%Ib-eE`|E?!Yn0f`KztTFzLRYcFg00t34Vw`hy| zE8BXa!hOSHWa@EqzajIcG`gK#9K7kNWu?HLZ0}7Xc^#X&M5tk6sT4cL@YRytF7n3# z&E~t1QkSHns2z@DGJ~2OJqz{{hgLWLZ+y~b!0+F^M)@I%weP@UDL|4J3WkiwuJX9KJ&j z(q(Yi8}talH42l#@d~R<#|%po=6u~pFkw)ONGRAcL~_aDK92R+>08oF)+`&(<7X=y z&z10S3HJD5Pm^oYj&J8fy|>>1=d-7R7akD3!R*a}@RsEnp`ZZPK+^am`E^C9HtPrW zC8Anod020dqVf3)9i+pU;L79~BI+1!_rJrBFzzM1Zaw|B2kDX<)I)KO5PbcXvU`Yt z!*K2w3Yms`Dz63_I&Eg^WLU~K*Ik*akW0Ql42i#a7E(a4`%|&wnUTI^4&6Q8|`qfR+^ZTfbUI~C3rb)KP^Q_#WDv* zCUPs=erQr?leJ<`6N(6)f056e4@3XdzQ(xa9ZsV1y$zCnz?%j$|M0mh;?tPcvZbxq z>m)oAb(&k;OZ#-8yFpeP4OY~0t|*S76!21NN;TIanBRm_GnWFxC5zu5*Yd|+P9<$+ z!cTU08#0@;h+MzWi})X4i#%}vJQLu?4k z*Yn4y^!VY;ZvaP?-idDAiNOOyEUw#;nJD_>sWyrB{WggRSd7fum$QW&;bY{x5DbqE z^7THPvq;aCGN%qkR2tltFys|*(1glwci3Bq7tD&yGZ!IzLzGw0*b05TkDs~CB6P^J z2SE$JJV(BS5v^UONyD$jFrd-hAMXSWazV*we|J|sy5BHuwWUFC311N6k*6Q=moN0* zY6@NXxoS5f>YRVWFh6Mn)&1vv0!c zzqdB*Yg&FymNkL17CpCK!rpk@>Bj6n!w+wl&u3Auw>v9+h66Nx+C;}RsnV{P6MXSU zvQkt0O0|_bs4d2yty>9kG=^jM!xoeHTwBi>o@eAeRMzTH5=3K*XUc!9*;dR%_Hk>h z6LN)@e}0Y#JYM1_5fa|{NJ+0ySMKZPv?7GC(C2)gGBy?kHc(E0qvHJ8Z0AYTV1Qf6 uzAHejY-6;4Mos}86#U)OfBY}(!2bVX$A3Rk{qMt$zp(*-vNwMU(*FTx@!9?W literal 0 HcmV?d00001 diff --git a/rootfiles/fonll/README.md b/rootfiles/fonll/README.md new file mode 100644 index 0000000..c8abbbb --- /dev/null +++ b/rootfiles/fonll/README.md @@ -0,0 +1,3 @@ +# FONLL Heavy Quark Production + +Files in this directory contain histograms for dsigma/dpT and dsigma/deta with data obtained from https://www.lpthe.jussieu.fr/~cacciari/fonll/fonllform.html. It uses CTEQ6.6 PDFs except 13.6 TeV which uses NNPDF30nlo_as0118 set. For b quark, it is cross-section to B hadron rather than bare quark while for c quark it uses mixture of 70% of D0 and 30% of D+. The cross-sections are done in the region 0.125 GeV <= pT <= 299.875 and -7.98 <= eta <= 7.98 with 1200 points for each variable. Other parameters are kept at default values. diff --git a/utils/makeFONLLHistograms.py b/utils/makeFONLLHistograms.py new file mode 100644 index 0000000..299e165 --- /dev/null +++ b/utils/makeFONLLHistograms.py @@ -0,0 +1,39 @@ +'''Script to create histograms for FONLL cross-sections used to generate +kinematics. Cross-sections are obtained from +https://www.lpthe.jussieu.fr/~cacciari/fonll/fonllform.html and this script +needs two of them, as a pT and eta which are inputs and script outputs root +file with necessary histograms. Details about individual settings are in +separate README.md in corresponding directory.''' + +import argparse +import numpy +import yaml +import ROOT + +parser = argparse.ArgumentParser( prog='FONLL histogram maker', + description='' ) + +parser.add_argument('-pt') +parser.add_argument('-eta') +parser.add_argument('-out') + +args = parser.parse_args() + +print(args) + +datapt = numpy.loadtxt( args.pt ) +dataeta = numpy.loadtxt( args.eta ) + +outFile = ROOT.TFile.Open( args.out, 'RECREATE' ) + +hheta = ROOT.TH1D( 'eta', '', 1200, -8, 8 ) +hhpt = ROOT.TH1D( 'pT', '', 1200, 0, 300 ) + +for ii in range(1, 1200): + xbin = hheta.FindBin(dataeta[ii-1][0]) + hheta.SetBinContent(xbin, dataeta[ii-1][1]) + xbin = hhpt.FindBin(datapt[ii-1][0]) + hhpt.SetBinContent(xbin, datapt[ii-1][1]) + +outFile.Write() +outFile.Close()