GSI Forum: Bugs, Fixes, Releases » negative masses in geometry (TG3)

Home » PANDA » PandaRoot » Bugs, Fixes, Releases » negative masses in geometry (TG3)

Show: Today's Messages :: Polls :: Message Navigator

Re: negative masses in geometry (TG3) [message #11691 is a reply to message #11688]

Wed, 20 April 2011 15:19

Mohammad Al-Turany
Messages: 518
Registered: April 2004
Location: GSI, Germany

first-grade participant

From: *gsi.de

Hi Olaf,

it could be a bug in TGeoXtru::Capacity:

Quote:

Double_t TGeoXtru::Capacity() const
{
// Compute capacity [length^3] of this shape.
Int_t iz;
Double_t capacity = 0;
Double_t area, dz, sc1, sc2;
TGeoXtru *xtru = (TGeoXtru*)this;
xtru->SetCurrentVertices(0.,0.,1.);
area = fPoly->Area();
for (iz=0; iz<fNz-1; iz++) {
dz = fZ[iz+1]-fZ[iz];
if (TGeoShape::IsSameWithinTolerance(dz,0)) continue;
sc1 = fScale[iz];
sc2 = fScale[iz+1];
capacity += (area*dz/3.)*(sc1*sc1+sc1*sc2+sc2*sc2);
}
return capacity;
}

in the line dz = fZ[iz+1]-fZ[iz]; if fZ[iz+1] is smaller than fZ[iz] then you get a negative value. Anyway I did not find the volume "Supporto2" in the geometry. can you please tell me where to find it or post a root file with a TGeoManager having this problem.

regards

Mohammad

Report message to a moderator

[Message index]

		negative masses in geometry (TG3) By: Olaf Hartmann on Wed, 20 April 2011 14:36
		Re: negative masses in geometry (TG3) By: Mohammad Al-Turany on Wed, 20 April 2011 15:19
		Re: negative masses in geometry (TG3) By: Olaf Hartmann on Wed, 20 April 2011 16:10
		Re: negative masses in geometry (TG3) By: Mohammad Al-Turany on Wed, 20 April 2011 22:56
		Re: negative masses in geometry (TG3) By: StefanoSpataro on Thu, 21 April 2011 08:00
		Re: negative masses in geometry (TG3) By: Olaf Hartmann on Thu, 21 April 2011 08:35
		Re: negative masses in geometry (TG3) By: Olaf Hartmann on Thu, 21 April 2011 10:26

Previous Topic:	TPC compilation error
Next Topic:	STT tracking crash during check getMCInfo.

Goto Forum:

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sun Oct 01 23:53:16 CEST 2023

Total time taken to generate the page: 0.00889 seconds