Changeset 6452 for trunk/MagicSoft/Mars/mtemp/mmpi

Timestamp:

02/14/05 11:06:10 (20 years ago)

Author:

rwagner

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mtemp/mmpi

Files:

• FluxCalcAndUnfold.C (added)
• MUnfoldSpectrum.cc (modified) (5 diffs)
• MUnfoldSpectrum.h (modified) (2 diffs)

trunk/MagicSoft/Mars/mtemp/mmpi/MUnfoldSpectrum.cc

@@ -26 +26 @@
 //  MUnfoldSpectrum
 //
+//  Unfolds a gamma spectrum using the algorithms given in the MUnfold class
+//
 /////////////////////////////////////////////////////////////////////////////   
 #include "MUnfoldSpectrum.h"
@@ -54 +56 @@
   //   TH3D* migration = migm->GetHist();
   
-  const Int_t energyBins = fToBeUnfolded->GetXaxis()->GetNbins();
-  const Int_t thetaBins =  fToBeUnfolded->GetYaxis()->GetNbins();
-  const TAxis* axisEnergy = fToBeUnfolded->GetXaxis();
-  const TAxis* axisTheta  = fToBeUnfolded->GetYaxis();
-  cout << "Distribution to be unfolded has " << energyBins << ", " << thetaBins << " bins" <<endl;
+  //const Int_t energyBins = fToBeUnfolded->GetXaxis()->GetNbins();
+  //const Int_t thetaBins =  fToBeUnfolded->GetYaxis()->GetNbins();
+  //const TAxis* axisEnergy = fToBeUnfolded->GetXaxis();
+  //const TAxis* axisTheta  = fToBeUnfolded->GetYaxis();
+  //cout << "Distribution to be unfolded has " << energyBins 
+  //     << ", " << thetaBins << " bins" <<endl;
   
   TAxis &taxis  = *fToBeUnfolded->GetYaxis();
@@ -140 +143 @@
     Double_t blow = baxismig.GetBinLowEdge(1);
     Double_t bup  = baxismig.GetBinLowEdge(nbmig+1);
-      
-    
+          
     // -----------------------------------------
     // dummy ideal distribution
@@ -260 +262 @@
       hb.SetBinError(a+1, sqrt(Vbcov(a, a)) );
     }
-
     
     for (Int_t k=1; k<=nb; k++) {
@@ -272 +273 @@
     }           
   }
-  
-        
-
-
-
 }
 

trunk/MagicSoft/Mars/mtemp/mmpi/MUnfoldSpectrum.h

@@ -13 +13 @@
 {
  public:
-    MUnfoldSpectrum();
-    ~MUnfoldSpectrum();
+  MUnfoldSpectrum();
+  ~MUnfoldSpectrum();
 
  private:
@@ -35 +35 @@
   void SetPriorInput(TH1D *hpr);
   void SetPriorPower(Double_t gamma);
+  void SetUnfoldingMethod(Short_t method) { fUnfoldingMethod = method; }
   
   ClassDef(MUnfoldSpectrum, 1)