D Hadronic Analysis - Log

03/03/2006

03/06/2006

03/07/2006

03/08/2006

03/10/2006

03/14/2006

03/15/2006

Find out the blocks for MC
Creat the beam constrained mass and output the ascii file

03/23/2006

Search through the Monte Carlo particle list
Start the systematic analysis of delta E

03/24/2006

Search through the Monte Carlo particle list
Systematic analysis of delta E
Create yield table

03/27/2006

Creat script to read the numbers from eps

03/28/2006

Calculate the difference

03/29/2006

script handle-table.py

03/30/2006

03/03/2006

Move the old fig to 20060302

Select the events based on the new ntuple.

Check the fitpars.C and TranslateFitParms.C before fit

/cdat/daf9/ponyisi/hadD/summerconf/fits/fitpars.C .  You may be able to
directly use the fitpars_C.so file in that directory.

link the fitpars_C.so.

Need to change the fit code to refresh the evt file!

03/06/2006

change the DHadFit code to refresh the evt file.

Signal MC

c3python DHad.py -t s -m 0 --tag s --sign 1
c3python DHad.py -t s -m 0 --tag s --sign -1
c3python DHadFit.py -t s -m 0 --tag s


c3python DHad.py -t s -m 1 --tag s --sign 1
c3python DHad.py -t s -m 1 --tag s --sign -1
c3python DHadFit.py -t s -m 1 --tag s

c3python DHad.py -t s -m 3 --tag s --sign 1
c3python DHad.py -t s -m 3 --tag s --sign -1
c3python DHadFit.py -t s -m 3 --tag s

c3python DHad.py -t s -m 200 --tag s --sign 1
c3python DHad.py -t s -m 200 --tag s --sign -1
c3python DHadFit.py -t s -m 200 --tag s

c3python DHad.py -t s -m 201 --tag s --sign 1
c3python DHad.py -t s -m 201 --tag s --sign -1
c3python DHadFit.py -t s -m 201 --tag s

c3python DHad.py -t s -m 202 --tag s --sign 1
c3python DHad.py -t s -m 202 --tag s --sign -1
c3python DHadFit.py -t s -m 202 --tag s


c3python DHad.py -t s -m 203 --tag s --sign 1
c3python DHad.py -t s -m 203 --tag s --sign 1
c3python DHadFit.py -t s -m 203 --tag s

c3python DHad.py -t s -m 204 --tag s --sign 1
c3python DHad.py -t s -m 204 --tag s --sign -1
c3python DHadFit.py -t s -m 204 --tag s

c3python DHad.py -t s -m 205 --tag s --sign 1
c3python DHad.py -t s -m 205 --tag s --sign -1
c3python DHadFit.py -t s -m 205 --tag s

Generic MC

c3python DHad.py -t g -m 0 --tag s --sign 1

events : 0

check the DHad.py. The data format of the root file is now:

"data31_mc_20051121_mc-ddbar-generic_10xlumi_200978_201345.root",

So, chage:

DD-> ddbar

c3python DHad.py -t g -m 0 --tag s --sign -1
c3python DHadFit.py -t g -m 0 --tag s


c3python DHad.py -t g -m 1 --tag s --sign 1
c3python DHad.py -t g -m 1 --tag s --sign -1
c3python DHadFit.py -t g -m 1 --tag s


c3python DHad.py -t g -m 3 --tag s --sign 1
c3python DHad.py -t g -m 3 --tag s --sign -1
c3python DHadFit.py -t g -m 3 --tag s

c3python DHad.py -t g -m 200 --tag s --sign 1
c3python DHad.py -t g -m 200 --tag s --sign -1
c3python DHadFit.py -t g -m 200 --tag s

c3python DHad.py -t g -m 201 --tag s --sign 1
c3python DHad.py -t g -m 201 --tag s --sign -1
c3python DHadFit.py -t g -m 201 --tag s

c3python DHad.py -t g -m 202 --tag s --sign 1
c3python DHad.py -t g -m 202 --tag s --sign -1
c3python DHadFit.py -t g -m 202 --tag s

c3python DHad.py -t g -m 203 --tag s --sign 1
c3python DHad.py -t g -m 203 --tag s --sign -1
c3python DHadFit.py -t g -m 203 --tag s

