Changeset 5955

Timestamp:

01/24/05 10:28:16 (20 years ago)

Author:

domingo

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mtemp/mifae

Files:

• Changelog (modified) (1 diff)
• library/MDispCalc.cc (modified) (2 diffs)
• library/MFDisp.cc (modified) (3 diffs)
• library/MFDisp.h (modified) (2 diffs)
• library/MHDisp.cc (modified) (3 diffs)

trunk/MagicSoft/Mars/mtemp/mifae/Changelog

@@ -18 +18 @@
 
                                                  -*-*- END OF LINE -*-*-
+
+ 2005/01/24 Eva Domingo
+    * library/MFDisp.[cc,h] 
+     - Number of islands read from MImagePar container
+       to be used as cut parameter if wanted
+    * library/MHDisp.cc 
+     - corrected to be able to run Disp method also over
+       data (MMcEvt container info just used if available)
+    * library/MDispCalc.cc
+     - introduced DISP = Xi*(1/(1+width/length)) as default
+       Disp parameterization for the moment
 
  2005/01/24 Eva Domingo

trunk/MagicSoft/Mars/mtemp/mifae/library/MDispCalc.cc

@@ -300 +300 @@
     //      + (p[2] + p[3]*(logsize-fLogSize0))*width/length;
 
-    Double_t disp = p[0] + p[1]*(logsize-fLogSize0) + p[4]*(length-fLength0) 
-      + (p[2] + p[3]*(logsize-fLogSize0))*width/length;
+//    Double_t disp = p[0] + p[1]*(logsize-fLogSize0) + p[4]*(length-fLength0) 
+//      + (p[2] + p[3]*(logsize-fLogSize0))*width/length;
 
     //    Double_t disp = p[0] + p[1]*(logsize-fLogSize0) + p[4]*(length-fLength0) 
@@ -327 +327 @@
     //                    *(1-width/length);
 
+    Double_t disp = ( p[0] + p[1]*(logsize) 
+                    + p[2]*pow(logsize,2)
+                    + p[3]*pow(logsize,3) 
+                    + p[4]*pow(logsize,4) )
+                        *( 1/(1+width/length) );
     return disp;
 }

trunk/MagicSoft/Mars/mtemp/mifae/library/MFDisp.cc

@@ -40 +40 @@
 #include "MHillas.h"
 #include "MHillasSrc.h"
+#include "MImagePar.h"
 #include "MNewImagePar.h"
 
@@ -130 +131 @@
     }
 
+    fImgPar = (MImagePar*)pList->FindObject("MImagePar");
+    if (!fImgPar)
+    {
+        *fLog << err << "MImagePar not found... aborting." << endl;
+        return kFALSE;
+    }
+
     fNewImgPar = (MNewImagePar*)pList->FindObject("MNewImagePar");
     if (!fNewImgPar)
@@ -200 +208 @@
     const Int_t numusedpixels = fNewImgPar->GetNumUsedPixels();
     const Int_t numcorepixels = fNewImgPar->GetNumCorePixels();
-    //const Int_t numislands    = fNewImgPar->GetNumIslands();
+    const Int_t numislands    = fImgPar->GetNumIslands();
 
     fResult = kFALSE;
 
-    //if  (numislands<fIslandsMin  || numislands>fIslandsMax )
-    //    return Set(1);
+    if  (numislands<fIslandsMin  || numislands>fIslandsMax )
+        return Set(1);
 
     if  (numusedpixels<fUsedPixelsMin  || numusedpixels>fUsedPixelsMax )

trunk/MagicSoft/Mars/mtemp/mifae/library/MFDisp.h

@@ -8 +8 @@
 class MHillas;
 class MHillasSrc;
+class MImagePar;
 class MNewImagePar;
 
@@ -16 +17 @@
     MHillas      *fHil;
     MHillasSrc   *fHilSrc;
+    MImagePar    *fImgPar;
     MNewImagePar *fNewImgPar;
 

trunk/MagicSoft/Mars/mtemp/mifae/library/MHDisp.cc

@@ -299 +299 @@
     if (!fMcEvt)
     {
-        *fLog << err << "MMcEvt not found... aborting." << endl;
-        return kFALSE;
+        *fLog << err << "MMcEvt not found... This is not a MC file,"
+              << " you are not trying to optimize Disp, just calculating it."
+              << endl;
+        //        return kFALSE;
     }
 
@@ -424 +426 @@
       xmax     = fMatrix ? GetVal(kLongitmax)  : fMcEvt->GetLongitmax();      
     }
-    else
-      *fLog << "No MMcEvt container available (no Energy,ImpactPar or Xmax)" << endl;
 
     Double_t theta      = fMatrix ? GetVal(kZd)  : fPointing->GetZd();
@@ -504 +504 @@
     
     // Fill Energy, Size and ZA distributions
-    fHistEnergy->Fill(log10(energy));
+    if (fMatrix || (!fMatrix && fMcEvt))
+      fHistEnergy->Fill(log10(energy));
     fHistSize->Fill(log10(size));
     fHistcosZA->Fill(cos(theta/kRad2Deg));
     
     // to check the size and energy dependence of the optimization
-    fHistMinParEnergy->Fill(log10(energy),d2);
+    if (fMatrix || (!fMatrix && fMcEvt))
+    {
+      fHistMinParEnergy->Fill(log10(energy),d2);
+      fHistDuEnergy->Fill(log10(energy),Du);
+      fHistDvEnergy->Fill(log10(energy),Dv);
+    }
     fHistMinParSize->Fill(log10(size),d2);
-    fHistDuEnergy->Fill(log10(energy),Du);
     fHistDuSize->Fill(log10(size),Du);
-    fHistDvEnergy->Fill(log10(energy),Dv);
     fHistDvSize->Fill(log10(size),Dv);