Changeset 5490 for trunk/MagicSoft/Mars/mtemp/mmpi

Timestamp:

11/28/04 15:46:12 (20 years ago)

Author:

mazin

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mtemp/mmpi

Files:

• MPointingPosCalcWA.cc (added)
• MPointingPosCalcWA.h (added)
• MSkyPlot.cc (modified) (14 diffs)
• MSkyPlot.h (modified) (3 diffs)

trunk/MagicSoft/Mars/mtemp/mmpi/MSkyPlot.cc

@@ -123 +123 @@
 #include <TLatex.h>
 #include <TCanvas.h>
+#include <TPaveLabel.h>
 #include <TPaveText.h>
 #include <TStopwatch.h>
@@ -281 +282 @@
     }
 
+    kSaveAlphaPlots=kTRUE;
+    kSaveSkyPlots=kTRUE;
+    kSaveNexPlot=kTRUE;
     fAlphaHName = "alphaplot.root";
     fSkyHName   = "skyplot.root";
+    fNexHName   = "Nexcess.gif";
 }
 
@@ -352 +357 @@
     fMinYGrid    = ymin;
     fMaxYGrid    = ymax;
+    
+    *fLog << endl <<  inf << " SetSkyPlot: fMinXGrid, fMaxXGrid, fMinYGrid, fMaxYGrid, fBinStepGrid: " << endl;
+    *fLog << inf << "           " << fMinXGrid << ", " << fMaxXGrid << ", " << fMinYGrid << ", " 
+                 << fMaxYGrid<< ", " << fBinStepGrid << endl;
 }
 
@@ -473 +482 @@
     // for the current pointing position and add a offset in the
     // Fill function!
-    kSaveAlphaPlots=kTRUE;
-    kSaveSkyPlots=kTRUE;
 
    // prepare skyplot
@@ -480 +487 @@
     fNumStepsY     =  (int) ((fMaxYGrid - fMinYGrid) / fBinStepGrid + 1.5);
     fNumalphahist  =  (int) (fNumStepsX * fNumStepsY  + 0.5);
+
+    *fLog << inf << "SetSkyPlot: fNumStepsX, fNumStepsY, fNumalphahist: "
+                 << fNumStepsX << ", " << fNumStepsY << ", " << fNumalphahist << endl;
 
    // fHistSignif.SetName("SPSignif2ndOrder");
@@ -519 +529 @@
         fHistAlpha[i].SetDirectory(NULL);
     }
+//cout  << " fHistAlpha[10].GetBinContent(5) " << fHistAlpha[10].GetBinContent(5) << endl;
 
     delete temp;
@@ -626 +637 @@
 */
 
+// make the center of the plot different from the center of the camera
 /*
     x_0 = fPntPosCam->GetX()*fMm2Deg; 
     y_0 = fPntPosCam->GetY()*fMm2Deg; 
 */      
-    x_0 = 0.;
+    x_0 = 0.;  
     y_0 = 0.;
 
