CMGWPrimeGroup/python/selectors_cff.py

import FWCore.ParameterSet.Config as cms

#####################
######  Muons  ######
#####################

#Note: These cuts are <= and >= 
muonSelectors = cms.PSet(
    VBTF = cms.PSet(
       minPt = cms.untracked.double(10.),
       maxEta = cms.untracked.double(2.4),
       minIsGlobal = cms.untracked.int32(1),
       maxDxy = cms.untracked.double(0.2),
       maxNormalizedChi2 = cms.untracked.double(10.0),
       minNTrackerHits = cms.untracked.int32(11),
       minNPixelHits = cms.untracked.int32(1),
       minNMuonHits = cms.untracked.int32(1),
       minNMatches = cms.untracked.int32(2),
    ),
    PFIsoHighPtBoostedZRelaxed = cms.PSet(
       minPt = cms.untracked.double(10.),
       maxEta = cms.untracked.double(2.4),
       minIsGblOrTrk = cms.untracked.int32(1),
       maxDxy = cms.untracked.double(0.2),
       maxDz = cms.untracked.double(0.5),
       maxRelPtErr = cms.untracked.double(0.3),
       minNPixelHits = cms.untracked.int32(1),
       minNTrackerLayers = cms.untracked.int32(6),
       minNMatches = cms.untracked.int32(2),
       #maxPFIso = cms.untracked.double(???),#Added below in clones
    ),
    PFIsoHighPtRelaxed = cms.PSet(
       minPt = cms.untracked.double(10.),
       maxEta = cms.untracked.double(2.4),
       minIsGlobal = cms.untracked.int32(1),
       maxDxy = cms.untracked.double(0.2),
       maxDz = cms.untracked.double(0.5),
       maxRelPtErr = cms.untracked.double(0.3),
       minNPixelHits = cms.untracked.int32(1),
       minNMuonHits = cms.untracked.int32(1), 
       minNTrackerLayers = cms.untracked.int32(6),
       minNMatches = cms.untracked.int32(2),
       #maxPFIso = cms.untracked.double(???),#Added below in clones
    ),
    PFIsoRelaxed = cms.PSet(
       minPt = cms.untracked.double(10.),
       maxEta = cms.untracked.double(2.4),
       minIsGlobal = cms.untracked.int32(1),
       minIsPF = cms.untracked.int32(1),
       maxDxy = cms.untracked.double(0.2),
       maxDz = cms.untracked.double(0.5),
       maxNormalizedChi2 = cms.untracked.double(10.0), 
       maxRelPtErr = cms.untracked.double(0.3),
       minNPixelHits = cms.untracked.int32(1),
       minNMuonHits = cms.untracked.int32(1), 
       minNTrackerLayers = cms.untracked.int32(6),
       minNMatches = cms.untracked.int32(2),
       #maxPFIso = cms.untracked.double(???),#Added below in clones
    ),
    Reco = cms.PSet(
       minPt = cms.untracked.double(10.),
    ),
    exotica = cms.PSet(
       maxEta = cms.untracked.double(2.4),
       minIsGlobal = cms.untracked.int32(1),
       minIsTracker = cms.untracked.int32(1),
       maxDxy = cms.untracked.double(0.2),
#       maxNormalizedChi2 = cms.untracked.double(10.), not needed for TeV muons
#       minNTrackerHits = cms.untracked.int32(11), deprecated
       minNPixelHits = cms.untracked.int32(1),
       minNMuonHits = cms.untracked.int32(1),
       minNMatches = cms.untracked.int32(2),
       minNTrackerLayers = cms.untracked.int32(9),
#       minTrackerValidFrac = cms.untracked.double(),
    ),
    )
###WprimeWZ
muonSelectors.WZRelaxed = muonSelectors.VBTF.clone(
    )
muonSelectors.WZRelaxedPt20 = muonSelectors.WZRelaxed.clone(
    minPt = 20.
    )
muonSelectors.WZLoose = muonSelectors.WZRelaxed.clone(
    maxIso03 = cms.untracked.double(0.15)
    )
muonSelectors.WZTight = muonSelectors.WZRelaxed.clone(
    maxIso03 = cms.untracked.double(0.1)
    )
muonSelectors.WZTightPt20 = muonSelectors.WZTight.clone(
    minPt = 20.,
    )
###EWKWZ
muonSelectors.EWKWZRelaxed = muonSelectors.PFIsoRelaxed.clone(
    )
muonSelectors.EWKWZRelaxedPt20 = muonSelectors.EWKWZRelaxed.clone(
    minPt = 20.
    )
muonSelectors.EWKWZLoose = muonSelectors.EWKWZRelaxed.clone(
    maxPFIso = cms.untracked.double(0.2)
    )
muonSelectors.EWKWZLoosePt20 = muonSelectors.EWKWZLoose.clone(
    minPt = 20.,
    )
muonSelectors.EWKWZTight = muonSelectors.EWKWZRelaxed.clone(
    maxPFIso = cms.untracked.double(0.12)
    )
muonSelectors.EWKWZTightPt20 = muonSelectors.EWKWZTight.clone(
    minPt = 20.,
    )

###High Pt PfIso Boosted Z
muonSelectors.PFIsoHighPtBoostedZRelaxedPt20 = muonSelectors.PFIsoHighPtBoostedZRelaxed.clone(
    minPt = 20.
    )
muonSelectors.PFIsoHighPtBoostedZLoose = muonSelectors.PFIsoHighPtBoostedZRelaxed.clone(
    maxPFIso = cms.untracked.double(0.2)
    )
muonSelectors.PFIsoHighPtBoostedZLoosePt20 = muonSelectors.PFIsoHighPtBoostedZLoose.clone(
    minPt = 20.,
    )
muonSelectors.PFIsoHighPtBoostedZTight = muonSelectors.PFIsoHighPtBoostedZRelaxed.clone(
    maxPFIso = cms.untracked.double(0.12)
    )
