bizhan bakhtiyari | 18 Apr 06:17 2015
Picon

help me PLZ

Hello to every user
i have  problems with my CT simulations
my spectrum is between 10-140 kev
and my activity is 16000000 bequerel
but when i want to get a image from VGate
Due to the high activity the simulation takes too much time to give a projection
and in high energy, projection is blurred and ambiguous
1-how can i  reduce the time?
2-how can i reduse the image blurring?
thanks for your help
<div><div dir="ltr">Hello to every user<div>i have &nbsp;problems with my CT simulations</div>
<div>my spectrum is between 10-140 kev</div>
<div>and my activity is 16000000 bequerel</div>
<div>but when i want to get a image from VGate</div>
<div>Due to the high activity the simulation takes too much time to give a projection</div>
<div>and in high energy, projection is blurred and ambiguous</div>
<div>1-how can i &nbsp;reduce the time?</div>
<div>2-how can i reduse the image blurring?</div>
<div>thanks for your help</div> 		 	   		  </div></div>
seta al-enazi | 18 Apr 00:40 2015
Picon

benchmarkPET

Dear Gaters
 i tried to execute benchmarkPET for Zr-89 only,
i wrote Zr-89 as an ion source as following:
/gate/source/addSource Zr89
/gate/source/Zr89/setActivity  1000 Bq
/gate/source/Zr89/gps/particle  ion
/gate/source/Zr89/gps/ion 40 89 0 0
/gate/source/Zr89/setForcedUnstableFlag true
/gate/source/Zr89/setForcedHalfLife 282240.0 s
/gate/source/Zr89/gps/energytype Mono
/gate/source/Zr89/gps/monoenergy 511. keV
/gate/source/Zr89/gps/type Volume
/gate/source/Zr89/gps/shape Cylinder
/gate/source/Zr89/gps/halfz 34.0 cm
/gate/source/Zr89/gps/radius .5 mm
/gate/source/Zr89/gps/angtype iso
/gate/source/Zr89/gps/centre 0. 0. 0. cm

and i didn't add O-15 source, then i wrote the benchPET.C file as following;
{
   gROOT->Reset();
   TFile f("test.root");

   TTree *Coincidences = (TTree*)gDirectory->Get("Coincidences");
   TTree *Hits = (TTree*)gDirectory->Get("Hits");
   TTree *Singles = (TTree*)gDirectory->Get("Singles");
   TTree *delay = (TTree*)gDirectory->Get("delay");

////Declaration of leaves types - TTree Coincidences
//  
   Int_t           RayleighCrystal1;
   Int_t           RayleighCrystal2;
   Int_t           RayleighPhantom1;
   Int_t           RayleighPhantom2;
   Char_t          RayleighVolName1[40];
   Char_t          RayleighVolName2[40];
   Float_t         axialPos;
   Char_t          comptVolName1[40];
   Char_t          comptVolName2[40];
   Int_t           compton1;
   Int_t           compton2;
   Int_t           crystalID1;
   Int_t           crystalID2;
   Int_t           comptonPhantom1;
   Int_t           comptonPhantom2;
   Float_t         energy1;
   Float_t         energy2;   
   Int_t           eventID1;
   Int_t           eventID2;
   Float_t         globalPosX1;
   Float_t         globalPosX2;
   Float_t         globalPosY1;
   Float_t         globalPosY2;      
   Float_t         globalPosZ1;
   Float_t         globalPosZ2;
   Int_t           layerID1;
   Int_t           layerID2;
   Int_t           moduleID1;
   Int_t           moduleID2;
   Float_t         rotationAngle;
   Int_t           rsectorID1;
   Int_t           rsectorID2;
   Int_t           runID;
   Float_t         sinogramS;
   Float_t         sinogramTheta;
   Int_t           sourceID1;
   Int_t           sourceID2;
   Float_t         sourcePosX1;
   Float_t         sourcePosX2;
   Float_t         sourcePosY1;
   Float_t         sourcePosY2;
   Float_t         sourcePosZ1;
   Float_t         sourcePosZ2;
   Int_t           submoduleID1;
   Int_t           submoduleID2;
   Double_t         time1;
   Double_t         time2;
   
   Float_t         zmin,zmax,z;
//   
//Set branch addresses - TTree Coincicences
//  
   Coincidences->SetBranchAddress("RayleighCrystal1",&RayleighCrystal1);
   Coincidences->SetBranchAddress("RayleighCrystal2",&RayleighCrystal2);
   Coincidences->SetBranchAddress("RayleighPhantom1",&RayleighPhantom1);
   Coincidences->SetBranchAddress("RayleighPhantom2",&RayleighPhantom2);
   Coincidences->SetBranchAddress("RayleighVolName1",&RayleighVolName1);
   Coincidences->SetBranchAddress("RayleighVolName2",&RayleighVolName2);
   Coincidences->SetBranchAddress("axialPos",&axialPos);
   Coincidences->SetBranchAddress("comptVolName1",&comptVolName1);
   Coincidences->SetBranchAddress("comptVolName2",&comptVolName2);
   Coincidences->SetBranchAddress("comptonCrystal1",&compton1);
   Coincidences->SetBranchAddress("comptonCrystal2",&compton2);
   Coincidences->SetBranchAddress("crystalID1",&crystalID1);
   Coincidences->SetBranchAddress("crystalID2",&crystalID2);
   Coincidences->SetBranchAddress("comptonPhantom1",&comptonPhantom1);
   Coincidences->SetBranchAddress("comptonPhantom2",&comptonPhantom2);
   Coincidences->SetBranchAddress("energy1",&energy1);
   Coincidences->SetBranchAddress("energy2",&energy2);   
   Coincidences->SetBranchAddress("eventID1",&eventID1);
   Coincidences->SetBranchAddress("eventID2",&eventID2);
   Coincidences->SetBranchAddress("globalPosX1",&globalPosX1);
   Coincidences->SetBranchAddress("globalPosX2",&globalPosX2);
   Coincidences->SetBranchAddress("globalPosY1",&globalPosY1);
   Coincidences->SetBranchAddress("globalPosY2",&globalPosY2);      
   Coincidences->SetBranchAddress("globalPosZ1",&globalPosZ1);
   Coincidences->SetBranchAddress("globalPosZ2",&globalPosZ2);
   Coincidences->SetBranchAddress("layerID1",&layerID1);
   Coincidences->SetBranchAddress("layerID2",&layerID2);
   Coincidences->SetBranchAddress("moduleID1",&moduleID1);
   Coincidences->SetBranchAddress("moduleID2",&moduleID2);
   Coincidences->SetBranchAddress("rotationAngle",&rotationAngle);
   Coincidences->SetBranchAddress("rsectorID1",&rsectorID1);
   Coincidences->SetBranchAddress("rsectorID2",&rsectorID2);
   Coincidences->SetBranchAddress("runID",&runID);
   Coincidences->SetBranchAddress("sinogramS",&sinogramS);
   Coincidences->SetBranchAddress("sinogramTheta",&sinogramTheta);
   Coincidences->SetBranchAddress("sourceID1",&sourceID1);
   Coincidences->SetBranchAddress("sourceID2",&sourceID2);
   Coincidences->SetBranchAddress("sourcePosX1",&sourcePosX1);
   Coincidences->SetBranchAddress("sourcePosX2",&sourcePosX2);
   Coincidences->SetBranchAddress("sourcePosY1",&sourcePosY1);
   Coincidences->SetBranchAddress("sourcePosY2",&sourcePosY2);
   Coincidences->SetBranchAddress("sourcePosZ1",&sourcePosZ1);
   Coincidences->SetBranchAddress("sourcePosZ2",&sourcePosZ2);
   Coincidences->SetBranchAddress("submoduleID1",&submoduleID1);
   Coincidences->SetBranchAddress("submoduleID2",&submoduleID2);
   Coincidences->SetBranchAddress("time1",&time1);
   Coincidences->SetBranchAddress("time2",&time2);
   
 

   /////////////STAT////////   
   gStyle -> SetStatW(0.28);
   gStyle -> SetStatH(0.3);
   gStyle -> SetStatColor(41);   
   gStyle -> SetStatX(1);
   gStyle -> SetStatY(1);   
   gStyle -> SetStatFont(42);
   gStyle->SetOptStat(0);
   gStyle->SetOptFit(0);
   /////////////////////////

   
   TH1F *DecayZr89 = new TH1F("Decroissance Zr-89","Zr-89 decay",60,0.,240.);
   TH1F *DetectPosAxial = new TH1F("DetectPosAxial","Axial detection position",52.,-208.,+208.);

   TH1F *AxialSensitivityDet = new TH1F("AxialSensitivityDet","Axial sensitivity",50,-200.,+200.);
   TH1F *AxialScattersDet = new TH1F("AxialScattersDet","Axial scatters distribution",50,-200.,+200.);
   TH1F *AxialProfileDet = new TH1F("AxialProfileDet","",50,-200.,+200.);
   TH1F *ScatterFractionAxialDet = new TH1F("ScatterFractionAxialDet","Axial scatter fraction",50,-200.,+200.);
   
   TH2F *DetectPos1 = new TH2F("DetectPos1","Transaxial detection position",252,-504.,+504.,252,-504.,+504.);
   TH2F *DetectPos2 = new TH2F("DetectPos2","Transaxial detection position",252,-504.,+504.,252,-504.,+504.);
   Int_t nentries = Coincidences->GetEntries();
   Int_t nbytes = 0, nbytesdelay = 0, nrandom = 0, nscatter = 0, ntrue = 0;
   
   
   Double_t Zr89Activity = 100000.;
   Double_t StartTime   = 0.;
   Double_t StopTime    = 240.;
   
   Double_t Zr89HalfLife = 282276;    
       
   Double_t Zr89DecayFactor = (1.0 - exp(-log(2.0)*(StopTime-StartTime)/Zr89HalfLife))/
                             (exp(+log(2.0)*StartTime/Zr89HalfLife)*log(2.0)/Zr89HalfLife*(StopTime-StartTime));
                    
   Double_t Zr89Decay = Zr89Activity * (StopTime-StartTime) * Zr89DecayFactor;                 


//
// Loop for each event in the TTree Coincidences
//
   for (Int_t i=0; i<nentries;i++) {
     if (fmod((double)i,10000.0) == 0.0) cout << ".";
     nbytes += Coincidences->GetEntry(i);
     if (eventID1 != eventID2) { // Random coincidence
       ++nrandom;
     } else {  // True coincidence
       if (runID == 0) { // First frame
         DetectPos1->Fill(globalPosX1,globalPosY1);
         DetectPos1->Fill(globalPosX2,globalPosY2);
         DetectPos1->Fill(sourcePosX1,sourcePosY1);
       }  else if (runID == 1) { // Second frame
         DetectPos2->Fill(globalPosX1,globalPosY1);
         DetectPos2->Fill(globalPosX2,globalPosY2);
         DetectPos2->Fill(sourcePosX1,sourcePosY1);
       }
       DetectPosAxial->Fill(globalPosZ1);
       DetectPosAxial->Fill(globalPosZ2);
       z = (globalPosZ1+globalPosZ2)/2.;
       if (comptonPhantom1 == 0 && comptonPhantom2 == 0 &&
           RayleighPhantom1 == 0 && RayleighPhantom2 == 0) {  // true unscattered coincidence
         AxialSensitivityDet->Fill(z);
     ntrue++;
       } else { // true scattered coincidence
         AxialScattersDet->Fill(z);
         nscatter++;
       }  
       AxialProfileDet->Fill(z);
     
        if ((sourceID1 == 0) && (sourceID2 == 0)) DecayZr89->Fill(time1);
     }  
   }
   cout << endl << endl;
   e1 = new TF1("e1","expo",0.,240.);
   //g1 = new TF1("g1","gaus",-5.,+5.);
   //Acolinea_Angle_Distribution_deg->Fit("g1");
 
   //Double_t ndecay = Ion_decay_time_s->GetEntries();
   cout << endl << endl;     
   cout << " *********************************************************************************** " << endl;
   cout << " *                                                                                 * " << endl;
   cout << " *   G A T E    P E T    B E N C H M A R K    R E S U L T S    A N A L Y S I S     * " << endl;
   cout << " *                                                                                 * " << endl;
   cout << " *********************************************************************************** " << endl;
   cout << endl << endl;     
   cout << " Acquisition start time = " << StartTime << " [s]" << endl;
   cout << " Acquisition stop time  = " << StopTime  << " [s]" << endl;
 
   cout << " Zr-89 decay factor = " << Zr89DecayFactor << endl;
   
   cout << " Zr-89 initial activity = " << Zr89Activity << " [Bq]" << endl;   
                       
   cout << " Zr-89 decays = " << Zr89Decay << endl;                    
   cout << " ==> Expected total number of decays during the acquisition is " << Zr89Decay << " +/- " << sqrt(Zr89Decay) << endl;   
  // cout << " There are " << ndecay << " recorded decays" << endl;
   cout << " There are " << ntrue << " true unscattered coincidences" << endl;
   cout << " There are " << nrandom << " random coincidences" << endl;
   cout << " There are " << nscatter << " scattered coincidences" << endl;
   cout << "  ==> there are " << nentries << " coincidences (true, scattered, and random)" << endl;   
   cout << "  ==> global scatter fraction = " << (float)nscatter/(float)(nentries-nrandom) << endl;
  // cout << "  ==> absolute sensitivity = " << 100.*(float)ntrue/ndecay << " % " << endl;
   Double_t p1 = e1->GetParameter(1);
   
   
   
   //Double_t p2 = g1->GetParameter(2);
   //cout << " Gamma acolinearity FWHM = " << p2*2.3548 << " degree (expected: 0.58)" << endl;
   c1 = new TCanvas("c1","GATE",3,28,970,632);
   Int_t pos=1;
   c1->Divide(3,2);
   c1->SetFillColor(0);
   c1->SetBorderMode(0);
   c1->cd(pos++);
   DetectPos1->Draw();
   //c1->cd(pos++);
   DetectPos2->SetMarkerColor(kRed);
   DetectPos2->Draw("same");
   tex = new TLatex(10.,200.,"Run 1");
   tex->SetTextColor(1);
   tex->SetLineWidth(2);
   tex->Draw();
   tex1 = new TLatex(10.,100.,"Run 2");
   tex1->SetTextColor(2);
   tex1->SetLineWidth(2);
   tex1->Draw();
   c1->cd(pos++);
   DetectPosAxial->Draw();
   c1->cd(pos++);
   DecayO15->SetLineColor(kRed);
   DecayO15->Draw();
   DecayZr89->SetLineColor(kBlue);
   DecayZr89->Draw("SAME");
   c1->cd(pos++);
   AxialSensitivityDet->Draw();
   c1->cd(pos++);
   ScatterFractionAxialDet->Divide(AxialScattersDet,AxialProfileDet,1.,1.,"");  
   ScatterFractionAxialDet->Draw();
   //c1->cd(pos++);
   //Acolinea_Angle_Distribution_deg->Draw();
   c1->cd(pos++);
   delay->Draw("time1");
   tex = new TLatex(30.,440.,"Delays: with coincidence sorter");
   tex->SetTextColor(1);
   tex->SetLineWidth(2);
   tex->Draw();
   tex1 = new TLatex(30.,410.,"Randoms: eventID1 != eventID2");
   tex1->SetTextColor(2);
   tex1->SetLineWidth(2);
   tex1->Draw();
   Coincidences->SetLineColor(2);
   Coincidences->Draw("time1","eventID1 != eventID2","same");

   
   c1->Update();   
   c1->SaveAs("benchmarkPET.gif");
}   
when i run benchmarkPET.C in root , i the output just includes tranaxial detection position and axial detection position
no output of axial sensitivity and axial scatter freaction

