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

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.

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

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

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

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
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

dhad table yields R_8.3 regular R_16.3 -t d

dhad table signal_line_shape_syst

dhad bf

dhad table data_results_pdg dhad table data_results_bf_ratio_pdg

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

ndm40000 regular

Test:

dhad fit ndm40000 M 3.7723 -t d -m 0 OK.

dhad proc fits ndm40000 M 3.7723 -t d

ndm40000 M 3.7723

dhad proc fits ndm40000 M 3.7713 -t d

ndm40000 M 3.7713

dhad proc fits ndm40000 Gamma 0.0311 -t d

ndm40000 Gamma 0.0311

dhad proc fits ndm40000 Gamma 0.0261 -t d

ndm40000 Gamma 0.0261

dhad proc fits ndm40000 R 16.3 -t d

ndm40000 R 16.3

dhad proc fits ndm40000 R 8.3 -t d

ndm40000 R 8.3

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.

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.

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

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

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:

Raw 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.)

dhad bf syst FSR 0

done.

dhad bf syst DCSD 0

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

log-2007-07-31 17:01:50