muonSelectors.PFIsoHighPtBoostedZTightPt20 = muonSelectors.PFIsoHighPtBoostedZTight.clone(
    minPt = 20.,
    )

#####High Pt PfIso
muonSelectors.PFIsoHighPtRelaxedPt20 = muonSelectors.PFIsoHighPtRelaxed.clone(
    minPt = 20.
    )
muonSelectors.PFIsoHighPtLoose = muonSelectors.PFIsoHighPtRelaxed.clone(
    maxPFIso = cms.untracked.double(0.2)
    )
muonSelectors.PFIsoHighPtLoosePt20 = muonSelectors.PFIsoHighPtLoose.clone(
    minPt = 20.,
    )
muonSelectors.PFIsoHighPtTight = muonSelectors.PFIsoHighPtRelaxed.clone(
    maxPFIso = cms.untracked.double(0.12)
    )
muonSelectors.PFIsoHighPtTightPt20 = muonSelectors.PFIsoHighPtTight.clone(
    minPt = 20.,
    )


###HadVZ
muonSelectors.HadVZLoose = muonSelectors.exotica.clone(
    minPt = cms.untracked.double(20.),
#    maxIso03 = cms.untracked.double(0.1)
    )
muonSelectors.HadVZTight = muonSelectors.exotica.clone(
    minPt = cms.untracked.double(35.),
#    maxIso03 = cms.untracked.double(0.1)
    )

#####################
####  Electrons  ####
#####################

cutsMissingHits = [0, 0, 0, 0, 0, 0]
cutsConvDist = [0., 0., 0.02, 0.02, 0.02, 0.02]
cutsConvDcot = [0., 0., 0.02, 0.02, 0.02, 0.02]
cutsEBCombRelIso = [0.150, 0.085, 0.053, 0.070, 0.030, 0.016]#2010
#cutsEBCombRelIso = [0.150, 0.085, 0.053, 0.040, 0.030, 0.016]#2011
cutsEBSigmaIEtaIEta = [0.012, 0.01, 0.01, 0.01, 0.01, 0.01]
cutsEBDeltaPhi = [0.800, 0.071, 0.039, 0.027, 0.020, 0.020]
cutsEBDeltaEta = [0.007, 0.007, 0.005, 0.005, 0.004, 0.004]
#cutsEECombRelIso = [0.100, 0.051, 0.042, 0.033, 0.016, 0.008]#2011
cutsEECombRelIso = [0.100, 0.051, 0.042, 0.060, 0.016, 0.008]#2010
cutsEESigmaIEtaIEta = [0.031, 0.031, 0.031, 0.031, 0.031, 0.031]
cutsEEDeltaPhi = [0.7, 0.047, 0.028, 0.021, 0.021, 0.021]
cutsEEDeltaEta = [0.011, 0.011, 0.007, 0.006, 0.005, 0.004]

