Changes between Version 19 and Version 20 of WorkBook/Tutorials/Mc4Bsm

Timestamp:

Jul 21, 2016, 4:29:38 AM (8 years ago)

Author:

Michele Selvaggi

Comment:

--

WorkBook/Tutorials/Mc4Bsm

@@ -169 +169 @@
 }}}
 
-6. Now improve the calorimeter resolution in the card, both energy and angular resolutions (by a factor 20) 
-{{{
+6. Now create a copy of the existing CMS card:
+{{{
+cp cards/delphes_card_CMS.tcl cards/delphes_card_CMS_nc.tcl
+}}}
+7. Now edit "cards/delphes_card_CMS_nc.tcl" edit the copy to improve the calorimeter resolution, both energy and angular resolutions (by a factor 20) :
+{{{
+#############
+# Calorimeter
+#############
+
+...
+
 set PhiBins {}
   for {set i -360} {$i <= 360} {incr i} {
@@ -182 +192 @@
   }
 
-
  ...
 
-  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
-  # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
+ # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
   set ECalResolutionFormula {                  (abs(eta) <= 1.5) *0.05*(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)}
 
+# set HCalResolutionFormula {resolution formula as a function of eta and energy}
+  set HCalResolutionFormula {                  (abs(eta) <= 3.0) *0.05 sqrt(energy^2*0.050^2 + energy*1.50^2) +
+                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * 0.05* sqrt(energy^2*0.130^2 + energy*2.70^2)}
 
 }}}