# 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() extractor1.ExtractAllValues() ImE0p.BuildGraph("energy")
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 Qsquared.SetBinValueSize(5.0,1.0) beta.SetBinValueSize(6.0,0.1) aDataPoint.SetValue(0.008) aDataPoint.GetSystError().SetError(0.004)
# 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() A12.BuildGraph("Qsquared") A32 = experiment.GetMeasurement(str('Ap(3/2)')+meas)
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 RExtractor.Initialize() RExtractor.ExtractAllValues() R.BuildGraph()