D Hadronic Analysis - Log - June 2006

<<< 05/31/2006

06/01/2006

Float the width of the 3770 in Single tag

For data

DHadFit.py -a 4

qsub gfloat.sh
Your job 394109 ("gfloat.log") has been submitted.

See plots

Single tag yields for generic MC

06/02/2006

P

Check the gfloat.log, only first one has the output. Need to be fixed. DHadronic_Log_200606: 05

• Write up the report

Select the widede data evt:

DHad.py  -t d -p widede

06/03/2006

• Write up the report

06/04/2006

• Write report

1. FSR event number

DHad.py -t s -p noFSR

DHadTable.py -a 5

Refresh the table, only concern the tex format.

06/05/2006

From DHadronic_Log_200606: 02 : Check the gfloat.log, only first one has the output. Need to be fixed.

Check: DHadFit.py -a 4

Using more intuitive options.

06/06/2006

Make the float run:

DHadFit.py -t d -p Peter_Gamma_float -m 1

Message :

RooRealVar::isValid(Ebeam): value 1.87679 out of range (1.8815 - 1.892)

In mbc_singletag_3s.C:

RooRealVar Ebeam("Ebeam","Ebeam",1.8815,1.892,"GeV");

Although it can generate pictures. Send mail to Peter.

Test all modes by printing the mode number:

DHadFit.py -t d -p Peter_Gamma_float

OK. Then submit the qsub job:

qsub gfloat.sh
Your job 397764 ("gfloat.log") has been submitted.

• Choose regular data

06/07/2006

Check the Float width fit

Look at gfloat.log397764, still the same problem.

Error: Ambiguous overload resolution (1000001,2) FILE:/home/xs32/work/CLEO/analysis/DHad/script/Roofit/mbc_singletag_3s.C LINE:599
Calling : basic_ofstream<char,char_traits<char> >::basic_ofstream<char,char_traits<char> >(TString);

Add "table.close();" in the .C file. Resubmit:

qsub gfloat.sh
Your job 397899 ("gfloat.log") has been submitted.

Choose regular data

edit : data_DHad.sh

qsub data_DHad.sh
Your job 397912 ("data_DHad.log") has been submitted.

Diagonal double tag fit with parameters print out

Table 3 in the CBX.

DHadFit.py -t d --tag d -p Peter_diag -m 0

qsub DHadFit_d_d_Peter_diag.sh
Your job 399765 ("DHadFit_d_d_Peter_diag.log") has been submitted.

06/08/2006

Check the fit result

Peter_Gamma_float

gfloat.log.o397899 , still the same symptom. So change the sh file and resubmit.

qsub DHadFit_d_Peter_Gamma_float.sh
Your job 400387 ("DHadFit_d_Peter_Gamma_float.log") has been submitted.

Peter_diag

Check DHadFit_d_d_Peter_diag.log.o399765

See plot:

• Diagnal Double Tag fit using Peter's Evt Data fit

Fit the widede data and signal MC for tab

DHadFit.py -t d -p widede -m 0

Create sh file and qsub

qsub DHadFit_d_widede.sh
Your job 400439 ("DHadFit_d_widede.log") has been submitted.

For the signal MC:

qsub DHadFit_widede.sh
Your job 400446 ("DHadFit_widede.log") has been submitted.

test the ofstream

DHadFit.py -t d

Error: Ambiguous overload resolution (1000001,2) FILE:/home/xs32/work/CLEO/analysis/DHad/script/Roofit/mbc_singletag_3s.C LINE:599
Calling : basic_ofstream<char,char_traits<char> >::basic_ofstream<char,char_traits<char> >(TString);
Match rank: file     line  signature
*  1000001 (compiled)   0 basic_ofstream<char,char_traits<char> > basic_ofstream<char,char_traits<char> >::basic_ofstream<char,char_traits<char> >(const char*,ios_base::openmode);
  ffffffff (compiled)   0 basic_ofstream<char,char_traits<char> > basic_ofstream<char,char_traits<char> >::basic_ofstream<char,char_traits<char> >();
*  1000001 (compiled)   0 basic_ofstream<char,char_traits<char> > basic_ofstream<char,char_traits<char> >::basic_ofstream<char,char_traits<char> >(const char*,ios_base::openmode);
  ffffffff (compiled)   0 basic_ofstream<char,char_traits<char> > basic_ofstream<char,char_traits<char> >::basic_ofstream<char,char_traits<char> >();
