Changeset 8957 for trunk


Ignore:
Timestamp:
06/14/08 16:55:58 (16 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars
Files:
24 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/Changelog

    r8956 r8957  
    8989       the default was not correctly treated because char* was
    9090       compared with char*
     91
     92   * mbase/MLog.cc, mbase/MMath.cc, mbase/MParContainer.cc,
     93     mcalib/MCalibrateData.cc, mcalib/MCalibrationPulseTimeCam.cc,
     94     mhbase/MBinning.h, mhbase/MH.cc, mhbase/MHMatrix.cc,
     95     mhbase/MHn.cc, mhcalib/MHCalibrationChargeCam.cc, mhflux/MHDisp.cc,
     96     mimage/MCameraSmooth.cc, mimage/MConcentration.cc,
     97     mjobs/MJCalibration.cc, mjoptim/MJOptimize.cc,
     98     mjtrain/MJTrainSeparation.cc, mmuon/MMuonSearchPar.cc,
     99     mreport/MReportCurrents.cc, mreport/MReportRun.cc,
     100     mtools/MChisqEval.cc:
     101     - fixed some cast-away constness warnings
     102     - fixed some shadow declarations
     103
     104   * mhbase/MH.[h,cc]:
     105     - removed an old obsolete workaround (ProjectionX/Y)
     106
     107   * mimage/Makefile, mimage/IMageLinkDef.h:
     108     - removed MCameraSmooth
    91109
    92110
  • trunk/MagicSoft/Mars/mbase/MLog.cc

    r8805 r8957  
    11/* ======================================================================== *\
    2 ! $Name: not supported by cvs2svn $:$Id: MLog.cc,v 1.59 2008-01-11 10:37:01 tbretz Exp $
     2! $Name: not supported by cvs2svn $:$Id: MLog.cc,v 1.60 2008-06-14 15:55:50 tbretz Exp $
    33! --------------------------------------------------------------------------
    44!
     
    663663{
    664664    // gcc 3.2:
    665     char *txt = (char*)"logXXXXXX";
     665    const char *txt = "logXXXXXX";
    666666
    667667    TString n(fname ? fname : txt);
  • trunk/MagicSoft/Mars/mbase/MMath.cc

    r8673 r8957  
    11/* ======================================================================== *\
    2 ! $Name: not supported by cvs2svn $:$Id: MMath.cc,v 1.37 2007-08-17 10:53:48 tbretz Exp $
     2! $Name: not supported by cvs2svn $:$Id: MMath.cc,v 1.38 2008-06-14 15:55:50 tbretz Exp $
    33! --------------------------------------------------------------------------
    44!
     
    128128}
    129129
     130/*
     131Double_t MMath::SignificanceLiMaErr(Double_t s, Double_t b, Double_t alpha)
     132{
     133    Double_t S = SignificanceLiMa(s, b, alpha);
     134    if (S<0)
     135        return -1;
     136
     137    const Double_t sum = s+b;
     138
     139
     140    Double_t l = TMath::Log(s/sum*(alpha+1)/alpha)/TMath::Sqrt(2*S);
     141    Double_t m = TMath::Log(s/sum*(alpha+1)/alpha)/TMath::Sqrt(2*S);
     142
     143
     144    const Double_t sum = s+b;
     145
     146    if (s<0 || b<0 || alpha<=0)
     147        return -1;
     148
     149    const Double_t l = s==0 ? 0 : s*TMath::Log(s/sum*(alpha+1)/alpha);
     150    const Double_t m = b==0 ? 0 : b*TMath::Log(b/sum*(alpha+1)      );
     151
     152    return l+m<0 ? -1 : TMath::Sqrt((l+m)*2);
     153}
     154*/
     155
    130156// --------------------------------------------------------------------------
    131157//
     
    752778        static const Double_t sqrt3 = TMath::Sqrt(3.);
    753779
    754         const Double_t s = TMath::Sin(t)*sqrt3;
    755         const Double_t c = TMath::Cos(t);
    756 
    757         x1 = 2*sqrtq *      c  - a/3;
    758         x2 =  -sqrtq * (s + c) - a/3;
    759         x3 =   sqrtq * (s - c) - a/3;
     780        const Double_t sn = TMath::Sin(t)*sqrt3;
     781        const Double_t cs = TMath::Cos(t);
     782
     783        x1 = 2*sqrtq *       cs  - a/3;
     784        x2 =  -sqrtq * (sn + cs) - a/3;
     785        x3 =   sqrtq * (sn - cs) - a/3;
    760786
    761787        /* --- Easier to understand but slower ---
  • trunk/MagicSoft/Mars/mbase/MParContainer.cc

    r8724 r8957  
    359359        return kFALSE;
    360360
    361     /*const*/ TMethodCall *call = ((TDataMember*)member)->GetterMethod(); //FIXME: Root
     361    /*const*/ TMethodCall *call = const_cast<TDataMember*>(member)->GetterMethod(); //FIXME: Root
    362362    if (!call)
    363363    {
  • trunk/MagicSoft/Mars/mcalib/MCalibrateData.cc

    r8642 r8957  
    11/* ======================================================================== *\
    2 ! $Name: not supported by cvs2svn $:$Id: MCalibrateData.cc,v 1.70 2007-07-18 19:30:12 tbretz Exp $
     2! $Name: not supported by cvs2svn $:$Id: MCalibrateData.cc,v 1.71 2008-06-14 15:55:51 tbretz Exp $
    33! --------------------------------------------------------------------------
    44!
     
    558558                // Correct for global shifts in light emission
    559559                //
    560                 const MCalibrationChargePix &ugpix = (MCalibrationChargePix&)updatecam->GetAverageArea(0);
     560                const MCalibrationChargePix &ugpix = static_cast<const MCalibrationChargePix&>(updatecam->GetAverageArea(0));
    561561
    562562                const Float_t globmean = avpix.GetConvertedMean();
    563563                calibUpdate = globmean==0 ? 1 : ugpix.GetConvertedMean()/globmean;
    564564
    565                 MBadPixelsPix &ubad  = (MBadPixelsPix&)updatecam->GetAverageBadArea(0);
     565                const MBadPixelsPix &ubad  = static_cast<const MBadPixelsPix&>(updatecam->GetAverageBadArea(0));
    566566                if (ubad.IsUncalibrated(MBadPixelsPix::kChargeIsPedestal))
    567567                  {
     
    698698// shall be calibrated, respectively.
    699699//
    700 Int_t MCalibrateData::Calibrate(Bool_t data, Bool_t pedestal) const
    701 {
    702     if (!data && !pedestal)
     700Int_t MCalibrateData::Calibrate(Bool_t data, Bool_t doped) const
     701{
     702    if (!data && !doped)
    703703        return kTRUE;
    704704
     
    763763        } /* if (data) */
    764764     
    765         if (pedestal)
     765        if (doped)
    766766        {
    767767            TIter NextPed(&fPedestalCams);
  • trunk/MagicSoft/Mars/mcalib/MCalibrationPulseTimeCam.cc

    r8506 r8957  
    199199  Float_t area = cam[idx].GetA();
    200200
    201  if (area == 0)
     201  if (area == 0)
    202202    return kFALSE;
    203203
    204  MCalibrationPix &pix = (MCalibrationPix&)(*this)[idx];
     204  const MCalibrationPix &pix = static_cast<const MCalibrationPix&>((*this)[idx]);
    205205
    206206  switch (type)
  • trunk/MagicSoft/Mars/mhbase/MBinning.h

    r8679 r8957  
    9494    Int_t GetNumBins() const   { return fEdges.GetSize()-1; }
    9595
    96     Double_t *GetEdges() const { return (Double_t*)fEdges.GetArray(); }
     96    const Double_t *GetEdges() const { return fEdges.GetArray(); }
    9797    const TArrayD &GetEdgesD() const { return fEdges; }
    9898
  • trunk/MagicSoft/Mars/mhbase/MH.cc

    r8916 r8957  
    11/* ======================================================================== *\
    2 ! $Name: not supported by cvs2svn $:$Id: MH.cc,v 1.38 2008-06-02 17:21:24 tbretz Exp $
     2! $Name: not supported by cvs2svn $:$Id: MH.cc,v 1.39 2008-06-14 15:55:52 tbretz Exp $
    33! --------------------------------------------------------------------------
    44!
     
    554554    if (h->InheritsFrom(TH3::Class()) && x->InheritsFrom(TH3::Class()))
    555555    {
    556         SetBinning((TH3*)h, ((TH1*)x)->GetXaxis(), ((TH1*)x)->GetYaxis(), ((TH1*)x)->GetZaxis());
     556        SetBinning((TH3*)h, x->GetXaxis(), x->GetYaxis(), x->GetZaxis());
    557557        return;
    558558    }
     
    561561    if (h->InheritsFrom(TH2::Class()) && x->InheritsFrom(TH2::Class()))
    562562    {
    563         SetBinning((TH2*)h, ((TH1*)x)->GetXaxis(), ((TH1*)x)->GetYaxis());
     563        SetBinning((TH2*)h, x->GetXaxis(), x->GetYaxis());
    564564        return;
    565565    }
     
    568568    if (h->InheritsFrom(TH1::Class()) && x->InheritsFrom(TH1::Class()))
    569569    {
    570         SetBinning(h, ((TH1*)x)->GetXaxis());
     570        SetBinning(h, x->GetXaxis());
    571571        return;
    572572    }
     
    869869    // Draw first histogram
    870870    //
    871     TH1 *h1 = ((TH1&)hist1).DrawCopy();
     871    TH1 *h1 = hist1.DrawCopy();
    872872    gPad->SetBorderMode(0);
    873873    gPad->Update();
     
    921921    // Draw second histogram
    922922    //
    923     TH1 *h2 = ((TH1&)hist2).DrawCopy("sames");
     923    TH1 *h2 = hist2.DrawCopy("sames");
    924924    gPad->Update();
    925925
     
    12081208}
    12091209
    1210 void MH::ProjectionX(TH1D &dest, const TH2 &src, Int_t firstybin, Int_t lastybin)
    1211 {
    1212     //*-*-*-*-*Project a 2-D histogram into a 1-D histogram along X*-*-*-*-*-*-*
    1213     //*-*      ====================================================
    1214     //
    1215     //   The projection dest is always of the type TH1D.
    1216     //   The projection is made from the channels along the Y axis
    1217     //   ranging from firstybin to lastybin included.
    1218     //   By default, bins 1 to ny are included
    1219     //   When all bins are included, the number of entries in the projection
    1220     //   is set to the number of entries of the 2-D histogram, otherwise
    1221     //   the number of entries is incremented by 1 for all non empty cells.
    1222     //
    1223     //   if Sumw2() was called for dest, the errors are computed.
    1224     //
    1225     TAxis &axex = *((TH2&)src).GetXaxis();
    1226     TAxis &axey = *((TH2&)src).GetYaxis();
    1227 
    1228     const Int_t nx = axex.GetNbins();
    1229     const Int_t ny = axey.GetNbins();
    1230     if (firstybin < 0)
    1231         firstybin = 1;
    1232     if (lastybin > ny)
    1233         lastybin = ny;
    1234 
    1235     dest.Reset();
    1236     SetBinning(&dest, &axex);
    1237 
    1238     // Create the projection histogram
    1239     const Bool_t computeErrors = dest.GetSumw2N() ? 1 : 0;
    1240 
    1241     // Fill the projected histogram
    1242     for (Int_t binx=0; binx<=nx+1; binx++)
    1243     {
    1244         Double_t err2 = 0;
    1245         for (Int_t biny=firstybin; biny<=lastybin; biny++)
    1246         {
    1247             const Double_t cont = src.GetCellContent(binx,biny);
    1248             const Double_t err1 = src.GetCellError(binx,biny);
    1249             err2 += err1*err1;
    1250             if (cont)
    1251                 dest.Fill(axex.GetBinCenter(binx), cont);
    1252         }
    1253         if (computeErrors)
    1254             dest.SetBinError(binx, TMath::Sqrt(err2));
    1255     }
    1256     if (firstybin <=1 && lastybin >= ny)
    1257         dest.SetEntries(src.GetEntries());
    1258 }
    1259 
    1260 void MH::ProjectionY(TH1D &dest, const TH2 &src, Int_t firstxbin, Int_t lastxbin)
    1261 {
    1262     //*-*-*-*-*Project a 2-D histogram into a 1-D histogram along X*-*-*-*-*-*-*
    1263     //*-*      ====================================================
    1264     //
    1265     //   The projection dest is always of the type TH1D.
    1266     //   The projection is made from the channels along the Y axis
    1267     //   ranging from firstybin to lastybin included.
    1268     //   By default, bins 1 to ny are included
    1269     //   When all bins are included, the number of entries in the projection
    1270     //   is set to the number of entries of the 2-D histogram, otherwise
    1271     //   the number of entries is incremented by 1 for all non empty cells.
    1272     //
    1273     //   if Sumw2() was called for dest, the errors are computed.
    1274     //
    1275     TAxis &axex = *((TH2&)src).GetXaxis();
    1276     TAxis &axey = *((TH2&)src).GetYaxis();
    1277 
    1278     const Int_t nx = axex.GetNbins();
    1279     const Int_t ny = axey.GetNbins();
    1280     if (firstxbin < 0)
    1281         firstxbin = 1;
    1282     if (lastxbin > nx)
    1283         lastxbin = nx;
    1284 
    1285     dest.Reset();
    1286     SetBinning(&dest, &axey);
    1287 
    1288     // Create the projection histogram
    1289     const Bool_t computeErrors = dest.GetSumw2N() ? 1 : 0;
    1290 
    1291     // Fill the projected histogram
    1292     for (Int_t biny=0; biny<=ny+1; biny++)
    1293     {
    1294         Double_t err2 = 0;
    1295         for (Int_t binx=firstxbin; binx<=lastxbin; binx++)
    1296         {
    1297             const Double_t cont = src.GetCellContent(binx,biny);
    1298             const Double_t err1 = src.GetCellError(binx,biny);
    1299             err2 += err1*err1;
    1300             if (cont)
    1301                 dest.Fill(axey.GetBinCenter(biny), cont);
    1302         }
    1303         if (computeErrors)
    1304             dest.SetBinError(biny, TMath::Sqrt(err2));
    1305     }
    1306     if (firstxbin <=1 && lastxbin >= nx)
    1307         dest.SetEntries(src.GetEntries());
    1308 }
    1309 
    13101210// --------------------------------------------------------------------------
    13111211//
  • trunk/MagicSoft/Mars/mhbase/MH.h

    r8892 r8957  
    119119                              const char* name="ProjectArray", const char* title="Projected Array");
    120120   
    121     static void ProjectionX(TH1D &dest, const TH2 &src, Int_t firstybin=-1, Int_t lastybin=9999);
    122     static void ProjectionY(TH1D &dest, const TH2 &src, Int_t firstxbin=-1, Int_t lastxbin=9999);
    123 
    124121    static void GetRangeUser(const TH1 &hist, Axis_t &lo, Axis_t &hi);
    125122    static void GetRangeUserX(const TH1 &hist, Axis_t &lo, Axis_t &hi);
  • trunk/MagicSoft/Mars/mhbase/MHMatrix.cc

    r8907 r8957  
    783783
    784784    th1->cd(1);
    785     ((TH1&)hth).DrawCopy();   // real histogram before
     785    hth.DrawCopy();   // real histogram before
    786786
    787787    th1->cd(2);
    788     ((TH1&)hta).DrawCopy();   // histogram after
     788    hta.DrawCopy();   // histogram after
    789789
    790790    th1->cd(3);
    791     ((TH1&)hthd).DrawCopy();  // correction factors
     791    hthd.DrawCopy();  // correction factors
    792792
    793793    th1->cd(4);
    794     ((TH1&)thsh).DrawCopy();  // target
     794    thsh.DrawCopy();  // target
    795795}
    796796
     
    895895    SetBinning(&hthd, &thsh);
    896896    hthd.SetDirectory(NULL);
    897     hthd.Divide((TH1F*)&thsh, &hth, 1, 1);
     897    hthd.Divide(&thsh, &hth, 1, 1);
    898898
    899899    if (hthd.GetMaximum() <= 0)
  • trunk/MagicSoft/Mars/mhbase/MHn.cc

    r8893 r8957  
    584584    for (int i=0; i<fNum; i++)
    585585    {
    586         TString opt(fDrawOption[i]);
     586        TString dopt(fDrawOption[i]);
    587587        if (same)
    588             opt += " same";
     588            dopt += " same";
    589589
    590590        pad->cd(i+1);
    591         fHist[i]->Draw(opt);
    592     }
    593 }
     591        fHist[i]->Draw(dopt);
     592    }
     593}
  • trunk/MagicSoft/Mars/mhcalib/MHCalibrationChargeCam.cc

    r8552 r8957  
    11/* ======================================================================== *\
    2 ! $Name: not supported by cvs2svn $:$Id: MHCalibrationChargeCam.cc,v 1.60 2007-06-12 20:50:30 tbretz Exp $
     2! $Name: not supported by cvs2svn $:$Id: MHCalibrationChargeCam.cc,v 1.61 2008-06-14 15:55:53 tbretz Exp $
    33! --------------------------------------------------------------------------
    44!
     
    743743  }
    744744
    745   const Int_t npix = fRawEvt->GetNumPixels();
    746   const Int_t nhi  = fRunHeader->GetNumSamplesHiGain();
    747   for (int pixid=0; pixid<npix; pixid++)
     745  const Int_t num = fRawEvt->GetNumPixels();
     746  const Int_t nhi = fRunHeader->GetNumSamplesHiGain();
     747  for (int pixid=0; pixid<num; pixid++)
    748748  {
    749749      MHCalibrationChargePix &histhi = (MHCalibrationChargePix&)(*this)[pixid];
  • trunk/MagicSoft/Mars/mhflux/MHDisp.cc

    r8719 r8957  
    243243    const Int_t    bz     = fHist.GetZaxis()->FindFixBin(0);
    244244
    245     TH2 &bg = fHalf ? fHistBg1 : fHistBg2;
     245    TH2 &hbg = fHalf ? fHistBg1 : fHistBg2;
    246246
    247247    const Bool_t smear = fSmearing>0;
     
    249249    {
    250250        if (!fHistOff)
    251             bg.Fill(pos1.X(), pos1.Y(), w*w0);
     251            hbg.Fill(pos1.X(), pos1.Y(), w*w0);
    252252        else
    253253            fHist.Fill(pos1.X(), pos1.Y(), 0.0, w*w0);
     
    276276                        const Double_t rc = TMath::Exp(-dp)/normg;
    277277                        if (!fHistOff)
    278                             bg.AddBinContent(bg.GetBin(x, y), rc);
     278                            hbg.AddBinContent(hbg.GetBin(x, y), rc);
    279279                        else
    280280                            fHist.AddBinContent(fHist.GetBin(x, y, bz), rc);
     
    310310
    311311    if (!fHistOff)
    312         bg.SetEntries(bg.GetEntries()+1);
     312        hbg.SetEntries(hbg.GetEntries()+1);
    313313    else
    314314        fHist.SetEntries(fHist.GetEntries()+1);
  • trunk/MagicSoft/Mars/mimage/ImageLinkDef.h

    r6890 r8957  
    66
    77#pragma link C++ class MImgCleanStd+;
    8 #pragma link C++ class MCameraSmooth+;
    98
    109#pragma link C++ class MHillas+;
  • trunk/MagicSoft/Mars/mimage/MCameraSmooth.cc

    r8489 r8957  
    110110            const MGeomPix &gpix = (*fGeomCam)[i];
    111111
    112             const Int_t n = gpix.GetNumNeighbors();
     112            const Int_t nn = gpix.GetNumNeighbors();
    113113
    114114            Int_t num  = fUseCentralPixel ? 1 : 0;
     
    116116            errors[i]  = fUseCentralPixel ? pix.GetErrorPhot()  : 0;
    117117
    118             for (int j=0; j<n; j++)
     118            for (int j=0; j<nn; j++)
    119119            {
    120120                const UShort_t nid = gpix.GetNeighbor(j);
  • trunk/MagicSoft/Mars/mimage/MConcentration.cc

    r6855 r8957  
    111111
    112112   const UInt_t npix = evt.GetNumPixels();
    113    for (UInt_t i=0; i<npix; i++)
     113   for (UInt_t j=0; j<npix; j++)
    114114   {
    115        const MSignalPix &pix = evt[i];
     115       const MSignalPix &pix = evt[j];
    116116       if (!pix.IsPixelUsed())
    117117           continue;
    118118
    119         const Double_t nphot = pix.GetNumPhotons()* geom.GetPixRatio(i);
     119        const Double_t nphot = pix.GetNumPhotons()* geom.GetPixRatio(j);
    120120
    121121        // Get number of photons in the 8 most populated pixels
     
    134134                continue;
    135135
    136             for(int j=0;j<i-1;j++)
    137                 maxpix[7-j]=maxpix[6-j];                 // [#phot]
     136            for(int k=0;k<i-1;k++)
     137                maxpix[7-k]=maxpix[6-k];                 // [#phot]
    138138
    139139            maxpix[8-i]=nphot;
  • trunk/MagicSoft/Mars/mimage/Makefile

    r6890 r8957  
    2727
    2828SRCFILES = MImgCleanStd.cc \
    29            MCameraSmooth.cc \
    3029           MHillas.cc \
    3130           MHillasSrc.cc \
  • trunk/MagicSoft/Mars/mjobs/MJCalibration.cc

    r8948 r8957  
    917917          gPad->SetBorderMode(0);
    918918          gPad->SetTicks();
    919           MHCamera *obj10=(MHCamera*)disp30.DrawCopy("hist");
     919          MHCamera *obj12=(MHCamera*)disp30.DrawCopy("hist");
    920920          //
    921921          // for the datacheck, fix the ranges!!
    922922          //
    923           obj10->SetMinimum(fTimeOffsetMin);
    924           obj10->SetMaximum(fTimeOffsetMax);
     923          obj12->SetMinimum(fTimeOffsetMin);
     924          obj12->SetMaximum(fTimeOffsetMax);
    925925          //
    926926          // Set the datacheck sizes:
    927927          //
    928           FixDataCheckHist((TH1D*)obj10);
    929           obj10->SetStats(kFALSE);
     928          FixDataCheckHist((TH1D*)obj12);
     929          obj12->SetStats(kFALSE);
    930930          //
    931931          // set reference lines
    932932          //
    933           DisplayReferenceLines(obj10,5);
     933          DisplayReferenceLines(obj12,5);
    934934     
    935935          c5.cd(3);
     
    949949          gPad->SetBorderMode(0);
    950950          gPad->SetTicks();
    951           MHCamera *obj11=(MHCamera*)disp31.DrawCopy("hist");
     951          MHCamera *obj13=(MHCamera*)disp31.DrawCopy("hist");
    952952          //
    953953          // for the datacheck, fix the ranges!!
    954954          //
    955           obj11->SetMinimum(fTimeResolutionMin);
    956           obj11->SetMaximum(fTimeResolutionMax);
     955          obj13->SetMinimum(fTimeResolutionMin);
     956          obj13->SetMaximum(fTimeResolutionMax);
    957957          //
    958958          // Set the datacheck sizes:
    959959          //
    960           FixDataCheckHist((TH1D*)obj11);
    961           obj11->SetStats(kFALSE);
     960          FixDataCheckHist((TH1D*)obj13);
     961          obj13->SetStats(kFALSE);
    962962          //
    963963          // set reference lines
    964964          //
    965           DisplayReferenceLines(obj11,6);
     965          DisplayReferenceLines(obj13,6);
    966966         
    967967          c5.cd(4);
    968968          gPad->SetBorderMode(0);
    969           obj11->SetPrettyPalette();
    970           obj11->Draw();
     969          obj13->SetPrettyPalette();
     970          obj13->Draw();
    971971         
    972972          c5.cd(6);
     
    11511151            disp34.CamDraw(c32, 2, 2, 0);
    11521152
    1153             MHCalibrationCam *cam = (MHCalibrationCam*)plist.FindObject("MHCalibrationRelTimeCam");
     1153            MHCalibrationCam *ccam = (MHCalibrationCam*)plist.FindObject("MHCalibrationRelTimeCam");
    11541154           
    1155             for (Int_t sector=1;sector<cam->GetAverageSectors();sector++)
     1155            for (Int_t sector=1;sector<ccam->GetAverageSectors();sector++)
    11561156            {
    1157                 cam->GetAverageHiGainSector(sector).DrawClone("fourierevents");
    1158                 cam->GetAverageLoGainSector(sector).DrawClone("fourierevents");
     1157                ccam->GetAverageHiGainSector(sector).DrawClone("fourierevents");
     1158                ccam->GetAverageLoGainSector(sector).DrawClone("fourierevents");
    11591159            }
    11601160
  • trunk/MagicSoft/Mars/mjoptim/MJOptimize.cc

    r8704 r8957  
    413413Int_t MJOptimize::Minuit(TMinuit &minuit, const char *cmd) const
    414414{
    415     Int_t err;
     415    Int_t er;
    416416    Double_t tmp[2] = { fNumMaxCalls, fTolerance };
    417     minuit.mnexcm(cmd, tmp, 2, err);
    418 
    419     switch (err)
     417    minuit.mnexcm(cmd, tmp, 2, er);
     418
     419    switch (er)
    420420    {
    421421    case 0:
     
    461461    }
    462462
    463     return err;
     463    return er;
    464464}
    465465
  • trunk/MagicSoft/Mars/mjtrain/MJTrainSeparation.cc

    r8888 r8957  
    227227        const Float_t s = gx->Integral(1, y+1);
    228228        const Float_t b = hx->Integral(1, y+1);
    229         const Float_t sig1 = MMath::SignificanceLiMa(s+b, b);
    230         const Float_t sig2 = s<1 ? 0 : MMath::SignificanceLiMa(s+b, b)*TMath::Log10(s);
    231 
    232         gr3.SetPoint(y, h.GetYaxis()->GetBinLowEdge(y+2), sig1);
    233         gr4.SetPoint(y, h.GetYaxis()->GetBinLowEdge(y+2), sig2);
     229        const Float_t sg1 = MMath::SignificanceLiMa(s+b, b);
     230        const Float_t sg2 = s<1 ? 0 : MMath::SignificanceLiMa(s+b, b)*TMath::Log10(s);
     231
     232        gr3.SetPoint(y, h.GetYaxis()->GetBinLowEdge(y+2), sg1);
     233        gr4.SetPoint(y, h.GetYaxis()->GetBinLowEdge(y+2), sg2);
    234234    }
    235235    delete hx;
  • trunk/MagicSoft/Mars/mmuon/MMuonSearchPar.cc

    r8911 r8957  
    228228
    229229    // Perform Simplex minimization
    230     Int_t err;
     230    Int_t er;
    231231    Double_t tmp[2] = { 0, 0 };
    232     minuit.mnexcm("simplex", tmp, 2, err);
     232    minuit.mnexcm("simplex", tmp, 2, er);
    233233
    234234    // Get resulting parameters
  • trunk/MagicSoft/Mars/mreport/MReportCurrents.cc

    r4575 r8957  
    9494        pos += 4;
    9595
    96         const Int_t n=sscanf(hex, "%4x", &c);
    97         if (n!=1)
     96        const Int_t nn=sscanf(hex, "%4x", &c);
     97        if (nn!=1)
    9898        {
    9999            *fLog << warn << "WARNING - Reading hexadecimal DC information." << endl;
  • trunk/MagicSoft/Mars/mreport/MReportRun.cc

    r8955 r8957  
    8585    if (ver>=200805190)
    8686    {
    87         Int_t len;
    88         const Int_t n=sscanf(str.Data(), "%d %n", &file, &len);
    89         if (n!=1)
     87        const Int_t nn=sscanf(str.Data(), "%d %n", &file, &len);
     88        if (nn!=1)
    9089        {
    9190            *fLog << warn << "WARNING - Wrong number of arguments." << endl;
  • trunk/MagicSoft/Mars/mtools/MChisqEval.cc

    r8679 r8957  
    146146    const Double_t y2 = fData1 ? fData1->GetValue() : 0;
    147147
    148     const Double_t dy  = y2-y1;
    149     const Double_t err = fData1 ? y1*y2 : 1;
     148    const Double_t dy = y2-y1;
     149    const Double_t er = fData1 ? y1*y2 : 1;
    150150
    151     const Double_t w   = fWeight ? fWeight->GetVal() : 1;
     151    const Double_t w  = fWeight ? fWeight->GetVal() : 1;
    152152
    153     fChisq += w*dy*dy/err;
     153    fChisq += w*dy*dy/er;
    154154    fSumW  += w;
    155155
Note: See TracChangeset for help on using the changeset viewer.