Changes in / [efc1546:c0756c4] in git

Location:

cards/CMS_PhaseII

Files:

• CMS_PhaseII_0PU.tcl (modified) (8 diffs)
• CMS_PhaseII_200PU.tcl (modified) (10 diffs)

cards/CMS_PhaseII/CMS_PhaseII_0PU.tcl

@@ -29 +29 @@
   HCal
 
+  PhotonEnergySmearing
+
   TowerMerger
   EFlowMerger
@@ -50 +52 @@
   GenJetFinder
   FastJetFinder
+  FatJetFinder
 
   ScalarHT
@@ -435 +438 @@
 }
 
+#################################
+# Energy resolution for electrons
+#################################
+
+module EnergySmearing PhotonEnergySmearing {
+  set InputArray ECal/eflowPhotons
+  set OutputArray eflowPhotons
+
+  # adding 1% extra photon smearing
+  set ResolutionFormula {energy*0.01}
+
+}
+
+
 ####################
 # Energy flow merger
@@ -442 +459 @@
 # add InputArray InputArray
   add InputArray HCal/eflowTracks
-  add InputArray ECal/eflowPhotons
+  add InputArray PhotonEnergySmearing/eflowPhotons
   add InputArray HCal/eflowNeutralHadrons
   set OutputArray eflow
@@ -537 +554 @@
 }
 
+
+
+
+############
+# Jet finder
+############
+
+module FastJetFinder FatJetFinder {
+#  set InputArray TowerMerger/towers
+  set InputArray EFlowMerger/eflow
+
+  set OutputArray jets
+
+  set JetAlgorithm 5
+  set ParameterR 0.8
+
+  set ComputeNsubjettiness 1
+  set Beta 1.0
+  set AxisMode 4
+ 
+  set ComputeTrimming 1
+  set RTrim 0.2
+  set PtFracTrim 0.05
+
+  set ComputePruning 1 
+  set ZcutPrun 0.1
+  set RcutPrun 0.5
+  set RPrun 0.8
+
+  set ComputeSoftDrop 1
+  set BetaSoftDrop 0.0
+  set SymmetryCutSoftDrop 0.1
+  set R0SoftDrop 0.8
+
+  set JetPTMin 200.0
+}
+
+
+
 ##################
 # Jet Energy Scale
@@ -556 +612 @@
   
   # particle for which calculate the isolation
-  set CandidateInputArray ECal/eflowPhotons
+  set CandidateInputArray PhotonEnergySmearing/eflowPhotons
   
   # isolation collection
@@ -801 +857 @@
 
   add Branch HCal/eflowTracks EFlowTrack Track
-  add Branch ECal/eflowPhotons EFlowPhoton Tower
+  add Branch PhotonEnergySmearing/eflowPhotons EFlowPhoton Tower
   add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
   
@@ -810 +866 @@
   
   add Branch JetEnergyScale/jets Jet Jet
+  add Branch FatJetFinder/jets FatJet Jet
   
   add Branch MissingET/momentum MissingET MissingET

cards/CMS_PhaseII/CMS_PhaseII_200PU.tcl

@@ -30 +30 @@
   ECal
   HCal
-
+  
+  PhotonEnergySmearing
   ElectronFilter
   TrackPileUpSubtractor
@@ -65 +66 @@
   GenJetFinder
   FastJetFinder
+  FatJetFinder
 
   ScalarHT
@@ -96 +98 @@
 
   # pre-generated minbias input file
-  set PileUpFile ../eos/cms/store/group/upgrade/delphes/PhaseII/MinBias_100k.pileup
+  #set PileUpFile ../eos/cms/store/group/upgrade/delphes/PhaseII/MinBias_100k.pileup
+  set PileUpFile MinBias.pileup
 
   # average expected pile up
@@ -464 +467 @@
 }
 
+#################################
+# Energy resolution for electrons
+#################################
+
+module EnergySmearing PhotonEnergySmearing {
+  set InputArray ECal/eflowPhotons
+  set OutputArray eflowPhotons
+
+  # adding 1% extra photon smearing
+  set ResolutionFormula {energy*0.01}
+
+}
+
+
+
 #################
 # Electron filter
@@ -512 +530 @@
 module Merger NeutralEFlowMerger {
 # add InputArray InputArray
-  add InputArray ECal/eflowPhotons
+  add InputArray PhotonEnergySmearing/eflowPhotons
   add InputArray HCal/eflowNeutralHadrons
   set OutputArray eflowTowers
@@ -525 +543 @@
 # add InputArray InputArray
   add InputArray HCal/eflowTracks
-  add InputArray ECal/eflowPhotons
+  add InputArray PhotonEnergySmearing/eflowPhotons
   add InputArray HCal/eflowNeutralHadrons
   set OutputArray eflow
@@ -537 +555 @@
 # add InputArray InputArray
   add InputArray TrackMerger/tracks
-  add InputArray ECal/eflowPhotons
+  add InputArray PhotonEnergySmearing/eflowPhotons
   add InputArray HCal/eflowNeutralHadrons
   set OutputArray eflow
@@ -708 +726 @@
 }
 
+
+
+############
+# Jet finder
+############
+
+module FastJetFinder FatJetFinder {
+#  set InputArray TowerMerger/towers
+  set InputArray RunPUPPI/PuppiParticles
+
+  set OutputArray jets
+
+  set JetAlgorithm 5
+  set ParameterR 0.8
+
+  set ComputeNsubjettiness 1
+  set Beta 1.0
+  set AxisMode 4
+
+  set ComputeTrimming 1
+  set RTrim 0.2
+  set PtFracTrim 0.05
+
+  set ComputePruning 1
+  set ZcutPrun 0.1
+  set RcutPrun 0.5
+  set RPrun 0.8
+
+  set ComputeSoftDrop 1
+  set BetaSoftDrop 0.0
+  set SymmetryCutSoftDrop 0.1
+  set R0SoftDrop 0.8
+
+  set JetPTMin 200.0
+}
+
+
 ##################
 # Jet Energy Scale
@@ -726 +781 @@
 
 module PdgCodeFilter PhotonFilter {
-  set InputArray ECal/eflowPhotons
+  set InputArray PhotonEnergySmearing/eflowPhotons
   set OutputArray photons
   set Invert true
@@ -1016 +1071 @@
 
   add Branch JetEnergyScale/jets Jet Jet
+  add Branch FatJetFinder/jets FatJet Jet
 
   add Branch MissingET/momentum MissingET MissingET