#include <TROOT.h>
#include <TTree.h>
#include <TProfile.h>
#include <TH1D.h>

void getampsum()
{
   TH1D *ampsum = new TH1D("ampsum","amp sum",512,-10,10);
   TProfile *amp = new TProfile("amp","amp",512,-10,10,-10,10);
   TTree *tree = gROOT->FindObject("run26_8e");
   tree->Project("amp","ampsub[0]:xcm-14.30");
   for (int i=1; i<=512; ++i) {
      ampsum->SetBinContent(i,amp->GetBinContent(i));
   }
   amp->Reset();
   tree->Project("amp","ampsub[1]:xcm-14.30");
   for (int i=1; i<=512; ++i) {
      ampsum->SetBinContent(i,ampsum->GetBinContent(i)+amp->GetBinContent(i));
   }
   amp->Reset();
   tree->Project("amp","ampsub[2]:xcm-14.30");
   for (int i=1; i<=512; ++i) {
      ampsum->SetBinContent(i,ampsum->GetBinContent(i)+amp->GetBinContent(i));
   }
   amp->Reset();
   tree->Project("amp","ampsub[3]:xcm-14.30");
   for (int i=1; i<=512; ++i) {
      ampsum->SetBinContent(i,ampsum->GetBinContent(i)+amp->GetBinContent(i));
   }
   delete amp;
}