is it necessary to add O-15 source to the macro file and benchmarkPET.C file?
<div><div><div>
<div dir="ltr">Dear Gaters</div>
<div dir="ltr">&nbsp;i tried to execute benchmarkPET for Zr-89 only,</div>
<div dir="ltr">i wrote Zr-89 as an ion source as following:</div>
<div dir="ltr">/gate/source/addSource Zr89<br class="">/gate/source/Zr89/setActivity&nbsp; 1000 Bq<br class="">/gate/source/Zr89/gps/particle&nbsp; ion<br class="">/gate/source/Zr89/gps/ion 40 89 0 0<br class="">/gate/source/Zr89/setForcedUnstableFlag true<br class="">/gate/source/Zr89/setForcedHalfLife 282240.0 s<br class="">/gate/source/Zr89/gps/energytype Mono<br class="">/gate/source/Zr89/gps/monoenergy 511. keV<br class="">/gate/source/Zr89/gps/type Volume<br class="">/gate/source/Zr89/gps/shape Cylinder<br class="">/gate/source/Zr89/gps/halfz 34.0 cm<br class="">/gate/source/Zr89/gps/radius .5 mm<br class="">/gate/source/Zr89/gps/angtype iso<br class="">/gate/source/Zr89/gps/centre 0. 0. 0. cm</div>
<div dir="ltr"><br></div>
<div dir="ltr">and i didn't add O-15 source, then i wrote the benchPET.C file as following;</div>
<div dir="ltr">{<br class="">&nbsp;&nbsp; gROOT-&gt;Reset();<br class="">&nbsp;&nbsp; TFile f("test.root");<br class=""><br class="">&nbsp;&nbsp; TTree *Coincidences = (TTree*)gDirectory-&gt;Get("Coincidences");<br class="">&nbsp;&nbsp; TTree *Hits = (TTree*)gDirectory-&gt;Get("Hits");<br class="">&nbsp;&nbsp; TTree *Singles = (TTree*)gDirectory-&gt;Get("Singles");<br class="">&nbsp;&nbsp; TTree *delay = (TTree*)gDirectory-&gt;Get("delay");<br class=""><br class="">////Declaration of leaves types - TTree Coincidences<br class="">// &nbsp;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; RayleighCrystal1;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; RayleighCrystal2;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; RayleighPhantom1;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; RayleighPhantom2;<br class="">&nbsp;&nbsp; Char_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; RayleighVolName1[40];<br class="">&nbsp;&nbsp; Char_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; RayleighVolName2[40];<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; axialPos;<br class="">&nbsp;&nbsp; Char_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; comptVolName1[40];<br class="">&nbsp;&nbsp; Char_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; comptVolName2[40];<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; compton1;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; compton2;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; crystalID1;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; crystalID2;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; comptonPhantom1;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; comptonPhantom2;<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; energy1;<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; energy2;&nbsp; &nbsp;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; eventID1;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; eventID2;<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; globalPosX1;<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; globalPosX2;<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; globalPosY1;<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; globalPosY2;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; globalPosZ1;<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; globalPosZ2;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; layerID1;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; layerID2;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; moduleID1;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; moduleID2;<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; rotationAngle;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; rsectorID1;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; rsectorID2;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; runID;<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; sinogramS;<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; sinogramTheta;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; sourceID1;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; sourceID2;<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; sourcePosX1;<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; sourcePosX2;<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; sourcePosY1;<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; sourcePosY2;<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; sourcePosZ1;<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; sourcePosZ2;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; submoduleID1;<br class="">&nbsp;&nbsp; Int_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; submoduleID2;<br class="">&nbsp;&nbsp; Double_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; time1;<br class="">&nbsp;&nbsp; Double_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; time2;<br class="">&nbsp; &nbsp;<br class="">&nbsp;&nbsp; Float_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; zmin,zmax,z;<br class="">//&nbsp; &nbsp;<br class="">//Set branch addresses - TTree Coincicences<br class="">// &nbsp;<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("RayleighCrystal1",&amp;RayleighCrystal1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("RayleighCrystal2",&amp;RayleighCrystal2);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("RayleighPhantom1",&amp;RayleighPhantom1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("RayleighPhantom2",&amp;RayleighPhantom2);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("RayleighVolName1",&amp;RayleighVolName1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("RayleighVolName2",&amp;RayleighVolName2);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("axialPos",&amp;axialPos);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("comptVolName1",&amp;comptVolName1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("comptVolName2",&amp;comptVolName2);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("comptonCrystal1",&amp;compton1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("comptonCrystal2",&amp;compton2);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("crystalID1",&amp;crystalID1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("crystalID2",&amp;crystalID2);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("comptonPhantom1",&amp;comptonPhantom1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("comptonPhantom2",&amp;comptonPhantom2);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("energy1",&amp;energy1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("energy2",&amp;energy2);&nbsp; &nbsp;<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("eventID1",&amp;eventID1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("eventID2",&amp;eventID2);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("globalPosX1",&amp;globalPosX1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("globalPosX2",&amp;globalPosX2);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("globalPosY1",&amp;globalPosY1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("globalPosY2",&amp;globalPosY2);&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("globalPosZ1",&amp;globalPosZ1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("globalPosZ2",&amp;globalPosZ2);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("layerID1",&amp;layerID1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("layerID2",&amp;layerID2);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("moduleID1",&amp;moduleID1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("moduleID2",&amp;moduleID2);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("rotationAngle",&amp;rotationAngle);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("rsectorID1",&amp;rsectorID1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("rsectorID2",&amp;rsectorID2);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("runID",&amp;runID);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("sinogramS",&amp;sinogramS);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("sinogramTheta",&amp;sinogramTheta);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("sourceID1",&amp;sourceID1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("sourceID2",&amp;sourceID2);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("sourcePosX1",&amp;sourcePosX1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("sourcePosX2",&amp;sourcePosX2);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("sourcePosY1",&amp;sourcePosY1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("sourcePosY2",&amp;sourcePosY2);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("sourcePosZ1",&amp;sourcePosZ1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("sourcePosZ2",&amp;sourcePosZ2);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("submoduleID1",&amp;submoduleID1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("submoduleID2",&amp;submoduleID2);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("time1",&amp;time1);<br class="">&nbsp;&nbsp; Coincidences-&gt;SetBranchAddress("time2",&amp;time2);<br class="">&nbsp; &nbsp;<br class="">&nbsp; <br class=""><br class="">&nbsp;&nbsp; /////////////STAT////////&nbsp; &nbsp;<br class="">&nbsp;&nbsp; gStyle -&gt; SetStatW(0.28);<br class="">&nbsp;&nbsp; gStyle -&gt; SetStatH(0.3);<br class="">&nbsp;&nbsp; gStyle -&gt; SetStatColor(41);&nbsp; &nbsp;<br class="">&nbsp;&nbsp; gStyle -&gt; SetStatX(1);<br class="">&nbsp;&nbsp; gStyle -&gt; SetStatY(1);&nbsp; &nbsp;<br class="">&nbsp;&nbsp; gStyle -&gt; SetStatFont(42);<br class="">&nbsp;&nbsp; gStyle-&gt;SetOptStat(0);<br class="">&nbsp;&nbsp; gStyle-&gt;SetOptFit(0);<br class="">&nbsp;&nbsp; /////////////////////////<br class=""><br class="">&nbsp; &nbsp;<br class="">&nbsp;&nbsp; TH1F *DecayZr89 = new TH1F("Decroissance Zr-89","Zr-89 decay",60,0.,240.);<br class="">&nbsp;&nbsp; TH1F *DetectPosAxial = new TH1F("DetectPosAxial","Axial detection position",52.,-208.,+208.);<br class=""><br class="">&nbsp;&nbsp; TH1F *AxialSensitivityDet = new TH1F("AxialSensitivityDet","Axial sensitivity",50,-200.,+200.);<br class="">&nbsp;&nbsp; TH1F *AxialScattersDet = new TH1F("AxialScattersDet","Axial scatters distribution",50,-200.,+200.);<br class="">&nbsp;&nbsp; TH1F *AxialProfileDet = new TH1F("AxialProfileDet","",50,-200.,+200.);<br class="">&nbsp;&nbsp; TH1F *ScatterFractionAxialDet = new TH1F("ScatterFractionAxialDet","Axial scatter fraction",50,-200.,+200.);<br class="">&nbsp; &nbsp;<br class="">&nbsp;&nbsp; TH2F *DetectPos1 = new TH2F("DetectPos1","Transaxial detection position",252,-504.,+504.,252,-504.,+504.);<br class="">&nbsp;&nbsp; TH2F *DetectPos2 = new TH2F("DetectPos2","Transaxial detection position",252,-504.,+504.,252,-504.,+504.);<br class="">&nbsp;&nbsp; Int_t nentries = Coincidences-&gt;GetEntries();<br class="">&nbsp;&nbsp; Int_t nbytes = 0, nbytesdelay = 0, nrandom = 0, nscatter = 0, ntrue = 0;<br class="">&nbsp; &nbsp;<br class="">&nbsp; &nbsp;<br class="">&nbsp;&nbsp; Double_t Zr89Activity = 100000.;<br class="">&nbsp;&nbsp; Double_t StartTime&nbsp;&nbsp; = 0.;<br class="">&nbsp;&nbsp; Double_t StopTime&nbsp;&nbsp;&nbsp; = 240.;<br class="">&nbsp; &nbsp;<br class="">&nbsp;&nbsp; Double_t Zr89HalfLife = 282276;&nbsp;&nbsp; &nbsp;<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<br class="">&nbsp;&nbsp; Double_t Zr89DecayFactor = (1.0 - exp(-log(2.0)*(StopTime-StartTime)/Zr89HalfLife))/<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (exp(+log(2.0)*StartTime/Zr89HalfLife)*log(2.0)/Zr89HalfLife*(StopTime-StartTime)); <br class="">&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<br class="">&nbsp;&nbsp; Double_t Zr89Decay = Zr89Activity * (StopTime-StartTime) * Zr89DecayFactor;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<br class=""><br class=""><br class="">//<br class="">// Loop for each event in the TTree Coincidences<br class="">//<br class="">&nbsp;&nbsp; for (Int_t i=0; i&lt;nentries;i++) {<br class="">&nbsp;&nbsp;&nbsp;&nbsp; if (fmod((double)i,10000.0) == 0.0) cout &lt;&lt; ".";<br class="">&nbsp;&nbsp;&nbsp;&nbsp; nbytes += Coincidences-&gt;GetEntry(i);<br class="">&nbsp;&nbsp;&nbsp;&nbsp; if (eventID1 != eventID2) { // Random coincidence<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ++nrandom;<br class="">&nbsp;&nbsp;&nbsp;&nbsp; } else {&nbsp; // True coincidence<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if (runID == 0) { // First frame<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; DetectPos1-&gt;Fill(globalPosX1,globalPosY1);<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; DetectPos1-&gt;Fill(globalPosX2,globalPosY2);<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; DetectPos1-&gt;Fill(sourcePosX1,sourcePosY1);<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; }&nbsp; else if (runID == 1) { // Second frame<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; DetectPos2-&gt;Fill(globalPosX1,globalPosY1);<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; DetectPos2-&gt;Fill(globalPosX2,globalPosY2);<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; DetectPos2-&gt;Fill(sourcePosX1,sourcePosY1);<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; }<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; DetectPosAxial-&gt;Fill(globalPosZ1);<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; DetectPosAxial-&gt;Fill(globalPosZ2);<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; z = (globalPosZ1+globalPosZ2)/2.;<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if (comptonPhantom1 == 0 &amp;&amp; comptonPhantom2 == 0 &amp;&amp;<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; RayleighPhantom1 == 0 &amp;&amp; RayleighPhantom2 == 0) {&nbsp; // true unscattered coincidence<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; AxialSensitivityDet-&gt;Fill(z);<br class="">&nbsp;&nbsp; &nbsp; ntrue++;<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; } else { // true scattered coincidence<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; AxialScattersDet-&gt;Fill(z);<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nscatter++;<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; } &nbsp;<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; AxialProfileDet-&gt;Fill(z);<br class="">&nbsp;&nbsp;&nbsp; &nbsp;<br class="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if ((sourceID1 == 0) &amp;&amp; (sourceID2 == 0)) DecayZr89-&gt;Fill(time1);<br class="">&nbsp;&nbsp;&nbsp;&nbsp; } &nbsp;<br class="">&nbsp;&nbsp; }<br class="">&nbsp;&nbsp; cout &lt;&lt; endl &lt;&lt; endl;<br class="">&nbsp;&nbsp; e1 = new TF1("e1","expo",0.,240.);<br class="">&nbsp;&nbsp; //g1 = new TF1("g1","gaus",-5.,+5.);<br class="">&nbsp;&nbsp; //Acolinea_Angle_Distribution_deg-&gt;Fit("g1");<br class="">&nbsp; <br class="">&nbsp;&nbsp; //Double_t ndecay = Ion_decay_time_s-&gt;GetEntries();<br class="">&nbsp;&nbsp; cout &lt;&lt; endl &lt;&lt; endl;&nbsp;&nbsp;&nbsp; &nbsp;<br class="">&nbsp;&nbsp; cout &lt;&lt; " *********************************************************************************** " &lt;&lt; endl;<br class="">&nbsp;&nbsp; cout &lt;&lt; " *&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; * " &lt;&lt; endl;<br class="">&nbsp;&nbsp; cout &lt;&lt; " *&nbsp;&nbsp; G A T E&nbsp;&nbsp;&nbsp; P E T&nbsp;&nbsp;&nbsp; B E N C H M A R K&nbsp;&nbsp;&nbsp; R E S U L T S&nbsp;&nbsp;&nbsp; A N A L Y S I S&nbsp;&nbsp;&nbsp;&nbsp; * " &lt;&lt; endl;<br class="">&nbsp;&nbsp; cout &lt;&lt; " *&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; * " &lt;&lt; endl;<br class="">&nbsp;&nbsp; cout &lt;&lt; " *********************************************************************************** " &lt;&lt; endl;<br class="">&nbsp;&nbsp; cout &lt;&lt; endl &lt;&lt; endl;&nbsp;&nbsp;&nbsp; &nbsp;<br class="">&nbsp;&nbsp; cout &lt;&lt; " Acquisition start time = " &lt;&lt; StartTime &lt;&lt; " [s]" &lt;&lt; endl;<br class="">&nbsp;&nbsp; cout &lt;&lt; " Acquisition stop time&nbsp; = " &lt;&lt; StopTime&nbsp; &lt;&lt; " [s]" &lt;&lt; endl;<br class="">&nbsp; <br class="">&nbsp;&nbsp; cout &lt;&lt; " Zr-89 decay factor = " &lt;&lt; Zr89DecayFactor &lt;&lt; endl;<br class="">&nbsp; &nbsp;<br class="">&nbsp;&nbsp; cout &lt;&lt; " Zr-89 initial activity = " &lt;&lt; Zr89Activity &lt;&lt; " [Bq]" &lt;&lt; endl;&nbsp; &nbsp;<br class="">&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<br class="">&nbsp;&nbsp; cout &lt;&lt; " Zr-89 decays = " &lt;&lt; Zr89Decay &lt;&lt; endl;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<br class="">&nbsp;&nbsp; cout &lt;&lt; " ==&gt; Expected total number of decays during the acquisition is " &lt;&lt; Zr89Decay &lt;&lt; " +/- " &lt;&lt; sqrt(Zr89Decay) &lt;&lt; endl;&nbsp; &nbsp;<br class="">&nbsp; // cout &lt;&lt; " There are " &lt;&lt; ndecay &lt;&lt; " recorded decays" &lt;&lt; endl;<br class="">&nbsp;&nbsp; cout &lt;&lt; " There are " &lt;&lt; ntrue &lt;&lt; " true unscattered coincidences" &lt;&lt; endl;<br class="">&nbsp;&nbsp; cout &lt;&lt; " There are " &lt;&lt; nrandom &lt;&lt; " random coincidences" &lt;&lt; endl;<br class="">&nbsp;&nbsp; cout &lt;&lt; " There are " &lt;&lt; nscatter &lt;&lt; " scattered coincidences" &lt;&lt; endl;<br class="">&nbsp;&nbsp; cout &lt;&lt; "&nbsp; ==&gt; there are " &lt;&lt; nentries &lt;&lt; " coincidences (true, scattered, and random)" &lt;&lt; endl;&nbsp; &nbsp;<br class="">&nbsp;&nbsp; cout &lt;&lt; "&nbsp; ==&gt; global scatter fraction = " &lt;&lt; (float)nscatter/(float)(nentries-nrandom) &lt;&lt; endl;<br class="">&nbsp; // cout &lt;&lt; "&nbsp; ==&gt; absolute sensitivity = " &lt;&lt; 100.*(float)ntrue/ndecay &lt;&lt; " % " &lt;&lt; endl;<br class="">&nbsp;&nbsp; Double_t p1 = e1-&gt;GetParameter(1);<br class="">&nbsp; &nbsp;<br class="">&nbsp; &nbsp;<br class="">&nbsp; &nbsp;<br class="">&nbsp;&nbsp; //Double_t p2 = g1-&gt;GetParameter(2);<br class="">&nbsp;&nbsp; //cout &lt;&lt; " Gamma acolinearity FWHM = " &lt;&lt; p2*2.3548 &lt;&lt; " degree (expected: 0.58)" &lt;&lt; endl;<br class="">&nbsp;&nbsp; c1 = new TCanvas("c1","GATE",3,28,970,632);<br class="">&nbsp;&nbsp; Int_t pos=1;<br class="">&nbsp;&nbsp; c1-&gt;Divide(3,2);<br class="">&nbsp;&nbsp; c1-&gt;SetFillColor(0);<br class="">&nbsp;&nbsp; c1-&gt;SetBorderMode(0);<br class="">&nbsp;&nbsp; c1-&gt;cd(pos++);<br class="">&nbsp;&nbsp; DetectPos1-&gt;Draw();<br class="">&nbsp;&nbsp; //c1-&gt;cd(pos++);<br class="">&nbsp;&nbsp; DetectPos2-&gt;SetMarkerColor(kRed);<br class="">&nbsp;&nbsp; DetectPos2-&gt;Draw("same");<br class="">&nbsp;&nbsp; tex = new TLatex(10.,200.,"Run 1");<br class="">&nbsp;&nbsp; tex-&gt;SetTextColor(1);<br class="">&nbsp;&nbsp; tex-&gt;SetLineWidth(2);<br class="">&nbsp;&nbsp; tex-&gt;Draw();<br class="">&nbsp;&nbsp; tex1 = new TLatex(10.,100.,"Run 2");<br class="">&nbsp;&nbsp; tex1-&gt;SetTextColor(2);<br class="">&nbsp;&nbsp; tex1-&gt;SetLineWidth(2);<br class="">&nbsp;&nbsp; tex1-&gt;Draw();<br class="">&nbsp;&nbsp; c1-&gt;cd(pos++);<br class="">&nbsp;&nbsp; DetectPosAxial-&gt;Draw();<br class="">&nbsp;&nbsp; c1-&gt;cd(pos++);<br class="">&nbsp;&nbsp; DecayO15-&gt;SetLineColor(kRed);<br class="">&nbsp;&nbsp; DecayO15-&gt;Draw();<br class="">&nbsp;&nbsp; DecayZr89-&gt;SetLineColor(kBlue);<br class="">&nbsp;&nbsp; DecayZr89-&gt;Draw("SAME");<br class="">&nbsp;&nbsp; c1-&gt;cd(pos++);<br class="">&nbsp;&nbsp; AxialSensitivityDet-&gt;Draw();<br class="">&nbsp;&nbsp; c1-&gt;cd(pos++);<br class="">&nbsp;&nbsp; ScatterFractionAxialDet-&gt;Divide(AxialScattersDet,AxialProfileDet,1.,1.,""); &nbsp;<br class="">&nbsp;&nbsp; ScatterFractionAxialDet-&gt;Draw();<br class="">&nbsp;&nbsp; //c1-&gt;cd(pos++);<br class="">&nbsp;&nbsp; //Acolinea_Angle_Distribution_deg-&gt;Draw();<br class="">&nbsp;&nbsp; c1-&gt;cd(pos++);<br class="">&nbsp;&nbsp; delay-&gt;Draw("time1");<br class="">&nbsp;&nbsp; tex = new TLatex(30.,440.,"Delays: with coincidence sorter");<br class="">&nbsp;&nbsp; tex-&gt;SetTextColor(1);<br class="">&nbsp;&nbsp; tex-&gt;SetLineWidth(2);<br class="">&nbsp;&nbsp; tex-&gt;Draw();<br class="">&nbsp;&nbsp; tex1 = new TLatex(30.,410.,"Randoms: eventID1 != eventID2");<br class="">&nbsp;&nbsp; tex1-&gt;SetTextColor(2);<br class="">&nbsp;&nbsp; tex1-&gt;SetLineWidth(2);<br class="">&nbsp;&nbsp; tex1-&gt;Draw();<br class="">&nbsp;&nbsp; Coincidences-&gt;SetLineColor(2);<br class="">&nbsp;&nbsp; Coincidences-&gt;Draw("time1","eventID1 != eventID2","same"); <br class=""><br class="">&nbsp; &nbsp;<br class="">&nbsp;&nbsp; c1-&gt;Update();&nbsp; &nbsp;<br class="">&nbsp;&nbsp; c1-&gt;SaveAs("benchmarkPET.gif");<br class="">}&nbsp;&nbsp;&nbsp; <br>
</div>
<div dir="ltr">when i run benchmarkPET.C in root , i the output just includes tranaxial detection position and axial detection position <br>
</div>
<div dir="ltr">no output of axial sensitivity and axial scatter freaction</div>
<div dir="ltr"><br></div>
<div dir="ltr">is it necessary to add O-15 source to the macro file and benchmarkPET.C file?<br>
</div>
</div></div></div>
Débora Salvado | 17 Apr 20:54 2015
Picon

Inconsistent energy spectra with multiple sources

Dear GATE users

I have a SPECT setup as exemplified in the attached file setup.png

I run different simulations:
1) I123 (159 and 529keV) 185MBq in the patient volume
2) Lu176 (88, 200 and 300keV) 5.6MBq in the LSO volume
3) I123 185MBq and Lu176 5.6MBq in the respective volumes

