Changeset 1668 for trunk/MagicSoft/Mars/mhist/MHAlphaEnergyTheta.cc

Timestamp:

11/25/02 09:12:47 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhist/MHAlphaEnergyTheta.cc (modified) (10 diffs)

trunk/MagicSoft/Mars/mhist/MHAlphaEnergyTheta.cc

@@ -64 +64 @@
     fHist.SetDirectory(NULL);
 
-    fHist.SetTitle("3D-plot of alpha, E-est, Theta");
+    fHist.SetTitle("3D-plot of alpha, E_{est}, Theta");
     fHist.SetXTitle("\\alpha [\\circ]");
-    fHist.SetYTitle("E-est [GeV]            ");
+    fHist.SetYTitle("E_{est} [GeV]");
     fHist.SetZTitle("\\Theta [\\circ]");
 }
@@ -79 +79 @@
    if (!fEnergy)
    {
-       *fLog << err << dbginf << "MHAlphaEnergyTheta : MEnergyEst not found... aborting." << endl;
+       *fLog << err << dbginf << "MEnergyEst not found... aborting." << endl;
        return kFALSE;
    }
@@ -86 +86 @@
    if (!fMcEvt)
    {
-       *fLog << err << dbginf << "MHAlphaEnergyTheta : MMcEvt not found... aborting." << endl;
+       *fLog << err << dbginf << "MMcEvt not found... aborting." << endl;
        return kFALSE;
    }
@@ -95 +95 @@
    if (!binsenergy || !binsalphaflux || !binstheta)
    {
-       *fLog << err << dbginf << "MHAlphaEnergyTheta : At least one MBinning not found... aborting." << endl;
+       *fLog << err << dbginf << "At least one MBinning not found... aborting." << endl;
        return kFALSE;      
    }
@@ -103 +103 @@
    fHist.Sumw2(); 
 
-   fHist.Sumw2(); 
-
    return kTRUE;
 }
@@ -116 +114 @@
     MHillasSrc &hil = *(MHillasSrc*)par;
 
-    fHist.Fill(fabs(hil.GetAlpha()), fEnergy->GetEnergy(), fMcEvt->GetTheta()*kRad2Deg);
+    fHist.Fill(fabs(hil.GetAlpha()), fEnergy->GetEnergy(),
+               fMcEvt->GetTelescopeTheta()*kRad2Deg);
 
     return kTRUE;
@@ -148 +147 @@
 
     h->SetTitle("Distribution of E-est [GeV]");
-    h->SetXTitle("E-est [GeV]            ");
+    h->SetXTitle("E_{est} [GeV]");
     h->SetYTitle("Counts");
 
@@ -186 +185 @@
 
     c.cd(1);
-    h = ((TH3*)(&fHist))->Project3D("expro");
+    h = ((TH3*)(&fHist))->Project3D(fName+"_expro");
 
     h->SetTitle("Distribution of \\alpha [\\circ]");
@@ -196 +195 @@
 
     c.cd(2);
-    h = ((TH3*)(&fHist))->Project3D("eypro");
+    h = ((TH3*)(&fHist))->Project3D(fName+"_eypro");
 
     h->SetTitle("Distribution of E-est [GeV]");
-    h->SetXTitle("E-est [GeV]            ");
+    h->SetXTitle("E_{est} [GeV]");
     h->SetYTitle("Counts");
 
@@ -207 +206 @@
 
     c.cd(3);
-    h = ((TH3*)(&fHist))->Project3D("ezpro");
+    h = ((TH3*)(&fHist))->Project3D(fName+"_ezpro");
 
     h->SetTitle("Distribution of \\Theta [\\circ]");