electronSelectors = cms.PSet(
    WZLoose = cms.PSet(),
    WZRelaxed95 = cms.PSet(),
    WZRelaxed80 = cms.PSet(),
    WZTight = cms.PSet(),
    exotica = cms.PSet(
       joint = cms.PSet(
          minPt = cms.untracked.double(10.),
          minIsEcalDriven  = cms.untracked.int32(1),
          maxDeltaPhi = cms.untracked.double(0.06),
          maxHoverE = cms.untracked.double(0.05),
          #maxTrackIso = cms.untracked.double(5.),
          maxPFIso = cms.untracked.double(0.15),          
          maxMissingHits = cms.untracked.int32(1),
       ),
       barrel = cms.PSet(
          maxDeltaEta = cms.untracked.double(0.005),
          #maxSigmaIEtaIEta = cms.untracked.double(),
          minPassEX5overE55 = cms.untracked.int32(1),
          #minPassEMHadDepth1Iso = cms.untracked.int32(1),
          maxd0 = cms.untracked.double(0.02),
       ),
       endcap = cms.PSet(
          maxDeltaEta = cms.untracked.double(0.007),
          maxSigmaIEtaIEta = cms.untracked.double(0.03),
          minPassEX5overE55 = cms.untracked.int32(1),
          #minPassEMHadDepth1Iso = cms.untracked.int32(1),
          maxd0 = cms.untracked.double(0.05),
       ),
    ),
    mva = cms.PSet(
       joint = cms.PSet(
          minPt = cms.untracked.double(10.),
          minpassMVAPresel = cms.untracked.int32(1),
          minpassMVATrig = cms.untracked.int32(1),
          #maxPFIso = cms.untracked.double(???),#Added below in clones
       ),
       barrel = cms.PSet(
       ),
       endcap = cms.PSet(
       ),
    ),
    EWKWZLoose = cms.PSet(),
    EWKWZTight = cms.PSet(),

    CiC2012VetoRelaxed = cms.PSet(
       joint = cms.PSet(
          minPt = cms.untracked.double(10.),
          maxIsGap = cms.untracked.int32(0),
          minPassConvVeto = cms.untracked.int32(1),
       ),
       barrel = cms.PSet(
           #maxMissingHits = cms.untracked.int32(1),
           #minVtxFitProb = cms.untracked.double(1e-6),
           maxSigmaIEtaIEta = cms.untracked.double(0.01),
           maxDeltaPhi = cms.untracked.double(0.8),
           maxDeltaEta = cms.untracked.double(0.007),
           maxHoverE = cms.untracked.double(0.15),
           maxd0 = cms.untracked.double(0.04),
           maxdz = cms.untracked.double(0.2),
           #maxfabsdiffEp = cms.untracked.double(0.05),
       ),
       endcap = cms.PSet(
           #maxMissingHits = cms.untracked.int32(1),
           #minVtxFitProb = cms.untracked.double(1e-6),
           maxSigmaIEtaIEta = cms.untracked.double(0.03),
           maxDeltaPhi = cms.untracked.double(0.7),
           maxDeltaEta = cms.untracked.double(0.01),
           #maxHoverE = cms.untracked.double(0.10),
           maxd0 = cms.untracked.double(0.04),
           maxdz = cms.untracked.double(0.2),
           #maxfabsdiffEp = cms.untracked.double(0.05),
       ),
    ),
    CiC2012LooseRelaxed = cms.PSet(
       joint = cms.PSet(
          minPt = cms.untracked.double(10.),
          maxIsGap = cms.untracked.int32(0),
          minPassConvVeto = cms.untracked.int32(1),
       ),
       barrel = cms.PSet(
           maxMissingHits = cms.untracked.int32(1),
           minVtxFitProb = cms.untracked.double(1e-6),
           maxSigmaIEtaIEta = cms.untracked.double(0.01),
           maxDeltaPhi = cms.untracked.double(0.15),
           maxDeltaEta = cms.untracked.double(0.007),
           maxHoverE = cms.untracked.double(0.12),
           maxd0 = cms.untracked.double(0.02),
           maxdz = cms.untracked.double(0.2),
           maxfabsdiffEp = cms.untracked.double(0.05),
       ),
       endcap = cms.PSet(
           maxMissingHits = cms.untracked.int32(1),
           minVtxFitProb = cms.untracked.double(1e-6),
           maxSigmaIEtaIEta = cms.untracked.double(0.03),
           maxDeltaPhi = cms.untracked.double(0.10),
           maxDeltaEta = cms.untracked.double(0.009),
           maxHoverE = cms.untracked.double(0.10),
           maxd0 = cms.untracked.double(0.02),
           maxdz = cms.untracked.double(0.2),
           maxfabsdiffEp = cms.untracked.double(0.05),
       ),
    ),
    Reco = cms.PSet(
       joint = cms.PSet(
          minPt = cms.untracked.double(10.),
       ),
       barrel = cms.PSet(),
       endcap = cms.PSet(),
    ),
    CiC2012Loose = cms.PSet(),
    CiC2012LoosePt20 = cms.PSet(),
    CiC2012MediumRelaxed = cms.PSet(
       joint = cms.PSet(
          minPt = cms.untracked.double(10.),
          maxIsGap = cms.untracked.int32(0),
          minPassConvVeto = cms.untracked.int32(1),
       ),
       barrel = cms.PSet(
           maxMissingHits = cms.untracked.int32(1),
           minVtxFitProb = cms.untracked.double(1e-6),
           maxSigmaIEtaIEta = cms.untracked.double(0.01),
           maxDeltaPhi = cms.untracked.double(0.06),
           maxDeltaEta = cms.untracked.double(0.004),
           maxHoverE = cms.untracked.double(0.12),
           maxd0 = cms.untracked.double(0.02),
           maxdz = cms.untracked.double(0.1),
           maxfabsdiffEp = cms.untracked.double(0.05),
       ),
       endcap = cms.PSet(
           maxMissingHits = cms.untracked.int32(1),
           minVtxFitProb = cms.untracked.double(1e-6),
           maxSigmaIEtaIEta = cms.untracked.double(0.03),
           maxDeltaPhi = cms.untracked.double(0.03),
           maxDeltaEta = cms.untracked.double(0.007),
           maxHoverE = cms.untracked.double(0.10),
           maxd0 = cms.untracked.double(0.02),
           maxdz = cms.untracked.double(0.1),
           maxfabsdiffEp = cms.untracked.double(0.05),
       ),
    ),
    CiC2012Medium = cms.PSet(),
    CiC2012MediumPt20 = cms.PSet(),
    )

for i, s in enumerate(["wp95", "wp90", "wp85", "wp80", "wp70", "wp60"]):
    pset = cms.PSet(
        joint = cms.PSet(
           minPt = cms.untracked.double(10.),
        ),
        barrel = cms.PSet(
           maxMissingHits = cms.untracked.int32(cutsMissingHits[i]),
           minConv = cms.untracked.double(cutsConvDist[i]),#Hack bc we need an OR of these two cuts below
           maxCombRelIso = cms.untracked.double(cutsEBCombRelIso[i]),
           maxSigmaIEtaIEta = cms.untracked.double(cutsEBSigmaIEtaIEta[i]),
           maxDeltaPhi = cms.untracked.double(cutsEBDeltaPhi[i]),
           maxDeltaEta = cms.untracked.double(cutsEBDeltaEta[i]),
           ),
        endcap = cms.PSet(
           maxMissingHits = cms.untracked.int32(cutsMissingHits[i]),
           minConv = cms.untracked.double(cutsConvDist[i]),#Hack bc we need an or of these two cuts below
           maxCombRelIso = cms.untracked.double(cutsEECombRelIso[i]),
           maxSigmaIEtaIEta = cms.untracked.double(cutsEESigmaIEtaIEta[i]),
           maxDeltaPhi = cms.untracked.double(cutsEEDeltaPhi[i]),
           maxDeltaEta = cms.untracked.double(cutsEEDeltaEta[i]),
           ),
        )
    setattr(electronSelectors, s, pset)
    
electronSelectors.WZLoose = electronSelectors.wp95.clone()

electronSelectors.WZRelaxed95 = electronSelectors.wp95.clone()
electronSelectors.WZRelaxed95.barrel.maxCombRelIso = cms.untracked.double(999999999.)
electronSelectors.WZRelaxed95.endcap.maxCombRelIso = cms.untracked.double(999999999.)

electronSelectors.WZRelaxed80 = electronSelectors.wp80.clone()
electronSelectors.WZRelaxed80.barrel.maxCombRelIso = cms.untracked.double(999999999.)
electronSelectors.WZRelaxed80.endcap.maxCombRelIso = cms.untracked.double(999999999.)

electronSelectors.WZRelaxed80Pt20 = electronSelectors.wp80.clone()
electronSelectors.WZRelaxed80Pt20.joint.minPt = cms.untracked.double(20.)
electronSelectors.WZRelaxed80Pt20.barrel.maxCombRelIso = cms.untracked.double(999999999.)
electronSelectors.WZRelaxed80Pt20.endcap.maxCombRelIso = cms.untracked.double(999999999.)

