Skip to content

Commit

Permalink
Run ID_V2 on 2017 data,MC
Browse files Browse the repository at this point in the history
  • Loading branch information
@swagata87
swagata87 committed Aug 30, 2018
1 parent 9ecdfd0 commit ec65ca6
Show file tree
Hide file tree
Showing 5 changed files with 71 additions and 28 deletions.
39 changes: 24 additions & 15 deletions python/TnPTreeProducer_cfg.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -8,8 +8,8 @@
## argument line options
###################################################################
varOptions = VarParsing('analysis')

varOptions.register(
#"isMC", True,
"isMC", False,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
Expand Down Expand Up @@ -39,7 +39,6 @@

varOptions.register(
"doRECO", False,
# "doRECO", True,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"Include tree for Reco SF"
Expand All @@ -49,18 +48,9 @@
"calibEn", False,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,

"use EGM smearer to calibrate photon and electron energy"
)

varOptions.register(
"isAOD", False,
# "isAOD", True,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"switch to run other AOD (for RECO SFs)"
)

#### HLTname is HLT2 in reHLT samples
varOptions.register(
"HLTname", "HLT",
Expand All @@ -71,20 +61,28 @@

varOptions.register(
#"GT","auto",
"GT","101X_dataRun2_Prompt_v9",
#"GT","101X_dataRun2_Prompt_v9",
"GT","94X_dataRun2_ReReco_EOY17_v6",
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Global Tag to be used"
)

varOptions.register(
"includeSUSY", True,
# "isAOD", True,
"includeSUSY", False,
# "isAOD", True,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"add also the variables used by SUSY"
)

varOptions.register(
"isAOD", True,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"use AOD"
)

varOptions.parseArguments()


Expand Down Expand Up @@ -162,7 +160,8 @@
###################################################################
from EgammaAnalysis.TnPTreeProducer.etc.tnpInputTestFiles_cff import filesMiniAOD_Preliminary2018 as inputs
#if options['useAOD'] : from EgammaAnalysis.TnPTreeProducer.etc.tnpInputTestFiles_cff import filesAOD_23Sep2016 as inputs #switch to 2017 samples if want to cmsRun on AOD
if options['useAOD'] : from EgammaAnalysis.TnPTreeProducer.etc.tnpInputTestFiles_cff import filesAOD_Preliminary2018 as inputs #switch to 2017 samples if want to cmsRun on AOD
#if options['useAOD'] :
from EgammaAnalysis.TnPTreeProducer.etc.tnpInputTestFiles_cff import filesAOD_Preliminary2017 as inputs #switch to 2017 samples if want to cmsRun on AOD
#if options['useAOD'] : from EgammaAnalysis.TnPTreeProducer.etc.tnpInputTestFiles_cff import filesAOD_empty as inputs #switch to 2017 samples if want to cmsRun on AOD

options['INPUT_FILE_NAME'] = inputs['data']
Expand Down Expand Up @@ -285,6 +284,11 @@
passingMedium94X = cms.InputTag("probeEleCutBasedMedium94X"),
passingTight94X = cms.InputTag("probeEleCutBasedTight94X" ),

passingVeto94XV2 = cms.InputTag("probeEleCutBasedVeto94XV2" ),
passingLoose94XV2 = cms.InputTag("probeEleCutBasedLoose94XV2" ),
passingMedium94XV2 = cms.InputTag("probeEleCutBasedMedium94XV2"),
passingTight94XV2 = cms.InputTag("probeEleCutBasedTight94XV2" ),

passingMVA94XwpLnoiso = cms.InputTag("probeEleMVA94XwpLnoiso" ),
passingMVA94Xwp90noiso = cms.InputTag("probeEleMVA94Xwp90noiso" ),
passingMVA94Xwp80noiso = cms.InputTag("probeEleMVA94Xwp80noiso" ),
Expand Down Expand Up @@ -321,6 +325,11 @@
passingLoose94X = cms.InputTag("probePhoCutBasedLoose94X"),
passingMedium94X = cms.InputTag("probePhoCutBasedMedium94X"),
passingTight94X = cms.InputTag("probePhoCutBasedTight94X"),

passingLoose100XV2 = cms.InputTag("probePhoCutBasedLoose100XV2"),
passingMedium100XV2 = cms.InputTag("probePhoCutBasedMedium100XV2"),
passingTight100XV2 = cms.InputTag("probePhoCutBasedTight100XV2"),

passingMVA94Xwp90 = cms.InputTag("probePhoMVA94Xwp90"),
passingMVA94Xwp80 = cms.InputTag("probePhoMVA94Xwp80"),
)
Expand Down
35 changes: 23 additions & 12 deletions python/egmElectronIDModules_cff.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -20,12 +20,13 @@ def setIDs(process, options):
my_id_modules = [
'RecoEgamma.ElectronIdentification.Identification.cutBasedElectronID_Summer16_80X_V1_cff',
'RecoEgamma.ElectronIdentification.Identification.heepElectronID_HEEPV60_cff',
'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_Spring16_GeneralPurpose_V1_cff',
'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_Spring16_GeneralPurpose_V1_cff',
'RecoEgamma.ElectronIdentification.Identification.cutBasedElectronID_Fall17_94X_V1_cff',
'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_Fall17_noIso_V1_cff',
'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_Fall17_iso_V1_cff',
# 'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_Fall17_noIso_V2_cff',
# 'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_Fall17_iso_V2_cff',
'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_Fall17_noIso_V2_cff',
'RecoEgamma.ElectronIdentification.Identification.mvaElectronID_Fall17_iso_V2_cff',
'RecoEgamma.ElectronIdentification.Identification.cutBasedElectronID_Fall17_94X_V2_cff'
]

### add only miniAOD supported IDs
Expand Down Expand Up @@ -120,23 +121,33 @@ def setIDs(process, options):
process.probeEleMVA94Xwp80isoV2 = process.probeEleCutBasedVeto.clone()
process.probeEleMVA94XwpHZZisoV2 = process.probeEleCutBasedVeto.clone()

# process.probeEleMVA94XwpLnoisoV2.selection = cms.InputTag("egmGsfElectronIDs:mvaEleID-Fall17-noIso-V2-wpLoose" )
# process.probeEleMVA94Xwp90noisoV2.selection = cms.InputTag("egmGsfElectronIDs:mvaEleID-Fall17-noIso-V2-wp90" )
# process.probeEleMVA94Xwp80noisoV2.selection = cms.InputTag("egmGsfElectronIDs:mvaEleID-Fall17-noIso-V2-wp80" )
# process.probeEleMVA94XwpLisoV2.selection = cms.InputTag("egmGsfElectronIDs:mvaEleID-Fall17-iso-V2-wpLoose" )
# process.probeEleMVA94Xwp90isoV2.selection = cms.InputTag("egmGsfElectronIDs:mvaEleID-Fall17-iso-V2-wp90" )
# process.probeEleMVA94Xwp80isoV2.selection = cms.InputTag("egmGsfElectronIDs:mvaEleID-Fall17-iso-V2-wp80" )
# process.probeEleMVA94XwpHZZisoV2.selection = cms.InputTag("egmGsfElectronIDs:mvaEleID-Fall17-iso-V2-wpHZZ" )
process.probeEleMVA94XwpLnoisoV2.selection = cms.InputTag("egmGsfElectronIDs:mvaEleID-Fall17-noIso-V2-wpLoose" )
process.probeEleMVA94Xwp90noisoV2.selection = cms.InputTag("egmGsfElectronIDs:mvaEleID-Fall17-noIso-V2-wp90" )
process.probeEleMVA94Xwp80noisoV2.selection = cms.InputTag("egmGsfElectronIDs:mvaEleID-Fall17-noIso-V2-wp80" )
process.probeEleMVA94XwpLisoV2.selection = cms.InputTag("egmGsfElectronIDs:mvaEleID-Fall17-iso-V2-wpLoose" )
process.probeEleMVA94Xwp90isoV2.selection = cms.InputTag("egmGsfElectronIDs:mvaEleID-Fall17-iso-V2-wp90" )
process.probeEleMVA94Xwp80isoV2.selection = cms.InputTag("egmGsfElectronIDs:mvaEleID-Fall17-iso-V2-wp80" )
process.probeEleMVA94XwpHZZisoV2.selection = cms.InputTag("egmGsfElectronIDs:mvaEleID-Fall17-iso-V2-wpHZZ" )


process.probeEleCutBasedVeto94XV2 = process.probeEleCutBasedVeto.clone()
process.probeEleCutBasedLoose94XV2 = process.probeEleCutBasedVeto.clone()
process.probeEleCutBasedMedium94XV2 = process.probeEleCutBasedVeto.clone()
process.probeEleCutBasedTight94XV2 = process.probeEleCutBasedVeto.clone()

process.probeEleCutBasedVeto94XV2.selection = cms.InputTag("egmGsfElectronIDs:cutBasedElectronID-Fall17-94X-V2-veto" )
process.probeEleCutBasedLoose94XV2.selection = cms.InputTag("egmGsfElectronIDs:cutBasedElectronID-Fall17-94X-V2-loose" )
process.probeEleCutBasedMedium94XV2.selection = cms.InputTag("egmGsfElectronIDs:cutBasedElectronID-Fall17-94X-V2-medium" )
process.probeEleCutBasedTight94XV2.selection = cms.InputTag("egmGsfElectronIDs:cutBasedElectronID-Fall17-94X-V2-tight" )


process.tagEleCutBasedTight = cms.EDProducer(eleProducer,
input = cms.InputTag("goodElectrons"),
cut = cms.string(options['ELECTRON_TAG_CUTS']),
selection = cms.InputTag("egmGsfElectronIDs:cutBasedElectronID-Fall17-94X-V1-tight"),
selection = cms.InputTag("egmGsfElectronIDs:cutBasedElectronID-Fall17-94X-V2-tight"),
id_cut = cms.bool(True)
)
process.tagEleCutBasedTight.selection = cms.InputTag("egmGsfElectronIDs:cutBasedElectronID-Fall17-94X-V1-tight")
process.tagEleCutBasedTight.selection = cms.InputTag("egmGsfElectronIDs:cutBasedElectronID-Fall17-94X-V2-tight")

if options['addSUSY'] :

Expand Down
7 changes: 7 additions & 0 deletions python/egmPhotonIDModules_cff.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -21,6 +21,7 @@ def setIDs(process, options):
'RecoEgamma.PhotonIdentification.Identification.mvaPhotonID_Spring16_nonTrig_V1_cff',
'RecoEgamma.PhotonIdentification.Identification.cutBasedPhotonID_Fall17_94X_V1_cff',
'RecoEgamma.PhotonIdentification.Identification.mvaPhotonID_Fall17_94X_V1_cff',
'RecoEgamma.PhotonIdentification.Identification.cutBasedPhotonID_Fall17_100X_V2_cff'
]

