Merge pull request cms-sw#213 from cbotta/CMGTools-from-CMSSW_7_2_3

xml for PHYS14 ele mva ID, xml for PHYS14 v2 susy and ttH lepMVA ID, cleaned 2lss cut files

CMGTools/TTHAnalysis/macros/leptons/trainLeptonID.cxx

@@ -2,8 +2,7 @@ void trainLeptonID(TString name, TString train="GoodvsBad") {
     TTree *dSig = (TTree*) _file0->Get("tree");
     TTree *dBg1 = (TTree*) _file1->Get("tree");
     TFile *fOut = new TFile(name+".root","RECREATE");
-    TMVA::Factory *factory = new TMVA::Factory(name, fOut, "!V:!Color");
-
+    TMVA::Factory *factory = new TMVA::Factory(name, fOut, "!V:!Color:Transformations=I");
 
     if (!name.Contains("pteta")) {
         factory->AddSpectator("LepGood_pt", 'D');

CMGTools/TTHAnalysis/python/analyzers/ntupleTypes.py

@@ -14,6 +14,7 @@
     NTupleVariable("mvaId",         lambda lepton : lepton.mvaNonTrigV0(full5x5=True) if abs(lepton.pdgId()) == 11 else lepton.mvaId(), help="EGamma POG MVA ID for non-triggering electrons (as HZZ); MVA Id for muons (BPH+Calo+Trk variables)"),
     NTupleVariable("mvaIdTrig",     lambda lepton : lepton.mvaTrigV0(full5x5=True)    if abs(lepton.pdgId()) == 11 else 1, help="EGamma POG MVA ID for triggering electrons; 1 for muons"),
    # Lepton MVA-id related variables
+    NTupleVariable("mvaTTH",    lambda lepton : getattr(lepton, 'mvaValueTTH', -1), help="Lepton MVA (TTH version)"),
     NTupleVariable("mvaSusy",    lambda lepton : getattr(lepton, 'mvaValueSusy', -1), help="Lepton MVA (SUSY version)"),
     NTupleVariable("jetPtRatio", lambda lepton : lepton.pt()/lepton.jet.pt() if hasattr(lepton,'jet') else -1, help="pt(lepton)/pt(nearest jet)"),
     NTupleVariable("jetBTagCSV", lambda lepton : lepton.jet.btag('combinedInclusiveSecondaryVertexV2BJetTags') if hasattr(lepton,'jet') and hasattr(lepton.jet, 'btag') else -99, help="CSV btag of nearest jet"),

CMGTools/TTHAnalysis/python/analyzers/ttHCoreEventAnalyzer.py

@@ -10,7 +10,7 @@ class ttHCoreEventAnalyzer( Analyzer ):
     def __init__(self, cfg_ana, cfg_comp, looperName ):
         super(ttHCoreEventAnalyzer,self).__init__(cfg_ana,cfg_comp,looperName)
         self.maxLeps = cfg_ana.maxLeps
-        #self.leptonMVATTH  = LeptonMVA("TTH", "%s/src/CMGTools/TTHAnalysis/data/leptonMVA/%%s_BDTG.weights.xml" % os.environ['CMSSW_BASE'], self.cfg_comp.isMC)
+        self.leptonMVATTH  = LeptonMVA("Susy", "%s/src/CMGTools/TTHAnalysis/data/leptonMVA/tth/%%s_BDTG.weights.xml" % os.environ['CMSSW_BASE'], self.cfg_comp.isMC)
         self.leptonMVASusy = LeptonMVA("Susy","%s/src/CMGTools/TTHAnalysis/data/leptonMVA/susy/%%s_BDTG.weights.xml" % os.environ['CMSSW_BASE'], self.cfg_comp.isMC)
 
     def declareHandles(self):
@@ -284,11 +284,11 @@ def process(self, event):
             if thisDeltaPhi < event.deltaPhiMin : event.deltaPhiMin = thisDeltaPhi
 
         for lep in event.selectedLeptons:
-            #lep.mvaValue     = self.leptonMVATTH(lep)
+            lep.mvaValueTTH     = self.leptonMVATTH(lep)
             lep.mvaValueSusy = self.leptonMVASusy(lep)
         for lep in event.inclusiveLeptons:
             if lep not in event.selectedLeptons:
-                #lep.mvaValue     = self.leptonMVATTH(lep)
+                lep.mvaValueTTH     = self.leptonMVATTH(lep)
                 lep.mvaValueSusy = self.leptonMVASusy(lep)
 
 

CMGTools/TTHAnalysis/python/plotter/bins/susy_2lss_sync.sh → CMGTools/TTHAnalysis/python/plotter/bins/susymultilepton/susy_2lss_sync.sh

@@ -1,12 +1,18 @@
 #!/bin/bash
-T="~/w/TREES_72X_SYNC"
+
+#T="~/w/TREES_72X_SYNC"
+#CORE="-P $T --s2v --tree treeProducerSusyMultilepton "
+#CORE="${CORE} -F sf/t $T/1_susyVars_2lssInc_v0/evVarFriend_{cname}.root -X exclusive --mcc bins/susy_2lssinc_lepchoice_sync.txt"
+T=" /data/g/gpetrucc/TREES_72X_040115/9_skim_2lssSR_sync_v1"
 CORE="-P $T --s2v --tree treeProducerSusyMultilepton "
-CORE="${CORE} -F sf/t $T/1_susyVars_2lssInc_v0/evVarFriend_{cname}.root -X exclusive --mcc bins/susy_2lssinc_lepchoice_sync.txt"
+CORE="${CORE} -F sf/t $T/0_allfriends/evVarFriend_{cname}.root -X exclusive --mcc bins/susy_2lssinc_lepchoice_sync.txt"
 
 POST="";
 if [[ "$1" == "mccounts" ]]; then
     GO="python mcAnalysis.py $CORE mca-Phys14.txt bins/susy_2lss_sync.txt -p T1tttt_HM -f -G -u "
     POST="| awk '/all/{print \$2}' "
+elif [[ "$1" == "mcyields" ]]; then
+    GO="python mcAnalysis.py $CORE mca-Phys14.txt bins/susy_2lss_sync.txt -p 'T1tttt_HM,T5tttt_MLDx,T5qqqqWW_HM,T5qqqqWWD,TTX,WZ,TT,WJets,DY' --pgroup TTX=TT[WZH] -f -G -l 4.0"
 elif [[ "$1" == "mcdumps" ]]; then
     FMT='{run:1d} {lumi:9d} {evt:12d}\t{nLepGood10:2d}\t{LepGood1_pdgId:+2d} {LepGood1_pt:5.1f}\t{LepGood2_pdgId:+2d} {LepGood2_pt:5.1f}\t{nJet40}\t{nBJetMedium40:2d}\t{met_pt:5.1f}\t{htJet40j:6.1f}'
     python mcDump.py $CORE  mca-Phys14.txt bins/susy_2lss_sync.txt -p T1tttt_HM -X lep1_pt25 -X lep2_pt25   | sort -n -k1 -k2 > 2lssInc_all.txt
@@ -20,14 +26,20 @@ fi
 
 SAVE="${GO}"
 for LL  in ee em mm; do 
-for SR  in 0 10 20 30; do 
-for LPt in hh hl; do
+for SR  in 00 10 20 30; do 
+for LPt in hh hl ll; do
 
 GO="${SAVE}"
 case $SR in
-1[0-9])  GO="${GO} -R nBjet nBjet1 'nBJetMedium40 == 1' " ;;
-2[0-9])  GO="${GO} -R nBjet nBjet2 'nBJetMedium40 == 2' " ;;
-3[0-9])  GO="${GO} -R nBjet nBjet3 'nBJetMedium40 >= 3' " ;;
+0)   GO="${GO} -R nBjet nBjet0 'nBJetMedium25 >= 0' " ;;
+00)  GO="${GO} -R nBjet nBjet0 'nBJetMedium25 == 0' " ;;
+10)  GO="${GO} -R nBjet nBjet1 'nBJetMedium25 == 1' " ;;
+20)  GO="${GO} -R nBjet nBjet2 'nBJetMedium25 == 2' " ;;
+30)  GO="${GO} -R nBjet nBjet3 'nBJetMedium25 >= 3' " ;;
+0[1-9X])  GO="${GO} -R nBjet nBjet0 'nBJetMedium25 == 0' -R met metSR 'met_pt > 50' -R ht htSR 'htJet40j > 200'" ;;
+1[1-9X])  GO="${GO} -R nBjet nBjet1 'nBJetMedium25 == 1' -R met metSR 'met_pt > 50' -R ht htSR 'htJet40j > 200'" ;;
+2[1-9X])  GO="${GO} -R nBjet nBjet2 'nBJetMedium25 == 2' -R met metSR 'met_pt > 50' -R ht htSR 'htJet40j > 200'" ;;
+3[1-9X])  GO="${GO} -R nBjet nBjet3 'nBJetMedium25 >= 3' -R met metSR 'met_pt > 50' -R ht htSR 'htJet40j > 200'" ;;
 esac;
 case $LL in
 ee)  GO="${GO} -R anyll ee 'abs(LepGood1_pdgId) == 11 && abs(LepGood2_pdgId) == 11' " ;;

CMGTools/TTHAnalysis/python/plotter/bins/susy_2lss_sync.txt → CMGTools/TTHAnalysis/python/plotter/bins/susymultilepton/susy_2lss_sync.txt

@@ -14,4 +14,4 @@ ele cuts: (abs(LepGood1_pdgId) == 13 || (LepGood1_convVeto && LepGood1_lostHits
 ht: htJet40j > 80
 met: met_pt > 30 || htJet40j > 500
 njet: nJet40 >=2
-nBjet: nBJetMedium40 >= 0
+nBjet: nBJetMedium25 >= 0