00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016 #include "TGo4FitDataHistogram.h"
00017
00018 #include <iostream.h>
00019
00020 #include "TArrayD.h"
00021 #include "TArrayF.h"
00022
00023 TGo4FitDataHistogram::TGo4FitDataHistogram() :
00024 TGo4FitData(), fxHistogram(this, TH1::Class()) {
00025 }
00026
00027 TGo4FitDataHistogram::TGo4FitDataHistogram(const char* iName, TH1 *iHistogram, Bool_t iHistogramOwned, Bool_t AddAmpl) :
00028 TGo4FitData(iName,"TH1 and derived classes", dtHistogram, AddAmpl),
00029 fxHistogram("Histogram","Histogram object of TH1 class ",this, TH1::Class(), kTRUE, iHistogram, iHistogramOwned) {
00030 }
00031
00032 TGo4FitDataHistogram::~TGo4FitDataHistogram() {
00033 }
00034
00035 TGo4FitDataIter* TGo4FitDataHistogram::MakeIter() {
00036 return new TGo4FitDataHistogramIter(this);
00037 }
00038
00039 void TGo4FitDataHistogram::FillSlotList(TSeqCollection* list) {
00040 TGo4FitData::FillSlotList(list);
00041 list->Add(&fxHistogram);
00042 }
00043
00044 void TGo4FitDataHistogram::Print(Option_t* option) const {
00045 TGo4FitData::Print(option);
00046 cout << " ";
00047 fxHistogram.Print(option);
00048 }
00049
00050 ClassImp(TGo4FitDataHistogram)
00051
00052
00053
00054
00055
00056 TGo4FitDataHistogramIter::TGo4FitDataHistogramIter() : TGo4FitDataIter(),
00057 fxData(0), fxLimits(), fxOwnScales(), fxOwnWidths() {
00058 }
00059
00060 TGo4FitDataHistogramIter::TGo4FitDataHistogramIter(TGo4FitDataHistogram* data) :
00061 TGo4FitDataIter(), fxData(data), fxLimits(), fxOwnScales(), fxOwnWidths() {
00062 }
00063
00064 TGo4FitDataHistogramIter::~TGo4FitDataHistogramIter() {
00065 }
00066
00067 Bool_t TGo4FitDataHistogramIter::StartReset() {
00068 if (fxData==0) return kFALSE;
00069 TH1* histo = fxData->GetHistogram();
00070 if (histo==0) return kFALSE;
00071
00072 Int_t NumDimen = histo->GetDimension();
00073 fxLimits.Set(NumDimen);
00074 fxLimits[0] = histo->GetNbinsX();
00075 if (NumDimen>1) fxLimits[1] = histo->GetNbinsY();
00076 if (NumDimen>2) fxLimits[2] = histo->GetNbinsZ();
00077
00078 fxOwnScales.Set(NumDimen); fxOwnScales.Reset(0.);
00079 fxOwnWidths.Set(NumDimen); fxOwnWidths.Reset(0.);
00080
00081 return ReserveArrays(NumDimen, NumDimen, kTRUE);
00082 }
00083
00084 Bool_t TGo4FitDataHistogramIter::ReadCurrentPoint() {
00085 if (fxData==0) return kFALSE;
00086 TH1* histo = fxData->GetHistogram();
00087 if (histo==0) return kFALSE;
00088
00089 switch (histo->GetDimension()) {
00090 case 1: fdValue = histo->GetBinContent(fxIndexes[0]+1); break;
00091 case 2: fdValue = histo->GetBinContent(fxIndexes[0]+1, fxIndexes[1]+1); break;
00092 case 3: fdValue = histo->GetBinContent(fxIndexes[0]+1, fxIndexes[1]+1, fxIndexes[2]+1); break;
00093 default: fdValue = 0;
00094 }
00095
00096 if (!GetDeviation()) fdStandardDeviation = fdValue;
00097
00098 for (Int_t n=0;n<fxOwnScales.GetSize();n++)
00099 switch (n) {
00100 case 0: fxOwnScales[0] = histo->GetXaxis()->GetBinCenter(fxIndexes[0]+1); break;
00101 case 1: fxOwnScales[1] = histo->GetYaxis()->GetBinCenter(fxIndexes[1]+1); break;
00102 case 2: fxOwnScales[2] = histo->GetZaxis()->GetBinCenter(fxIndexes[2]+1); break;
00103 }
00104
00105 for (Int_t n=0;n<fxOwnWidths.GetSize();n++)
00106 switch (n) {
00107 case 0: fxOwnWidths[0] = histo->GetXaxis()->GetBinWidth(fxIndexes[0]+1); break;
00108 case 1: fxOwnWidths[1] = histo->GetYaxis()->GetBinWidth(fxIndexes[1]+1); break;
00109 case 2: fxOwnWidths[2] = histo->GetZaxis()->GetBinWidth(fxIndexes[2]+1); break;
00110 }
00111
00112 return ProduceScales(fxIndexes.GetArray(), fxOwnScales.GetArray(), fxOwnWidths.GetArray());
00113 }
00114
00115 Bool_t TGo4FitDataHistogramIter::ShiftToNextPoint() {
00116 return NextIndex(fxIndexes,fxLimits);
00117 }
00118
00119 ClassImp(TGo4FitDataHistogramIter)
00120
00121
00122
00123