Changeset 7691 for trunk/MagicSoft/Mars/mhflux
- Timestamp:
- 05/05/06 07:30:44 (19 years ago)
- Location:
- trunk/MagicSoft/Mars/mhflux
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/mhflux/MHEnergyEst.cc
r7409 r7691 71 71 fTitle = title ? title : "3-D histogram E-true E-est Theta"; 72 72 73 fNameEnergy = "MEnergyEst"; 74 fNameResult = "MinimizationValue"; 75 73 76 fHEnergy.SetDirectory(NULL); 74 77 fHEnergy.SetName("EnergyEst"); … … 184 187 } 185 188 189 // -------------------------------------------------------------------------- 190 // 191 // Divide chisq and bias by number of executions 192 // Print result 193 // 186 194 Bool_t MHEnergyEst::Finalize() 187 195 { … … 191 199 fResult->SetVal(fChisq); 192 200 201 *fLog << all << endl; 193 202 Print(); 194 203 … … 196 205 } 197 206 207 // -------------------------------------------------------------------------- 208 // 209 // Print result 210 // 198 211 void MHEnergyEst::Print(Option_t *o) const 199 212 { 200 // const Double_t resp = TMath::Power(10, res)-1;201 // const Double_t resm = 1-TMath::Power(10, -res);202 203 213 const Double_t res = TMath::Sqrt(fChisq-fBias*fBias); 204 214 if (TMath::IsNaN(res)) … … 209 219 210 220 TH1D *h = (TH1D*)fHResolution.ProjectionZ("Resolution"); 211 h->Fit("gaus", "Q0" );221 h->Fit("gaus", "Q0", "", -1.0, 0.25); 212 222 213 223 TF1 *f = h->GetFunction("gaus"); … … 225 235 } 226 236 237 // -------------------------------------------------------------------------- 238 // 239 // Return Correction Coefficients (weights) 240 // hist = E_mc/E_est 241 // 227 242 void MHEnergyEst::GetWeights(TH1D &hist) const 228 243 { … … 301 316 } 302 317 303 hx->Fit("gaus", "Q" );318 hx->Fit("gaus", "Q", "", -1.0, 0.25); 304 319 gPad=NULL; 305 320 gStyle->SetOptFit(101); … … 545 560 fMatrix = NULL; 546 561 } 547 -
trunk/MagicSoft/Mars/mhflux/MHEnergyEst.h
r6954 r7691 22 22 { 23 23 private: 24 TString fNameEnergy; 25 TString fNameResult; 26 24 27 MMcEvt *fMcEvt; //! 25 28 MParameterD *fEnergy; //! 26 29 MParameterD *fResult; //! 27 30 28 Int_t fMap[100]; // FIXME!31 Int_t fMap[100]; // FIXME! 29 32 MHMatrix *fMatrix; //! 30 33 … … 40 43 41 44 Double_t GetVal(Int_t i) const; 45 void CalcChisq(Double_t &chisq, Double_t &prob) const; 42 46 43 47 Bool_t SetupFill(const MParList *pList); … … 58 62 void Print(Option_t *o="") const; 59 63 60 ClassDef(MHEnergyEst, 1) //64 ClassDef(MHEnergyEst, 2) // 61 65 62 66 };
Note:
See TracChangeset
for help on using the changeset viewer.