Hi,
the weights that can be used for this scan are in
/u/ec/rfaccini/dev/gen12.4.0g/workdir/weights/*
the files are called
b<i>_l<j>
where the algorithm to get mb and a is
my $bmin=4.55;
my $bstep=0.05;
my $l1min=-1.1;
my $l1step=0.1;
my $b=$bmin+$i*$bstep;
my $a=-1-3*(5.3-$b)*(5.3-$b)/($l1min+$j*$l1step);
I produced the files for the 1D reweighting with the 14 bins
double histbA[14] = {0.0, 1.2, 1.341, 1.414, 1.549, 1.612, 1.732,
1.843, 1.949, 2.5, 3.0, 3.5, 4.0, 4.5};
there is a fraction of low a runs that failed because the genarator is
unable to generate them!!!!
The points that failed are shown in the bottom plot of
http://www.slac.stanford.edu/~rfaccini/lost+found/weight_mba.ps
[x axis: Mb, y axis 'a')
(top plot is all generated points)
I would say that these points (which have empty weight file)
committed suicide and should not be considered
ciao
ric
On Fri, 27 Jun 2003, Riccardo Faccini wrote:
> ups
> >
> > I have then related a with mb and l1 with
> > a=1-3*(5.3-mb)*(5.3-mb)/l1
>
> I meant
> a=-1-3*(5.3-mb)*(5.3-mb)/l1
>
> > please let me know if you disagree with this.
> bye
> ric
>
|