def getit(fn): f = ROOT.TFile(fn) t = f.Get('mfvMiniTree/t') hr = draw_hist_register(t, True) h,n = hr.draw('weight', 'nvtx>=2', binning='1,0,1', get_n=True, goff=True) i = get_integral(h)[0] d = Samples.norm_from_file(f) e,l,u = clopper_pearson(i, d) ee = (u-l)/2 return e,ee,i,d
def getit(fn): f = ROOT.TFile(fn) t = f.Get('mfvMiniTree/t') hr = draw_hist_register(t, True) h, n = hr.draw('weight', 'nvtx>=2' + ('&& genmatch0 && genmatch1' if req_genmatch else ''), binning='1,0,1', get_n=True, goff=True) #i = get_integral(h)[0] #c = i/n d = Samples.norm_from_file(f) e, l, u = clopper_pearson(n, d) ee = (u - l) / 2 return e, ee, n, d
'drop *', 'keep MFVEvent_mfvEvent__*', 'keep MFVVertexAuxs_mfvSelectedVerticesTight__*', 'keep *_mfvSampleInfo_*_*', ) ) process.outp = cms.EndPath(process.out) process.options.emptyRunLumiMode = cms.untracked.string('doNotHandleEmptyRunsAndLumis') process.source.inputCommands = cms.untracked.vstring('keep *', 'drop *_MEtoEDMConverter_*_*') process.out.outputCommands += ['drop LumiDetails_lumiProducer_*_*', 'drop LumiSummary_lumiProducer_*_*', 'drop RunSummary_lumiProducer_*_*'] process.out.dropMetaData = cms.untracked.string('ALL') if __name__ == '__main__' and hasattr(sys, 'argv') and 'submit' in sys.argv: samples = Samples.from_argv([Samples.mfv_neutralino_tau0100um_M0400, Samples.mfv_neutralino_tau1000um_M0400, Samples.mfv_neutralino_tau0300um_M0400, Samples.mfv_neutralino_tau9900um_M0400] + Samples.ttbar_samples + Samples.qcd_samples) run_half = True def modify(sample): to_add = [] to_replace = [] to_add.append(''' process.mfvSampleInfo.sample = '%s' process.mfvSampleInfo.num_events = %s process.mfvSampleInfo.cross_section = %g ''' % (sample.name, 'int(%i*0.5)' % sample.nevents if (run_half and sample in Samples.ttbar_samples + Samples.qcd_samples) else sample.nevents, sample.cross_section)
#!/usr/bin/env python import sys from JMTucker.Tools.Merge_cfg import cms, process if __name__ == '__main__' and hasattr(sys, 'argv') and 'submit' in sys.argv: import JMTucker.Tools.Samples as Samples samples = Samples.from_argv(Samples.mfv_signal_samples + [Samples.qcdht0100]) from JMTucker.Tools.CRABSubmitter import CRABSubmitter cs = CRABSubmitter('MergeNtupleV18_2', use_ana_dataset = True, total_number_of_events = -1, events_per_job = 50000, get_edm_output = True, data_retrieval = 'fnal', max_threads = 3, publish_data_name = 'mfvmergentuple_v18', ) cs.submit_all(samples)
size=(600, 600), pdf=True, log=False) multijet = Samples.mfv_signal_samples_2017 dijet = Samples.mfv_stopdbardbar_samples_2017 for sample in multijet + dijet: fn = os.path.join('/uscms_data/d2/tucker/crab_dirs/MiniTree%s' % version, sample.name + '.root') if not os.path.exists(fn): print 'no', sample.name continue f = ROOT.TFile(fn) t = f.Get('mfvMiniTree/t') hr = draw_hist_register(t, True) cut = 'nvtx>=2' # && svdist > 0.04' h = hr.draw('weight', cut, binning='1,0,1', goff=True) num, _ = get_integral(h) den = Samples.norm_from_file(f) sample.y, sample.yl, sample.yh = clopper_pearson( num, den) # ignore integral != entries, just get central value right print '%26s: efficiency = %.3f (%.3f, %.3f)' % (sample.name, sample.y, sample.yl, sample.yh) per = PerSignal('efficiency', y_range=(0., 1.05)) per.add(multijet, title='#tilde{N} #rightarrow tbs') per.add(dijet, title='#tilde{t} #rightarrow #bar{d}#bar{d}', color=ROOT.kBlue) per.draw(canvas=ps.c) ps.save('sigeff')
import FWCore.ParameterSet.Config as cms process = cms.Process("PVAnalyzer") process.load("FWCore.MessageService.MessageLogger_cfi") process.MessageLogger.cerr.INFO = cms.untracked.PSet(limit = cms.untracked.int32(-1)) process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(-1)) process.source = cms.Source("PoolSource", fileNames = cms.untracked.vstring('file:/store/mc/Summer12_DR53X/QCD_HT-1000ToInf_TuneZ2star_8TeV-madgraph-pythia6/AODSIM/PU_S10_START53_V7A-v1/00000/16716FE6-440E-E211-A50B-0024E8767D86.root')) process.Analyzer = cms.EDAnalyzer('PVAnalyzer',primary_vertices_src = cms.InputTag('offlinePrimaryVertices'), track_src = cms.InputTag('generalTracks'), beamspot_src = cms.InputTag('offlineBeamSpot'), use_only_pv_tracks = cms.bool(True), use_only_pvs_tracks = cms.bool(True), maxNormChi2 = cms.double(20.0), minPxLayer = cms.int32(2), minSilLayer = cms.int32(5) ) process.TFileService = cms.Service("TFileService", fileName = cms.string('PVAnalyzer_histo.root')) process.p = cms.Path(process.Analyzer) if __name__ == '__main__' and hasattr(sys, 'argv') and 'submit' in sys.argv: import JMTucker.Tools.Samples as Samples samples = Samples.from_argv([Samples.qcdht1000, Samples.mfv_neutralino_tau1000um_M0400, Samples.mfv_neutralino_tau0100um_M0400]) from JMTucker.Tools.CRABSubmitter import CRABSubmitter cs = CRABSubmitter('PVAnalyzer', total_number_of_events = -1, events_per_job = 100000, ) cs.submit_all(samples)
import sys from JMTucker.Tools.BasicAnalyzer_cfg import * file_event_from_argv(process) process.TFileService.fileName = cms.string('evids.root') process.evids = cms.EDAnalyzer('EventIdRecorder') process.p = cms.Path(process.evids) process.source.duplicateCheckMode = cms.untracked.string('noDuplicateCheck') if __name__ == '__main__' and hasattr(sys, 'argv') and 'submit' in sys.argv: import JMTucker.Tools.Samples as Samples samples = Samples.from_argv(Samples.ttbar_samples + Samples.qcd_samples + Samples.mfv_signal_samples) from JMTucker.Tools.CRABSubmitter import CRABSubmitter cs = CRABSubmitter('EventIdsV18', total_number_of_events = -1, events_per_job = 1000000, use_ana_dataset = True, ) cs.submit_all(samples)
process.source.secondaryFileNames = cms.untracked.vstring(*'''/store/mc/Summer12_DR53X/TTJets_HadronicMGDecays_8TeV-madgraph/AODSIM/PU_S10_START53_V7A-v1/00000/0E9DCF01-0216-E211-934F-20CF3019DF0F.root /store/mc/Summer12_DR53X/TTJets_HadronicMGDecays_8TeV-madgraph/AODSIM/PU_S10_START53_V7A-v1/00000/E0F2FA59-1016-E211-A573-00259073E3A8.root /store/mc/Summer12_DR53X/TTJets_HadronicMGDecays_8TeV-madgraph/AODSIM/PU_S10_START53_V7A-v1/00000/5205EEA5-1016-E211-BBE1-90E6BA442F2B.root'''.split('\n')) process.load('JMTucker.MFVNeutralino.VertexSelector_cfi') process.load('JMTucker.MFVNeutralino.AnalysisCuts_cfi') process.mfvAnalysisCutsOneVtx = process.mfvAnalysisCuts.clone(min_nvertex = 1) process.WhichGenParticleOneVtx = cms.EDAnalyzer('MFVWhichGenParticle', gen_particles_src = cms.InputTag('genParticles'), mevent_src = cms.InputTag('mfvEvent'), vertices_src = cms.InputTag('mfvSelectedVerticesTight'), ) process.WhichGenParticleTwoVtx = process.WhichGenParticleOneVtx.clone() process.p = cms.Path(process.mfvSelectedVerticesSeq * process.mfvAnalysisCutsOneVtx * process.WhichGenParticleOneVtx * process.mfvAnalysisCuts * process.WhichGenParticleTwoVtx) if __name__ == '__main__' and hasattr(sys, 'argv') and 'submit' in sys.argv: import JMTucker.Tools.Samples as Samples samples = Samples.from_argv(Samples.ttbar_samples + Samples.qcd_samples) from JMTucker.Tools.CRABSubmitter import CRABSubmitter cs = CRABSubmitter('WhichGenParticle', job_control_from_sample = True, use_ana_dataset = True, use_parent = True, run_half_mc = True, ) cs.submit_all(samples)
err = 'trying to submit on data, and tuple template does not contain the magic string "%s"' % magic to_replace.append((magic, 'runOnMC = False', err)) to_add.append('process.dummyToMakeDiffHash = cms.PSet(sampleName = cms.string("%s"))' % sample.name) to_add.append('process.patJetCorrFactors.levels.append("L2L3Residual")') return to_add, to_replace cs = CRABSubmitter('QuadJetTrigEff', pset_modifier = modify, job_control_from_sample = True, max_threads = 2, ) mc_samples = [Samples.qcdmupt15] + Samples.ttbar_samples + Samples.leptonic_background_samples data_samples = [ Samples.DataSample('SingleMu2012A', '/SingleMu/Run2012A-22Jan2013-v1/AOD'), Samples.DataSample('SingleMu2012B', '/SingleMu/Run2012B-22Jan2013-v1/AOD'), Samples.DataSample('SingleMu2012C', '/SingleMu/Run2012C-22Jan2013-v1/AOD'), Samples.DataSample('SingleMu2012D', '/SingleMu/Run2012D-22Jan2013-v1/AOD'), ] for sample in mc_samples: sample.events_per = 20000 for sample in data_samples: sample.lumis_per = 50 samples = Samples.from_argv(mc_samples + data_samples) cs.submit_all(samples)
process.p *= obj vtx_srcs.append(cms.InputTag(name)) process.nmx = cms.EDAnalyzer( 'NmxHistos', weight_src=cms.InputTag('mfvWeight'), use_weight=cms.bool(False), vertex_srcs=cms.VInputTag(*vtx_srcs), ) process.p *= process.nmx if __name__ == '__main__' and hasattr(sys, 'argv') and 'submit' in sys.argv: import JMTucker.Tools.Samples as Samples samples = Samples.from_argv([ Samples.mfv_neutralino_tau0100um_M0400, Samples. mfv_neutralino_tau1000um_M0400, Samples.mfv_neutralino_tau0300um_M0400, Samples.mfv_neutralino_tau9900um_M0400 ] + Samples.ttbar_samples + Samples.qcd_samples + Samples.data_samples) for s in Samples.data_samples: s.json = 'ana_all.json' from JMTucker.Tools.CRABSubmitter import CRABSubmitter cs = CRABSubmitter( 'NmxV20', job_control_from_sample=True, use_ana_dataset=True, ) cs.submit_all(samples)