SaturationFit Class Reference

#include <SaturationFit.h>

Public Member Functions
	SaturationFit (int N, float *charge_ptr, float *adc_ptr, float *u0_ptr, float *u1_ptr, float *u2_ptr, float *u3_ptr)
	~SaturationFit ()

Detailed Description

Definition at line 9 of file SaturationFit.h.

Constructor & Destructor Documentation

SaturationFit::SaturationFit ( int  N, float *  charge_ptr, float *  adc_ptr, float *  u0_ptr, float *  u1_ptr, float *  u2_ptr, float *  u3_ptr  )

inline

Definition at line 13 of file SaturationFit.h.

References gather_cfg::cout, SaturationFcn::set_data(), ctppsDiamondLocalTracks_cfi::sigma, SaturationFcn::x0start, SaturationFcn::x1start, SaturationFcn::x2start, and SaturationFcn::x3start.

                                                                                                                  {
   double u[3],sigma[3],chisq;
   VariableMetricMinimizer* pMinimizer=new VariableMetricMinimizer() ; 
   SaturationFcn* pFcn=new SaturationFcn();
   pFcn->set_data(N,charge_ptr,adc_ptr);
   std::vector<double> par(4,0);
   std::vector<double> err(4,0);
   //printf(" about to fill pars \n");
   par[0]=pFcn->x0start;
   par[1]=pFcn->x1start;
   par[2]=pFcn->x2start;
   par[3]=pFcn->x3start;
   //printf(" pars: %f %f %f %f \n",par[0],par[1],par[2],par[3]);
   err[0]=20.0;
   err[1]=0.0001;
   err[2]=1.0;
   err[3]=20.0;
   FunctionMinimum fmin = pMinimizer->Minimize(*pFcn, par, err, 1, 5000, 0.01);
   if( !fmin.IsValid()){
     printf(" minuit did not converge \n");
   } else {
     std::cout<<"fit succeeded... results: "<<fmin<<std::endl;
    chisq = fmin.Fval();
    u[0]     = fmin.UserParameters().Value( static_cast<unsigned int>(0) );
    *u0_ptr=u[0];
    sigma[0] = fmin.UserParameters().Error( static_cast<unsigned int>(0) );
    u[1]     = fmin.UserParameters().Value( static_cast<unsigned int>(1) );
    *u1_ptr=u[1];
    sigma[1] = fmin.UserParameters().Error( static_cast<unsigned int>(1) );
    u[2]     = fmin.UserParameters().Value( static_cast<unsigned int>(2) );
    *u2_ptr=u[2];
    sigma[2] = fmin.UserParameters().Error( static_cast<unsigned int>(2) );
    u[3]     = fmin.UserParameters().Value( static_cast<unsigned int>(3) );
    *u3_ptr=u[3];
    sigma[3] = fmin.UserParameters().Error( static_cast<unsigned int>(3) );
    //std::cout<<"fitresults: "<<u[0]<<" "<<u[1]<<" "<<u[2]<<" "<<u[3]<<std::endl;
    
   }
 }

SaturationFit::~SaturationFit ( )

inline

Definition at line 53 of file SaturationFit.h.

{}