In the file energySpectra.png, I show you the energy spectra obtained when running every simulation for 1 minute. Can someone explain me why the dashed and the black energy spectra are different?
The dashed spectrum is just a sum of the red and green data series.

In my opinion, the proportion of the peaks in the black spectrum seems correct (simulation with both sources), but then why doesn't it correspond to the sum of the counts obtained in separate simulations with only one of the sources?

Any suggestions?

Regards,
Débora

Debora Salvado, MSc
PhD Student
University College London
INM, UCLH,
(T-5) 235 Euston Road
London NW1 2BU
UK
<div><div dir="ltr">
<div>
<div>Dear GATE users<br><br>
</div>
<div>I have a SPECT setup as exemplified in the attached file setup.png<br><br>
</div>
<div>I run different simulations:<br>1) I123 (159 and 529keV) 185MBq in the patient volume<br>
</div>
<div>2) Lu176 (88, 200 and 300keV) 5.6MBq in the LSO volume<br>
</div>
<div>3) I123 185MBq and Lu176 5.6MBq in the respective volumes<br><br>
</div>
<div>In the file energySpectra.png, I show you the energy spectra obtained when running every simulation for 1 minute. Can someone explain me why the dashed and the black energy spectra are different?<br>
</div>The dashed spectrum is just a sum of the red and green data series.<br><br>
</div>
<div>In my opinion, the proportion of the peaks in the black spectrum seems correct (simulation with both sources), but then why doesn't it correspond to the sum of the counts obtained in separate simulations with only one of the sources?<br><br>
</div>
<div>Any suggestions?<br><br>
</div>
<div>Regards,<br>D&eacute;bora<br><br>
</div>
<div><div><div>
<div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>
<span></span><span><span></span><span>Debora Salvado, MSc</span></span><br><div><span><span>PhD Student<br></span></span></div>
<span><span>University College London<br>INM, UCLH, </span></span><span><span><span><span><span> (T-5) </span></span></span>235 Euston Road<br></span></span><span><span>London NW1 2BU<br></span></span><span><span>UK</span></span><span><span></span></span><br>
</div></div></div></div></div></div>
</div></div></div>
</div></div>
Tommaso Stefano Carzaniga | 17 Apr 17:59 2015
Picon

G8 PET-CT from Perkin Elmer

Hello everyone,

I'm a newbie user of OpenGate and I'm setting up a simulation of a small animal PET system for 43Sc.
Before starting to model the G8 PET-CT from Perkin Elmer, I would like to ask if somebody already worked with this system.

Thanks a lot for your help!

Tommaso

--
Tommaso Stefano CARZANIGA

PhD candidate at LHEP - Universität Bern
Student at MAS in Medical Physics - ETHZ
Alumnus UNITECH International Society - ETHZ
Alumnus T.I.M.E. - Top Industrial Manager for Europe - ECLille
Alumnus Politecnico di Milano - Physics Engineering -

<div><div dir="ltr">Hello everyone,<div><br></div>
<div>I'm a newbie user of OpenGate and I'm setting up a simulation of a small animal PET system for 43Sc.</div>
<div>Before starting to model the&nbsp;G8 PET-CT from Perkin Elmer, I would like to ask if somebody already worked with this system.</div>
<div><br></div>
<div>Thanks a lot for your help!</div>
<div><br></div>
<div>Tommaso<br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr">Tommaso Stefano CARZANIGA<br><br>PhD candidate at LHEP - Universit&auml;t Bern<br><div>Student at MAS in Medical Physics -&nbsp;<span>ETHZ</span>
</div>
<div><div><div>Alumnus UNITECH International Society - ETHZ</div></div></div>
<div>
<span>Alumnus T.I.M.E. - Top Industrial Manager for Europe - ECLille</span><br>
</div>
<div><div></div></div>
<div>
<span>Alumnus Politecnico di Milano - Physics Engineering -</span><br>
</div>
<div><span><br></span></div>
</div></div></div></div>
</div>
</div></div>
abdallah smida | 14 Apr 19:09 2015
Picon

Re: Gate-users Digest, Vol 107, Issue 32

Please don't send me your mails in future




Le Mardi 14 avril 2015 17h19, "gate-users-request-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhz460WcJ8gAc@public.gmane.orgg" <gate-users-request-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhSoIsB4E12gc@public.gmane.org> a écrit :


Send Gate-users mailing list submissions to
    gate-users-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhSoIsB4E12gc@public.gmane.org

To subscribe or unsubscribe via the World Wide Web, visit
    http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
or, via email, send a message with subject or body 'help' to
    gate-users-request-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhFnp2FX2i8LQ@public.gmane.orgn.org

You can reach the person managing the list at
    gate-users-owner-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhSoIsB4E12gc@public.gmane.org

When replying, please edit your Subject line so it is more specific
than "Re: Contents of Gate-users digest..."


Today's Topics:

  1. Re: ion source timestamp bug, Gate 7.1 (Gergely Patay)
  2. Re: Gate 7.1 build broken after updating to    Geant4.10.01.p01
      (Alex Vergara Gil)
  3. Re: Gate 7.1 build broken after updating to    Geant4.10.01.p01
      (Marc Verderi)
  4. Geant4 and Eclipse (imen elhamrouni)
  5. Re: Root not working in VGate 3.0 (Mathieu DUPONT)


----------------------------------------------------------------------

Message: 1
Date: Tue, 14 Apr 2015 15:13:32 +0200
From: Gergely Patay <gergely.patay-wwqpr9UM3B/QT0dZR+AlfA@public.gmane.org>
To: Gate User <gate-users-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEh9Sl0XMgJPXI@public.gmane.orgon.org>
Subject: Re: [Gate-users] ion source timestamp bug, Gate 7.1
Message-ID: <552D127C.7040800-wwqpr9UM3B/QT0dZR+AlfA@public.gmane.org>
Content-Type: text/plain; charset="utf-8"; format=flowed

Dear Gaters,

to answer my own question: the problem is indeed deep inside Geant4, and as
usual, originally this is "not a bug, but a feature".

In GateVSource.cc:268 we (try to) force life time of the source particle to zero:
// Force life time to 0, time is managed by GATE not G4
GetParticleDefinition()->SetPDGLifeTime(0);

However, from v10.0+ when Geant4 reads back the life time of a general ion
source, it does *not* retrieve the previously set value, but an (otherwise
physically correct) value from an external data file
(G4ENSDFSTATE1.0/ENSDFSTATE.dat) :

G4ParticleDefinition.icc:45
inline
  G4double G4ParticleDefinition::GetPDGLifeTime() const
{
  if(IsGeneralIon())
  { return GetIonLifeTime(); }
  else
  { return thePDGLifeTime; }
}

Unfortunately, the underlying isGeneralIon variable is not directly accessible
from outside (i.e. from Gate).

A minimal workaround might be as follows:

1) change visibility of isGeneralIon from protected to public:
geant4/particles/management/include/G4ParticleDefinition.hh:402
-  protected:
+  public:
        G4bool isGeneralIon;

(this might be modified *without* recompiling Geant4)

2) add the following line to gate/source/physics/src/GateVSource.cc:269
          GetParticleDefinition()->SetPDGLifeTime(0);
+        GetParticleDefinition()->isGeneralIon=false;

3) recompile Gate with the modified header file.

Since this modification affects Geant4 code too, I cannot provide a simple patch
file yet.

best regards,
  Gergely

  --
Gergely Patay
Mediso Ltd.



