Les Houches

This is an old revision of the document!

---

Part I: Preamble (define detector objects, variables)

The idea is to write a plain text file that can be read easily into xml if needed (and hence universally readable). The name of the block become sthe xml tag. Anything that follows is turned into attributes. e.g. “name ATLAS” inside the block becomes <experiment name=“ATLAS”>. Anything inside the block starting with “#” is added as a comment, i.e. the first block below automatically converts to

<analysis name="ATLAS-2014-XXX" publication="JHEP11(2014)118" cme="8.0" lumi="20.0" description="Details about analysis">

There are five types of pre-defined blocks – detectorobj (defines reconstructed objects), function (defines functions that act on these), cut (blocks of cuts), SR (contains all cuts for a particular signal region) and results (for giving numbers). Besides this, the user can define whatever they want – e.g. here, “analysis” gives extra info, and “table” gives a way to provide tables or digitised histograms.

For any block, a link to code or external documentation can be provided. We can also include a MCsamples block to document the programs used to generate signal and background samples.

When a block inherits the properties of other blocks (e.g. medium electron inherits from loose), the names of those can be added after the name of the current block.

detectorobj elec_m elec_l

analysis
  # Details about experiment 
  name 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
 
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
 
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

PART 2: uses the objects defined above to define the signal regions

cut cleaning
  # List of preselection cuts: Give human readable descriptions (for
  # your experimental collaborators) and efficiency for theorists (if applicable)
  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
 
 
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
end SR

Part 3: Present results of the analysis. Give numbers, exclusions, digitised histograms, link to plots, etc.

result tN_diag
  exp_sig
  stat_sig
  sys_sig
  exp_bg
  stat_bg
  sys_bg
  obs
  s95
  hist data1
end result tN_diag
 
table data1
  # Not real data in example; could be used to provide histograms/tables/likelihoods
  bins nn
  data xlabel ylabel yerr 
   1.0 2.0 1.0
   2.0 1.0 1.0
   3.0 10.0 1.0
   4.0 15.0 1.0
end table data1