Exemplo n.º 1
0
#!/usr/bin/env python

import FWCore.ParameterSet.Config as cms

from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cfg import process
from SUSYBSMAnalysis.Zprime2muAnalysis.HistosFromPAT_cfi import HistosFromPAT
from SUSYBSMAnalysis.Zprime2muAnalysis.ResolutionUsingMC_cfi import ResolutionUsingMC

levels = ['Global', 'TkOnly', 'TPFMS', 'Picky', 'TMR', 'TuneP', 'SigmaSwitch']

from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cff import rec_levels
rec_levels(process, levels)

useMC = True
process.p = cms.Path(process.Zprime2muAnalysisSequence)
for level in levels:
    lepton_src = cms.InputTag('leptons', level)
    dilepton_src = cms.InputTag('dimuons' + level)

    h = HistosFromPAT.clone(lepton_src=lepton_src, dilepton_src=dilepton_src)
    setattr(process, 'HistosFromPAT' + level, h)

    if useMC:
        r = ResolutionUsingMC.clone(lepton_src=lepton_src,
                                    dilepton_src=dilepton_src)
        setattr(process, 'ResolutionUsingMC' + level, r)

    process.p *= h * r
Exemplo n.º 2
0
HistosFromPAT.leptonsFromDileptons = True
#HistosFromPAT.useMadgraphWeight= True #if we want to re-wigth madgraph samples. For other samples the weight is always 1. Does not work with data.
process.load('SUSYBSMAnalysis.Zprime2muAnalysis.HistosFromPAT_cfi')
process.HistosFromPAT.leptonsFromDileptons = True
# These modules define the basic selection cuts. For the monitoring
# sets below, we don't need to define a whole new module, since they
# just change one or two cuts -- see below.

import SUSYBSMAnalysis.Zprime2muAnalysis.OurSelectionDec2012_cff as OurSelectionDec2012

from SUSYBSMAnalysis.Zprime2muAnalysis.hltTriggerMatch_cfi import trigger_match, prescaled_trigger_match, trigger_paths, prescaled_trigger_paths, overall_prescale, offline_pt_threshold, prescaled_offline_pt_threshold

##################
from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cff import rec_levels, rec_level_module
tracks = ['inner', 'tunepnew', 'startup']
rec_levels(process, tracks)
##################
#Zprime2muAnalysis_cfg imports #Zprime2muAnalysis_cfg where is defined Zprime2muAnalysisSequence = cms.Sequence(muonPhotonMatch * leptons * allDimuons * dimuons)
# full sequence process.p = cms.Path(process.goodDataFilter*process.muonPhotonMatch * process.leptons * process.allDimuons * process.dimuons * process.HistosFromPAT)

process.p = cms.Path(process.goodDataFilter *
                     process.Zprime2muAnalysisSequence)

process.p *= rec_level_module(process, process.HistosFromPAT, 'Histos', tracks)


def rec_level_module_ntupler(process, module, name, tracks):
    p = []
    for t in tracks:
        h = module.clone()
        if hasattr(h, 'dimu_src'):
process.maxEvents.input = -1
#process.source.fileNames = ['file:crab1/crab_effres_zp1000/res/pat_2_1_375.root']
process.source.fileNames = ['/store/user/slava/DYToMuMu_M_20_TuneZ2star_8TeV_pythia6/effres_dy20/20941d9c676d6826327c8223aa3d20e0/pat_9_1_5zY.root']
process.options.wantSummary = True

ex = ''

if use_old_selection:
    ex += 'oldsel'
    from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cff import switch_to_old_selection
    switch_to_old_selection(process)

from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cff import rec_levels, rec_level_module
tracks = ['global', 'inner', 'tpfms', 'picky', 'tunep', 'tmr', 'tunepnew']
rec_levels(process, tracks)

process.load('SUSYBSMAnalysis.Zprime2muAnalysis.HardInteractionFilter_cfi')
process.HardInteractionFilterRes = process.HardInteractionFilter.clone(use_resonance_mass=True)

process.load('SUSYBSMAnalysis.Zprime2muAnalysis.EfficiencyFromMC_cfi')

#process.HardInteractionFilter.use_resonance_mass = True
#process.EfficiencyFromMC.use_resonance_mass_denom = True

# Since LooseTightPairSelector ignores the cutFor that
# Zprime2muLeptonProducer sets, don't need to redo leptons for the
# VBTF path.
import SUSYBSMAnalysis.Zprime2muAnalysis.VBTFSelection_cff as VBTFSelection
process.allDimuonsVBTF = VBTFSelection.allDimuons.clone()
process.dimuonsVBTF = VBTFSelection.dimuons.clone(src = 'allDimuonsVBTF')
#!/usr/bin/env python

import FWCore.ParameterSet.Config as cms

from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cfg import process
from SUSYBSMAnalysis.Zprime2muAnalysis.HistosFromPAT_cfi import HistosFromPAT
from SUSYBSMAnalysis.Zprime2muAnalysis.ResolutionUsingMC_cfi import ResolutionUsingMC

levels = ['Global', 'TkOnly', 'TPFMS', 'Picky', 'TMR', 'TuneP', 'SigmaSwitch']

from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cff import rec_levels
rec_levels(process, levels)

useMC = True
process.p = cms.Path(process.Zprime2muAnalysisSequence)
for level in levels:
    lepton_src = cms.InputTag('leptons', level)
    dilepton_src = cms.InputTag('dimuons' + level)
    
    h = HistosFromPAT.clone(lepton_src=lepton_src, dilepton_src=dilepton_src)
    setattr(process, 'HistosFromPAT' + level, h)

    if useMC:
        r = ResolutionUsingMC.clone(lepton_src=lepton_src, dilepton_src=dilepton_src)
        setattr(process, 'ResolutionUsingMC' + level, r)

    process.p *= h*r