From 6a90570870c0d17eca5e3bcba6cd3e13cf61526b Mon Sep 17 00:00:00 2001 From: kowal2 Date: Thu, 2 Nov 2000 07:24:11 +0000 Subject: [PATCH] Correction in the TPC geometry. Changes due to the new hit structure. --- TPC/AliTPCv0.cxx | 51 ++++++++++++-------------------- TPC/AliTPCv1.cxx | 76 +++++++++++++++++++++++++----------------------- TPC/AliTPCv1.h | 2 +- TPC/AliTPCv2.cxx | 72 ++++++++++++++++++++++----------------------- TPC/AliTPCv3.cxx | 53 ++++++++++++--------------------- 5 files changed, 113 insertions(+), 141 deletions(-) diff --git a/TPC/AliTPCv0.cxx b/TPC/AliTPCv0.cxx index 86fb39a4aa9..606594598b6 100644 --- a/TPC/AliTPCv0.cxx +++ b/TPC/AliTPCv0.cxx @@ -15,6 +15,9 @@ /* $Log$ +Revision 1.17 2000/07/10 20:57:39 hristov +Update of TPC code and macros by M.Kowalski + Revision 1.16 2000/06/30 12:07:50 kowal2 Updated from the TPC-PreRelease branch @@ -423,38 +426,31 @@ void AliTPCv0::CreateGeometry() // Daughter volumes - // Tpc SAndwich 5 - Al + // Tpc SAndwich 5 - Tedlar dm[0]= 258.; dm[1]= 260.05; dm[2]= 251.7; - gMC->Gsvolu("TSA5","TUBE",idtmed[4],dm,3); + gMC->Gsvolu("TSA5","TUBE",idtmed[9],dm,3); - // Tpc SAndwich 6 - Tedlar + // Tpc SAndwich 6 - Kevlar dm[0] += 5.e-3; dm[1] -= 5.e-3; - gMC->Gsvolu("TSA6","TUBE",idtmed[9],dm,3); - - // Tpc SAndwich 7 - Kevlar - - dm[0] += 5.e-3; - dm[1] -= 5.e-3; + gMC->Gsvolu("TSA6","TUBE",idtmed[5],dm,3); - gMC->Gsvolu("TSA7","TUBE",idtmed[5],dm,3); - // Tpc SAndwich 8 - NOMEX + // Tpc SAndwich 7 - NOMEX dm[0] += 0.02; dm[1] -= 0.02; - gMC->Gsvolu("TSA8","TUBE",idtmed[6],dm,3); + gMC->Gsvolu("TSA7","TUBE",idtmed[6],dm,3); - // 8->7->6->5->TOFC + // 7->6->5->TOFC - gMC->Gspos("TSA8",1,"TSA7",0.,0.,0.,0,"ONLY"); gMC->Gspos("TSA7",1,"TSA6",0.,0.,0.,0,"ONLY"); gMC->Gspos("TSA6",1,"TSA5",0.,0.,0.,0,"ONLY"); @@ -892,38 +888,31 @@ void AliTPCv0::CreateGeometry() // Daughter volumes - // Tpc Sandwich 17 - Al + // Tpc Sandwich 17 - Tedlar dm[0]= 77.15; dm[1]= 79.2; dm[2]= 251.7; - gMC->Gsvolu("TS17","TUBE",idtmed[4],dm,3); + gMC->Gsvolu("TS17","TUBE",idtmed[9],dm,3); - // Tpc Sandwich 18 - Tedlar + // Tpc Sandwich 18 - Kevlar dm[0]+= 5.e-3; dm[1]-= 5.e-3; - gMC->Gsvolu("TS18","TUBE",idtmed[9],dm,3); + gMC->Gsvolu("TS18","TUBE",idtmed[5],dm,3); - // Tpc Sandwich 19 - Kevlar - dm[0]+= 5.e-3; - dm[1]-= 5.e-3; - - gMC->Gsvolu("TS19","TUBE",idtmed[5],dm,3); - - // Tpc Sandwich 20 - NOMEX + // Tpc Sandwich 19 - NOMEX - dm[0]+= 0.06; - dm[1]-= 0.06; + dm[0]+= 0.02; + dm[1]-= 0.02; - gMC->Gsvolu("TS20","TUBE",idtmed[6],dm,3); + gMC->Gsvolu("TS19","TUBE",idtmed[6],dm,3); - // 20->19->18->17 + // 19->18->17 - gMC->Gspos("TS20",1,"TS19",0.,0.,0.,0,"ONLY"); gMC->Gspos("TS19",1,"TS18",0.,0.,0.,0,"ONLY"); gMC->Gspos("TS18",1,"TS17",0.,0.,0.,0,"ONLY"); @@ -1714,7 +1703,6 @@ void AliTPCv0::DrawDetector() gMC->Gsatt("TSA5","SEEN",0); gMC->Gsatt("TSA6","SEEN",0); gMC->Gsatt("TSA7","SEEN",0); - gMC->Gsatt("TSA8","SEEN",0); gMC->Gsatt("TIIN","COLO",7); gMC->Gsatt("TIIN","SEEN",1); gMC->Gsatt("TICL","SEEN",0); @@ -1732,7 +1720,6 @@ void AliTPCv0::DrawDetector() gMC->Gsatt("TS17","SEEN",0); gMC->Gsatt("TS18","SEEN",0); gMC->Gsatt("TS19","SEEN",0); - gMC->Gsatt("TS20","SEEN",0); gMC->Gsatt("TS21","SEEN",0); gMC->Gsatt("TS22","SEEN",0); gMC->Gsatt("TS23","SEEN",0); diff --git a/TPC/AliTPCv1.cxx b/TPC/AliTPCv1.cxx index d5b4e5c5638..36c662aaa22 100644 --- a/TPC/AliTPCv1.cxx +++ b/TPC/AliTPCv1.cxx @@ -15,6 +15,9 @@ /* $Log$ +Revision 1.19 2000/07/10 20:57:39 hristov +Update of TPC code and macros by M.Kowalski + Revision 1.18 2000/06/30 12:07:50 kowal2 Updated from the TPC-PreRelease branch @@ -93,6 +96,14 @@ Introduction of the Copyright and cvs Log #include ClassImp(AliTPCv1) + +//_____________________________________________________________________________ + + AliTPCv1::AliTPCv1(){ + + fHitType = 1; + +} //_____________________________________________________________________________ AliTPCv1::AliTPCv1(const char *name, const char *title) @@ -106,6 +117,8 @@ AliTPCv1::AliTPCv1(const char *name, const char *title) if (fTPCParam) fTPCParam->Write(fTPCParam->GetTitle()); + + fHitType=1; } //_____________________________________________________________________________ @@ -423,38 +436,31 @@ void AliTPCv1::CreateGeometry() // Daughter volumes - // Tpc SAndwich 5 - Al + // Tpc SAndwich 5 - Tedlar dm[0]= 258.; dm[1]= 260.05; dm[2]= 251.7; - gMC->Gsvolu("TSA5","TUBE",idtmed[4],dm,3); + gMC->Gsvolu("TSA5","TUBE",idtmed[9],dm,3); - // Tpc SAndwich 6 - Tedlar + // Tpc SAndwich 6 - Kevlar dm[0] += 5.e-3; dm[1] -= 5.e-3; - gMC->Gsvolu("TSA6","TUBE",idtmed[9],dm,3); - - // Tpc SAndwich 7 - Kevlar - - dm[0] += 5.e-3; - dm[1] -= 5.e-3; + gMC->Gsvolu("TSA6","TUBE",idtmed[5],dm,3); - gMC->Gsvolu("TSA7","TUBE",idtmed[5],dm,3); - // Tpc SAndwich 8 - NOMEX + // Tpc SAndwich 7 - NOMEX dm[0] += 0.02; dm[1] -= 0.02; - gMC->Gsvolu("TSA8","TUBE",idtmed[6],dm,3); + gMC->Gsvolu("TSA7","TUBE",idtmed[6],dm,3); - // 8->7->6->5->TOFC + // 7->6->5->TOFC - gMC->Gspos("TSA8",1,"TSA7",0.,0.,0.,0,"ONLY"); gMC->Gspos("TSA7",1,"TSA6",0.,0.,0.,0,"ONLY"); gMC->Gspos("TSA6",1,"TSA5",0.,0.,0.,0,"ONLY"); @@ -892,38 +898,31 @@ void AliTPCv1::CreateGeometry() // Daughter volumes - // Tpc Sandwich 17 - Al + // Tpc Sandwich 17 - Tedlar dm[0]= 77.15; dm[1]= 79.2; dm[2]= 251.7; - gMC->Gsvolu("TS17","TUBE",idtmed[4],dm,3); + gMC->Gsvolu("TS17","TUBE",idtmed[9],dm,3); - // Tpc Sandwich 18 - Tedlar + // Tpc Sandwich 18 - Kevlar dm[0]+= 5.e-3; dm[1]-= 5.e-3; - gMC->Gsvolu("TS18","TUBE",idtmed[9],dm,3); - - // Tpc Sandwich 19 - Kevlar - - dm[0]+= 5.e-3; - dm[1]-= 5.e-3; + gMC->Gsvolu("TS18","TUBE",idtmed[5],dm,3); - gMC->Gsvolu("TS19","TUBE",idtmed[5],dm,3); - // Tpc Sandwich 20 - NOMEX + // Tpc Sandwich 19 - NOMEX - dm[0]+= 0.06; - dm[1]-= 0.06; + dm[0]+= 0.02; + dm[1]-= 0.02; - gMC->Gsvolu("TS20","TUBE",idtmed[6],dm,3); + gMC->Gsvolu("TS19","TUBE",idtmed[6],dm,3); - // 20->19->18->17 + // 19->18->17 - gMC->Gspos("TS20",1,"TS19",0.,0.,0.,0,"ONLY"); gMC->Gspos("TS19",1,"TS18",0.,0.,0.,0,"ONLY"); gMC->Gspos("TS18",1,"TS17",0.,0.,0.,0,"ONLY"); @@ -1809,7 +1808,6 @@ void AliTPCv1::DrawDetector() gMC->Gsatt("TSA5","SEEN",0); gMC->Gsatt("TSA6","SEEN",0); gMC->Gsatt("TSA7","SEEN",0); - gMC->Gsatt("TSA8","SEEN",0); gMC->Gsatt("TIIN","COLO",7); gMC->Gsatt("TIIN","SEEN",1); gMC->Gsatt("TICL","SEEN",0); @@ -1827,7 +1825,6 @@ void AliTPCv1::DrawDetector() gMC->Gsatt("TS17","SEEN",0); gMC->Gsatt("TS18","SEEN",0); gMC->Gsatt("TS19","SEEN",0); - gMC->Gsatt("TS20","SEEN",0); gMC->Gsatt("TS21","SEEN",0); gMC->Gsatt("TS22","SEEN",0); gMC->Gsatt("TS23","SEEN",0); @@ -1922,8 +1919,7 @@ void AliTPCv1::StepManager() Float_t hits[4]; Int_t vol[2]; TLorentzVector p; - TClonesArray &lhits = *fHits; - // + // if(gMC->TrackCharge() && gMC->IsTrackEntering()) { @@ -1968,21 +1964,27 @@ void AliTPCv1::StepManager() hits[1]=p[1]; hits[2]=p[2]; hits[3]=0.; // this hit has no energy loss - new(lhits[fNhits++]) AliTPChit(fIshunt,gAlice->CurrentTrack(),vol,hits); + // new(lhits[fNhits++]) AliTPChit(fIshunt,gAlice->CurrentTrack(),vol,hits); + + AddHit(gAlice->CurrentTrack(), vol,hits); // M.I. gMC->TrackPosition(p); hits[0]=p[0]; hits[1]=p[1]; hits[2]=p[2]; hits[3]=0.; // this hit has no energy loss - new(lhits[fNhits++]) AliTPChit(fIshunt,gAlice->CurrentTrack(),vol,hits); + // new(lhits[fNhits++]) AliTPChit(fIshunt,gAlice->CurrentTrack(),vol,hits); + + AddHit(gAlice->CurrentTrack(), vol,hits); // M.I. } gMC->TrackPosition(p); for(i=0;i<3;++i) hits[i]=p[i]; hits[3]=1; //I'd like to have something positive here (I.Belikov) - new(lhits[fNhits++]) AliTPChit(fIshunt,gAlice->CurrentTrack(),vol,hits); + // new(lhits[fNhits++]) AliTPChit(fIshunt,gAlice->CurrentTrack(),vol,hits); + + AddHit(gAlice->CurrentTrack(), vol,hits); // M.I. } diff --git a/TPC/AliTPCv1.h b/TPC/AliTPCv1.h index ca45d42069a..6d2b9801e3d 100644 --- a/TPC/AliTPCv1.h +++ b/TPC/AliTPCv1.h @@ -15,7 +15,7 @@ class AliTPCv1 : public AliTPC { public: - AliTPCv1() {} + AliTPCv1(); AliTPCv1(const char *name, const char *title); virtual ~AliTPCv1() {} virtual void CreateGeometry(); diff --git a/TPC/AliTPCv2.cxx b/TPC/AliTPCv2.cxx index 26bf561fe6e..670818912df 100644 --- a/TPC/AliTPCv2.cxx +++ b/TPC/AliTPCv2.cxx @@ -15,6 +15,9 @@ /* $Log$ +Revision 1.26 2000/10/05 16:16:29 kowal2 +Corrections of the hit recording algorithm. + Revision 1.25 2000/10/02 21:28:18 fca Removal of useless dependecies via forward declarations @@ -98,6 +101,8 @@ Introduction of the Copyright and cvs Log #include "AliPDG.h" #include "AliTPCParam.h" #include "AliTPCParamSR.h" +#include "AliTPCTrackHits.h" + ClassImp(AliTPCv2) @@ -487,38 +492,31 @@ void AliTPCv2::CreateGeometry() // Daughter volumes - // Tpc SAndwich 5 - Al + // Tpc SAndwich 5 - Tedlar dm[0]= 258.; dm[1]= 260.05; dm[2]= 251.7; - gMC->Gsvolu("TSA5","TUBE",idtmed[4],dm,3); + gMC->Gsvolu("TSA5","TUBE",idtmed[9],dm,3); - // Tpc SAndwich 6 - Tedlar + // Tpc SAndwich 6 - Kevlar dm[0] += 5.e-3; dm[1] -= 5.e-3; - gMC->Gsvolu("TSA6","TUBE",idtmed[9],dm,3); - - // Tpc SAndwich 7 - Kevlar + gMC->Gsvolu("TSA6","TUBE",idtmed[5],dm,3); - dm[0] += 5.e-3; - dm[1] -= 5.e-3; - - gMC->Gsvolu("TSA7","TUBE",idtmed[5],dm,3); // Tpc SAndwich 8 - NOMEX dm[0] += 0.02; dm[1] -= 0.02; - gMC->Gsvolu("TSA8","TUBE",idtmed[6],dm,3); + gMC->Gsvolu("TSA7","TUBE",idtmed[6],dm,3); - // 8->7->6->5->TOFC + // 7->6->5->TOFC - gMC->Gspos("TSA8",1,"TSA7",0.,0.,0.,0,"ONLY"); gMC->Gspos("TSA7",1,"TSA6",0.,0.,0.,0,"ONLY"); gMC->Gspos("TSA6",1,"TSA5",0.,0.,0.,0,"ONLY"); @@ -956,38 +954,31 @@ void AliTPCv2::CreateGeometry() // Daughter volumes - // Tpc Sandwich 17 - Al + // Tpc Sandwich 17 - Tedlar dm[0]= 77.15; dm[1]= 79.2; dm[2]= 251.7; - gMC->Gsvolu("TS17","TUBE",idtmed[4],dm,3); + gMC->Gsvolu("TS17","TUBE",idtmed[9],dm,3); - // Tpc Sandwich 18 - Tedlar + // Tpc Sandwich 18 - Kevlar dm[0]+= 5.e-3; dm[1]-= 5.e-3; - gMC->Gsvolu("TS18","TUBE",idtmed[9],dm,3); - - // Tpc Sandwich 19 - Kevlar + gMC->Gsvolu("TS18","TUBE",idtmed[5],dm,3); - dm[0]+= 5.e-3; - dm[1]-= 5.e-3; - gMC->Gsvolu("TS19","TUBE",idtmed[5],dm,3); + // Tpc Sandwich 19 - NOMEX - // Tpc Sandwich 20 - NOMEX - - dm[0]+= 0.06; - dm[1]-= 0.06; + dm[0]+= 0.02; + dm[1]-= 0.02; - gMC->Gsvolu("TS20","TUBE",idtmed[6],dm,3); + gMC->Gsvolu("TS19","TUBE",idtmed[6],dm,3); - // 20->19->18->17 + // 19->18->17 - gMC->Gspos("TS20",1,"TS19",0.,0.,0.,0,"ONLY"); gMC->Gspos("TS19",1,"TS18",0.,0.,0.,0,"ONLY"); gMC->Gspos("TS18",1,"TS17",0.,0.,0.,0,"ONLY"); @@ -1915,7 +1906,6 @@ void AliTPCv2::DrawDetector() gMC->Gsatt("TSA5","SEEN",0); gMC->Gsatt("TSA6","SEEN",0); gMC->Gsatt("TSA7","SEEN",0); - gMC->Gsatt("TSA8","SEEN",0); gMC->Gsatt("TIIN","COLO",7); gMC->Gsatt("TIIN","SEEN",1); gMC->Gsatt("TICL","SEEN",0); @@ -1933,7 +1923,6 @@ void AliTPCv2::DrawDetector() gMC->Gsatt("TS17","SEEN",0); gMC->Gsatt("TS18","SEEN",0); gMC->Gsatt("TS19","SEEN",0); - gMC->Gsatt("TS20","SEEN",0); gMC->Gsatt("TS21","SEEN",0); gMC->Gsatt("TS22","SEEN",0); gMC->Gsatt("TS23","SEEN",0); @@ -2055,7 +2044,6 @@ void AliTPCv2::StepManager() Int_t id,copy; Float_t hits[4]; Int_t vol[2]; - TClonesArray &lhits = *fHits; TLorentzVector p; vol[1]=0; // preset row number to 0 @@ -2102,7 +2090,7 @@ void AliTPCv2::StepManager() // upper sector vol[0] = copy-1+fTPCParam->GetNInnerSector(); // sector number - vol[1] -= fTPCParam->GetNRowLow(); // row number + vol[1] -= fTPCParam->GetNRowLow(); // row number (starts also from 0) } @@ -2115,7 +2103,9 @@ void AliTPCv2::StepManager() hits[1]=p[1]; hits[2]=p[2]; hits[3]=0.; // this hit has no energy loss - new(lhits[fNhits++]) AliTPChit(fIshunt,gAlice->CurrentTrack(),vol,hits); + // new(lhits[fNhits++]) AliTPChit(fIshunt,gAlice->CurrentTrack(),vol,hits); + + AddHit(gAlice->CurrentTrack(), vol,hits); //MI change } @@ -2124,9 +2114,10 @@ void AliTPCv2::StepManager() hits[1]=p[1]; hits[2]=p[2]; hits[3]=0.; // this hit has no energy loss - new(lhits[fNhits++]) AliTPChit(fIshunt,gAlice->CurrentTrack(),vol,hits); + // new(lhits[fNhits++]) AliTPChit(fIshunt,gAlice->CurrentTrack(),vol,hits); + + AddHit(gAlice->CurrentTrack(), vol,hits); //MI change - } else return; @@ -2147,7 +2138,14 @@ void AliTPCv2::StepManager() // Add this hit - new(lhits[fNhits++]) AliTPChit(fIshunt,gAlice->CurrentTrack(),vol,hits); + // new(lhits[fNhits++]) AliTPChit(fIshunt,gAlice->CurrentTrack(),vol,hits); + if (fHitType&&2){ + gMC->TrackMomentum(p); + Float_t momentum = TMath::Sqrt(p[0]*p[0]+p[1]*p[1]); + Float_t precision = (momentum>0.1) ? 0.002 :0.01; + fTrackHits->SetHitPrecision(precision); + } + AddHit(gAlice->CurrentTrack(), vol,hits); //MI change } diff --git a/TPC/AliTPCv3.cxx b/TPC/AliTPCv3.cxx index 0fec01b0c02..4190b85fdef 100644 --- a/TPC/AliTPCv3.cxx +++ b/TPC/AliTPCv3.cxx @@ -15,6 +15,9 @@ /* $Log$ +Revision 1.16 2000/10/02 21:28:18 fca +Removal of useless dependecies via forward declarations + Revision 1.15 2000/07/10 20:57:39 hristov Update of TPC code and macros by M.Kowalski @@ -424,38 +427,31 @@ void AliTPCv3::CreateGeometry() // Daughter volumes - // Tpc SAndwich 5 - Al + // Tpc SAndwich 5 - Tedlar dm[0]= 258.; dm[1]= 260.05; dm[2]= 251.7; - gMC->Gsvolu("TSA5","TUBE",idtmed[4],dm,3); + gMC->Gsvolu("TSA5","TUBE",idtmed[9],dm,3); - // Tpc SAndwich 6 - Tedlar + // Tpc SAndwich 6 - Kevlar dm[0] += 5.e-3; dm[1] -= 5.e-3; - gMC->Gsvolu("TSA6","TUBE",idtmed[9],dm,3); - - // Tpc SAndwich 7 - Kevlar + gMC->Gsvolu("TSA6","TUBE",idtmed[5],dm,3); - dm[0] += 5.e-3; - dm[1] -= 5.e-3; - gMC->Gsvolu("TSA7","TUBE",idtmed[5],dm,3); - - // Tpc SAndwich 8 - NOMEX + // Tpc SAndwich 7 - NOMEX dm[0] += 0.02; dm[1] -= 0.02; - gMC->Gsvolu("TSA8","TUBE",idtmed[6],dm,3); + gMC->Gsvolu("TSA7","TUBE",idtmed[6],dm,3); - // 8->7->6->5->TOFC + // 7->6->5->TOFC - gMC->Gspos("TSA8",1,"TSA7",0.,0.,0.,0,"ONLY"); gMC->Gspos("TSA7",1,"TSA6",0.,0.,0.,0,"ONLY"); gMC->Gspos("TSA6",1,"TSA5",0.,0.,0.,0,"ONLY"); @@ -893,38 +889,31 @@ void AliTPCv3::CreateGeometry() // Daughter volumes - // Tpc Sandwich 17 - Al + // Tpc Sandwich 17 - Tedlar dm[0]= 77.15; dm[1]= 79.2; dm[2]= 251.7; - gMC->Gsvolu("TS17","TUBE",idtmed[4],dm,3); + gMC->Gsvolu("TS17","TUBE",idtmed[9],dm,3); - // Tpc Sandwich 18 - Tedlar + // Tpc Sandwich 18 - Kevlar dm[0]+= 5.e-3; dm[1]-= 5.e-3; - gMC->Gsvolu("TS18","TUBE",idtmed[9],dm,3); + gMC->Gsvolu("TS18","TUBE",idtmed[5],dm,3); - // Tpc Sandwich 19 - Kevlar - dm[0]+= 5.e-3; - dm[1]-= 5.e-3; + // Tpc Sandwich 19 - NOMEX - gMC->Gsvolu("TS19","TUBE",idtmed[5],dm,3); - - // Tpc Sandwich 20 - NOMEX - - dm[0]+= 0.06; - dm[1]-= 0.06; + dm[0]+= 0.02; + dm[1]-= 0.02; - gMC->Gsvolu("TS20","TUBE",idtmed[6],dm,3); + gMC->Gsvolu("TS19","TUBE",idtmed[6],dm,3); - // 20->19->18->17 + // 19->18->17 - gMC->Gspos("TS20",1,"TS19",0.,0.,0.,0,"ONLY"); gMC->Gspos("TS19",1,"TS18",0.,0.,0.,0,"ONLY"); gMC->Gspos("TS18",1,"TS17",0.,0.,0.,0,"ONLY"); @@ -1678,8 +1667,6 @@ void AliTPCv3::CreateGeometry() } // end of function - - //_____________________________________________________________________________ void AliTPCv3::DrawDetector() @@ -1713,7 +1700,6 @@ void AliTPCv3::DrawDetector() gMC->Gsatt("TSA5","SEEN",0); gMC->Gsatt("TSA6","SEEN",0); gMC->Gsatt("TSA7","SEEN",0); - gMC->Gsatt("TSA8","SEEN",0); gMC->Gsatt("TIIN","COLO",7); gMC->Gsatt("TIIN","SEEN",1); gMC->Gsatt("TICL","SEEN",0); @@ -1731,7 +1717,6 @@ void AliTPCv3::DrawDetector() gMC->Gsatt("TS17","SEEN",0); gMC->Gsatt("TS18","SEEN",0); gMC->Gsatt("TS19","SEEN",0); - gMC->Gsatt("TS20","SEEN",0); gMC->Gsatt("TS21","SEEN",0); gMC->Gsatt("TS22","SEEN",0); gMC->Gsatt("TS23","SEEN",0); -- 2.43.0