GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Bugs, Fixes, Releases » negative masses in geometry (TG3)
Re: negative masses in geometry (TG3) [message #11691 is a reply to message #11688] Wed, 20 April 2011 15:19 Go to previous messageGo to previous message
Mohammad Al-Turany is currently offline  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




 
Read Message icon5.gif
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message icon14.gif
Previous Topic: TPC compilation error
Next Topic: STT tracking crash during check getMCInfo.
Goto Forum:
  


Current Time: Fri Dec 03 01:47:11 CET 2021

Total time taken to generate the page: 0.02250 seconds