@@ -641 +653 @@
         for (Int_t gridx = 0; gridx < fNumStepsX; gridx++)
         {
+
            xsource = fMinXGrid + gridx * fBinStepGrid;
+
  /*     derotation    : rotate  into camera coordinates */
  /*     formel: (x_cam)    (cos(gam)  -sin(gam))   (xtilde)   (x_0)
@@ -649 +663 @@
            xsourcam = cosgam * xsource - singam * ysource + x_0;
            ysourcam = singam * xsource + cosgam * ysource + y_0;
+
+
  /*    end derotatiom    */
+//           xsourcam = xsource;
+//           ysourcam = ysource;
 
        /* calculate ALPHA und DIST according to the source position */
@@ -698 +716 @@
     const Int_t onbinsalpha = TMath::Nint(fAlphaONMax/fHistAlphaBinWidth);
 
-// fit with gaus 
-    fHistSignifGaus.Fit("gaus");
 
 // fit function for the background
@@ -713 +729 @@
     alpha_iterator = fHistAlpha.begin();
 
+    fHistNexcess.Reset();
+    fHistOn.Reset();
+    fHistSignif.Reset();   
+    fHistSignifGaus.Reset();
+    
     while( alpha_iterator != fHistAlpha.end() ) {
          // find the bin in the 2dim histogram
@@ -751 +772 @@
     }
 
+// fit with gaus 
+    fHistSignifGaus.Fit("gaus","N");
+
+
 //temp
 /*
@@ -770 +795 @@
 //  TString stri2 = "skyplots.root";
   if(kSaveSkyPlots==kTRUE) SaveSkyPlots(fSkyHName);
+
+// save nex plot:
+  if(kSaveNexPlot==kTRUE) SaveNexPlot(fNexHName);
 
   fHistAlpha.clear();
@@ -800 +828 @@
 }
 
+void MSkyPlot::SaveNexPlot(TString nexplotname)
+{
+    gStyle->SetCanvasBorderMode(0);
+    gStyle->SetCanvasBorderSize(0);
+    gStyle->SetCanvasColor(10);
+//    gStyle->SetPadBorderMode(0);
+//    gStyle->SetPadBorderSize(0);
+    gStyle->SetPadColor(10);
+    gStyle->SetOptFit(0);
+    gStyle->SetOptStat(0);
+    gStyle->SetStatColor(10);
+    gStyle->SetPalette(1);
+    gStyle->SetPadRightMargin(0.16);
+    gStyle->SetPadLeftMargin(0.13);
+
+    TH2D tmp = fHistNexcess;
+    tmp.SetMaximum(470); 
+    tmp.SetMinimum(-90); 
+    TCanvas can("SkyPlot","SkyPlot", 0, 0, 800, 400);
+    can.Divide(2,1);
+    can.cd(1);
+    tmp.GetYaxis()->SetTitleOffset(1.3);
+    tmp.Draw("colz"); 
+    TPaveLabel myname(fMinXGrid-0.5,fMinYGrid-0.4,fMinXGrid+0.3,fMinYGrid-0.2,"by D. Mazin");
+    myname.Draw();
+
+    can.cd(2);
+    fHistNexcess.SetMaximum(470);    
+    fHistNexcess.SetMinimum(-40);    
+    fHistNexcess.GetXaxis()->SetTitleOffset(1.3);
+    fHistNexcess.GetYaxis()->SetTitleOffset(1.6);
+    gPad->SetTheta(20);
+    fHistNexcess.Draw("lego2");
+    can.Print(nexplotname);  
+//    can.Print("test.root");
+}
 
 void MSkyPlot::SaveSkyPlots(TString skyplotfilename)

trunk/MagicSoft/Mars/mtemp/mmpi/MSkyPlot.h

@@ -80 +80 @@
     Bool_t kSaveAlphaPlots;
     Bool_t kSaveSkyPlots;
+    Bool_t kSaveNexPlot;
 
     Float_t fMinXGrid;                  //  [deg] , left edge of the skyplot
@@ -105 +106 @@
     TString fAlphaHName;          // name for histogram with alpha plots
     TString fSkyHName;            // name for histogram with sky plots
+    TString fNexHName;            // name for canvas with Nex plot
 
     Int_t DistancetoPrimitive(Int_t px, Int_t py);
@@ -128 +130 @@
     void SetSizeMin(Float_t size) { fSizeMin = size; } // Absolute minimum Size
     void SetSizeMax(Float_t size) { fSizeMax = size; } // Absolute maximum Size
-//    void SetMinLD(Float_t ratio)  { fMinLD = ratio; }  // Minimum ratio between length/dist
-//    void SetMaxLD(Float_t ratio)  { fMaxLD = ratio; }  // Maximum ratio between length/dist
     void ReadCuts(const TString parSCinit);
     void SetSkyPlot(Float_t xmin, Float_t xmax, Float_t ymin, Float_t ymax, Float_t step);
     Double_t CalcLimit(Double_t *a, Double_t ls, Double_t ls2, Double_t dd2);
-    //void SaveSkyPlots(const TString skyplotfilename="skyplots.root");
+
+    void SetOutputSkyName(TString outname2)     { fSkyHName = outname2; }
+    void SaveSkyPlots(TString stri);
+    void SetNotSaveSkyPlots()                   { kSaveSkyPlots = kFALSE; }
+
 
     void SetOutputAlphaName(TString outname1)   { fAlphaHName = outname1; }
-    void SaveSkyPlots(TString stri);
+    void SaveAlphaPlots(const TString stri2);
+    void SetNotSaveAlphaPlots()                 { kSaveAlphaPlots = kFALSE; }
 
+    void SetOutputNexName(TString outname3)     { fNexHName = outname3; }
+    void SaveNexPlot(const TString stri3);
+    void SetNotSaveNexPlot()                    { kSaveNexPlot = kFALSE; }
 
-    //void SaveAlphaPlots(const TString alphaplotfilename="alphaplots.root");
-    void SetOutputSkyName(TString outname2)     { fSkyHName = outname2; }
-    void SaveAlphaPlots(const TString stri2);
 
     void SetAlphaCut(Float_t alpha);