c3python DHad.py -t g -m 204 --tag s --sign 1
c3python DHad.py -t g -m 204 --tag s --sign -1
c3python DHadFit.py -t g -m 204 --tag s

c3python DHad.py -t g -m 205 --tag s --sign 1
c3python DHad.py -t g -m 205 --tag s --sign -1
c3python DHadFit.py -t g -m 205 --tag s

03/07/2006

Questions:

There are messages after each fit.

 *** Break *** segmentation violation
 Generating stack trace...
/usr/bin/addr2line: python: No such file or directory
 0xb6f1259d in THashList::~THashList() + 0x2f from /nfs/cleo3/Offline/rel/20050417_FULL/other_sources/Root/lib/libCore.so
 0xb55a3b1c in <unknown> from /nfs/cleo3/Offline/other_sources/pkg/Root/Root-4.03.02/root/install/lib/libHtml.so
 0xb78c76c3 in exit + 0x63 from /lib/tls/libc.so.6
 0xb78b27ec in <unknown> from /lib/tls/libc.so.6
 0x08048625 in Py_Main + 0x31 from python
/home/xs32/local/bin/c3python: line 19:  7121 Aborted                 python $@
x

Read out the masses

Mode	md (Signal MC)	md (Generic MC)
0 : 'kpi'	1.8645941 +/- 0.0000077	1.8646088 +/- 0.0000022
1 : 'kpipi0'	1.8646603 +/- 0.0000056	1.8646645 +/- 0.0000022
3 : 'kpipipi'	1.8645991 +/- 0.0000052	1.8646057 +/- 0.0000019
200: 'kpipi'	1.8694081 +/- 0.0000052	1.8694070 +/- 0.0000018
201: 'kpipipi0'	1.8694413 +/- 0.0000094	1.8694266 +/- 0.0000039
202: 'kspi'	1.8694112 +/- 0.0000057	1.8694118 +/- 0.0000049
203: 'kspipi0'	1.869604 +/- 0.000012	1.8694323 +/- 0.0000049
204: 'kspipipi'	1.869388 +/- 0.000010	1.8693989 +/- 0.0000048
205: 'kkpi'	1.869420 +/- 0.000010	1.8694064 +/- 0.0000052

Look up the $C3_DATA/evt.pdl

evt.pdl

 *      type        name     id  mass/GeV width/GeV max_Dm/GeV 3*charge 2*spin lifetime*c/mm Lund-KC


add  p Meson       D+      411  1.8693          0       0        3       0  0.315 104
add  p Meson       D-     -411  1.8693          0       0       -3       0  0.315  0
add  p Meson       D0      421  1.8645          0       0        0       0  0.1234 105

Mass difference M(fit)- M(input)

Mode	md (Signal MC)	md (Generic MC)
0 : 'kpi'	0.0000941 +/- 0.0000077	0.0001088 +/- 0.0000022
1 : 'kpipi0'	0.0001603 +/- 0.0000056	0.0001645 +/- 0.0000022
3 : 'kpipipi'	0.0000991 +/- 0.0000052	0.0001057 +/- 0.0000019
200: 'kpipi' \|-0.0002081 +/- 0.0000052	0.0001070 +/- 0.0000018
201: 'kpipipi0'	0.0001413 +/- 0.0000094	0.0001266 +/- 0.0000039
202: 'kspi'	0.0001112 +/- 0.0000057	0.0001118 +/- 0.0000049
203: 'kspipi0'	0.000304 +/- 0.000012	0.0001323 +/- 0.0000049
204: 'kspipipi'	0.000088 +/- 0.000010	0.0000989 +/- 0.0000048
205: 'kkpi'	0.000220 +/- 0.000010	0.0001064 +/- 0.0000052

03/08/2006

Correct the typos

Mass difference M(fit)- M(input)

