Changeset 1370 in svn for trunk/examples/delphes_card_TLEP.tcl

Timestamp:

Apr 17, 2014, 10:23:04 AM (10 years ago)

Author:

Michele Selvaggi

Message:

committed delphes_card_FCC_basic.tcl

File:

• trunk/examples/delphes_card_TLEP.tcl (modified) (7 diffs)

trunk/examples/delphes_card_TLEP.tcl

@@ -15 +15 @@
 
   TrackMerger
-  Calorimeter
+  ImpactParameterSmearing
+  
+  Ecal
+  Hcal
+  
+  TowerMerger
   EFlowMerger
 
@@ -24 +29 @@
 
   JetEnergyScale
-
-  BTagging
+  
+  TrackCountingBTagging
   TauTagging
 
@@ -187 +192 @@
 }
 
+################################
+# Track impact parameter smearing
+################################
+
+module ImpactParameterSmearing ImpactParameterSmearing {
+  set InputArray TrackMerger/tracks
+  set OutputArray tracks
+
+
+# absolute impact parameter smearing formula (in mm) as a function of pt and eta
+set ResolutionFormula {(pt > 0.1  && pt <= 5.0)   * (0.010) + \
+                       (pt > 5.0)                 * (0.005)}
+
+}
+
 #############
-# Calorimeter
+#   ECAL
 #############
 
-module Calorimeter Calorimeter {
+module SimpleCalorimeter Ecal {
   set ParticleInputArray ParticlePropagator/stableParticles
-  set TrackInputArray TrackMerger/tracks
-
-  set TowerOutputArray towers
-  set PhotonOutputArray photons
-
-  set EFlowTrackOutputArray eflowTracks
-  set EFlowPhotonOutputArray eflowPhotons
-  set EFlowNeutralHadronOutputArray eflowNeutralHadrons
+  set TrackInputArray ImpactParameterSmearing/tracks
+
+  set TowerOutputArray ecalTowers
+  set EFlowTowerOutputArray eflowPhotons
   
   set pi [expr {acos(-1)}]
@@ -234 +250 @@
     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) <= 4.0) * sqrt(energy^2*0.005^2 + energy*0.02^2) + \
+                            (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.05^2  + energy*1.00^2)}
+
+ 
+}
+
+#############
+#   HCAL
+#############
+
+module SimpleCalorimeter Hcal {
+  set ParticleInputArray ParticlePropagator/stableParticles
+  set TrackInputArray ImpactParameterSmearing/tracks
+
+  set TowerOutputArray hcalTowers
+  set EFlowTowerOutputArray eflowNeutralHadrons
+  
+  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
+
+  # 5 degrees towers
+  set PhiBins {}
+  for {set i -36} {$i <= 36} {incr i} {
+    add PhiBins [expr {$i * $pi/36.0}]
+  }
+  foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
+    add EtaPhiBins $eta $PhiBins
+  }
+
+  # 10 degrees towers
+  set PhiBins {}
+  for {set i -18} {$i <= 18} {incr i} {
+    add PhiBins [expr {$i * $pi/18.0}]
+  }
+  foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
+    add EtaPhiBins $eta $PhiBins
+  }
+
+  # 20 degrees towers
+  set PhiBins {}
+  for {set i -9} {$i <= 9} {incr i} {
+    add PhiBins [expr {$i * $pi/9.0}]
+  }
+  foreach eta {6 -5.6 -5.3 -5 -4.7 -4.525 4.7 5 5.3 5.6 6} {
+    add EtaPhiBins $eta $PhiBins
+  }
+  
   # default energy fractions {abs(PDG code)} {Fecal Fhcal}
-  add EnergyFraction {0} {0.0 1.0}
+  add EnergyFraction {0} {1.0}
   # energy fractions for e, gamma and pi0
