Changes between Version 4 and Version 5 of WorkBook/Tutorials/Mc4Bsm

Timestamp:

Jul 20, 2016, 5:42:07 AM (8 years ago)

Author:

Michele Selvaggi

Comment:

--

WorkBook/Tutorials/Mc4Bsm

@@ -76 +76 @@
 To run on our your input file, type:
 {{{
-./DelphesSTDHEP cards/delphes_card_CMS.tcl delphes_output_pp_zp_ww.root pp_zp_ww.hep
+./DelphesSTDHEP cards/delphes_card_CMS.tcl out_pp_zp_ww.root pp_zp_ww.hep
 }}}
 5) Open freshly produced Delphes output with ROOT, and explore it.
 {{{
-root -l  delphes_output_pp_zp_ww.root
+root -l  out_pp_zp_ww.root
 TBrowser t;
 }}}
-In the browser, double click on the " delphes_output_pp_zp_ww.root", and then on the "Delphes" tree. Play around by double clicking on the various branches/observables.
+In the browser, double click on the "out_pp_zp_ww.root", and then on the "Delphes" tree. Play around by double clicking on the various branches/observables.
 
 You can then play plot important observable with a simple selection with the following syntax:
@@ -104 +104 @@
 
 http://cp3.irmp.ucl.ac.be/downloads/RootTreeDescription.html
+
+=== Part II - Understand and modify the configuration file ===
+
+The delphes card can be schematically divided in three parts:
+The !ExecutionPath is where the simulation/reconstruction sequence of modules is defined
+The list of modules configurations.
+The TreeWriter, where the user defines which objects he stores in the output tree.
+
+You can find an explanation for most Delphes modules here:
+https://cp3.irmp.ucl.ac.be/projects/delphes/wiki/WorkBook/Modules
+
+1. Open the card cards/delphes_card_CMS.tcl with your favorite editor and try to make sense of it. 2. 
+
+3. The present version of the card produces a Jet collection without computing and storing in the output substructure variable. Add a new collection of fat-jets called FatJet in the card. In order to do this you will need to first call the process in the ExecutionPath:
+{{{
+set ExecutionPath {
+ ...
+ FastJetFinder
+ FatJetFinder
+ ...  
+}
+}}}
+4. Then configure the FastJetFinder module by switching on the options for substructure:
+{{{
+############
+# Fat Jet finder
+############
+
+module FastJetFinder FatJetFinder {
+  set InputArray EFlowMerger/eflow
+
+  set OutputArray jets
+
+  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 
+  set JetAlgorithm 5 
+  set ParameterR 0.8
+
+  set ComputeNsubjettiness 1
+  set Beta 1.0
+  set AxisMode 4
+
+  set ComputeTrimming 1
+  set RTrim 0.2
+  set PtFracTrim 0.05
+
+  set ComputePruning 1
+  set ZcutPrun 0.1
+  set RcutPrun 0.5
+  set RPrun 0.8
+
+  set ComputeSoftDrop 1
+  set BetaSoftDrop 0.0
+  set SymmetryCutSoftDrop 0.1
+  set R0SoftDrop 0.8
+
+  set JetPTMin 200.0
+}
+}}}
+5. Finally, write the fat-jet collection in your output tree by adding it to the !TreeWriter:
+
+{{{
+module TreeWriter TreeWriter {
+  ...
+  add Branch JetEnergyScale/jets Jet Jet
+  add Branch FatJetFinder/jets FatJet Jet
+  ...
+}
+}}}
+
+6. Now run Delphes with the newly modified card with jet substructure, on both the dijet and Z' samples:
+{{{
+./DelphesSTDHEP cards/delphes_card_CMS.tcl out_pp_zp_ww_js.root pp_zp_ww.hep
+./DelphesSTDHEP cards/delphes_card_CMS.tcl out_pp_jj_js.root pp_jj.hep
+}}}
+
+7. Open the files with ROOT TBrowser as in the Getting Started section, and make sure that the new jet substructure variables are properly stored.