Mode	md (Signal MC)	md (Generic MC)
0 : 'kpi'	0.0000941 +/- 0.0000077	0.0001088 +/- 0.0000022
1 : 'kpipi0'	0.0001603 +/- 0.0000056	0.0001645 +/- 0.0000022
3 : 'kpipipi'	0.0000991 +/- 0.0000052	0.0001057 +/- 0.0000019
200: 'kpipi'	0.0001081 +/- 0.0000052	0.0001070 +/- 0.0000018
201: 'kpipipi0'	0.0001413 +/- 0.0000094	0.0001266 +/- 0.0000039
202: 'kspi'	0.0001112 +/- 0.0000057	0.0001118 +/- 0.0000049
203: 'kspipi0'	0.000304 +/- 0.000012	0.0001323 +/- 0.0000049
204: 'kspipipi'	0.000088 +/- 0.000010	0.0000989 +/- 0.0000048
205: 'kkpi'	0.000120 +/- 0.000010	0.0001064 +/- 0.0000052

Friday March 10, 2006 - 09:44

tuple.Draw("ecm-mce[0]")

Start from DHad.py

root [] TFile f("Single_Dp_to_Kspi.root")
root [] f.ls()
  OBJ: TTree    dnt     Hadronic D Ntuple : 0 at: 0x90c2990

root [] dnt->Draw("ecm-mce[0]")

c3python DHad.py -a 1 -t s -m 0 --tag s --sign 1

Error message:

Traceback (most recent call last):
  File "DHad.py", line 326, in ?
    main()
  File "DHad.py", line 320, in main
    Evt.drawTree("ecm")
  File "DHad.py", line 172, in drawTree
    self.pt.Draw(x)
AttributeError: 'PyTTree' object has no attribute 'Draw'

Add method at PyROOTExt.py

class PyTTree(object):
    def Draw(self, entry):
        self.rtree.Draw(entry)

Error:

Error in <TTreeFormula::DefinedVariable>: the branch "ecm" has to be
enabled to be used

 ERROR 31 :
 The Variable :  "ecm" exists but is not accessible

Message from Peter:

It does not work because those branches are disabled.  If you want to use
TTree::Draw in python, do _not_ do it from a PyTTree object, in fact never
make the PyTTree object.

03/10/2006

tuple.Draw("ecm-mce[0]")

Refine DHad.py

c3python DHad.py -a 1 -t s --tag s

c3python DHad.py -a 1 -t g --tag s

Pass the test.

03/14/2006

Several pages need to be bookmarked:

CLEO III Software Cross-Reference

Pages for the TUPLE_DEFINITION

03/15/2006

Find out the blocks for MC

In the Pages for the TUPLE_DEFINITION, we find:

424 BLOCK MC OPTIONAL
425 INTEGER MCDMODE
426 DOC "Decay mode of D0 or D+.  NOTE: *NOT* the same as dmode!!!"
427 INTEGER MCDBMODE
428 DOC "Decay mode of D0bar or D-.  NOTE: *NOT* the same as dmode!!!"
429 INTEGER MCDMATCH NDCAND
430 DOC "Index to MC Block if the DTag has a MC match"
431 INDEX NMCPART 0 50
432 DOC "Number of Monte Carlo particles in this block"
433 INTEGER MCIDENT NMCPART
434 DOC "Internal node identifier"
435 INTEGER MCPDGID NMCPART
436 DOC "PDG ID of particle"
437 INTEGER MCPARENT NMCPART
438 DOC "Index of parent in block"
439 FLOAT MCE NMCPART
440 DOC "Energy of MC particle"
441 FLOAT MCPX NMCPART
442 DOC "Px of MC particle"
443 FLOAT MCPY NMCPART
444 DOC "Py of MC particle"
445 FLOAT MCPZ NMCPART
446 DOC "Pz of MC particle"
447 INTEGER MCTRUTH NMCPART
448 DOC "Index of the reconstructed object corresponding to this particle"

Creat the beam constrained mass and output the ascii file

Look at the DHad.py

Need:

1. pte.dmbc[d]

2. pte.ecm/2

3. pte.dcsign[d]

Look at the root file

root []:

TFile f("/nfs/cor/an1/ponyisi/dhad/summerconf/signal/dtuple_c_0525_p2_1104/Single_Dp_to_Kspi.root")
f.ls()
dnt->Show(10)

...

