///2012-09
///Source:http://root.cern.ch/phpBB3/viewtopic.php?f=3&t=12888
void GetHistFromFile ( TString A_file, TString A_outputDir ="./" ) {
  gStyle->SetPalette(1);
  string str(A_file.Data());

  Double_t energy = atof( (str.substr( str.find("Ekin") + 4, 4 )).c_str() );

  TString fParticleStringID = "unknown";
  if( fParticleStringID.CompareTo("unknown") == 0 ){
     Int_t particleStringLength =  (str.find("_Ekin") - (str.find("DeltaE-E_mcs_") + 13));
     fParticleStringID  = str.substr( str.find("DeltaE-E_mcs_") + 13, particleStringLength );
  }
  std::cout << "Particle: " << fParticleStringID.Data() << ", with energy " << energy << endl;

  TCanvas * c1 = new TCanvas("c1", Form("\\Delta E - E for %s energy = %f MeV", fParticleStringID.Data(), energy), 100, 100, 800, 600);

  TFile file0;
  file0.Open( A_file );

  hist1->SetTitle(Form ( "\\Delta E-E for %s energy = %4.0f MeV",fParticleStringID.Data(), energy ));
  hist1->Draw("COLSCAT");
  c1->SaveAs( Form("%s/ForGif_DeltaE-E_%s_Ekin%4.0f.png", A_outputDir.Data(), fParticleStringID.Data(), energy) );

}