from JMTucker.Tools.BasicAnalyzer_cfg import cms, process
process.source.fileNames = [
'/store/user/tucker/sstop_genfsimreco_test/sstop_genfsimreco_test//15c4250952b10a469cc6da8beaecd65e/fastsim_9_1_Ziq.root'
]
process.MessageLogger.cerr.FwkReport.reportEvery = 1
process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
process.ParticleListDrawer = cms.EDAnalyzer(
'ParticleListDrawer',
maxEventsToPrint=cms.untracked.int32(1000000),
src=cms.InputTag('genParticles'),
printOnlyHardInteraction=cms.untracked.bool(False),
useMessageLogger=cms.untracked.bool(False))
process.a = cms.EDAnalyzer('MCInteractionTest')
process.p = cms.Path(process.ParticleListDrawer * process.a)
process.trignjetsht750 = cuts.clone(min_ht=750)
process.trigmu = cuts.clone(trigger_bit=3, min_4th_jet_pt=20)
process.trigmuht500 = cuts.clone(trigger_bit=3, min_4th_jet_pt=20, min_ht=500)
process.trigtopmu = cuts.clone(trigger_bit=4, min_4th_jet_pt=20)
process.trigtopmunmu = cuts.clone(trigger_bit=4,
min_4th_jet_pt=20,
min_nsemilepmuons=1)
process.trigtopmunmuht500 = cuts.clone(trigger_bit=4,
min_4th_jet_pt=20,
min_nsemilepmuons=1,
min_ht=500)
process.trignjetsslep = cuts.clone(min_nsemileptons=1)
process.trigonly = cuts.clone(min_njets=0, min_4th_jet_pt=0)
for name in process.filters.keys():
setattr(process, 'p' + name, cms.Path(getattr(process, name)))
if use_weights:
process.load('JMTucker.MFVNeutralino.WeightProducer_cfi')
import JMTucker.Tools.SimpleTriggerEfficiency_cfi as SimpleTriggerEfficiency
SimpleTriggerEfficiency.setup_endpath(
process, weight_src='mfvWeight' if use_weights else '')
if __name__ == '__main__' and hasattr(sys, 'argv') and 'submit' in sys.argv:
import JMTucker.Tools.Samples as Samples
samples = Samples.qcd_samples + Samples.ttbar_samples + Samples.mfv_signal_samples + Samples.leptonic_background_samples
from JMTucker.Tools.CRABSubmitter import CRABSubmitter
from JMTucker.Tools.SampleFiles import SampleFiles
from stopBL_samples import *
inputFiles = stopBL_M400_CTau_0p1mm
process.source = cms.Source(
"PoolSource",
fileNames=cms.untracked.vstring(*inputFiles),
)
process.TFileService.fileName = 'stopBL_M400_CTau_0p1mm.root'
process.RandomNumberGeneratorService = cms.Service(
'RandomNumberGeneratorService')
process.RandomNumberGeneratorService.SimpleTriggerEfficiency = cms.PSet(
initialSeed=cms.untracked.uint32(1219))
#import prescales
process.load('JMTucker.Tools.SimpleTriggerEfficiency_cfi')
process.SimpleTriggerEfficiency.prescale_paths = cms.vstring(
) #*prescales.prescales.keys()),
process.SimpleTriggerEfficiency.prescale_values = cms.vuint32(
) #*[o for l,h,o in prescales.prescales.itervalues()]),
process.p = cms.Path(process.SimpleTriggerEfficiency)
for x in sys.argv:
if x.startswith('process='):
process_name = x.replace('process=', '')
process.SimpleTriggerEfficiency.trigger_results_src = cms.InputTag(
'TriggerResults', '', process_name)
'E': 'max_geo2ddist = 2.5',
'F': 'max_bs2derr = 0.0025',
'G': 'min_njetsntks = 1',
'H': 'min_ntracksptgt3 = 3',
#'I': 'max_sumnhitsbehind = 0',
}
all_cuts = ''.join(sorted(cuts.keys()))
process.load('JMTucker.MFVNeutralino.AnalysisCuts_cfi')
process.ana = process.mfvAnalysisCuts.clone(apply_vertex_cuts=False)
process.load('JMTucker.MFVNeutralino.VertexSelector_cfi')
vtx_sel = process.mfvSelectedVertices
process.p = cms.Path(process.ana)
process.load('JMTucker.MFVNeutralino.WeightProducer_cfi')
process.p *= process.mfvWeight
vtx_srcs = []
for i in xrange(len(cuts) + 1):
for cut_combs in combinations(all_cuts, i):
name = ''.join(cut_combs)
if name == '':
name = 'none'
cut = ', '.join(cuts[x] for x in cut_combs)
obj = eval('vtx_sel.clone(vertex_src = "%s", %s)' % (name, cut))
print name, cut
setattr(process, name, obj)