{
#include <stdlib.h>
#include <stdio.h>
#include <fcntl.h>

  FILE *fp, *fp1, *fp2;
  float ptmp,p[30];
  int i, iline=0;
  ntuple = new TNtuple("ntuple","NTUPLE","depth:gam:ep:em:mup:mum:had:charg:nuc:cer"); 
 
  fp = fopen("dataFile","r"); 
  while ( fscanf(fp,"%f",&ptmp) != EOF ){
  p[i++]=ptmp;
  if (i==10)
  {
   i=0;   iline++; 
   ntuple->Fill(p[0],p[1],p[2],p[3],p[4],p[5],p[6],p[7],p[8],p[9]); 
  } 
  gROOT->SetStyle("Plain");  
  gStyle->SetTitleBorderSize(0);
  gStyle->SetTitleX(0.0f);
  gStyle->SetTitleW(1.4f); 
  gStyle->SetTitleOffset(1.3,"y"); 
  gStyle->SetTitleOffset(1.0,"x"); 
  } 

 leg = new TLegend(0.2,0.6,0.3,0.8); 
 leg->SetFillStyle(0); 
 leg->SetBorderSize(0); 

 ntuple->SetMarkerColor(2); 
 ntuple->SetMarkerStyle(20); 
 ntuple->SetMarkerSize(1);
 ntuple.Draw("gam:depth","","lp");

 leg->AddEntry(ntuple,"","p"); 
 leg->Draw();
 gPad->Update();

c1->SetCanvasSize(600,400);  
c1->SetTitle("component");
htemp->GetXaxis()->SetTitle("g/cm^{2}"); 
htemp->GetYaxis()->SetTitle("N_{p}/m^{2}"); 
htemp->SetTitle("component");
gPad->Modified();
c1->Print("macroAlma.png");
}