mu2 = 'mu2Pt,mu2Eta,mu2Phi,0.1056'
mu1gen = 'mu1GenPt,mu1GenEta,mu1GenPhi,0.1056'
mu2gen = 'mu2GenPt,mu2GenEta,mu2GenPhi,0.1056'
phogen = 'phoGenEt,phoGenEta,phoGenPhi,0'
mmMassGen = 'twoBodyMass({mu1}, {mu2})'.format(mu1=mu1gen, mu2=mu2gen)
mmgMassGen = 'threeBodyMass({mu1}, {mu2}, {pho})'.format(mu1=mu1gen,
mu2=mu2gen,
pho=phogen)
kRatioGen = 'kRatio({mmgMass}, {mmMass})'.format(mmgMass=mmgMassGen,
mmMass=mmMassGen)
## ----------------------------------------------------------------------------
## Customize below
#srecofit.applyDefinitions([Model('cbShape')])
srecofit.applyDefinitions([DimuonMassMax(80)])
_fits = []
for subdet_r9_cat in list(subdet_r9_categories)[:]:
for lo, hi in list(BinEdges([10, 12, 15, 20, 25, 30, 100]))[:]:
for model in 'gauss cbShape bifurGauss cruijff'.split():
print '+++++', model
fit = srecofit.clone()
fit.applyDefinitions([subdet_r9_cat, PhoEtBin(lo, hi),
Model(model)])
_fits.append(fit)
maxIterations = 1
fSigma = 1.5
pullEpsilon = 0.1
mu1 = 'mu1Pt,mu1Eta,mu1Phi,0.1056'
mu2 = 'mu2Pt,mu2Eta,mu2Phi,0.1056'
mu1gen = 'mu1GenPt,mu1GenEta,mu1GenPhi,0.1056'
mu2gen = 'mu2GenPt,mu2GenEta,mu2GenPhi,0.1056'
phogen = 'phoGenEt,phoGenEta,phoGenPhi,0'
mmMassGen = 'twoBodyMass({mu1}, {mu2})'.format(mu1=mu1gen, mu2=mu2gen)
mmgMassGen = 'threeBodyMass({mu1}, {mu2}, {pho})'.format(mu1=mu1gen,
mu2=mu2gen,
pho=phogen)
kRatioGen = 'kRatio({mmgMass}, {mmMass})'.format(mmgMass=mmgMassGen,
mmMass=mmMassGen)
## ----------------------------------------------------------------------------
## Customize below
struefit.applyDefinitions([DimuonMassMax(80)])
struefits =[]
# models = ('gauss lognormal bifurGauss cbShape gamma'.split() +
# 'cruijff gsh bifurGsh bw sumGaussGauss'.split())
models = 'bifurGauss'.split()
# models += 'sumGaussGauss sumGauss3 sumCruijffGauss sumBwGauss'.split()
for subdet_r9_cat in subdet_r9_categories:
for lo, hi in BinEdges([10, 12, 15, 20, 25, 30, 100]):
## for lo, hi in BinEdges([10, 12, 15]):
for model in models:
fit = struefit.clone().applyDefinitions([subdet_r9_cat,
PhoEtBin(lo, hi),
Model(model)])
if ('EB_lowR9_PhoEt10-12' in fit.name or
mu2 = 'mu2Pt,mu2Eta,mu2Phi,0.1056'
mu1gen = 'mu1GenPt,mu1GenEta,mu1GenPhi,0.1056'
mu2gen = 'mu2GenPt,mu2GenEta,mu2GenPhi,0.1056'
phogen = 'phoGenEt,phoGenEta,phoGenPhi,0'
mmMassGen = 'twoBodyMass({mu1}, {mu2})'.format(mu1=mu1gen, mu2=mu2gen)
mmgMassGen = 'threeBodyMass({mu1}, {mu2}, {pho})'.format(mu1=mu1gen,
mu2=mu2gen,
pho=phogen)
kRatioGen = 'kRatio({mmgMass}, {mmMass})'.format(mmgMass=mmgMassGen,
mmMass=mmMassGen)
## ----------------------------------------------------------------------------
## Customize below
struefit.applyDefinitions([Model('bifurGauss')])
srecofit.applyDefinitions([Model('cbShape')])
defaultfits = [struefit, srecofit]
for fit in defaultfits:
fit.applyDefinitions([DimuonMassMax(90)])
srecofits, struefits =[], []
for subdet_r9_cat in subdet_r9_categories:
for lo, hi in BinEdges([10, 12, 15, 20, 25, 30, 100]):
## for lo, hi in BinEdges([10, 12, 15]):
srecofit_cat = srecofit.clone().applyDefinitions([subdet_r9_cat,
PhoEtBin(lo, hi)])
struefit_cat = struefit.clone().applyDefinitions([subdet_r9_cat,
PhoEtBin(lo, hi)])
