Changes between Version 1 and Version 2 of WorkBook/Tutorials/Prefit

Timestamp:

Mar 5, 2020, 2:14:30 PM (5 years ago)

Author:

Michele Selvaggi

Comment:

--

WorkBook/Tutorials/Prefit

@@ -9 +9 @@
 https://twiki.cern.ch/twiki/bin/view/VBSCan/PREFIT20
 
+If you have successfully produced Pythia8 event files in the MadGraph and Parton Shower lecture (Admir Greljo, Ilaria Brivio), copy (or move) the pp->ll (SM and BSM) into the Delphes directory Delphes-3.4.2
+
+FIXME: find out how these files where named during tutorial
+
+If not you can download the files from here (but it might take a while):
+
+{{{
+cd Delphes-3.4.2
+curl -o pp_ll_sm.hepmc  https://cernbox.cern.ch/index.php/s/4rnlQ0bybBYvLTS/download
+curl -o pp_ll_bsm.hepmc  https://cernbox.cern.ch/index.php/s/DfDDS5QixZaBsYO/download
+}}}
 
 
-- gcc/tcl:
 
-For linux users gcc/tcl should be already installed. For Mac users you should install XCode.
+== I) Produce Delphes root files ==
 
-- ROOT:
+In this section you will run Delphes Fast Detector simulation using the default CMS card. Delphes will produce a ROOT file starting from the hepmc format produced by MG5+Pythia8.
 
-can be downloaded from https://root.cern.ch/downloading-root
-Go on latest release, and download a version under "Binary distributions".
+./DelphesHepMC cards/delphes_card_CMS.tcl pp_ll_sm.root pp_ll_sm.hepmc
+./DelphesHepMC cards/delphes_card_CMS.tcl pp_ll_bsm.root pp_ll_bsm.hepmc
 
-- Pythia8:
 
-following instructions from here (or using the Pythia8 installation in !MadGraph):
-
-https://cp3.irmp.ucl.ac.be/projects/delphes/wiki/WorkBook/Pythia8
-
-The solutions for all the exercises can be found in the attachment file (suggestion: download file locally). 
-
-== I) Event generation with Pythia8 + Delphes sample ==
-
-This exercise will teach how to configure the Pythia8 event generator for a simple production of e+e- -> ZH events. Next, you will generate events and simulate the detector with the DelphesPythia8 executable.
-
-1) Create a Pythia8 configuration card that generates N=10k events of ee->Zh->mumu at sqrt(s)=240 GeV (call it "examples/Pythia8/config_ee_zh_zmumu.cmd").
-
-{{{
-Main:numberOfEvents = 10000         ! number of events to generate
-
-Beams:idA = 11                   ! first beam, e- = -11
-Beams:idB = -11                  ! second beam, e+ = 11
-Beams:eCM = 240.                 ! CM energy of collision
-
-! Higgsstrahlung process
-HiggsSM:ffbar2HZ = on
-
-! 5) Force the Z decays to muons
-23:onMode = off
-23:onIfAny = 13 -13
-}}}
-
-2) Produce Delphes events using the above Pythia8 configuration (this command should run Pythia and Delphes on the fly!), using the CEPC detector card "cards/delphes_card_CEPC.tcl"
-
-{{{
-./DelphesPythia8 cards/delphes_card_CEPC.tcl examples/Pythia8/config_ee_zh_zmumu.cmd delphes_ee_zh_zmumu.root
-}}}
-
-== II) Simple Tree analysis ==
+== II) Simple Interactive Tree analysis ==
 
 
@@ -61 +37 @@
 
 {{{
-root -l delphes_ee_zh_zmumu.root
+root -l pp_ll_sm.root
 gSystem->Load("libDelphes");
 TBrowser t;
@@ -67 +43 @@
 
 
-Note: Most objects are described in terms of pp specific variables (PT, Eta, Phi).
-This is simply for historical reasons (and makes of course no difference whatsoever) since Delphes was developed originally as a tool for LHC physics. To plot ee-like variables, one needs to write the translation (or make use of the very useful TLorentzVector of ROOT, see part III).
-
-2) Interactively draw the muon multiplicity and the jet multiplicity. You first have to double-click on the root file icon in the TBrowser. Do you understand these distributions?
+2) Interactively draw the muon and electron multiplicity. You first have to double-click on the root file icon in the TBrowser. Do you understand these distributions?
 
 {{{
 Delphes->Draw("Muon_size");
-Delphes->Draw("Jet_size");
+Delphes->Draw("Electron_size");
 }}}
+
+3) Compare the leading muon transverse momentum and pseudo-rapidity in the SM and the BSM sample.
+
+
+
+{{{
+Delphes->Draw("Muon_size");
+Delphes->Draw("Electron_size");
+}}}
+
+
 
 == III) Write a simple analysis macro ==