Double_t *in = new Double_t[nppw];
//TH1 *hb = NULL;
TH1 *hm = 0;
int g;
void LPF()
{
  //set in, do smoothing, (store new arr values each fft/ifft in filter script), repeat
  for (g=1; g<= 4; g++) {
    cout << "point 5" << endl;
    for (int h=0; h<nppw; h++) {
      in[h] = data[(g-1)*(nppw+1)+h];
    }
    cout << "point 4" << endl;
    gROOT->ProcessLine(".x DoLPFilterNewest.C");
    cout << "point 6" << endl;
  }
}