TGraphAligned.h

#ifndef _TGRAPH_ALIGNED_H
#define _TGRAPH_ALIGNED_H

#include "TGraph.h" 

class TGraphAligned : public TGraph {


#define TGRAPH_ALIGNED_ALIGNMENT 32 
  typedef Double_t * aligned_double_v __attribute__((aligned (TGRAPH_ALIGNED_ALIGNMENT))); 
  //Have to reimplement all interesting constructors, unfortunately, since we
  //unfortunately cannot call the TGraph constructor and have it call our
  //implementation of CtorAllocate 

  public: 
    TGraphAligned(); 

    TGraphAligned(Int_t n); 
    TGraphAligned & operator=(const TGraphAligned&);
    TGraphAligned(Int_t n, const Double_t * x, const Double_t * y); 

    TGraphAligned(const TGraph  &); 


    double * getMoments(int i, double origin, double * moments = 0) const; 

    void dBize(double mindB=-200); 
    void undBize(); 


    void zeroMean(); 

    /* just copies the numbers from g*/ 
    void adopt(const TGraphAligned * g); 



    //Shift the values by a number of samples. 
    //by default will set any new values to zero, unless wrap is true
    void shift(int nsamples, bool zero = true); 

    void getMeanAndRMS(Double_t * mean, Double_t *  rms, Int_t istart = 0, Int_t iend = -1) const; 
    Double_t getSumV2(Int_t istart = 0, Int_t iend = -1) const; 

    Double_t peakVal(Int_t * location = 0, Int_t istart= 0, Int_t iend = -1, bool abs = false) const; 
    Double_t pk2pk(Int_t nthMax = 0, Int_t nthMin = 0,Int_t * location_max = 0, Int_t * location_min = 0, Int_t istart= 0, Int_t iend = -1) const; 

    TH1 * valueHist(int nbins = 100, const double * weights =0,  TH1 * out = 0 ) const; 

    void setPlottingLimits(double mult = 1.1, bool sym = true, double dt_from_peak =-1); 

    aligned_double_v GetX() const { return fX; } 
    aligned_double_v GetY() const { return fY; } 

    void setBelow(double val, bool abso=true, double to=0); 

    virtual ~TGraphAligned();

  protected: 

    virtual Double_t **AllocateAlignedArrays(Int_t Narrays, Int_t arraySize); 
    inline Double_t **Allocate(Int_t newsize) { return AllocateAlignedArrays(2, newsize); }
    virtual Bool_t CtorAllocate (void); 

    //needs to be reimplement to avoid delete[] on memalign-allocated pointer, which Valgrind will perhaps unjustly complain about
    virtual void CopyAndRelease(Double_t ** newarrays, Int_t ibegin, Int_t iend, Int_t obegin); 
    
    ClassDef(TGraphAligned,1); 

}; 


#endif