print ' output parameter file: %s' % parFileOut
print ' output fit result file: %s' % resultFileOut
print ' time acceptance file 2011: %s' % timeAccFile2011
print ' time acceptance file 2012: %s' % timeAccFile2012
print ' angular acceptance type: %s' % args.angAccType
print ' angular acceptance file: %s' % angAccFile
print ' constant angular acceptance parameters: %s' % ( 'true' if constAngAcc else 'false' )
# clear command-line options
import sys
sys.argv = sys.argv[ 0 : 1 ]
# PDF options
from math import pi
from P2VV.Parameterizations.FullPDFs import Bs2Jpsiphi_RunIAnalysis as PdfConfig
pdfConfig = PdfConfig( RunPeriods = args.runPeriod )
pdfConfig['blind'] = { }
pdfConfig['splitParams'] = { }
pdfConfig['lambdaCPParam'] = 'observables_CPVDecay' if args.model == 'polarDep' else 'lambPhi'
pdfConfig['paramKKMass'] = ''
pdfConfig['KKMassBinBounds'] = [ 990., 1050. ]
pdfConfig['CSPValues'] = [ 1. ]
pdfConfig['timeEffParameters'] = dict( Cache = False )
pdfConfig['timeEffHistFiles']['file'] = timeAccFile2011 if args.runPeriod == '2011' else timeAccFile2012
if args.runPeriod == '2011' :
#timeAccHistHLT1UB2011 = 'Bs_HltPropertimeAcceptance_Data_2011_40bins_Hlt1DiMuon_Hlt2DiMuonDetached_Reweighted'
timeAccHistHLT1UB = 'Bs_HltPropertimeAcceptance_Data_2011_40bins_Hlt1DiMuon_Hlt2DiMuonDetached'
pdfConfig['timeEffHistFiles']['hlt1UB'] = timeAccHistHLT1UB
print ' output parameter file: %s' % parFileOut
print ' output fit result file: %s' % resultFileOut
print ' time acceptance file 2011: %s' % timeAccFile2011
print ' time acceptance file 2012: %s' % timeAccFile2012
print ' angular acceptance type: %s' % args.angAccType
print ' angular acceptance file: %s' % angAccFile
print ' constant angular acceptance parameters: %s' % ( 'true' if constAngAcc else 'false' )
# clear command-line options
import sys
sys.argv = sys.argv[ 0 : 1 ]
# PDF options
from math import pi
from P2VV.Parameterizations.FullPDFs import Bs2Jpsiphi_RunIAnalysis as PdfConfig
pdfConfig = PdfConfig()
if parFileIn :
# read parameters from file
pdfConfig.readParametersFromFile( filePath = parFileIn )
if not blind :
pdfConfig['blind'] = { }
#pdfConfig['KKMassBinBounds'] = [ 1020. - 12., 1020. - 4., 1020., 1020. + 4., 1020. + 12. ]
#pdfConfig['CSPValues'] = [ 0.9022, 0.8619, 0.8875, 0.9360 ]
#pdfConfig['KKMassBinBounds'] = [ 1020. - 12., 1020. + 12. ]
#pdfConfig['CSPValues'] = [ 0.498 ]
pdfConfig['tagCatsType'] = args.tagCatsType
pdfConfig['contEstWTag'] = contEstWTag
# data
dataPath = '/project/bfys/jleerdam/data/Bs2Jpsiphi/'
myPath = '/project/bfys/vsyropou/data/Bs2JpsiPhi/'
dataSetFile = options.FitData if options.FitData else \
dataPath + 'angEff/P2VVDataSets20112012Reco14_I2Mass_6KKMassBins_2TagCats_HLT2B.root'
dataSetName = 'JpsiKK_sigSWeight'
# read / write fited parameters from file
parFileIn = options.ParFileIn if options.ParFileIn \
else myPath + 'nominalFitResults/corrAngAccDEC/20112012Reco14DataFitValues_6KKMassBins.root'
parFileOut = options.ParFileOut if options.ParFileOut else '20112012Reco14DataFitValues_6KKMassBins.par'
# PDF configuration
from P2VV.Parameterizations.FullPDFs import Bs2Jpsiphi_RunIAnalysis as PdfConfig
pdfConfig = PdfConfig()
# swich blidn on/off
if 'True' in options.unblind: pdfConfig['blind'] = {}
# cpv parametrization
pdfConfig['lambdaCPParam'] = 'lambPhi'
# acceptances
# decay time
pdfConfig['timeEffType'] = 'paper2012'
pdfConfig['timeEffHistFiles'].getSettings( [ ( 'runPeriod', 'p2011' ) ] )['file']\
= dataPath + 'Reco14/timeAcceptanceFit_2011.root' # Bs_HltPropertimeAcceptance_Data_2011_40bins_TOS.root'
pdfConfig['timeEffHistFiles'].getSettings( [ ( 'runPeriod', 'p2012' ) ] )['file']\
= dataPath + 'Reco14/timeAcceptanceFit_2012.root' # Bs_HltPropertimeAcceptance_Data_2012_40bins_TOS.root'
# decay angles
parser.add_argument('plotComp')
args = parser.parse_args()
plotVar = args.plotVar
plotComp = args.plotComp
assert plotVar in [ 'timeLin', 'timeLog', 'ctk', 'ctl', 'phi' ]
assert plotComp in [ 'total', 'phi', 'even', 'odd', 'S' ]
print 'job parameters:'
print ' plot variable: %s' % plotVar
print ' plot component: %s' % plotComp
from math import pi
from P2VV.Parameterizations.FullPDFs import Bs2Jpsiphi_RunIAnalysis as PdfConfig
pdfConfig = PdfConfig()
# job parameters
#parFileIn = '/project/bfys/jleerdam/softDevel/P2VV2/test/plots/Reco14/parVals.par'
parFileIn = '/project/bfys/jleerdam/softDevel/P2VV2/test/plots/Reco14/LHCb/parVals.par'
dataPath = '/project/bfys/jleerdam/data/Bs2Jpsiphi/Reco14/'
dataSetName = 'JpsiKK_sigSWeight'
#dataSetFile = dataPath + 'P2VVDataSets20112012Reco14_I2Mass_6KKMassBins_2TagCats_HLT2B_20140309.root'
dataSetFile = dataPath + 'P2VVDataSets20112012Reco14_I2Mass_6KKMassBins_2TagCats_HLT2B_20140822_fromNTuple_v2_corrWeights.root'
#plotsFilePath = '/project/bfys/jleerdam/softDevel/P2VV2/test/plots/Reco14/projPlot_%s_%s.ps' % ( plotVar, plotComp )
#ROOTFilePath = '/project/bfys/jleerdam/softDevel/P2VV2/test/plots/Reco14/projPlot_%s_%s.root' % ( plotVar, plotComp )
plotsFilePath = '/project/bfys/jleerdam/softDevel/P2VV2/test/plots/Reco14/LHCb/projPlot_%s_%s.ps' % ( plotVar, plotComp )
ROOTFilePath = '/project/bfys/jleerdam/softDevel/P2VV2/test/plots/Reco14/LHCb/projPlot_%s_%s.root' % ( plotVar, plotComp )
# PDF options
#pdfConfig['timeEffHistFiles'].getSettings( [ ( 'runPeriod', 'p2011' ) ] )['file'] = dataPath + 'timeAcceptanceFit_2011.root'
args = parser.parse_args()
jobID = args.jobID
startIndex = args.startIndex
testParNames = [ args.parName ]
fitParVals = [ [ val ] for val in args.parVals ]
print 'job parameters:'
print ' ID:', jobID
print ' start index:', startIndex
print ' index width:', indexWidth
print ' parameter names:', testParNames
print ' parameter values:', fitParVals
from math import pi
from P2VV.Parameterizations.FullPDFs import Bs2Jpsiphi_RunIAnalysis as PdfConfig
pdfConfig = PdfConfig()
# job parameters
parFileIn = '20112012Reco14DataPlotValues_6KKMassBins_CPVDecay.par'
dataPath = '/project/bfys/jleerdam/data/Bs2Jpsiphi/Reco14/'
dataSetName = 'JpsiKK_sigSWeight'
dataSetFile = dataPath + 'P2VVDataSets20112012Reco14_I2DiegoMass_6KKMassBins_2TagCats.root'
# fit options
fitOpts = dict( NumCPU = 1
, Optimize = 2
, Minimizer = 'Minuit2'
)
pdfConfig['fitOptions'] = fitOpts
corrSFitErrCats = [ 'runPeriod', 'KKMassCat' ]
randomParVals = ( ) # ( 1., 12345 )
