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 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():
## 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()