07/01/2007
Check the previous fit
dhad fit regular test1 -t d -m 0 —qsub log-2007-06-30 23:14:22 FAILED.
Uncomment and fit again:
RooNumIntConfig bkgd1config(*Bkgd1.getIntegratorConfig()); bkgd1config.setEpsAbs(1e-8); bkgd1config.setEpsRel(1e-8); Bkgd1.setIntegratorConfig(bkgd1config); RooNumIntConfig bkgd2config(*Bkgd2.getIntegratorConfig()); bkgd2config.setEpsAbs(1e-8); bkgd2config.setEpsRel(1e-8); Bkgd2.setIntegratorConfig(bkgd2config);
dhad fit regular test2 -t d -m 0 —qsub log-2007-07-01 09:15:48
BF fitting
dhad bf
- data results pdg
- data results bf ratio pdg
- Fitting result for stat errors RAW format
- Fitting result for stat and syst errors RAW format
- By mode syst error
BF noFSR fitting
dhad bf noFSR
dhad table compare bf noFSR
07/02/2007
Fix the fitting
Message from Peter:
It turns out that the "analytic" integral was never getting called. I think the version of RooFit that was being used must have changed; the expected function signature changed (RooFit now expects getAnalyticalIntegral() to take an extra rangeName parameter) and so the function was no longer being overridden. [Ah, things that strong static typing can't save you from.] The code in my area (/nfs/cor/user/ponyisi/xin/RooFitModels) contains new RooArgusBG.cc and RooArgusBG.rdl files which do the following: 1. Fix the override problem [always returns 1 for normalization] 2. Use Anders's normalization except now with 10000 bins in mbc
Update RooArgusBG.cc to v1.3
Update RooArgusBG.rdl to v1.2
compile gmake -f GNUmakefile.standalone
dhad fit test3 -t d -m 0 —qsub log-2007-07-02 09:54:14 OK.
Comment and fit again:
RooNumIntConfig bkgd1config(*Bkgd1.getIntegratorConfig()); bkgd1config.setEpsAbs(1e-8); bkgd1config.setEpsRel(1e-8); Bkgd1.setIntegratorConfig(bkgd1config); RooNumIntConfig bkgd2config(*Bkgd2.getIntegratorConfig()); bkgd2config.setEpsAbs(1e-8); bkgd2config.setEpsRel(1e-8); Bkgd2.setIntegratorConfig(bkgd2config);
dhad fit test4 -t d -m 0 —qsub log-2007-07-02 10:51:34 OK.
- Submit the previous failed jobs
dhad fit M 3.7713 -t d -m 201 —qsub log-2007-07-02 13:37:06 OK.
dhad fit M 3.7713 -t d -m 204 —qsub log-2007-07-02 13:37:15 OK.
dhad fit M 3.7723 -t d -m 204 —qsub log-2007-07-02 13:37:24 OK.
dhad fit M 3.7723 -t d -m 205 —qsub log-2007-07-02 13:37:33 OK.
dhad fit Gamma 0.0311 -t d -m 1 —qsub log-2007-07-02 13:37:42 OK.
dhad fit Gamma 0.0311 -t d -m 201 —qsub log-2007-07-02 13:37:50 OK.
dhad fit regular -t d -m 0 —qsub log-2007-07-02 13:38:09 OK.
dhad fit Gamma 0.0261 -t d -m 1 —qsub log-2007-07-02 13:38:17 OK.
dhad fit Gamma 0.0261 -t d -m 205 —qsub log-2007-07-02 13:38:24 OK.
dhad fit R 16.3 -t d -m 1 —qsub log-2007-07-02 13:38:31 OK.
dhad fit R 16.3 -t d -m 201 —qsub log-2007-07-02 13:38:39 OK.
dhad fit R 16.3 -t d -m 204 —qsub log-2007-07-02 13:38:47 OK
dhad fit R 8.3 -t d -m 0 —qsub log-2007-07-02 13:39:12 FAILED.
dhad fit R 8.3 -t d -m 1 —qsub log-2007-07-02 13:39:20 OK
dhad fit R 8.3 -t d -m 3 —qsub log-2007-07-02 13:39:28 OK
dhad fit R 8.3 -t d -m 202 —qsub log-2007-07-02 13:39:37 OK
dhad fit regular -t s -m 203 —qsub log-2007-07-02 13:39:45 OK
dhad fit noFSR -t s -m 3 —qsub log-2007-07-02 13:39:55 OK
Only one failed mode.
Update Systmatic table
'Fractional uncertainty per double tag' : 0.002,
Summit the other regular jobs
dhad fit regular -t d -m 1 —qsub log-2007-07-02 15:07:23
dhad fit regular -t d -m 3 —qsub log-2007-07-02 15:07:33
dhad fit regular -t d -m 200 —qsub log-2007-07-02 15:07:46
dhad fit regular -t d -m 201 —qsub log-2007-07-02 15:07:56
dhad fit regular -t d -m 202 —qsub log-2007-07-02 15:08:07
dhad fit regular -t d -m 203 —qsub log-2007-07-02 15:08:17
dhad fit regular -t d -m 204 —qsub log-2007-07-02 15:08:30
dhad fit regular -t d -m 205 —qsub log-2007-07-02 15:08:42
dhad fit regular -t s -m 0 —qsub log-2007-07-02 15:09:04
dhad fit regular -t s -m 1 —qsub log-2007-07-02 15:09:14
dhad fit regular -t s -m 3 —qsub log-2007-07-02 15:09:24
dhad fit regular -t s -m 200 —qsub log-2007-07-02 15:09:35
dhad fit regular -t s -m 201 —qsub log-2007-07-02 15:09:44 Failed.
dhad fit regular -t s -m 202 —qsub log-2007-07-02 15:09:55
dhad fit regular -t s -m 203 —qsub log-2007-07-02 15:10:05
dhad fit regular -t s -m 204 —qsub log-2007-07-02 15:10:15
dhad fit regular -t s -m 205 —qsub log-2007-07-02 15:10:25
All the others are OK.
07/03/2007
Check the previous fits
- Failed fits
dhad fit R 8.3 -t d -m 0 —qsub log-2007-07-02 13:39:12
dhad fit regular -t s -m 201 —qsub log-2007-07-02 15:09:44
- Suggestion from Peter:
Change the RooArgusBG.cc (v1.3) to
L91: int ndm=40000;
Make snapshot for the previous code to: ndm10000
Backup the lib file as:
/home/xs32/work/CLEO/analysis/DHad/lib/RooFitModels/libRooFitModels_0.so
Check-in the ndm=40000 code, RooArgusBG.cc to (v1.4). Make snampshot name ndm40000
Compile and copy
gmake -f GNUmakefile.standalone cp $dhad/src/RooFitModels/tmp/libRooFitModels.so $dhad/lib/RooFitModels/libRooFitModels_1.so
dhad fit ndm40000 -t s -m 201 -l libRooFitModels_1.so —qsub log-2007-07-03 12:08:41 OK.
dhad fit R 8.3 -t d -m 0 -l libRooFitModels_1.so —qsub log-2007-07-03 12:09:10 OK.
Compare regular and ndm40000
/home/xs32/work/CLEO/analysis/DHad/7.06/dat/fit/regular/signal_Single_Dp_to_Kpipipi0__Dm_to_Kpipipi0.txt
Name | Value | Error |
---|---|---|
N1 | 20430.7061708558 | 152.5610860071 |
N2 | 20700.0144614164 | 153.6087891129 |
Nbkgd1 | 2853.4521428793 | 75.4643757358 |
Nbkgd2 | 2788.4466858067 | 75.3862148874 |
md | 1.8694285157 | 0.0000093903 |
p | 0.7237705190 | 0.0548593063 |
sigmap1 | 0.0041640505 | 0.0000522708 |
xi | -33.3319410287 | 2.8439010925 |
chisq1 | 39.4343574977 | 0.0000000000 |
chisq2 | 63.3017876614 | 0.0000000000 |
/home/xs32/work/CLEO/analysis/DHad/7.06/dat/fit/ndm40000/signal_Single_Dp_to_Kpipipi0__Dm_to_Kpipipi0.txt
Name | Value | Error |
---|---|---|
N1 | 20433.2424540185 | 152.2822892205 |
N2 | 20704.9605774617 | 153.3136035901 |
Nbkgd1 | 2855.7456309701 | 74.9571157617 |
Nbkgd2 | 2790.0200359321 | 74.8093280231 |
md | 1.8694286504 | 0.0000092604 |
p | 0.7188053268 | 0.0500125990 |
sigmap1 | 0.0041684345 | 0.0000509769 |
xi | -33.0744468683 | 2.5997952332 |
chisq1 | 39.3554922614 | 0.0000000000 |
chisq2 | 63.4520087330 | 0.0000000000 |
Update the BF
- Backup the regular file:
cp /home/xs32/work/CLEO/analysis/DHad/7.06/dat/fit/regular/signal_Single_Dp_to_Kpipipi0__Dm_to_Kpipipi0.txt /home/xs32/work/CLEO/analysis/DHad/7.06/dat/fit/regular/signal_Single_Dp_to_Kpipipi0__Dm_to_Kpipipi0.txt.bak
- Use the ndm40000 cover the regular file:
cp /home/xs32/work/CLEO/analysis/DHad/7.06/dat/fit/ndm40000/signal_Single_Dp_to_Kpipipi0__Dm_to_Kpipipi0.txt /home/xs32/work/CLEO/analysis/DHad/7.06/dat/fit/regular/signal_Single_Dp_to_Kpipipi0__Dm_to_Kpipipi0.txt
- Update the R syst
dhad table yields R_8.3 regular R_16.3 -t d
- Update the signal line shape syst
dhad table signal_line_shape_syst
- Get the BF
dhad bf
dhad table data_results_pdg dhad table data_results_bf_ratio_pdg
- Compare with noFSR
dhad bf noFSR
dhad table compare bf noFSR
Fit ndm40000
dhad proc fits ndm40000 regular -t d
dhad fit ndm40000 regular -t d -m 0 —qsub
07/04/2007
Fit ndm40000
dhad proc fits ndm40000 regular -t d
Test:
dhad fit ndm40000 M 3.7723 -t d -m 0 OK.
dhad proc fits ndm40000 M 3.7723 -t d
dhad proc fits ndm40000 M 3.7713 -t d
dhad proc fits ndm40000 Gamma 0.0311 -t d
dhad proc fits ndm40000 Gamma 0.0261 -t d
dhad proc fits ndm40000 R 16.3 -t d
dhad proc fits ndm40000 R 8.3 -t d
07/05/2007
Check the fitting results
Failed:
dhad fit ndm40000 M 3.7723 -t d -m 0 log-2007-07-05 00:45:52
All the others are OK.
Update the PRD draft
dhad table prd_line_shape_paras
dhad table prd_dt_dz_eff_yield
dhad table prd_dt_dp_eff_yield
dhad table prd_st_eff_yield
dhad table prd_data_results
dhad table prd_data_correlation_matrix
07/06/2007
Run MINOS for the failed fit
dhad fit ndm40000 M 3.7723 -t d -m 0 —set MINUIT=er4 —qsub log-2007-07-06 09:37:05
Problems.
Increase to 160,000 bins to fit
int ndm=160000; cd $dhad/src/RooFitModels gmake -f GNUmakefile.standalone cp tmp/libRooFitModels.so $dhad/lib/RooFitModels/libRooFitModels_2.so
dhad fit ndm160000 M 3.7723 -t d -m 0 —qsub log-2007-07-06 15:38:18 FAILED.
dhad fit ndm160000 M 3.7723 -t d -m 0 —set MINUIT=er4 —qsub log-2007-07-07 10:40:32
Remove err in noFSR comparison table
dhad table compare bf noFSR
Widede fit
dhad fit widede -t d OK.
dhad proc fits widede -t d widede
dhad proc fits widede -t s widede
dhad fit widede -t d —tag d —qsub log-2007-07-06 11:56:28 OK.
dhad fit widede -t s —tag d —qsub log-2007-07-06 11:56:47 OK.
07/07/2007
Check widede fits
Failed:
dhad fit widede -t d -m 0
Try:
dhad fit ndm40000 widede -t d -m 0 —qsub log-2007-07-07 10:46:28 OK.
Compare yields with ndm40000
dhad table yields regular ndm40000_regular -t d —set rnd=.001
07/08/2007
Update CBX
dhad table residual_double
dhad update cbx table
Widede BF fit
07/09/2007
fa plus minus fit
Message from Peter:
The RooFit files which I fixed the analytical integral for are: RooErf.rdl, RooErf.cc, RooDLineShape.rdl, RooDLineShape.cc all in /nfs/cor/user/ponyisi/xin/RooFitModels/.
RooErf.rdl to 1.2
RooErf.cc to 1.2
RooDLineShape.rdl to 1.2
RooDLineShape.cc 1.2
gmake -f GNUmakefile.standalone cp $dhad/src/RooFitModels/tmp/libRooFitModels.so $dhad/lib/RooFitModels/libRooFitModels_3.so
Use the ordinary code to do the fits for consistency
dhad fit faplus -t d -m 0 OK.
dhad proc fits faplus -t d faplus
dhad proc fits faplus -t s faplus
dhad fit faminus -t d -m 0 OK.
dhad proc fits faminus -t d faminus
dhad proc fits faminus -t s faminus
dhad fit saplus -t d -m 0 OK.
dhad proc fits saplus -t d saplus
dhad proc fits saplus -t s saplus
dhad fit saminus -t d -m 0 OK.
dhad proc fits saminus -t d saminus
dhad proc fits saminus -t s saminus
Failed fits:
dhad fit faplus -t d -m 0
dhad fit faplus -t s -m 0
dhad fit faplus -t s -m 205
dhad fit saplus -t s -m 0
dhad fit saminus -t s -m 200
dhad fit saminus -t s -m 200
Formating the Table in PRD
Message from Peter:
You may already know about it, but here is some code in Python to put commas in numbers with more than three digits: import locale locale.setlocale(locale.LC_ALL, '') locale.format('%.2f', 722291.432, True) # replace 722291.432 with your number
dhad table prd_dt_dz_eff_yield
dhad table prd_dt_dp_eff_yield
dhad table prd_st_eff_yield
OK.
07/10/2007
Check the fits
dhad check fits faplus -t d OK.
dhad check fits faplus -t s
-------------------------------------------------- dhad fit faplus -t s -m 205 -------------------------------------------------- FCN=-320296 FROM MIGRAD STATUS=FAILED 309 CALLS 310 TOTAL -------------------------------------------------- dhad fit faplus -t s -m 0 -------------------------------------------------- FCN=-1.51474e+06 FROM MIGRAD STATUS=FAILED 310 CALLS 311 TOTAL
dhad check fits faminus -t d OK.
dhad check fits faminus -t s
-------------------------------------------------- dhad fit faminus -t s -m 200 -------------------------------------------------- FCN=-1.64174e+06 FROM MIGRAD STATUS=FAILED 334 CALLS 335 TOTAL
dhad check fits saplus -t d OK.
dhad check fits saplus -t s
-------------------------------------------------- dhad fit saplus -t s -m 200 -------------------------------------------------- FCN=-1.64173e+06 FROM MIGRAD STATUS=FAILED 304 CALLS 305 TOTAL -------------------------------------------------- dhad fit saplus -t s -m 203 -------------------------------------------------- FCN=-510657 FROM MIGRAD STATUS=FAILED 480 CALLS 481 TOTAL -------------------------------------------------- dhad fit saplus -t s -m 202 -------------------------------------------------- FCN=-1.34838e+06 FROM MIGRAD STATUS=FAILED 371 CALLS 372 TOTAL
dhad check fits saminus -t d
-------------------------------------------------- dhad fit saminus -t d -m 204 -------------------------------------------------- FCN=-1.63936e+06 FROM MIGRAD STATUS=FAILED 361 CALLS 362 TOTAL
dhad check fits saminus -t s
-------------------------------------------------- dhad fit saminus -t s -m 205 -------------------------------------------------- FCN=-320297 FROM MIGRAD STATUS=FAILED 381 CALLS 382 TOTAL -------------------------------------------------- dhad fit saminus -t s -m 203 -------------------------------------------------- FCN=-510653 FROM MIGRAD STATUS=FAILED 554 CALLS 555 TOTAL -------------------------------------------------- dhad fit saminus -t s -m 0 -------------------------------------------------- FCN=-1.51474e+06 FROM MIGRAD STATUS=FAILED 286 CALLS 287 TOTAL
Submit (fixed the analytical integral) jobs for the failed fits
dhad fit faplus -t s -m 205 -l libRooFitModels_3.so —qsub log-2007-07-10 16:07:38 FAILED.
dhad fit faplus -t s -m 0 -l libRooFitModels_3.so —qsub log-2007-07-10 16:08:08 FAILED.
dhad fit faminus -t s -m 200 -l libRooFitModels_3.so —qsub log-2007-07-10 16:08:13 OK.
dhad fit saplus -t s -m 200 -l libRooFitModels_3.so —qsub log-2007-07-10 16:08:17 OK.
dhad fit saplus -t s -m 203 -l libRooFitModels_3.so —qsub log-2007-07-10 16:08:23 OK.
dhad fit saplus -t s -m 202 -l libRooFitModels_3.so —qsub log-2007-07-10 16:08:29 OK.
dhad fit saminus -t d -m 204 -l libRooFitModels_3.so —qsub log-2007-07-10 16:08:33 OK.
dhad fit saminus -t s -m 205 -l libRooFitModels_3.so —qsub log-2007-07-10 16:08:41 FAILED.
dhad fit saminus -t s -m 203 -l libRooFitModels_3.so —qsub log-2007-07-10 16:08:47 OK.
dhad fit saminus -t s -m 0 -l libRooFitModels_3.so —qsub log-2007-07-10 16:08:52 OK.
- try ndm 160000 for the failed modes:
dhad fit faplus -t s -m 205 -l libRooFitModels_2.so —qsub log-2007-07-10 16:58:24 OK.
dhad fit faplus -t s -m 0 -l libRooFitModels_2.so —qsub log-2007-07-10 16:58:41 OK.
dhad fit saminus -t s -m 205 -l libRooFitModels_2.so —qsub log-2007-07-10 16:58:59 FAILED.
07/11/2007
Add FSR table into the PRD draft
dhad table prd_data_results
Done.
Check the fitting result
One failed (fixed the analytical integral, ndm160000):
dhad fit saminus -t s -m 205
07/12/2007
Update PRD results
07/13/2007
Submit fixed the analytical integral + ndm640000
RooArgusBG.cc:
int ndm=640000;
cd $dhad/src/RooFitModels gmake -f GNUmakefile.standalone cp tmp/libRooFitModels.so $dhad/lib/RooFitModels/libRooFitModels_4.so
dhad fit saminus -t s -m 205 -l libRooFitModels_4.so —qsub
log-2007-07-13 11:06:50 FAILED.
07/16/2007
Try more precision in the integrator
In mbc_sigletag_3s.C add:
RooNumIntConfig fcn1_1config(*fcn1_1.getIntegratorConfig()); fcn1_1config.setEpsAbs(1e-8); fcn1_1config.setEpsRel(1e-8); fcn1_1.setIntegratorConfig(fcn1_1config); RooNumIntConfig fcn1_2config(*fcn1_2.getIntegratorConfig()); fcn1_2config.setEpsAbs(1e-8); fcn1_2config.setEpsRel(1e-8); fcn1_2.setIntegratorConfig(fcn1_2config); RooNumIntConfig fcn1_3config(*fcn1_3.getIntegratorConfig()); fcn1_3config.setEpsAbs(1e-8); fcn1_3config.setEpsRel(1e-8); fcn1_3.setIntegratorConfig(fcn1_3config); RooNumIntConfig fnc2_1config(*fnc2_1.getIntegratorConfig()); fnc2_1config.setEpsAbs(1e-8); fnc2_1config.setEpsRel(1e-8); fnc2_1.setIntegratorConfig(fnc2_1config); RooNumIntConfig fnc2_2config(*fnc2_2.getIntegratorConfig()); fnc2_2config.setEpsAbs(1e-8); fnc2_2config.setEpsRel(1e-8); fnc2_2.setIntegratorConfig(fnc2_2config); RooNumIntConfig fnc2_3config(*fnc2_3.getIntegratorConfig()); fnc2_3config.setEpsAbs(1e-8); fnc2_3config.setEpsRel(1e-8); fnc2_3.setIntegratorConfig(fnc2_3config);
Save as:
mbc_singletag_3s_1.C
dhad fit saminus -t s -m 205 -l mbc_singletag_3s_1.C —qsub log-2007-07-16 14:22:13 FAILED.
07/18/2007
Compare yields between regular and faminus
dhad table yields regular faminus -t d
dhad table yields regular faplus -t d
dhad table yields regular saminus -t d
dhad table yields regular saplus -t d
07/19/2007
Update CBX
dhad update cbx table
07/20/2007
Update CBX
dhad update cbx table
Ready for updating the tables
07/23/2007
Update CBX
dhad update cbx table
07/24/2007
Mordify the tables in PRD draft
Message from Steve:
Numbers: [page #]-[comment # on page] -------- XIN --------- 10-1: more space between columns 11-2, 12-1, & 13-1: align $\pm$'s within each column
Done.
07/25/2007
Update CBX
dhad update cbx table
07/27/2007
Update CBX
dhad update cbx table
Done with table 3 and 4.
07/29/2007
Update CBX
dhad table cbx_double_data_gamma
dhad table cbx_single_signal_eff
dhad table cbx_single_data_gamma
OK.
- cbx_single_widede
dhad fit widede noFSR -t s -m 0 OK. dhad proc fits widede noFSR -t s Save plots as: widede noFSR
dhad table cbx_single_widede
- cbx_single_fix_p
dhad proc fits P 0.5 -t d Save plots as: P 0.5
dhad proc fits P 0.5 -t s Save plots as: P 0.5
dhad table cbx_single_fix_p
- cbx_compare_yields_DDbar
dhad table cbx_compare_yields_DDbar
OK.
dhad table cbx_noFSR_single_signal_eff OK.
dhad table cbx_data_results_pdg OK.
dhad table cbx_data_results_bf_ratio_pdg OK.
dhad table cbx_data_correlation_matrix OK.
dhad table cbx_residual_single OK.
dhad table cbx_residual_double OK.
dhad update cbx table
Need to be done:
dhad table cbx_data_results_variations
BF fits
dhad bf widede
dhad bf Gamma 0.0261
dhad bf Gamma 0.0311
07/30/2007
Check fitting results
dhad check fits widede noFSR -t s
-------------------------------------------------- dhad fit widede noFSR -t s -m 204 -------------------------------------------------- FCN=-470584 FROM MIGRAD STATUS=FAILED 406 CALLS 407 TOTAL
Move on.
dhad check fits P 0.5 -t d OK.
dhad check fits P 0.5 -t s OK.
Update CBX table
dhad table cbx_single_widede OK.
dhad table cbx_single_fix_p OK.
dhad update cbx table
Various BF fits
P = 0.5 double tag fit
dhad fit P 0.5 -t s —tag d —qsub log-2007-07-30 10:20:40
dhad fit P 0.5 -t d —tag d —qsub log-2007-07-30 10:20:57
Make plots page. Done.
dhad check fits P 0.5 -t s —tag d dhad check fits P 0.5 -t d —tag d OK.
dhad bf P 0.5
dhad table cbx_data_results_variations
dhad update cbx table
Modify lineshapepara table
dhad table cbx_line_shape_paras
dhad table cbx_double_signal_eff
edit: dhad table cbx_double_generic_eff
dhad update cbx table
Update CBX plots
dhad update cbx plots done.
Update CBX results
From the RAW file of the bf results:
sigma(D0D0bar) = 3660.64 +- 28.7876 (stat) +- 59.7543 (syst) sigma(D+D-) = 2910.02 +- 29.5951 (stat) +- 46.2134 (syst) Correlation coeff between sigma(D0D0bar) and sigma(D+D-): 0.570614 sigma(DDbar) = 6570.66 +- 41.2868 (stat) +- 99.2985 (syst) chg/neu = 0.794948 +- 0.0103848 (stat) +- 0.0088281 (syst)
07/31/2007
Alternative BF fits
Message from Anders:
We should do two more alternative branching fraction fits: 1) Setting the systematic uncertainty due to FSR to zero. 2) Setting the systematic uncertainty due to the double DCSD interference to zero. (Actually I guess we might need to do one fit where we set them both to zero.) I'd like to do these two fits to see what the remaining systematic uncertainty is due to the pure experimental effects. Lawrence and I discussed that we should consider quoting the results for the branching fractions as brfr. +- stat +- exp syst +- FSR syst. +- double DCSD syst. (The last only applies for D0 channels.)
- FSR syst. zero fit
dhad bf syst FSR 0
done.
- Double DCSD interference syst.zero fit
dhad bf syst DCSD 0
- Both Zero fit
dhad bf syst FSR 0 DCSD 0
Comparison table FSR DCSD
dhad table compare bf syst/FSR/0 —set rnd=.001
dhad table compare bf syst/DCSD/0 —set rnd=.001
dhad table compare bf syst/FSR/0/DCSD/0 —set rnd=.001
Rerun double tag fits
Message from Anders:
I looked at the logfile for your doubletag fits and noticed that we are rejecting 35 out of 25131 events. This is a 0.1% effect, and not a very large effect. However, we should probably try to correct this. As a first step, could you rerun the fits for the double tag yields? In the script you use for the DT fits you should have a line like: RooRealVar Ebeam("Ebeam","Ebeam",1.8854,1.8911,"GeV"); can you change this to: RooRealVar Ebeam("Ebeam","Ebeam",1.882,1.8911,"GeV"); and rerun the fits.
Save to lineshapefit2d_1.C
dhad fit regular1 -t d —tag d -l lineshapefit2d_1.C —qsub