Esempio n. 1
0
    ak8JetExtender,
    goodCA15Jets,
    puppiAK8CorrectionMod,
    puppiCA15CorrectionMod,
    goodAK8PuppiJets,
    goodCA15PuppiJets,
])

#############################
### DATA/MC CUSTOMIZATION ###
#############################

if analysis.isRealData:
    badEventsFilterMod = mithep.BadEventsFilterMod(
        'BadEventsFilterMod',
        FillHist=True,
        TaggingMode=False  # filter out events
    )
    badEventsFilterMod.SetFilter('HBHENoiseFilter')
    badEventsFilterMod.SetFilter('EEBadScFilter')
    if int(analysis.book.split('/')[2]) < 43:  # e.g. t2mit/filefi/042
        badEventsFilterMod.AddEventList(
            'CSCTightHaloFilter', mitdata + '/eventlist/csc2015_Dec01.txt')
        badEventsFilterMod.AddEventList(
            'EEBadScFilter', mitdata + '/eventlist/ecalscn1043093_Dec01.txt')
        badEventsFilterMod.AddEventList(
            'CHTrackResolutionFilter',
            mitdata + '/eventlist/badResolutionTrack_Jan13.txt')
        badEventsFilterMod.AddEventList(
            'MuBadTrackFilter', mitdata + '/eventlist/muonBadTrack_Jan13.txt')
        badEventsFilterMod.AddEventList(
Esempio n. 2
0
if analysis.isRealData:
    jecPattern = mitdata + '/JEC/Summer15_' + jecVersion + '/Summer15_' + jecVersion + '_DATA_{level}_{jettype}.txt'
    jecLevels = ['L1FastJet', 'L2Relative', 'L3Absolute', 'L2L3Residual']

else:
    jecPattern = mitdata + '/JEC/Summer15_' + jecVersion + '/Summer15_' + jecVersion + '_MC_{level}_{jettype}.txt'
    jecLevels = ['L1FastJet', 'L2Relative', 'L3Absolute']

#########################################
### MODULES RUN WITH DEFAULT SETTINGS ###
#########################################

badEventsFilterMod = mithep.BadEventsFilterMod(
    'BadEventsFilterMod',
    EEBadScFilter=True,
    CSCTightHaloFilter=True,
    HBHENoiseFilter=analysis.isRealData,
    FillHist=True)
from MitPhysics.Mods.GoodPVFilterMod import goodPVFilterMod
from MitPhysics.Mods.SeparatePileUpMod import separatePileUpMod
from MitPhysics.Mods.PuppiMod import puppiMod
from MitPhysics.Mods.PuppiPFJetMod import puppiPFJetMod

################################
### JET/MET ID & CORRECTIONS ###
################################

jetCorrection = mithep.JetCorrectionMod(
    'JetCorrection',
    InputName='AKt4PFJetsCHS',
    CorrectedJetsName='CorrectedJets',