Changeset b24f05a in git

Timestamp:

Jun 6, 2018, 7:34:16 PM (6 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

8a58fff

Parents:

ec37bc3

Message:

push PhaseII final card

Files:

• cards/CMS_PhaseII/CMS_PhaseII_200PU_v03.tcl (modified) (12 diffs)
• modules/DenseTrackFilter.cc (modified) (3 diffs)

cards/CMS_PhaseII/CMS_PhaseII_200PU_v03.tcl

@@ -1 @@
-set RandomSeed 123
-set MaxEvents 10
-
 ########################################
 #
@@ -38 +35 @@
   HCal
 
-  MuonFilter
-
   PhotonEnergySmearing
   ElectronFilter
@@ -167 +162 @@
 
 
-#################################
-# Dense Track propagation calo
-#################################
+#####################################
+# Track propagation to calorimeters
+#####################################
 
 module ParticlePropagator ParticlePropagator {
@@ -203 +198 @@
 
 
-#################################
-# Dense Track propagation
-#################################
+####################################
+# Track propagation to pseudo-pixel
+####################################
 
 module ParticlePropagator DenseProp {
 
   set InputArray TrackMergerProp/tracks
-
-  set OutputArray stableParticles
-  set NeutralOutputArray neutralParticles
-  set ChargedHadronOutputArray chargedHadrons
-  set ElectronOutputArray electrons
-  set MuonOutputArray muons
 
   # radius of the first pixel layer
@@ -246 +235 @@
 
 module DenseTrackFilter DenseTrackFilter {
+
   set TrackInputArray DenseMergeTracks/tracks
+
   set TrackOutputArray tracks
-
   set ChargedHadronOutputArray chargedHadrons
   set ElectronOutputArray electrons
@@ -446 +436 @@
 
 
+#################################
+# Energy resolution for electrons
+#################################
+
+module EnergySmearing ElectronEnergySmearing {
+  set InputArray ElectronTrackingEfficiency/electrons
+  set OutputArray electrons
+
+  # set ResolutionFormula {resolution formula as a function of eta and energy}
+
+  # resolution formula for electrons
+
+  # taking something flat in energy for now, ECAL will take over at high energy anyway.
+  # inferred from hep-ex/1306.2016 and 1502.02701
+  set ResolutionFormula {
+
+                        (abs(eta) <= 1.5)  * (energy*0.028) +
+    (abs(eta) > 1.5  && abs(eta) <= 1.75)  * (energy*0.037) +
+    (abs(eta) > 1.75  && abs(eta) <= 2.15) * (energy*0.038) +
+    (abs(eta) > 2.15  && abs(eta) <= 3.00) * (energy*0.044) +
+    (abs(eta) > 3.00  && abs(eta) <= 4.00) * (energy*0.10)}
+
+}
+
+###############################
+# Momentum resolution for muons
+###############################
+
+module MomentumSmearing MuonMomentumSmearing {
+  set InputArray MuonTrackingEfficiency/muons
+  set OutputArray muons
+  # resolution formula for muons
+
+  # up to |eta| < 2.8 take measurement from tracking + muon chambers
+  # for |eta| > 2.8 and pT < 5.0 take measurement from tracking alone taken from
+  # http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
+  source muonMomentumResolution.tcl
+}
+
+
+
 ##############
 # Track merger
@@ -457 +488 @@
   set OutputArray tracks
 }
-
-#################################
-# Energy resolution for electrons
-#################################
-
-module EnergySmearing ElectronEnergySmearing {
-  set InputArray ElectronTrackingEfficiency/electrons
-  set OutputArray electrons
-
-  # set ResolutionFormula {resolution formula as a function of eta and energy}
-
-  # resolution formula for electrons
-
-  # taking something flat in energy for now, ECAL will take over at high energy anyway.
-  # inferred from hep-ex/1306.2016 and 1502.02701
-  set ResolutionFormula {
-
-                        (abs(eta) <= 1.5)  * (energy*0.028) +
-    (abs(eta) > 1.5  && abs(eta) <= 1.75)  * (energy*0.037) +
-    (abs(eta) > 1.75  && abs(eta) <= 2.15) * (energy*0.038) +
-    (abs(eta) > 2.15  && abs(eta) <= 3.00) * (energy*0.044) +
-    (abs(eta) > 3.00  && abs(eta) <= 4.00) * (energy*0.10)}
-
-}
-
-###############################
-# Momentum resolution for muons
-###############################
-
-module MomentumSmearing MuonMomentumSmearing {
-  set InputArray MuonTrackingEfficiency/muons
-  set OutputArray muons
-  # resolution formula for muons
-
-  # up to |eta| < 2.8 take measurement from tracking + muon chambers
-  # for |eta| > 2.8 and pT < 5.0 take measurement from tracking alone taken from
-  # http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
-  source muonMomentumResolution.tcl
-}
-
-
 
 
@@ -721 +711 @@
 }
 
-#################
-# Muon filter
-#################
-
-module PdgCodeFilter MuonFilter {
-  set InputArray HCal/eflowTracks
-  set OutputArray muons
-  set Invert true
-  add PdgCode {13}
-  add PdgCode {-13}
-}
 
 
@@ -742 +721 @@
   add InputArray HCal/eflowTracks eflowTracks
   add InputArray ElectronFilter/electrons electrons
-  add InputArray MuonFilter/muons muons
+  add InputArray MuonMomentumSmearing/muons muons
 
   set VertexInputArray PileUpMerger/vertices
@@ -1200 +1179 @@
 
 module Cloner MuonCloner {
-  set InputArray MuonFilter/muons
+  set InputArray MuonMomentumSmearing/muons
   set OutputArray muons
 }
@@ -2207 +2186 @@
 
 module Isolation MuonIsolation {
-  set CandidateInputArray MuonFilter/muons
+  set CandidateInputArray MuonMomentumSmearing/muons
 
   # isolation collection
@@ -4122 +4101 @@
   add Branch GenMissingET/momentum GenMissingET MissingET
 
-  add Branch TrackMerger/tracks TrackMerger Track
-  add Branch HCal/eflowTracks EFlowTrack Track
 #  add Branch HCal/eflowTracks EFlowTrack Track
 #  add Branch ECal/eflowPhotons EFlowPhoton Tower

modules/DenseTrackFilter.cc

@@ -243 +243 @@
   Int_t numberOfCandidates;
 
-  
-  cout<<fBestTrack<<endl;
-  
   // saving track with highest pT that hit the tower
   if(fTowerTrackHits > 0) 
@@ -251 +248 @@
     
     numberOfCandidates = fBestTrack->GetCandidates()->GetEntriesFast();
-    if (numberOfCandidates < 2)
-       cout<<fBestTrack<<endl;
-
     if (numberOfCandidates > 1)
     {  
@@ -278 +272 @@
           fChargedHadronOutputArray->Add(candidate);
       }
- 
-    }
-  }
-
-}
+
+    }
+  }
+
+}