D Had Log - November 2006

DHad.py fit fa 0.142 -m 0 --qsub
job 53353

DHad.py fit fa 0.152 -m 0 --qsub
job 53354

DHad.py fit fa 0.162 -m 0 --qsub
job 53355

DHad.py fit fa 0.172 -m 0 --qsub
job 53356

DHad.py fit fa 0.182 -m 0 --qsub
job 53357

DHad.py fit fa 0.192 -m 0 --qsub
job 53358

DHad.py fit fa 0.202 -m 0 --qsub
job 53359

DHad.py fit fa 0.212 -m 0 --qsub
job 53360

DHad.py fit fa 0.222 -m 0 --qsub
job 53361

Plots list

For Kpipi mode, fa = 0.196 +/- 0.036,

DHad.py fit fa 0.196+/-0.036,10 -m 1 --qsub --commands

/home/xs32/local/bin/DHad.py fit fa 0.16 -m 1 --qsub
Your job 53364 ("DHad_fit_fa_0.16_-m_1.log") has been submitted

/home/xs32/local/bin/DHad.py fit fa 0.1672 -m 1 --qsub
Your job 53365 ("DHad_fit_fa_0.1672_-m_1.log") has been submitted

/home/xs32/local/bin/DHad.py fit fa 0.1744 -m 1 --qsub
Your job 53366 ("DHad_fit_fa_0.1744_-m_1.log") has been submitted

/home/xs32/local/bin/DHad.py fit fa 0.1816 -m 1 --qsub
Your job 53367 ("DHad_fit_fa_0.1816_-m_1.log") has been submitted

/home/xs32/local/bin/DHad.py fit fa 0.1888 -m 1 --qsub
Your job 53368 ("DHad_fit_fa_0.1888_-m_1.log") has been submitted

/home/xs32/local/bin/DHad.py fit fa 0.196 -m 1 --qsub
Your job 53369 ("DHad_fit_fa_0.196_-m_1.log") has been submitted

/home/xs32/local/bin/DHad.py fit fa 0.2032 -m 1 --qsub
Your job 53370 ("DHad_fit_fa_0.2032_-m_1.log") has been submitted

/home/xs32/local/bin/DHad.py fit fa 0.2104 -m 1 --qsub
Your job 53371 ("DHad_fit_fa_0.2104_-m_1.log") has been submitted

/home/xs32/local/bin/DHad.py fit fa 0.2176 -m 1 --qsub
Your job 53372 ("DHad_fit_fa_0.2176_-m_1.log") has been submitted

/home/xs32/local/bin/DHad.py fit fa 0.2248 -m 1 --qsub
Your job 53373 ("DHad_fit_fa_0.2248_-m_1.log") has been submitted

Consider using the Roofit inner function to finish this job

source p.55

RooPlot* frame = sin2b.frame(0,1,20) ;
cpmixPdf->plotNLLOn(frame,cpmixDada,1.0,kTRUE) ;

DHad.py fit fa plotNLL -m 0

12/04/2006

Using the exsiting one and make the plot

12/05/2006

plot fa vs. logL

For mode 0:

Get the table of minNll and fa

DHad.py table fa,minNll -m 0

Consider using the Roofit inner function for Data

source p.55

RooPlot* frame = sin2b.frame(0,1,20) ;
cpmixPdf->plotNLLOn(frame,cpmixDada,1.0,kTRUE) ;

DHad.py fit fa plotNll -t d -m 0

pi0 mass in Kpipi0 for Data

DHad.py plot mass pi0  -t d -m 1  --qsub

job 54211

12/06/2006

plot fa vs. logL

DHad.py fit fa plotNll -t d -m 0

Error message

Hit fraction shrink to 60%(70,80,90) in the Ntuple DT/MC

variable in the ntuple: trhitf[ntrack]
Select events for 60% print it out

DHad.py yield trhitf 0.6 —qsub job 54715

DHad.py fit trhitf 0.6

DHad.py yield trhitf 0.6 -t d —qsub job 54716

DHad.py yield trhitf 0.6 -t g —qsub job 54717

DHad.py yield trhitf 0.7 -t s —qsub job 54718

DHad.py yield trhitf 0.7 -t d —qsub job 54719

DHad.py yield trhitf 0.7 -t g —qsub job 54720

DHad.py yield trhitf 0.8 -t s —qsub job 54721

DHad.py yield trhitf 0.8 -t d —qsub job 54722

DHad.py yield trhitf 0.8 -t g —qsub job 54723

DHad.py yield trhitf 0.9 -t s —qsub job 54724

DHad.py yield trhitf 0.9 -t d —qsub job 54725

DHad.py yield trhitf 0.9 -t g —qsub job 54726

Error happend:

in mode 1:

index = 3, means 4 tracks, but only 3 tracks has trhitf info in the ntuple. Need to redo.

12/07/2006

fa loglikelihood

float fa plot

DHad.py fit fa float -t d --qsub
job 54867

