Changeset 2762 for trunk/MagicSoft/Mars/mhist

Timestamp:

01/09/04 22:47:53 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhist/MHCamera.cc (modified) (9 diffs)

trunk/MagicSoft/Mars/mhist/MHCamera.cc

@@ -103 +103 @@
 //  Default Constructor. To be used by the root system ONLY.
 //
-MHCamera::MHCamera() : TH1D(), fGeomCam(NULL), fColors(kItemsLegend)
+MHCamera::MHCamera() : TH1D(), fGeomCam(NULL), fColors(kItemsLegend), fYProj(NULL)
 {
     Init();
@@ -116 +116 @@
 MHCamera::MHCamera(const MGeomCam &geom, const char *name, const char *title)
 : fGeomCam(NULL), /*TH1D(name, title, geom.GetNumPixels(), -0.5, geom.GetNumPixels()-0.5),
-fUsed(geom.GetNumPixels()),*/ fColors(kItemsLegend)
+                    fUsed(geom.GetNumPixels()),*/ fColors(kItemsLegend), fYProj(NULL)
 {
     //fGeomCam = (MGeomCam*)geom.Clone();
@@ -160 +160 @@
 {
     if (fGeomCam)
-        delete fGeomCam;
+      delete fGeomCam;
     if (fNotify)
-        delete fNotify;
+      delete fNotify;
+    if (fYProj)
+      delete fYProj;
+
 }
 
@@ -399 +402 @@
     gPad->SetBorderMode(0);
 
+    TString opt(option);
+    opt.ToLower();
+
     AppendPad(option);
 
@@ -414 +420 @@
 // If this makes trouble please write a bug report.
 //
-TObject *MHCamera::DrawClone(Option_t *option) const
+TObject *MHCamera::DrawClone(Option_t *option) const 
 {
    // Draw a clone of this object in the current pad
@@ -423 +429 @@
 
    TObject *newobj = Clone();
+
    if (!newobj)
        return 0;
@@ -437 +444 @@
    */
 
-   const TString opt(option);
+   TString opt(option);
+   opt.ToLower();
+
    newobj->Draw(opt.IsNull() ? GetDrawOption() : option);
 
@@ -445 +454 @@
    return newobj;
 }
+
+void MHCamera::CreateProjection()
+{
+  
+  Int_t nbins = 50;
+
+  // Create the projection histogram
+  TString ytitle(GetYaxis()->GetTitle());
+  fYProj = new TH1D(ytitle.Data(),GetTitle(),nbins,GetMinimum()-0.1,GetMaximum()+0.1);
+  fYProj->SetXTitle(ytitle.Data());
+  fYProj->SetYTitle("Nr. of pixels");
+  fYProj->Sumw2();
+  fYProj->SetDirectory(NULL);
+
+  // Fill the projected histogram
+  Double_t cont;
+  for (Int_t binx =0;binx<=GetNbinsX()+1;binx++) {
+    cont  = GetBinContent(binx);
+    if (cont) 
+      fYProj->Fill(cont);
+  }
+}
+
 
 // ------------------------------------------------------------------------
@@ -592 +624 @@
     }
 
+    if (opt.Contains("proj"))
+    {
+
+      CreateProjection();
+      opt.ReplaceAll("proj", "");
+      fYProj->Paint(opt);
+      return;
+    }
+    
     gPad->Clear();