|Message: RE: Brachytherapy example: TG43_relative_dose.C issue||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
Dear Lara, I executed the example and I did not find any problem. I list here the steps to undertake: 1) compile the brachytherapy with the analysis component on, e.g.: > cmake -DWITH_ANALYSIS_USE=ON -DGeant4_DIR=PATH/TO/geant4.10.03 PATH_TO/brachytherapy > make 2) run : ./Brachy FlexiSourceMacro.mac (I suggest to increase the number of events, e.g. 10^6 at least) The output of the simulation is brachytherapy.root. Check you have it. 3) root TG43_relative_dose.C Execute the root macro in the same directory where you have brachytherapy.root. This will generate geant4_dose.txt from reading brachytherapy.root. 4) Finally you can use the macro compare.C in the directory "comparison" to compare your simulation results to reference data (granero.txt). Just copy your geant4_dose.txt in the comparison directory and execute: root comparison.C To have a statistically meaningful comparison, you need to generate at least 300 million events in the simulation. I hope this helps you, Kind Regards, Susanna -----Original Message----- From: Lara Mason <firstname.lastname@example.org> Sent: Tuesday, 10 April 2018 10:30 PM To: email@example.com Subject: Brachytherapy example: TG43_relative_dose.C issue *** Discussion title: Medical Applications Hi experts, I'm attempting to use the ROOT macro TG43_relative_dose.C to compare to reference data. When I try to run it out of the box I get the error Processing TG43_relative_dose.C... input_line_13:2:6: error: member reference type 'TH2F *' is a pointer; did you mean to use '->'? (h20.GetXaxis()->GetBinCenter(((*(int*)0x7fff503611d0)))) ~~~^ -> libc++abi.dylib: terminating with uncaught exception of type libc++cling::CompilationException: Error evaluating expression libc++(h20.GetXaxis()->GetBinCenter(((*(int*)0x7fff503611d0)))) This is complaining about the lines Double_t xx_histo = h20.GetXaxis()->GetBinCenter(k); Double_t yy_histo = h20.GetYaxis()->GetBinCenter(m); Double_t edep_histo=h20.GetBinContent(k, m); However, when I switch them to Double_t xx_histo = h20->GetXaxis()->GetBinCenter(k); Double_t yy_histo = h20->GetYaxis()->GetBinCenter(m); Double_t edep_histo=h20->GetBinContent(k, m); And try root -l TG43_relative_dose.C again, it just runs indefinitely, returning no output and no errors and never completing. What am I doing wrong here? Thanks in advance for your help. Lara ------------------------------------------------------------- Visit this GEANT4 at hypernews.slac.stanford.edu message (to reply or unsubscribe) at: http://hypernews.slac.stanford.edu/HyperNews/geant4/get/medical_app/409.html
|Inline Depth:||Outline Depth:||Add message:|