Merge pull request #1 from jsturdy/simplified-workflow

Simplified workflow

Author: jsturdy

.gitignore

@@ -0,0 +1,21 @@
+# Ignore everything
+*
+# Don't ignore subdirectories
+!*/
+# Don't ignore special files
+!*.gitignore
+!*.md
+!*.rst
+!LICENSE
+# Don't ignore source files, etc
+!*.C
+!*.cc
+!*.h
+!*.py
+!*.sh
+# Don't ignore build files
+!BuildFile.xml
+# Ignore __init__.py files
+__init__.py
+*.#*
+# Others to ignore

Analyzer/plugins/MakeZprimeMiniAodTreeHEEP.cc

@@ -1,18 +1,14 @@
-#==============================================================
-#          Analysis code for Z' boson to dilepton analysis    =  
-#          This python to run the tree producer code          = 
-#          To run over MINIAOD MC or Data with fixed trigger  = 
-#                  Author:  Sherif Elgammal                   = 
-#                       01/01/2016                            = 
-#==============================================================
 import FWCore.ParameterSet.Config as cms
+from FWCore.ParameterSet.VarParsing import VarParsing
 from RecoMuon.TrackingTools.MuonServiceProxy_cff import *
 from PhysicsTools.PatAlgos.triggerLayer1.triggerProducer_cfi import *
 #from TrackingTools.TransientTrack.TransientTrackBuilder_cfi import *
+from PhysicsTools.SelectorUtils.tools.vid_id_tools import switchOnVIDElectronIdProducer
+from PhysicsTools.SelectorUtils.tools.vid_id_tools import setupVIDElectronSelection
+from PhysicsTools.SelectorUtils.tools.vid_id_tools import setupAllVIDIdsInModule
+from PhysicsTools.SelectorUtils.tools.vid_id_tools import DataFormat
 
 process = cms.Process("tree")
-process.load("RecoEgamma.ElectronIdentification.heepIdVarValueMapProducer_cfi")
-process.load("TrackingTools/TransientTrack/TransientTrackBuilder_cfi")
 process.load("Configuration.StandardSequences.MagneticField_38T_cff")
 process.load('Configuration/Geometry/GeometryRecoDB_cff')
 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff')
@@ -24,17 +20,25 @@
 process.load("FWCore.MessageLogger.MessageLogger_cfi")
 #Track isolation correction
 process.load("TrackingTools.TransientTrack.TransientTrackBuilder_cfi")
-process.load("PhysicsTools.PatAlgos.slimming.packedCandidatesForTrkIso_cfi")
 process.load("PhysicsTools.PatAlgos.slimming.primaryVertexAssociation_cfi")
 process.load("RecoEgamma.ElectronIdentification.heepIdVarValueMapProducer_cfi")
 
+
+#EGamma VID
+from PhysicsTools.SelectorUtils.tools.vid_id_tools import *
+dataFormat = DataFormat.MiniAOD
+switchOnVIDElectronIdProducer(process, dataFormat)
+
+
+
+
 process.options = cms.untracked.PSet(
     Rethrow = cms.untracked.vstring('ProductNotFound'),
     wantSummary = cms.untracked.bool(True) 
 )
 
 process.maxEvents = cms.untracked.PSet(
-    input = cms.untracked.int32(-1),
+    input = cms.untracked.int32(10000),
     SkipEvent = cms.untracked.vstring('ProductNotFound')
 )
 
@@ -43,7 +47,15 @@
 
 process.source = cms.Source("PoolSource",
     fileNames = cms.untracked.vstring(
-'/store/mc/RunIISpring16MiniAODv1/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/MINIAODSIM/PUSpring16_80X_mcRun2_asymptotic_2016_v3_ext1-v1/20000/0017320C-7BFC-E511-9B2D-0CC47A4C8E34.root'
+
+'/store/mc/RunIISpring16MiniAODv2/ZToEE_NNPDF30_13TeV-powheg_M_200_400/MINIAODSIM/PUSpring16RAWAODSIM_reHLT_80X_mcRun2_asymptotic_v14-v1/20000/0C242E56-BA3A-E611-9B2D-0242AC130004.root',
+
+
+
+#'/store/data/Run2016B/DoubleEG/MINIAOD/03Feb2017_ver2-v2/50000/00054938-CEEA-E611-889E-0CC47A4D7650.root',
+#'/store/data/Run2016B/DoubleEG/MINIAOD/03Feb2017_ver2-v2/50000/001BC16F-90EA-E611-87E2-F04DA27540BB.root',
+#'/store/data/Run2016B/DoubleEG/MINIAOD/03Feb2017_ver2-v2/50000/001EB4EF-D3EA-E611-B94E-0CC47A4C8F26.root'
+
 #'/store/data/Run2016G/SingleMuon/MINIAOD/23Sep2016-v1/90000/9695BF0B-8E97-E611-A1A4-00259044051C.root',
 #'/store/data/Run2016G/SingleMuon/MINIAOD/23Sep2016-v1/90000/1C603D88-C597-E611-A0BE-008CFAF2931E.root'
 	)
@@ -64,22 +76,29 @@
                 thresh = cms.untracked.double(0.25)
                 )
 ##################################################################
+#from HEEP.VID.tools import addHEEPV70ElesMiniAOD
+#addHEEPV70ElesMiniAOD(process,useStdName=True)
+
+from HEEP.VID.tools import addHEEPV70ElesMiniAOD
+addHEEPV70ElesMiniAOD(process,useStdName=True)
+
 #####################################################################
 # Global tag (data)
