void Core2(int ip=0) {
  TCanvas *c1(new TCanvas());

  TGraph *tg;
  tg=new TGraph(16);

  Double_t p[6],q[6];

  int e[16]={1,2,3,4,6,10,20,30,40,60,100,200,300,400,600,1000};

  for(unsigned i(0);i<11;i++) {
    if(i== 0) TFile::Open("SimAnalysisLayerTracks900066.root");
    if(i== 1) TFile::Open("SimAnalysisLayerTracks900067.root");
    if(i== 2) TFile::Open("SimAnalysisLayerTracks900068.root");
    if(i== 3) TFile::Open("SimAnalysisLayerTracks900069.root");
    if(i== 4) TFile::Open("SimAnalysisLayerTracks900070.root");
    if(i== 5) TFile::Open("SimAnalysisLayerTracks900071.root");
    if(i== 6) TFile::Open("SimAnalysisLayerTracks900057.root");
    if(i== 7) TFile::Open("SimAnalysisLayerTracks900054.root");
    if(i== 8) TFile::Open("SimAnalysisLayerTracks900051.root");
    if(i== 9) TFile::Open("SimAnalysisLayerTracks900048.root");
    if(i==10) TFile::Open("SimAnalysisLayerTracks900045.root");
    if(i==11) TFile::Open("SimAnalysisLayerTracks900046.root");
    if(i==12) TFile::Open("SimAnalysisLayerTracks900049.root");
    if(i==13) TFile::Open("SimAnalysisLayerTracks900052.root");
    if(i==14) TFile::Open("SimAnalysisLayerTracks900055.root");
    if(i==15) TFile::Open("SimAnalysisLayerTracks900058.root");

    double maxDensity(0.0);
    int maxLayer(-1);
    for(int j(0);j<30;j++) {
      double density(0.0);
      double density(TLayerAll->GetBinContent(j+1));
      if(density>maxDensity) {
      	maxDensity=density;
	maxLayer=j;
      }
    } 

    if(i<6)  maxDensity/=100000.0;
    if(i>=6 && i<11)  maxDensity/=10000.0;
    if(i>=11) maxDensity/=1000.0;

    if(ip==0) tg->SetPoint(i,e[i],maxDensity);
    else      tg->SetPoint(i,e[i],maxLayer);
  }

  if(ip==0) tg->SetTitle(";Photon energy (GeV);Maximum core density (particles/mm^{2})");
  else tg->SetTitle(";Photon energy (GeV);Layer with maximum core density");
  tg->SetMarkerStyle(20);
  tg->SetMarkerColor(2);
  tg->SetLineColor(2);

  tg->GetYaxis()->SetTitleOffset(1.32);


  tg->Draw("AP");
  
  c1->Update();
  if(ip==0) c1->Print("Core0.png");
  else      c1->Print("Core1.png");
}