Merge pull request cms-sw#200 from mariadalfonso/MT272Xconfig

MT2 72x config

CMGTools/TTHAnalysis/cfg/run_susyMT2_cfg.py

@@ -20,6 +20,7 @@
 lepAna.loose_electron_relIso = 0.15
 lepAna.loose_electron_lostHits = 999 # no cut
 lepAna.inclusive_electron_lostHits = 999 # no cut
+lepAna.mu_isoCorr = "deltaBeta"
 lepAna.ele_isoCorr = "deltaBeta"
 lepAna.ele_tightId = "Cuts_2012"
 lepAna.notCleaningElectrons = True
@@ -30,7 +31,7 @@
 jetAna.jetEta = 5.2
 jetAna.jetEtaCentral = 2.5
 jetAna.jetPt = 10.
-jetAna.recalibrateJets = False
+jetAna.recalibrateJets = True
 jetAna.jetLepDR = 0.4
 jetAna.smearJets = False
 jetAna.jetGammaDR = 0.4
@@ -52,6 +53,9 @@
 # Isolated Track
 isoTrackAna.setOff=False
 
+# recalibrate MET
+metAna.recalibrate = False
+
 ##------------------------------------------ 
 ##  CONTROL VARIABLES
 ##------------------------------------------ 
@@ -78,13 +82,13 @@
 ##  PRODUCER
 ##------------------------------------------
 
-from CMGTools.TTHAnalysis.samples.samples_13TeV_PHYS14 import triggers_HT900, triggers_MET170, triggers_HTMET, triggers_MT2_mumu, triggers_MT2_ee, triggers_MT2_mue, triggers_1mu, triggers_photon155
+from CMGTools.TTHAnalysis.samples.samples_13TeV_PHYS14 import triggers_HT900, triggers_MET170, triggers_HTMET, triggers_MT2_mumu, triggers_MT2_ee, triggers_MT2_mue, triggers_1mu, triggers_photon155,triggers_1mu_isolow
 
 triggerFlagsAna.triggerBits = {
-            'HT650' : triggers_HT900,
-            'MET150' : triggers_MET170,
-            'ht350met100' : triggers_HTMET,
-            'SingleMu' : triggers_1mu,
+            'HT900' : triggers_HT900,
+            'MET170' : triggers_MET170,
+            'ht350met120' : triggers_HTMET,
+            'SingleMu' : triggers_1mu_isolow,
             'DoubleMu' : triggers_MT2_mumu,
             'DoubleEl' : triggers_MT2_ee,
             'MuEG'     : triggers_MT2_mue,
@@ -107,7 +111,11 @@
      collections = susyFullHad_collections,
 )
 
