meas = str(arg)
elif opt in ("-t", "--title"):
mtitle = str(arg)
# else:
# file = 'analyses/MAID2007/data/mult1p/Maid2007_E0p_1p.dat'
# elif opt in ('-v', '--verbose'):
# verbose = True
experiment = manager.GetExperiment("MAID2007")
if not experiment:
experiment = NucDBExperiment("MAID2007", "MAID2007")
# proton Im and Re parts of E0(proton,I=1/2)
ImE0p = experiment.GetMeasurement(str("Im_") + meas)
if not ImE0p:
ImE0p = NucDBMeasurement(str("Im_") + meas, str("Im(") + meas + str(")"))
experiment.AddMeasurement(ImE0p)
ImE0p.ClearDataPoints()
ImE0p.fColor = 8
extractor1 = MAID2007Extractor()
extractor1.iValueRow = 3
extractor1.iStatErr = 4
extractor1.SetMeasurement(ImE0p)
# extractor1.SetInputFile("analyses/MAID2007/data/mult1p/Maid2007_E0p_1p.dat")
extractor1.SetInputFile(file)
extractor1.linestoskip = 4
energy = NucDBBinnedVariable("energy", "E_{#gamma}")
extractor1.fCurrentDataPoint.AddBinVariable(energy)
energyunit = NucDBEnergyMeV()
energy.SetUnit(energyunit)
extractor1.Initialize()
manager = NucDBManager.GetManager(1)
paper = manager.GetPaper("LATTICE2001")
if not paper :
paper = NucDBPaper("LATTICE2001","LATTICE2001")
# d2p - From Phys.Rev.D 63 074506, page 8
d2Calc = paper.GetCalculation('d2')
if not d2Calc :
d2Calc = NucDBCalculation('d2','d_{2}')
paper.AddCalculation(d2Calc)
d2p = d2Calc.GetMeasurement('d2p')
if not d2p :
d2p = NucDBMeasurement('d2p','d_{2}^{p}')
d2Calc.AddMeasurement(d2p)
d2p.ClearDataPoints()
d2p.fColor=4
aDataPoint = NucDBDataPoint()
# Q^2 = \mu^2 = 5.0 GeV^2 From table 5
Qsquared = NucDBBinnedVariable("Qsquared","Q^{2}")
aDataPoint.AddBinVariable(Qsquared)
beta = NucDBBinnedVariable("beta","#beta")
aDataPoint.AddBinVariable(beta)
# 1
meas = str(arg)
elif opt in ('-t','--title'):
mtitle = str(arg)
# else:
# file = 'analyses/MAID2007/data/mult1p/Maid2007_E0p_1p.dat'
# elif opt in ('-v', '--verbose'):
# verbose = True
experiment = manager.GetExperiment("MAID2007")
if not experiment :
experiment = NucDBExperiment("MAID2007","MAID2007")
#proton Im and Re parts of E0(proton,I=1/2)
A12 = experiment.GetMeasurement(str('Ap(1/2)')+meas)
if not A12 :
A12 = NucDBMeasurement(str('Ap(1/2)')+meas,str('Ap(1/2)')+meas)
experiment.AddMeasurement(A12)
A12.ClearDataPoints()
A12.fColor=1
extractor1 = MAID2007Extractor()
extractor1.iValueRow=1
extractor1.iStatErr=2
extractor1.SetMeasurement(A12)
extractor1.SetInputFile(file)
extractor1.linestoskip=6
Qsq = NucDBBinnedVariable("Qsquared","Q^{2}")
extractor1.fCurrentDataPoint.AddBinVariable(Qsq)
#energyunit = NucDBEnergyMeV()
#energy.SetUnit(energyunit)
extractor1.Initialize()
extractor1.ExtractAllValues()
self.fCurrentDataPoint.fValue=float(values[self.iValueRow])
self.fCurrentDataPoint.fStatisticalError.SetError(float(values[self.istatErr]))
self.fCurrentDataPoint.fSystematicError.SetError(float(values[self.isysErr]))
self.fCurrentDataPoint.CalculateTotalError()
self.fCurrentDataPoint.Print()
# get the NucDB manager
manager = NucDBManager.GetManager()
# get/create the experiment
experiment = manager.GetExperiment("SLAC_E143")
if not experiment :
experiment = NucDBExperiment("SLAC_E143","SLAC_E143")
# get/create the measurement
R = experiment.GetMeasurement("R") # get measurement if it already exists
if not R :
R = NucDBMeasurement("R","R")
experiment.AddMeasurement(R)
R.ClearDataPoints() # in case datapoints already exist
R.fColor=2 # default color used for plotting
# configure and run the extractor
RExtractor = SLACE143ExtractorR()
RExtractor.SetMeasurement(R)
RExtractor.fCurrentDataPoint.fName = "R"
Xbjorken = NucDBBinnedVariable("x","x")
Qsq = NucDBBinnedVariable("Qsquared","Q^{2}")
Epsilon = NucDBBinnedVariable("epsilon","#epsilon")
RExtractor.fCurrentDataPoint.fBinnedVariables.Add(Xbjorken)
RExtractor.fCurrentDataPoint.fBinnedVariables.Add(Qsq)
RExtractor.fCurrentDataPoint.fBinnedVariables.Add(Epsilon)
RExtractor.SetInputFile("experiments/SLAC-E143/R.dat") # relative to the NucDB directory
RExtractor.linestoskip=15 # data starts on line 15