-#from Configuration.AlCa.GlobalTag_condDBv2 import GlobalTag
-#process.GlobalTag = GlobalTag(process.GlobalTag, '80X_dataRun2_2016SeptRepro_v4', '')
+from Configuration.AlCa.GlobalTag_condDBv2 import GlobalTag
+process.GlobalTag = GlobalTag(process.GlobalTag, '80X_dataRun2_2016SeptRepro_v7', '')
 
 # Global tag (MC)
-from Configuration.AlCa.GlobalTag_condDBv2 import GlobalTag
-process.GlobalTag = GlobalTag(process.GlobalTag, '80X_mcRun2_asymptotic_2016_v3', '')
+#from Configuration.AlCa.GlobalTag_condDBv2 import GlobalTag
+#process.GlobalTag = GlobalTag(process.GlobalTag, '80X_mcRun2_asymptotic_2016_v3', '')
 
-process.demo = cms.EDAnalyzer("MakeZprimeMiniAodTree",
-    outputFile = cms.string('MC.root'),
+process.demo = cms.EDAnalyzer("MakeZprimeMiniAodTreeHEEPData",
+    outputFile = cms.string('Data.root'),
     eleTrkPtIsoLabel = cms.InputTag("heepIDVarValueMaps","eleTrkPtIso"),
     scProducer = cms.InputTag("reducedEgamma:reducedSuperClusters"),
     vertices   = cms.InputTag("offlineSlimmedPrimaryVertices"),
     muons      = cms.InputTag("slimmedMuons"),
     electrons  = cms.InputTag("slimmedElectrons"),
+    eles       = cms.InputTag("slimmedElectrons"),
     taus       = cms.InputTag("slimmedTaus"),
     photons    = cms.InputTag("slimmedPhotons"),
     #jets       = cms.InputTag("slimmedJetsCMSTopTagCHSPacked:SubJets"),
@@ -123,6 +142,7 @@
         #'pfCombinedMVABJetTags'
     ),
 
+    #Analysis = cms.string('ZprimeToEE')
     Analysis = cms.string('ZprimeToMuMu')
 )
 

Analyzer/plugins/MakeZprimeMiniAodTreeMC.cc

@@ -0,0 +1,13 @@
+# `ZprimeDiLeptons`
+This package contains several pieces of code
+* `Analyzer`: Contains the code for creating the `ntuples` used for making plots
+
+## Working with the `ZprimeDiLeptons` `ntuple` making code
+
+### Obtaining the code and setting up a working release area
+(Currently tested and working with `CMSSW_8_0_25`, which is checked out as part of the script)
+```
+wget https://raw.githubusercontent.com/wsu-fnal-ci-analysis/ZprimeDiLeptons/<tag>/scripts/setup.sh
+cd <location where you want to set up the release>
+source setup.sh
+```

Analyzer/plugins/MakeZprimeMiniAodTreeMuon.cc

@@ -0,0 +1,62 @@
+#!/bin/sh
+
+# usage() {
+#     echo 'Usage : ./setup.sh -r<ROOT version> [-s<sub version for ROOT6>] -c<compiler>'
+#     echo '         -r must be either 5 or 6'
+#     echo '         -s must be either 6, 8, or 10 if -r is 6'
+#     echo '         -c must be either gcc or clang'
+#     kill -INT $$;
+# }
+
+# while getopts "r:c:s:h" o; do
+#     case "${o}" in
+#         r)
+#             local rootver=${OPTARG}
+#             if [ "${rootver}" != "5" ] && [ "${rootver}" != "6" ]
+#             then
+#                 echo "Invalid ROOT version specified\n"
+#                 usage
+#             fi
+#             ;;
+#         s)
+#             local subver=${OPTARG}
+#             ;;
+#         c)
+#             local compiler=${OPTARG}
+#             if [ "${compiler}" != "gcc" ] && [ "${compiler}" != "clang" ]
+#             then
+#                 echo "Invalid compiler option\n"
+#                 usage
+#             fi
+#             ;;
+#         h)
+#             echo "Printing help"
+#             usage
+#             ;;
+#         *)
+#             echo "Unknown option"
+#             usage
+#             ;;
+#     esac
+# done
+
+# shift $((OPTIND-1))
+
+echo cmsrel CMSSW_8_0_25
+echo cd CMSSW_8_0_25/src
+echo cmsenv
+echo git cms-init
+# brings in HEEP V70 into VID
+echo git cms-merge-topic Sam-Harper:HEEPV70VID_8010_ReducedCheckout 
+# for other E/gamma IDs in VID if you wish to have them
+echo git cms-merge-topic ikrav:egm_id_80X_v3
+# only necessary to run HEEP V70 on AOD (it will crash if this is not present looking for puppi candidates
+echo git cms-merge-topic Sam-Harper:PackedCandNoPuppi
+# we need this for the mva weights which runs in VID regardless if you need it or not
+echo mkdir -p ../external/slc6_amd64_gcc530/data/RecoEgamma/ElectronIdentification/
+echo git clone git@github.com:cms-data/RecoEgamma-ElectronIdentification ../external/slc6_amd64_gcc530/data/RecoEgamma/ElectronIdentification/data
+# needed for HEEP modified electrons with value maps
+echo git clone -b HEEPV70 git@github.com:Sam-Harper/HEEP.git
+echo git cms-merge-topic Sam-Harper:NewEGModifiers_8010
+# the zprime selection code
+echo git clone https://github.com/cms-analysis/ZprimeDiLeptons