Message: 3D histogram Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question 3D histogram 

Forum: Analysis
Date: 14 Sep, 2010
From: Ernesto Amato <Ernesto Amato>

Dear Geant4 experts, I am trying to introduce a three-dimensional histogram of deposited energy into the example TestEm12. It would be acceptable for me to have the same limits and number of bins along the three axes; so I modified in HistoManager.hh:

 class IHistogram3D;
...
  private:
...
    AIDA::IHistogram3D*      histo[MaxHisto];

  and, in HistoManager.cc:

void HistoManager::book()
...
    histo[1] = hf->createHistogram3D("EDEPOSIT", Nbins[1], Vmin[1],
Vmax[1], Nbins[1], Vmin[1], Vmax[1], Nbins[1], Vmin[1], Vmax[1]);
...

void HistoManager::FillHisto(G4int ih, G4double e, G4double weight)
...
  histo[1]->fill(1.,1.,1.,1.);

Obviously, this is only a first attempt, in which at every call of Fillhisto, the (1,1,1) voxel is incremented of weight=1.

Compilation and linking is successfull, but the execution is aborted with the messages:

# An unexpected error has been detected by HotSpot Virtual Machine:
#
#  SIGSEGV (0xb) at pc=0x08056d82, pid=16342, tid=3075435296
#
# Java VM: Java HotSpot(TM) Client VM (1.4.2_16-b05 mixed mode)
# Problematic frame:
# C  [voxel+0xed82]

I use geant4 version 9.1.p02 on Linux Ubuntu 8.04.

java version "1.6.0_0" OpenJDK Runtime Environment (build 1.6.0_0-b11) OpenJDK Server VM (build 1.6.0_0-b11, mixed mode)

If necessary, I can send also the error log file.

Thank you in advance for your precious help.

  Ernesto Amato (Italy)

 Add Message Add Message
to: "3D histogram"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]