mcpdgid         = 91,
                  30443, 22, 411, -411, -311, 211, 310, 111, 111, 22,
                  22, 22, 22, 311, -211, 111, 310, 211, -211
 mcparent        = -1,
                  0, 0, 1, 1, 3, 3, 5, 7, 7, 8,
                  8, 9, 9, 4, 4, 4, 14, 17, 17
 mce             = 3.77045,
                  3.76173, 0.0087195, 1.88047, 1.88126, 0.912706,
                  0.967767, 0.912706, 0.193089, 0.719617, 0.14091,
                  0.0521915, 0.625231, 0.0943865, 0.882455, 0.286929,
                  0.71187, 0.882455, 0.257815, 0.624641
 mcpx            = -0.0118842,
                  -0.0118572, -3.52188e-05, 0.17902, -0.190884, -0.478865,
                  0.657878, -0.478865, -0.114716, -0.364157, -0.104544,
                  -0.0101802, -0.306937, -0.0572272, -0.601427, -0.161637,
                  0.572164, -0.601427, -0.0281432, -0.573291
 mcpy            = 0,
                  0, 0, -0.0186152, 0.0185998, 0.130876,
                  -0.149499, 0.130876, 0.0739338, 0.0569348, 0.0903948,
                  -0.0164692, 0.00993078, 0.0469968, 0.0131778, -0.142675,
                  0.148082, 0.0131778, 0.0486508, -0.0354802
 mcpz            = 0.000997781,
                  -0.00772122, 0.00871178, -0.0974952, 0.0897658, 0.582163,
                  -0.679665, 0.582163, -0.0209822, 0.603137, 0.0274768,
                  -0.0484662, 0.544611, 0.0585178, -0.411315, 0.127939,
                  0.373128, -0.411315, -0.209364, -0.201959
 mctruth         = -1,
                  -1, -1, -1, -1, -1, 0, -1, -1, 0, 6,
                  8, 0, 7, -1, 1, 1, -1, -1, 2

Find out the defination of the beam constrained mass

In Peter's talk on CLEO101:

Two pages explaining this: [ 1 | 2 ]

The equation for m_BC is :

dnt->Draw("ecm**2-(mcpx[0]**2+mcpy[0]**2+mcpz[0]**2)")

Message from Peter

To first order, it should be
"sqrt((ecm/2)**2-(mcpx[0]**2+mcpy[0]**2+mcpz[0]**2))"
(there are small corrections due to crossing angle)

Apply it in the DHad.py to output ascii file

Ask Peter :

        for pte in self.pt:
            if pte.ecm < 3.7:
                continue
            d = PueoUtils.chooseD(self.mode, pte, self.sign)
            if d != None:
                fo.write('%.5f %.5f %d\n' % \
   (sqrt((pte.ecm/2)**2-(pte.mcpx[0]**2+pte.mcpy[0]**2+pte.mcpz[0]**2)),\
    pte.ecm/2, pte.dcsign[d]))

 or [0] -> [d] ?

Go on to test.

c3python DHad.py --test -a 2  -t s -m 0 --tag s --sign 1

Reply from Peter:

You do not want to choose the D.  The whole point is to avoid reconstruction.
You want to search through the Monte Carlo particle list -- not the
reconstructed particle list.

03/23/2006

Search through the Monte Carlo particle list

Dig in the chooseD in PueoUtils.py

Ask Peter for how to search MC list.

Start the systematic analysis of delta E

Formula

      ( data yield wide Delta E)          (MC yield wide Delta E)
 ( ------------------------------- ) * ( -------------------------- ) ^ -1
     ( data yiels regular Delta E)       (MC yield regular Delta E)

Here MC stands for signal MC.

Location of the evt file

/cdat/daf9/ponyisi/hadD/summerconf/txt/

03/24/2006

Search through the Monte Carlo particle list

Reply from Peter

You do not care about which mode it is, because you aren't reconstructing
it from its decay products.  There is only one D for each sign in the MC
list.

"one D for each sign in the MC list." Not quite understand. Let's try to dump out the selected D+ by pid first.

 D+  pid   411
 D-  pid  -411
 D0  pid   421

In the DHad.py:

In the 'chooseD' of PueoUtils.py, the returned value of d is pte.ndcand. Which means: "Index to MC Block if the DTag has a MC match" refer to Tuple_Definition Page.

for i in range(pte.ndcand):

Try d=0 at this time, and look at the output.

c3python DHad.py --test -a 2  -t s --tag s

Can not use the "mcpdgid" as the selection criteria.

Suggestion from Peter: Move on to the systematic delta E first.

