VUB-RECOIL Archives

Vub measurement using recoil of fully reconstructed Bs

VUB-RECOIL@LISTSERV.SLAC.STANFORD.EDU

Options: Use Forum View

Use Monospaced Font
Show Text Part by Default
Show All Mail Headers

Message: [<< First] [< Prev] [Next >] [Last >>]
Topic: [<< First] [< Prev] [Next >] [Last >>]
Author: [<< First] [< Prev] [Next >] [Last >>]

Print Reply
Subject:
From:
Daniele del Re <[log in to unmask]>
Date:
08 Oct 2002 15:10:54 -0700 (PDT)Tue, 08 Oct 2002 15:10:54 -0700 (PDT)
Content-Type:
TEXT/PLAIN
Parts/Attachments:
TEXT/PLAIN (57 lines)

Yes, this is the wrong piece of code.

You should look at the method Fitmessub() and not Fitmes() in fitUtil.icc

Daniele


> Hi all,
> After a very quick look at the code in fitNtp.cc I see that the b0 mixing
> correction is still done at mes level (maybe I'm wrong but I see that
> piece of code in :
> theFit() :
>
>      for(int k=1;k<40;k++){
>        tempbinchb = ((TH1D*)gDirectory->Get(namebch))->GetBinContent(k);
>        tempbinb0os = ((TH1D*)gDirectory->Get(nameb0os))->GetBinContent(k);
>        tempbinb0ss = ((TH1D*)gDirectory->Get(nameb0ss))->GetBinContent(k);
>        temperrchb = ((TH1D*)gDirectory->Get(namebch))->GetBinError(k);
>        temperrb0os = ((TH1D*)gDirectory->Get(nameb0os))->GetBinError(k);
>        temperrb0ss = ((TH1D*)gDirectory->Get(nameb0ss))->GetBinError(k);
>        tempbin = tempbinchb + ((1-chid)/(1-2*chid)) * tempbinb0os -
> (chid/(1-2*chid)) * tempbinb0ss;
>        if(MIXCORR==0){
> 	 tempbin = tempbinchb + tempbinb0os;
> 	 temperr = sqrt(temperrchb * temperrchb + temperrb0os *
> temperrb0os);
>        }else{
> 	 if(tempbin<0) tempbin=0;
> 	 temperr = sqrt(temperrchb*temperrchb + ((1-chid)/(1-2*chid)) *
> ((1-chid)/(1-2*chid)) * temperrb0os * temperrb0os + (chid/(1-2*chid)) *
> (chid/(1-2*chid)) * temperrb0ss * temperrb0ss);
>        }
>        ((TH1D*)gDirectory->Get(name))->SetBinContent(k, tempbin);
>        ((TH1D*)gDirectory->Get(name))->SetBinError(k, temperr);
>     }
>
> This is the same code as that i'm using right now to make comparisons...
>
> Maybe I'm lloking at the wrong piece of code..
> Can please Daniele point me to the right one?
> Alessio
>
> ______________________________________________________
> Alessio Sarti     Universita' & I.N.F.N. Ferrara
>  tel  +39-0532-781928  Ferrara
> roma  +39-06-49914338
> SLAC +001-650-926-2972
>
> "... e a un Dio 'fatti il culo' non credere mai..."
> (F. De Andre')
>
>




ATOM RSS1 RSS2