Ignore:
Timestamp:
02/14/05 11:06:10 (20 years ago)
Author:
rwagner
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/mtemp/mmpi/MUnfoldSpectrum.cc

    r6384 r6452  
    2626//  MUnfoldSpectrum
    2727//
     28//  Unfolds a gamma spectrum using the algorithms given in the MUnfold class
     29//
    2830/////////////////////////////////////////////////////////////////////////////   
    2931#include "MUnfoldSpectrum.h"
     
    5456  //   TH3D* migration = migm->GetHist();
    5557 
    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;
    6164 
    6265  TAxis &taxis  = *fToBeUnfolded->GetYaxis();
     
    140143    Double_t blow = baxismig.GetBinLowEdge(1);
    141144    Double_t bup  = baxismig.GetBinLowEdge(nbmig+1);
    142      
    143    
     145         
    144146    // -----------------------------------------
    145147    // dummy ideal distribution
     
    260262      hb.SetBinError(a+1, sqrt(Vbcov(a, a)) );
    261263    }
    262 
    263264   
    264265    for (Int_t k=1; k<=nb; k++) {
     
    272273    }           
    273274  }
    274  
    275        
    276 
    277 
    278 
    279275}
    280276
Note: See TracChangeset for help on using the changeset viewer.