Systematic analysis of delta E

Signal MC

Mode 0 :'kpi'

cdfit
c3python DHadFit.py -a 1 -t s -m 0 --tag s

Mode 1 : 'kpipi0'

c3python DHadFit.py -a 1 -t s -m 1 --tag s

Mode 3 : 'kpipipi'

c3python DHadFit.py -a 1 -t s -m 3 --tag s

Mode 200: 'kpipi'

c3python DHadFit.py -a 1 -t s -m 200 --tag s

Mode 201: 'kpipipi0'

c3python DHadFit.py -a 1 -t s -m 201 --tag s

Mode 202: 'kspi'

c3python DHadFit.py -a 1 -t s -m 202 --tag s

Mode 203: 'kspipi0'

c3python DHadFit.py -a 1 -t s -m 203 --tag s

Mode 204: 'kspipipi'

c3python DHadFit.py -a 1 -t s -m 204 --tag s

Mode 205: 'kkpi'

c3python DHadFit.py -a 1 -t s -m 205 --tag s

Data

Mode 0 :'kpi'

c3python DHadFit.py -a 1 -t d -m 0 --tag s

Mode 1 : 'kpipi0'

c3python DHadFit.py -a 1 -t d -m 1 --tag s

Mode 3 : 'kpipipi'

c3python DHadFit.py -a 1 -t d -m 3 --tag s

Mode 200: 'kpipi'

c3python DHadFit.py -a 1 -t d -m 200 --tag s

Mode 201: 'kpipipi0'

c3python DHadFit.py -a 1 -t d -m 201 --tag s

Mode 202: 'kspi'

c3python DHadFit.py -a 1 -t d -m 202 --tag s

Mode 203: 'kspipi0'

c3python DHadFit.py -a 1 -t d -m 203 --tag s

Mode 204: 'kspipipi'

c3python DHadFit.py -a 1 -t d -m 204 --tag s

Mode 205: 'kkpi'

c3python DHadFit.py -a 1 -t d -m 205 --tag s

Regular Data Fit

Mode 0 :'kpi'

c3python DHadFit.py -t d -m 0 --tag s

Mode 1 : 'kpipi0'

c3python DHadFit.py -t d -m 1 --tag s

Mode 3 : 'kpipipi'

c3python DHadFit.py  -t d -m 3 --tag s

Mode 200: 'kpipi'

c3python DHadFit.py -t d -m 200 --tag s

Mode 201: 'kpipipi0'

c3python DHadFit.py  -t d -m 201 --tag s

Mode 202: 'kspi'

c3python DHadFit.py  -t d -m 202 --tag s

Mode 203: 'kspipi0'

c3python DHadFit.py  -t d -m 203 --tag s

Mode 204: 'kspipipi'

c3python DHadFit.py -t d -m 204 --tag s

Mode 205: 'kkpi'

c3python DHadFit.py -t d -m 205 --tag s

Create yield table

Signal MC

Mode	MC regular -N1	N2	MC wide N1	N2
0	20460	20809	20901	21217
1	71301	71668	70205	71368
3	46329	46903	47619	47735
200	43972	44131	44726	45168
201	20965	21005	20623	21207
202	9076	9118	9341	9190
203	13600	13551	13379	13409
204	12920	12847	12785	12646
205	9495	9398	9512	9509

Data

Create a script to read out the number from eps file

03/27/2006

Creat script to read the numbers from eps

python readeps.py -f ../../doc/fig/data_Single_D0AndD0B_to_Kpi.eps

N1 = 25757
N2 = 26254

python readeps.py -a 1 -t d --tag s


------   Processing more files. ------
mode     N1      N2
0        25757   26254
1        50351   50621
3        39782   39682
200      40188   40677
201      13071   12980
202      5770    5848
203      13332   13184
204      8281    8290
205      3527    3508

python readeps.py -a 1 -t wd --tag s

Mode	data N1	data N2	wide N1	wide N2
0	25757	26254	26187	26743
1	50351	50621	51082	51195
3	39782	39682	40790	40769
200	40188	40677	41162	41565
201	13071	12980	13114	13217
202	5770	5848	5877	6012
203	13332	13184	13469	13511
204	8281	8290	8372	8289
205	3527	3508	3622	3599

