#include "TROOT.h"
#include "TBuffer.h"
#include "TClass.h"
#include "THashList.h"
#include "TH2.h"
#include "TVirtualPad.h"
#include "TF2.h"
#include "TProfile.h"
#include "TRandom.h"
#include "TMatrixFBase.h"
#include "TMatrixDBase.h"
#include "THLimitsFinder.h"
#include "TError.h"
#include "TMath.h"
#include "TObjString.h"
#include "TObjArray.h"
#include "TVirtualHistPainter.h"
#include "snprintf.h"
 
//#include<TSpectrum>
//%using namespace std;

void M1back(){			
	char name[20]="cut3823";
	int maxnum = 3;
	int maxeventnum=10000;
	TChain tr("ge"); // treeNameという名のTTreeのchain
	tr.Add(Form("%s.root",name)); // treeNameなTTreeを含むROOTファイル
	for(int i=1;i<=maxnum;i++){ // 連番ファイルの連結の例
		if(i<=9){
			tr.Add(Form("%s_00%d.root",name,i));}
		else{
			tr.Add(Form("%s_0%d.root",name,i));}
	}
	tr.GetEntries();

	TH1D *h =new TH1D("h","h",16000,-0.5,15999.5);
//	ge->Draw("c22.E>>h","c29.E<100&&c17.E<100&&c18.E<100&&c19.E<100&&c20.E<100&&c21.E<100&&c16.E<100&&c23.E<100&&c24.E<100&&c25.E<100&&c26.E<100&&c27.E<100&&c28.E<100");
	ge->Draw("c16.E>>h","c29.E>1310&&c29.E<1350");
/*
	
	int tag_energy = 1173;
	int tag_crystal = 22;
	double total;
	double C_factor[14];
	double Num[14];
	TH1D *h[14];
	total = ge->GetEntries(Form("abs(c%d.E-%d)<10",tag_crystal,tag_energy));
	
	cout<<"Tag Crystal = c"<<tag_crystal<<" : Tag Energy ="<<tag_energy<<endl;
	cout<<"Total Numebr ="<<total<<endl;
	
	for(int i=0;i<14;i++){
		h[i] = new TH1D(Form("h%d",i+1),Form("c%d",i+1),10000,-0.5,9999.5);
		ge->Project(Form("h%d",i+1),Form("c%d.E",i+16),Form("abs(c%d.E-%d)<10",tag_crystal,tag_energy));
		Num[i] = h[i]->Integral(110,1500);
		C_factor[i] = Num[i]/total;
		cout<<" c"<<i+1<<" : number ="<<Num[i]<<", C factor="<<C_factor[i]<<endl;
	}
			
*/
	
	
	TSpectrum *s = new TSpectrum;
	s->Background(h,20,"");

	TH1D *h1 =new TH1D("h1","h1",16000,-0.5,15999.5);
	h1->Add(h);
	h1->Add(h_background,-1.);
}