Ticket #993: delphes_card_nodet.tcl

#######################################
# Order of execution of various modules
#######################################

set ExecutionPath {

  NeutrinoFilter

  ElectronFilter
  MuonFilter
  PhotonFilter

  ElectronIsolation
  MuonIsolation
  PhotonIsolation

  MissingET
  FastJetFinder
  UniqueObjectFinder

  TreeWriter
}


#####################
# Neutrino Filter
#####################

module PdgCodeFilter NeutrinoFilter {

  set InputArray Delphes/stableParticles
  set OutputArray filteredParticles

  set PTMin 0.0

  add PdgCode {12}
  add PdgCode {14}
  add PdgCode {16}
  add PdgCode {-12}
  add PdgCode {-14}
  add PdgCode {-16}

}

#################
# Electron filter
#################

module PdgCodeFilter ElectronFilter {
  set InputArray NeutrinoFilter/filteredParticles
  set OutputArray electrons
  set Invert true
  add PdgCode {11}
  add PdgCode {-11}
}

#################
# Muon filter
#################

module PdgCodeFilter MuonFilter {
  set InputArray NeutrinoFilter/filteredParticles
  set OutputArray muons
  set Invert true
  add PdgCode {13}
  add PdgCode {-13}
}

#################
# Photon filter
#################

module PdgCodeFilter PhotonFilter {
  set InputArray NeutrinoFilter/filteredParticles
  set OutputArray photons
  set Invert true
  add PdgCode {22}
}


####################
# Electron isolation
####################

module Isolation ElectronIsolation {
  set CandidateInputArray ElectronFilter/electrons
  set IsolationInputArray NeutrinoFilter/filteredParticles

  set OutputArray electrons

  set DeltaRMax 0.5

  set PTMin 0.5

  set PTRatioMax 0.12
}

################
# Muon isolation
################

module Isolation MuonIsolation {
  set CandidateInputArray MuonFilter/muons
  set IsolationInputArray NeutrinoFilter/filteredParticles

  set OutputArray muons

  set DeltaRMax 0.5

  set PTMin 0.5

  set PTRatioMax 0.25
}

##################
# Photon isolation
##################

module Isolation PhotonIsolation {
  set CandidateInputArray PhotonFilter/photons
  set IsolationInputArray NeutrinoFilter/filteredParticles

  set OutputArray photons

  set DeltaRMax 0.5

  set PTMin 0.5

  set PTRatioMax 0.12
}



#####################
# MC truth jet finder
#####################

module FastJetFinder FastJetFinder {
  set InputArray NeutrinoFilter/filteredParticles

  set OutputArray jets

  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
  set JetAlgorithm 6
  set ParameterR 0.5

  set JetPTMin 20.0
}


#########################
# Gen Missing ET merger
########################

module Merger MissingET {
# add InputArray InputArray
  add InputArray NeutrinoFilter/filteredParticles
  set MomentumOutputArray momentum
}

#####################
# MC truth jet finder
#####################

module FastJetFinder GenJetFinder {
  set InputArray NeutrinoFilter/filteredParticles

  set OutputArray jets

  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
  set JetAlgorithm 6
  set ParameterR 0.5

  set JetPTMin 20.0
}


#####################################################
# Find uniquely identified photons/electrons/tau/jets
#####################################################

module UniqueObjectFinder UniqueObjectFinder {
# earlier arrays take precedence over later ones
# add InputArray InputArray OutputArray
  add InputArray PhotonIsolation/photons photons
  add InputArray ElectronIsolation/electrons electrons
  add InputArray MuonIsolation/muons muons
  add InputArray FastJetFinder/jets jets
}

##################
# ROOT tree writer
##################

# tracks, towers and eflow objects are not stored by default in the output.
# if needed (for jet constituent or other studies), uncomment the relevant
# "add Branch ..." lines.

module TreeWriter TreeWriter {
# add Branch InputArray BranchName BranchClass
  add Branch Delphes/allParticles Particle GenParticle

  add Branch UniqueObjectFinder/jets Jet Jet
  add Branch UniqueObjectFinder/electrons Electron Electron
  add Branch UniqueObjectFinder/photons Photon Photon
  add Branch UniqueObjectFinder/muons Muon Muon
  add Branch MissingET/momentum MissingET MissingET
}