process.trignjets = cuts.clone()
process.trignjetsht500 = cuts.clone(min_ht = 500)
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
cs = CRABSubmitter('EventsCutplay',
total_number_of_events = 1000000,
events_per_job = 500000,
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)
import os, sys
from JMTucker.Tools.BasicAnalyzer_cfg import cms, process, add_analyzer, geometry_etc
process.MessageLogger.cerr.FwkReport.reportEvery = 1
geometry_etc(process, 'START53_V27::All')
process.source.fileNames = ['file:reco.root']
process.source.secondaryFileNames = cms.untracked.vstring('file:gensimhlt.root')
process.load('SimTracker.TrackAssociation.TrackAssociatorByHits_cfi')
process.TFileService.fileName = 'TrackAssociator.root'
process.TrackAssociator = cms.EDAnalyzer('MFVTrackAssociator',
tracks_src = cms.InputTag('generalTracks'),
sim_tracks_src = cms.InputTag('g4SimHits'),
sim_vertices_src = cms.InputTag('g4SimHits'),
tracking_particles_src = cms.InputTag('mergedtruth','MergedTrackTruth'),
gen_particles_src = cms.InputTag('genParticles'),
rec_vertices_src = cms.InputTag('mfvInclusiveVertexFinder'),
do_checks = cms.bool(True),
)
process.load('JMTucker.MFVNeutralino.GenParticles_cff')
process.load('JMTucker.Tools.ParticleListDrawer_cff')
process.ParticleListDrawer2 = process.ParticleListDrawer.clone(src = cms.InputTag('mfvGenParticles', 'All'))
process.ParticleListDrawer3 = process.ParticleListDrawer.clone(src = cms.InputTag('mfvGenParticles', 'Visible'))
process.load('CommonTools.ParticleFlow.goodOfflinePrimaryVertices_cfi')
process.goodOfflinePrimaryVertices.filter = cms.bool(True)
process.load('JMTucker.MFVNeutralino.VertexReco_cff')
import os, sys
from JMTucker.Tools.BasicAnalyzer_cfg import cms, process
from JMTucker.Tools import SampleFiles
use_weights = False
SampleFiles.setup(process, 'MFVNtupleV15', 'qcdht1000', 100000)
process.TFileService.fileName = 'pileup.root'
process.load('JMTucker.MFVNeutralino.VertexSelector_cfi')
vtx_sel = process.mfvSelectedVerticesTight
process.load('JMTucker.MFVNeutralino.AnalysisCuts_cfi')
ana_sel = process.mfvAnalysisCuts
def pize(f,sz):
fmt = '%.' + str(sz) + 'f'
return (fmt % f).replace('.','p').replace('-','n')
changes = []
changes.append(('nm1', '', ''))
changes.append(('tightX3', '', 'max_npv = 5'))
changes.append(('ABCDX3', 'min_ntracks = 5, min_maxtrackpt = 0', 'min_ntracks01 = 0, min_maxtrackpt01 = 0, max_npv = 5'))
changes.append(('triggerX3', '', 'apply_vertex_cuts = False, max_npv = 5'))
changes.append(('nocutsX3', '', 'min_4th_jet_pt = 0, min_njets = 0, apply_vertex_cuts = False, max_npv = 5, min_ht = 0'))
for i in xrange(6,30,3):
changes.append(('tightX%i'%i, '', 'min_npv = %i, max_npv = %i+2'%(i,i)))
changes.append(('ABCDX%i'%i, 'min_ntracks = 5, min_maxtrackpt = 0', 'min_ntracks01 = 0, min_maxtrackpt01 = 0, min_npv = %i, max_npv = %i+2'%(i,i)))
changes.append(('triggerX%i'%i, '', 'apply_vertex_cuts = False, min_npv = %i, max_npv = %i+2'%(i,i)))
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)
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)
########################################################################
from HLTrigger.HLTfilters.hltHighLevel_cfi import hltHighLevel
process.goodDataFilter = hltHighLevel.clone()
process.goodDataFilter.TriggerResultsTag = cms.InputTag('TriggerResults', '', 'PAT')
process.goodDataFilter.HLTPaths = ['eventCleaningAll'] # can set to just 'goodOfflinePrimaryVertices', for example
process.goodDataFilter.andOr = False # = AND
process.triggerFilter = hltHighLevel.clone()
process.triggerFilter.HLTPaths = ['HLT_QuadJet50_v*']
process.triggerFilter.andOr = True # = OR
########################################################################
process.load('JMTucker.Tools.ResolutionsHistogrammer_cfi')
for x in ['WithTrigger']:
setattr(process, 'histos' + x, process.histos.clone())
########################################################################
process.p0 = cms.Path( process.histos)
process.p1 = cms.Path(process.triggerFilter * process.histosWithTrigger)
if 'debug' in sys.argv:
from JMTucker.Tools.CMSSWTools import file_event_from_argv
file_event_from_argv(process)
process.MessageLogger.cerr.FwkReport.reportEvery = 1
process.histos.print_info = True
process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
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
cs = CRABSubmitter(
'EventsCutplay',
total_number_of_events=1000000,
nonpucands_src = cms.InputTag('pfNoPileUpPF'),
pf_jets_src = cms.InputTag('ak4PFJets'),
pat_jets_src = cms.InputTag('selectedPatJetsPF'),
pv_src = cms.InputTag('goodOfflinePrimaryVertices'),
ltmm_src = cms.InputTag('mfvTrackMatches'),
pt_cut = cms.double(0),
verbose = cms.bool(True),
)
process.PileupRemovalPlayNoZCheck = process.PileupRemovalPlay.clone(pucands_src = 'pfPileUpNoClosestZVertexPF',
nonpucands_src = 'pfNoPileUpNoClosestZVertexPF')
process.PileupRemovalPlayPt5 = process.PileupRemovalPlay.clone(pt_cut = 5)
process.PileupRemovalPlayNoZCheckPt5 = process.PileupRemovalPlayNoZCheck.clone(pt_cut = 5)
#process.pPileupRemovalPlay *= process.PileupRemovalPlayNoZCheck * process.PileupRemovalPlayPt5 * process.PileupRemovalPlayNoZCheckPt5
process.load('JMTucker.MFVNeutralino.RedoPURemoval_cff')
process.pPileupRemovalPlay.replace(process.PileupRemovalPlay, process.mfvRedoPURemoval * process.PileupRemovalPlay)
if __name__ == '__main__' and hasattr(sys, 'argv') and 'submit' in sys.argv:
from JMTucker.Tools.Samples import *
mfv_neutralino_tau0000um_M0400.ana_dataset_override = '/mfv_neutralino_tau0000um_M0400/jchu-jtuple_pileupremovalstudies_v6-5c257295836d18fbf020b73449ad9b5f/USER'
mfv_neutralino_tau0100um_M0400.ana_dataset_override = '/mfv_neutralino_tau0100um_M0400/jchu-jtuple_pileupremovalstudies_v6-5c257295836d18fbf020b73449ad9b5f/USER'
mfv_neutralino_tau1000um_M0400.ana_dataset_override = '/mfv_neutralino_tau1000um_M0400/jchu-jtuple_pileupremovalstudies_v6-5c257295836d18fbf020b73449ad9b5f/USER'
mfv_neutralino_tau1000um_M1000.ana_dataset_override = '/mfv_neutralino_tau1000um_M1000/jchu-jtuple_pileupremovalstudies_v6-5c257295836d18fbf020b73449ad9b5f/USER'
mfv_neutralino_tau9900um_M0400.ana_dataset_override = '/mfv_neutralino_tau9900um_M0400/jchu-jtuple_pileupremovalstudies_v6-5c257295836d18fbf020b73449ad9b5f/USER'
samples = [mfv_neutralino_tau0000um_M0400, mfv_neutralino_tau0100um_M0400, mfv_neutralino_tau1000um_M0400, mfv_neutralino_tau1000um_M1000, mfv_neutralino_tau9900um_M0400]
from JMTucker.Tools.CRABSubmitter import CRABSubmitter
cs = CRABSubmitter('PileupPlay',
import os, sys
from JMTucker.Tools.BasicAnalyzer_cfg import cms, process
process.source.fileNames = ['/store/user/jchu/mfv_neutralino_tau1000um_M0400/jtuple_v7/5d4c2a74c85834550d3f9609274e8548/pat_1_1_hdB.root']
process.source.secondaryFileNames = cms.untracked.vstring('/store/user/tucker/mfv_neutralino_tau1000um_M0400/mfv_neutralino_tau1000um_M0400/a6ab3419cb64660d6c68351b3cff9fb0/aodpat_891_1_sZ9.root','/store/user/tucker/mfv_neutralino_tau1000um_M0400/mfv_neutralino_tau1000um_M0400/a6ab3419cb64660d6c68351b3cff9fb0/aodpat_948_2_lgB.root')
process.TFileService.fileName = 'leptonic_nb_analysis.root'
process.load('Configuration.Geometry.GeometryIdeal_cff')
process.load('Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
process.GlobalTag.globaltag = 'START53_V21::All'
process.load('TrackingTools.TransientTrack.TransientTrackBuilder_cfi')
process.load('CommonTools.ParticleFlow.goodOfflinePrimaryVertices_cfi')
process.goodOfflinePrimaryVertices.filter = cms.bool(True)
process.load('JMTucker.MFVNeutralino.Vertexer_cff')
process.p = cms.EDAnalyzer('LeptonicNBAnalysis',
trigger_results_src = cms.InputTag('TriggerResults', '', 'HLT'),
primary_vertex_src = cms.InputTag('goodOfflinePrimaryVertices'),
secondary_vertex_src = cms.InputTag('mfvSelectedVerticesTight'),
muon_src = cms.InputTag('selectedPatMuonsPF'),
min_nmuons = cms.int32(1),
jet_src = cms.InputTag('selectedPatJetsPF'),
min_jet_pt = cms.double(30.0),
max_jet_eta = cms.double(2.4),
min_njets = cms.int32(6),
b_discriminator_name = cms.string('combinedSecondaryVertexBJetTags'),
bdisc_min = cms.double(0.679),
min_nbtags = cms.int32(1),
muon_stations_min = cms.bool(True),
import sys
from JMTucker.Tools.BasicAnalyzer_cfg import cms, process, geometry_etc
from JMTucker.Tools import SampleFiles
#SampleFiles.setup(process, 'MFVNtupleV18', 'mfv_neutralino_tau1000um_M0400', 100)
process.source.fileNames = ['/store/user/tucker/TTJets_HadronicMGDecays_8TeV-madgraph/mfvntuple_v18/c761ddfa7f093d8f86a338439e06a1d4/ntuple_100_1_NHs.root']
process.TFileService.fileName = 'histos.root'
process.load('JMTucker.MFVNeutralino.Histos_cff')
import JMTucker.Tools.SimpleTriggerResults_cfi as SimpleTriggerResults
SimpleTriggerResults.setup_endpath(process, weight_src='mfvWeight')
nm1s = [
('Ntracks', 'min_ntracks = 0, min_njetsntks = 0'),
('Drmin', 'max_drmin = 1e9'),
('Drmax', 'max_drmax = 1e9'),
('Mindrmax','min_drmax = 0'),
('Bs2derr', 'max_bs2derr = 1e9'),
('Njets', 'min_njetsntks = 0'),
('Ntracksptgt3', 'min_ntracksptgt3 = 0'),
('Sumnhitsbehind', 'max_sumnhitsbehind = 1000000'),
('ButNtracksAndGt3', 'max_drmin = 1e9, max_drmax = 1e9, min_drmax = 0, max_bs2derr = 1e9, min_njetsntks = 0'),
]
for name, cut in nm1s:
evt_cut = ''
if type(cut) == tuple:
cut, evt_cut = cut
vtx = eval('process.mfvSelectedVerticesTight.clone(%s)' % cut)
import sys
from JMTucker.Tools.BasicAnalyzer_cfg import cms, process, file_event_from_argv
file_event_from_argv(process)
process.TFileService.fileName = 'simple_trigger_efficiency.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)
cuts = {
'A': 'min_ntracks = 5',
'B': 'max_drmin = 0.4',
'C': 'max_drmax = 4',
'D': 'min_drmax = 1.2',
'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)