Message: Re: Problem with G4Polycone Not Logged In (login)

## Re: Problem with G4Polycone

Forum: Geometry
Re: Problem with G4Polycone (Filippo Ambroglini)
Re: Re: Problem with G4Polycone (Filippo Ambroglini)
Re: Re: Problem with G4Polycone (John Apostolakis)
Re: Re: Problem with G4Polycone (Filippo Ambroglini)
Date: 14 Jan, 2013
From: Filippo Ambroglini <Filippo Ambroglini>

 ```Dear John, any news on this. Best regards, Filippo > *** Discussion title: Geometry > > Dear John, > > the order of the Z was descendig then I have reverse to ascending but I > still continue to have the problem. > > I have tested the solid with these: > if(physiVanAllenBelt->IsRegularStructure()) > G4cout<<"The solid is not a regular structure"< if(physiVanAllenBelt->CheckOverlaps()) > G4cout<<"The solid had some overlaps"< > but no error reported and trying to visualize it it seems ok. > > In attach you can find the list of z Rmin and Rmax and also an image > from visualization. > > Thanks for the help. > > Filippo > > P.S. > > The value fo Z Rmin and Rmax are listed normalized to rplanet ( earth > radius ) > >> Dear Filippo, >> >> Can you double check that your points are in ascending order - i.e. z0 < z1 < z2 < z3 < z4 < z5 < z6 < z7 ... >> >> r_i < R_i ( where r is r_min and R is r_max for location 'i' ) >> >> If this is ok, can you give a well-formatted printout of the values of < z, Rmin, Rmax > ? >> >> Also, have you tried to test the solid that you created ? >> >> Best regards, >> John >> =================================================== >> John Apostolakis, SFT (SoFTware) Group, PH Department, CERN >> Email: john.apostolakis@cern.ch Office Tel: +41-22-767-7239 >> ------------------------------------------------------------------------------------------ >> >> On Dec 5, 2012, at 3:21 PM, Filippo Ambroglini wrote: >> >>> *** Discussion title: Geometry >>> >>> Dear all, >>> >>> please some one could try to help me on understand why I had the below >>> mentioned problem with G4Polycone? >>> >>> Best regards, Filippo On Thu, 15 Nov 2012 10:10:47 GMT, Filippo >>> Ambroglini wrote: >>> >>>> Dear G4 expert, >>>> >>>> I had a problem with a G4Polycone object that I create in order to >>>> simulate a complex volume like an elliptical torus with an elliptical >>>> hole in the central part. >>>> >>>> This is the implementation: >>>> >>>> std::vector tempZ, tempInner,tempOuter; >>>> >>>> for(int j=0; j<9;j++){ >>>> tempZ.push_back((2-0.1*j)*rplanet); >>>> tempInner.push_back((2.5-sqrt(1-((1.9-0.1*j)*(1.9-0.1*j))/4)*2.5)*rplanet); >>>> tempOuter.push_back((2.5+sqrt(1-((1.9-0.1*j)*(1.9-0.1*j))/4)*2.5)*rplanet); >>>> } >>>> for(int j=0;j<23;j++){ >>>> tempZ.push_back((1.1-0.1*j)*rplanet); >>>> tempInner.push_back((sqrt(1-((1.1-0.1*j)*(1.1-0.1*j))/(1.2*1.2))*1.2)*rplanet); >>>> tempOuter.push_back((2.5+sqrt(1-((1.1-0.1*j)*(1.1-0.1*j))/4)*2.5)*rplanet); >>>> } >>>> for(int j=0; j<9;j++){ >>>> tempZ.push_back((-1.1-0.1*j)*rplanet); >>>> tempInner.push_back((2.5-sqrt(1-((1.1+0.1*j)*(1.1+0.1*j))/4)*2.5)*rplanet); >>>> tempOuter.push_back((2.5+sqrt(1-((1.1+0.1*j)*(1.1+0.1*j))/4)*2.5)*rplanet); >>>> } >>>> G4double rInner[41],rOuter[41],zPlane[41]; >>>> for(int i=0;i<41;i++){ >>>> rInner[i]=tempInner[i]; >>>> rOuter[i]=tempOuter[i]; >>>> zPlane[i]=tempZ[i]; >>>> } >>>> G4Polycone* VanAllenBelt = new >>>> G4Polycone("VanAllenBelt",0,2*pi,41,zPlane,rInner,rOuter); >>>> logicVanAllenBelt = new >>>> G4LogicalVolume(VanAllenBelt,Vacuum,"VanAllenBelt",0,0,0); >>>> >>>> logicVanAllenBelt->SetVisAttributes(VisAttInvisible); >>>> logicVanAllenBelt->SetUserLimits(theMagnetosphereUserLimits); >>>> >>>> G4VPhysicalVolume* physiVanAllenBelt= new >>>> G4PVPlacement(0,G4ThreeVector(),"VanAllenBeltPV",logicVanAllenBelt,physiMagnetosphere,false,0); >>>> >>>> If I try to visualize it everything seams ok but when I try to use it in >>>> my program and I will try to transport a particle inside this volume >>>> accessing the volume name with >>>> aStep->GetPreStepPoint()->GetPhysicalVolume()->GetName(); I always get >>>> the name of the mother volume. >>>> >>>> If instead of G4Polycone I will use this definition >>>> >>>> G4Torus* outerInnerVAB = new >>>> G4Torus("ExternInnerVAB",0,1.25*rplanet,0.75*rplanet+ZpositionForAtmosphereBottom,0,2*pi); >>>> G4Orb* innerInnerVAB = new >>>> G4Orb("InnerInnerVAB",0.2*rplanet+ZpositionForAtmosphereBottom); >>>> G4SubtractionSolid* VanAllenBelt = new >>>> G4SubtractionSolid("VanAllenBelt",outerInnerVAB,innerInnerVAB); >>>> >>>> everything is ok. >>>> >>>> Someone could explain me why and where is my mistake? >>>> >>>> Best regards, >>>> Filippo >>>> >>> ------------------------------------------------------------- >>> Visit this GEANT4 at hypernews.slac.stanford.edu message (to reply or unsubscribe) at: >>> http://hypernews.slac.stanford.edu/HyperNews/geant4/get/geometry/1242/1.html > Z[0] = -2; rInner[0] = 1.71938; rOuter[0] = 3.28062 > Z[1] = -1.9; rInner[1] = 1.41028; rOuter[1] = 3.58972 > Z[2] = -1.8; rInner[2] = 1.18304; rOuter[2] = 3.81696 > Z[3] = -1.7; rInner[3] = 1; rOuter[3] = 4 > Z[4] = -1.6; rInner[4] = 0.846405; rOuter[4] = 4.15359 > Z[5] = -1.5; rInner[5] = 0.714643; rOuter[5] = 4.28536 > Z[6] = -1.4; rInner[6] = 0.600164; rOuter[6] = 4.39984 > Z[7] = -1.3; rInner[7] = 0.5; rOuter[7] = 4.5 > Z[8] = -1.2; rInner[8] = 0.412088; rOuter[8] = 4.58791 > Z[9] = -1.1; rInner[9] = 0.479583; rOuter[9] = 4.58791 > Z[10] = -1; rInner[10] = 0.663325; rOuter[10] = 4.66506 > Z[11] = -0.9; rInner[11] = 0.793725; rOuter[11] = 4.73257 > Z[12] = -0.8; rInner[12] = 0.894427; rOuter[12] = 4.79129 > Z[13] = -0.7; rInner[13] = 0.974679; rOuter[13] = 4.84187 > Z[14] = -0.6; rInner[14] = 1.03923; rOuter[14] = 4.88485 > Z[15] = -0.5; rInner[15] = 1.09087; rOuter[15] = 4.92061 > Z[16] = -0.4; rInner[16] = 1.13137; rOuter[16] = 4.94949 > Z[17] = -0.3; rInner[17] = 1.1619; rOuter[17] = 4.97171 > Z[18] = -0.2; rInner[18] = 1.18322; rOuter[18] = 4.98747 > Z[19] = -0.1; rInner[19] = 1.19583; rOuter[19] = 4.99687 > Z[20] = 0; rInner[20] = 1.2; rOuter[20] = 5 > Z[21] = 0.1; rInner[21] = 1.19583; rOuter[21] = 4.99687 > Z[22] = 0.2; rInner[22] = 1.18322; rOuter[22] = 4.98747 > Z[23] = 0.3; rInner[23] = 1.1619; rOuter[23] = 4.97171 > Z[24] = 0.4; rInner[24] = 1.13137; rOuter[24] = 4.94949 > Z[25] = 0.5; rInner[25] = 1.09087; rOuter[25] = 4.92061 > Z[26] = 0.6; rInner[26] = 1.03923; rOuter[26] = 4.88485 > Z[27] = 0.7; rInner[27] = 0.974679; rOuter[27] = 4.84187 > Z[28] = 0.8; rInner[28] = 0.894427; rOuter[28] = 4.79129 > Z[29] = 0.9; rInner[29] = 0.793725; rOuter[29] = 4.73257 > Z[30] = 1; rInner[30] = 0.663325; rOuter[30] = 4.66506 > Z[31] = 1.1; rInner[31] = 0.479583; rOuter[31] = 4.58791 > Z[32] = 1.2; rInner[32] = 0.412088; rOuter[32] = 4.58791 > Z[33] = 1.3; rInner[33] = 0.5; rOuter[33] = 4.5 > Z[34] = 1.4; rInner[34] = 0.600164; rOuter[34] = 4.39984 > Z[35] = 1.5; rInner[35] = 0.714643; rOuter[35] = 4.28536 > Z[36] = 1.6; rInner[36] = 0.846405; rOuter[36] = 4.15359 > Z[37] = 1.7; rInner[37] = 1; rOuter[37] = 4 > Z[38] = 1.8; rInner[38] = 1.18304; rOuter[38] = 3.81696 > Z[39] = 1.9; rInner[39] = 1.41028; rOuter[39] = 3.58972 > Z[40] = 2; rInner[40] = 1.71938; rOuter[40] = 3.28062 > > Attachment: > http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2012/12/05/07.57-51810-VolumeView.jpg > > ------------------------------------------------------------- > Visit this GEANT4 at hypernews.slac.stanford.edu message (to reply or unsubscribe) at: > http://hypernews.slac.stanford.edu/HyperNews/geant4/get/geometry/1242/1/1/1.html ```

 Inline Depth: Outline Depth: Add message:

1 Re: Problem with G4Polycone   (John Allison - 14 Jan, 2013)
Re: Problem with G4Polycone   (Filippo Ambroglini - 15 Jan, 2013)
Re: Problem with G4Polycone   (John Allison - 15 Jan, 2013)
 to: "Re: Problem with G4Polycone"
This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews

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