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')
