Les Houches

This is an old revision of the document!

---

# The idea is to write a plain text file that can be read easily into
# xml (and hence universally readable).  The name of the block becomes
# the xml tag.  Anything that follows is turned into
# attibutes. e.g. "name ATLAS" inside the block becomes <experiment
# name="ATLAS">.  Anything inside the block starting with "#" becomes
# a comment, i.e. the first block below becomes <analysis
# name="ATLAS" analysis="ATLAS-2014-XXX" cme="8.0" lumi="20.0"
# description="Details about analysis">
# DETECTOROBJ: JET, ELEC, MUON, MET 
# BASE-VARIABLES: PX, PY, PZ, E, PT, ETA, PHI, BETA
# USER-DEFINED VARS: MT(), MT2(), DPHIJMET(J, MET)
 
analysis
  # Details about experiment 
  name ATLAS
  analysis ATLAS-2014-XXX
  publication JHEP11(2014)118
  cme 8.0
  lumi 20.0
end anslysis
 
detectorobj jet
  # detectorobj are lists of objects of same type
  # <detectorobj type="jet" ... >
  # Some human readable definition goes here; we assume sorted by pt
  algorithm anti-kt
  R 0.4
  ptmin 20
  etamax 2.5
  pileup_sub _link_
  # Can provide code and/or link to writeup with details of implementation
  code = _code_link
  doc _link_to_file_with_details	
end detectorobj jet
 
function isolation 2
  # Human readable description: takes 2 arguments; returns (sum of pT in cone of 0.1 around
  # electron)/(pt of lepton) takes in lepton and conesize
  code = _code_link
  doc _link_to_file_with_details	
end function
 
function overlap_removal 2
  # takes in jet list and lepton, calculates based on recipe
  code = _code_link
  doc _link_to_file_with_details	
end function ovelap_removal
 
detectorobj elec_l
  # Can provide code and/or link to writeup with details of implementation
  cut isolation<0.1
  cut overlap_removal=true
  code = _code_link
  doc _link_to_file_with_details	
end detectorobj elec_l
 
detectorobj elec_m elec_l
  # electron-medium inherits from loose.
  # Can provide code and/or link to writeup with details of implementation
  code = _code_link
  doc _link_to_file_with_details	
end detoctorobj elect_t
 
detectorobj lep elec_m mu_m
  # lep inherits from medium electrons or muons
end detectorobj
 
detectorobj met
  # Can provide code and/or link to writeup with details of implementation
  code _code_link
  doc _link_to_file_with_details	
end detectorobj met
 
trigger e_trig
  # Human readable info
  # Each signal analysis will need atleast one trigger to be satisfied
  # Can provide code and/or link to writeup with details of implementation
  eff 0.8
  code _code_link
  doc _link_to_file_with_details	
end trigger e_trig
 
trigger m_trig
  # human readable info
  eff 0.8
  code _code_link
  doc _link_to_file_with_details
end trigger m_trig
 
trigger met_trig
  # human readable info
  eff 0.8	   
  code _code_link
  doc _link_to_file_with_details
end trigger met_trig
 
function pt 2
  # p_T is a function (wrapper) that takes two inputs (1) name of list
  # of detectorobj and (2) the index in that list.  E.g. pt(jet2) =
  # pt(jets,2); pt(elec_tight,1) = pt(elec_t,1) In principle, all
  # common p(), eta(), phi() etc. need not be explicitly defined, this
  # is just an example
  code _code_link
  doc _link_to_file_with_details
end function pt
 
function mt 2
  # mT requires one lepton and met i.e. mt(lep(1),met)
  code _code_link
  doc _link_to_file_with_details
end function mt
 
function mt2 3
  code _code_link
  doc _link_to_file_with_details
end function mt2
 
function b-tag 1
  # takes in a jet and returns whether it is b-tagged
  eff 0.7
  code _code_link
  doc _link_to_file_with_details
end function b-tag  
 
function select 2
  # select(listname, criterion)
  # e.g. select(jet, b-tag) returns number of objects (i.e. jets) give
  # true for the test (i.e. b-tag).  Should be made a global function like pT?
  code _code_link
  doc _link_to_file_with_details
end function select
 
MCsamples
  # Human readable info regarding what programs were used to make what
  # sample suggested format: Signal program extra details (extras are
  # stored in description for xml
  signal Sherpa1.4.1 stop-pair
  background POWHEG ttbar mtop = 172.5 NLO
  doc _link_to_file_with_details
end MCsamples
 
cut cleaning
  # List of preselection cuts: Give human readable descriptions (for
  # your collaborators) and efficiency for theorists
  code _link_to_code
  doc _link_to_doc
end cut
 
cut preselect
  # Pre-selection cuts
  eff 0.95
  trigger trigger_e.OR.trigger_mu.OR.trigger_met
  cut pt(lep(1))>25
  cut size(lep)=1
  cut size(jet)>2
  cut met>100
end preselect
 
detectorobj hadtop
  # Define hadtop 
end detectorobj hadtop
 
function deltaR 2
  # takes in two objects; returns deltaR
end function deltaR
 
detectorobj bjet
  #
  eff 0.7
  code _code_
  doc _doc_link_
end detectorobj bjet
 
function HT
  code _code_
  doc _doc_link_
end function HT
 
function metBYsqrtHT
  # returns MET/sqrt(HT) 
end function metBYsqrtHT
 
function delPhi 2
  # returns minimum delta phi between objects
end function delPhi
 
 
SR tN_diag
  cut preselect
  cut size(jet)>4
  cut pt(jet,1)>60
  cut pt(jet,2)>60
  cut pt(jet,3)>40
  cut pt(jet,4)>25
  cut mt(lep(1),met)>=60
  cut m(hadtop)>=130
  cut m(hadtop)<=205
  cut size(tau)=0
  cut deltaR(bjet(1),lep(1))<2.5
  cut metBYsqrtHT>5
  cut delPhi(jet(1),met)>0.8
  cut delPhi(jet(2),met)>0.8
  hist SHAPE_SR1
end SR
 
result tN_diag
  exp_sig
  stat_sig
  sys_sig
  exp_bg
  stat_bg
  sys_bg
  obs
  s95
  hist SHAPE1
  hist SHAPE2
end result tN_diag