Hi Sheila,
is the solution of general interest?
Heiko
On Thu, 30 Mar 2006, Sheila Mclachlin wrote:
>
> Hello,
>
> Just FYI:
>
> I found the problem, so I no longer need help
> with this.
>
> sheila
>
> ---------- Forwarded message ----------
> Date: Tue, 28 Mar 2006 18:27:37 -0500 (EST)
> From: Sheila Mclachlin <[log in to unmask]>
> To: vub-recoil <[log in to unmask]>
> Subject: ntuple help
>
>
> Hello!
>
> I am trying to solve a problem in my ntuple-making
> code, and I thought maybe someone could help.
>
> The main code is in:
>
> ~penguin/vubrecoil/vxb30/VxbRecoil/VxbXlnu.cc
>
> (VxbRecoil is a package I made, based on
> VubRecoilUser and VcbRecoil.)
>
> Here is the problem: I run 30 events with the command:
>
> bin/$BFARCH/VxbRecoilApp VxbXlnu.tcl SP-1235-BSemiExcl-R18b-2 MC 1235 2
>
> which produces the following output files:
>
> workdir/log/SP-1235-BSemiExcl-R18b-2.log
> workdir/log/SP-1235-BSemiExcl-R18b-2.root
>
> Then I look at my ntuple in ROOT.
> Here is the confusing part. For example:
>
> root [0] TFile f("data/SP-1235-BSemiExcl-R18b-2.root");
> root [1] ntp1->Scan("eUps:eBlab:decayMode");
> ************************************************
> * Row * eUps * eBlab * decayMode *
> ************************************************
> * 0 * 12.101385 * 6.0293622 * 11101 *
> * 1 * 12.101385 * -10 * 0 *
> * 2 * 12.101385 * -10 * 0 *
> * 3 * 12.101385 * 6.3097918 * 12201 *
> * 4 * 12.101385 * -10 * 0 *
> * 5 * 12.101448 * 5.9822311 * 15123 *
> * 6 * 12.101448 * 6.0727898 * 15236 *
> * 7 * 12.101371 * 6.3998302 * 11236 *
> * 8 * 12.101428 * 6.1657033 * 14203 *
> * 9 * 12.101428 * 6.3021577 * 14323 *
> * 10 * 12.101420 * 5.7747154 * 14141 *
> * 11 * 12.101420 * -10 * 0 *
> * 12 * 12.101420 * 5.8976044 * 11109 *
> * 13 * 12.101369 * -10 * 0 *
> * 14 * 12.101369 * 6.1250192 * 14322 *
>
> The first column has a value of eUps in every row.
> However, for some reason, the values of eBlab have
> been filled with their default values (-10.0) in
> some of the rows.
>
> There are exactly 30 rows (ie, one per event),
> and each has exactly one B meson (I checked).
>
> According to HepTuple/Tuple.h:
>
> "On a particular "event", if no call is made to "column"
> for a particular column, that column's
> default value is used when filling the ntuple."
>
> This tells me that for some reason, in some events
> there is no call made to "column". But as far
> as I can see, there is nothing in VxbXlnu.cc that would
> cause that to happen. The dumpBreco function should
> be called every time the dumpEventInfo function is
> called.
>
> I checked the Breco pointers, and they were all nonzero.
>
> Does anyone know what is causing this?
>
> Thanks,
>
> sheila
>
>
>
>
>
|