import ROOT,math,sys from DataFormats.FWLite import Events,Handle from ROOT import TLorentzVector from selectionFunction import myClass #---------------Information--From--TTTTMCFile------------------------------------------ events = Events('/user/Samples/TTTTJetsMC/TTTT_TuneZ2star_8TeV-madgraph-tauola_patTuple.root') test = myClass(events,"TTTTHist") test.selectionloop() print "-----------------------------Selection-----------------------------------------" print "Total Events ",test.numbersignalevents print "triggeredEvts ",test.numbertrigger print "muonsEvts ",test.numbertotalmuons print '---6jets Criteria------' print "jetsEvts ",test.numbertotalsixjets print "wmasscriteria ",test.wmassselected print '---7jets Criteria------' print "jetsEvts ",test.numbertotalsevenjets print "jetsEvts+highenergetic",test.numbertotalsevenjetsmomentum print "jetsEvts + Wmass criteria ",test.wmassselected7 print "Number Comb Jets ",test.totaljetcombevents test.jetsmomentumbefore.SaveAs(test.jetsmomentumbefore.GetName()+'_mc.root') test.jetsmomentumafter.SaveAs(test.jetsmomentumafter.GetName()+'_mc.root') test.numberjetsBeforeDataHist.SaveAs(test.numberjetsBeforeDataHist.GetName()+'_mc.root') test.numberjetsAfterDataHist.SaveAs(test.numberjetsAfterDataHist.GetName()+'_mc.root') test.wmassdistribution.SaveAs(test.wmassdistribution.GetName()+'.root') test.topmassdistribution.SaveAs(test.topmassdistribution.GetName()+'.root') test.cutprocess.SaveAs(test.cutprocess.GetName()+'.root')
import ROOT,math,sys from DataFormats.FWLite import Events,Handle from ROOT import TLorentzVector from selectionFunction import myClass multiJetFiles = open('/home/home2/institut_3b/almazan/MultiJetfiles.txt') files=multiJetFiles.readlines() multiJetFiles=[] for f in files: multiJetFiles.append ('dcap://grid-dcap-extern.physik.rwth-aachen.de/pnfs/physik.rwth-aachen.de/cms'+f.rstrip('\n')) events = Events(multiJetFiles) test = myClass(events,"DataHist") test.selectionloop() print 'Number of selected events is ',test.numberselectedevents, 'in a total of ',test.numbersignalevents,' events' print "triggeredEvts ",test.numbertrigger print "muonsEvts ",test.numbertotalmuons print "jetsEvts ",test.numbertotaljets print "btagEvts ",test.numbertotalbjets
#multiJetFiles = open('/user/almazan/SingleMu__Run2012A-22Jan2013-v1__AOD_gridOutputFiles.txt') #multiJetFiles = open('/home/home2/institut_3b/almazan/SingleMu.txt') #---------------Information--From--DataRunA------------------------------------------ multiJetFiles = open('/user/almazan/CMSSW_5_3_16/FWLitePython/SingleMu_SoftJetSkim.txt') files=multiJetFiles.readlines() multiJetFiles=[] for f in files: multiJetFiles.append ('dcap://grid-dcap-extern.physik.rwth-aachen.de/pnfs/physik.rwth-aachen.de/cms'+f.rstrip('\n')) print "before creating Events" events = Events(multiJetFiles) print "after creating Events" test = myClass(events,"DataHistRunA") test.selectionloop() print "-----------------------------Selection-----------------------------------------" print "triggeredEvts ",test.numbertrigger print "muonsEvts ",test.numbertotalmuons print '---6jets Criteria------' print "jetsEvts ",test.numbertotalsixjets print "wmasscriteria ",test.wmassselected print '---7jets Criteria------' print "jetsEvts ",test.numbertotalsevenjets print "jetsEvts+highenergetic",test.numbertotalsevenjetsmomentum print "jetsEvts + Wmass criteria ",test.wmassselected7 print "Number Comb Jets ",test.totaljetcombevents test.jetsmomentumbefore.SaveAs(test.jetsmomentumbefore.GetName()+'_mc.root') test.jetsmomentumafter.SaveAs(test.jetsmomentumafter.GetName()+'_mc.root')