00001 #ifndef __INCLUDE_SFIT_BLAST_CC__
00002 #define __INCLUDE_SFIT_BLAST_CC__
00003 #include "sfit.h"
00004
00005 CCF2SFit_Blast::CCF2SFit_Blast(CSourceCalc *scset,C3DArray *cexp3Dset,C3DArray *cerror3Dset, C3DArray *ctheory3Dset,CCHArray *ctheoryset,CCHArray *sourceset,CKernel *kernelset){
00006 int i,j;
00007
00008
00009
00010 sourcecalc=scset;
00011 cexp3D=cexp3Dset;
00012 cerror3D=cerror3Dset;
00013 ctheoryCH=ctheoryset;
00014 ctheory3D=ctheory3Dset;
00015 sourceCH=sourceset;
00016 kernel=kernelset;
00017 Init();
00018
00019
00020
00021 AddPar("lambda",parameter::getD(sourcecalc->spars,"lambda",1),
00022 0.01,0.0,1.5);
00023 AddPar("R",parameter::getD(sourcecalc->spars,"R",13),
00024 0.05,5.0,20.0);
00025 AddPar("Tau",parameter::getD(sourcecalc->spars,"Tau",12),
00026 0.05,5.0,20.0);
00027 AddPar("DelTau",parameter::getD(sourcecalc->spars,"DelTau",5),
00028 0.05,0.0,20.0);
00029
00030
00031 for(i=0;i<nfreepars;i++){
00032 for(j=0;j<nfreepars;j++){
00033 StepMatrix[i][j]=0.0;
00034 ErrorMatrix[i][j]=0.0;
00035 if(i==j){
00036 StepMatrix[i][j]=par[i]->error;
00037 ErrorMatrix[i][j]=par[i]->error*par[i]->error;
00038 }
00039 }
00040 }
00041 }
00042 #endif
00043