def doF0Fit(inFileName, sectors, startBin, stopBin, tBins, sectorRangeMap = {}, referenceWave = "", deg = 0, writeResultToFile = None):
fixedAMPD_sub = pc.fixedParameterization("/nfs/freenas/tuph/e18/project/compass/analysis/fkrinner/fkrinner/trunk/massDependentFit/scripts/anything/zeroModes/bwAmplitudes_noBF/amp_1pp0pSigmaPiP")
f0Mass = ptc.parameter(mF0,"f0Mass")
f0g1 = ptc.parameter(g1, "g1")
f0g2 = ptc.parameter(g2, "g2")
# f0g1.lock = True
# f0g2.lock = True
f0 = ptc.flatte([f0Mass, f0g1, f0g2], mPi, mPi, mPi, mK, mK, 0, 1, False)
polyPars = []
for i in range(deg):
polyPars.append(ptc.parameter(0.1, "reC"+str(i+1)))
polyPars.append(ptc.parameter(0.1, "imC"+str(i+1)))
poly = ptc.complexPolynomial(deg, polyPars)
fitF0 = amplitudeAnalysis(inFileName, sectors, {"1++0+[pi,pi]0++PiP":[f0,fixedAMPD_sub]}, startBin, stopBin, tBins, sectorRangeMap = sectorRangeMap)
fitF0.loadData(referenceWave = referenceWave)
fitF0.finishModelSetup()
# fitF0.unifyComa()
fitF0.fitShapeParameters()
fitF0.calculateNonShapeParameters()
fitF0.mode = AMPL
if writeResultToFile:
with open(writeResultToFile, 'a') as outFile:
if len(tBins) > 1:
raise ValueError("More than one t' bin not supported")
resultString = str(tBins[0])+ " 666. " + str(f0Mass.value) + ' ' + str(f0Mass.error) + ' ' + str(f0g1.value) + ' ' + str(f0g1.error) + ' ' + str(f0g2.value) + ' ' + str(f0g2.error)
for i in range(deg):
resultString += ' ' + str(polyPars[i].value) + ' ' + str(polyPars[i].error)
resultString += "\n"
outFile.write(resultString)
return fitF0
def doFitRho(inFileName, sectors, startBin, stopBin, tBins, sectorRangeMap = {}, referenceWave = "", writeResultToFile = None, acv = None, polDeg = None):
rhoMass = ptc.parameter( mRho, "rhoMass" )
rhoWidth = ptc.parameter( Grho , "rhoWidth")
rho = ptc.relativisticBreitWigner([rhoMass,rhoWidth], mPi, mPi, mPi, 1, 1, False)
model = [rho]
if polDeg is not None:
params = []
sVal = 1.
for i in range(polDeg):
params.append(ptc.parameter( sVal, "c"+str(i+1)+"_re" ))
sVal = 0. # Make the start value to constant 1
params.append(ptc.parameter( 0., "c"+str(i+1)+"_im" ))
polynome = ptc.complexPolynomial(polDeg,params)
model.append(polynome)
fitRho = amplitudeAnalysis(inFileName, sectors, {"1-+1+[pi,pi]1--PiP":model}, startBin, stopBin, tBins, sectorRangeMap = sectorRangeMap)
fitRho.loadData(loadIntegrals = True, referenceWave = referenceWave)
fitRho.finishModelSetup()
fitRho.fitShapeParameters()
fitRho.calculateNonShapeParameters()
fitRho.mode = AMPL
# fitRho.removeGlobalPhaseFromComa()
if writeResultToFile:
with open(writeResultToFile, 'a') as outFile:
if len(tBins) > 1:
raise ValueError("More than one t' bin not supported")
resultString = str(tBins[0])+ " 666. " + str(rhoMass.value) + ' ' + str(rhoMass.error) + ' ' + str(rhoWidth.value) + ' ' + str(rhoWidth.error) + "\n"
outFile.write(resultString)
return fitRho