*  1000001 (compiled)   0 basic_ofstream<char,char_traits<char> > basic_ofstream<char,char_traits<char> >::basic_ofstream<char,char_traits<char> >(const char*,SZ_T);
  ffffffff (compiled)   0 basic_ofstream<char,char_traits<char> > basic_ofstream<char,char_traits<char> >::basic_ofstream<char,char_traits<char> >();
*** Interpreter error recovered ***

code:

599:  ofstream table(txtfile);
615:  table.precision(10);

06/12/2006

Check the previous fits

• DHadFit_d_widede.log.o400439 On mode 205:

  File "/home/xs32/local/bin/DHadFit.py", line 218
    diag_double_tag_fit(prefix, dt_type, tag, mode)
                                                  ^
IndentationError: unindent does not match any outer indentation level

Check the other modes:

Gamma should not appear in the parameters. So, fix:

1. Make the mode 205 work in data 2. Make sure no gamma in the parameter list

1. > DHadFit.py -t d -p widede -m 205 2. > Add lines in the mbc_singletag_3s.C:

  if (floatwidth != 0){
  table << "Gamma"  <<"\t|  "<< Gamma    << "\t|  " << Gamma_err  << endl;
  }

06/14/2006

Check the regular data fit

qsub DHadFit_d.sh
Your job 405208 ("DHadFit_d.log") has been submitted.

OK.

Check the regular single tag fit

Reselect the evt

qsub DHad.sh
Your job 405253 ("DHad.log") has been submitted.

log at dat/evt/log/

OK.

qsub DHadFit.sh
Your job 405284 ("DHadFit.log") has been submitted.

&gt;&gt;&gt;

Count the Generated number of Generic MC

How to count?

DHad.py -t g -m 0

&gt;&gt;&gt;

Fix the table problem in the modes loop

Use the test version during the qsub period to avoid messing things up.

test-DHadFit.py -p Peter_evt

Currently is OK for 2 modes without plotting.

Try plot with all the command:

test-DHadFit.py -p Peter_evt

Then the error occured again as below:

Error: Ambiguous overload resolution (1000001,2) FILE:/home/xs32/work/CLEO/analysis/DHad/script/Roofit/test_mbc_singletag_3s.C LINE:473
Calling : basic_ofstream<char,char_traits<char> >::basic_ofstream<char,char_traits<char> >(TString);
Match rank: file     line  signature
*  1000001 (compiled)   0 basic_ofstream<char,char_traits<char> > basic_ofstream<char,char_traits<char> >::basic_ofstream<char,char_traits<char> >(const char*,ios_base::openmode);
  ffffffff (compiled)   0 basic_ofstream<char,char_traits<char> > basic_ofstream<char,char_traits<char> >::basic_ofstream<char,char_traits<char> >();
*  1000001 (compiled)   0 basic_ofstream<char,char_traits<char> > basic_ofstream<char,char_traits<char> >::basic_ofstream<char,char_traits<char> >(const char*,SZ_T);
  ffffffff (compiled)   0 basic_ofstream<char,char_traits<char> > basic_ofstream<char,char_traits<char> >::basic_ofstream<char,char_traits<char> >();
*** Interpreter error recovered ***
S

So, the conflict should be between the plot and writing out txt file.

Comment out the fitting part and see. From this plot starting, the conflict starts:

TCanvas *canvas= new TCanvas("canvas","mbc",900,300);

&gt;&gt;&gt;

No FSR Event Check

test-DHad.py

&gt;&gt;&gt;

06/15/2006

Check the regular single tag fit

Check the DHadFit.log.o40528 ... OK.

Check the Regular fit- signal MC ... OK.

Debug the ofstream in more modes

Use the C style way:

#include <stdio.h>
...

FILE* table = fopen(txtfile.Data(), "w+");
fprintf(table, "Name\t||Value\t|| Error\n");
fprintf(table, "N1\t|%10f\t|%10f\n", yield1, yield1_err);

fclose(table);

Apply the change to the mbc_singletag_3s.C. Test:

DHadFit.py -p Peter_evt

Check the digits of decimal ... OK.

Now, ready for more modes fit.

Count the Generated number of Generic MC

Wait for Peter's code.

&gt;&gt;&gt;

No FSR Event Check

Seclect the event

Make sure no qsub job is working. Use the DHad.py directly for testing.

DHad.py -p noFSR

List the raw efficiency.

DHadTable.py -a 7

Input from Anders

FSR events should be selected only for the interested mode.

Look at the mcdmode doc.