for idmod in my_id_modules:
Expand Down Expand Up @@ -58,6 +59,12 @@ def setIDs(process, options):
process.probePhoCutBasedMedium94X.selection = cms.InputTag("egmPhotonIDs:cutBasedPhotonID-Fall17-94X-V1-medium" )
process.probePhoCutBasedTight94X.selection = cms.InputTag("egmPhotonIDs:cutBasedPhotonID-Fall17-94X-V1-tight" )

process.probePhoCutBasedLoose100XV2 = process.probePhoCutBasedLoose.clone()
process.probePhoCutBasedMedium100XV2 = process.probePhoCutBasedLoose.clone()
process.probePhoCutBasedTight100XV2 = process.probePhoCutBasedLoose.clone()
process.probePhoCutBasedLoose100XV2.selection = cms.InputTag("egmPhotonIDs:cutBasedPhotonID-Fall17-100X-V2-loose" )
process.probePhoCutBasedMedium100XV2.selection = cms.InputTag("egmPhotonIDs:cutBasedPhotonID-Fall17-100X-V2-medium" )
process.probePhoCutBasedTight100XV2.selection = cms.InputTag("egmPhotonIDs:cutBasedPhotonID-Fall17-100X-V2-tight" )

process.probePhoMVA = process.probePhoCutBasedLoose.clone()
process.probePhoMVA.selection = cms.InputTag("egmPhotonIDs:mvaPhoID-Spring16-nonTrig-V1-wp90")
Expand Down
7 changes: 7 additions & 0 deletions python/egmTreesSetup_cff.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -184,6 +184,10 @@ def setSequences(process, options):
process.probeEleCutBasedLoose94X +
process.probeEleCutBasedMedium94X +
process.probeEleCutBasedTight94X +
process.probeEleCutBasedVeto94XV2 +
process.probeEleCutBasedLoose94XV2 +
process.probeEleCutBasedMedium94XV2 +
process.probeEleCutBasedTight94XV2 +
process.probeEleMVA94XwpLnoiso +
process.probeEleMVA94Xwp90noiso +
process.probeEleMVA94Xwp80noiso +
Expand Down Expand Up @@ -216,6 +220,9 @@ def setSequences(process, options):
process.probePhoCutBasedLoose94X +
process.probePhoCutBasedMedium94X +
process.probePhoCutBasedTight94X +
process.probePhoCutBasedLoose100XV2 +
process.probePhoCutBasedMedium100XV2 +
process.probePhoCutBasedTight100XV2 +
process.probePhoMVA94Xwp90 +
process.probePhoMVA94Xwp80 +
process.probePho
Expand Down
11 changes: 10 additions & 1 deletion python/etc/tnpInputTestFiles_cff.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -15,13 +15,22 @@

