Changeset 6452
- Timestamp:
- 02/14/05 11:06:10 (20 years ago)
- Location:
- trunk/MagicSoft/Mars/mtemp
- Files:
-
- 1 added
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/mtemp/Changelog
r6447 r6452 18 18 19 19 -*-*- END OF LINE -*-*- 20 20 21 21 22 2005/02/14: Hendrik Bartko … … 25 26 26 27 28 29 2005/02/14: Robert Wagner 30 * mtemp/mmpi/FluxCalcAndUnfold.C 31 - Example macro for calculating a flux and unfolding it using 32 MUnfoldSpectrum and MUnfold 33 34 35 36 27 37 2005/02/11: Robert Wagner 28 38 * mtemp/mmpi/MUnfoldSpectrum.[cc,h] 29 39 - A wrapper class which applies the general unfolding algorithms 30 40 of MUnfold to a gamma spectrum 41 31 42 32 43 -
trunk/MagicSoft/Mars/mtemp/mmpi/MUnfoldSpectrum.cc
r6384 r6452 26 26 // MUnfoldSpectrum 27 27 // 28 // Unfolds a gamma spectrum using the algorithms given in the MUnfold class 29 // 28 30 ///////////////////////////////////////////////////////////////////////////// 29 31 #include "MUnfoldSpectrum.h" … … 54 56 // TH3D* migration = migm->GetHist(); 55 57 56 const Int_t energyBins = fToBeUnfolded->GetXaxis()->GetNbins(); 57 const Int_t thetaBins = fToBeUnfolded->GetYaxis()->GetNbins(); 58 const TAxis* axisEnergy = fToBeUnfolded->GetXaxis(); 59 const TAxis* axisTheta = fToBeUnfolded->GetYaxis(); 60 cout << "Distribution to be unfolded has " << energyBins << ", " << thetaBins << " bins" <<endl; 58 //const Int_t energyBins = fToBeUnfolded->GetXaxis()->GetNbins(); 59 //const Int_t thetaBins = fToBeUnfolded->GetYaxis()->GetNbins(); 60 //const TAxis* axisEnergy = fToBeUnfolded->GetXaxis(); 61 //const TAxis* axisTheta = fToBeUnfolded->GetYaxis(); 62 //cout << "Distribution to be unfolded has " << energyBins 63 // << ", " << thetaBins << " bins" <<endl; 61 64 62 65 TAxis &taxis = *fToBeUnfolded->GetYaxis(); … … 140 143 Double_t blow = baxismig.GetBinLowEdge(1); 141 144 Double_t bup = baxismig.GetBinLowEdge(nbmig+1); 142 143 145 144 146 // ----------------------------------------- 145 147 // dummy ideal distribution … … 260 262 hb.SetBinError(a+1, sqrt(Vbcov(a, a)) ); 261 263 } 262 263 264 264 265 for (Int_t k=1; k<=nb; k++) { … … 272 273 } 273 274 } 274 275 276 277 278 279 275 } 280 276 -
trunk/MagicSoft/Mars/mtemp/mmpi/MUnfoldSpectrum.h
r6384 r6452 13 13 { 14 14 public: 15 16 15 MUnfoldSpectrum(); 16 ~MUnfoldSpectrum(); 17 17 18 18 private: … … 35 35 void SetPriorInput(TH1D *hpr); 36 36 void SetPriorPower(Double_t gamma); 37 void SetUnfoldingMethod(Short_t method) { fUnfoldingMethod = method; } 37 38 38 39 ClassDef(MUnfoldSpectrum, 1)
Note:
See TracChangeset
for help on using the changeset viewer.