electronSelectors.WZTight = electronSelectors.wp80.clone()

electronSelectors.WZTightPt20 = electronSelectors.wp80.clone()
electronSelectors.WZTightPt20.joint.minPt = cms.untracked.double(20.)

##EWKWZ
electronSelectors.EWKWZLoose = electronSelectors.mva.clone()
electronSelectors.EWKWZLoose.joint.maxPFIso = cms.untracked.double(0.1)

electronSelectors.EWKWZRelaxed = electronSelectors.mva.clone()

electronSelectors.EWKWZRelaxedPt20 = electronSelectors.mva.clone()
electronSelectors.EWKWZRelaxedPt20.joint.minPt = 20.

electronSelectors.EWKWZTight = electronSelectors.mva.clone()
electronSelectors.EWKWZTight.joint.maxPFIso = cms.untracked.double(0.1)

electronSelectors.EWKWZTightPt20 = electronSelectors.mva.clone()
electronSelectors.EWKWZTightPt20.joint.minPt = 20.
electronSelectors.EWKWZTightPt20.joint.maxPFIso = cms.untracked.double(0.1)

#CiC2012

electronSelectors.CiC2012Veto = electronSelectors.CiC2012VetoRelaxed.clone()
electronSelectors.CiC2012Veto.joint.maxPFIso = cms.untracked.double(0.15)

electronSelectors.CiC2012VetoPt20 = electronSelectors.CiC2012Veto.clone()
electronSelectors.CiC2012VetoPt20.joint.minPt = 20.

electronSelectors.CiC2012Loose = electronSelectors.CiC2012LooseRelaxed.clone()
electronSelectors.CiC2012Loose.joint.maxPFIso = cms.untracked.double(0.15)

electronSelectors.CiC2012LoosePt20 = electronSelectors.CiC2012Loose.clone()
electronSelectors.CiC2012LoosePt20.joint.minPt = 20.
    
electronSelectors.CiC2012Medium = electronSelectors.CiC2012MediumRelaxed.clone()
electronSelectors.CiC2012Medium.joint.maxPFIso = cms.untracked.double(0.15)

electronSelectors.CiC2012MediumPt20 = electronSelectors.CiC2012Medium.clone()
electronSelectors.CiC2012MediumPt20.joint.minPt = 20.

##Exotica
electronSelectors.exoticaPt20 = electronSelectors.exotica.clone()
electronSelectors.exoticaPt20.joint.minPt = cms.untracked.double(20.)

####################
#####  Jets  #######
####################
jetSelectors = cms.PSet(
    Base = cms.PSet(
       minPt = cms.untracked.double(30.),
       maxEta = cms.untracked.double(2.4),
       maxNHF = cms.untracked.double(0.99),
       maxNEF = cms.untracked.double(0.99),
       minNDaughters = cms.untracked.int32(1),
       minCHF = cms.untracked.double(0.0),
       maxCEF = cms.untracked.double(0.99),
       minCMult = cms.untracked.int32(0),
       ),
    Pat = cms.PSet(
       minPt = cms.untracked.double(30.),
       maxEta = cms.untracked.double(3.0),
       ),
    )

####################
#### Photons #######
####################
photonSelectors = cms.PSet(
    Base = cms.PSet(#Made up values
       minPt = cms.untracked.double(10.),
       maxEta = cms.untracked.double(5),
       maxECalIso = cms.untracked.double(0.1),
       maxHCalIso = cms.untracked.double(0.1),
       maxTrkIso = cms.untracked.double(0.1),
       maxHoE = cms.untracked.double(0.05),
       maxSigmaee = cms.untracked.double(9e9),
       minHasSeed = cms.untracked.bool(True),       
       ),
    )

