Changes between Version 10 and Version 11 of WorkBook/Tutorials/Mc4Bsm

Timestamp:

Jul 20, 2016, 9:58:24 AM (8 years ago)

Author:

Michele Selvaggi

Comment:

--

WorkBook/Tutorials/Mc4Bsm

@@ -165 +165 @@
 5. Now [https://cp3.irmp.ucl.ac.be/projects/delphes/raw-attachment/wiki/WorkBook/Tutorials/Mc4Bsm/macro.C download] and run this simple event selection macro. It will plot two crucial jet observables that can help in discriminating between a w-jet and a light jet: the jet mass, and the N-subjettiness ratio tau2/tau1. 
 {{{
+root examples/macro.C'("out_pp_jj_js.root","out_pp_zp_ww_js.root")' 
+}}}
+
+6. Now improve the calorimeter resolution in the card, both energy and angular resolutions (roughly by a factor 20) 
+{{{
+set PhiBins {}
+  for {set i -360} {$i <= 360} {incr i} {
+    add PhiBins [expr {$i * $pi/360.0}]
+  }
+
+  # 0.01 unit in eta up to eta = 2.5
+  for {set i -1000} {$i <= 1000} {incr i} {
+    set eta [expr {$i * 0.005}]
+    add EtaPhiBins $eta $PhiBins
+  }
+
+
+ ...
+
+  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
+  # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
+  set ECalResolutionFormula {                  (abs(eta) <= 1.5) *0.01*(1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
+                             (abs(eta) > 1.5 && abs(eta) <= 2.5) *0.05* (2.16 + 5.6*(abs(eta)-2)^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
+                             (abs(eta) > 2.5 && abs(eta) <= 5.0) *0.05* sqrt(energy^2*0.107^2 + energy*2.08^2)}
+
+
+}}}
+
+7. Re-run Delphes simulation and macro.C with the new configuration and appreciate the effect of improved resolution on the final distributions:
+{{{
+./DelphesSTDHEP cards/delphes_card_CMS.tcl out_pp_zp_ww_js_nc.root pp_zp_ww.hep
+./DelphesSTDHEP cards/delphes_card_CMS.tcl out_pp_jj_js_nc.root pp_jj.hep
 root examples/macro.C'("out_pp_jj_js_bc.root","out_pp_zp_ww_js_bc.root")' 
 }}}