Changes between Version 3 and Version 4 of WorkBook/Tutorials/Hefei

Timestamp:

Nov 20, 2018, 11:42:59 AM (6 years ago)

Author:

Michele Selvaggi

Comment:

--

WorkBook/Tutorials/Hefei

@@ -68 +68 @@
 2) Interactively draw the muon multiplicity and the jet multiplicity. Do you understand these distributions?
 
+{{{
+Delphes->Draw("Muon_size");
+Delphes->Draw("Jet_size");
+}}}
+
 == III) Write a simple analysis macro ==
 
@@ -80 +85 @@
 This Example1.py macro does not produce anything interesting here (it most likely produce an empty plot). The above command is simply shown as an example for how to run a macro. You should open Example1.py with a text editor, and write a small analysis that first selects events with two muons and then reconstructs and plot the recoil Higgs mass using the formula found in III.1)
 
+
+Solution:
+
+{{{
+# Recoil Mass macro
+#!/usr/bin/env python
+
+import sys
+import ROOT
+
+if len(sys.argv) < 2:
+  print " Usage: python examples/MissingMass.py delphes_ee_zh_zmumu.root hist_mrec.root"
+  sys.exit(1)
+
+ROOT.gSystem.Load("libDelphes")
+
+try:
+  ROOT.gInterpreter.Declare('#include "classes/SortableObject.h"')
+  ROOT.gInterpreter.Declare('#include "classes/DelphesClasses.h"')
+  ROOT.gInterpreter.Declare('#include "external/ExRootAnalysis/ExRootTreeReader.h"')
+except:
+  pass
+
+inputFile = sys.argv[1]
+outputFile = sys.argv[2]
+
+# Create chain of root trees
+chain = ROOT.TChain("Delphes")
+chain.Add(inputFile)
+
+# Create object of class ExRootTreeReader
+treeReader = ROOT.ExRootTreeReader(chain)
+numberOfEntries = treeReader.GetEntries()
+
+# Get pointers to branches used in this analysis
+branchMuon = treeReader.UseBranch("Muon")
+
+# Book histograms
+histMass = ROOT.TH1F("mass", "M_{recoil} [GeV]", 60, 100.0, 160.0)
+
+ptot = ROOT.TLorentzVector(0.,0.,0.,240.)
+# Loop over all events
+for entry in range(0, numberOfEntries):
+  # Load selected branches with data from specified event
+  treeReader.ReadEntry(entry)
+
+  # If event contains at least 2 muons
+  if branchMuon.GetEntries() > 1:
+
+    mu1 = branchMuon.At(0)
+    mu2 = branchMuon.At(1)
+
+    pll = mu1.P4() + mu2.P4()
+    ph = ptot - pll
+
+    histMass.Fill(ph.M())
+
+out_root = ROOT.TFile(outputFile,"RECREATE")
+histMass.Write()
+}}}
 
 == IV) Modify the Delphes detector card ==