Creat the MC table

python readeps.py -a 1 -t s --tag s
python readeps.py -a 1 -t ws --tag s

Mode	MC N1	MC N2	wide N1	wide N2
0	20617	20950	20901	21217
1	69752	70842	70205	71368
3	46661	46771	47619	47735
200	43894	44276	44726	45168
201	20766	20814	20623	21207
202	36537	36171	9341	9190
203	13304	13322	13379	13409
204	12582	12440	12785	12646
205	9399	9409	9512	9509

CBX 05-08
Redo the fit for signal MC, using Peter's Before that, mv the original figure to fig/20060327/
Change the command line as "-t ws" for "wide". So, now -a 1 only stands for Peter's evt file locations.

cdfit
c3python DHadFit.py -a 1 -t s -m 0 --tag s

c3python DHadFit.py -a 1 -t s -m 1 --tag s

c3python DHadFit.py -a 1 -t s -m 3 --tag s

c3python DHadFit.py -a 1 -t s -m 200 --tag s

c3python DHadFit.py -a 1 -t s -m 201 --tag s

c3python DHadFit.py -a 1 -t s -m 202 --tag s

c3python DHadFit.py -a 1 -t s -m 203 --tag s

c3python DHadFit.py -a 1 -t s -m 204 --tag s

c3python DHadFit.py -a 1 -t s -m 205 --tag s

03/28/2006

Update the new yield list for MC

python readeps.py -a 1 -t s --tag s
python readeps.py -a 1 -t ws --tag s

Still get the same table.

202     | 36537 | 36171

Compare the original evt file with Peters:

DpAndDm_to_Kspi

/cdat/daf9/ponyisi/hadD/summerconf/txt/signal_Single_Dp_to_Kspi.evt

dat/evt/signal_Single_Dp_to_Kspi.evt

They are the same file.

So, there might be something strange in the widede file.

>| Pause here. (Peter regenerate the evt file later after the meeting)

Calculate the difference

Save the data talbe below as dat/table/deltae.txt

Mode	data N1	data N2	wide N1	wide N2	MC N1	MC N2	wide N1	wide N2
0	25757	26254	26187	26743	20617	20950	20901	21217
1	50351	50621	51082	51195	69752	70842	70205	71368
3	39782	39682	40790	40769	46661	46771	47619	47735
200	40188	40677	41162	41565	43894	44276	44726	45168
201	13071	12980	13114	13217	20766	20814	20623	21207
202	5770	5848	5877	6012	36537	36171	9341	9190
203	13332	13184	13469	13511	13304	13322	13379	13409
204	8281	8290	8372	8289	12582	12440	12785	12646
205	3527	3508	3622	3599	9399	9409	9512	9509

Create a script to get the

handle-table.py

fix "p"

03/29/2006

script handle-table.py

03/30/2006

Fit regular MC:

cdfit
c3python DHadFit.py -a 1 -t s -m 0 --tag s
c3python DHadFit.py -a 1 -t s -m 1 --tag s
c3python DHadFit.py -a 1 -t s -m 3 --tag s
c3python DHadFit.py -a 1 -t s -m 200 --tag s
c3python DHadFit.py -a 1 -t s -m 201 --tag s
c3python DHadFit.py -a 1 -t s -m 202 --tag s
c3python DHadFit.py -a 1 -t s -m 203 --tag s
c3python DHadFit.py -a 1 -t s -m 204 --tag s
c3python DHadFit.py -a 1 -t s -m 205 --tag s

Continue work on script ... done
Read out the siganl MC yield table

python readeps.py -a 1 -t s --tag s

Same as before except one mode 203:

before	13304	13322
after	13270	13269

Fit widede signal MC:

cdfit
c3python DHadFit.py -a 1 -t ws -m 0 --tag s
c3python DHadFit.py -a 1 -t ws -m 1 --tag s
c3python DHadFit.py -a 1 -t ws -m 3 --tag s
c3python DHadFit.py -a 1 -t ws -m 200 --tag s
c3python DHadFit.py -a 1 -t ws -m 201 --tag s
c3python DHadFit.py -a 1 -t ws -m 202 --tag s
c3python DHadFit.py -a 1 -t ws -m 203 --tag s
c3python DHadFit.py -a 1 -t ws -m 204 --tag s
c3python DHadFit.py -a 1 -t ws -m 205 --tag s

