[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: fits with full kinetics function



You mean lambda_of, right?
I had a factor of 2 mistake in the x-axis scaling.  Also, the
parameter for lambda_op below should be this (in 1/ns)
   func->FixParameter(TLifetimeFit::kParKin_lambda_op, 0.08e-3);
not 0.08e-6 as I had before.  The corrected files are attached.
It looks like about 22 Hz for l_of 0->2.3, and I used a density
of 0.0116.  (I'm not sure why the fit doesn't give exactly 664 Hz
for l_of = 0...maybe because I didn't set lambda_pf to 0.)

Sorry for the mistake earlier.  Now the code takes the parameter
value from the x-axis of the results histogram, not scaled by
hand as before.

Steven


On Wed, Oct 04, 2006 at 09:15:41PM -0500, Peter Kammel wrote:
> Dear Steve,
> 
> You get a 66 Hz change for l_op 0->2.3, I get only 40ppm*0.45=19 Hz.
> That is a major discrepancy. Did you account for the fact that
> the effective l_op has to be multiplied by the 0.01 density
> relative to LH2 (l_op=2.3 for LH2).
> 
> Best regards
> 
> Peter
> 
> 
> On Wed, 4 Oct 2006, Steven Clayton wrote:
> 
> > ROOT files are attached.  To save as a GIF, EPS or whatever,
> > start a ROOT session and do this:
> > root [0] TFile f("LifetimeStudyFullKin_vs_lambda_op.root")
> > root [1] LifetimeStudyFullKin->Draw()
> > (save via the file menu of the canvas)
> > and similarly for the other file, LifetimeStudyFullKin_vs_lambda_of.root.
> > 
> > The fixed parameters are:
> >   // All rates specified should be in 1/ns.
> >   // Result histograms are scaled to give 1/s.
> >   //   lambda0 = lambda_fit - 664e-9 1/ns
> >   //   Lambda_s = 664e-9 1/ns
> >   //   where lambda_fit is from the 3-parameter fit.
> >   func->FixParameter(TLifetimeFit::kParKin_lambda0, ltfit->GetR());
> >   func->FixParameter(TLifetimeFit::kParKin_Lambda_t, 12.0e-9);
> >   func->FixParameter(TLifetimeFit::kParKin_gamma_o, 1.009/2.);
> >   func->FixParameter(TLifetimeFit::kParKin_gamma_p, 1.143/2.);
> >   func->FixParameter(TLifetimeFit::kParKin_phi, 0.0116);
> >   func->FixParameter(TLifetimeFit::kParKin_lambda_pf, 0.0074e-3);
> >   if(fVaryLambdaOP) {
> >     func->FixParameter(TLifetimeFit::kParKin_lambda_of, 1.80e-3);
> >     func->FixParameter(TLifetimeFit::kParKin_lambda_op, val);
> >   } else {
> >     func->FixParameter(TLifetimeFit::kParKin_lambda_of, val);
> >     func->FixParameter(TLifetimeFit::kParKin_lambda_op, 0.08e-6);
> >   }
> >   func->FixParameter(TLifetimeFit::kParKin_cZ, 0);
> >   func->FixParameter(TLifetimeFit::kParKin_lambda_Ztr, 0);
> >   func->FixParameter(TLifetimeFit::kParKin_Lambda_Z, 0);
> > 
> > 

Attachment: LifetimeStudyFullKin_vs_lambda_of.root
Description: Binary data

Attachment: LifetimeStudyFullKin_vs_lambda_op.root
Description: Binary data