On 2015-04-01 13:38, Gergely Patay wrote:
> Dear Gaters,
>
> I installed the latest version of Gate from git. Everything went fine
> (benchmarks etc), until I tried to use the 'generic ion' source for a PET
> simulation.
>
> The timestamps of the hits (and consequently, of the singles and coincidences)
> are weird. They should be distributed uniformly in the [0,1] s interval (as seen
> on the ion_timestamp_gate7.0.png figure), but in 7.1 they are distributed
> exponentially over a much-much larger time interval with an exponent which
> equals to the decay rate of the ion, in this case Na-22.
>
> The macro files were identical, and the corresponding software versions were
> - Gate7.0, Geant4.9.6.p03, root5.34.18
> - Gate7.1, Geant4.10.01, root5.34.26
>
> This bug is actually the very same to that was reported on March 31, 2011 for
> Gate 6.0.1
> (http://lists.opengatecollaboration.org/mailman/private/gate-users/2011-March/004895.html)
> which somehow disappeared from 7.0 but re-emerged in 7.1.
>
> The bug does not affect the backtoback or e+ sources, and energy spectra are fine.
> It might be buried somewhere deep in Geant4.
>
> Do you have any idea what went wrong?
>
> thanks,
>  Gergely
>
> --
> Gergely Patay
> Mediso Ltd.
>
>
> _______________________________________________
> Gate-users mailing list
> Gate-users-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhSoIsB4E12gc@public.gmane.org
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>


------------------------------

Message: 2
Date: Tue, 14 Apr 2015 09:16:01 -0400
From: Alex Vergara Gil <alexvergaragil-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org>
To: REUZE SYLVAIN <SYLVAIN.REUZE-BQgGim59Xlv5SJMB2oCDiQ@public.gmane.org>
Cc: Gate Users <gate-users-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhSoIsB4E12gc@public.gmane.org>
Subject: Re: [Gate-users] Gate 7.1 build broken after updating to
    Geant4.10.01.p01
Message-ID:
    <CAMF2x=nDWVfRCzt5sn0RuhkwBpEme594sZpB70+Z_N5pbjtO4Q-JsoAwUIsXosN+BqQ9rBEUg@public.gmane.org>
Content-Type: text/plain; charset=UTF-8

Indeed! In the G4 page there is no link to the source code of 10.1
version without the patch01.

This is dangerous for Gate, shall it be created a mirror image of the
current working G4 source code in the
http://www.opengatecollaboration.org page??? I think this would be a
temporary solution at least for now.

Regards
Alex


2015-04-14 9:06 GMT-04:00, REUZE SYLVAIN <SYLVAIN.REUZE-BQgGim59Xlv5SJMB2oCDiQ@public.gmane.org>:
> Hello,
>
>
>
> I have the same compilation errors with Geant4 10.01 p01. I'd like to
> install Gate 7.1 in order to solve another problem, but I can't find the
> link to download Geant4 10.01, without the patch.
>
> Could you help me ?
>
>
>
> Best regards,
>
> Sylvain Reuz?
> Stagiaire physicien
> D?partement de Radioth?rapie
> Service de physique m?dicale
>
> T?l. : 01 42 11 54 06
> www.gustaveroussy.fr<http://www.gustaveroussy.fr>
> [cid:image001.gif <at> 01D076C4.8D1D9C60]
>
>
>
>
>
>
>
> -----Message d'origine-----
> De : gate-users-bounces-ZwoEplunGu2uWGBtOU8ma2ZaI6lBa+6R@public.gmane.orgboration.org
> [mailto:gate-users-bounces-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhSoIsB4E12gc@public.gmane.org] De la part de
> Alex Vergara Gil
> Envoy? : lundi 13 avril 2015 13:12
> ? : Marc Chamberland
> Cc : Gate Users
> Objet : Re: [Gate-users] Gate 7.1 build broken after updating to
> Geant4.10.01.p01
>
>
>
> For what it's worth, I can reproduce the compilation errors, too.
>
> There might be a compatibility issue with the
>
> G4SPSRandomGenerator.hh:156:7:
>
> error: field of type 'G4Cache<G4SPSRandomGenerator::a_check>' has
>
> private copy constructor
>
> class G4SPSRandomGenerator {
>
>
>
> I'm compiling with:
>
>
>
> GCC version 4.8
>
>
>
> Alex
>
>
>
> 2015-04-10 11:53 GMT-04:00, Marc Chamberland
>
> <MarcChamberland <at> cmail.carleton.ca>:
>
>> For what it's worth, I can reproduce the compilation errors, too. I'm
>
>> compiling with:
>
>>
>
>> Apple LLVM version 6.0 (clang-600.0.54) (based on LLVM 3.5svn)
>
>>
>
>>
>
>> Marc
>
>>
>
>>
>
>>
>
>> __________________________
>
>>
>
>> Marc Chamberland, MSc
>
>> PhD candidate
>
>> Department of Physics
>
>> Carleton University
>
>> Ottawa (ON)
>
>>
>
>>> On Apr 9, 2015, at 17:17, Mojca Miklavec
>>> <mojca.miklavec.lists-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org>
>
>>> wrote:
>
>>>
>
>>> Hi,
>
>>>
>
>>> Today I updated Geant4 from version 4.10.01 to 4.10.01.p01 and Gate no
>
>>> longer compiles for me. Below is the first error (followed by many
>
>>> others). I was surprised because Geant4 developers aren't really known
>
>>> for breaking backward compatibility with patch releases.
>
>>>
>
>>> I'm using clang 3.4 which is not fully C++11 compatible (at least not
>
>>> without switching to a different stdlib) in case that this was the
>
>>> reason for the failure.
>
>>>
>
>>> I will try to investigate where the problem could be, but I wanted to
>
>>> report an issue. Maybe others will be able to reproduce it and
>
>>> diagnose it faster than I can. (There is a chance that it's not
>
>>> related to the code in Gate and that Geant4 is to blame.)
>
>>>
>
>>> Mojca
>
>>>
>
>>>
>
>>>
>
>>> [ 87%] Building CXX object
>
>>> CMakeFiles/Gate.dir/source/physics/src/GateVSourceMessenger.cc.o
>
>>> /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/clang++
>
>>>  -DG4INTY_USE_QT -DG4UI_USE -DG4UI_USE_QT -DG4UI_USE_TCSH -DG4VERBOSE
>
>>> -DG4VIS_USE -DG4VIS_USE_OPENGL -DG4VIS_USE_OPENGLQT
>
>>> -DG4_STORE_TRAJECTORY -pipe -Os -arch x86_64 -stdlib=libstdc++ -W
>
>>> -Wall -pedantic -Wno-non-virtual-dtor -Wno-long-long -Wwrite-strings
>
>>> -Wpointer-arith -Woverloaded-virtual -Wno-variadic-macros -Wshadow
>
>>> -pipe -Qunused-arguments -std=c++98 -pipe -Os -arch x86_64
>
>>> -stdlib=libstdc++  -Wno-shadow -DNDEBUG -DNDEBUG -arch x86_64
>
>>> -mmacosx-version-min=10.7 -F//System/Library/Frameworks
>
>>> -I/opt/local/libexec/root5/include/root -I/path/to/build -isystem
>
>>> /opt/local/include/Geant4/Geant4.10.1/Geant4 -isystem
>
>>> /opt/local/include -isystem /opt/local/include/QtCore -isystem
>
>>> /opt/local/include/QtGui -isystem /opt/local/include/QtOpenGL
>
>>> -I/path/to/Gate-a20d5f0/source/arf/include
>
>>> -I/path/to/Gate-a20d5f0/source/geometry/include
>
>>> -I/path/to/Gate-a20d5f0/source/physics/include
>
>>> -I/path/to/Gate-a20d5f0/source/digits_hits/include
>
>>> -I/path/to/Gate-a20d5f0/source/general/include
>
>>> -I/path/to/Gate-a20d5f0/source/externals/clhep/include
>
>>> -I/path/to/Gate-a20d5f0/source/externals/itk-mhd
>
>>> -I/path/to/Gate-a20d5f0/source/externals/itk-mhd/itkzlib
>
>>> -I/path/to/Gate-a20d5f0/source/externals/itk-mhd/MetaIO
>
>>> -I/path/to/build/itk-mhd    -o
>
>>> CMakeFiles/Gate.dir/source/physics/src/GateVSourceMessenger.cc.o -c
>
>>> /path/to/Gate-a20d5f0/source/physics/src/GateVSourceMessenger.cc
>
>>> In file included from
>
>>> /path/to/Gate-a20d5f0/source/physics/src/GateVSource.cc:16:
>
>>> In file included from
>
>>> /path/to/Gate-a20d5f0/source/physics/include/GateBackToBack.hh:14:
>
>>> In file included from
>
>>> /path/to/Gate-a20d5f0/source/physics/include/GateVSource.hh:16:
>
>>> In file included from
>
>>> /opt/local/include/Geant4/Geant4.10.1/Geant4/G4SingleParticleSource.hh:124:
>
>>> In file included from
>
>>> /opt/local/include/Geant4/Geant4.10.1/Geant4/G4SPSPosDistribution.hh:141:
>
>>> /opt/local/include/Geant4/Geant4.10.1/Geant4/G4SPSRandomGenerator.hh:156:7:
>
>>> error: field of type 'G4Cache<G4SPSRandomGenerator::a_check>' has
>
>>> private copy constructor
>
>>> class G4SPSRandomGenerator {
>
>>>      ^
>
>>> _______________________________________________
>
>>> Gate-users mailing list
>
>>> Gate-users-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhSoIsB4E12gc@public.gmane.org
>
>>> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>
>>
>
>> _______________________________________________
>
>> Gate-users mailing list
>
>> Gate-users-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhSoIsB4E12gc@public.gmane.org
>
>> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>
> _______________________________________________
>
> Gate-users mailing list
>
> Gate-users-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhSoIsB4E12gc@public.gmane.org
>
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>


------------------------------

Message: 3
Date: Tue, 14 Apr 2015 15:21:17 +0200
From: Marc Verderi <verderi-wBnxOHnizbs@public.gmane.org>
To: gate-users-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhSoIsB4E12gc@public.gmane.org
Subject: Re: [Gate-users] Gate 7.1 build broken after updating to
    Geant4.10.01.p01
Message-ID: <552D144D.5050808-wBnxOHnizbs@public.gmane.org>
Content-Type: text/plain; charset=utf-8; format=flowed

Dear all,

Andrea Dotti <adotti <at> slac.stanford.edu>, from Geant4, has been informed
of the problem. Could involved Gate developers contact him for providing
a solution ?

Cheers,
Marc


On 04/14/2015 03:16 PM, Alex Vergara Gil wrote:
> Indeed! In the G4 page there is no link to the source code of 10.1
> version without the patch01.
>
> This is dangerous for Gate, shall it be created a mirror image of the
> current working G4 source code in the
> http://www.opengatecollaboration.org page??? I think this would be a
> temporary solution at least for now.
>
> Regards
> Alex
>
>
> 2015-04-14 9:06 GMT-04:00, REUZE SYLVAIN <SYLVAIN.REUZE <at> gustaveroussy.fr>:
>> Hello,
>>
>>
>>
>> I have the same compilation errors with Geant4 10.01 p01. I'd like to
>> install Gate 7.1 in order to solve another problem, but I can't find the
>> link to download Geant4 10.01, without the patch.
>>
>> Could you help me ?
>>
>>
>>
>> Best regards,
>>
>> Sylvain Reuz?
>> Stagiaire physicien
>> D?partement de Radioth?rapie
>> Service de physique m?dicale
>>
>> T?l. : 01 42 11 54 06
>> www.gustaveroussy.fr<http://www.gustaveroussy.fr>
>> [cid:image001.gif-WvmHr/fbjPMkNSLZUKRfcjRmByFHzeGd@public.gmane.org]
>>
>>
>>
>>
>>
>>
>>
>> -----Message d'origine-----
>> De : gate-users-bounces-ZwoEplunGu0bJ6lf+0Zptg@public.gmane.orgcollaboration.org
>> [mailto:gate-users-bounces-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEh@public.gmane.orgion.org] De la part de
>> Alex Vergara Gil
>> Envoy? : lundi 13 avril 2015 13:12
>> ? : Marc Chamberland
>> Cc : Gate Users
>> Objet : Re: [Gate-users] Gate 7.1 build broken after updating to
>> Geant4.10.01.p01
>>
>>
>>
>> For what it's worth, I can reproduce the compilation errors, too.
>>
>> There might be a compatibility issue with the
>>
>> G4SPSRandomGenerator.hh:156:7:
>>
>> error: field of type 'G4Cache<G4SPSRandomGenerator::a_check>' has
>>
>> private copy constructor
>>
>> class G4SPSRandomGenerator {
>>
>>
>>
>> I'm compiling with:
>>
>>
>>
>> GCC version 4.8
>>
>>
>>
>> Alex
>>
>>
>>
>> 2015-04-10 11:53 GMT-04:00, Marc Chamberland
>>
>> <MarcChamberland-binDkDF82+79HEwTC5yFM/d9D2ou9A/h@public.gmane.org>:
>>
>>> For what it's worth, I can reproduce the compilation errors, too. I'm
>>> compiling with:
>>> Apple LLVM version 6.0 (clang-600.0.54) (based on LLVM 3.5svn)
>>> Marc
>>> __________________________
>>> Marc Chamberland, MSc
>>> PhD candidate
>>> Department of Physics
>>> Carleton University
>>> Ottawa (ON)
>>>> On Apr 9, 2015, at 17:17, Mojca Miklavec
>>>> <mojca.miklavec.lists-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org>
>>>> wrote:
>>>> Hi,
>>>> Today I updated Geant4 from version 4.10.01 to 4.10.01.p01 and Gate no
>>>> longer compiles for me. Below is the first error (followed by many
>>>> others). I was surprised because Geant4 developers aren't really known
>>>> for breaking backward compatibility with patch releases.
>>>> I'm using clang 3.4 which is not fully C++11 compatible (at least not
>>>> without switching to a different stdlib) in case that this was the
>>>> reason for the failure.
>>>> I will try to investigate where the problem could be, but I wanted to
>>>> report an issue. Maybe others will be able to reproduce it and
>>>> diagnose it faster than I can. (There is a chance that it's not
>>>> related to the code in Gate and that Geant4 is to blame.)
>>>> Mojca
>>>> [ 87%] Building CXX object
>>>> CMakeFiles/Gate.dir/source/physics/src/GateVSourceMessenger.cc.o
>>>> /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/clang++
>>>>  -DG4INTY_USE_QT -DG4UI_USE -DG4UI_USE_QT -DG4UI_USE_TCSH -DG4VERBOSE
>>>> -DG4VIS_USE -DG4VIS_USE_OPENGL -DG4VIS_USE_OPENGLQT
>>>> -DG4_STORE_TRAJECTORY -pipe -Os -arch x86_64 -stdlib=libstdc++ -W
>>>> -Wall -pedantic -Wno-non-virtual-dtor -Wno-long-long -Wwrite-strings
>>>> -Wpointer-arith -Woverloaded-virtual -Wno-variadic-macros -Wshadow
>>>> -pipe -Qunused-arguments -std=c++98 -pipe -Os -arch x86_64
>>>> -stdlib=libstdc++  -Wno-shadow -DNDEBUG -DNDEBUG -arch x86_64
>>>> -mmacosx-version-min=10.7 -F//System/Library/Frameworks
>>>> -I/opt/local/libexec/root5/include/root -I/path/to/build -isystem
>>>> /opt/local/include/Geant4/Geant4.10.1/Geant4 -isystem
>>>> /opt/local/include -isystem /opt/local/include/QtCore -isystem
>>>> /opt/local/include/QtGui -isystem /opt/local/include/QtOpenGL
>>>> -I/path/to/Gate-a20d5f0/source/arf/include
>>>> -I/path/to/Gate-a20d5f0/source/geometry/include
>>>> -I/path/to/Gate-a20d5f0/source/physics/include
>>>> -I/path/to/Gate-a20d5f0/source/digits_hits/include
>>>> -I/path/to/Gate-a20d5f0/source/general/include
>>>> -I/path/to/Gate-a20d5f0/source/externals/clhep/include
>>>> -I/path/to/Gate-a20d5f0/source/externals/itk-mhd
>>>> -I/path/to/Gate-a20d5f0/source/externals/itk-mhd/itkzlib
>>>> -I/path/to/Gate-a20d5f0/source/externals/itk-mhd/MetaIO
>>>> -I/path/to/build/itk-mhd    -o
>>>> CMakeFiles/Gate.dir/source/physics/src/GateVSourceMessenger.cc.o -c
>>>> /path/to/Gate-a20d5f0/source/physics/src/GateVSourceMessenger.cc
>>>> In file included from
>>>> /path/to/Gate-a20d5f0/source/physics/src/GateVSource.cc:16:
>>>> In file included from
>>>> /path/to/Gate-a20d5f0/source/physics/include/GateBackToBack.hh:14:
>>>> In file included from
>>>> /path/to/Gate-a20d5f0/source/physics/include/GateVSource.hh:16:
>>>> In file included from
>>>> /opt/local/include/Geant4/Geant4.10.1/Geant4/G4SingleParticleSource.hh:124:
>>>> In file included from
>>>> /opt/local/include/Geant4/Geant4.10.1/Geant4/G4SPSPosDistribution.hh:141:
>>>> /opt/local/include/Geant4/Geant4.10.1/Geant4/G4SPSRandomGenerator.hh:156:7:
>>>> error: field of type 'G4Cache<G4SPSRandomGenerator::a_check>' has
>>>> private copy constructor
>>>> class G4SPSRandomGenerator {
>>>>      ^
>>>> _______________________________________________
>>>> Gate-users mailing list
>>>> Gate-users-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhSoIsB4E12gc@public.gmane.org
>>>> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>>> _______________________________________________
>>> Gate-users mailing list
>>> Gate-users-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEh9Sl0XMgJPXI@public.gmane.orgon.org
>>> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>> _______________________________________________
>>
>> Gate-users mailing list
>>
>> Gate-users-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhSoIsB4E12gc@public.gmane.org
>>
>> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>>
> _______________________________________________
> Gate-users mailing list
> Gate-users-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhSoIsB4E12gc@public.gmane.org
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users



------------------------------

Message: 4
Date: Tue, 14 Apr 2015 16:47:15 +0200
From: imen elhamrouni <imenelhamrouni-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org>
To: gate-users-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhSoIsB4E12gc@public.gmane.org
Subject: [Gate-users] Geant4 and Eclipse
Message-ID:
    <CAMaL-We2FxJYTeK_jYhpLu3F9cw9Qm5HQA_sGfohhLs-yscG5g-JsoAwUIsXosN+BqQ9rBEUg@public.gmane.org>
Content-Type: text/plain; charset="utf-8"

Hello all,

Is There any way to run an example of Geant4 ,B1 for example, from a new
project in  eclipse , which are the libraries i should add to this
project? Is there any specific configuration for eclipse?What should i do
?..

Thanks!
Imen
?
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20150414/699be61a/attachment-0001.html>

------------------------------

Message: 5
Date: Tue, 14 Apr 2015 18:19:09 +0200 (CEST)
From: Mathieu DUPONT <mdupont-/G0kxDag6VrIdshrf1d+Pg@public.gmane.org>
To: BERGAMASCHI Antoine <antoine.bergamaschi-eE6tJkUsrWpDGTAf/f00U7HwLv5LXcrC@public.gmane.org>
Cc: gate-users-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEh4KqKFgQWhpY@public.gmane.orgorg
Subject: Re: [Gate-users] Root not working in VGate 3.0
Message-ID:
    <1799207786.5910617.1429028349637.JavaMail.zimbra-/G0kxDag6VrIdshrf1d+Pg@public.gmane.org>
Content-Type: text/plain; charset=utf-8

Hi Antoine,

you are right. It's not working.

For the libjpeg missing, simply do :

sudo apt-get install libjpeg62

And after that, you will miss some ?libtiff stuff?.

To fix this, I did

sudo apt-get install libtiff4
sudo ln -s /usr/lib/x86_64-linux-gnu/libtiff.so.4 /usr/lib/libtiff.so.3

And it's working now (at least, the Draw() command in root). It's not a clean solution, but I hope it will help you.


Mathieu Dupont



----- Original Message -----
From: "BERGAMASCHI Antoine" <antoine.bergamaschi-eE6tJkUsrWpDGTAf/f00U7HwLv5LXcrC@public.gmane.org>
To: gate-users <at> lists.opengatecollaboration.org
Sent: Tuesday, April 14, 2015 10:05:51 AM
Subject: [Gate-users] Root not working in VGate 3.0

Hi all,

Everything is in the title, i am using the Gate V delivered as the version 3.0 and the root programme crash when using the "draw()" method due to the lack or a corrupted "libjpeg.so.62" and the fail to load libASImage.so ( which is present in the directory ).

I did not change a single bit of configuration on the VM.

Is a working version of the very usefull vGate disponible anywhere to download ?

Best regards,

Antoine Bergamaschi
_______________________________________________
Gate-users mailing list
Gate-users-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhSoIsB4E12gc@public.gmane.org
http://lists.opengatecollaboration.org/mailman/listinfo/gate-users


------------------------------

_______________________________________________
Gate-users mailing list
Gate-users-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhSoIsB4E12gc@public.gmane.org
http://lists.opengatecollaboration.org/mailman/listinfo/gate-users

End of Gate-users Digest, Vol 107, Issue 32
*******************************************


<div><div>
<div dir="ltr"><span>Please don't send me your mails in future</span></div>
<div dir="ltr"><br></div>  <br><div class="qtdSeparateBR">
<br><br>
</div>
<div class="yahoo_quoted"> <div> <div> <div dir="ltr">  Le Mardi 14 avril 2015 17h19, "gate-users-request@...g" &lt;gate-users-request@...&gt; a &eacute;crit :<br> </div>  <br><br><div class="y_msg_container">Send Gate-users mailing list submissions to<br>&nbsp;&nbsp;&nbsp; <a ymailto="mailto:gate-users@..." href="mailto:gate-users <at> lists.opengatecollaboration.org">gate-users@...</a><br><br>To subscribe or unsubscribe via the World Wide Web, visit<br>&nbsp;&nbsp;&nbsp; <a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><br>or, via email, send a message with subject or body 'help' to<br>&nbsp;&nbsp;&nbsp; <a ymailto="mailto:gate-users-request@..." href="mailto:gate-users-request <at> lists.opengatecollaboration.org">gate-users-request@...n.org</a><br><br>You can reach the person managing the list at<br>&nbsp;&nbsp;&nbsp; <a ymailto="mailto:gate-users-owner@....org" href="mailto:gate-users-owner@...">gate-users-owner@...</a><br><br>When replying, please edit your Subject line so it is more specific<br>than "Re: Contents of Gate-users digest..."<br><br><br>Today's Topics:<br><br>&nbsp;  1. Re: ion source timestamp bug, Gate 7.1 (Gergely Patay)<br>&nbsp;  2. Re: Gate 7.1 build broken after updating to&nbsp;&nbsp;&nbsp; Geant4.10.01.p01<br>&nbsp; &nbsp; &nbsp; (Alex Vergara Gil)<br>&nbsp;  3. Re: Gate 7.1 build broken after updating to&nbsp;&nbsp;&nbsp; Geant4.10.01.p01<br>&nbsp; &nbsp; &nbsp; (Marc Verderi)<br>&nbsp;  4. Geant4 and Eclipse (imen elhamrouni)<br>&nbsp;  5. Re: Root not working in VGate 3.0 (Mathieu DUPONT)<br><br><br>----------------------------------------------------------------------<br><br>Message: 1<br>Date: Tue, 14 Apr 2015 15:13:32 +0200<br>From: Gergely Patay &lt;<a ymailto="mailto:gergely.patay@..." href="mailto:gergely.patay <at> mediso.com">gergely.patay@...</a>&gt;<br>To: Gate User &lt;<a ymailto="mailto:gate-users@..." href="mailto:gate-users@...">gate-users@...on.org</a>&gt;<br>Subject: Re: [Gate-users] ion source timestamp bug, Gate 7.1<br>Message-ID: &lt;<a ymailto="mailto:552D127C.7040800@..." href="mailto:552D127C.7040800@...">552D127C.7040800@...</a>&gt;<br>Content-Type: text/plain; charset="utf-8"; format=flowed<br><br>Dear Gaters,<br><br>to answer my own question: the problem is indeed deep inside Geant4, and as <br>usual, originally this is "not a bug, but a feature".<br><br>In GateVSource.cc:268 we (try to) force life time of the source particle to zero:<br>// Force life time to 0, time is managed by GATE not G4<br>GetParticleDefinition()-&gt;SetPDGLifeTime(0);<br><br>However, from v10.0+ when Geant4 reads back the life time of a general ion <br>source, it does *not* retrieve the previously set value, but an (otherwise <br>physically correct) value from an external data file <br>(G4ENSDFSTATE1.0/ENSDFSTATE.dat) :<br><br>G4ParticleDefinition.icc:45<br>inline<br>&nbsp; G4double G4ParticleDefinition::GetPDGLifeTime() const<br>{<br>&nbsp;  if(IsGeneralIon())<br>&nbsp;  { return GetIonLifeTime(); }<br>&nbsp;  else<br>&nbsp;  { return thePDGLifeTime; }<br>}<br><br>Unfortunately, the underlying isGeneralIon variable is not directly accessible <br>from outside (i.e. from Gate).<br><br>A minimal workaround might be as follows:<br><br>1) change visibility of isGeneralIon from protected to public:<br>geant4/particles/management/include/G4ParticleDefinition.hh:402<br>-&nbsp;  protected:<br>+&nbsp;  public:<br>&nbsp; &nbsp; &nbsp; &nbsp; G4bool isGeneralIon;<br><br>(this might be modified *without* recompiling Geant4)<br><br>2) add the following line to gate/source/physics/src/GateVSource.cc:269<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  GetParticleDefinition()-&gt;SetPDGLifeTime(0);<br>+&nbsp; &nbsp; &nbsp; &nbsp;  GetParticleDefinition()-&gt;isGeneralIon=false;<br><br>3) recompile Gate with the modified header file.<br><br>Since this modification affects Geant4 code too, I cannot provide a simple patch <br>file yet.<br><br>best regards,<br>&nbsp; Gergely<br><br>&nbsp; --<br>Gergely Patay<br>Mediso Ltd.<br><br><br><br>On 2015-04-01 13:38, Gergely Patay wrote:<br>&gt; Dear Gaters,<br>&gt;<br>&gt; I installed the latest version of Gate from git. Everything went fine<br>&gt; (benchmarks etc), until I tried to use the 'generic ion' source for a PET<br>&gt; simulation.<br>&gt;<br>&gt; The timestamps of the hits (and consequently, of the singles and coincidences)<br>&gt; are weird. They should be distributed uniformly in the [0,1] s interval (as seen<br>&gt; on the ion_timestamp_gate7.0.png figure), but in 7.1 they are distributed<br>&gt; exponentially over a much-much larger time interval with an exponent which<br>&gt; equals to the decay rate of the ion, in this case Na-22.<br>&gt;<br>&gt; The macro files were identical, and the corresponding software versions were<br>&gt; - Gate7.0, Geant4.9.6.p03, root5.34.18<br>&gt; - Gate7.1, Geant4.10.01, root5.34.26<br>&gt;<br>&gt; This bug is actually the very same to that was reported on March 31, 2011 for<br>&gt; Gate 6.0.1<br>&gt; (<a href="http://lists.opengatecollaboration.org/mailman/private/gate-users/2011-March/004895.html" target="_blank">http://lists.opengatecollaboration.org/mailman/private/gate-users/2011-March/004895.html</a>)<br>&gt; which somehow disappeared from 7.0 but re-emerged in 7.1.<br>&gt;<br>&gt; The bug does not affect the backtoback or e+ sources, and energy spectra are fine.<br>&gt; It might be buried somewhere deep in Geant4.<br>&gt;<br>&gt; Do you have any idea what went wrong?<br>&gt;<br>&gt; thanks,<br>&gt;&nbsp;  Gergely<br>&gt;<br>&gt; --<br>&gt; Gergely Patay<br>&gt; Mediso Ltd.<br>&gt;<br>&gt;<br>&gt; _______________________________________________<br>&gt; Gate-users mailing list<br>&gt; <a ymailto="mailto:Gate-users@..." href="mailto:Gate-users@...ngatecollaboration.org">Gate-users@...</a><br>&gt; <a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><br>&gt;<br><br><br>------------------------------<br><br>Message: 2<br>Date: Tue, 14 Apr 2015 09:16:01 -0400<br>From: Alex Vergara Gil &lt;<a ymailto="mailto:alexvergaragil@..." href="mailto:alexvergaragil@...">alexvergaragil@...</a>&gt;<br>To: REUZE SYLVAIN &lt;<a ymailto="mailto:SYLVAIN.REUZE@..." href="mailto:SYLVAIN.REUZE@...">SYLVAIN.REUZE@...</a>&gt;<br>Cc: Gate Users &lt;<a ymailto="mailto:gate-users@...on.org" href="mailto:gate-users@...">gate-users@...</a>&gt;<br>Subject: Re: [Gate-users] Gate 7.1 build broken after updating to<br>&nbsp;&nbsp;&nbsp; Geant4.10.01.p01<br>Message-ID:<br>&nbsp;&nbsp;&nbsp; &lt;CAMF2x=nDWVfRCzt5sn0RuhkwBpEme594sZpB70+<a ymailto="mailto:Z_N5pbjtO4Q@..." href="mailto:Z_N5pbjtO4Q@...">Z_N5pbjtO4Q@...</a>&gt;<br>Content-Type: text/plain; charset=UTF-8<br><br>Indeed! In the G4 page there is no link to the source code of 10.1<br>version without the patch01.<br><br>This is dangerous for Gate, shall it be created a mirror image of the<br>current working G4 source code in the<br><a href="http://www.opengatecollaboration.org/" target="_blank">http://www.opengatecollaboration.org </a>page??? I think this would be a<br>temporary solution at least for now.<br><br>Regards<br>Alex<br><br><br>2015-04-14 9:06 GMT-04:00, REUZE SYLVAIN &lt;<a ymailto="mailto:SYLVAIN.REUZE@..." href="mailto:SYLVAIN.REUZE <at> gustaveroussy.fr">SYLVAIN.REUZE@...</a>&gt;:<br>&gt; Hello,<br>&gt;<br>&gt;<br>&gt;<br>&gt; I have the same compilation errors with Geant4 10.01 p01. I'd like to<br>&gt; install Gate 7.1 in order to solve another problem, but I can't find the<br>&gt; link to download Geant4 10.01, without the patch.<br>&gt;<br>&gt; Could you help me ?<br>&gt;<br>&gt;<br>&gt;<br>&gt; Best regards,<br>&gt;<br>&gt; Sylvain Reuz?<br>&gt; Stagiaire physicien<br>&gt; D?partement de Radioth?rapie<br>&gt; Service de physique m?dicale<br>&gt;<br>&gt; T?l. : 01 42 11 54 06<br>&gt; www.gustaveroussy.fr&lt;<a href="http://www.gustaveroussy.fr/" target="_blank">http://www.gustaveroussy.fr</a>&gt;<br>&gt; [cid:<a ymailto="mailto:image001.gif <at> 01D076C4.8D1D9C60" href="mailto:image001.gif@...">image001.gif <at> 01D076C4.8D1D9C60</a>]<br>&gt;<br>&gt;<br>&gt;<br>&gt;<br>&gt;<br>&gt;<br>&gt;<br>&gt; -----Message d'origine-----<br>&gt; De : <a ymailto="mailto:gate-users-bounces@..." href="mailto:gate-users-bounces@...">gate-users-bounces@...boration.org</a><br>&gt; [mailto:<a ymailto="mailto:gate-users-bounces <at> lists.opengatecollaboration.org" href="mailto:gate-users-bounces@...gatecollaboration.org">gate-users-bounces@...</a>] De la part de<br>&gt; Alex Vergara Gil<br>&gt; Envoy? : lundi 13 avril 2015 13:12<br>&gt; ? : Marc Chamberland<br>&gt; Cc : Gate Users<br>&gt; Objet : Re: [Gate-users] Gate 7.1 build broken after updating to<br>&gt; Geant4.10.01.p01<br>&gt;<br>&gt;<br>&gt;<br>&gt; For what it's worth, I can reproduce the compilation errors, too.<br>&gt;<br>&gt; There might be a compatibility issue with the<br>&gt;<br>&gt; G4SPSRandomGenerator.hh:156:7:<br>&gt;<br>&gt; error: field of type 'G4Cache&lt;G4SPSRandomGenerator::a_check&gt;' has<br>&gt;<br>&gt; private copy constructor<br>&gt;<br>&gt; class G4SPSRandomGenerator {<br>&gt;<br>&gt;<br>&gt;<br>&gt; I'm compiling with:<br>&gt;<br>&gt;<br>&gt;<br>&gt; GCC version 4.8<br>&gt;<br>&gt;<br>&gt;<br>&gt; Alex<br>&gt;<br>&gt;<br>&gt;<br>&gt; 2015-04-10 11:53 GMT-04:00, Marc Chamberland<br>&gt;<br>&gt; &lt;<a ymailto="mailto:MarcChamberland@...ton.ca" href="mailto:MarcChamberland@...">MarcChamberland <at> cmail.carleton.ca</a>&gt;:<br>&gt;<br>&gt;&gt; For what it's worth, I can reproduce the compilation errors, too. I'm<br>&gt;<br>&gt;&gt; compiling with:<br>&gt;<br>&gt;&gt;<br>&gt;<br>&gt;&gt; Apple LLVM version 6.0 (clang-600.0.54) (based on LLVM 3.5svn)<br>&gt;<br>&gt;&gt;<br>&gt;<br>&gt;&gt;<br>&gt;<br>&gt;&gt; Marc<br>&gt;<br>&gt;&gt;<br>&gt;<br>&gt;&gt;<br>&gt;<br>&gt;&gt;<br>&gt;<br>&gt;&gt; __________________________<br>&gt;<br>&gt;&gt;<br>&gt;<br>&gt;&gt; Marc Chamberland, MSc<br>&gt;<br>&gt;&gt; PhD candidate<br>&gt;<br>&gt;&gt; Department of Physics<br>&gt;<br>&gt;&gt; Carleton University<br>&gt;<br>&gt;&gt; Ottawa (ON)<br>&gt;<br>&gt;&gt;<br>&gt;<br>&gt;&gt;&gt; On Apr 9, 2015, at 17:17, Mojca Miklavec<br>&gt;&gt;&gt; &lt;<a ymailto="mailto:mojca.miklavec.lists@..." href="mailto:mojca.miklavec.lists@...">mojca.miklavec.lists@...</a>&gt;<br>&gt;<br>&gt;&gt;&gt; wrote:<br>&gt;<br>&gt;&gt;&gt;<br>&gt;<br>&gt;&gt;&gt; Hi,<br>&gt;<br>&gt;&gt;&gt;<br>&gt;<br>&gt;&gt;&gt; Today I updated Geant4 from version 4.10.01 to 4.10.01.p01 and Gate no<br>&gt;<br>&gt;&gt;&gt; longer compiles for me. Below is the first error (followed by many<br>&gt;<br>&gt;&gt;&gt; others). I was surprised because Geant4 developers aren't really known<br>&gt;<br>&gt;&gt;&gt; for breaking backward compatibility with patch releases.<br>&gt;<br>&gt;&gt;&gt;<br>&gt;<br>&gt;&gt;&gt; I'm using clang 3.4 which is not fully C++11 compatible (at least not<br>&gt;<br>&gt;&gt;&gt; without switching to a different stdlib) in case that this was the<br>&gt;<br>&gt;&gt;&gt; reason for the failure.<br>&gt;<br>&gt;&gt;&gt;<br>&gt;<br>&gt;&gt;&gt; I will try to investigate where the problem could be, but I wanted to<br>&gt;<br>&gt;&gt;&gt; report an issue. Maybe others will be able to reproduce it and<br>&gt;<br>&gt;&gt;&gt; diagnose it faster than I can. (There is a chance that it's not<br>&gt;<br>&gt;&gt;&gt; related to the code in Gate and that Geant4 is to blame.)<br>&gt;<br>&gt;&gt;&gt;<br>&gt;<br>&gt;&gt;&gt; Mojca<br>&gt;<br>&gt;&gt;&gt;<br>&gt;<br>&gt;&gt;&gt;<br>&gt;<br>&gt;&gt;&gt;<br>&gt;<br>&gt;&gt;&gt; [ 87%] Building CXX object<br>&gt;<br>&gt;&gt;&gt; CMakeFiles/Gate.dir/source/physics/src/GateVSourceMessenger.cc.o<br>&gt;<br>&gt;&gt;&gt; /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/clang++<br>&gt;<br>&gt;&gt;&gt;&nbsp; -DG4INTY_USE_QT -DG4UI_USE -DG4UI_USE_QT -DG4UI_USE_TCSH -DG4VERBOSE<br>&gt;<br>&gt;&gt;&gt; -DG4VIS_USE -DG4VIS_USE_OPENGL -DG4VIS_USE_OPENGLQT<br>&gt;<br>&gt;&gt;&gt; -DG4_STORE_TRAJECTORY -pipe -Os -arch x86_64 -stdlib=libstdc++ -W<br>&gt;<br>&gt;&gt;&gt; -Wall -pedantic -Wno-non-virtual-dtor -Wno-long-long -Wwrite-strings<br>&gt;<br>&gt;&gt;&gt; -Wpointer-arith -Woverloaded-virtual -Wno-variadic-macros -Wshadow<br>&gt;<br>&gt;&gt;&gt; -pipe -Qunused-arguments -std=c++98 -pipe -Os -arch x86_64<br>&gt;<br>&gt;&gt;&gt; -stdlib=libstdc++&nbsp; -Wno-shadow -DNDEBUG -DNDEBUG -arch x86_64<br>&gt;<br>&gt;&gt;&gt; -mmacosx-version-min=10.7 -F//System/Library/Frameworks<br>&gt;<br>&gt;&gt;&gt; -I/opt/local/libexec/root5/include/root -I/path/to/build -isystem<br>&gt;<br>&gt;&gt;&gt; /opt/local/include/Geant4/Geant4.10.1/Geant4 -isystem<br>&gt;<br>&gt;&gt;&gt; /opt/local/include -isystem /opt/local/include/QtCore -isystem<br>&gt;<br>&gt;&gt;&gt; /opt/local/include/QtGui -isystem /opt/local/include/QtOpenGL<br>&gt;<br>&gt;&gt;&gt; -I/path/to/Gate-a20d5f0/source/arf/include<br>&gt;<br>&gt;&gt;&gt; -I/path/to/Gate-a20d5f0/source/geometry/include<br>&gt;<br>&gt;&gt;&gt; -I/path/to/Gate-a20d5f0/source/physics/include<br>&gt;<br>&gt;&gt;&gt; -I/path/to/Gate-a20d5f0/source/digits_hits/include<br>&gt;<br>&gt;&gt;&gt; -I/path/to/Gate-a20d5f0/source/general/include<br>&gt;<br>&gt;&gt;&gt; -I/path/to/Gate-a20d5f0/source/externals/clhep/include<br>&gt;<br>&gt;&gt;&gt; -I/path/to/Gate-a20d5f0/source/externals/itk-mhd<br>&gt;<br>&gt;&gt;&gt; -I/path/to/Gate-a20d5f0/source/externals/itk-mhd/itkzlib<br>&gt;<br>&gt;&gt;&gt; -I/path/to/Gate-a20d5f0/source/externals/itk-mhd/MetaIO<br>&gt;<br>&gt;&gt;&gt; -I/path/to/build/itk-mhd&nbsp; &nbsp; -o<br>&gt;<br>&gt;&gt;&gt; CMakeFiles/Gate.dir/source/physics/src/GateVSourceMessenger.cc.o -c<br>&gt;<br>&gt;&gt;&gt; /path/to/Gate-a20d5f0/source/physics/src/GateVSourceMessenger.cc<br>&gt;<br>&gt;&gt;&gt; In file included from<br>&gt;<br>&gt;&gt;&gt; /path/to/Gate-a20d5f0/source/physics/src/GateVSource.cc:16:<br>&gt;<br>&gt;&gt;&gt; In file included from<br>&gt;<br>&gt;&gt;&gt; /path/to/Gate-a20d5f0/source/physics/include/GateBackToBack.hh:14:<br>&gt;<br>&gt;&gt;&gt; In file included from<br>&gt;<br>&gt;&gt;&gt; /path/to/Gate-a20d5f0/source/physics/include/GateVSource.hh:16:<br>&gt;<br>&gt;&gt;&gt; In file included from<br>&gt;<br>&gt;&gt;&gt; /opt/local/include/Geant4/Geant4.10.1/Geant4/G4SingleParticleSource.hh:124:<br>&gt;<br>&gt;&gt;&gt; In file included from<br>&gt;<br>&gt;&gt;&gt; /opt/local/include/Geant4/Geant4.10.1/Geant4/G4SPSPosDistribution.hh:141:<br>&gt;<br>&gt;&gt;&gt; /opt/local/include/Geant4/Geant4.10.1/Geant4/G4SPSRandomGenerator.hh:156:7:<br>&gt;<br>&gt;&gt;&gt; error: field of type 'G4Cache&lt;G4SPSRandomGenerator::a_check&gt;' has<br>&gt;<br>&gt;&gt;&gt; private copy constructor<br>&gt;<br>&gt;&gt;&gt; class G4SPSRandomGenerator {<br>&gt;<br>&gt;&gt;&gt;&nbsp; &nbsp; &nbsp; ^<br>&gt;<br>&gt;&gt;&gt; _______________________________________________<br>&gt;<br>&gt;&gt;&gt; Gate-users mailing list<br>&gt;<br>&gt;&gt;&gt; <a ymailto="mailto:Gate-users <at> lists.opengatecollaboration.org" href="mailto:Gate-users@...ration.org">Gate-users@...</a><br>&gt;<br>&gt;&gt;&gt; <a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><br>&gt;<br>&gt;&gt;<br>&gt;<br>&gt;&gt; _______________________________________________<br>&gt;<br>&gt;&gt; Gate-users mailing list<br>&gt;<br>&gt;&gt; <a ymailto="mailto:Gate-users@...ollaboration.org" href="mailto:Gate-users@...">Gate-users@...</a><br>&gt;<br>&gt;&gt; <a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><br>&gt;<br>&gt; _______________________________________________<br>&gt;<br>&gt; Gate-users mailing list<br>&gt;<br>&gt; <a ymailto="mailto:Gate-users@..." href="mailto:Gate-users <at> lists.opengatecollaboration.org">Gate-users@...</a><br>&gt;<br>&gt; <a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><br>&gt;<br><br><br>------------------------------<br><br>Message: 3<br>Date: Tue, 14 Apr 2015 15:21:17 +0200<br>From: Marc Verderi &lt;<a ymailto="mailto:verderi@..." href="mailto:verderi@...">verderi@...</a>&gt;<br>To: <a ymailto="mailto:gate-users@..." href="mailto:gate-users@...engatecollaboration.org">gate-users@...</a><br>Subject: Re: [Gate-users] Gate 7.1 build broken after updating to<br>&nbsp;&nbsp;&nbsp; Geant4.10.01.p01<br>Message-ID: &lt;<a ymailto="mailto:552D144D.5050808@..." href="mailto:552D144D.5050808@...">552D144D.5050808@...</a>&gt;<br>Content-Type: text/plain; charset=utf-8; format=flowed<br><br>Dear all,<br><br>Andrea Dotti &lt;<a ymailto="mailto:adotti@..." href="mailto:adotti@...">adotti <at> slac.stanford.edu</a>&gt;, from Geant4, has been informed <br>of the problem. Could involved Gate developers contact him for providing <br>a solution ?<br><br>Cheers,<br>Marc<br><br><br>On 04/14/2015 03:16 PM, Alex Vergara Gil wrote:<br>&gt; Indeed! In the G4 page there is no link to the source code of 10.1<br>&gt; version without the patch01.<br>&gt;<br>&gt; This is dangerous for Gate, shall it be created a mirror image of the<br>&gt; current working G4 source code in the<br>&gt; <a href="http://www.opengatecollaboration.org/" target="_blank">http://www.opengatecollaboration.org </a>page??? I think this would be a<br>&gt; temporary solution at least for now.<br>&gt;<br>&gt; Regards<br>&gt; Alex<br>&gt;<br>&gt;<br>&gt; 2015-04-14 9:06 GMT-04:00, REUZE SYLVAIN &lt;<a ymailto="mailto:SYLVAIN.REUZE@..." href="mailto:SYLVAIN.REUZE@...">SYLVAIN.REUZE <at> gustaveroussy.fr</a>&gt;:<br>&gt;&gt; Hello,<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt; I have the same compilation errors with Geant4 10.01 p01. I'd like to<br>&gt;&gt; install Gate 7.1 in order to solve another problem, but I can't find the<br>&gt;&gt; link to download Geant4 10.01, without the patch.<br>&gt;&gt;<br>&gt;&gt; Could you help me ?<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt; Best regards,<br>&gt;&gt;<br>&gt;&gt; Sylvain Reuz?<br>&gt;&gt; Stagiaire physicien<br>&gt;&gt; D?partement de Radioth?rapie<br>&gt;&gt; Service de physique m?dicale<br>&gt;&gt;<br>&gt;&gt; T?l. : 01 42 11 54 06<br>&gt;&gt; www.gustaveroussy.fr&lt;<a href="http://www.gustaveroussy.fr/" target="_blank">http://www.gustaveroussy.fr</a>&gt;<br>&gt;&gt; [cid:<a ymailto="mailto:image001.gif@..." href="mailto:image001.gif@...">image001.gif@...</a>]<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt; -----Message d'origine-----<br>&gt;&gt; De : <a ymailto="mailto:gate-users-bounces@..." href="mailto:gate-users-bounces@...">gate-users-bounces@...collaboration.org</a><br>&gt;&gt; [mailto:<a ymailto="mailto:gate-users-bounces@..." href="mailto:gate-users-bounces <at> lists.opengatecollaboration.org">gate-users-bounces@...ion.org</a>] De la part de<br>&gt;&gt; Alex Vergara Gil<br>&gt;&gt; Envoy? : lundi 13 avril 2015 13:12<br>&gt;&gt; ? : Marc Chamberland<br>&gt;&gt; Cc : Gate Users<br>&gt;&gt; Objet : Re: [Gate-users] Gate 7.1 build broken after updating to<br>&gt;&gt; Geant4.10.01.p01<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt; For what it's worth, I can reproduce the compilation errors, too.<br>&gt;&gt;<br>&gt;&gt; There might be a compatibility issue with the<br>&gt;&gt;<br>&gt;&gt; G4SPSRandomGenerator.hh:156:7:<br>&gt;&gt;<br>&gt;&gt; error: field of type 'G4Cache&lt;G4SPSRandomGenerator::a_check&gt;' has<br>&gt;&gt;<br>&gt;&gt; private copy constructor<br>&gt;&gt;<br>&gt;&gt; class G4SPSRandomGenerator {<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt; I'm compiling with:<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt; GCC version 4.8<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt; Alex<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt; 2015-04-10 11:53 GMT-04:00, Marc Chamberland<br>&gt;&gt;<br>&gt;&gt; &lt;<a ymailto="mailto:MarcChamberland <at> cmail.carleton.ca" href="mailto:MarcChamberland@...">MarcChamberland@...</a>&gt;:<br>&gt;&gt;<br>&gt;&gt;&gt; For what it's worth, I can reproduce the compilation errors, too. I'm<br>&gt;&gt;&gt; compiling with:<br>&gt;&gt;&gt; Apple LLVM version 6.0 (clang-600.0.54) (based on LLVM 3.5svn)<br>&gt;&gt;&gt; Marc<br>&gt;&gt;&gt; __________________________<br>&gt;&gt;&gt; Marc Chamberland, MSc<br>&gt;&gt;&gt; PhD candidate<br>&gt;&gt;&gt; Department of Physics<br>&gt;&gt;&gt; Carleton University<br>&gt;&gt;&gt; Ottawa (ON)<br>&gt;&gt;&gt;&gt; On Apr 9, 2015, at 17:17, Mojca Miklavec<br>&gt;&gt;&gt;&gt; &lt;<a ymailto="mailto:mojca.miklavec.lists@..." href="mailto:mojca.miklavec.lists@...">mojca.miklavec.lists@...</a>&gt;<br>&gt;&gt;&gt;&gt; wrote:<br>&gt;&gt;&gt;&gt; Hi,<br>&gt;&gt;&gt;&gt; Today I updated Geant4 from version 4.10.01 to 4.10.01.p01 and Gate no<br>&gt;&gt;&gt;&gt; longer compiles for me. Below is the first error (followed by many<br>&gt;&gt;&gt;&gt; others). I was surprised because Geant4 developers aren't really known<br>&gt;&gt;&gt;&gt; for breaking backward compatibility with patch releases.<br>&gt;&gt;&gt;&gt; I'm using clang 3.4 which is not fully C++11 compatible (at least not<br>&gt;&gt;&gt;&gt; without switching to a different stdlib) in case that this was the<br>&gt;&gt;&gt;&gt; reason for the failure.<br>&gt;&gt;&gt;&gt; I will try to investigate where the problem could be, but I wanted to<br>&gt;&gt;&gt;&gt; report an issue. Maybe others will be able to reproduce it and<br>&gt;&gt;&gt;&gt; diagnose it faster than I can. (There is a chance that it's not<br>&gt;&gt;&gt;&gt; related to the code in Gate and that Geant4 is to blame.)<br>&gt;&gt;&gt;&gt; Mojca<br>&gt;&gt;&gt;&gt; [ 87%] Building CXX object<br>&gt;&gt;&gt;&gt; CMakeFiles/Gate.dir/source/physics/src/GateVSourceMessenger.cc.o<br>&gt;&gt;&gt;&gt; /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/clang++<br>&gt;&gt;&gt;&gt;&nbsp;  -DG4INTY_USE_QT -DG4UI_USE -DG4UI_USE_QT -DG4UI_USE_TCSH -DG4VERBOSE<br>&gt;&gt;&gt;&gt; -DG4VIS_USE -DG4VIS_USE_OPENGL -DG4VIS_USE_OPENGLQT<br>&gt;&gt;&gt;&gt; -DG4_STORE_TRAJECTORY -pipe -Os -arch x86_64 -stdlib=libstdc++ -W<br>&gt;&gt;&gt;&gt; -Wall -pedantic -Wno-non-virtual-dtor -Wno-long-long -Wwrite-strings<br>&gt;&gt;&gt;&gt; -Wpointer-arith -Woverloaded-virtual -Wno-variadic-macros -Wshadow<br>&gt;&gt;&gt;&gt; -pipe -Qunused-arguments -std=c++98 -pipe -Os -arch x86_64<br>&gt;&gt;&gt;&gt; -stdlib=libstdc++&nbsp; -Wno-shadow -DNDEBUG -DNDEBUG -arch x86_64<br>&gt;&gt;&gt;&gt; -mmacosx-version-min=10.7 -F//System/Library/Frameworks<br>&gt;&gt;&gt;&gt; -I/opt/local/libexec/root5/include/root -I/path/to/build -isystem<br>&gt;&gt;&gt;&gt; /opt/local/include/Geant4/Geant4.10.1/Geant4 -isystem<br>&gt;&gt;&gt;&gt; /opt/local/include -isystem /opt/local/include/QtCore -isystem<br>&gt;&gt;&gt;&gt; /opt/local/include/QtGui -isystem /opt/local/include/QtOpenGL<br>&gt;&gt;&gt;&gt; -I/path/to/Gate-a20d5f0/source/arf/include<br>&gt;&gt;&gt;&gt; -I/path/to/Gate-a20d5f0/source/geometry/include<br>&gt;&gt;&gt;&gt; -I/path/to/Gate-a20d5f0/source/physics/include<br>&gt;&gt;&gt;&gt; -I/path/to/Gate-a20d5f0/source/digits_hits/include<br>&gt;&gt;&gt;&gt; -I/path/to/Gate-a20d5f0/source/general/include<br>&gt;&gt;&gt;&gt; -I/path/to/Gate-a20d5f0/source/externals/clhep/include<br>&gt;&gt;&gt;&gt; -I/path/to/Gate-a20d5f0/source/externals/itk-mhd<br>&gt;&gt;&gt;&gt; -I/path/to/Gate-a20d5f0/source/externals/itk-mhd/itkzlib<br>&gt;&gt;&gt;&gt; -I/path/to/Gate-a20d5f0/source/externals/itk-mhd/MetaIO<br>&gt;&gt;&gt;&gt; -I/path/to/build/itk-mhd&nbsp; &nbsp; -o<br>&gt;&gt;&gt;&gt; CMakeFiles/Gate.dir/source/physics/src/GateVSourceMessenger.cc.o -c<br>&gt;&gt;&gt;&gt; /path/to/Gate-a20d5f0/source/physics/src/GateVSourceMessenger.cc<br>&gt;&gt;&gt;&gt; In file included from<br>&gt;&gt;&gt;&gt; /path/to/Gate-a20d5f0/source/physics/src/GateVSource.cc:16:<br>&gt;&gt;&gt;&gt; In file included from<br>&gt;&gt;&gt;&gt; /path/to/Gate-a20d5f0/source/physics/include/GateBackToBack.hh:14:<br>&gt;&gt;&gt;&gt; In file included from<br>&gt;&gt;&gt;&gt; /path/to/Gate-a20d5f0/source/physics/include/GateVSource.hh:16:<br>&gt;&gt;&gt;&gt; In file included from<br>&gt;&gt;&gt;&gt; /opt/local/include/Geant4/Geant4.10.1/Geant4/G4SingleParticleSource.hh:124:<br>&gt;&gt;&gt;&gt; In file included from<br>&gt;&gt;&gt;&gt; /opt/local/include/Geant4/Geant4.10.1/Geant4/G4SPSPosDistribution.hh:141:<br>&gt;&gt;&gt;&gt; /opt/local/include/Geant4/Geant4.10.1/Geant4/G4SPSRandomGenerator.hh:156:7:<br>&gt;&gt;&gt;&gt; error: field of type 'G4Cache&lt;G4SPSRandomGenerator::a_check&gt;' has<br>&gt;&gt;&gt;&gt; private copy constructor<br>&gt;&gt;&gt;&gt; class G4SPSRandomGenerator {<br>&gt;&gt;&gt;&gt;&nbsp; &nbsp; &nbsp;  ^<br>&gt;&gt;&gt;&gt; _______________________________________________<br>&gt;&gt;&gt;&gt; Gate-users mailing list<br>&gt;&gt;&gt;&gt; <a ymailto="mailto:Gate-users@...engatecollaboration.org" href="mailto:Gate-users@...ion.org">Gate-users@...</a><br>&gt;&gt;&gt;&gt; <a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><br>&gt;&gt;&gt; _______________________________________________<br>&gt;&gt;&gt; Gate-users mailing list<br>&gt;&gt;&gt; <a ymailto="mailto:Gate-users@..." href="mailto:Gate-users@...">Gate-users@...on.org</a><br>&gt;&gt;&gt; <a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><br>&gt;&gt; _______________________________________________<br>&gt;&gt;<br>&gt;&gt; Gate-users mailing list<br>&gt;&gt;<br>&gt;&gt; <a ymailto="mailto:Gate-users@...llaboration.org" href="mailto:Gate-users@...">Gate-users@...</a><br>&gt;&gt;<br>&gt;&gt; <a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><br>&gt;&gt;<br>&gt; _______________________________________________<br>&gt; Gate-users mailing list<br>&gt; <a ymailto="mailto:Gate-users@..." href="mailto:Gate-users@...atecollaboration.org">Gate-users@...</a><br>&gt; <a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><br><br><br><br>------------------------------<br><br>Message: 4<br>Date: Tue, 14 Apr 2015 16:47:15 +0200<br>From: imen elhamrouni &lt;<a ymailto="mailto:imenelhamrouni@..." href="mailto:imenelhamrouni@...">imenelhamrouni@...</a>&gt;<br>To: <a ymailto="mailto:gate-users@..." href="mailto:gate-users <at> lists.opengatecollaboration.org">gate-users@...</a><br>Subject: [Gate-users] Geant4 and Eclipse<br>Message-ID:<br>&nbsp;&nbsp;&nbsp; &lt;<a ymailto="mailto:CAMaL-We2FxJYTeK_jYhpLu3F9cw9Qm5HQA_sGfohhLs-yscG5g@..." href="mailto:CAMaL-We2FxJYTeK_jYhpLu3F9cw9Qm5HQA_sGfohhLs-yscG5g@...">CAMaL-We2FxJYTeK_jYhpLu3F9cw9Qm5HQA_sGfohhLs-yscG5g@...</a>&gt;<br>Content-Type: text/plain; charset="utf-8"<br><br>Hello all,<br><br>Is There any way to run an example of Geant4 ,B1 for example, from a new<br>project in&nbsp; eclipse , which are the libraries i should add to this<br> project? Is there any specific configuration for eclipse?What should i do<br>?..<br><br>Thanks!<br>Imen<br>?<br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: &lt;<a href="http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20150414/699be61a/attachment-0001.html" target="_blank">http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20150414/699be61a/attachment-0001.html</a>&gt;<br><br>------------------------------<br><br>Message: 5<br>Date: Tue, 14 Apr 2015 18:19:09 +0200 (CEST)<br>From: Mathieu DUPONT &lt;<a ymailto="mailto:mdupont <at> cppm.in2p3.fr" href="mailto:mdupont@...">mdupont@...</a>&gt;<br>To: BERGAMASCHI Antoine &lt;<a ymailto="mailto:antoine.bergamaschi <at> synchrotron-soleil.fr" href="mailto:antoine.bergamaschi <at> synchrotron-soleil.fr">antoine.bergamaschi@...</a>&gt;<br>Cc: <a ymailto="mailto:gate-users@..." href="mailto:gate-users@...">gate-users@...org</a><br>Subject: Re: [Gate-users] Root not working in VGate 3.0<br>Message-ID:<br>&nbsp;&nbsp;&nbsp; &lt;<a ymailto="mailto:1799207786.5910617.1429028349637.JavaMail.zimbra@..." href="mailto:1799207786.5910617.1429028349637.JavaMail.zimbra@...">1799207786.5910617.1429028349637.JavaMail.zimbra@...</a>&gt;<br>Content-Type: text/plain; charset=utf-8<br><br>Hi Antoine, <br><br>you are right. It's not working. <br><br>For the libjpeg missing, simply do : <br><br>sudo apt-get install libjpeg62<br><br>And after that, you will miss some ?libtiff stuff?. <br><br>To fix this, I did <br><br>sudo apt-get install libtiff4<br>sudo ln -s /usr/lib/x86_64-linux-gnu/libtiff.so.4 /usr/lib/libtiff.so.3<br><br>And it's working now (at least, the Draw() command in root). It's not a clean solution, but I hope it will help you. <br><br><br>Mathieu Dupont<br><br><br><br>----- Original Message -----<br>From: "BERGAMASCHI Antoine" &lt;<a ymailto="mailto:antoine.bergamaschi@..." href="mailto:antoine.bergamaschi@...">antoine.bergamaschi@...</a>&gt;<br>To: <a ymailto="mailto:gate-users@...org" href="mailto:gate-users@...">gate-users <at> lists.opengatecollaboration.org</a><br>Sent: Tuesday, April 14, 2015 10:05:51 AM<br>Subject: [Gate-users] Root not working in VGate 3.0<br><br>Hi all,<br><br>Everything is in the title, i am using the Gate V delivered as the version 3.0 and the root programme crash when using the "draw()" method due to the lack or a corrupted "libjpeg.so.62" and the fail to load libASImage.so ( which is present in the directory ).<br><br>I did not change a single bit of configuration on the VM.<br><br>Is a working version of the very usefull vGate disponible anywhere to download ?<br><br>Best regards,<br><br>Antoine Bergamaschi<br>_______________________________________________<br>Gate-users mailing list<br><a ymailto="mailto:Gate-users@...aboration.org" href="mailto:Gate-users@...">Gate-users@...</a><br><a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><br><br><br>------------------------------<br><br>_______________________________________________<br>Gate-users mailing list<br><a ymailto="mailto:Gate-users <at> lists.opengatecollaboration.org" href="mailto:Gate-users@...ollaboration.org">Gate-users@...</a><br><a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><br><br>End of Gate-users Digest, Vol 107, Issue 32<br>*******************************************<br><br><br>
</div>  </div> </div>  </div>
</div></div>
imen elhamrouni | 14 Apr 16:47 2015
Picon

Geant4 and Eclipse

Hello all,

Is There any way to run an example of Geant4 ,B1 for example, from a new project in  eclipse , which are the libraries i should add to this  project? Is there any specific configuration for eclipse?What should i do ?..

Thanks!
Imen
<div><div dir="ltr">Hello all,<div><br></div>
<div>
<span>Is There any way to run an example of Geant4 ,B1 for example, from a new project in &nbsp;eclipse , which are the libraries i should add to this &nbsp;project? Is there any specific configuration for eclipse?What should i do ?..</span><br>
</div>
<div><span><br></span></div>
<div><span>Thanks!</span></div>
<div><span>Imen</span></div>
<div hspace="streak-pt-mark">&#5159;</div>
</div></div>
Florian Bretin | 14 Apr 10:55 2015

Dose uncertainty

Dear Gaters, 

I sent an email already a few month back, but I still couldn't solve my
issue (although a friendly GATE user replied!). 

I am running a CT simulation of multiple field of views with a phantom. I
shift the phantom to simulate the translational source movement and
therefore have several dose maps and dose uncertainties per scan. I would
like to combine the uncertainty maps to form one single file expressing the
overall uncertainty of the simulation. All files are Analyze75 images (with
header .hdr and image .img).

I read in all files (dose squared map, "squared" in the equation; dose map,
"mean" in the equation and N being the overall emitted photons) in MATLAB
and fail to compute the final uncertainty or even reproduce the single
uncertainty.
I am using the following equation given in the manual (and in the file:
GateImageWithStatistic.cc, line 346, after Chetty et al. from 2006):

    // Chetty2006 p1250 : relative statistical uncertainty
    // exactly same than Ma2002
    if (mean != 0.0 && N != 1 && squared != 0.0){
      *po = sqrt( (1.0/(N-1))*(squared/N - pow(mean/N, 2)))/(mean/N);
    }
    else *po = 1;

All I get is a uncertainty map with complex numbers in MATLAB.
Did anyone face a similar problem?
Does anyone have a tool to compute the uncertainty or the final uncertainty
of several dose maps? 

Any help would be much appreciated, deadlines are approaching! ;) 

Kind regards,
Florian

---
This email is free from viruses and malware because avast! Antivirus protection is active.
http://www.avast.com

BERGAMASCHI Antoine | 14 Apr 10:05 2015
Picon

Root not working in VGate 3.0

Hi all,

Everything is in the title, i am using the Gate V delivered as the version 3.0 and the root programme crash
when using the "draw()" method due to the lack or a corrupted "libjpeg.so.62" and the fail to load
libASImage.so ( which is present in the directory ).

I did not change a single bit of configuration on the VM.

Is a working version of the very usefull vGate disponible anywhere to download ?

Best regards,

Antoine Bergamaschi
AYCIN BORA | 13 Apr 16:57 2015
Picon

Dose Actor Problem

Dear Gate Users;

I'm new at Gate.I'm trying to set up a simulation about radiotheraphy
application.I confine my source distribution to my volume-unitvolume-
and attach the dose actor to that volume. But I get no output. Would
you please help me?

Kind Regards

Aycin BORA
Matthieu Dondey | 10 Apr 16:01 2015
Picon

Parallel simulations and hardware specifications

​​

Dear Gaters,

 

We expect to soon have some funds in our lab to buy a simulation workstation and one of its purposes will be Gate Monte Carlo simulations.

I’ve tried to gather what the important hardware specifications are but I’d like to have advice from people who are already performing parallel simulations on multi-thread workstation or clusters.

My first question is about methodology:

To split a simulation into several jobs when using a multi-core computer, is the process the same as the one described in the Users guide - Parallel Computing -> How to use Gate on a cluster?

Next, is it really beneficial to use GPU + CPU instead of CPUs alone? Is the Nvidia CUDA technology well adapted to Gate simulations and what are its benefits and downsides?

Lastly, I’ve read in the Users Guide:

In planning simulation time, it is important to be aware that Gate simulations currently seem to benefit only from the addition of physical CPUs. A computer with 8 hyper-threaded physical CPU cores (16 logical CPUs) will have the same computational efficiency if 8 processes are run simultaneously as it would with 16 simultaneous processes.

 

Can someone confirm and explain why it is impossible to use logical threads with Gate (“seem to benefit” doesn’t seem realy categorical)? Is it an issue which will be fixed in next versions of GATE? Is it a Geant4 specific problem? Or is it possible to exploit logical threads on Geant4?


Any help will be greatly appreciated and I'd be very grateful about any other advice on hardware specifications (RAM usage, important spec of processors, benefit of a SSD, etc.) 


Cheers,


Matthieu

<div><div dir="ltr">
<div class="gmail_default">&#8203;&#8203;</div>
<div class="gmail_default">

<p class="MsoNormal"><span lang="EN-US">Dear Gaters,</span></p>

<p class="MsoNormal"><span lang="EN-US">&nbsp;</span></p>

<p class="MsoNormal"><span lang="EN-US">We expect
to soon have some funds in our lab to buy a simulation workstation and one of
its purposes will be Gate Monte Carlo simulations.</span></p>

<p class="MsoNormal"><span lang="EN-US">I&rsquo;ve tried
to gather what the important hardware specifications are but I&rsquo;d like to have advice
from people who are already performing parallel simulations on multi-thread
workstation or clusters.</span></p>

<p class="MsoNormal"><span lang="EN-US">My first
question is about methodology:</span></p>

<p class="MsoNormal"><span lang="EN-US">To split a
simulation into several jobs when using a multi-core computer, is the process
the same as the one described in the Users guide - Parallel Computing -&gt; How
to use Gate on a cluster?</span></p>

<p class="MsoNormal"><span lang="EN-US">Next, is it
really beneficial to use GPU + CPU instead of CPUs alone? Is the Nvidia CUDA
technology well adapted to Gate simulations and what are its benefits and
downsides?</span></p>

<p class="MsoNormal"><span lang="EN-US">Lastly, I&rsquo;ve
read in the Users Guide: </span></p>

<p class="MsoNormal" align="left"><span lang="EN-US">In planning
simulation time, it is important to be aware that Gate simulations currently
seem to benefit only from the addition of physical CPUs. A computer with 8
hyper-threaded physical CPU cores (16 logical CPUs) will have the same
computational efficiency if 8 processes are run simultaneously as it would with
16 simultaneous processes.</span></p>

<p class="MsoNormal"><span lang="EN-US">&nbsp;</span></p>

<p class="MsoNormal"><span lang="EN-US">Can someone
confirm and explain why it is impossible to use logical threads with Gate (&ldquo;seem
to benefit&rdquo; doesn&rsquo;t seem realy categorical)? Is it an issue which will be fixed
in next versions of GATE? Is it a Geant4 specific problem? Or is it possible to
exploit logical threads on Geant4?</span></p>
<p class="MsoNormal"><br><span lang="EN-US"></span></p>
<p class="MsoNormal"><span lang="EN-US"></span></p>

<p class="MsoNormal"><span lang="EN-US">Any help
will be greatly appreciated and I'd be very grateful about any other advice on hardware specifications (RAM usage, important spec of processors, benefit of a SSD, etc.)&nbsp; <br></span></p>
<p class="MsoNormal"><br></p>
<p class="MsoNormal">Cheers,<br></p>
<p class="MsoNormal"><span lang="EN-US"></span></p>
<p class="MsoNormal"><span lang="EN-US"><br></span></p>
<p class="MsoNormal"><span lang="EN-US">Matthieu</span></p>

</div>
</div></div>
Grevillot Loic | 10 Apr 15:57 2015
Picon

Re: Retrieving Fluorescence photon position of production

Hello,

Please always answer to the whole mailing list to facilitate tracking.
I would advise you to use root output format instead of IAEA.
You can find a lot of related matters in the attached e-mail, see below copy/pasted relevant matter from
this e-mail.

> If you still need IAEA, from what I remember there is an issue with
> the IAEA format (corrupted Phase space or else). Sometimes ago,
> "Stahis" (look on the mailing-list exchanges) made a lot of posts
> about that and proposed a solution. I attached one e-mail if that can
> help. In future someone should implement correctly a patch for IAEA
> PhS in Gate - if you can do it that would be a very good contribution!

Hope this help.

Sincerely,
Loïc

-----Original Message-----
From: BERGAMASCHI Antoine [mailto:antoine.bergamaschi@...]
Sent: Friday, April 10, 2015 2:21 PM
To: Grevillot Loic
Subject: RE:[Gate-users] Retrieving Fluorescence photon position of production

Hi Loic thx for the answer.

A new problem appears, what to do with this AIEAphsp file ? how can i open it and read it ?

Regards,

Antoine
________________________________________
De : Grevillot Loic [loic.grevillot@...] Envoyé : vendredi 10
avril 2015 12:39 À : BERGAMASCHI Antoine; gate-users@...ration.org
Objet : RE: [Gate-users] Retrieving Fluorescence photon position of     production

Hello,

Phase space actor allows retrieving particle production volume. It might corresponds to you needs,
please have a look at the documentation.

Sincerely,
Loïc

-----Original Message-----
From: gate-users-bounces@...
[mailto:gate-users-bounces@...]
On Behalf Of BERGAMASCHI Antoine
Sent: Friday, April 10, 2015 10:52 AM
To: gate-users@...
Subject: Re: [Gate-users] Retrieving Fluorescence photon position of production

up,

Have anybody any hints ?

Regards,

Antoine Bergamaschi
________________________________________
De : gate-users-bounces@...
[gate-users-bounces <at> lists.opengatecollaboration.org] de la part de BERGAMASCHI Antoine Envoyé :
mercredi 8 avril 2015 15:06 À : gate-users@...ion.org
Objet : [Gate-users] Retrieving Fluorescence photon position of production

Hi all Gaters,

I am a beginner in Gate. I want to simulate self absorption phenomena in phantoms.
My goal is to produce a 3D map representing the position of production of each fluorescence photon hitting
my photon detector.
I already tested using Actors and Filters on the phantom itself, which do work has a first "filter" but the
result is acquiring the "total" fluorescence signal without self absorption.

My questions are:
- Is it possible to retrieve the position of creation of a given particle while only using the scripting
tools available in gate. ( for example given a photon detected in a Volume "A", tell if its come from the
primary source or any other Volume "X" )

- And if it is, what kind of strategy should i use to perform this.

Kind regards

Antoine Bergamaschi
_______________________________________________
Gate-users mailing list
Gate-users@...
http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
_______________________________________________
Gate-users mailing list
Gate-users@...
http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
Disclaimer:
Please notice our E-Mail Disclaimer http://www.medaustron.at/email-disclaimer/

Disclaimer:
Please notice our E-Mail Disclaimer http://www.medaustron.at/email-disclaimer/

Picon
From: Stathis Kamperis <ekamperi@...>
Subject: Re: [Gate-users] Can't use IAEA phase space
Date: 2014-10-28 00:57:43 GMT
Hi Davis,

I also suggest that you just use .root. The IAEA code has many issues...

That said, I think that the problem you are encountering isn't related
to the errors I analyzed in my previous emails. I think that your
problem might be just a typographic error in the code (I'm using Gate
6.2):

stathis:Downloads/gate_v6.2% cat -n
./source/digits_hits/src/GatePhaseSpaceActor.cc | grep -A5 '// gamma'
   279      if( pdg == 22) pIAEARecordType->particle = 1; // gamma
   280      else if( pdg == 11) pIAEARecordType->particle = 2; // electron
   281      else if( pdg == -11) pIAEARecordType->particle = 3; // positron
   282      else if( pdg == 2112) pIAEARecordType->particle = 4; // neutron
   283      else if( pdg == 2122) pIAEARecordType->particle = 5; // proton
   284      else GateError("Actor phase space: particle not available
in IAEA format." );

At line 283 the code checks whether the particle is a proton by
comparing its PDG code (Particle Data Group, IIRC)  against the
constant 2122, which supposedly corresponds to protons. If you check
the Geant4 documentation though[1], you will see that 2122 corresponds
to some other baryon, not proton. Proton has a PDG code 2212, so I
suggest that you replace line 283:

283      else if( pdg == 2122) pIAEARecordType->particle = 5; // proton

with this:

283      else if( pdg == 2212) pIAEARecordType->particle = 5; // proton

rebuild Gate and try again. This is why when you are using an e-
source, everything works well. Because the check for electrons is
correct since electrons do have a PDG code 11.

Disclaimer: I haven't read the code thoroughly, so I might be wrong or
there might be other issues as well. The lack of a consistent coding
style makes the code review a tiresome process (an issue I also had
with Geant4 as well).

Good luck.

Best regards,
Stathis Kamperis


____
[1] http://geant4.cern.ch/G4UsersDocuments/UsersGuides/ForApplicationDeveloper/html/Appendix/particleList.src/baryons/index.html



On Mon, Oct 27, 2014 at 10:15 AM, Grevillot Loic
<loic.grevillot <at> medaustron.at> wrote:
> Hi David,
>
>
>
> For phase spaces, I would suggest to use root format instead of IAEA (just
> change the output to .root).
>
> If you still need IAEA, from what I remember there is an issue with the IAEA
> format (corrupted Phase space or else). Sometimes ago, “Stahis” (look on the
> mailing-list exchanges) made a lot of posts about that and proposed a
> solution. I attached one e-mail if that can help. In future someone should
> implement correctly a patch for IAEA PhS in Gate – if you can do it that
> would be a very good contribution!
>
>
>
> Sincerely,
>
> Loïc
>
>
>
>
>
>
>
>
>
> From: gate-users-bounces <at> lists.opengatecollaboration.org
> [mailto:gate-users-bounces <at> lists.opengatecollaboration.org] On Behalf Of
> David Watts
> Sent: Mittwoch, 22. Oktober 2014 11:53
> To: gate-users <at> lists.opengatecollaboration.org
> Subject: [Gate-users] Can't use IAEA phase space
>
>
>
> Hi Gate community,
>
>
>
> I’m doing some tests recently with Phase Spaces. In particular, I’m trying
> to generate a IAEA Phase Space.
>
>
>
> I’ve just had the error, “src/GatePhaseSpaceActor.cc (l.267): Actor phase
> space: particle not available in IAEA format.”
>
>
>
> This is odd because I am generating a beam of protons using the gps source…
> like this:
>
>                 /gate/source/addSource PBS-source gps
>
>                 /gate/source/PBS-source/gps/particle proton
>
>                 /gate/source/PBS-source/gps/pos/type Beam
>
>                 /gate/source/PBS-source/gps/pos/rot1 0 1 0
>
>                 /gate/source/PBS-source/gps/pos/rot2 1 0 0
>
>                 /gate/source/PBS-source/gps/pos/shape Circle
>
>                 /gate/source/PBS-source/gps/pos/centre 0 0 0 mm
>
>                 /gate/source/PBS-source/gps/pos/sigma_x 0.5 mm
>
>                 /gate/source/PBS-source/gps/pos/sigma_y 0.5 mm
>
>                 /gate/source/PBS-source/gps/ene/mono 70 MeV
>
>                 /gate/source/PBS-source/gps/ene/sigma 0 MeV
>
>                 /gate/source/PBS-source/gps/ene/type Gauss
>
>                 /gate/source/PBS-source/gps/direction 0 0 1
>
>                 /gate/application/setTotalNumberOfPrimaries 10000
>
>
>
> If I change the second line to ‘e-‘ , the simulation runs fine, and the
> output IAEA PS is generated (along with the header file).
>
>
>
> Can someone help please? I’m running the GATE virtual image “vGATE_2.0.vdi”
> on Virtual Box
>
>
>
> In order to avoid starting a new thread on a similar topic, I will include
> another question. What I’m trying to do is develop a method for importing
> particle distributions from beam dynamics codes (such as LINAC… Parmela
> format) as a PS source in GATE. If someone has already done this and is
> willing to share the solution that would be great!
>
>
>
>
>
> Cheers,
>
> Dave.
>
> Disclaimer:
> Please notice our E-Mail Disclaimer
> http://www.medaustron.at/email-disclaimer/

>
>
>
> _______________________________________________
> Gate-users mailing list
> Gate-users <at> lists.opengatecollaboration.org
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users

Hi Davis,

I also suggest that you just use .root. The IAEA code has many issues...

That said, I think that the problem you are encountering isn't related
to the errors I analyzed in my previous emails. I think that your
problem might be just a typographic error in the code (I'm using Gate
6.2):

stathis:Downloads/gate_v6.2% cat -n
./source/digits_hits/src/GatePhaseSpaceActor.cc | grep -A5 '// gamma'
   279      if( pdg == 22) pIAEARecordType->particle = 1; // gamma
   280      else if( pdg == 11) pIAEARecordType->particle = 2; // electron
   281      else if( pdg == -11) pIAEARecordType->particle = 3; // positron
   282      else if( pdg == 2112) pIAEARecordType->particle = 4; // neutron
   283      else if( pdg == 2122) pIAEARecordType->particle = 5; // proton
   284      else GateError("Actor phase space: particle not available
in IAEA format." );

At line 283 the code checks whether the particle is a proton by
comparing its PDG code (Particle Data Group, IIRC)  against the
constant 2122, which supposedly corresponds to protons. If you check
the Geant4 documentation though[1], you will see that 2122 corresponds
to some other baryon, not proton. Proton has a PDG code 2212, so I
suggest that you replace line 283:

283      else if( pdg == 2122) pIAEARecordType->particle = 5; // proton

with this:

283      else if( pdg == 2212) pIAEARecordType->particle = 5; // proton

rebuild Gate and try again. This is why when you are using an e-
source, everything works well. Because the check for electrons is
correct since electrons do have a PDG code 11.

Disclaimer: I haven't read the code thoroughly, so I might be wrong or
there might be other issues as well. The lack of a consistent coding
style makes the code review a tiresome process (an issue I also had
with Geant4 as well).

Good luck.

Best regards,
Stathis Kamperis


____
[1] http://geant4.cern.ch/G4UsersDocuments/UsersGuides/ForApplicationDeveloper/html/Appendix/particleList.src/baryons/index.html



On Mon, Oct 27, 2014 at 10:15 AM, Grevillot Loic
<loic.grevillot <at> medaustron.at> wrote:
> Hi David,
>
>
>
> For phase spaces, I would suggest to use root format instead of IAEA (just
> change the output to .root).
>
> If you still need IAEA, from what I remember there is an issue with the IAEA
> format (corrupted Phase space or else). Sometimes ago, “Stahis” (look on the
> mailing-list exchanges) made a lot of posts about that and proposed a
> solution. I attached one e-mail if that can help. In future someone should
> implement correctly a patch for IAEA PhS in Gate – if you can do it that
> would be a very good contribution!
>
>
>
> Sincerely,
>
> Loïc
>
>
>
>
>
>
>
>
>
> From: gate-users-bounces <at> lists.opengatecollaboration.org
> [mailto:gate-users-bounces <at> lists.opengatecollaboration.org] On Behalf Of
> David Watts
> Sent: Mittwoch, 22. Oktober 2014 11:53
> To: gate-users <at> lists.opengatecollaboration.org
> Subject: [Gate-users] Can't use IAEA phase space
>
>
>
> Hi Gate community,
>
>
>
> I’m doing some tests recently with Phase Spaces. In particular, I’m trying
> to generate a IAEA Phase Space.
>
>
>
> I’ve just had the error, “src/GatePhaseSpaceActor.cc (l.267): Actor phase
> space: particle not available in IAEA format.”
>
>
>
> This is odd because I am generating a beam of protons using the gps source…
> like this:
>
>                 /gate/source/addSource PBS-source gps
>
>                 /gate/source/PBS-source/gps/particle proton
>
>                 /gate/source/PBS-source/gps/pos/type Beam
>
>                 /gate/source/PBS-source/gps/pos/rot1 0 1 0
>
>                 /gate/source/PBS-source/gps/pos/rot2 1 0 0
>
>                 /gate/source/PBS-source/gps/pos/shape Circle
>
>                 /gate/source/PBS-source/gps/pos/centre 0 0 0 mm
>
>                 /gate/source/PBS-source/gps/pos/sigma_x 0.5 mm
>
>                 /gate/source/PBS-source/gps/pos/sigma_y 0.5 mm
>
>                 /gate/source/PBS-source/gps/ene/mono 70 MeV
>
>                 /gate/source/PBS-source/gps/ene/sigma 0 MeV
>
>                 /gate/source/PBS-source/gps/ene/type Gauss
>
>                 /gate/source/PBS-source/gps/direction 0 0 1
>
>                 /gate/application/setTotalNumberOfPrimaries 10000
>
>
>
> If I change the second line to ‘e-‘ , the simulation runs fine, and the
> output IAEA PS is generated (along with the header file).
>
>
>
> Can someone help please? I’m running the GATE virtual image “vGATE_2.0.vdi”
> on Virtual Box
>
>
>
> In order to avoid starting a new thread on a similar topic, I will include
> another question. What I’m trying to do is develop a method for importing
> particle distributions from beam dynamics codes (such as LINAC… Parmela
> format) as a PS source in GATE. If someone has already done this and is
> willing to share the solution that would be great!
>
>
>
>
>
> Cheers,
>
> Dave.
>
> Disclaimer:
> Please notice our E-Mail Disclaimer
> http://www.medaustron.at/email-disclaimer/

>
>
>
> _______________________________________________
> Gate-users mailing list
> Gate-users <at> lists.opengatecollaboration.org
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users


Gmane