Read out the widede siganl MC yield table

python readeps.py -a 1 -t ws --tag s

mode	N1	N2
0	20901	21217
1	70205	71368
3	47619	47735
200	44726	45168

Two modes changed

mode ||N1 ||N2 ||N1 before||N2 before

201	20969	21053	20623	21207
202	37074	36704	9341	9190
203	13379	13409
204	12785	12646
205	9512	9509

Update the deltae.txt:

Signal MC

203  | 13304  | 13322   -> | 13270 | 13269
201  | 20623  | 21207   -> | 20969 | 21053
202  | 9341   | 9190    -> | 37074 | 36704

Mode	data N1	data_N2	wide_N1	wide N2	MC N1	MC N2\|\| wide N1	wide N2
0	25757 \|26254\| 26187	26743	20617	20950	20901	21217
1	50351 \|50621\| 51082	51195	69752	70842	70205	71368
3	39782 \|39682\| 40790	40769	46661	46771	47619	47735
200	40188 \|40677\| 41162	41565	43894	44276	44726	45168
201	13071 \|12980\| 13114	13217	20766	20814	20969	21053
202	5770 \|5848	5877	6012	36537	36171	37074	36704
203	13332 \|13184\| 13469	13511	13270	13269	13379	13409
204	8281 \|8290	8372	8289	12582	12440	12785	12646
205	3527 \|3508	3622	3599	9399	9409	9512	9509

Calculate the values

python handle-table.py -f ../../dat/table/deltae.txt

Mode	N1	N2
0	1.003	1.006
1	1.008	1.004
3	1.005	1.007
200	1.005	1.002
201	0.994	1.007
202	1.004	1.013
203	1.002	1.014
204	0.995	0.984
205	1.015	1.015

Fit the fixed "p" value Look at the .c file

/cdat/daf9/ponyisi/hadD/summerconf/fits/mbc_singletag_3s.c

Variable 'double p_side', the No.10th parameter in DHadFit.py

('', '', evtfile, epsfile, outfile, mc, fp.getSigma(), fp.getSa(),
fp.getSb(), fp.getFa(), fp.getFb(), optstr)

c3python DHadFit.py -a 1 -t s -m 0 --tag s -p 0.5

[ Front | MyDHad-archives | Index ]

Updated January 11, 2007

D Hadronic Analysis - Log - March 2006

Home > D-Hadronic > Archives

Archives

Links

03/03/2006

03/06/2006

03/07/2006

Read out the masses

Look up the $C3_DATA/evt.pdl

03/08/2006

03/10/2006

03/14/2006

03/15/2006

Find out the blocks for MC

Creat the beam constrained mass and output the ascii file

Look at the DHad.py

Look at the root file

Find out the defination of the beam constrained mass

Apply it in the DHad.py to output ascii file

03/23/2006

Search through the Monte Carlo particle list

Start the systematic analysis of delta E

03/24/2006

Search through the Monte Carlo particle list

Systematic analysis of delta E

Signal MC

Mode 0 :'kpi'

Mode 1 : 'kpipi0'

Mode 3 : 'kpipipi'

Mode 200: 'kpipi'

Mode 201: 'kpipipi0'

Mode 202: 'kspi'

Mode 203: 'kspipi0'

Mode 204: 'kspipipi'

Mode 205: 'kkpi'

Data

Mode 0 :'kpi'

Mode 1 : 'kpipi0'

Mode 3 : 'kpipipi'

Mode 200: 'kpipi'

Mode 201: 'kpipipi0'

Mode 202: 'kspi'

Mode 203: 'kspipi0'

Mode 204: 'kspipipi'

Mode 205: 'kkpi'

Regular Data Fit

Mode 0 :'kpi'

Mode 1 : 'kpipi0'

Mode 3 : 'kpipipi'

Mode 200: 'kpipi'

Mode 201: 'kpipipi0'

Mode 202: 'kspi'

Mode 203: 'kspipi0'

Mode 204: 'kspipipi'

Mode 205: 'kkpi'

Create yield table

03/27/2006

Creat script to read the numbers from eps

03/28/2006

Calculate the difference

03/29/2006

script handle-table.py

03/30/2006