-  add EnergyFraction {11} {1.0 0.0}
-  add EnergyFraction {22} {1.0 0.0}
-  add EnergyFraction {111} {1.0 0.0}
+  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 0.0}
-  add EnergyFraction {13} {0.0 0.0}
-  add EnergyFraction {14} {0.0 0.0}
-  add EnergyFraction {16} {0.0 0.0}
-  add EnergyFraction {1000022} {0.0 0.0}
-  add EnergyFraction {1000023} {0.0 0.0}
-  add EnergyFraction {1000025} {0.0 0.0}
-  add EnergyFraction {1000035} {0.0 0.0}
-  add EnergyFraction {1000045} {0.0 0.0}
+  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 0.7}
-  add EnergyFraction {3122} {0.3 0.7}
-
-  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
-  set ECalResolutionFormula {                 (abs(eta) <= 4.0) * sqrt(energy^2*0.005^2 + energy*0.02^2) + \
-                            (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.05^2  + energy*1.00^2)}
-
-  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
-  set HCalResolutionFormula {                 (abs(eta) <= 4.0) * sqrt(energy^2*0.03^2 + energy*0.50^2) + \
+  add EnergyFraction {310} {0.7}
+  add EnergyFraction {3122} {0.7}
+
+   # set HCalResolutionFormula {resolution formula as a function of eta and energy}
+  set ResolutionFormula {                     (abs(eta) <= 4.0) * sqrt(energy^2*0.03^2 + energy*0.50^2) + \
                             (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.05^2 + energy*1.00^2)}
+}
+
+
+####################
+# Tower Merger (in case not using e-flow algorithm)
+####################
+
+module Merger TowerMerger {
+# add InputArray InputArray
+  add InputArray Ecal/ecalTowers
+  add InputArray Hcal/hcalTowers
+  set OutputArray towers
 }
 
@@ -270 +365 @@
 module Merger EFlowMerger {
 # add InputArray InputArray
-  add InputArray Calorimeter/eflowTracks
-  add InputArray Calorimeter/eflowPhotons
-  add InputArray Calorimeter/eflowNeutralHadrons
+  add InputArray ImpactParameterSmearing/tracks
+  add InputArray Ecal/eflowPhotons
+  add InputArray Hcal/eflowNeutralHadrons
   set OutputArray eflow
 }
+
 
 ###################
@@ -342 +438 @@
 }
 
-###########
-# b-tagging
-###########
-
-module BTagging BTagging {
-  set PartonInputArray Delphes/partons
+##########################
+# Track Counting b-tagging
+##########################
+
+module TrackCountingBTagging TrackCountingBTagging {
+  set TrackInputArray ImpactParameterSmearing/tracks
   set JetInputArray JetEnergyScale/jets
 
   set BitNumber 0
 
-  set DeltaR 0.5
-
-  set PartonPTMin 1.0
-
-  set PartonEtaMax 4.0
-
-  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
-  # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
-  # gluon's PDG code has the lowest priority
-
-  # https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsBTV
-  # default efficiency formula (misidentification rate)
-  add EfficiencyFormula {0} {0.001}
-
-  # efficiency formula for c-jets (misidentification rate)
-  add EfficiencyFormula {4} {                                      (pt <= 15.0) * (0.000) + \
-                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
-                              (abs(eta) > 1.2 && abs(eta) <= 4.0) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
-                              (abs(eta) > 4.0)                                  * (0.000)}
-
-  # efficiency formula for b-jets
-  add EfficiencyFormula {5} {                                      (pt <= 15.0) * (0.000) + \
-                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
-                              (abs(eta) > 1.2 && abs(eta) <= 4.0) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
-                              (abs(eta) > 4.0)                                  * (0.000)}
-}
+  # maximum distance between jet and track
+  set DeltaR 0.3
+  
+  # minimum pt of tracks 
+  set TrackPTMin 1.0
+  
+  # minimum transverse impact parameter (in mm)
+  set TrackIPMax 2.0
+  
+  # minimum ip significance for the track to be counted
+  set SigMin 6.5
+  
+  # minimum number of tracks (high efficiency n=2, high purity n=3)
+  set Ntracks 3
+}
+
+
+##########################
+# tau-tagging
+##########################
+
 
 module TauTagging TauTagging {
@@ -406 +497 @@
   add Branch Delphes/allParticles Particle GenParticle
   add Branch GenJetFinder/jets GenJet Jet
-  add Branch Calorimeter/eflowNeutralHadrons NeutralHadron Tower
+  
+  add Branch Hcal/eflowNeutralHadrons NeutralHadron Tower
+  add Branch Ecal/eflowPhotons Photon Photon
   add Branch ChargedHadronMomentumSmearing/chargedHadrons ChargedHadron Track
-  add Branch Calorimeter/eflowPhotons Photon Photon
   add Branch ElectronEnergySmearing/electrons Electron Electron
   add Branch MuonMomentumSmearing/muons Muon Muon