def theConfig(Nj, mcdir = '', initFile = '', btag = False, toydataFile=""):
fitterPars = RooWjjFitterParams()
fitterPars.smoothingOrder = 1
# fitterPars.MCDirectory = "/uscms_data/d2/kalanand/WjjTrees/Full2011DataFall11MC/ReducedTree_PAT/RD_"
fitterPars.MCDirectory = "/uscms_data/d2/kalanand/WjjTrees/Full2011DataFall11MC/OriginalTree_PAT/"
fitterPars.WpJDirectory = fitterPars.MCDirectory
if (len(mcdir) > 0):
fitterPars.WpJDirectory = mcdir
fitterPars.toyWpJ = False
##fitterPars.QCDDirectory = fitterPars.MCDirectory[:-3]
fitterPars.QCDDirectory = "/uscms_data/d2/kalanand/WjjTrees/Full2011DataFall11MC/ReducedTree_PAT/"
fitterPars.initParamsFile = initFile
fitterPars.constraintParamsFile = initFile
fitterPars.DataDirectory = fitterPars.MCDirectory
fitterPars.muonData = 'WmunuJets_DataAll_GoldenJSON_4p7invfb.root'
fitterPars.includeMuons = True
fitterPars.electronData = 'WenuJets_DataAllSingleElectronTrigger_GoldenJSON_4p7invfb.root'
fitterPars.includeElectrons = True
fitterPars.NewPhysicsDirectory = fitterPars.MCDirectory
fitterPars.toydataFile = toydataFile
fitterPars.minMass = 40.
fitterPars.maxMass = 150.
fitterPars.nbins = 32
fitterPars.intLumi = 5020.
binEdge = fitterPars.minMass
fitterPars.binEdges.push_back(binEdge)
print binEdge,' ',
while (binEdge < fitterPars.maxMass):
binEdge += round(0.10*binEdge)
fitterPars.binEdges.push_back(binEdge)
print binEdge,' ',
## fitterPars.binEdges[fitterPars.binEdges.size()-1] = fitterPars.maxMass
fitterPars.maxMass = fitterPars.binEdges[fitterPars.binEdges.size()-1]
print "mass range:",fitterPars.minMass,'-',fitterPars.maxMass
fitterPars.truncRange = False
fitterPars.blind = False
fitterPars.minTrunc = 65.
binMin = fitterPars.minMass
for binEdge in fitterPars.binEdges:
if binEdge > fitterPars.minTrunc:
fitterPars.minTrunc = binMin
break
binMin = binEdge
fitterPars.maxTrunc = 94.
for binEdge in fitterPars.binEdges:
if binEdge > fitterPars.maxTrunc:
fitterPars.maxTrunc = binEdge
break
fitterPars.minFit = fitterPars.minMass
binMin = fitterPars.minMass
for binEdge in fitterPars.binEdges:
if binEdge > fitterPars.minFit:
fitterPars.minFit = binMin
break
binMin = binEdge
fitterPars.maxFit = fitterPars.maxMass
for binEdge in fitterPars.binEdges:
if binEdge > fitterPars.maxFit:
fitterPars.maxFit = binEdge
break
fitterPars.minWmass = fitterPars.minTrunc
fitterPars.maxWmass = fitterPars.maxTrunc
print 'truncate:',fitterPars.truncRange,':', fitterPars.minTrunc,'-',\
fitterPars.maxTrunc
print 'fit range:',fitterPars.minFit,'-',fitterPars.maxFit
fitterPars.njets = Nj
fitterPars.constrainDiboson = False
fitterPars.constrainDibosonShape = False
fitterPars.constrainWpJShape = False
fitterPars.constrainWpJ = True
fitterPars.doEffCorrections = True
fitterPars.muIdEffFiles.push_back("EffTableDir/muonEffsRecoToIso_ScaleFactors.txt")
fitterPars.muHLTEffFiles.push_back("EffTableDir/muonEffsIsoToHLT_data_LP_LWA.txt")
fitterPars.lumiPerEpochMuon.push_back(fitterPars.intLumi)
fitterPars.eleIdEffFiles.push_back("EffTableDir/eleEffsRecoToWP80_ScaleFactors.txt")
fitterPars.eleRecoEffFiles.push_back("EffTableDir/eleEffsSCToReco_ScaleFactors.txt")
fitterPars.eleHLTEffFiles.push_back("EffTableDir/eleEffsSingleElectron.txt")
fitterPars.eleJ30EffFiles.push_back("EffTableDir/FullyEfficient.txt")
fitterPars.eleJ25NoJ30EffFiles.push_back("EffTableDir/FullyEfficient_Jet2NoJet1.txt")
fitterPars.eleMHTEffFiles.push_back("EffTableDir/FullyEfficient_MHT.txt")
fitterPars.eleWMtEffFiles.push_back("EffTableDir/WMt50TriggerEfficiency.txt")
fitterPars.lumiPerEpochElectron.push_back(fitterPars.intLumi)
fitterPars.useExternalMorphingPars = False
## jetCut = '(jetPFCor_Pt[%i] < 2) ' % (Nj)
## for j in range(0, Nj):
## jetCut += '&& (JetPFCor_Pt[%i] > 30) ' % (j)
## jetCut = '(ggdevt == %i)' % (Nj)
## if (Nj < 2):
## jetCut = '((ggdevt == 2) || (ggdevt == 3))'
btagCut = 1.74
fitterPars.cuts = '(sqrt(JetPFCor_Pt[0]**2+JetPFCor_Pt[1]**2+2*JetPFCor_Pt[0]*JetPFCor_Pt[1]*cos(JetPFCor_Phi[0]-JetPFCor_Phi[1]))>20.) ' + \
'&& (abs(JetPFCor_Eta[0]-JetPFCor_Eta[1])<1.5) ' + \
'&& (abs(JetPFCor_dphiMET[0])>0.4) '
fitterPars.cuts += '&& (JetPFCor_Pt[2]<30) '
#fitterPars.cuts += '&& (fit_status==0) '
if btag:
fitterPars.useWbbPDF = False
fitterPars.cuts += '&& ( (JetPFCor_bDiscriminator[0]>%.2f) || (JetPFCor_bDiscriminator[1]>%.2f) ) ' % (btagCut,btagCut)
fitterPars.cuts += '&& (JetPFCor_Pt[1]>35) '
else:
fitterPars.useWbbPDF = False
fitterPars.cuts += '&& (JetPFCor_bDiscriminator[0]<%.2f) && (JetPFCor_bDiscriminator[1]<%.2f) ' % (btagCut,btagCut)
fitterPars.cuts += '&& (JetPFCor_Pt[1]>35.) '
return fitterPars
def theConfig(Nj, mcdir = '', initFile = ''):
fitterPars = RooWjjFitterParams()
fitterPars.smoothingOrder = 1
fitterPars.MCDirectory = "/uscms_data/d2/kalanand/WjjTrees/Full2011DataFall11MC/ReducedTree_VBF/RD_"
## fitterPars.WpJDirectory = str(fitterPars.MCDirectory).replace("PAT", "AOD")
fitterPars.WpJDirectory = fitterPars.MCDirectory
if (len(mcdir) > 0):
fitterPars.WpJDirectory = mcdir
fitterPars.toyWpJ = False
fitterPars.QCDDirectory = fitterPars.MCDirectory[:-3]
fitterPars.initParamsFile = initFile
fitterPars.constraintParamsFile = initFile;
fitterPars.DataDirectory = fitterPars.MCDirectory
fitterPars.muonData = 'WmunuJets_DataAll_GoldenJSON_4p7invfb.root'
fitterPars.includeMuons = True
fitterPars.electronData = 'WenuJets_DataAllSingleElectronTrigger_GoldenJSON_4p7invfb.root'
fitterPars.includeElectrons = True
fitterPars.overwriteDibosonPDFWithVBFWW = True
fitterPars.NewPhysicsDirectory = fitterPars.MCDirectory
fitterPars.minMass = 40.
fitterPars.maxMass = 200.
fitterPars.nbins = 30
fitterPars.intLumi = 5020.
fitterPars.binEdges.push_back(fitterPars.minMass)
#fitterPars.binEdges.push_back(40.)
fitterPars.binEdges.push_back(45.)
fitterPars.binEdges.push_back(50.)
fitterPars.binEdges.push_back(55.)
fitterPars.binEdges.push_back(60.)
fitterPars.binEdges.push_back(65.)
fitterPars.binEdges.push_back(70.)
fitterPars.binEdges.push_back(75.)
fitterPars.binEdges.push_back(80.)
fitterPars.binEdges.push_back(85.)
fitterPars.binEdges.push_back(95.)
fitterPars.binEdges.push_back(105.)
fitterPars.binEdges.push_back(115.)
fitterPars.binEdges.push_back(125.)
fitterPars.binEdges.push_back(135.)
fitterPars.binEdges.push_back(145.)
fitterPars.binEdges.push_back(160.)
fitterPars.binEdges.push_back(175.)
fitterPars.binEdges.push_back(200.)
## binEdge = fitterPars.minMass
## print binEdge,' ',
## while (binEdge <= fitterPars.maxMass):
## binEdge += round(0.14*binEdge)
## fitterPars.binEdges.push_back(binEdge)
## print binEdge,' ',
## fitterPars.binEdges[fitterPars.binEdges.size()-1] = fitterPars.maxMass
print "mass range:",fitterPars.minMass,'-',fitterPars.maxMass
fitterPars.truncRange = True
fitterPars.blind = False
fitterPars.minTrunc = 65.
binMin = fitterPars.minMass
for binEdge in fitterPars.binEdges:
if binEdge > fitterPars.minTrunc:
fitterPars.minTrunc = binMin
break
binMin = binEdge
fitterPars.maxTrunc = 94.
for binEdge in fitterPars.binEdges:
if binEdge > fitterPars.maxTrunc:
fitterPars.maxTrunc = binEdge
break
fitterPars.minFit = fitterPars.minMass
binMin = fitterPars.minMass
for binEdge in fitterPars.binEdges:
if binEdge > fitterPars.minFit:
fitterPars.minFit = binMin
break
binMin = binEdge
fitterPars.maxFit = fitterPars.maxMass
for binEdge in fitterPars.binEdges:
if binEdge > fitterPars.maxFit:
fitterPars.maxFit = binEdge
break
fitterPars.minWmass = fitterPars.minTrunc
fitterPars.maxWmass = fitterPars.maxTrunc
print 'truncate:',fitterPars.truncRange,':', fitterPars.minTrunc,'-',\
fitterPars.maxTrunc
print 'fit range:',fitterPars.minFit,'-',fitterPars.maxFit
fitterPars.njets = Nj
fitterPars.constrainDiboson = True
fitterPars.doEffCorrections = True
fitterPars.muIdEffFiles.push_back("EffTableDir/muonEffsRecoToIso_ScaleFactors.txt")
fitterPars.muHLTEffFiles.push_back("EffTableDir/muonEffsIsoToHLT_data_LP_LWA.txt")
fitterPars.lumiPerEpochMuon.push_back(fitterPars.intLumi)
fitterPars.eleIdEffFiles.push_back("EffTableDir/eleEffsRecoToWP80_ScaleFactors.txt")
fitterPars.eleRecoEffFiles.push_back("EffTableDir/eleEffsSCToReco_ScaleFactors.txt")
fitterPars.eleHLTEffFiles.push_back("EffTableDir/eleEffsSingleElectron.txt")
fitterPars.eleJ30EffFiles.push_back("EffTableDir/FullyEfficient.txt")
fitterPars.eleJ25NoJ30EffFiles.push_back("EffTableDir/FullyEfficient_Jet2NoJet1.txt")
fitterPars.eleMHTEffFiles.push_back("EffTableDir/FullyEfficient_MHT.txt")
fitterPars.eleWMtEffFiles.push_back("EffTableDir/WMt50TriggerEfficiency.txt")
fitterPars.lumiPerEpochElectron.push_back(fitterPars.intLumi)
fitterPars.constrainDiboson = True
fitterPars.useExternalMorphingPars = False
fitterPars.var = "vbf_wjj_m"
fitterPars.cuts = '(vbf_event==1) ' + \
'&& (JetPFCor_bDiscriminator[vbf_waj_id] < 1.74) ' + \
'&& (JetPFCor_bDiscriminator[vbf_wbj_id] < 1.74) ' + \
'&& (vbf_lvjj_m > %0.1f) && (vbf_lvjj_m < %0.1f) ' % \
(minMlvjj, maxMlvjj)
return fitterPars
def theConfig(Nj, mcdir='', initFile='', btag=False, toydataFile=""):
fitterPars = RooWjjFitterParams()
fitterPars.smoothingOrder = 1
# fitterPars.MCDirectory = "/uscms_data/d2/kalanand/WjjTrees/Full2011DataFall11MC/ReducedTree_PAT/RD_"
fitterPars.MCDirectory = "/uscms_data/d2/kalanand/WjjTrees/Full2011DataFall11MC/OriginalTree_PAT/"
fitterPars.WpJDirectory = fitterPars.MCDirectory
if (len(mcdir) > 0):
fitterPars.WpJDirectory = mcdir
fitterPars.toyWpJ = False
##fitterPars.QCDDirectory = fitterPars.MCDirectory[:-3]
fitterPars.QCDDirectory = "/uscms_data/d2/kalanand/WjjTrees/Full2011DataFall11MC/ReducedTree_PAT/"
fitterPars.initParamsFile = initFile
fitterPars.constraintParamsFile = initFile
fitterPars.DataDirectory = fitterPars.MCDirectory
fitterPars.muonData = 'WmunuJets_DataAll_GoldenJSON_4p7invfb.root'
fitterPars.includeMuons = True
fitterPars.electronData = 'WenuJets_DataAllSingleElectronTrigger_GoldenJSON_4p7invfb.root'
fitterPars.includeElectrons = True
fitterPars.NewPhysicsDirectory = fitterPars.MCDirectory
fitterPars.toydataFile = toydataFile
fitterPars.minMass = 40.
fitterPars.maxMass = 150.
fitterPars.nbins = 32
fitterPars.intLumi = 5020.
binEdge = fitterPars.minMass
fitterPars.binEdges.push_back(binEdge)
print binEdge, ' ',
while (binEdge < fitterPars.maxMass):
binEdge += round(0.10 * binEdge)
fitterPars.binEdges.push_back(binEdge)
print binEdge, ' ',
## fitterPars.binEdges[fitterPars.binEdges.size()-1] = fitterPars.maxMass
fitterPars.maxMass = fitterPars.binEdges[fitterPars.binEdges.size() - 1]
print "mass range:", fitterPars.minMass, '-', fitterPars.maxMass
fitterPars.truncRange = False
fitterPars.blind = False
fitterPars.minTrunc = 65.
binMin = fitterPars.minMass
for binEdge in fitterPars.binEdges:
if binEdge > fitterPars.minTrunc:
fitterPars.minTrunc = binMin
break
binMin = binEdge
fitterPars.maxTrunc = 94.
for binEdge in fitterPars.binEdges:
if binEdge > fitterPars.maxTrunc:
fitterPars.maxTrunc = binEdge
break
fitterPars.minFit = fitterPars.minMass
binMin = fitterPars.minMass
for binEdge in fitterPars.binEdges:
if binEdge > fitterPars.minFit:
fitterPars.minFit = binMin
break
binMin = binEdge
fitterPars.maxFit = fitterPars.maxMass
for binEdge in fitterPars.binEdges:
if binEdge > fitterPars.maxFit:
fitterPars.maxFit = binEdge
break
fitterPars.minWmass = fitterPars.minTrunc
fitterPars.maxWmass = fitterPars.maxTrunc
print 'truncate:',fitterPars.truncRange,':', fitterPars.minTrunc,'-',\
fitterPars.maxTrunc
print 'fit range:', fitterPars.minFit, '-', fitterPars.maxFit
fitterPars.njets = Nj
fitterPars.constrainDiboson = False
fitterPars.constrainDibosonShape = False
fitterPars.constrainWpJShape = False
fitterPars.constrainWpJ = True
fitterPars.doEffCorrections = True
fitterPars.muIdEffFiles.push_back(
"EffTableDir/muonEffsRecoToIso_ScaleFactors.txt")
fitterPars.muHLTEffFiles.push_back(
"EffTableDir/muonEffsIsoToHLT_data_LP_LWA.txt")
fitterPars.lumiPerEpochMuon.push_back(fitterPars.intLumi)
fitterPars.eleIdEffFiles.push_back(
"EffTableDir/eleEffsRecoToWP80_ScaleFactors.txt")
fitterPars.eleRecoEffFiles.push_back(
"EffTableDir/eleEffsSCToReco_ScaleFactors.txt")
fitterPars.eleHLTEffFiles.push_back(
"EffTableDir/eleEffsSingleElectron.txt")
fitterPars.eleJ30EffFiles.push_back("EffTableDir/FullyEfficient.txt")
fitterPars.eleJ25NoJ30EffFiles.push_back(
"EffTableDir/FullyEfficient_Jet2NoJet1.txt")
fitterPars.eleMHTEffFiles.push_back("EffTableDir/FullyEfficient_MHT.txt")
fitterPars.eleWMtEffFiles.push_back(
"EffTableDir/WMt50TriggerEfficiency.txt")
fitterPars.lumiPerEpochElectron.push_back(fitterPars.intLumi)
fitterPars.useExternalMorphingPars = False
## jetCut = '(jetPFCor_Pt[%i] < 2) ' % (Nj)
## for j in range(0, Nj):
## jetCut += '&& (JetPFCor_Pt[%i] > 30) ' % (j)
## jetCut = '(ggdevt == %i)' % (Nj)
## if (Nj < 2):
## jetCut = '((ggdevt == 2) || (ggdevt == 3))'
btagCut = 1.74
fitterPars.cuts = '(sqrt(JetPFCor_Pt[0]**2+JetPFCor_Pt[1]**2+2*JetPFCor_Pt[0]*JetPFCor_Pt[1]*cos(JetPFCor_Phi[0]-JetPFCor_Phi[1]))>20.) ' + \
'&& (abs(JetPFCor_Eta[0]-JetPFCor_Eta[1])<1.5) ' + \
'&& (abs(JetPFCor_dphiMET[0])>0.4) '
fitterPars.cuts += '&& (JetPFCor_Pt[2]<30) '
#fitterPars.cuts += '&& (fit_status==0) '
if btag:
fitterPars.useWbbPDF = False
fitterPars.cuts += '&& ( (JetPFCor_bDiscriminator[0]>%.2f) || (JetPFCor_bDiscriminator[1]>%.2f) ) ' % (
btagCut, btagCut)
fitterPars.cuts += '&& (JetPFCor_Pt[1]>35) '
else:
fitterPars.useWbbPDF = False
fitterPars.cuts += '&& (JetPFCor_bDiscriminator[0]<%.2f) && (JetPFCor_bDiscriminator[1]<%.2f) ' % (
btagCut, btagCut)
fitterPars.cuts += '&& (JetPFCor_Pt[1]>35.) '
return fitterPars
def theConfig(Nj, mcdir='', initFile=''):
fitterPars = RooWjjFitterParams()
fitterPars.smoothingOrder = 1
fitterPars.MCDirectory = "/uscms_data/d2/kalanand/WjjTrees/Full2011DataFall11MC/ReducedTree_VBF/RD_"
## fitterPars.WpJDirectory = str(fitterPars.MCDirectory).replace("PAT", "AOD")
fitterPars.WpJDirectory = fitterPars.MCDirectory
if (len(mcdir) > 0):
fitterPars.WpJDirectory = mcdir
fitterPars.toyWpJ = False
fitterPars.QCDDirectory = fitterPars.MCDirectory[:-3]
fitterPars.initParamsFile = initFile
fitterPars.constraintParamsFile = initFile
fitterPars.DataDirectory = fitterPars.MCDirectory
fitterPars.muonData = 'WmunuJets_DataAll_GoldenJSON_4p7invfb.root'
fitterPars.includeMuons = True
fitterPars.electronData = 'WenuJets_DataAllSingleElectronTrigger_GoldenJSON_4p7invfb.root'
fitterPars.includeElectrons = True
fitterPars.overwriteDibosonPDFWithVBFWW = True
fitterPars.NewPhysicsDirectory = fitterPars.MCDirectory
fitterPars.minMass = 40.
fitterPars.maxMass = 200.
fitterPars.nbins = 30
fitterPars.intLumi = 5020.
fitterPars.binEdges.push_back(fitterPars.minMass)
#fitterPars.binEdges.push_back(40.)
fitterPars.binEdges.push_back(45.)
fitterPars.binEdges.push_back(50.)
fitterPars.binEdges.push_back(55.)
fitterPars.binEdges.push_back(60.)
fitterPars.binEdges.push_back(65.)
fitterPars.binEdges.push_back(70.)
fitterPars.binEdges.push_back(75.)
fitterPars.binEdges.push_back(80.)
fitterPars.binEdges.push_back(85.)
fitterPars.binEdges.push_back(95.)
fitterPars.binEdges.push_back(105.)
fitterPars.binEdges.push_back(115.)
fitterPars.binEdges.push_back(125.)
fitterPars.binEdges.push_back(135.)
fitterPars.binEdges.push_back(145.)
fitterPars.binEdges.push_back(160.)
fitterPars.binEdges.push_back(175.)
fitterPars.binEdges.push_back(200.)
## binEdge = fitterPars.minMass
## print binEdge,' ',
## while (binEdge <= fitterPars.maxMass):
## binEdge += round(0.14*binEdge)
## fitterPars.binEdges.push_back(binEdge)
## print binEdge,' ',
## fitterPars.binEdges[fitterPars.binEdges.size()-1] = fitterPars.maxMass
print "mass range:", fitterPars.minMass, '-', fitterPars.maxMass
fitterPars.truncRange = True
fitterPars.blind = False
fitterPars.minTrunc = 65.
binMin = fitterPars.minMass
for binEdge in fitterPars.binEdges:
if binEdge > fitterPars.minTrunc:
fitterPars.minTrunc = binMin
break
binMin = binEdge
fitterPars.maxTrunc = 94.
for binEdge in fitterPars.binEdges:
if binEdge > fitterPars.maxTrunc:
fitterPars.maxTrunc = binEdge
break
fitterPars.minFit = fitterPars.minMass
binMin = fitterPars.minMass
for binEdge in fitterPars.binEdges:
if binEdge > fitterPars.minFit:
fitterPars.minFit = binMin
break
binMin = binEdge
fitterPars.maxFit = fitterPars.maxMass
for binEdge in fitterPars.binEdges:
if binEdge > fitterPars.maxFit:
fitterPars.maxFit = binEdge
break
fitterPars.minWmass = fitterPars.minTrunc
fitterPars.maxWmass = fitterPars.maxTrunc
print 'truncate:',fitterPars.truncRange,':', fitterPars.minTrunc,'-',\
fitterPars.maxTrunc
print 'fit range:', fitterPars.minFit, '-', fitterPars.maxFit
fitterPars.njets = Nj
fitterPars.constrainDiboson = True
fitterPars.doEffCorrections = True
fitterPars.muIdEffFiles.push_back(
"EffTableDir/muonEffsRecoToIso_ScaleFactors.txt")
fitterPars.muHLTEffFiles.push_back(
"EffTableDir/muonEffsIsoToHLT_data_LP_LWA.txt")
fitterPars.lumiPerEpochMuon.push_back(fitterPars.intLumi)
fitterPars.eleIdEffFiles.push_back(
"EffTableDir/eleEffsRecoToWP80_ScaleFactors.txt")
fitterPars.eleRecoEffFiles.push_back(
"EffTableDir/eleEffsSCToReco_ScaleFactors.txt")
fitterPars.eleHLTEffFiles.push_back(
"EffTableDir/eleEffsSingleElectron.txt")
fitterPars.eleJ30EffFiles.push_back("EffTableDir/FullyEfficient.txt")
fitterPars.eleJ25NoJ30EffFiles.push_back(
"EffTableDir/FullyEfficient_Jet2NoJet1.txt")
fitterPars.eleMHTEffFiles.push_back("EffTableDir/FullyEfficient_MHT.txt")
fitterPars.eleWMtEffFiles.push_back(
"EffTableDir/WMt50TriggerEfficiency.txt")
fitterPars.lumiPerEpochElectron.push_back(fitterPars.intLumi)
fitterPars.constrainDiboson = True
fitterPars.useExternalMorphingPars = False
fitterPars.var = "vbf_wjj_m"
fitterPars.cuts = '(vbf_event==1) ' + \
'&& (JetPFCor_bDiscriminator[vbf_waj_id] < 1.74) ' + \
'&& (JetPFCor_bDiscriminator[vbf_wbj_id] < 1.74) ' + \
'&& (vbf_lvjj_m > %0.1f) && (vbf_lvjj_m < %0.1f) ' % \
(minMlvjj, maxMlvjj)
return fitterPars