Changeset 13001 for trunk


Ignore:
Timestamp:
03/07/12 16:25:53 (13 years ago)
Author:
tbretz
Message:
Added AddCamDifference and weights to some Add-functions.
Location:
trunk/Mars/mhist
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/Mars/mhist/MHCamera.cc

    r12921 r13001  
    13591359// Call this function to add a MCamEvent on top of the present contents.
    13601360//
    1361 void MHCamera::AddMeanShift(const MCamEvent &event, Int_t type)
     1361void MHCamera::AddMeanShift(const MCamEvent &event, Int_t type, Stat_t w)
    13621362{
    13631363    if (fNcells<=1 || IsFreezed())
     
    13731373        {
    13741374            SetUsed(idx);
    1375             Fill(idx, val-mean); // FIXME: Slow!
     1375            Fill(idx, (val-mean)*w); // FIXME: Slow!
    13761376        }
    13771377    }
     
    13831383// Call this function to add a MCamEvent on top of the present contents.
    13841384//
    1385 void MHCamera::AddMedianShift(const MCamEvent &event, Int_t type)
     1385void MHCamera::AddMedianShift(const MCamEvent &event, Int_t type, Stat_t w)
    13861386{
    13871387    if (fNcells<=1 || IsFreezed())
     
    13971397        {
    13981398            SetUsed(idx);
    1399             Fill(idx, val-median); // FIXME: Slow!
     1399            Fill(idx, (val-median)*w); // FIXME: Slow!
    14001400        }
    14011401    }
     
    14071407// Call this function to add a MCamEvent on top of the present contents.
    14081408//
    1409 void MHCamera::AddCamContent(const MCamEvent &event, Int_t type)
     1409void MHCamera::AddCamContent(const MCamEvent &event, Int_t type, Stat_t w)
    14101410{
    14111411    if (fNcells<=1 || IsFreezed())
     
    14191419        {
    14201420            SetUsed(idx);
    1421             Fill(idx, val); // FIXME: Slow!
     1421            Fill(idx, val*w); // FIXME: Slow!
    14221422        }
     1423    }
     1424    fEntries++;
     1425}
     1426
     1427// ------------------------------------------------------------------------
     1428//
     1429void MHCamera::AddCamDifference(const MCamEvent &event1, const MCamEvent &event2, Int_t type)
     1430{
     1431    if (fNcells<=1 || IsFreezed())
     1432        return;
     1433
     1434    // FIXME: Security check missing!
     1435    for (Int_t idx=0; idx<fNcells-2; idx++)
     1436    {
     1437        Double_t val1=0;
     1438        if (!event1.GetPixelContent(val1, idx, *fGeomCam, type))
     1439            continue;
     1440
     1441        Double_t val2=0;
     1442        if (!event2.GetPixelContent(val2, idx, *fGeomCam, type))
     1443            continue;
     1444
     1445        SetUsed(idx);
     1446        Fill(idx, val1-val2); // FIXME: Slow!
    14231447    }
    14241448    fEntries++;
  • trunk/Mars/mhist/MHCamera.h

    r12921 r13001  
    163163
    164164    virtual Double_t GetPixContent(Int_t idx) const { return GetBinContent(idx+1); }
    165     virtual void     AddCamContent(const MCamEvent &evt, Int_t type=0);
     165    virtual void     AddCamDifference(const MCamEvent &evt, const MCamEvent &evt2, Int_t type=0);
     166    virtual void     AddCamContent(const MCamEvent &evt, Int_t type=0, Stat_t weight=1);
    166167    virtual void     AddCamContent(const MHCamera &evt, Int_t type=0);
    167168    virtual void     AddCamContent(const TArrayD &arr, const TArrayC *used=NULL);
     
    175176    virtual void     SetUsed(const TArrayC &arr);
    176177
    177     virtual void     AddMedianShift(const MCamEvent &evt, Int_t type=0);
    178     virtual void     AddMeanShift(const MCamEvent &evt, Int_t type=0);
     178    virtual void     AddMedianShift(const MCamEvent &evt, Int_t type=0, Stat_t weight=1);
     179    virtual void     AddMeanShift(const MCamEvent &evt, Int_t type=0, Stat_t weight=1);
    179180
    180181    virtual void     CntCamContent(const MCamEvent &evt, Double_t threshold, Int_t type=0, Bool_t isabove=kTRUE);
Note: See TracChangeset for help on using the changeset viewer.