#include <Riostream>

int Round(double x)	 
{
 int i;
 i=(int) floor(x);
 if ((x-i) < 0.5) return i; 
 return i+1;
}


void caen()
{
const int VectSize =256;

FILE *fpIn, *fpOut, *fpText;
char FileName[100]=" ";
char FileExt[12]="001";
char InFile [60];
double W1[VectSize],W2[VectSize];
double I1[VectSize],I2[VectSize];
long PHa[VectSize]={0};long QHa[VectSize]={0};long QHb[VectSize]={0};

int Ext,i,s1, s2;

int ;
double d1,d2,ds1,ds2;
float f1,f2,Qslow,Qtot,x,y;
char line[128];

TNtuple *ntuple=new TNtuple("ntuple"," ","x:y");
TFile *f=new TFile("tree1.root","RECREATE");
TH1F *h1=new TH1F("h1","x distribution",VectSize,0,VectSize);
TH1F *h2=new TH1F("h2","x distribution",VectSize,0,VectSize);
TH1F *h3=new TH1F("h3","x distribution",VectSize,0,VectSize*200);
TH1F *h4=new TH1F("h4","x distribution",VectSize,0,VectSize*200);
TH2F *h5=new TH2F("h5","x distribution",VectSize,0,VectSize*200,VectSize,0,VectSize*200);

	fpIn=fopen("Co938.txt","r");
	if(fpIn==NULL) return;
	Ext=0;

	while(!feof(fpIn)) {
	for(i=0;i<VectSize;i++) {
			fscanf(fpIn,"%f",&f1);
			W1[i]=f1;
	};


	
		ds1=0;
		for(i=0;i<70;i++) {
		ds1+=W1[i];
				
		};

		for(i=0;i<VectSize;i++) {
		W1[i]-=ds1/70.0;
		}								

		for(i=0;i<VectSize;i++) W2[i]=(-1)*W1[i];

		for(i=0;i<VectSize;i++){
		h1->Fill(i,W2[i]);
		h2->Fill(W2[i]);
		}

		Qtot=Qslow=0;
		
		for(i=80;i<180;i++){ 
		Qtot+=W2[i];}
		s1=Round(Qtot/100);
		if(s1>15&&s1<VectSize) QHa[s1]++;
		
		for(i=96;i<180;i++) {
		Qslow+=W2[i];}
		s2=Round(Qslow/100);
		if(s2>15&&s2<VectSize) QHb[s2]++;

		fpOut=fopen("PH.dat","w+");
		for(i=0;i<VectSize;i++) fprintf(fpOut,"%d %d %d\n",i,QHa[i],QHb[i]);
		fclose(fpOut);
	       
               for(i=0;i<VectSize;i++){
                ntuple->Fill(QHa[i],QHb[i]);
		h3->Fill(QHa[i]);
		h4->Fill(i,QHa[i]);
		h5->Fill(QHa[i],QHb[i]);
                }
		Ext++;
		if(Ext%100==0) {
		printf("Ext = %d\n",Ext);
		};
};
fclose(fpIn);
ntuple->Write();
}