Changes between Version 3 and Version 4 of WorkBook/Tutorials/Mc4Bsm
- Timestamp:
- Jul 20, 2016, 5:26:20 AM (8 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
WorkBook/Tutorials/Mc4Bsm
v3 v4 56 56 make -j 4 57 57 }}} 58 3) Download Z' to WWevents in stdhep format58 3) Download Z' (m= 2 TeV) to WW and dijet (pT > 1 TeV) events in stdhep format 59 59 {{{ 60 gunzip ../MG5_aMC_v2_4_0/pp_h_eemm_res/Events/run_01/tag_1_pythia_events.hep.gz 61 gunzip ../MG5_aMC_v2_4_0/pp_h_eemm_res/Events/run_02/tag_1_pythia_events.hep.gz 62 }}} 63 {{{ 64 mkdir -p input 65 mv ../MG5_aMC_v2_4_0/pp_h_eemm_res/Events/run_01/tag_1_pythia_events.hep input/pp_h_eemm_m125.hep 66 mv ../MG5_aMC_v2_4_0/pp_h_eemm_res/Events/run_02/tag_1_pythia_events.hep input/pp_h_eemm_m750.hep 60 wget 61 wget 62 gunzip pp_zp_ww.hep.gz 63 gunzip pp_jj.hep.gz 67 64 }}} 68 65 … … 79 76 To run on our your input file, type: 80 77 {{{ 81 ./DelphesSTDHEP cards/delphes_card_CMS.tcl delphes_output_ m125.root input/pp_h_eemm_m125.hep78 ./DelphesSTDHEP cards/delphes_card_CMS.tcl delphes_output_pp_zp_ww.root pp_zp_ww.hep 82 79 }}} 83 80 5) Open freshly produced Delphes output with ROOT, and explore it. 84 81 {{{ 85 root -l delphes_output_m125.root82 root -l delphes_output_pp_zp_ww.root 86 83 TBrowser t; 87 84 }}} 88 In the browser, double click on the " delphes_output_m125.root", and then on the "Delphes" tree. Play around by double clicking on the various branches/observables.85 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. 89 86 90 87 You can then play plot important observable with a simple selection with the following syntax: … … 96 93 - Note 2: !Muon/Electron - branch name; PT - variable (leaf) of this branch 97 94 {{{ 98 Delphes->Draw("Muon.Eta", "Muon.PT > 20"); 99 Delphes->Draw("Electron.Eta", "Electron.PT > 20"); 100 Delphes->Draw("Jet_size","Electron_size > 1 && Muon_size > 1"); 95 Delphes->Draw("Jet.Mass","Electron_size + Muon_size == 1"); 96 Delphes->Draw("Jet.Mass","Electron_size + Muon_size == 1 && Jet.PT > 500"); 101 97 }}} 102 Objects are already ordered in PT, you can then plot leading ele/muin this way:98 Objects are already ordered in PT, you can then plot the leading jet observables in this way: 103 99 {{{ 104 Delphes->Draw(" Muon[0].PT", "Muon.PT > 20");105 Delphes->Draw(" Electron[0].PT", "Electron.PT > 20");100 Delphes->Draw("Jet[0].Mass"); 101 Delphes->Draw("Jet[0].PT"); 106 102 }}} 107 103 For more information on ROOT trees: 108 104 109 105 http://cp3.irmp.ucl.ac.be/downloads/RootTreeDescription.html 110