RecoEgamma/python/ElectronMcSignalValidatorPt1000_gedGsfElectrons_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 electronMcSignalHistosCfg = cms.PSet(
0004   Nbinxyz = cms.int32(50),
0005   Nbinp = cms.int32(50), Nbinp2D = cms.int32(50),
0006   Pmax = cms.double(4500),
0007   Nbinpt = cms.int32(50), Nbinpt2D = cms.int32(50), Nbinpteff = cms.int32(19),Ptmax = cms.double(1500.0),
0008   Nbinfhits = cms.int32(30), Fhitsmax = cms.double(30.0),
0009   Nbinlhits = cms.int32(5), Lhitsmax = cms.double(10.0),
0010   Nbineta = cms.int32(50), Nbineta2D = cms.int32(50),Etamin = cms.double(-2.5), Etamax = cms.double(2.5),
0011   NbinetaExtended = cms.int32(60), Nbineta2DExtended = cms.int32(60),EtaminExtended = cms.double(-3.0), EtamaxExtended = cms.double(3.0),
0012   Nbindeta = cms.int32(100), Detamin = cms.double(-0.005), Detamax = cms.double(0.005),
0013   Nbindetamatch = cms.int32(100), Nbindetamatch2D = cms.int32(50), Detamatchmin = cms.double(-0.05), Detamatchmax = cms.double(0.05),
0014   Nbinphi = cms.int32(64), Nbinphi2D = cms.int32(32), Phimin = cms.double(-3.2), Phimax = cms.double(3.2),
0015   Nbindphi = cms.int32(100), Dphimin = cms.double(-0.01), Dphimax = cms.double(0.01),
0016   Nbindphimatch = cms.int32(100), Nbindphimatch2D = cms.int32(50), Dphimatchmin = cms.double(-0.2), Dphimatchmax = cms.double(0.2),
0017   Nbineop = cms.int32(50), Nbineop2D = cms.int32(30),
0018   Eopmax = cms.double(10),
0019   Eopmaxsht = cms.double(3.0),
0020   Nbinmee = cms.int32(100), Meemin = cms.double(0.0), Meemax = cms.double(150.),
0021   Nbinhoe = cms.int32(100), Hoemin = cms.double(0.0), Hoemax = cms.double(0.5),
0022   Nbinpoptrue = cms.int32(75), Poptruemin = cms.double(0.0), Poptruemax = cms.double(1.5),
0023   Nbinerror = cms.int32(30), Energyerrormax = cms.double(150.0),
0024   EfficiencyFlag = cms.bool(True), StatOverflowFlag = cms.bool(False),
0025   NbinOPV = cms.int32(80), OPV_min = cms.double(-0.5), OPV_max = cms.double(79.5), # OPV : Offline Primary Vertices
0026   NbinELE = cms.int32(11), ELE_min = cms.double(-0.5), ELE_max = cms.double(10.5), # ELE : recEleNum
0027   NbinCORE = cms.int32(21), CORE_min = cms.double(-0.5), CORE_max = cms.double(20.5), # CORE : recCoreNum
0028   NbinTRACK = cms.int32(41), TRACK_min = cms.double(-0.5), TRACK_max = cms.double(40.5), # TRACK : recTrackNum
0029   NbinSEED = cms.int32(101), SEED_min = cms.double(-0.5), SEED_max = cms.double(100.5), # SEED : recSeedNum
0030 )
0031
0032 from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
0033 electronMcSignalValidatorPt1000 = DQMEDAnalyzer('ElectronMcSignalValidator',
0034
0035   Verbosity = cms.untracked.int32(0),
0036   FinalStep = cms.string("AtJobEnd"),
0037   InputFile = cms.string(""),
0038   OutputFile = cms.string(""),
0039   InputFolderName = cms.string("EgammaV/ElectronMcSignalValidatorPt1000"),
0040   OutputFolderName = cms.string("EgammaV/ElectronMcSignalValidatorPt1000"),
0041
0042   mcTruthCollection = cms.InputTag("genParticles"),
0043   electronCollection = cms.InputTag("gedGsfElectrons"),
0044   electronCollectionEndcaps = cms.InputTag("gedGsfElectrons"),
0045   electronCoreCollection = cms.InputTag("gedGsfElectronCores"),
0046   electronCoreCollectionEndcaps = cms.InputTag("gedGsfElectronCores"),
0047   electronTrackCollection = cms.InputTag("electronGsfTracks"),
0048   electronSeedCollection = cms.InputTag("electronMergedSeeds"),
0049   offlinePrimaryVertices = cms.InputTag("offlinePrimaryVertices"),
0050
0051   beamSpot = cms.InputTag("offlineBeamSpot"),
0052   readAOD = cms.bool(False),
0053
0054   isoFromDepsTk03            = cms.InputTag(""),
0055   isoFromDepsTk04            = cms.InputTag(""),
0056   isoFromDepsEcalFull03      = cms.InputTag(""),
0057   isoFromDepsEcalFull04      = cms.InputTag(""),
0058   isoFromDepsEcalReduced03   = cms.InputTag(""),
0059   isoFromDepsEcalReduced04   = cms.InputTag(""),
0060   isoFromDepsHcal03          = cms.InputTag(""),
0061   isoFromDepsHcal04          = cms.InputTag(""),
0062
0063   MaxPt = cms.double(1500.0),
0064   DeltaR = cms.double(0.05),
0065   MaxAbsEta = cms.double(2.5),
0066   MaxAbsEtaExtended = cms.double(3.0),
0067   MatchingID = cms.vint32(11,-11),
0068   MatchingMotherID = cms.vint32(23,24,-24,32,990),
0069   histosCfg = cms.PSet(electronMcSignalHistosCfg)
0070 )
0071
0072 from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
0073 phase2_hgcal.toModify(
0074     electronMcSignalValidatorPt1000,
0075     electronCollectionEndcaps = 'ecalDrivenGsfElectronsHGC',
0076     electronCoreCollectionEndcaps = 'ecalDrivenGsfElectronCoresHGC',
0077     MaxAbsEta = 3.0,
0078     histosCfg = dict(
0079         Nbineta = 60 ,
0080         Nbineta2D = 60 ,
0081         Etamin = -3.0 ,
0082         Etamax = 3.0 ,
0083
0084         NbinOPV = 125, OPV_min = -0.5, OPV_max = 249.5,
0085         NbinELE = 100, ELE_min = -0.5, ELE_max = 99.5,
0086         NbinCORE = 100, CORE_min = -0.5, CORE_max = 499.5,
0087         NbinTRACK = 100, TRACK_min = -0.5, TRACK_max = 999.5,
0088         NbinSEED = 100, SEED_min = -0.5, SEED_max = 9999.5,
0089    ),
0090 )