Use the "mcdmode" and "mcdbmode" in the dictionary. &gt;&gt;&gt;

Compare the evt file with Peter's

Data

/nfs/cor/user/ponyisi/hadD/summerconf/txt/data_Single_D0_to_Kpi.evt /home/xs32/work/CLEO/analysis/DHad/dat/evt/data_Single_D0_to_Kpi.evt

They are different! There must be something wrong with the root file of the data. &gt;&gt;&gt;

Signal

/nfs/cor/user/ponyisi/hadD/summerconf/txt/signal_Single_D0_to_Kpi.evt /home/xs32/work/CLEO/analysis/DHad/dat/evt/signal_Single_D0_to_Kpi.evt

They are the same.

Fit Peter's data evt

Test DHadFit.py -t d -p Peter_evt ... OK. Submit qsub:

qsub DHadFit_d_Peter_evt.sh
Your job 405708 ("DHadFit_d_Peter_evt.log") has been submitted.

&gt;&gt;&gt;

Fit the yields

Test DHadFit.py -p noFSR ... OK.

Send to the qsub job:

qsub DHadFit_noFSR.sh
Your job 405651 ("DHadFit_noFSR.log") has been submitted.

As the event selection is not correct, this fit is not very useful. &gt;&gt;&gt;

Next

Change the r to 12.7 in some modes to compare the difference.

mbc_singletag_3s.C

r=new RooRealVar("r","r",15.0);

&gt;&gt;&gt;

06/16/2006

Create Reference Page

&lt;&lt;&lt;

Hold the numbers of the generated MC

Tables from Peter etc.

&gt;&gt;&gt;

No FSR Event Selection

&lt;&lt;&lt;

Use the "mcdmode" and "mcdbmode" in the dictionary

DHad.py -p noFSR

The efficiencies looks OK.

Data NTuple difference

&lt;&lt;&lt; &gt;&gt;&gt;

Check Peter's data evt fit

&lt;&lt;&lt;

Check DHadFit_d_Peter_evt.log.o405708 ... OK. Make plots Need to change the table format. ... Done. The plots are at:

• Regular fit- data

Fit the new FSR signal MC

&lt;&lt;&lt;

qsub DHadFit_noFSR.sh
Your job 406110 ("DHadFit_noFSR.log") has been submitted.

&gt;&gt;&gt;

Use r = 12.7 to fit certain modes

&lt;&lt;&lt;

test-DHadFit.py -p R_12.7 -m 0 ... OK.

See plot: R = 12.7 Signal MC

Compare the yield:

Name	Regular Signal MC	R = 12.7
N1	20617.326 &plusmn 143.834	20617.326 &plusmn 143.834
N2	20950.098 &plusmn 144.993	20950.098 &plusmn 144.993

Check the R value is really 12.7 ... True.

Hardcode that and try again.

test-DHadFit.py -p R_12.7 -m 0

&gt;&gt;&gt;

06/19/2006

Create table for generated generic MC

&lt;&lt;&lt;

See Table.

Data NTuple Difference

&lt;&lt;&lt; Send email to Peter.

It turns out that I didn't inclue the /cdat/dafc/dhad/summerconf/data/data37_dskim_evtstore__1.root

Polish the script ... OK

Check the no FSR signal MC Fit

&lt;&lt;&lt;

Check DHadFit_noFSR.log.o406110 ... OK Outputfile, table is not good. Need to change the format ... OK

06/20/2006

Select the Data evt

/nfs/cor/an1/ponyisi/dhad/summerconf/data_srs/

Before selection, backup the old one first. evt, fig, dat/fit, : old_data

Test DHad.py -t d -m 0 ... OK.

qsub DHad_d.sh
Your job 412149 ("DHad_d.log") has been submitted.

Compare the evt with Peter's

/nfs/cor/user/ponyisi/hadD/summerconf/txt/data_Single_D0_to_Kpi.evt /home/xs32/work/CLEO/analysis/DHad/dat/evt/data_Single_D0_to_Kpi.evt

They are the same.

So, the difference comes out from not including the rootfile:

/cdat/dafc/dhad/summerconf/data/data37_dskim_evtstore__1.root

• Learn how to get email notice from qsub

#$ -m bes

Submit the r = 12.7 job for data

Use Peter's data evt first

DHadFit.py -t d -p Peter_evt_R_12.7 -m 0

See plot R = 12.7 Data

Use my evt file

DHadFit.py -t d -p R_12.7 -m 0

Refresh data related fit

Regular fit

