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')
>
>
|