Exemple #1
0
    'JetCorrection',
    InputName='AKt4PFJetsCHS',
    CorrectedJetsName='CorrectedJets',
    RhoAlgo=mithep.PileupEnergyDensity.kFixedGridFastjetAll)

puppiJetCorrection = mithep.JetCorrectionMod(
    'PuppiJetCorrection',
    InputName='AKt4PuppiJets',
    CorrectedJetsName='CorrectedPuppiJets',
    RhoAlgo=mithep.PileupEnergyDensity.kFixedGridFastjetAll)

metCorrection = mithep.MetCorrectionMod(
    'MetCorrection',
    InputName='PFMet',
    OutputName='PFType1CorrectedMet',
    JetsName='AKt4PFJetsCHS',
    RhoAlgo=mithep.PileupEnergyDensity.kFixedGridFastjetAll,
    MaxEMFraction=0.9,
    SkipMuons=True,
    MinJetPt=15.,
    MuonGeometricMatch=True)
metCorrection.ApplyType0(False)
metCorrection.ApplyType1(True)
metCorrection.ApplyShift(False)
metCorrection.IsData(analysis.isRealData)

metCorrectionJESUp = metCorrection.clone('MetCorrectionJESUp',
                                         OutputName='PFType1CorrectedMetJESUp',
                                         JESUncertaintySigma=1.)
metCorrectionJESDown = metCorrection.clone(
    'MetCorrectionJESDown',
    OutputName='PFType1CorrectedMetJESDown',
jetId = mithep.JetIdMod(
    InputName = jetCorr.GetOutputName(),
    OutputName = 'CleanedJets',
    UseL1Correction = False,
    UseL2Correction = False,
    UseL3Correction = False,
    ApplyPFLooseId = True,
    ApplyMVACut = True,
    PtMin = 30.,
    EtaMax = 2.5
)

metCorr = mithep.MetCorrectionMod(
    CorrectedName = 'PFType1CorrectedMet',
    JetsName = 'AKt4PFJetsCHS',
    CorrectedJetsName = jetCorr.GetOutputName()
)
metCorr.ApplyType1(True)
metCorr.IsData(False)

preRun2Sync = mithep.PreRun2SynchExercise(
    VerticesName = mithep.Names.gkPVBrn,
    MetName = metCorr.GetOutputName(),
    JetsName = jetId.GetOutputName(),
    ElectronsName = eleId.GetOutputName(),
    MuonsName = muId.GetOutputName(),
    TausName = tauId.GetOutputName(),
    PhotonsName = photonId.GetOutputName()
)