-sequence = cfg.Sequence(susyCoreSequence+[
+
+susyCoreSequence.insert(susyCoreSequence.index(skimAnalyzer),
+                        susyCounter)
+sequence = cfg.Sequence(
+    susyCoreSequence+[
     ttHMT2Control,
     ttHTopoJetAna,
     ttHFatJetAna,
@@ -161,8 +169,13 @@
     ## 50 ns ttbar CSAv2
     #    comp=TTJets
     ## 25 ns ttbar PHYS14
-    comp = TTJets
-    comp.files = comp.files[:1]
+#    comp = TTJets
+#    comp.files = comp.files[:1]
+
+#    comp=TTJets
+#    comp.files = ['/afs/cern.ch/work/d/dalfonso/public/ttjets_miniaodsim_00C90EFC-3074-E411-A845-002590DB9262.root']
+    comp=GJets_HT200to400
+    comp.files = ['/afs/cern.ch/work/d/dalfonso/public/gjets_ht200to400_miniaodsim_060B8ED3-8571-E411-A2CD-002590D0AFEA.root']
     selectedComponents = [comp]
     comp.splitFactor = 10
 elif test==2:
@@ -173,11 +186,20 @@
         comp.files = comp.files[:]
         #comp.files = comp.files[:1]
 
+from PhysicsTools.HeppyCore.framework.services.tfile import TFileService 
+output_service = cfg.Service(
+      TFileService,
+      'outputfile',
+      name="outputfile",
+      fname='mt2.root',
+      option='recreate'
+    )
+
 # the following is declared in case this cfg is used in input to the heppy.py script                                                                                                                   
 from PhysicsTools.HeppyCore.framework.eventsfwlite import Events
 config = cfg.Config( components = selectedComponents,
                      sequence = sequence,
-                     services = [],
+                     services = [output_service],
                      events_class = Events)
 
 #printComps(config.components, True)

CMGTools/TTHAnalysis/python/analyzers/susyCore_modules_cff.py

@@ -9,6 +9,11 @@
 from PhysicsTools.Heppy.analyzers.gen.all import *
 import os
 
+from CMGTools.TTHAnalysis.analyzers.ttHhistoCounterAnalyzer import ttHhistoCounterAnalyzer
+susyCounter = cfg.Analyzer(
+    ttHhistoCounterAnalyzer, name="ttHhistoCounterAnalyzer",
+    )
+
 PDFWeights = []
 #PDFWeights = [ ("CT10",53), ("MSTW2008lo68cl",41), ("NNPDF21_100",101) ]
 
@@ -60,6 +65,7 @@
         "trkPOG_manystripclus53X" : [ "Flag_trkPOG_manystripclus53X" ],
         "trkPOG_toomanystripclus53X" : [ "Flag_trkPOG_toomanystripclus53X" ],
         "trkPOG_logErrorTooManyClusters" : [ "Flag_trkPOG_logErrorTooManyClusters" ],
+        "METFilters" : [ "Flag_METFilters" ],
     }
     )
 
@@ -294,6 +300,7 @@
     METAnalyzer, name="metAnalyzer",
     doTkMet = False,
     doMetNoMu = False,
+    recalibrate = True,
     candidates='packedPFCandidates',
     candidatesTypes='std::vector<pat::PackedCandidate>',
     dzMax = 0.1,

CMGTools/TTHAnalysis/python/analyzers/treeProducerSusyFullHad.py

@@ -115,13 +115,13 @@
 susyFullHad_collections.update({
         # put more here
 ##        "gennus"         : NTupleCollection("genNu",     genParticleWithSourceType, 10, help="Generated neutrinos (nue/numu/nutau) from W/Z decays"),
-        "selectedLeptons" : NTupleCollection("lep", leptonTypeSusy, 50, help="Leptons after the preselection", filter=lambda l : l.pt()>10 ),
+        "selectedLeptons" : NTupleCollection("lep", leptonType, 50, help="Leptons after the preselection", filter=lambda l : l.pt()>10 ),
         "selectedTaus"    : NTupleCollection("tau", tauTypeSusy, 50, help="Taus after the preselection"),
         "cleanJetsAll"       : NTupleCollection("jet",     jetTypeSusy, 100, help="all jets (w/ x-cleaning, w/ ID applied w/o PUID applied pt>10 |eta|<5.2) , sorted by pt", filter=lambda l : l.pt()>10  ),
         "fatJets"         : NTupleCollection("fatJet", fatJetType, 15, help="Cental jets after full selection and cleaning, sorted by pt"),
         "selectedPhotons"    : NTupleCollection("gamma", photonTypeSusy, 50, help="photons with pt>20 and loose cut based ID"),
         "selectedIsoTrack"    : NTupleCollection("isoTrack", isoTrackType, 50, help="isoTrack, sorted by pt"),
-        "genParticles" : NTupleCollection("genPart", genParticleWithMotherId, 200, help="all pruned genparticles"),
+        "genParticles" : NTupleCollection("genPart", genParticleWithMotherId, 300, help="all pruned genparticles"),
 })
 
 

CMGTools/TTHAnalysis/python/analyzers/ttHhistoCounterAnalyzer.py

@@ -0,0 +1,37 @@
+import operator 
+import itertools
+import copy
+from math import *
+
+from ROOT import std 
+from ROOT import TLorentzVector, TVectorD
+
+from PhysicsTools.Heppy.analyzers.core.Analyzer import Analyzer
+from PhysicsTools.HeppyCore.framework.event import Event
+from PhysicsTools.HeppyCore.statistics.counter import Counter, Counters
+from PhysicsTools.Heppy.analyzers.core.AutoHandle import AutoHandle
+
+import ROOT
+import os
+
+class ttHhistoCounterAnalyzer( Analyzer ):
+    def __init__(self, cfg_ana, cfg_comp, looperName ):
+        super(ttHhistoCounterAnalyzer,self).__init__(cfg_ana,cfg_comp,looperName) 
+
+    def declareHandles(self):
+        super(ttHhistoCounterAnalyzer, self).declareHandles()
+
+    def beginLoop(self, setup):
+        super(ttHhistoCounterAnalyzer,self).beginLoop(setup)
+        self.counters.addCounter('pairs')
+        count = self.counters.counter('pairs')
+        count.register('all events')
+        if "outputfile" in setup.services :
+            setup.services["outputfile"].file.cd()
+            self.inputCounter = ROOT.TH1F("Count","Count",1,0,2)
+
+    def process(self, event):
+        self.readCollections( event.input )
+        self.inputCounter.Fill(1)
+
+        return True

CMGTools/TTHAnalysis/python/plotter/bins/validation_plots_MT2.txt

@@ -28,18 +28,34 @@ count_nGammas20 : nGammas20 : [-0.5,0.5,1.5,2.5,3.5,4.5,5.5]; XTitle="N(#gamma\,
 count_nJet40 : nJet40 : [-0.5,0.5,1.5,2.5,3.5,4.5,5.5,6.5,7.5,8.5,9.5,10.5] ; XTitle='N(jet\, p_{T} > 40)'
 count_nBJet40 : nBJet40 : [-0.5,0.5,1.5,2.5,3.5,4.5,5.5,6.5] ; XTitle='N(b-jet\, p_{T} > 40\, CSVM)'
 
+count_ngenParticles: ngenPart : 500,0,500  ; XTitle="nGenParticles"
+
 ##### TRIGGERS
 
-HLT_ht350met100: HLT_ht350met100: [-0.5,0.5,1.5]; XTitle="HLT"
-HLT_MET150     : HLT_MET150     : [-0.5,0.5,1.5]; XTitle="HLT" 
-HLT_HT650      : HLT_HT650      : [-0.5,0.5,1.5]; XTitle="HLT"
+HLT_ht350met120: HLT_ht350met120: [-0.5,0.5,1.5]; XTitle="HLT"
+HLT_MET170     : HLT_MET170     : [-0.5,0.5,1.5]; XTitle="HLT" 
+HLT_HT900      : HLT_HT900      : [-0.5,0.5,1.5]; XTitle="HLT"
 HLT_SingleMu   : HLT_SingleMu   : [-0.5,0.5,1.5]; XTitle="HLT"
 HLT_DoubleMu   : HLT_DoubleMu   : [-0.5,0.5,1.5]; XTitle="HLT"
 HLT_DoubleEl   : HLT_DoubleEl   : [-0.5,0.5,1.5]; XTitle="HLT"
 HLT_MuEG       : HLT_MuEG       : [-0.5,0.5,1.5]; XTitle="HLT"
 HLT_Photons    : HLT_Photons    : [-0.5,0.5,1.5]; XTitle="HLT"
 ###HLT_htXprescale: HLT_htXprescale: [-0.5,0.5,1.5]; XTitle="HLT"
 
+Flag_EcalDeadCellTriggerPrimitiveFilter : Flag_EcalDeadCellTriggerPrimitiveFilter: [-0.5,0.5,1.5]; XTitle="metFilter"
+Flag_trkPOG_manystripclus53X : Flag_trkPOG_manystripclus53X: [-0.5,0.5,1.5]; XTitle="metFilter"
+Flag_ecalLaserCorrFilter : Flag_ecalLaserCorrFilter: [-0.5,0.5,1.5]; XTitle="metFilter"
+Flag_trkPOG_toomanystripclus53X : Flag_trkPOG_toomanystripclus53X : [-0.5,0.5,1.5]; XTitle="metFilter"
+Flag_hcalLaserEventFilter : Flag_hcalLaserEventFilter : [-0.5,0.5,1.5]; XTitle="metFilter"
+Flag_trkPOG_logErrorTooManyClusters : Flag_trkPOG_logErrorTooManyClusters : [-0.5,0.5,1.5]; XTitle="metFilter"
+Flag_trkPOGFilters : Flag_trkPOGFilters : [-0.5,0.5,1.5]; XTitle="metFilter"
+Flag_trackingFailureFilter : Flag_trackingFailureFilter: [-0.5,0.5,1.5]; XTitle="metFilter"
+Flag_CSCTightHaloFilter : Flag_CSCTightHaloFilter: [-0.5,0.5,1.5]; XTitle="metFilter"
+Flag_HBHENoiseFilter : Flag_HBHENoiseFilter : [-0.5,0.5,1.5]; XTitle="metFilter"
+Flag_goodVertices : Flag_goodVertices: [-0.5,0.5,1.5]; XTitle="metFilter"
+Flag_eeBadScFilter : Flag_eeBadScFilter: [-0.5,0.5,1.5]; XTitle="metFilter"
+Flag_METFilters : Flag_METFilters: [-0.5,0.5,1.5]; XTitle="metFilter"
+
 ##### MT2
 
 mt2:   mt2: 16,0,800    ; XTitle="MT2  [GeV]"
@@ -139,9 +155,9 @@ gamma_chHadIso: gamma_chHadIso: 100,0,0.5    ; XTitle="ch gamma ISO  [GeV]"
 gamma_neuHadIso: gamma_neuHadIso: 100,0,0.5    ; XTitle="neu gamma ISO  [GeV]"   
 gamma_phIso: gamma_phIso: 100,0,0.5    ; XTitle="ph gamma ISO  [GeV]"   
 
-gamma_sigmaIetaIeta: gamma_sigmaIetaIeta: 100,0,5.    ; XTitle="gamma sigmaIetaIeta"
-gamma_r9: gamma_r9: 100,0,5.    ; XTitle="gamma r9"
-gamma_hOverE: gamma_hOverE: 100,0,5.    ; XTitle="gamma hOverE"
+gamma_sigmaIetaIeta: gamma_sigmaIetaIeta: 100,0,0.1    ; XTitle="gamma sigmaIetaIeta"
+gamma_r9: gamma_r9: 100,0,1.    ; XTitle="gamma r9"
+gamma_hOverE: gamma_hOverE: 100,0,0.1    ; XTitle="gamma hOverE"
 gamma_idCutBased: gamma_idCutBased: 10,0,10.    ; XTitle="gamma ID Cut based"
 gamma_mcMatchId:  gamma_mcMatchId: 28,-1.5,26.5;  Legend='TR'
 

CMGTools/TTHAnalysis/python/samples/samples_13TeV_CSA14v2.py

@@ -61,10 +61,12 @@
 ZJetsToNuNu_HT100to200 = kreator.makeMCComponent("ZJetsToNuNu_HT100to200", "/ZJetsToNuNu_HT-100to200_Tune4C_13TeV-madgraph-tauola/Spring14miniaod-141029_PU40bx50_PLS170_V6AN2-v1/MINIAODSIM", "CMS", ".*root",372.6*1.27)
 ZJetsToNuNu_HT200to400 = kreator.makeMCComponent("ZJetsToNuNu_HT200to400", "/ZJetsToNuNu_HT-200to400_Tune4C_13TeV-madgraph-tauola/Spring14miniaod-141029_PU40bx50_PLS170_V6AN2-v1/MINIAODSIM", "CMS", ".*root",100.8*1.27)
 ZJetsToNuNu_HT400to600 = kreator.makeMCComponent("ZJetsToNuNu_HT400to600", "/ZJetsToNuNu_HT-400to600_Tune4C_13TeV-madgraph-tauola/Spring14miniaod-141029_PU40bx50_PLS170_V6AN2-v1/MINIAODSIM", "CMS", ".*root",11.99*1.27)
+ZJetsToNuNu_HT600toInf = kreator.makeMCComponent("ZJetsToNuNu_HT600toInf", "/ZJetsToNuNu_HT-600toInf_Tune4C_13TeV-madgraph-tauola/Spring14miniaod-141029_PU40bx50_PLS170_V6AN2-v1/MINIAODSIM", "CMS", ".*root",4.113*1.27)
 ZJetsToNuNuHT = [
 ZJetsToNuNu_HT100to200,
 ZJetsToNuNu_HT200to400,
 ZJetsToNuNu_HT400to600,
+ZJetsToNuNu_HT600toInf,
 ]
 WJetsToLNu_HT100to200 = kreator.makeMCComponent("WJetsToLNu_HT100to200", "/WJetsToLNu_HT-100to200_Tune4C_13TeV-madgraph-tauola/Spring14miniaod-141029_PU40bx50_PLS170_V6AN2-v1/MINIAODSIM", "CMS", ".*root",1817*1.23)
 WJetsToLNu_HT200to400 = kreator.makeMCComponent("WJetsToLNu_HT200to400", "/WJetsToLNu_HT-200to400_Tune4C_13TeV-madgraph-tauola/Spring14miniaod-141029_PU40bx50_PLS170_V6AN2-v1/MINIAODSIM", "CMS", ".*root",471.6*1.23)

CMGTools/TTHAnalysis/python/samples/samples_13TeV_PHYS14.py

@@ -48,7 +48,7 @@
 
 triggers_photon155 = ["HLT_Photon155_v*"]
 
-triggers_MT2_mumu = triggers_mumu_iso + triggers_mumu_noniso
+triggers_MT2_mumu = triggers_mumu_iso
 triggers_MT2_ee   = triggers_ee
 
 triggers_MT2_mue = triggers_mue
@@ -131,10 +131,12 @@
 ]
 ZJetsToNuNu_HT100to200 = kreator.makeMCComponent("ZJetsToNuNu_HT100to200", "/ZJetsToNuNu_HT-100to200_Tune4C_13TeV-madgraph-tauola/Phys14DR-PU20bx25_PHYS14_25_V1-v1/MINIAODSIM", "CMS", ".*root",372.6*1.27)
 ZJetsToNuNu_HT200to400 = kreator.makeMCComponent("ZJetsToNuNu_HT200to400", "/ZJetsToNuNu_HT-200to400_Tune4C_13TeV-madgraph-tauola/Phys14DR-PU20bx25_PHYS14_25_V1-v1/MINIAODSIM", "CMS", ".*root",100.8*1.27)
+ZJetsToNuNu_HT400to600 = kreator.makeMCComponent("ZJetsToNuNu_HT400to600", "/ZJetsToNuNu_HT-400to600_Tune4C_13TeV-madgraph-tauola/Phys14DR-PU20bx25_PHYS14_25_V1-v2/MINIAODSIM", "CMS", ".*root",11.99*1.27)
 ZJetsToNuNu_HT600toInf = kreator.makeMCComponent("ZJetsToNuNu_HT600toInf", "/ZJetsToNuNu_HT-600toInf_Tune4C_13TeV-madgraph-tauola/Phys14DR-PU20bx25_PHYS14_25_V1-v1/MINIAODSIM", "CMS", ".*root",4.113*1.27)
 ZJetsToNuNuHT = [
 ZJetsToNuNu_HT100to200,
 ZJetsToNuNu_HT200to400,
+ZJetsToNuNu_HT400to600,
 ZJetsToNuNu_HT600toInf,
 ]