Changeset 26a21bd in git for cards/delphes_card_IDEA.tcl

Timestamp:

Oct 15, 2021, 5:23:17 PM (3 years ago)

Author:

michele <michele.selvaggi@…>

Branches:

master

Children:

78e1fb3

Parents:

a3784f9

Message:

added missing loopers to eflowtracks, fixed missing bin in dr calo

File:

• cards/delphes_card_IDEA.tcl (modified) (17 diffs)

cards/delphes_card_IDEA.tcl

@@ +1 @@
+set RandomSeed 123
+
 ####################################################################                                l
 # FCC-ee IDEA detector model
@@ -40 +42 @@
 
   TrackMerger
+  ForwardLooperTracks
   Calorimeter
 
@@ -45 +48 @@
   TimeOfFlightNeutralHadron
 
+  EFlowTrackMerger
   EFlowMerger
 
@@ -121 +125 @@
     set UseMomentumVector true
 
-    # We use only one efficiency, we set only 0 effincency out of eta bounds:
     set EfficiencyFormula {
         (abs(eta) > 2.56)                                  * (0.000) +
-        (pt < 0.1) * (abs(eta) <= 2.56)                    * (0.000)
+        (pt < 0.1) * (abs(eta) <= 2.56)                    * (0.000) +
         (pt >= 0.1) * (abs(eta) <= 2.56)                   * (1.000)
     }
@@ -140 +143 @@
     set UseMomentumVector true
 
-    # Current full simulation with CLICdet provides for electrons:
     set EfficiencyFormula {
         (abs(eta) > 2.56)                                  * (0.000) +
-        (pt < 0.1) * (abs(eta) <= 2.56)                    * (0.000)
+        (pt < 0.1) * (abs(eta) <= 2.56)                    * (0.000) +
         (pt >= 0.1) * (abs(eta) <= 2.56)                   * (1.000)
     }
@@ -158 +160 @@
     set UseMomentumVector true
 
-    # Current full simulation with CLICdet provides for muons:
     set EfficiencyFormula {
         (abs(eta) > 2.56)                                  * (0.000) +
-        (pt < 0.1) * (abs(eta) <= 2.56)                    * (0.000)
+        (pt < 0.1) * (abs(eta) <= 2.56)                    * (0.000) +
         (pt >= 0.1) * (abs(eta) <= 2.56)                   * (1.000)
     }
@@ -177 +178 @@
   set OutputArray tracks
 }
+
 
 
@@ -437 +439 @@
 
 
+######################
+# Looper Selection
+######################
+
+module Efficiency ForwardLooperTracks  {
+  set InputArray TrackMerger/tracks
+  set OutputArray tracks
+  set UseMomentumVector False
+
+  ## select looping tracks that end up in position |eta| > 3.142 (lost by calo)
+  set EfficiencyFormula {
+    (abs(eta) > 3.0 )                                 * (1.000) +
+    (abs(eta) <= 3.0 )                                * (0.000)
+  }
+
+}
+
+
 #############
 # Calorimeter
@@ -488 +508 @@
     #deta=0.02 units for 0.88 < |eta| <= 3.0
     #first, from -3.0 to -0.88
-    for {set i 1} {$i <=106} {incr i} {
+    for {set i 0} {$i <=106} {incr i} {
         set eta [expr {-3.00 + $i * 0.02}]
         add EtaPhiBins $eta $PhiBins
@@ -565 +585 @@
 }
 
+
+############################
+# Energy flow track merger
+############################
+
+module Merger EFlowTrackMerger {
+# add InputArray InputArray
+  add InputArray Calorimeter/eflowTracks
+  add InputArray ForwardLooperTracks/tracks
+  set OutputArray eflowTracks
+}
+
+
+
 ####################
 # Energy flow merger
@@ -571 +605 @@
 module Merger EFlowMerger {
 # add InputArray InputArray
-  add InputArray Calorimeter/eflowTracks
+  add InputArray EFlowTrackMerger/eflowTracks
   add InputArray Calorimeter/eflowPhotons
   add InputArray TimeOfFlightNeutralHadron/eflowNeutralHadrons
@@ -585 +619 @@
   set InputArray Calorimeter/eflowPhotons
   set OutputArray photons
-  set UseMomentumVector true
 
   # set EfficiencyFormula {efficiency formula as a function of eta and pt}
@@ -619 +652 @@
 
 module PdgCodeFilter ElectronFilter {
-  set InputArray Calorimeter/eflowTracks
+  set InputArray EFlowTrackMerger/eflowTracks
   set OutputArray electrons
   set Invert true
@@ -631 +664 @@
 
 module PdgCodeFilter MuonFilter {
-  set InputArray Calorimeter/eflowTracks
+  set InputArray EFlowTrackMerger/eflowTracks
   set OutputArray muons
   set Invert true
@@ -646 +679 @@
   set InputArray ElectronFilter/electrons
   set OutputArray electrons
-  set UseMomentumVector true
 
   # set EfficiencyFormula {efficiency formula as a function of eta and pt}
@@ -683 +715 @@
   set InputArray MuonFilter/muons
   set OutputArray muons
-  set UseMomentumVector true
 
   # set EfficiencyFormula {efficiency as a function of eta and pt}
@@ -885 +916 @@
     add Branch TruthVertexFinder/vertices GenVertex Vertex
 
-    add Branch TrackMerger/tracks Track Track
-    add Branch Calorimeter/towers Tower Tower
-
-    add Branch Calorimeter/eflowTracks EFlowTrack Track
+    add Branch EFlowTrackMerger/eflowTracks EFlowTrack Track
     add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
     add Branch TimeOfFlightNeutralHadron/eflowNeutralHadrons EFlowNeutralHadron Tower