Changeset 3551 for trunk/MagicSoft/Mars/mhist

Timestamp:

03/19/04 16:56:19 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhist

Files:

• MHCamera.cc (modified) (2 diffs)
• MHCamera.h (modified) (1 diff)

trunk/MagicSoft/Mars/mhist/MHCamera.cc

@@ -594 +594 @@
 TProfile *MHCamera::RadialProfileS(const TArrayI &sector, const TArrayI &aidx, const char *name, const Int_t nbins) const
 {
-
-    // Create the projection histogram
-    TString pname(name);
-    if (name=="_rad")
-    {
-        pname.Prepend(GetName());
-        if (sector.GetSize()>0)
+  
+  // Create the projection histogram
+  TString pname(name);
+  if (name=="_rad")
+    {
+      pname.Prepend(GetName());
+      if (sector.GetSize()>0)
         {
-            pname += ";";
-            for (int i=0; i<sector.GetSize(); i++)
-                pname += sector[i];
+          pname += ";";
+          for (int i=0; i<sector.GetSize(); i++)
+            pname += sector[i];
         }
-        if (aidx.GetSize()>0)
+      if (aidx.GetSize()>0)
         {
-            pname += ";";
-            for (int i=0; i<aidx.GetSize(); i++)
-                pname += aidx[i];
+          pname += ";";
+          for (int i=0; i<aidx.GetSize(); i++)
+            pname += aidx[i];
         }
     }
-
-    TProfile *h1=0;
-
-    //check if histogram with identical name exist
-    TObject *h1obj = gROOT->FindObject(pname);
-    if (h1obj && h1obj->InheritsFrom("TProfile")) {
-        h1 = (TProfile*)h1obj;
-        h1->Reset();
-    }
-
-    if (!h1)
-    {
-
+  
+  TProfile *h1=0;
+  
+  //check if histogram with identical name exist
+  TObject *h1obj = gROOT->FindObject(pname);
+  if (h1obj && h1obj->InheritsFrom("TProfile")) {
+    h1 = (TProfile*)h1obj;
+    h1->Reset();
+  }
+  
+  if (!h1)
+    {
+      
       Double_t min = 0.;
       Double_t max = fGeomCam->GetMaxRadius();
@@ -638 +638 @@
       h1->SetYTitle(GetYaxis()->GetTitle());
     }
-    
-    // Fill the projected histogram
-    for (Int_t idx=0; idx<fNcells-2; idx++)
-      if (IsUsed(idx) && MatchSector(idx, sector, aidx))
-        h1->Fill(TMath::Hypot((*fGeomCam)[idx].GetX(),(*fGeomCam)[idx].GetY()),
-                 GetBinContent(idx+1));
-    return h1;
+  
+  // Fill the projected histogram
+  for (Int_t idx=0; idx<fNcells-2; idx++)
+    if (IsUsed(idx) && MatchSector(idx, sector, aidx))
+      h1->Fill(TMath::Hypot((*fGeomCam)[idx].GetX(),(*fGeomCam)[idx].GetY()),
+               GetBinContent(idx+1));
+  return h1;
 }
 

trunk/MagicSoft/Mars/mhist/MHCamera.h

@@ -236 +236 @@
 
     TProfile *RadialProfile(const char *name="_rad") const { return  RadialProfileS(TArrayI(), TArrayI(), name);}
-    TProfile *RadialProfileS(Int_t sector, Int_t aidx, const char *name="_rad", const Int_t nbins=20) const
+    TProfile *RadialProfileS(Int_t sector, Int_t aidx, const char *name="_rad", const Int_t nbins=25) const
     {
         return RadialProfileS(TArrayI(1, &sector), TArrayI(1, &aidx), name, nbins);
     }
-    TProfile *RadialProfileS(const TArrayI &sector, const TArrayI &aidx, const char *name="_rad", const Int_t nbins=20) const;
+    TProfile *RadialProfileS(const TArrayI &sector, const TArrayI &aidx, const char *name="_rad", const Int_t nbins=25) const;
     
     const MGeomCam &GetGeomCam() const { return *fGeomCam; }