Ejemplo n.º 1
0
        "leadPixSeed>>leadPixSeed(2,-0.5,1.5)",
        "subleadPassEleVeto>>subleadPassEleVeto(2,-0.5,1.5)",
        "subleadBlockChIso>>subleadBlockChIso(120,-10,50)",
        "subleadBlockPhoIso>>subleadBlockPhoIso(120,-10,50)",
        "subleadChIso>>subleadChIso(120,-10,50)",
        "subleadPhoIso>>subleadPhoIso(120,-10,50)",
        "subleadNeutIso>>subleadNeutIso(120,-10,50)",
        "subleadHoE>>subleadHoE(40,0,0.2)",
        "subleadSigmaIeIe>>subleadSigmaIeIe(50,0,5.e-2)",
        "subleadPixSeed>>subleadPixSeed(2,-0.5,1.5)",
        "subleadPassEleVeto>>subleadPassEleVeto(2,-0.5,1.5)",
        "subleadPt:leadPt>>ptSubVsLead(145,100,3000:145,100,3000)",
        "minR9>>minR9(110,0,1.1)", "maxEta>>maxEta(250,0,2.5)"
    ])

minimalDumper = diphotonDumper.clone()
cfgTools.dumpOnly(minimalDumper, [
    "mass",
    "pt",
    "leadPt",
    "leadEta",
    "leadScEta",
    "leadPhi",
    "subleadPt",
    "subleadEta",
    "subleadScEta",
    "subleadPhi",
    "leadBlockPhoIso",
    "subleadBlockPhoIso",
    "leadBlockChIso",
    "subleadBlockChIso",
genDiphotonDumper.dumpGlobalVariables = True
genDiphotonDumper.globalVariables = globalVariables
cfgTools.addCategories(genDiphotonDumper,
                       [("EB","max(abs(leadingPhoton.eta),abs(subLeadingPhoton.eta))<1.4442",0),
                        ("EE","1",0)
                        ],
                       variables=["genMass := mass","pt",
                                  "genLeadPt := leadingPhoton.pt",
                                  "genSubeadPt := subLeadingPhoton.pt",
                                  ],
                       histograms=["genMass>>genmass(1500,0,15000)",
                                   ]
                       )

# restricted set of variables to be dumped for ancillary event selections
minimalDumper = diphotonDumper.clone()
cfgTools.dumpOnly(minimalDumper,
                  ["mass","pt","genMass","satRegressedMass","regressedMass",
                   "leadEnergy","subLeadEnergy",
                   "leadSatRegressedEnergy","subLeadSatRegressedEnergy",
                   "leadRegressedEnergy","subLeadRegressedEnergy",
                   "leadIsSat","subLeadIsSat","leadIsWeird","subLeadIsWeird",
                   "leadPt","leadEta","leadScEta","leadPhi",
                   "subleadPt","subleadEta","subleadScEta","subleadPhi",
                   "leadBlockPhoIso","subleadBlockPhoIso",
                   "leadBlockChIso","subleadBlockChIso",
                   "leadRndConePhoIso","leadRndConeChIso",
                   "subleadRndConePhoIso","subleadRndConeChIso",
                   "leadRndConeChIso0","leadRndConeChIso1","leadRndConeChIso2","leadRndConeChIso3",
                   "leadRndConeChIso4","leadRndConeChIso5","leadRndConeChIso6","leadRndConeChIso7","leadRndConeChIso8",
                   "subleadRndConeChIso0","subleadRndConeChIso1","subleadRndConeChIso2","subleadRndConeChIso3",
Ejemplo n.º 3
0
process.GlobalTag.globaltag = 'POSTLS170_V5::All'
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
process.MessageLogger.cerr.FwkReport.reportEvery = cms.untracked.int32( 1000 )



process.source = cms.Source("PoolSource",
                            fileNames=cms.untracked.vstring(
        )
)

process.TFileService = cms.Service("TFileService",
                                   fileName = cms.string("test.root")
)

process.trigger=diphotonDumper.clone()
process.id=diphotonDumper.clone()

process.triggerMva=diphotonDumper.clone(src=cms.InputTag("hmvaDiPhotons"))
process.mva=diphotonDumper.clone(src=cms.InputTag("hmvaDiPhotons"))

process.egid=diphotonDumper.clone(src=cms.InputTag("egLooseDiPhotons"))

process.kin=diphotonDumper.clone(src=cms.InputTag("kinDiPhotons"))
process.kin.dumpTrees = True

process.isoKinDiphotons = process.tmpKinDiPhotons.clone(src="kinDiPhotons",
                                                     cut="leadingPhoton.userFloat('genIso') < 10. && subLeadingPhoton.userFloat('genIso') < 10.")
process.isohCic4Diphotons = process.tmpKinDiPhotons.clone(src="hcic4DiPhotons",
                                                       cut="leadingPhoton.userFloat('genIso') < 10. && subLeadingPhoton.userFloat('genIso') < 10.")
process.isohmvaDiphotons = process.tmpKinDiPhotons.clone(src="hmvaDiPhotons",