DHad.py fit fa float -t s --qsub
job 54868

DHad.py fit fa float -t g --qsub
job 54869

../plots/data_fa_float

Organize the previous fit

The table location for mode 0:

~/work/CLEO/analysis/DHad/dat/fit/fa/0.142/signal_Single_D0_to_Kpi__D0B_to_Kpi.txt

fa from 0.142 to 0.222

Fit the data:

DHad.py fit fa 0.182+/-0.045,10 -m 1 -t d --qsub --commands
/home/xs32/local/bin/DHad.py fit fa 0.137 -m 1 -t d --qsub
Your job 54847 ("DHad_fit_fa_0.137_-m_1_-t_d.log") has been submitted
/home/xs32/local/bin/DHad.py fit fa 0.146 -m 1 -t d --qsub
Your job 54848 ("DHad_fit_fa_0.146_-m_1_-t_d.log") has been submitted
/home/xs32/local/bin/DHad.py fit fa 0.155 -m 1 -t d --qsub
Your job 54849 ("DHad_fit_fa_0.155_-m_1_-t_d.log") has been submitted
/home/xs32/local/bin/DHad.py fit fa 0.164 -m 1 -t d --qsub
Your job 54850 ("DHad_fit_fa_0.164_-m_1_-t_d.log") has been submitted
/home/xs32/local/bin/DHad.py fit fa 0.173 -m 1 -t d --qsub
Your job 54851 ("DHad_fit_fa_0.173_-m_1_-t_d.log") has been submitted
/home/xs32/local/bin/DHad.py fit fa 0.182 -m 1 -t d --qsub
Your job 54852 ("DHad_fit_fa_0.182_-m_1_-t_d.log") has been submitted
/home/xs32/local/bin/DHad.py fit fa 0.191 -m 1 -t d --qsub
Your job 54853 ("DHad_fit_fa_0.191_-m_1_-t_d.log") has been submitted
/home/xs32/local/bin/DHad.py fit fa 0.2 -m 1 -t d --qsub
Your job 54854 ("DHad_fit_fa_0.2_-m_1_-t_d.log") has been submitted
/home/xs32/local/bin/DHad.py fit fa 0.209 -m 1 -t d --qsub
Your job 54855 ("DHad_fit_fa_0.209_-m_1_-t_d.log") has been submitted
/home/xs32/local/bin/DHad.py fit fa 0.218 -m 1 -t d --qsub
Your job 54856 ("DHad_fit_fa_0.218_-m_1_-t_d.log") has been submitted

DHad.py fit fa 0.196+/-0.036,10 -m 1 -t d --qsub --commands
/home/xs32/local/bin/DHad.py fit fa 0.16 -m 1 -t d --qsub
Your job 54857 ("DHad_fit_fa_0.16_-m_1_-t_d.log") has been submitted
/home/xs32/local/bin/DHad.py fit fa 0.1672 -m 1 -t d --qsub
Your job 54858 ("DHad_fit_fa_0.1672_-m_1_-t_d.log") has been submitted
/home/xs32/local/bin/DHad.py fit fa 0.1744 -m 1 -t d --qsub
Your job 54859 ("DHad_fit_fa_0.1744_-m_1_-t_d.log") has been submitted
/home/xs32/local/bin/DHad.py fit fa 0.1816 -m 1 -t d --qsub
Your job 54860 ("DHad_fit_fa_0.1816_-m_1_-t_d.log") has been submitted
/home/xs32/local/bin/DHad.py fit fa 0.1888 -m 1 -t d --qsub
Your job 54861 ("DHad_fit_fa_0.1888_-m_1_-t_d.log") has been submitted
/home/xs32/local/bin/DHad.py fit fa 0.196 -m 1 -t d --qsub
Your job 54862 ("DHad_fit_fa_0.196_-m_1_-t_d.log") has been submitted
/home/xs32/local/bin/DHad.py fit fa 0.2032 -m 1 -t d --qsub
Your job 54863 ("DHad_fit_fa_0.2032_-m_1_-t_d.log") has been submitted
/home/xs32/local/bin/DHad.py fit fa 0.2104 -m 1 -t d --qsub
Your job 54864 ("DHad_fit_fa_0.2104_-m_1_-t_d.log") has been submitted
/home/xs32/local/bin/DHad.py fit fa 0.2176 -m 1 -t d --qsub
Your job 54865 ("DHad_fit_fa_0.2176_-m_1_-t_d.log") has been submitted
/home/xs32/local/bin/DHad.py fit fa 0.2248 -m 1 -t d --qsub
Your job 54866 ("DHad_fit_fa_0.2248_-m_1_-t_d.log") has been submitted

DHad.py table fa,minNll -m 0

fa	minNll
0.142	-737707.854
0.152	-737707.667
0.162	-737707.468
0.172	-737707.270
0.182	-737707.073
0.192	-737706.884
0.202	-737706.701
0.212	-737706.529
0.222	-737706.363

