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