Changeset 8746 for trunk/MagicSoft/Mars/mjobs
- Timestamp:
- 09/17/07 16:56:20 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/mjobs/MJSpectrum.cc
r8744 r8746 333 333 { 334 334 // Some debug output 335 *fLog << all << endl; 335 336 fLog->Separator("Compiling original MC distribution"); 336 337 … … 455 456 if (!weight.Set(*head)) 456 457 { 457 *fLog << err << "ERROR - Initializing MMcSpectrumWeight f ailed." << endl;458 *fLog << err << "ERROR - Initializing MMcSpectrumWeight from " << fname << " failed." << endl; 458 459 return kFALSE; 459 460 } … … 463 464 464 465 // No we found everything... go on reading contents 465 *fLog << inf << "Reading OriginalMC of " << fname << endl;466 *fLog << inf2 << "Reading OriginalMC of " << fname << endl; 466 467 467 468 // Fill histogram from tree … … 948 949 gPad->SetGridy(); 949 950 TH1D *spec=(TH1D*)spectrum.DrawCopy(); 950 FitSpectrum(*spec);951 //FitSpectrum(*spec); 951 952 952 953 c1.cd(3); … … 985 986 spec = (TH1D*)spectrum.DrawCopy(); 986 987 988 TLatex tex; 989 987 990 TF1 fc("Crab", "7.0e-6*(x/300)^(-2.31-0.26*log10(x/300))", 100, 6000); 988 991 fc.SetLineStyle(7); … … 991 994 fc.DrawCopy("same"); 992 995 996 tex.DrawText(90, fc.Eval(100), "Crab (\\Gamma=-2.31)"); 997 993 998 TF1 fa("PG1553", "1.8e-6*(x/200)^-4.21", 90, 600); 994 999 static_cast<const TAttLine&>(fc).Copy(fa); 995 1000 fa.DrawCopy("same"); 1001 1002 tex.DrawText(90, fa.Eval(90), "PG1553 (\\Gamma=-4.21)"); 996 1003 997 1004 // Display dN/dE*E^2 for conveinience … … 1329 1336 *fLog << inf; 1330 1337 fLog->Separator(GetDescriptor()); 1331 *fLog << "Compile Monte Carlo Sample (dataset " << set.GetName() << ")" << endl;1338 *fLog << "Compile Monte Carlo sample (dataset " << set.GetName() << ")" << endl; 1332 1339 *fLog << endl; 1333 1340 … … 1696 1703 cont.Add(&area1); 1697 1704 cont.Add(&hest); 1705 cont.Add(&set); 1698 1706 1699 1707 if (fDisplay)
Note:
See TracChangeset
for help on using the changeset viewer.