#CiC2012Relaxed
Revision:	1.28
Committed:	Fri Mar 29 20:12:51 2013 UTC (12 years, 1 month ago) by fantasia
Content type:	text/x-python
Branch:	MAIN
Changes since 1.27:	+35 -0 lines
Log Message:	added veto id
#	User	Rev	Content
1	fantasia	1.1	import FWCore.ParameterSet.Config as cms
2	fantasia	1.5
3			#####################
4			###### Muons ######
5			#####################
6
7	fantasia	1.1	#Note: These cuts are <= and >=
8			muonSelectors = cms.PSet(
9	fantasia	1.5	VBTF = cms.PSet(
10	fantasia	1.1	minPt = cms.untracked.double(10.),
11			maxEta = cms.untracked.double(2.4),
12			minIsGlobal = cms.untracked.int32(1),
13			maxDxy = cms.untracked.double(0.2),
14			maxNormalizedChi2 = cms.untracked.double(10.0),
15			minNTrackerHits = cms.untracked.int32(11),
16			minNPixelHits = cms.untracked.int32(1),
17			minNMuonHits = cms.untracked.int32(1),
18			minNMatches = cms.untracked.int32(2),
19			),
20	fantasia	1.26	PFIsoHighPtBoostedZRelaxed = cms.PSet(
21	fantasia	1.17	minPt = cms.untracked.double(10.),
22	fantasia	1.22	maxEta = cms.untracked.double(2.4),
23	fantasia	1.20	minIsGblOrTrk = cms.untracked.int32(1),
24	fantasia	1.24	maxDxy = cms.untracked.double(0.2),
25			maxDz = cms.untracked.double(0.5),
26			maxRelPtErr = cms.untracked.double(0.3),
27			minNPixelHits = cms.untracked.int32(1),
28			minNTrackerLayers = cms.untracked.int32(6),
29			minNMatches = cms.untracked.int32(2),
30			#maxPFIso = cms.untracked.double(???),#Added below in clones
31			),
32			PFIsoHighPtRelaxed = cms.PSet(
33			minPt = cms.untracked.double(10.),
34			maxEta = cms.untracked.double(2.4),
35			minIsGlobal = cms.untracked.int32(1),
36	fantasia	1.17	maxDxy = cms.untracked.double(0.2),
37	fantasia	1.19	maxDz = cms.untracked.double(0.5),
38	fantasia	1.26	maxRelPtErr = cms.untracked.double(0.3),
39			minNPixelHits = cms.untracked.int32(1),
40			minNMuonHits = cms.untracked.int32(1),
41			minNTrackerLayers = cms.untracked.int32(6),
42			minNMatches = cms.untracked.int32(2),
43			#maxPFIso = cms.untracked.double(???),#Added below in clones
44			),
45			PFIsoRelaxed = cms.PSet(
46			minPt = cms.untracked.double(10.),
47			maxEta = cms.untracked.double(2.4),
48			minIsGlobal = cms.untracked.int32(1),
49	fantasia	1.27	minIsPF = cms.untracked.int32(1),
50	fantasia	1.26	maxDxy = cms.untracked.double(0.2),
51			maxDz = cms.untracked.double(0.5),
52			maxNormalizedChi2 = cms.untracked.double(10.0),
53	fantasia	1.21	maxRelPtErr = cms.untracked.double(0.3),
54	fantasia	1.17	minNPixelHits = cms.untracked.int32(1),
55	fantasia	1.24	minNMuonHits = cms.untracked.int32(1),
56	fantasia	1.19	minNTrackerLayers = cms.untracked.int32(6),
57	fantasia	1.17	minNMatches = cms.untracked.int32(2),
58			#maxPFIso = cms.untracked.double(???),#Added below in clones
59			),
60	clint	1.25	Reco = cms.PSet(
61			minPt = cms.untracked.double(10.),
62			),
63	fantasia	1.1	exotica = cms.PSet(
64	fantasia	1.9	maxEta = cms.untracked.double(2.4),
65			minIsGlobal = cms.untracked.int32(1),
66			minIsTracker = cms.untracked.int32(1),
67			maxDxy = cms.untracked.double(0.2),
68	cleonido	1.10	# maxNormalizedChi2 = cms.untracked.double(10.), not needed for TeV muons
69			# minNTrackerHits = cms.untracked.int32(11), deprecated
70	fantasia	1.9	minNPixelHits = cms.untracked.int32(1),
71			minNMuonHits = cms.untracked.int32(1),
72			minNMatches = cms.untracked.int32(2),
73	cleonido	1.10	minNTrackerLayers = cms.untracked.int32(9),
74	fantasia	1.9	# minTrackerValidFrac = cms.untracked.double(),
75	fantasia	1.1	),
76			)
77	fantasia	1.17	###WprimeWZ
78	fantasia	1.9	muonSelectors.WZRelaxed = muonSelectors.VBTF.clone(
79	fantasia	1.18	)
80	fantasia	1.19	muonSelectors.WZRelaxedPt20 = muonSelectors.WZRelaxed.clone(
81	fantasia	1.9	minPt = 20.
82			)
83	fantasia	1.19	muonSelectors.WZLoose = muonSelectors.WZRelaxed.clone(
84			maxIso03 = cms.untracked.double(0.15)
85			)
86			muonSelectors.WZTight = muonSelectors.WZRelaxed.clone(
87	fantasia	1.18	maxIso03 = cms.untracked.double(0.1)
88			)
89	fantasia	1.19	muonSelectors.WZTightPt20 = muonSelectors.WZTight.clone(
90	fantasia	1.9	minPt = 20.,
91			)
92	fantasia	1.17	###EWKWZ
93	fantasia	1.26	muonSelectors.EWKWZRelaxed = muonSelectors.PFIsoRelaxed.clone(
94	fantasia	1.18	)
95	fantasia	1.19	muonSelectors.EWKWZRelaxedPt20 = muonSelectors.EWKWZRelaxed.clone(
96	fantasia	1.17	minPt = 20.
97			)
98	fantasia	1.19	muonSelectors.EWKWZLoose = muonSelectors.EWKWZRelaxed.clone(
99			maxPFIso = cms.untracked.double(0.2)
100	fantasia	1.20	)
101			muonSelectors.EWKWZLoosePt20 = muonSelectors.EWKWZLoose.clone(
102			minPt = 20.,
103	fantasia	1.19	)
104			muonSelectors.EWKWZTight = muonSelectors.EWKWZRelaxed.clone(
105	fantasia	1.18	maxPFIso = cms.untracked.double(0.12)
106			)
107	fantasia	1.19	muonSelectors.EWKWZTightPt20 = muonSelectors.EWKWZTight.clone(
108	fantasia	1.17	minPt = 20.,
109			)
110	fantasia	1.18
111	fantasia	1.26	###High Pt PfIso Boosted Z
112			muonSelectors.PFIsoHighPtBoostedZRelaxedPt20 = muonSelectors.PFIsoHighPtBoostedZRelaxed.clone(
113			minPt = 20.
114			)
115			muonSelectors.PFIsoHighPtBoostedZLoose = muonSelectors.PFIsoHighPtBoostedZRelaxed.clone(
116			maxPFIso = cms.untracked.double(0.2)
117			)
118			muonSelectors.PFIsoHighPtBoostedZLoosePt20 = muonSelectors.PFIsoHighPtBoostedZLoose.clone(
119			minPt = 20.,
120			)
121			muonSelectors.PFIsoHighPtBoostedZTight = muonSelectors.PFIsoHighPtBoostedZRelaxed.clone(
122			maxPFIso = cms.untracked.double(0.12)
123			)
124			muonSelectors.PFIsoHighPtBoostedZTightPt20 = muonSelectors.PFIsoHighPtBoostedZTight.clone(
125			minPt = 20.,
126			)
127
128	fantasia	1.24	#####High Pt PfIso
129			muonSelectors.PFIsoHighPtRelaxedPt20 = muonSelectors.PFIsoHighPtRelaxed.clone(
130			minPt = 20.
131			)
132			muonSelectors.PFIsoHighPtLoose = muonSelectors.PFIsoHighPtRelaxed.clone(
133			maxPFIso = cms.untracked.double(0.2)
134			)
135			muonSelectors.PFIsoHighPtLoosePt20 = muonSelectors.PFIsoHighPtLoose.clone(
136			minPt = 20.,
137			)
138			muonSelectors.PFIsoHighPtTight = muonSelectors.PFIsoHighPtRelaxed.clone(
139			maxPFIso = cms.untracked.double(0.12)
140			)
141			muonSelectors.PFIsoHighPtTightPt20 = muonSelectors.PFIsoHighPtTight.clone(
142			minPt = 20.,
143			)
144
145
146	fantasia	1.17	###HadVZ
147	fantasia	1.11	muonSelectors.HadVZLoose = muonSelectors.exotica.clone(
148	fantasia	1.12	minPt = cms.untracked.double(20.),
149	fladias	1.16	# maxIso03 = cms.untracked.double(0.1)
150	fantasia	1.9	)
151	fantasia	1.11	muonSelectors.HadVZTight = muonSelectors.exotica.clone(
152	fantasia	1.12	minPt = cms.untracked.double(35.),
153	fladias	1.16	# maxIso03 = cms.untracked.double(0.1)
154	fantasia	1.9	)
155	fantasia	1.5
156			#####################
157			#### Electrons ####
158			#####################
159
160	fantasia	1.1	cutsMissingHits = [0, 0, 0, 0, 0, 0]
161			cutsConvDist = [0., 0., 0.02, 0.02, 0.02, 0.02]
162			cutsConvDcot = [0., 0., 0.02, 0.02, 0.02, 0.02]
163	fantasia	1.2	cutsEBCombRelIso = [0.150, 0.085, 0.053, 0.070, 0.030, 0.016]#2010
164			#cutsEBCombRelIso = [0.150, 0.085, 0.053, 0.040, 0.030, 0.016]#2011
165	fantasia	1.1	cutsEBSigmaIEtaIEta = [0.012, 0.01, 0.01, 0.01, 0.01, 0.01]
166			cutsEBDeltaPhi = [0.800, 0.071, 0.039, 0.027, 0.020, 0.020]
167			cutsEBDeltaEta = [0.007, 0.007, 0.005, 0.005, 0.004, 0.004]
168			#cutsEECombRelIso = [0.100, 0.051, 0.042, 0.033, 0.016, 0.008]#2011
169			cutsEECombRelIso = [0.100, 0.051, 0.042, 0.060, 0.016, 0.008]#2010
170			cutsEESigmaIEtaIEta = [0.031, 0.031, 0.031, 0.031, 0.031, 0.031]
171			cutsEEDeltaPhi = [0.7, 0.047, 0.028, 0.021, 0.021, 0.021]
172			cutsEEDeltaEta = [0.011, 0.011, 0.007, 0.006, 0.005, 0.004]
173
174			electronSelectors = cms.PSet(
175			WZLoose = cms.PSet(),
176	fantasia	1.18	WZRelaxed95 = cms.PSet(),
177			WZRelaxed80 = cms.PSet(),
178	fantasia	1.1	WZTight = cms.PSet(),
179	fantasia	1.19	exotica = cms.PSet(
180			joint = cms.PSet(
181			minPt = cms.untracked.double(10.),
182			minIsEcalDriven = cms.untracked.int32(1),
183			maxDeltaPhi = cms.untracked.double(0.06),
184			maxHoverE = cms.untracked.double(0.05),
185			#maxTrackIso = cms.untracked.double(5.),
186			maxPFIso = cms.untracked.double(0.15),
187			maxMissingHits = cms.untracked.int32(1),
188			),
189			barrel = cms.PSet(
190			maxDeltaEta = cms.untracked.double(0.005),
191			#maxSigmaIEtaIEta = cms.untracked.double(),
192			minPassEX5overE55 = cms.untracked.int32(1),
193			#minPassEMHadDepth1Iso = cms.untracked.int32(1),
194			maxd0 = cms.untracked.double(0.02),
195			),
196			endcap = cms.PSet(
197			maxDeltaEta = cms.untracked.double(0.007),
198			maxSigmaIEtaIEta = cms.untracked.double(0.03),
199			minPassEX5overE55 = cms.untracked.int32(1),
200			#minPassEMHadDepth1Iso = cms.untracked.int32(1),
201			maxd0 = cms.untracked.double(0.05),
202			),
203			),
204	fantasia	1.18	mva = cms.PSet(
205			joint = cms.PSet(
206			minPt = cms.untracked.double(10.),
207			minpassMVAPresel = cms.untracked.int32(1),
208			minpassMVATrig = cms.untracked.int32(1),
209			#maxPFIso = cms.untracked.double(???),#Added below in clones
210			),
211			barrel = cms.PSet(
212			),
213			endcap = cms.PSet(
214			),
215			),
216			EWKWZLoose = cms.PSet(),
217			EWKWZTight = cms.PSet(),
218	fantasia	1.19
219	fantasia	1.28	CiC2012VetoRelaxed = cms.PSet(
220			joint = cms.PSet(
221			minPt = cms.untracked.double(10.),
222			maxIsGap = cms.untracked.int32(0),
223			minPassConvVeto = cms.untracked.int32(1),
224			),
225			barrel = cms.PSet(
226			#maxMissingHits = cms.untracked.int32(1),
227			#minVtxFitProb = cms.untracked.double(1e-6),
228			maxSigmaIEtaIEta = cms.untracked.double(0.01),
229			maxDeltaPhi = cms.untracked.double(0.8),
230			maxDeltaEta = cms.untracked.double(0.007),
231			maxHoverE = cms.untracked.double(0.15),
232			maxd0 = cms.untracked.double(0.04),
233			maxdz = cms.untracked.double(0.2),
234			#maxfabsdiffEp = cms.untracked.double(0.05),
235			),
236			endcap = cms.PSet(
237			#maxMissingHits = cms.untracked.int32(1),
238			#minVtxFitProb = cms.untracked.double(1e-6),
239			maxSigmaIEtaIEta = cms.untracked.double(0.03),
240			maxDeltaPhi = cms.untracked.double(0.7),
241			maxDeltaEta = cms.untracked.double(0.01),
242			#maxHoverE = cms.untracked.double(0.10),
243			maxd0 = cms.untracked.double(0.04),
244			maxdz = cms.untracked.double(0.2),
245			#maxfabsdiffEp = cms.untracked.double(0.05),
246			),
247			),
248	fantasia	1.23	CiC2012LooseRelaxed = cms.PSet(
249	fantasia	1.22	joint = cms.PSet(
250			minPt = cms.untracked.double(10.),
251			maxIsGap = cms.untracked.int32(0),
252	fantasia	1.24	minPassConvVeto = cms.untracked.int32(1),
253	fantasia	1.22	),
254			barrel = cms.PSet(
255			maxMissingHits = cms.untracked.int32(1),
256			minVtxFitProb = cms.untracked.double(1e-6),
257			maxSigmaIEtaIEta = cms.untracked.double(0.01),
258			maxDeltaPhi = cms.untracked.double(0.15),
259			maxDeltaEta = cms.untracked.double(0.007),
260			maxHoverE = cms.untracked.double(0.12),
261			maxd0 = cms.untracked.double(0.02),
262			maxdz = cms.untracked.double(0.2),
263			maxfabsdiffEp = cms.untracked.double(0.05),
264			),
265			endcap = cms.PSet(
266			maxMissingHits = cms.untracked.int32(1),
267			minVtxFitProb = cms.untracked.double(1e-6),
268			maxSigmaIEtaIEta = cms.untracked.double(0.03),
269			maxDeltaPhi = cms.untracked.double(0.10),
270			maxDeltaEta = cms.untracked.double(0.009),
271			maxHoverE = cms.untracked.double(0.10),
272			maxd0 = cms.untracked.double(0.02),
273			maxdz = cms.untracked.double(0.2),
274			maxfabsdiffEp = cms.untracked.double(0.05),
275			),
276			),
277	clint	1.25	Reco = cms.PSet(
278			joint = cms.PSet(
279			minPt = cms.untracked.double(10.),
280			),
281			barrel = cms.PSet(),
282			endcap = cms.PSet(),
283			),
284	fantasia	1.23	CiC2012Loose = cms.PSet(),
285	fantasia	1.21	CiC2012LoosePt20 = cms.PSet(),
286	fantasia	1.23	CiC2012MediumRelaxed = cms.PSet(
287	fantasia	1.19	joint = cms.PSet(
288	fantasia	1.21	minPt = cms.untracked.double(10.),
289	fantasia	1.22	maxIsGap = cms.untracked.int32(0),
290	fantasia	1.24	minPassConvVeto = cms.untracked.int32(1),
291	fantasia	1.19	),
292			barrel = cms.PSet(
293			maxMissingHits = cms.untracked.int32(1),
294			minVtxFitProb = cms.untracked.double(1e-6),
295			maxSigmaIEtaIEta = cms.untracked.double(0.01),
296			maxDeltaPhi = cms.untracked.double(0.06),
297			maxDeltaEta = cms.untracked.double(0.004),
298			maxHoverE = cms.untracked.double(0.12),
299			maxd0 = cms.untracked.double(0.02),
300	fantasia	1.22	maxdz = cms.untracked.double(0.1),
301	fantasia	1.19	maxfabsdiffEp = cms.untracked.double(0.05),
302			),
303			endcap = cms.PSet(
304			maxMissingHits = cms.untracked.int32(1),
305			minVtxFitProb = cms.untracked.double(1e-6),
306			maxSigmaIEtaIEta = cms.untracked.double(0.03),
307			maxDeltaPhi = cms.untracked.double(0.03),
308			maxDeltaEta = cms.untracked.double(0.007),
309			maxHoverE = cms.untracked.double(0.10),
310			maxd0 = cms.untracked.double(0.02),
311	fantasia	1.22	maxdz = cms.untracked.double(0.1),
312	fantasia	1.19	maxfabsdiffEp = cms.untracked.double(0.05),
313			),
314			),
315	fantasia	1.23	CiC2012Medium = cms.PSet(),
316	fantasia	1.21	CiC2012MediumPt20 = cms.PSet(),
317	fantasia	1.1	)
318
319			for i, s in enumerate(["wp95", "wp90", "wp85", "wp80", "wp70", "wp60"]):
320			pset = cms.PSet(
321	fantasia	1.18	joint = cms.PSet(
322			minPt = cms.untracked.double(10.),
323			),
324	fantasia	1.1	barrel = cms.PSet(
325			maxMissingHits = cms.untracked.int32(cutsMissingHits[i]),
326			minConv = cms.untracked.double(cutsConvDist[i]),#Hack bc we need an OR of these two cuts below
327			maxCombRelIso = cms.untracked.double(cutsEBCombRelIso[i]),
328			maxSigmaIEtaIEta = cms.untracked.double(cutsEBSigmaIEtaIEta[i]),
329			maxDeltaPhi = cms.untracked.double(cutsEBDeltaPhi[i]),
330			maxDeltaEta = cms.untracked.double(cutsEBDeltaEta[i]),
331			),
332			endcap = cms.PSet(
333			maxMissingHits = cms.untracked.int32(cutsMissingHits[i]),
334			minConv = cms.untracked.double(cutsConvDist[i]),#Hack bc we need an or of these two cuts below
335			maxCombRelIso = cms.untracked.double(cutsEECombRelIso[i]),
336			maxSigmaIEtaIEta = cms.untracked.double(cutsEESigmaIEtaIEta[i]),
337			maxDeltaPhi = cms.untracked.double(cutsEEDeltaPhi[i]),
338			maxDeltaEta = cms.untracked.double(cutsEEDeltaEta[i]),
339			),
340			)
341			setattr(electronSelectors, s, pset)
342
343			electronSelectors.WZLoose = electronSelectors.wp95.clone()
344	fantasia	1.9
345	fantasia	1.18	electronSelectors.WZRelaxed95 = electronSelectors.wp95.clone()
346			electronSelectors.WZRelaxed95.barrel.maxCombRelIso = cms.untracked.double(999999999.)
347			electronSelectors.WZRelaxed95.endcap.maxCombRelIso = cms.untracked.double(999999999.)
348
349			electronSelectors.WZRelaxed80 = electronSelectors.wp80.clone()
350			electronSelectors.WZRelaxed80.barrel.maxCombRelIso = cms.untracked.double(999999999.)
351			electronSelectors.WZRelaxed80.endcap.maxCombRelIso = cms.untracked.double(999999999.)
352
353			electronSelectors.WZRelaxed80Pt20 = electronSelectors.wp80.clone()
354			electronSelectors.WZRelaxed80Pt20.joint.minPt = cms.untracked.double(20.)
355			electronSelectors.WZRelaxed80Pt20.barrel.maxCombRelIso = cms.untracked.double(999999999.)
356			electronSelectors.WZRelaxed80Pt20.endcap.maxCombRelIso = cms.untracked.double(999999999.)
357	fantasia	1.1
358			electronSelectors.WZTight = electronSelectors.wp80.clone()
359
360	fantasia	1.18	electronSelectors.WZTightPt20 = electronSelectors.wp80.clone()
361			electronSelectors.WZTightPt20.joint.minPt = cms.untracked.double(20.)
362
363			##EWKWZ
364			electronSelectors.EWKWZLoose = electronSelectors.mva.clone()
365			electronSelectors.EWKWZLoose.joint.maxPFIso = cms.untracked.double(0.1)
366
367			electronSelectors.EWKWZRelaxed = electronSelectors.mva.clone()
368
369			electronSelectors.EWKWZRelaxedPt20 = electronSelectors.mva.clone()
370			electronSelectors.EWKWZRelaxedPt20.joint.minPt = 20.
371
372			electronSelectors.EWKWZTight = electronSelectors.mva.clone()
373			electronSelectors.EWKWZTight.joint.maxPFIso = cms.untracked.double(0.1)
374
375			electronSelectors.EWKWZTightPt20 = electronSelectors.mva.clone()
376			electronSelectors.EWKWZTightPt20.joint.minPt = 20.
377			electronSelectors.EWKWZTightPt20.joint.maxPFIso = cms.untracked.double(0.1)
378	fantasia	1.9
379	fantasia	1.21	#CiC2012
380	fantasia	1.23
381	fantasia	1.28	electronSelectors.CiC2012Veto = electronSelectors.CiC2012VetoRelaxed.clone()
382			electronSelectors.CiC2012Veto.joint.maxPFIso = cms.untracked.double(0.15)
383
384			electronSelectors.CiC2012VetoPt20 = electronSelectors.CiC2012Veto.clone()
385			electronSelectors.CiC2012VetoPt20.joint.minPt = 20.
386
387	fantasia	1.23	electronSelectors.CiC2012Loose = electronSelectors.CiC2012LooseRelaxed.clone()
388			electronSelectors.CiC2012Loose.joint.maxPFIso = cms.untracked.double(0.15)
389
390	fantasia	1.21	electronSelectors.CiC2012LoosePt20 = electronSelectors.CiC2012Loose.clone()
391			electronSelectors.CiC2012LoosePt20.joint.minPt = 20.
392
393	fantasia	1.23	electronSelectors.CiC2012Medium = electronSelectors.CiC2012MediumRelaxed.clone()
394			electronSelectors.CiC2012Medium.joint.maxPFIso = cms.untracked.double(0.15)
395
396	fantasia	1.21	electronSelectors.CiC2012MediumPt20 = electronSelectors.CiC2012Medium.clone()
397			electronSelectors.CiC2012MediumPt20.joint.minPt = 20.
398	fantasia	1.19
399			##Exotica
400			electronSelectors.exoticaPt20 = electronSelectors.exotica.clone()
401			electronSelectors.exoticaPt20.joint.minPt = cms.untracked.double(20.)
402
403	fantasia	1.4	####################
404			##### Jets #######
405			####################
406			jetSelectors = cms.PSet(
407	fantasia	1.8	Base = cms.PSet(
408	fladias	1.14	minPt = cms.untracked.double(30.),
409	fantasia	1.4	maxEta = cms.untracked.double(2.4),
410			maxNHF = cms.untracked.double(0.99),
411			maxNEF = cms.untracked.double(0.99),
412	fantasia	1.6	minNDaughters = cms.untracked.int32(1),
413			minCHF = cms.untracked.double(0.0),
414	fantasia	1.4	maxCEF = cms.untracked.double(0.99),
415	fantasia	1.6	minCMult = cms.untracked.int32(0),
416	fantasia	1.4	),
417	fantasia	1.17	Pat = cms.PSet(
418			minPt = cms.untracked.double(30.),
419			maxEta = cms.untracked.double(3.0),
420			),
421	fantasia	1.4	)
422	fantasia	1.8
423			####################
424			#### Photons #######
425			####################
426			photonSelectors = cms.PSet(
427			Base = cms.PSet(#Made up values
428			minPt = cms.untracked.double(10.),
429			maxEta = cms.untracked.double(5),
430			maxECalIso = cms.untracked.double(0.1),
431			maxHCalIso = cms.untracked.double(0.1),
432			maxTrkIso = cms.untracked.double(0.1),
433			maxHoE = cms.untracked.double(0.05),
434			maxSigmaee = cms.untracked.double(9e9),
435			minHasSeed = cms.untracked.bool(True),
436			),
437			)
438
439	fantasia	1.23	#CiC2012Relaxed