struefit = ScaleFitter( name = 'strue_mc_NominalFitRange68', ## name = 'strue_mc_FitRangePositive', title = 'strue-Fit, Powheg S4', labels = ['Powheg S4 Summer11 MC'], source = _chains['z'], xName = 's', xTitle = 's_{true} = E^{#gamma}_{reco}/E^{#gamma}_{gen} - 1', xExpression = '100 * (phoE/phoGenE - 1)', cuts = ['isFSR', 'phoGenE > 0'], xRange = (-50, 100), xUnit = '%', nBins = 150, pdf = 'gauss', graphicsExtensions = ['png'], massWindowScale = 1.5, massWindow = (87.2, 95.2), fitScale = 1.2, fitRange = (-50,50), doAutoBinning = True, binContentMax = 200, binContentMin = 35, canvasStyle = 'landscape', doAutoXRange = True, doAutoXRangeZoom = True, doAutoFitRange = True, xRangeSigmaLevel = 5, xRangeSigmaLevelZoom = 2, fitRangeMode = 'Fraction', fitRangeSigmaLevel = 2.0, fitRangeNumberOfEntries = 3000, fitRangeFraction = 0.68, paramLayout = (0.57, 0.92, 0.92), useCustomChi2Calculator = True, )
## Default fit of s = Ereco / Ekin - 1 srecofit = ScaleFitter( name = 'sreco_mc', title = 'sreco-Fit, Powheg S4', labels = ['Powheg S4'], cuts = [], source = _chains['z'], xName = 's', xTitle = 's_{reco} = E^{#gamma}_{reco}/E^{kin}_{reco} - 1', xExpression = '100 * (1/kRatio - 1)', xRange = (-50, 100), xUnit = '%', nBins = 150, pdf = 'cbShape', graphicsExtensions = ['png'], massWindowScale = 1.5, massWindow = (87.2, 95.2), fitScale = 1.2, fitRange = (-50,100), doAutoBinning = True, binContentMax = 200, binContentMin = 35, canvasStyle = 'landscape', doAutoXRange = False, doAutoXRangeZoom = True, xRangeSigmaLevelZoom = 5, paramLayout = (.45, 0.75, 0.5), useCustomChi2Calculator = True, ) # srecofit.name += ('_FitRange%d' % int(sreco_fitrange))
srecofit = ScaleFitter( name = 'sreco_mc', title = 'sreco-Fit, Powheg S4', labels = ['Powheg S4'], cuts = [], source = _chains['z'], xName = 's', xTitle = 's_{reco} = E^{#gamma}_{reco}/E^{kin}_{reco} - 1', xExpression = '100 * (1/kRatio - 1)', xRange = (-50, 100), xUnit = '%', nBins = 150, pdf = 'cbShape', graphicsExtensions = ['png'], massWindowScale = 1.5, massWindow = (87.2, 95.2), fitScale = 1.2, fitRange = (-50,100), doAutoBinning = True, binContentMax = 200, binContentMin = 35, canvasStyle = 'landscape', doAutoXRange = True, doAutoXRangeZoom = True, doAutoFitRange = True, xRangeFraction = 1, xRangeMode = 'Fraction', xRangeSigmaLevelZoom = 3, paramLayout = (.45, 0.75, 0.5), useCustomChi2Calculator = True, )
## 715/pb for Vg Summer conferences # _chains = esChains.getChains('v7') ## 2/fb of LP11 dataset _chains = esChains.getChains("v11") ## Default fit of s = Ereco / Ekin - 1 srecofit = ScaleFitter( name="sreco_mc", title="sreco-Fit, Powheg S4", labels=["Powheg S4"], cuts=[], source=_chains["z"], xName="s", xTitle="s_{reco} = E^{#gamma}_{reco}/E^{kin}_{reco} - 1", xExpression="100 * (1/kRatio - 1)", xRange=(-50, 100), xUnit="%", nBins=150, # fitRange = (-50, 100), pdf="gauss", graphicsExtensions=["png"], # graphicsExtensions = [], massWindowScale=1.5, massWindow=(87.2, 95.2), fitScale=1.2, ) ## Default fit of strue = Ereco / Egen - 1 struefit = srecofit.clone( name="strue_mc", title="strue-Fit, Powheg S4", xTitle="s_{true} = E^{#gamma}_{reco}/E^{#gamma}_{gen} - 1",
## end of class BremBin ## Default fit of sraw = Eraw / Egen - 1 srawfit = ScaleFitter( name='sraw_mc', title='sraw-Fit, Powheg S4', labels=[ 'Flat-pt #gamma gun', '#eta_{SC} #in [1.16,1.44]', '#phi cracks removed', 'sech fit' ], cuts=['!isEBEtaGap & !isEBPhiGap & 1.16 < abs(scEta) & abs(scEta) < 1.44'], source=_chains['g93p01'], xName='s', xTitle='s_{true} = E^{SC}_{raw}/E^{#gamma}_{gen} - 1', xExpression='100 * (scRawE/genE - 1)', xRange=(-10, 5), xUnit='%', nBins=120, fitRange=(-10, 5), pdf='sech', # graphicsExtensions = ['png'], graphicsExtensions=[], paramLayout=(.25, 0.55, 0.92), # x1, x2, y1 labelsLayout=(.25, 0.6), # x1, y1 ) ## ---------------------------------------------------------------------------- ## Customize below _fits = [] for geant in 'g93p01 g94p02 g94cms'.split():
return daughter ## end of SubdetFitters if __name__ == '__main__': _chains = esChains.getChains('v10') test_fitter = ScaleFitter( name = 's', title = 's-Fit', source = '_chains["z"]', expression = '100 * (1/kRatio - 1)', cuts = ['mmMass < 80'], # 'phoR9 < 0.94', '30 < phoPt', 'phoPt < 9999'], labels = [#'R_{9}^{#gamma} < 0.94', #'E_{T}^{#gamma} #in [X,Y] GeV', 'Powheg S4', 'Test Model'], xRange = (-20, 40), nBins = 120, fitRange = (-100, 100), pdf = 'lognormal', # graphicsExtensions = ['png', 'eps'], graphicsExtensions = [], massWindowScale = 1.5, fitScale = 2.0 ) print test_fitter.pydump(), ',' for fit in SubdetR9Fitters(test_fitter): print fit.pydump() import user
## Default fit of s = Ereco / Ekin - 1 srecofit = ScaleFitter( name="sreco_mc", title="sreco-Fit, Powheg S4", labels=["Powheg S4"], cuts=[], source=_chains["z"], xName="s", xTitle="s_{reco} = E^{#gamma}_{reco}/E^{kin}_{reco} - 1", xExpression="100 * (1/kRatio - 1)", xRange=(-50, 100), xUnit="%", nBins=150, pdf="cbShape", graphicsExtensions=["png"], massWindowScale=1.5, massWindow=(87.2, 95.2), fitScale=1.2, fitRange=(-50, 100), doAutoBinning=True, binContentMax=200, binContentMin=35, canvasStyle="landscape", doAutoXRange=False, doAutoXRangeZoom=True, xRangeSigmaLevelZoom=5, paramLayout=(0.45, 0.75, 0.5), useCustomChi2Calculator=True, ) ## Default fit of strue = Ereco / Egen - 1
def __str__(self): return self.__class__.__name__ + '(%g, %g)' % self.bin_range ## end of class BremBin ## Default fit of sraw = Eraw / Egen - 1 srawfit = ScaleFitter( name = 'sraw_mc', title = 'sraw-Fit, Powheg S4', labels = ['Flat-pt #gamma gun', '#eta_{SC} #in [1.16,1.44]', '#phi cracks removed', 'sech fit'], cuts = ['!isEBEtaGap & !isEBPhiGap & 1.16 < abs(scEta) & abs(scEta) < 1.44'], source = _chains['g93p01'], xName = 's', xTitle = 's_{true} = E^{SC}_{raw}/E^{#gamma}_{gen} - 1', xExpression = '100 * (scRawE/genE - 1)', xRange = (-10, 5), xUnit = '%', nBins = 120, fitRange = (-10, 5), pdf = 'sech', # graphicsExtensions = ['png'], graphicsExtensions = [], paramLayout = (.25, 0.55, 0.92), # x1, x2, y1 labelsLayout = (.25, 0.6), # x1, y1 ) ## ---------------------------------------------------------------------------- ## Customize below _fits = [] for geant in 'g93p01 g94p02 g94cms'.split():
sfit = ScaleFitter( name = 'strue_mc_modulep4', title = 'strue-Fit, Powheg S4', labels = ['Flat-pt #gamma gun', '#eta_{SC} #in [1.16,1.44]', '#eta/#phi cracks removed',], cuts = ['!isEBEtaGap & !isEBPhiGap', '1.16 < abs(scEta) & abs(scEta) < 1.44', 'scE > 0', 'genE > 0',], source = _chains['g93p01'], xName = 's', xTitle = 's_{true} = E^{SC}_{corr}/E^{#gamma}_{gen} - 1', xExpression = '100 * (scE/genE - 1)', xRange = (-50, 50), xUnit = '%', nBins = 120, fitRange = (-49, 49), pdf = 'cruijff', graphicsExtensions = ['png'], paramLayout = (.2, 0.52, 0.92), # x1, x2, y1 labelsLayout = (.55, 0.55), # x1, y1 # In order to define chi2 well binContentMin = 10, binContentMax = 100, canvasStyle = 'extended', doAutoXRange = True, doAutoXRangeZoom = True, doAutoFitRange = True, xRangeSigmaLevel = 5, xRangeSigmaLevelZoom = 2, fitRangeSigmaLevel = 1.5, useCustomChi2Calculator = False, )
## end of SubdetFitters if __name__ == '__main__': _chains = esChains.getChains('v10') test_fitter = ScaleFitter( name='s', title='s-Fit', source='_chains["z"]', expression='100 * (1/kRatio - 1)', cuts=['mmMass < 80'], # 'phoR9 < 0.94', '30 < phoPt', 'phoPt < 9999'], labels=[ #'R_{9}^{#gamma} < 0.94', #'E_{T}^{#gamma} #in [X,Y] GeV', 'Powheg S4', 'Test Model' ], xRange=(-20, 40), nBins=120, fitRange=(-100, 100), pdf='lognormal', # graphicsExtensions = ['png', 'eps'], graphicsExtensions=[], massWindowScale=1.5, fitScale=2.0) print test_fitter.pydump(), ',' for fit in SubdetR9Fitters(test_fitter): print fit.pydump()