Changeset 861 for trunk/MagicSoft/Mars/mhist

Timestamp:

07/10/01 12:41:13 (23 years ago)

Author:

jlopez

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhist

Files:

• MHMcEnergies.cc (modified) (3 diffs)
• MHMcEnergies.h (modified) (2 diffs)
• MHMcEnergy.cc (modified) (11 diffs)
• MHMcEnergy.h (modified) (1 diff)

trunk/MagicSoft/Mars/mhist/MHMcEnergies.cc

@@ -23 +23 @@
 \* ======================================================================== */
 
+/////////////////////////////////////////////////////////////////////////////
+//
+//  MHMcEnergies
+//
+// This class holds and array of MHMcEnergy objects in order to be able
+// to compute the threshold for each of the different trigger conditions
+// in a root file.
+//
+////////////////////////////////////////////////////////////////////////////
 #include "MHMcEnergies.h" 
 
@@ -30 +39 @@
 ClassImp(MHMcEnergies);
 
+// --------------------------------------------------------------------------
+//
+// Default Constructor.
+//
 MHMcEnergies::MHMcEnergies(const UInt_t count, const char *name, const char *title)
     : fNum(count)
 { 
-    //
-    //   default constructor
-    //
     char aux[25]="MHMcEnergies";
     sprintf(aux+12, "[%i]", fNum);
@@ -48 +58 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Default Destructor.
+//
 MHMcEnergies::~MHMcEnergies()
 {
     delete fHists;
 }
-
+// --------------------------------------------------------------------------
+//
+// Add to the parameter list all the MHMcEnergy objects contained in the array.
+//
 void MHMcEnergies::AddEntriesToList(MParList *plist)
 {

trunk/MagicSoft/Mars/mhist/MHMcEnergies.h

@@ -20 +20 @@
 {
 private:
-    UInt_t        fNum;
-    TClonesArray *fHists; // histograms
+    UInt_t        fNum;   // Num of histograms
+    TClonesArray *fHists; // Array with the energy histograms
 
 public:
@@ -36 +36 @@
 
 #endif
+

trunk/MagicSoft/Mars/mhist/MHMcEnergy.cc

@@ -23 +23 @@
 \* ======================================================================== */
 
+/////////////////////////////////////////////////////////////////////////////
+//
+//  MHMcEnergy
+//
+// This class holds the information ( histogram and fit function )  
+// about the energy threshold for a particular trigger condition.
+//
+////////////////////////////////////////////////////////////////////////////
 #include "MHMcEnergy.h" 
 
@@ -34 +42 @@
 ClassImp(MHMcEnergy);
 
+// -------------------------------------------------------------------------
+//
+//  Default Constructor.
+//
 MHMcEnergy::MHMcEnergy(const UInt_t idx, const char *name, const char *title)
 { 
-    //
-    //   default constructor
-    //
     char aux[15]="MHMcEnergy";
 
@@ -58 +67 @@
     if (idx>0)
         sprintf(aux+10, ";%i", idx);
-    hLogEner = new TH1F(aux, "", 100, 0.5, 4.5);
+    hLogEner = new TH1F(aux, "", 40, 0.5, 4.5);
     hLogEner->SetXTitle("log(E) [GeV]");
     hLogEner->SetYTitle("dN/dE");
-    //hLogEner->SetBins(60);
 }
 
+//-------------------------------------------------------------------------
+//
+//  Defualt Destructor
+//
 MHMcEnergy::~MHMcEnergy()
 {
@@ -70 +82 @@
 }
 
+//--------------------------------------------------------------------------
+//
+//  Fill the histogram with the log10 of the energy for triggered events.
+//
 void MHMcEnergy::Fill(Float_t log10E, Float_t w)
 {
@@ -75 +91 @@
 }
 
+// -------------------------------------------------------------------------
+//
+// Fitting function
+//
 void MHMcEnergy::Fit(Axis_t xxmin, Axis_t xxmax)
 {
     //
     // 0: don't draw the function (it is drawn together with the histogram)
-    // +: add these function to the list of fits. Don't delete the last fit.
+    // Q: quiet mode
     //
-    // FIXME: R means: use the range specified in the function (xxmin, xxmax are ignored!)
-    //        Q means: quiet (why?)
-    //
-    //
-    hLogEner->Fit(fLogEner->GetName(), "Q0+", "", xxmin, xxmax);
+    hLogEner->Fit(fLogEner->GetName(), "Q0", "", xxmin, xxmax);
 }
 
+// ------------------------------------------------------------------------
+// 
+// Drawing function. It creates its own canvas.
+//
 void MHMcEnergy::Draw(Option_t *option)
 {
@@ -111 +131 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Set the number of bins in the histogran.
+//
+void MHMcEnergy::SetBins(Int_t nbins)
+{
+  hLogEner->SetBins(nbins,0.5,4.5);
+}
+// --------------------------------------------------------------------------
+//
+// Write the threshold and its error in the standard output
+//
 void MHMcEnergy::Print(Option_t*)
 {
@@ -116 +148 @@
 }
 
+// -------------------------------------------------------------------------
+//
+//  Return the threshold
+//
 Float_t MHMcEnergy::GetThreshold() const
 {
@@ -123 +159 @@
 }
 
+// -------------------------------------------------------------------------
+//
+// Return the error of the threshold.
+//
 Float_t MHMcEnergy::GetThresholdErr() const
 {
@@ -129 +169 @@
     const Float_t p1err = fLogEner->GetParError(1);
 
+    // The error has into accuont the error in the fit
     return pow(10, p1) * p1err * lg10;
 }
 
+// -------------------------------------------------------------------------
+//
+// Return the peak of the fitted gaussan function.
+//
 Float_t MHMcEnergy::GetGaussPeak() const
 {
@@ -137 +182 @@
 }
 
+// -------------------------------------------------------------------------
+//
+// Return the sigma of the fitted gaussan function.
+//
 Float_t MHMcEnergy::GetGaussSigma() const
 {
@@ -142 +191 @@
 }
 
+
+
+

trunk/MagicSoft/Mars/mhist/MHMcEnergy.h

@@ -32 +32 @@
     void Fill(Float_t log10E, Float_t w);
     void Fit(Axis_t xxmin, Axis_t xxmax);
+    void SetBins(Int_t nbins = 100);
 
     void Draw(Option_t* option = "");
     void Print(Option_t* option = NULL);
 
-    ClassDef(MHMcEnergy, 1)  // Histogram container for montecarlo energy
+    ClassDef(MHMcEnergy, 1)  // Histogram container for montecarlo energy threshold
 };