12/15/2006

plot fa vs. logL

DHad.py fit fa plotNll -t d -m 0

f2 == 0.000000 ---- New Nll == -981466.782695
f2 == 0.025000 ---- New Nll == -981534.372474
f2 == 0.050000 ---- New Nll == -981589.252634
f2 == 0.075000 ---- New Nll == -981632.945147
f2 == 0.100000 ---- New Nll == -981666.622045
f2 == 0.125000 ---- New Nll == -981691.210110
f2 == 0.150000 ---- New Nll == -981707.458049
f2 == 0.175000 ---- New Nll == -981715.981441
f2 == 0.200000 ---- New Nll == -981717.293867
f2 == 0.225000 ---- New Nll == -981711.829112
f2 == 0.250000 ---- New Nll == -981699.957393
f2 == 0.275000 ---- New Nll == -981681.997477
f2 == 0.300000 ---- New Nll == -981658.225911
f2 == 0.325000 ---- New Nll == -981628.884150
f2 == 0.350000 ---- New Nll == -981594.184157
f2 == 0.375000 ---- New Nll == -981554.312860
f2 == 0.400000 ---- New Nll == -981509.435728
f2 == 0.425000 ---- New Nll == -981459.699685
f2 == 0.450000 ---- New Nll == -981405.235497
f2 == 0.475000 ---- New Nll == -981346.159748

12/21/2006

Test the RooDLineShape

cd /home/xs32/local/app
cp -rp  /nfs/cor/user/ponyisi/daf9/RooFitCore .
cp -rp  /nfs/cor/user/ponyisi/daf9/RooFitModels .
cd RooFitCore
gmake -f GNUmakefile.standalone clean
gmake -f GNUmakefile.standalone

cd RooFitModels
gmake -f GNUmakefile.standalone clean
gmake -f GNUmakefile.standalone

Add print message in RooFitModels/RooDLineShapeImp.cc
Then recompile the code
Make the code to load file from the app/ area

Send out a test job DHad.py fit test -m 0 —qsub job 73806

Log file

12/22/2006

Make the fitting faster

List the available fitting parameters

To make the fitting debug faster, fix most of the parameters, only left N1, N2, fa:

DHad.py fit fa fixALL plotNll -t d -m 0

Get the Nll+/-0.5 value and plot:

DHad.py fit fa  plotNll -t d -m 0

Mode 0: 'kpi'

12/26/2006

calculate the Nll + 0.5 value

Test

DHad.py fit fa fixALL plotNll -t d -m 0

OK.

DHad.py fit fa plotNll -t d -m 0

Mode 0: 'kpi'

Get the f2 value: 0.192 +/- 0.010

Compare with the value from Signal MC: 0.182 +/- 0.045, it's in side the one sigma area.

Compare the difference of fa=0.182 and fa=0.202

fa = 0.182 happened to be the normal value. DHad.py fit fa 0.182 -t d -m 0 —qsub job 75943

DHad.py fit fa 0.202 -t d -m 0 —qsub job 75944

For mode 1

DHad.py fit fa plotNll -t d -m 1 —qsub job 75946

12/29/2006

Check Mode 1

From the logfile (75946), !!! User command too long (1024>1019)!!!

Remove the longer .C pic name for a moment. Then redo the job

DHad.py fit fa plotNll -t d -m 1 —qsub job 77440 done.

Log file

STATUS=NOT POSDEF

Continue the Mode 0

Check the previous jobs, normal done.

Layout the plots

Mode 0: 'kpi'

Compare fa 202 with normal fa

Before doing that, submit the fit for the S/G MC first:

DHad.py fit fa 0.202 -t s -m 0 —qsub job 77441 done

DHad.py fit fa 0.202 -t g -m 0 —qsub job 77442 done

Then coding for the compare:

ref: DHadTable.py -a compare -o yields_faplus_cleo200609

DHad.py table compare yields fa/0.202 cleo200609  -m 0

Compare yields between fa 0.202 and cleo200609

[ Front | MyDHad-archives | Index ]

Updated January 16, 2007

D Had Log - November 2006

Home > D-Hadronic > Archives

Archives

Links

12/03/2006

plot the kpi, kpipi mode to see the logL distribution of fa

Consider using the Roofit inner function to finish this job

12/04/2006

Using the exsiting one and make the plot

12/05/2006

plot fa vs. logL

Consider using the Roofit inner function for Data

pi0 mass in Kpipi0 for Data

12/06/2006

plot fa vs. logL

Hit fraction shrink to 60%(70,80,90) in the Ntuple DT/MC

12/07/2006

fa loglikelihood

Organize the previous fit

12/15/2006

plot fa vs. logL

12/21/2006

Test the RooDLineShape

12/22/2006

Make the fitting faster

12/26/2006

calculate the Nll + 0.5 value

12/29/2006

Check Mode 1

Continue the Mode 0