####################################### # CLICdet_fcal model for 3TeV CLIC stage # # Based on CLICdet model # by Ulrike Schnoor ulrike.schnoor@cern.ch # # Forward calorimeters implemented # by Aleksander Filip Zarnecki zarnecki@fuw.edu.pl # based on ILCDelphes model and # full simulation results for CLICdet # described in arXiv:1812.07337 ####################################### ####################################### # Order of execution of various modules ####################################### set ExecutionPath { ParticlePropagator ChargedHadronTrackingEfficiency ElectronTrackingEfficiency MuonTrackingEfficiency ChargedHadronMomentumSmearing ElectronMomentumSmearing MuonMomentumSmearing TrackMerger ECal HCal LumiCalF LumiCalR BeamCalF BeamCalR Calorimeter EFlowMerger LumiTowers BCalTowers LumiMerger BCalMerger LumiEfficiency BCalEfficiency PhotonEfficiency PhotonIsolation ElectronFilter ElectronEfficiency ElectronIsolation ChargedHadronFilter MuonEfficiency MuonIsolation EFlowFilter NeutrinoFilter GenJetFinder FastJetFinderKt FastJetFinderVLC_R05_N2 FastJetFinderVLC_R05_N3 FastJetFinderVLC_R05_N4 FastJetFinderVLC_R05_N5 FastJetFinderVLC_R05_N6 FastJetFinderVLC_R07_N2 FastJetFinderVLC_R07_N3 FastJetFinderVLC_R07_N4 FastJetFinderVLC_R07_N5 FastJetFinderVLC_R07_N6 FastJetFinderVLC_R10_N2 FastJetFinderVLC_R10_N3 FastJetFinderVLC_R10_N4 FastJetFinderVLC_R10_N5 FastJetFinderVLC_R10_N6 FastJetFinderVLC_R12_N2 FastJetFinderVLC_R12_N3 FastJetFinderVLC_R12_N4 FastJetFinderVLC_R12_N5 FastJetFinderVLC_R12_N6 FastJetFinderVLC_R15_N2 FastJetFinderVLC_R15_N3 FastJetFinderVLC_R15_N4 FastJetFinderVLC_R15_N5 FastJetFinderVLC_R15_N6 FastJetFinderVLC_R05_inclusive FastJetFinderVLC_R07_inclusive FastJetFinderVLC_R10_inclusive FastJetFinderVLC_R12_inclusive FastJetFinderVLC_R15_inclusive MissingET GenMissingET JetMomentumSmearing_VLCR05N2 JetMomentumSmearing_VLCR05N3 JetMomentumSmearing_VLCR05N4 JetMomentumSmearing_VLCR05N5 JetMomentumSmearing_VLCR05N6 JetMomentumSmearing_VLCR05_inclusive JetMomentumSmearing_VLCR07N2 JetMomentumSmearing_VLCR07N3 JetMomentumSmearing_VLCR07N4 JetMomentumSmearing_VLCR07N5 JetMomentumSmearing_VLCR07N6 JetMomentumSmearing_VLCR07_inclusive JetMomentumSmearing_VLCR10N2 JetMomentumSmearing_VLCR10N3 JetMomentumSmearing_VLCR10N4 JetMomentumSmearing_VLCR10N5 JetMomentumSmearing_VLCR10N6 JetMomentumSmearing_VLCR10_inclusive JetMomentumSmearing_VLCR12N2 JetMomentumSmearing_VLCR12N3 JetMomentumSmearing_VLCR12N4 JetMomentumSmearing_VLCR12N5 JetMomentumSmearing_VLCR12N6 JetMomentumSmearing_VLCR12_inclusive JetMomentumSmearing_VLCR15N2 JetMomentumSmearing_VLCR15N3 JetMomentumSmearing_VLCR15N4 JetMomentumSmearing_VLCR15N5 JetMomentumSmearing_VLCR15N6 JetMomentumSmearing_VLCR15_inclusive JetFlavorAssociation_R05N2 JetFlavorAssociation_R05N3 JetFlavorAssociation_R05N4 JetFlavorAssociation_R05N5 JetFlavorAssociation_R05N6 JetFlavorAssociation_R07N2 JetFlavorAssociation_R07N3 JetFlavorAssociation_R07N4 JetFlavorAssociation_R07N5 JetFlavorAssociation_R07N6 JetFlavorAssociation_R10N2 JetFlavorAssociation_R10N3 JetFlavorAssociation_R10N4 JetFlavorAssociation_R10N5 JetFlavorAssociation_R10N6 JetFlavorAssociation_R12N2 JetFlavorAssociation_R12N3 JetFlavorAssociation_R12N4 JetFlavorAssociation_R12N5 JetFlavorAssociation_R12N6 JetFlavorAssociation_R15N2 JetFlavorAssociation_R15N3 JetFlavorAssociation_R15N4 JetFlavorAssociation_R15N5 JetFlavorAssociation_R15N6 JetFlavorAssociation_R05_inclusive JetFlavorAssociation_R07_inclusive JetFlavorAssociation_R10_inclusive JetFlavorAssociation_R12_inclusive JetFlavorAssociation_R15_inclusive JetFlavorAssociation_JER_R05N2 JetFlavorAssociation_JER_R05N3 JetFlavorAssociation_JER_R05N4 JetFlavorAssociation_JER_R05N5 JetFlavorAssociation_JER_R05N6 JetFlavorAssociation_JER_R07N2 JetFlavorAssociation_JER_R07N3 JetFlavorAssociation_JER_R07N4 JetFlavorAssociation_JER_R07N5 JetFlavorAssociation_JER_R07N6 JetFlavorAssociation_JER_R10N2 JetFlavorAssociation_JER_R10N3 JetFlavorAssociation_JER_R10N4 JetFlavorAssociation_JER_R10N5 JetFlavorAssociation_JER_R10N6 JetFlavorAssociation_JER_R12N2 JetFlavorAssociation_JER_R12N3 JetFlavorAssociation_JER_R12N4 JetFlavorAssociation_JER_R12N5 JetFlavorAssociation_JER_R12N6 JetFlavorAssociation_JER_R15N2 JetFlavorAssociation_JER_R15N3 JetFlavorAssociation_JER_R15N4 JetFlavorAssociation_JER_R15N5 JetFlavorAssociation_JER_R15N6 JetFlavorAssociation_JER_R05_inclusive JetFlavorAssociation_JER_R07_inclusive JetFlavorAssociation_JER_R10_inclusive JetFlavorAssociation_JER_R12_inclusive JetFlavorAssociation_JER_R15_inclusive BTaggingWP50_R05N2 BTaggingWP70_R05N2 BTaggingWP90_R05N2 BTaggingWP50_R05N3 BTaggingWP70_R05N3 BTaggingWP90_R05N3 BTaggingWP50_R05N4 BTaggingWP70_R05N4 BTaggingWP90_R05N4 BTaggingWP50_R05N5 BTaggingWP70_R05N5 BTaggingWP90_R05N5 BTaggingWP50_R05N6 BTaggingWP70_R05N6 BTaggingWP90_R05N6 BTaggingWP50_R07N2 BTaggingWP70_R07N2 BTaggingWP90_R07N2 BTaggingWP50_R07N3 BTaggingWP70_R07N3 BTaggingWP90_R07N3 BTaggingWP50_R07N4 BTaggingWP70_R07N4 BTaggingWP90_R07N4 BTaggingWP50_R07N5 BTaggingWP70_R07N5 BTaggingWP90_R07N5 BTaggingWP50_R07N6 BTaggingWP70_R07N6 BTaggingWP90_R07N6 BTaggingWP50_R10N2 BTaggingWP70_R10N2 BTaggingWP90_R10N2 BTaggingWP50_R10N3 BTaggingWP70_R10N3 BTaggingWP90_R10N3 BTaggingWP50_R10N4 BTaggingWP70_R10N4 BTaggingWP90_R10N4 BTaggingWP50_R10N5 BTaggingWP70_R10N5 BTaggingWP90_R10N5 BTaggingWP50_R10N6 BTaggingWP70_R10N6 BTaggingWP90_R10N6 BTaggingWP50_R12N2 BTaggingWP70_R12N2 BTaggingWP90_R12N2 BTaggingWP50_R12N3 BTaggingWP70_R12N3 BTaggingWP90_R12N3 BTaggingWP50_R12N4 BTaggingWP70_R12N4 BTaggingWP90_R12N4 BTaggingWP50_R12N5 BTaggingWP70_R12N5 BTaggingWP90_R12N5 BTaggingWP50_R12N6 BTaggingWP70_R12N6 BTaggingWP90_R12N6 BTaggingWP50_R15N2 BTaggingWP70_R15N2 BTaggingWP90_R15N2 BTaggingWP50_R15N3 BTaggingWP70_R15N3 BTaggingWP90_R15N3 BTaggingWP50_R15N4 BTaggingWP70_R15N4 BTaggingWP90_R15N4 BTaggingWP50_R15N5 BTaggingWP70_R15N5 BTaggingWP90_R15N5 BTaggingWP50_R15N6 BTaggingWP70_R15N6 BTaggingWP90_R15N6 BTaggingWP50_R05_inclusive BTaggingWP70_R05_inclusive BTaggingWP90_R05_inclusive BTaggingWP50_R07_inclusive BTaggingWP70_R07_inclusive BTaggingWP90_R07_inclusive BTaggingWP50_R10_inclusive BTaggingWP70_R10_inclusive BTaggingWP90_R10_inclusive BTaggingWP50_R12_inclusive BTaggingWP70_R12_inclusive BTaggingWP90_R12_inclusive BTaggingWP50_R15_inclusive BTaggingWP70_R15_inclusive BTaggingWP90_R15_inclusive BTagging_JER_WP50_R05N2 BTagging_JER_WP70_R05N2 BTagging_JER_WP90_R05N2 BTagging_JER_WP50_R05N3 BTagging_JER_WP70_R05N3 BTagging_JER_WP90_R05N3 BTagging_JER_WP50_R05N4 BTagging_JER_WP70_R05N4 BTagging_JER_WP90_R05N4 BTagging_JER_WP50_R05N5 BTagging_JER_WP70_R05N5 BTagging_JER_WP90_R05N5 BTagging_JER_WP50_R05N6 BTagging_JER_WP70_R05N6 BTagging_JER_WP90_R05N6 BTagging_JER_WP50_R07N2 BTagging_JER_WP70_R07N2 BTagging_JER_WP90_R07N2 BTagging_JER_WP50_R07N3 BTagging_JER_WP70_R07N3 BTagging_JER_WP90_R07N3 BTagging_JER_WP50_R07N4 BTagging_JER_WP70_R07N4 BTagging_JER_WP90_R07N4 BTagging_JER_WP50_R07N5 BTagging_JER_WP70_R07N5 BTagging_JER_WP90_R07N5 BTagging_JER_WP50_R07N6 BTagging_JER_WP70_R07N6 BTagging_JER_WP90_R07N6 BTagging_JER_WP50_R10N2 BTagging_JER_WP70_R10N2 BTagging_JER_WP90_R10N2 BTagging_JER_WP50_R10N3 BTagging_JER_WP70_R10N3 BTagging_JER_WP90_R10N3 BTagging_JER_WP50_R10N4 BTagging_JER_WP70_R10N4 BTagging_JER_WP90_R10N4 BTagging_JER_WP50_R10N5 BTagging_JER_WP70_R10N5 BTagging_JER_WP90_R10N5 BTagging_JER_WP50_R10N6 BTagging_JER_WP70_R10N6 BTagging_JER_WP90_R10N6 BTagging_JER_WP50_R12N2 BTagging_JER_WP70_R12N2 BTagging_JER_WP90_R12N2 BTagging_JER_WP50_R12N3 BTagging_JER_WP70_R12N3 BTagging_JER_WP90_R12N3 BTagging_JER_WP50_R12N4 BTagging_JER_WP70_R12N4 BTagging_JER_WP90_R12N4 BTagging_JER_WP50_R12N5 BTagging_JER_WP70_R12N5 BTagging_JER_WP90_R12N5 BTagging_JER_WP50_R12N6 BTagging_JER_WP70_R12N6 BTagging_JER_WP90_R12N6 BTagging_JER_WP50_R15N2 BTagging_JER_WP70_R15N2 BTagging_JER_WP90_R15N2 BTagging_JER_WP50_R15N3 BTagging_JER_WP70_R15N3 BTagging_JER_WP90_R15N3 BTagging_JER_WP50_R15N4 BTagging_JER_WP70_R15N4 BTagging_JER_WP90_R15N4 BTagging_JER_WP50_R15N5 BTagging_JER_WP70_R15N5 BTagging_JER_WP90_R15N5 BTagging_JER_WP50_R15N6 BTagging_JER_WP70_R15N6 BTagging_JER_WP90_R15N6 BTagging_JER_WP50_R05_inclusive BTagging_JER_WP70_R05_inclusive BTagging_JER_WP90_R05_inclusive BTagging_JER_WP50_R07_inclusive BTagging_JER_WP70_R07_inclusive BTagging_JER_WP90_R07_inclusive BTagging_JER_WP50_R10_inclusive BTagging_JER_WP70_R10_inclusive BTagging_JER_WP90_R10_inclusive BTagging_JER_WP50_R12_inclusive BTagging_JER_WP70_R12_inclusive BTagging_JER_WP90_R12_inclusive BTagging_JER_WP50_R15_inclusive BTagging_JER_WP70_R15_inclusive BTagging_JER_WP90_R15_inclusive TauTagging_R05N2 TauTagging_R05N3 TauTagging_R05N4 TauTagging_R05N5 TauTagging_R05N6 TauTagging_R07N2 TauTagging_R07N3 TauTagging_R07N4 TauTagging_R07N5 TauTagging_R07N6 TauTagging_R10N2 TauTagging_R10N3 TauTagging_R10N4 TauTagging_R10N5 TauTagging_R10N6 TauTagging_R12N2 TauTagging_R12N3 TauTagging_R12N4 TauTagging_R12N5 TauTagging_R12N6 TauTagging_R15N2 TauTagging_R15N3 TauTagging_R15N4 TauTagging_R15N5 TauTagging_R15N6 TauTagging_R05_inclusive TauTagging_R07_inclusive TauTagging_R10_inclusive TauTagging_R12_inclusive TauTagging_R15_inclusive TauTagging_JER_R05N2 TauTagging_JER_R05N3 TauTagging_JER_R05N4 TauTagging_JER_R05N5 TauTagging_JER_R05N6 TauTagging_JER_R07N2 TauTagging_JER_R07N3 TauTagging_JER_R07N4 TauTagging_JER_R07N5 TauTagging_JER_R07N6 TauTagging_JER_R10N2 TauTagging_JER_R10N3 TauTagging_JER_R10N4 TauTagging_JER_R10N5 TauTagging_JER_R10N6 TauTagging_JER_R12N2 TauTagging_JER_R12N3 TauTagging_JER_R12N4 TauTagging_JER_R12N5 TauTagging_JER_R12N6 TauTagging_JER_R15N2 TauTagging_JER_R15N3 TauTagging_JER_R15N4 TauTagging_JER_R15N5 TauTagging_JER_R15N6 TauTagging_JER_R05_inclusive TauTagging_JER_R07_inclusive TauTagging_JER_R10_inclusive TauTagging_JER_R12_inclusive TauTagging_JER_R15_inclusive ScalarHT TreeWriter } ################################# # Propagate particles in cylinder ################################# module ParticlePropagator ParticlePropagator { set InputArray Delphes/stableParticles set OutputArray stableParticles set ChargedHadronOutputArray chargedHadrons set ElectronOutputArray electrons set MuonOutputArray muons # radius of the magnetic field coverage in the calorimeter, in m set Radius 1.5 # half-length of the magnetic field coverage in the calorimeter, in m set HalfLength 2.31 # magnetic field, in T set Bz 4.0 } #################################### # Charged hadron tracking efficiency #################################### module Efficiency ChargedHadronTrackingEfficiency { set InputArray ParticlePropagator/chargedHadrons set OutputArray chargedHadrons # Current full simulation with CLICdet provides for pions: set EfficiencyFormula { (abs(eta) > 2.54) * (0.000) + (energy >= 80) * (abs(eta) < 2.54) * (1.000) + (energy < 80 && energy >= 3) * (abs(eta) <=2.54 && abs(eta) > 2.34) * (0.994) + (energy < 80 && energy >= 3) * (abs(eta) <= 2.34) * (1.000) + (energy < 3) * (abs(eta) <= 2.54 && abs(eta) > 0.55 ) * (0.000) + (energy < 3) * (abs(eta) <= 0.55 ) * (1.000) } } ############################## # Electron tracking efficiency ############################## module Efficiency ElectronTrackingEfficiency { set InputArray ParticlePropagator/electrons set OutputArray electrons # Current full simulation with CLICdet provides for electrons: set EfficiencyFormula { (pt <= 1) * (0.000) + (abs(eta) > 2.54) * (0.000) + (energy >= 80) * (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (0.993) + (energy >= 80) * (abs(eta) <= 2.44 && abs(eta) > 2.34 ) * (0.997) + (energy >= 80) * (abs(eta) <= 2.34 ) * (1.000) + (energy < 80 && energy >= 5) * (abs(eta) <= 2.54 && abs(eta) > 2.17 ) * (0.998) + (energy < 80 && energy >= 5) * (abs(eta) <= 2.17) * (1.000) + (energy < 5) * (abs(eta) > 2.34 ) * (0.000) + (energy < 5) * (abs(eta) <= 2.34 && abs(eta) > 0.76 ) * (0.997) + (energy < 5) * (abs(eta) <= 0.76) * (0.999) } } ########################## # Muon tracking efficiency ########################## module Efficiency MuonTrackingEfficiency { set InputArray ParticlePropagator/muons set OutputArray muons # Current full simulation with CLICdet provides for muons: set EfficiencyFormula { (pt < 1) * (0.000) + (abs(eta) > 2.54) * (0.000) + (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 80) * (0.994) + (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 5 && energy < 80) * (0.996) + (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy < 5 ) * (0.996) + (abs(eta) <= 2.44 ) * (energy >= 5 ) * (1.000) + (abs(eta) <= 2.44 && abs(eta) > 2.25 ) * (energy < 5 && pt >=1 ) * (0.999) + (abs(eta) <= 2.25 ) * (energy >= 1) * (1.000) } } ######################################## # Momentum resolution for charged tracks ######################################## module MomentumSmearing ChargedHadronMomentumSmearing { set InputArray ChargedHadronTrackingEfficiency/chargedHadrons set OutputArray chargedHadrons # Resolution given in dpT/pT. # CLICdet internal studies set ResolutionFormula { (abs(eta) < 2.66 && abs(eta) >= 1.74 ) * 2 * sqrt( 8.56036e-05^2 * pt^2 +0.0148987^2 ) + (abs(eta) < 1.74 && abs(eta) >= 1.32 ) * 1.5 * sqrt( 8.56036e-05^2 * pt^2 +0.0148987^2 ) + (abs(eta) < 1.32 && abs(eta) >= 0.76 ) * sqrt( 1.12382e-05^2 * pt^2 +0.00391722^2 ) + (abs(eta) < 0.76 && abs(eta) >= 0.36 ) * sqrt( 1.16768e-05^2 * pt^2 +0.00255204^2 ) + (abs(eta) < 0.36 && abs(eta) >= 0.18 ) * sqrt( 1.28327e-05^2 * pt^2 +0.00220587^2 ) + (abs(eta) < 0.18) * sqrt( 1.32845e-05^2 * pt^2 +0.00209325^2 ) } } ################################### # Momentum resolution for electrons ################################### module MomentumSmearing ElectronMomentumSmearing { set InputArray ElectronTrackingEfficiency/electrons set OutputArray electrons # Resolution given in dpT/pT. # CLICdet internal studies set ResolutionFormula { (abs(eta) < 2.66 && abs(eta) >= 1.74 ) * 2 * sqrt( 8.62283e-05^2 * pt^2 + 0.0177556^2 ) + (abs(eta) < 1.74 && abs(eta) >= 1.32 ) * 1.5 * sqrt( 8.62283e-05^2 * pt^2 + 0.0177556^2 ) + (abs(eta) < 1.32 && abs(eta) >= 0.76 ) * sqrt( 1.0915e-05 ^2 * pt^2 + 0.00663766^2 ) + (abs(eta) < 0.76 && abs(eta) >= 0.36 ) * sqrt( 1.15518e-05^2 * pt^2 + 0.00398644^2 ) + (abs(eta) < 0.36 && abs(eta) >= 0.18 ) * sqrt( 1.3307e-05 ^2 * pt^2 + 0.00317807^2 ) + (abs(eta) < 0.18) * sqrt( 1.40722e-05^2 * pt^2 + 0.00292138^2 ) } } ############################### # Momentum resolution for muons ############################### module MomentumSmearing MuonMomentumSmearing { set InputArray MuonTrackingEfficiency/muons set OutputArray muons # Resolution given in dpT/pT. # CLICdet internal studies set ResolutionFormula { (abs(eta) < 2.66 && abs(eta) >= 1.74 ) * 2 * sqrt(4.57439e-05^2 * pt^2 + 0.0149328^2 ) + (abs(eta) < 1.74 && abs(eta) >= 1.32 ) * 1.5 * sqrt(4.57439e-05^2 * pt^2 + 0.0149328^2 ) + (abs(eta) < 1.32 && abs(eta) >= 0.76 ) * sqrt(9.81626e-06^2 * pt^2 + 0.00379895^2 ) + (abs(eta) < 0.76 && abs(eta) >= 0.36 ) * sqrt(1.1959e-05^2 * pt^2 + 0.00242417^2 ) + (abs(eta) < 0.36 && abs(eta) >= 0.18 ) * sqrt(1.20149e-05^2 * pt^2 + 0.00219291^2 ) + (abs(eta) < 0.18) * sqrt(1.29686e-05^2 * pt^2 + 0.0020392^2 ) } } ############## # Track merger ############## module Merger TrackMerger { # add InputArray InputArray add InputArray ChargedHadronMomentumSmearing/chargedHadrons add InputArray ElectronMomentumSmearing/electrons add InputArray MuonMomentumSmearing/muons set OutputArray tracks } ############# # ECAL ############# module SimpleCalorimeter ECal { set ParticleInputArray ParticlePropagator/stableParticles set TrackInputArray TrackMerger/tracks set TowerOutputArray ecalTowers set EFlowTrackOutputArray eflowTracks set EFlowTowerOutputArray eflowPhotons set IsEcal true set EnergyMin 0.5 set EnergySignificanceMin 1.0 set SmearTowerCenter true set pi [expr {acos(-1)}] # lists of the edges of each tower in eta and phi # each list starts with the lower edge of the first tower # the list ends with the higher edged of the last tower #ECAL barrel: dphi = 0.2 degree, deta=0.003 towers up to |eta| <=1.2 #ECAL endcaps: dphi = 0.8 degree, deta=0.02 towers up to |eta| <=2.5 #ECAL plug: dphi = 1 degree, deta = 0.02 up to |eta| <=3 #ECAL cell sizes always 5x5 mm^2 #barrel: #dphi = 0.2 degree towers up to eta <=1.2 set PhiBins {} for {set i -900} {$i <= 900} {incr i} { add PhiBins [expr {$i * $pi/900.0 }] } # 0.003 unit (5x5 mm^2) in eta up to eta <=1.2 for {set i -400} {$i <=400} {incr i} { set eta [expr {$i * 0.003}] add EtaPhiBins $eta $PhiBins } #endcaps: #dphi = 0.8 degree towers for 1.2 < eta <=2.5 set PhiBins {} for {set i -225} {$i <= 225} {incr i} { add PhiBins [expr {$i * $pi/225.}] } #deta=0.02 units for 1.2 < |eta| <=2.5 #first, from -2.5 to -1.2, there will be (1.3/0.02=)65 segments for {set i 1} {$i <=66} {incr i} { set eta [expr {-2.52 + $i * 0.02}] add EtaPhiBins $eta $PhiBins } #same for 1.2 to 2.5 for {set i 1} {$i <=66} {incr i} { set eta [expr {1.18 + $i*0.02}] add EtaPhiBins $eta $PhiBins } #plug: #dphi = 1 degree for 2.5 < eta <=3 set PhiBins {} for {set i -180} {$i <= 180} {incr i} { add PhiBins [expr {$i * $pi/180.}] } # deta = 0.02 for 2.5 < |eta| <=3 # from -3 to -2.5, there will be 25 segments for {set i 1} {$i <= 26} {incr i} { set eta [expr {-3.02 + $i * 0.02}] add EtaPhiBins $eta $PhiBins } #same for 2.5 to 3 for {set i 1} {$i <= 26} {incr i} { set eta [expr {2.48 + $i*0.02}] add EtaPhiBins $eta $PhiBins } # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL} add EnergyFraction {0} {0.0} # energy fractions for e, gamma and pi0 add EnergyFraction {11} {1.0} add EnergyFraction {22} {1.0} add EnergyFraction {111} {1.0} # energy fractions for muon, neutrinos and neutralinos add EnergyFraction {12} {0.0} add EnergyFraction {13} {0.0} add EnergyFraction {14} {0.0} add EnergyFraction {16} {0.0} add EnergyFraction {1000022} {0.0} add EnergyFraction {1000023} {0.0} add EnergyFraction {1000025} {0.0} add EnergyFraction {1000035} {0.0} add EnergyFraction {1000045} {0.0} # energy fractions for K0short and Lambda add EnergyFraction {310} {0.3} add EnergyFraction {3122} {0.3} # set ECalResolutionFormula {resolution formula as a function of eta and energy} set ResolutionFormula { (abs(eta) <= 0.78 ) * sqrt(energy^2*0.01^2 + energy*0.156^2)+ (abs(eta) > 0.78 && abs(eta) <=0.83 ) * sqrt( energy^0.01^2 + energy*0.175^2 ) + (abs(eta) <= 3 && abs(eta) > 0.83) * sqrt( energy^2*0.01^2 + energy*0.151^2 )} } ############# # HCAL ############# module SimpleCalorimeter HCal { set ParticleInputArray ParticlePropagator/stableParticles set TrackInputArray ECal/eflowTracks set TowerOutputArray hcalTowers set EFlowTrackOutputArray eflowTracks set EFlowTowerOutputArray eflowNeutralHadrons set IsEcal false set EnergyMin 1.0 set EnergySignificanceMin 1.0 set SmearTowerCenter true set pi [expr {acos(-1)}] # lists of the edges of each tower in eta and phi # each list starts with the lower edge of the first tower # the list ends with the higher edged of the last tower #HCAL barrel: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.8 #HCAL ring: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.9 #HCAL endcaps: dphi = 6 degree, deta = 0.1 up to |eta| <=3.5 #HCAL cell sizes always 30x30 mm^2 #barrel and ring: #dphi = 1 degree up to |eta| <=0.9 set PhiBins {} for {set i -180} {$i <=180} {incr i} { add PhiBins [expr {$i * $pi/180.0}] } #deta= 0.02 towers up to |eta| <=0.9 for {set i -45} {$i <=45} {incr i} { set eta [expr {$i * 0.02}] add EtaPhiBins $eta $PhiBins } #endcaps: # dphi = 6 degree set PhiBins {} for {set i -30} {$i <=30} {incr i} { add PhiBins [expr {$i * $pi/30.0}] } # deta =0.1 for 0.9 < |eta| <=3.0 #for -3.0 to -0.9, 21 segments for {set i 1} {$i <=22} {incr i} { set eta [expr {-3.1 + $i * 0.1}] add EtaPhiBins $eta $PhiBins } #same for 0.9 to 3.0 for {set i 1} {$i <=22} {incr i} { set eta [expr {0.8 + $i * 0.1 }] add EtaPhiBins $eta $PhiBins } # default energy fractions {abs(PDG code)} {Fecal Fhcal} add EnergyFraction {0} {1.0} # energy fractions for e, gamma and pi0 add EnergyFraction {11} {0.0} add EnergyFraction {22} {0.0} add EnergyFraction {111} {0.0} # energy fractions for muon, neutrinos and neutralinos add EnergyFraction {12} {0.0} add EnergyFraction {13} {0.0} add EnergyFraction {14} {0.0} add EnergyFraction {16} {0.0} add EnergyFraction {1000022} {0.0} add EnergyFraction {1000023} {0.0} add EnergyFraction {1000025} {0.0} add EnergyFraction {1000035} {0.0} add EnergyFraction {1000045} {0.0} # energy fractions for K0short and Lambda add EnergyFraction {310} {0.7} add EnergyFraction {3122} {0.7} # set HCalResolutionFormula {resolution formula as a function of eta and energy} #CLICdet internal studies set ResolutionFormula { (abs(eta)<= 0.3) * sqrt(1.38^2 + energy*0.308^2 + energy^2*0.050^2) + (abs(eta)<= 0.78 && abs(eta) > 0.3) * sqrt(1.25^2 + energy*0.322^2 + energy^2*0.048^2) + (abs(eta)<=1.099 && abs(eta) > 0.78) * sqrt( 1.159^2 + energy*0.341^2 + energy^2*0.049^2 ) + (abs(eta)<=3 && abs(eta)> 1.099) * sqrt( 1.09^2 + energy*0.319^2 + energy^2*0.052^2 ) } } ############## # LumiCal ############## module SimpleCalorimeter LumiCalF { set ParticleInputArray ParticlePropagator/stableParticles set TrackInputArray TrackMerger/tracks set TowerOutputArray lumicalTowers set EFlowTowerOutputArray lumicalPhotons set IsEcal true set EnergyMin 2.0 set EnergySignificanceMin 1.0 set SmearTowerCenter true source CLIC/CLICdet_LumiCalF_Binning.tcl source CLIC/CLICdet_LumiCal_EnergyFractions.tcl source CLIC/CLICdet_LumiCal_Resolution_Stage3.tcl } module SimpleCalorimeter LumiCalR { set ParticleInputArray ParticlePropagator/stableParticles set TrackInputArray TrackMerger/tracks set TowerOutputArray lumicalTowers set EFlowTowerOutputArray lumicalPhotons set IsEcal true set EnergyMin 2.0 set EnergySignificanceMin 1.0 set SmearTowerCenter true source CLIC/CLICdet_LumiCalR_Binning.tcl source CLIC/CLICdet_LumiCal_EnergyFractions.tcl source CLIC/CLICdet_LumiCal_Resolution_Stage3.tcl } ############## # BeamCal ############## module SimpleCalorimeter BeamCalR { set ParticleInputArray ParticlePropagator/stableParticles set TrackInputArray TrackMerger/tracks set TowerOutputArray bcalTowers set EFlowTowerOutputArray bcalPhotons set IsEcal true set EnergyMin 5.0 set EnergySignificanceMin 1.0 set SmearTowerCenter true source CLIC/CLICdet_BeamCalR_Binning.tcl source CLIC/CLICdet_BeamCal_EnergyFractions.tcl source CLIC/CLICdet_BeamCal_Resolution_Stage3.tcl } module SimpleCalorimeter BeamCalF { set ParticleInputArray ParticlePropagator/stableParticles set TrackInputArray TrackMerger/tracks set TowerOutputArray bcalTowers set EFlowTowerOutputArray bcalPhotons set IsEcal true set EnergyMin 5.0 set EnergySignificanceMin 1.0 set SmearTowerCenter true source CLIC/CLICdet_BeamCalF_Binning.tcl source CLIC/CLICdet_BeamCal_EnergyFractions.tcl source CLIC/CLICdet_BeamCal_Resolution_Stage3.tcl } ################# # Electron filter ################# module PdgCodeFilter ElectronFilter { set InputArray HCal/eflowTracks set OutputArray electrons set Invert true add PdgCode {11} add PdgCode {-11} } ###################### # ChargedHadronFilter ###################### module PdgCodeFilter ChargedHadronFilter { set InputArray HCal/eflowTracks set OutputArray chargedHadrons add PdgCode {11} add PdgCode {-11} add PdgCode {13} add PdgCode {-13} } ################################################### # Tower Merger (in case not using e-flow algorithm) ################################################### module Merger Calorimeter { # add InputArray InputArray add InputArray ECal/ecalTowers add InputArray HCal/hcalTowers set OutputArray towers } ############################### # LumiCal tower merger ############################### module Merger LumiTowers { # add InputArray InputArray add InputArray LumiCalF/lumicalTowers add InputArray LumiCalR/lumicalTowers set OutputArray lumicalTowers } ############################### # BeamCal tower merger ############################### module Merger BCalTowers { # add InputArray InputArray add InputArray BeamCalF/bcalTowers add InputArray BeamCalR/bcalTowers set OutputArray bcalTowers } #################### # Energy flow merger #################### module Merger EFlowMerger { # add InputArray InputArray add InputArray HCal/eflowTracks add InputArray ECal/eflowPhotons add InputArray HCal/eflowNeutralHadrons set OutputArray eflow } ############################### # LumiCal energy flow merger ############################### module Merger LumiMerger { # add InputArray InputArray add InputArray LumiCalF/lumicalPhotons add InputArray LumiCalR/lumicalPhotons set OutputArray lumicalPhotons } ############################### # BeamCal energy flow merger ############################### module Merger BCalMerger { # add InputArray InputArray add InputArray BeamCalF/bcalPhotons add InputArray BeamCalR/bcalPhotons set OutputArray bcalPhotons } ################### # Photon efficiency ################### module Efficiency PhotonEfficiency { set InputArray ECal/eflowPhotons set OutputArray photons # set EfficiencyFormula {efficiency formula as a function of eta and pt} # efficiency formula for photons # current full simulation of CLICdet yields: set EfficiencyFormula { (energy < 2.0 ) * (0.000) + (energy >= 2.0) * (abs(eta) < 0.7)*(0.94) + (energy >= 2.0) * (abs(eta) >=0.7 && abs(eta) <=3.0) * (0.9) } } ################## # Photon isolation ################## module Isolation PhotonIsolation { set CandidateInputArray PhotonEfficiency/photons set IsolationInputArray EFlowMerger/eflow set OutputArray photons set DeltaRMax 0.1 set PTMin 0.5 set PTRatioMax 0.2 } ############################## # LumiCal photon efficiency ############################## module Efficiency LumiEfficiency { set InputArray LumiMerger/lumicalPhotons set OutputArray lumicalPhotons source CLIC/CLICdet_LumiCal_Efficiency_Stage3.tcl } ############################## # BeamCal photon efficiency ############################## module Efficiency BCalEfficiency { set InputArray BCalMerger/bcalPhotons set OutputArray bcalPhotons source CLIC/CLICdet_BeamCal_Efficiency_Stage3.tcl } ##################### # Electron efficiency ##################### module Efficiency ElectronEfficiency { set InputArray ElectronFilter/electrons set OutputArray electrons # set EfficiencyFormula {efficiency formula as a function of eta and pt} set EfficiencyFormula { (energy < 3.0) * ( 0.00 ) + ( energy >=3 && energy < 8 ) * (abs(eta) > 1.95) * (0.58 ) + ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.7 ) + ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.6 ) + ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.7 ) + ( energy >=3 && energy < 8 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.8 ) + ( energy >=3 && energy < 8 ) * (abs(eta) <= 0.69) * (0.84 ) + ( energy >=8 && energy < 13 ) * (abs(eta) > 1.95) * ( 0.6 ) + ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.76 ) + ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.67 ) + ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.78 ) + ( energy >=8 && energy < 13 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.86 ) + ( energy >=8 && energy < 13 ) * (abs(eta) <= 0.69) * ( 0.88 ) + ( energy >=13 && energy < 18 ) * (abs(eta) > 1.95) * ( 0.6 ) + ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.8 ) + ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.68 ) + ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.84 ) + ( energy >=13 && energy < 18 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.88 ) + ( energy >=13 && energy < 18 ) * (abs(eta) <= 0.69) * ( 0.9 ) + ( energy >=18 && energy < 23 ) * (abs(eta) > 1.95) * (0.64 ) + ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.82 ) + ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.7 ) + ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.84 ) + ( energy >=18 && energy < 23 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.9 ) + ( energy >=18 && energy < 23 ) * (abs(eta) <= 0.69) * (0.92 ) + ( energy >= 23 && energy < 28 ) * (abs(eta) > 1.95) * (0.64 ) + ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.86 ) + ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.74 ) + ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.87 ) + ( energy >= 23 && energy < 28 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.91 ) + ( energy >= 23 && energy < 28 ) * (abs(eta) <= 0.69) * (0.94 ) + ( energy >=28 && energy < 35 ) * (abs(eta) > 1.95) * (0.67 ) + ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.88 ) + ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.78 ) + ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.9 ) + ( energy >=28 && energy < 35 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.94 ) + ( energy >=28 && energy < 35 ) * (abs(eta) <= 0.69) * (0.94 ) + ( energy >=35 && energy < 45 ) * (abs(eta) > 1.95) * (0.68 ) + ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.9 ) + ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.86 ) + ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.92 ) + ( energy >=35 && energy < 45 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.94 ) + ( energy >=35 && energy < 45 ) * (abs(eta) <= 0.69) * (0.96 ) + ( energy >=45 && energy < 80 ) * (abs(eta) > 1.95) * ( 0.7 ) + ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.92 ) + ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.8 ) + ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.94 ) + ( energy >=45 && energy < 80 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.96 ) + ( energy >=45 && energy < 80 ) * (abs(eta) <= 0.69) * ( 0.97 ) + ( energy >=80 && energy < 200 ) * (abs(eta) > 1.95) * (0.68 ) + ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.96 ) + ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.84 ) + ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.94 ) + ( energy >=80 && energy < 200 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.98 ) + ( energy >=80 && energy < 200 ) * (abs(eta) <= 0.69) * (0.98 ) + ( energy >=200 && energy < 400 ) * (abs(eta) > 1.95) * ( 0.68 ) + ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.97 ) + ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.86 ) + ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.96 ) + ( energy >=200 && energy < 400 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.98 ) + ( energy >=200 && energy < 400 ) * (abs(eta) <= 0.69) * ( 0.98 ) + ( energy >=400 ) * (abs(eta) > 1.95) * (0.68 ) + ( energy >=400 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.96 ) + ( energy >=400 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.82 ) + ( energy >=400 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.96 ) + ( energy >=400 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.98 ) + ( energy >=400 ) * (abs(eta) <= 0.69) * (0.98 ) } } #################### # Electron isolation #################### module Isolation ElectronIsolation { set CandidateInputArray ElectronEfficiency/electrons set IsolationInputArray EFlowMerger/eflow set OutputArray electrons set DeltaRMax 0.1 set PTMin 0.5 set PTRatioMax 0.2 } ################# # Muon efficiency ################# module Efficiency MuonEfficiency { set InputArray MuonMomentumSmearing/muons set OutputArray muons # set EfficiencyFormula {efficiency as a function of eta and pt} # efficiency formula for muons # current full simulation of CLICdet yields: set EfficiencyFormula { (energy < 2.0 ) * (0.00) + (energy>=2.0) * (0.999) } } ################ # Muon isolation ################ module Isolation MuonIsolation { set CandidateInputArray MuonEfficiency/muons set IsolationInputArray EFlowMerger/eflow set OutputArray muons set DeltaRMax 0.1 set PTMin 0.5 set PTRatioMax 0.2 } ################### # Missing ET merger ################### module Merger MissingET { # add InputArray InputArray add InputArray EFlowMerger/eflow set MomentumOutputArray momentum } ################## # Scalar HT merger ################## module Merger ScalarHT { # add InputArray InputArray add InputArray EFlowMerger/eflow set EnergyOutputArray energy } ###################### # EFlowFilter (UniqueObjectFinder) ###################### module UniqueObjectFinder EFlowFilter { add InputArray PhotonIsolation/photons photons add InputArray ElectronIsolation/electrons electrons add InputArray MuonIsolation/muons muons add InputArray EFlowMerger/eflow eflow } ################# # 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} } ##################### # 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, 7 anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness, 9 Valencia set JetAlgorithm 9 set ParameterR 0.5 set JetPTMin 20.0 } ######################### # Gen Missing ET merger ######################## module Merger GenMissingET { # add InputArray InputArray add InputArray NeutrinoFilter/filteredParticles set MomentumOutputArray momentum } ############ # Jet finder ############ module FastJetFinder FastJetFinderKt { # set InputArray Calorimeter/towers set InputArray EFlowMerger/eflow set OutputArray KTjets # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt, 7 anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness, 9 Valencia set JetAlgorithm 4 set ParameterR 0.5 set JetPTMin 20.0 } ################ # Jet finder VLC ################ source CLIC/CLICdet_JetReco.tcl ######################################### # Jet Momentum Smearing to mimick overlay ######################################### source CLIC/CLICdet_JetSmearing.tcl ######################## # Jet Flavor Association ######################## source CLIC/CLICdet_JetFlavorAssociation.tcl source CLIC/CLICdet_JetFlavorAssociation_JER.tcl ########### # b-tagging ########### # based on CLICdp-Note-2014-002 source CLIC/CLICdet_BTagging.tcl source CLIC/CLICdet_BTagging_JER.tcl ############# # tau-tagging ############# # based on LCD-2010-009 source CLIC/CLICdet_TauTagging.tcl source CLIC/CLICdet_TauTagging_JER.tcl ################## # ROOT tree writer ################## module TreeWriter TreeWriter { # add Branch InputArray BranchName BranchClass add Branch Delphes/allParticles Particle GenParticle add Branch GenJetFinder/jets GenJet Jet add Branch FastJetFinderKt/KTjets KTjet Jet add Branch FastJetFinderVLC_R05_N2/VLCjetsR05N2 VLCjetR05N2 Jet add Branch FastJetFinderVLC_R05_N3/VLCjetsR05N3 VLCjetR05N3 Jet add Branch FastJetFinderVLC_R05_N4/VLCjetsR05N4 VLCjetR05N4 Jet add Branch FastJetFinderVLC_R05_N5/VLCjetsR05N5 VLCjetR05N5 Jet add Branch FastJetFinderVLC_R05_N6/VLCjetsR05N6 VLCjetR05N6 Jet add Branch FastJetFinderVLC_R07_N2/VLCjetsR07N2 VLCjetR07N2 Jet add Branch FastJetFinderVLC_R07_N3/VLCjetsR07N3 VLCjetR07N3 Jet add Branch FastJetFinderVLC_R07_N4/VLCjetsR07N4 VLCjetR07N4 Jet add Branch FastJetFinderVLC_R07_N5/VLCjetsR07N5 VLCjetR07N5 Jet add Branch FastJetFinderVLC_R07_N6/VLCjetsR07N6 VLCjetR07N6 Jet add Branch FastJetFinderVLC_R10_N2/VLCjetsR10N2 VLCjetR10N2 Jet add Branch FastJetFinderVLC_R10_N3/VLCjetsR10N3 VLCjetR10N3 Jet add Branch FastJetFinderVLC_R10_N4/VLCjetsR10N4 VLCjetR10N4 Jet add Branch FastJetFinderVLC_R10_N5/VLCjetsR10N5 VLCjetR10N5 Jet add Branch FastJetFinderVLC_R10_N6/VLCjetsR10N6 VLCjetR10N6 Jet add Branch FastJetFinderVLC_R12_N2/VLCjetsR12N2 VLCjetR12N2 Jet add Branch FastJetFinderVLC_R12_N3/VLCjetsR12N3 VLCjetR12N3 Jet add Branch FastJetFinderVLC_R12_N4/VLCjetsR12N4 VLCjetR12N4 Jet add Branch FastJetFinderVLC_R12_N5/VLCjetsR12N5 VLCjetR12N5 Jet add Branch FastJetFinderVLC_R12_N6/VLCjetsR12N6 VLCjetR12N6 Jet add Branch FastJetFinderVLC_R15_N2/VLCjetsR15N2 VLCjetR15N2 Jet add Branch FastJetFinderVLC_R15_N3/VLCjetsR15N3 VLCjetR15N3 Jet add Branch FastJetFinderVLC_R15_N4/VLCjetsR15N4 VLCjetR15N4 Jet add Branch FastJetFinderVLC_R15_N5/VLCjetsR15N5 VLCjetR15N5 Jet add Branch FastJetFinderVLC_R15_N6/VLCjetsR15N6 VLCjetR15N6 Jet add Branch FastJetFinderVLC_R05_inclusive/VLCjetsR05_inclusive VLCjetR05_inclusive Jet add Branch FastJetFinderVLC_R07_inclusive/VLCjetsR07_inclusive VLCjetR07_inclusive Jet add Branch FastJetFinderVLC_R10_inclusive/VLCjetsR10_inclusive VLCjetR10_inclusive Jet add Branch FastJetFinderVLC_R12_inclusive/VLCjetsR12_inclusive VLCjetR12_inclusive Jet add Branch FastJetFinderVLC_R15_inclusive/VLCjetsR15_inclusive VLCjetR15_inclusive Jet ###with jet energy smearing add Branch JetMomentumSmearing_VLCR05N2/JER_VLCjetsR05N2 JER_VLCjetR05N2 Jet add Branch JetMomentumSmearing_VLCR05N3/JER_VLCjetsR05N3 JER_VLCjetR05N3 Jet add Branch JetMomentumSmearing_VLCR05N4/JER_VLCjetsR05N4 JER_VLCjetR05N4 Jet add Branch JetMomentumSmearing_VLCR05N5/JER_VLCjetsR05N5 JER_VLCjetR05N5 Jet add Branch JetMomentumSmearing_VLCR05N6/JER_VLCjetsR05N6 JER_VLCjetR05N6 Jet add Branch JetMomentumSmearing_VLCR07N2/JER_VLCjetsR07N2 JER_VLCjetR07N2 Jet add Branch JetMomentumSmearing_VLCR07N3/JER_VLCjetsR07N3 JER_VLCjetR07N3 Jet add Branch JetMomentumSmearing_VLCR07N4/JER_VLCjetsR07N4 JER_VLCjetR07N4 Jet add Branch JetMomentumSmearing_VLCR07N5/JER_VLCjetsR07N5 JER_VLCjetR07N5 Jet add Branch JetMomentumSmearing_VLCR07N6/JER_VLCjetsR07N6 JER_VLCjetR07N6 Jet add Branch JetMomentumSmearing_VLCR10N2/JER_VLCjetsR10N2 JER_VLCjetR10N2 Jet add Branch JetMomentumSmearing_VLCR10N3/JER_VLCjetsR10N3 JER_VLCjetR10N3 Jet add Branch JetMomentumSmearing_VLCR10N4/JER_VLCjetsR10N4 JER_VLCjetR10N4 Jet add Branch JetMomentumSmearing_VLCR10N5/JER_VLCjetsR10N5 JER_VLCjetR10N5 Jet add Branch JetMomentumSmearing_VLCR10N6/JER_VLCjetsR10N6 JER_VLCjetR10N6 Jet add Branch JetMomentumSmearing_VLCR12N2/JER_VLCjetsR12N2 JER_VLCjetR12N2 Jet add Branch JetMomentumSmearing_VLCR12N3/JER_VLCjetsR12N3 JER_VLCjetR12N3 Jet add Branch JetMomentumSmearing_VLCR12N4/JER_VLCjetsR12N4 JER_VLCjetR12N4 Jet add Branch JetMomentumSmearing_VLCR12N5/JER_VLCjetsR12N5 JER_VLCjetR12N5 Jet add Branch JetMomentumSmearing_VLCR12N6/JER_VLCjetsR12N6 JER_VLCjetR12N6 Jet add Branch JetMomentumSmearing_VLCR15N2/JER_VLCjetsR15N2 JER_VLCjetR15N2 Jet add Branch JetMomentumSmearing_VLCR15N3/JER_VLCjetsR15N3 JER_VLCjetR15N3 Jet add Branch JetMomentumSmearing_VLCR15N4/JER_VLCjetsR15N4 JER_VLCjetR15N4 Jet add Branch JetMomentumSmearing_VLCR15N5/JER_VLCjetsR15N5 JER_VLCjetR15N5 Jet add Branch JetMomentumSmearing_VLCR15N6/JER_VLCjetsR15N6 JER_VLCjetR15N6 Jet add Branch JetMomentumSmearing_VLCR05_inclusive/JER_VLCjetsR05_inclusive JER_VLCjetR05_inclusive Jet add Branch JetMomentumSmearing_VLCR07_inclusive/JER_VLCjetsR07_inclusive JER_VLCjetR07_inclusive Jet add Branch JetMomentumSmearing_VLCR10_inclusive/JER_VLCjetsR10_inclusive JER_VLCjetR10_inclusive Jet add Branch JetMomentumSmearing_VLCR12_inclusive/JER_VLCjetsR12_inclusive JER_VLCjetR12_inclusive Jet add Branch JetMomentumSmearing_VLCR15_inclusive/JER_VLCjetsR15_inclusive JER_VLCjetR15_inclusive Jet #### add Branch GenMissingET/momentum GenMissingET MissingET add Branch TrackMerger/tracks Track Track add Branch Calorimeter/towers Tower Tower add Branch HCal/eflowTracks EFlowTrack Track add Branch ECal/eflowPhotons EFlowPhoton Tower add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower add Branch EFlowFilter/photons Photon Photon add Branch EFlowFilter/electrons Electron Electron add Branch EFlowFilter/muons Muon Muon add Branch MissingET/momentum MissingET MissingET add Branch ScalarHT/energy ScalarHT ScalarHT add Branch BCalTowers/bcalTowers BCalTower Tower add Branch LumiTowers/lumicalTowers LumiTower Tower add Branch LumiEfficiency/lumicalPhotons LumiPhoton Photon add Branch BCalEfficiency/bcalPhotons BCalPhoton Photon }