Test DHadFit.py -t d ... OK.

qsub DHadFit_d.sh ("DHadFit_d.log.o412714") ... OK.

Wide Delta E

Test DHad.py -t d -p widede OK.

qsub DHad_d_widede.sh ("DHad_d_widede.log.o412715") ... OK.

Gamma = 35.0 MeV

Test DHadFit.py -t d -p Gamma_0.035 ... OK.

qsub DHadFit_d_Gamma_0.035.sh

("DHadFit_d_Gamma_0.035.log.o412716") ... OK.

Gamma = 23.6 MeV

Test DHadFit.py -t d -p Gamma_0.0236 ... OK.

qsub DHadFit_d_Gamma_0.0236.sh ("DHadFit_d_Gamma_0.0236.log.o412717") ... OK.

Gamma = 27 MeV

Test DHadFit.py -t d -p Gamma_0.027 -m 0,1,3 ... OK.

qsub DHadFit_d_Gamma_0.027_013.sh ("DHadFit_d_Gamma_0.027_013.log.o412718") ... OK.

Gamma = 20 MeV

Test DHadFit.py -t d -p Gamma_0.02 -m 200,201,202,203,204,205 ... OK.

qsub DHadFit_d_Gamma_0.02_200_205.sh ("DHadFit_d_Gamma_0.02_200_205.log.o412732") ... OK.

Float the width of Psi(3770)

Test DHadFit.py -t d -p Gamma_float ... OK.

qsub DHadFit_d_Gamma_float.sh ("DHadFit_d_Gamma_float.log.o412733") ... OK.

06/21/2006

Check the data related fits

Creating pages

• Regular Fit
• Gamma = 35.0 MeV
• Gamma = 23.6 MeV
• D0D0B : Gamma 27 MeV
• D+ D- : Gamma 20 MeV
• Float the width of Psi(3770)
• R = 12.7

Wide Delta E

Test DHadFit.py -t d -p widede ... OK.

qsub DHadFit_d_widede.sh ("DHadFit_d_widede.log.o415049") ... OK.

Refresh the fits for the signal MC

p = 0.5 fit

Test DHadFit.py -p p0.5 ... OK. qsub DHadFit_p0.5.sh - job 415079 ("DHadFit_p0.5.log") For mode 0, failed to fit.

One of the evt file missed.

DHad.py ... OK.

qsub DHadFit_p0.5.sh - job 415180 ("DHadFit_p0.5.log") ... OK.

Wide Delta E

DHad.py

Test DHadFit.py -p widede ... OK. qsub DHadFit_widede.sh job 415181 ("DHadFit_widede.log") ... OK.

06/22/2006

Check the previous fits

OK.

FSR table

Redo the fits for the noFSR

Test DHadFit.py -p noFSR ... OK. qsub DHadFit_noFSR.sh job 416099 ("DHadFit_noFSR.log") ... Event number is wrong

06/26/2006

Fix the FSR fit

Check the DHadFit.py ... OK. Test DHadFit.py -p noFSR ... OK.

qsub DHadFit_noFSR.sh Your job 416337 ("DHadFit_noFSR.log") has been submitted. ... OK.

06/27/2006

Work on the table

DHadTable.py -a 5

Got the table

Handle the ratio error

        eff_noFSR     Y_noFSR/N_noFSR     Y_noFSR       N
ratio = ---------  =  ----------------  = ------- * --------
           eff            Y / N              Y       N_noFSR


let
             Y_noFSR                 N_noFSR
  ratio_Y = ---------  , ratio_N =  ---------
	        Y    	                N

                  ------------------------
  ratio_Y_err = \/ Y * ratio_Y (1-ratio_Y)

                  ------------------------
  ratio_N_err = \/ Y * ratio_N (1-ratio_N)

                ------------------------------
  ratio_err = \/ ratio_Y_err^2 + ratio_N_err^2

06/28/2006

Implement the error calculation

Use the Psi(3770) mass of 3.774 GeV to fit the width

Test DHadFit.py -t d -p Gamma_float_M_3.774 ... OK.

qsub DHadFit_d_Gamma_float_M_3.774.sh Your job 418731 ("DHadFit_d_Gamma_float_M_3.774.log") ... OK.

06/29/2006

Generate the plot page for the M 3.774 GeV

• Float the width of Psi(3770), with mass = 3.774 GeV

06/30/2006

Moving on the table calculation

[

Base Input

Auxiliary

Final ]

Update the CBX

>>> 07/03/2006