filesAOD_Preliminary2018 = {
'mc' : cms.untracked.vstring(
'/store/mc/RunIISpring18DRPremix/DYJetsToLL_M-50_TuneCP5_13TeV-madgraphMLM-pythia8/AODSIM/100X_upgrade2018_realistic_v10-v2/100000/968B1BAD-332F-E811-8178-FA163E9F7460.root',
# '/store/mc/RunIISpring18DRPremix/DYJetsToLL_M-50_TuneCP5_13TeV-madgraphMLM-pythia8/AODSIM/100X_upgrade2018_realistic_v10-v2/90001/FE8F7D45-133E-E811-891E-FA163EA4957D.root',
),
'data' : cms.untracked.vstring(
'/store/data/Run2018B/EGamma/AOD/PromptReco-v1/000/317/864/00000/C269C719-EC71-E811-9C7E-FA163EF55202.root',
)
}

filesAOD_Preliminary2017 = {
'mc' : cms.untracked.vstring(
'/store/mc/RunIIFall17DRPremix/DYJetsToLL_M-50_TuneCP5_13TeV-amcatnloFXFX-pythia8/AODSIM/PU2017RECOPF_94X_mc2017_realistic_v11-v1/50000/D42B8057-9F67-E811-9656-549F3525C4EC.root'
),
'data' : cms.untracked.vstring(
'/store/data/Run2017F/SingleElectron/AOD/17Nov2017-v1/50000/005B2A56-96E0-E711-B727-0CC47A4D7690.root',
)
}


filesMiniAOD_Preliminary2017 = {
'mc' : cms.untracked.vstring(
Expand Down

0 comments on commit ec65ca6

Please sign in to comment.