Skip to content

Cavepipe update #14714

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
sawenzel merged 2 commits into from
Oct 3, 2025
Merged

Cavepipe update #14714

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
f9d0156
extension of barrel volume to avoid extrusion of FOCAL
amorsch Oct 2, 2025
3edc0f9
beampipes adapted to new barrel and corrections: overlaps, gaps, Alu …
amorsch Oct 2, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion Detectors/Passive/include/DetectorsPassive/PipeRun4.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -46,7 +46,7 @@ class PipeRun4 : public PassiveBase
TGeoPcon* makeMotherFromTemplate(const TGeoPcon* shape, int imin = -1, int imax = -1, float r0 = 0.,
int nz = -1);
TGeoPcon* makeInsulationFromTemplate(TGeoPcon* shape);
TGeoVolume* makeBellow(const char* ext, int nc, float rMin, float rMax, float dU, float rPlie,
TGeoVolume* makeBellow(const char* ext, int nc, float rMin, float rMax, float rPlie,
float dPlie);
TGeoVolume* makeBellowCside(const char* ext, int nc, float rMin, float rMax, float rPlie, float dPlie);

Expand Down
16 changes: 10 additions & 6 deletions Detectors/Passive/src/Cave.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -85,23 +85,27 @@ void Cave::ConstructGeometry()
shCaveTR1->DefineSection(0, -706. - 8.6, 0., 790.5);
shCaveTR1->DefineSection(1, 707. + 7.6, 0., 790.5);
TGeoTube* shCaveTR2 = new TGeoTube("shCaveTR2", 0., 150., 110.);
TGeoTube* shCaveTR3 = new TGeoTube("shCaveTR3", 0., 80., 75.);

TGeoTranslation* transCaveTR2 = new TGeoTranslation("transTR2", 0, 30., -505. - 110.);
TGeoTranslation* transCaveTR3 = new TGeoTranslation("transTR3", 0, 30., 714.6 + 75.);
transCaveTR2->RegisterYourself();
TGeoCompositeShape* shCaveTR = new TGeoCompositeShape("shCaveTR", "shCaveTR1-shCaveTR2:transTR2");
transCaveTR3->RegisterYourself();

TGeoCompositeShape* shCaveTR = new TGeoCompositeShape("shCaveTR", "shCaveTR1-shCaveTR2:transTR2+shCaveTR3:transTR3");
TGeoVolume* voBarrel = new TGeoVolume("barrel", shCaveTR, kMedAir);
cavevol->AddNode(voBarrel, 1, new TGeoTranslation(0., -30., 0.));
if (mHasRB24) { // should be not true only for alice 3
// mother volume for RB24 side (FDD, Compensator)
const Float_t kRB24CL = 2. * 597.9;
const Float_t kRB24CL = 2. * 597.9 - 150.;
auto shCaveRB24 = new TGeoPcon(0., 360., 6);
Float_t z0 = kRB24CL / 2 + 714.6;
Float_t z0 = kRB24CL / 2 + 714.6 + 150.;
shCaveRB24->DefineSection(0, -kRB24CL / 2., 0., 105.);
shCaveRB24->DefineSection(1, -z0 + 1705., 0., 105.);
shCaveRB24->DefineSection(2, -z0 + 1705., 0., 14.5);
shCaveRB24->DefineSection(3, -z0 + 1880., 0., 14.5);
shCaveRB24->DefineSection(4, -z0 + 1880., 0., 40.0);
shCaveRB24->DefineSection(5, kRB24CL / 2, 0., 40.0);
shCaveRB24->DefineSection(3, -z0 + 1878, 0., 14.5);
shCaveRB24->DefineSection(4, -z0 + 1878., 0., 40.0);
shCaveRB24->DefineSection(5, kRB24CL / 2., 0., 40.0);

TGeoVolume* caveRB24 = new TGeoVolume("caveRB24", shCaveRB24, kMedAir);
caveRB24->SetVisibility(0);
Expand Down
6 changes: 3 additions & 3 deletions Detectors/Passive/src/Pipe.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -786,13 +786,13 @@ void Pipe::ConstructGeometry()

// Copper Tube RB24/1
const Float_t kRB24CuTubeL = 381.5;
const Float_t kRB24cCuTubeL = 155.775;
const Float_t kRB24cCuTubeL = 155.775 - 150.;
const Float_t kRB24bCuTubeL = kRB24CuTubeL - kRB24cCuTubeL;
const Float_t kRB24CuTubeRi = 8.0 / 2.;
const Float_t kRB24CuTubeRo = 8.4 / 2.;
const Float_t kRB24CuTubeFRo = 7.6;
const Float_t kRB24CuTubeFL = 1.86;
const Float_t kRB24CL = 2. * 597.9;
const Float_t kRB24CL = 2. * 597.9 - 150.;
//
// introduce cut at end of barrel 714.6m
//
Expand All @@ -812,7 +812,7 @@ void Pipe::ConstructGeometry()
voRB24cCuTubeM->AddNode(voRB24cCuTube, 1, gGeoIdentity);

// Air outside tube with higher transport cuts
TGeoVolume* voRB24CuTubeA = new TGeoVolume("voRB24CuTubeA", new TGeoTube(80., 81., kRB24bCuTubeL / 2.), kMedAirHigh);
TGeoVolume* voRB24CuTubeA = new TGeoVolume("voRB24CuTubeA", new TGeoTube(79., 80., kRB24bCuTubeL / 2.), kMedAirHigh);
voRB24CuTubeA->SetVisibility(0);

// Simplified DN 100 Flange
Expand Down
117 changes: 52 additions & 65 deletions Detectors/Passive/src/PipeRun4.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -642,7 +642,7 @@ void PipeRun4::ConstructGeometry()
// Drawings from C. Gargiulo :
// \\cern.ch\dfs\Workspaces\c\cgargiul\EXPERIMENT\ALICE\ALICE_MECHANICS\ALICE_DATA_PACKAGE\IN\DETECTORS\ITS_UPGRADE\1-DESIGN\0-IF_Control_Drawing\20140207_ICD_ITS_MFT_BP
/////////////////////////////////////////////////////////////////////

TGeoVolumeAssembly* beamPipeAsideSection = new TGeoVolumeAssembly("BeamPipeAsideSection");
float kConicalBerilliumMinThickness = 0.08;
float kConicalBerilliumMaxThickness = 0.1;
float kFlangeZ = 483.75;
Expand All @@ -656,33 +656,34 @@ void PipeRun4::ConstructGeometry()
float kConicalBePipeEndOuterRadius = 3.0;

TGeoPcon* tube0 = new TGeoPcon(0., 360., 5);
tube0->DefineSection(0, kFlangeZ - kFlangeWidth / 2, kConicalBePipeEndOuterRadius - kConicalBerilliumMaxThickness, kConicalBePipeEndOuterRadius);
tube0->DefineSection(1, kConicalBerylliumEnd, kConicalBePipeEndOuterRadius - kConicalBerilliumMaxThickness, kConicalBePipeEndOuterRadius);
tube0->DefineSection(0, kFlangeZ - kFlangeWidth / 2, kConicalBePipeEndOuterRadius - kConicalBerilliumMinThickness, kConicalBePipeEndOuterRadius);
tube0->DefineSection(1, kConicalBerylliumEnd, kConicalBePipeEndOuterRadius - kConicalBerilliumMinThickness, kConicalBePipeEndOuterRadius);
tube0->DefineSection(2, kSupport1 + kSupportWidth, kPipeRadiusAtSupport1 - kConicalBerilliumMinThickness, kPipeRadiusAtSupport1);
tube0->DefineSection(3, kSupport1, kPipeRadiusAtSupport1 - kConicalBerilliumMinThickness, kPipeRadiusAtSupport1);
tube0->DefineSection(4, kBeryliumSectionZmax, kBeryliumSectionOuterRadius - kConicalBerilliumMinThickness, kBeryliumSectionOuterRadius); // need a transition to kConicalBerilliumMaxThickness

TGeoPcon* tube0vide = new TGeoPcon(0., 360., 5);
tube0vide->DefineSection(0, kFlangeZ - kFlangeWidth / 2, 0., kConicalBePipeEndOuterRadius - kConicalBerilliumMaxThickness - 0.01);
tube0vide->DefineSection(1, kConicalBerylliumEnd, 0., kConicalBePipeEndOuterRadius - kConicalBerilliumMaxThickness - 0.01);
tube0vide->DefineSection(2, kSupport1 + kSupportWidth, 0, kPipeRadiusAtSupport1 - kConicalBerilliumMinThickness - 0.01);
tube0vide->DefineSection(3, kSupport1, 0, kPipeRadiusAtSupport1 - kConicalBerilliumMinThickness - 0.01);
tube0vide->DefineSection(4, kBeryliumSectionZmax, 0., kBeryliumSectionOuterRadius - kConicalBerilliumMinThickness - 0.01);
TGeoPcon* tube0Mo = new TGeoPcon(0., 360., 5);
tube0Mo->DefineSection(0, kFlangeZ - kFlangeWidth / 2, 0., kConicalBePipeEndOuterRadius);
tube0Mo->DefineSection(1, kConicalBerylliumEnd, 0., kConicalBePipeEndOuterRadius);
tube0Mo->DefineSection(2, kSupport1 + kSupportWidth, 0, kPipeRadiusAtSupport1);
tube0Mo->DefineSection(3, kSupport1, 0, kPipeRadiusAtSupport1);
tube0Mo->DefineSection(4, kBeryliumSectionZmax, 0., kBeryliumSectionOuterRadius);

TGeoVolume* votube0 = new TGeoVolume("votube0", tube0, kMedBe);
votube0->SetLineColor(kRed);
TGeoVolume* votube0vide = new TGeoVolume("votube0vide", tube0vide, kMedVac);
votube0vide->SetLineColor(kGreen);
TGeoVolume* votube0Mo = new TGeoVolume("votube0Mo", tube0Mo, kMedVac);
votube0Mo->AddNode(votube0, 1, gGeoIdentity);
votube0Mo->SetLineColor(kGreen);

barrel->AddNode(votube0, 1, new TGeoTranslation(0., 30., 0.));
barrel->AddNode(votube0vide, 1, new TGeoTranslation(0., 30., 0.));
beamPipeAsideSection->AddNode(votube0Mo, 1, gGeoIdentity);

TGeoVolume* beampipeSupportA1 = makeSupportBar("A1", kPipeRadiusAtSupport1 + 0.01, kPipeRadiusAtSupport1 + 0.38, 20.67, 14.25);
barrel->AddNode(beampipeSupportA1, 1, new TGeoTranslation(0., 30, kSupport1 + kSupportWidth / 2.));
// already defined in IT3
// TGeoVolume* beampipeSupportA1 = makeSupportBar("A1", kPipeRadiusAtSupport1 + 0.01, kPipeRadiusAtSupport1 + 0.38, 20.67, 14.25);
// beamPipeAsideSection->AddNode(beampipeSupportA1, 1, new TGeoTranslation(0., 0., kSupport1 + kSupportWidth / 2.));

// Length is approximate
TGeoVolume* beampipeSupportA2 = makeSupportBar("A2", kConicalBePipeEndOuterRadius, kConicalBePipeEndOuterRadius + 0.38, 44, 37.5);
barrel->AddNode(beampipeSupportA2, 1, new TGeoTranslation(0., 30, kConicalBerylliumEnd + kSupportWidth / 2.));
beamPipeAsideSection->AddNode(beampipeSupportA2, 1, new TGeoTranslation(0., 0., kConicalBerylliumEnd + kSupportWidth / 2.));

TGeoPcon* Bolt1 = new TGeoPcon(0., 360, 8);
Bolt1->DefineSection(0, 0, 0, 0.5);
Expand Down Expand Up @@ -734,7 +735,7 @@ void PipeRun4::ConstructGeometry()
Bolts->AddNode(volBolt1, 7, t7);
Bolts->AddNode(volBolt1, 8, t8);

barrel->AddNode(Bolts, 1, new TGeoTranslation(0., 30., 0.));
beamPipeAsideSection->AddNode(Bolts, 1, gGeoIdentity);

TGeoTranslation* Tflange = new TGeoTranslation(0, 0, kFlangeZ);
Tflange->SetName("Tflange");
Expand All @@ -753,53 +754,40 @@ void PipeRun4::ConstructGeometry()
TGeoVolume* volflange = new TGeoVolume("voFlangeHoles", FlangeWithHoles, kMedAlBe);
volflange->SetLineWidth(2);
volflange->SetLineColor(kGray);

barrel->AddNode(volflange, 1, new TGeoTranslation(0., 30., 0.));

TGeoPcon* pipeSamell = new TGeoPcon(0., 360., 2);
pipeSamell->DefineSection(0, kFlangeZ + kFlangeWidth / 2, kConicalBePipeEndOuterRadius - kConicalBerilliumMaxThickness, kConicalBePipeEndOuterRadius);
pipeSamell->DefineSection(1, kFlangeZ + 5.13 + 0.435 + 0.4 + 0.08, kConicalBePipeEndOuterRadius - kConicalBerilliumMaxThickness, kConicalBePipeEndOuterRadius);
pipeSamell->SetName("pipeSamell");

TGeoVolume* VolpipeSmall = new TGeoVolume("voPipeSmallVac", pipeSamell, kMedAlu2219);
VolpipeSmall->SetLineWidth(2);
barrel->AddNode(VolpipeSmall, 1, new TGeoTranslation(0., 30., 0.));

TGeoPcon* pipeSmallVac = new TGeoPcon(0., 360., 2);
pipeSmallVac->DefineSection(0, kFlangeZ + kFlangeWidth / 2, 0, kConicalBePipeEndOuterRadius - kConicalBerilliumMaxThickness - 0.01);
pipeSmallVac->DefineSection(1, kFlangeZ + 5.13 + 0.435 + 0.4 + 0.08, 0, kConicalBePipeEndOuterRadius - kConicalBerilliumMaxThickness - 0.01);
TGeoVolume* vopipeSmallVac = new TGeoVolume("voPipeSmallVac", pipeSmallVac, kMedVac);
vopipeSmallVac->SetLineColor(kGreen);

barrel->AddNode(vopipeSmallVac, 1, new TGeoTranslation(0., 30., 0.));
beamPipeAsideSection->AddNode(volflange, 1, gGeoIdentity);

// -- Bellows on A side
// float plieradius = (3.72 + (2. * 7 - 2.) * 0.03) / (4. * 7); // radius of bellows "plis"
float plieradiusA = 0.2; // radius of bellow plies

// ------------------ First Bellow --------------------
// Inner: 3.0 cm, outer 3.97 cm length 8.47 cm with 10 wiggles
// check meaning of dU ; it is probably the total length, see also below
TGeoVolume* vobellows1A = makeBellow("bellows1A", 10, 3.0, 3.97, 8.47, plieradiusA, 0.03);
// Z position is rough for now.
barrel->AddNode(vobellows1A, 1, new TGeoTranslation(0., 30., kFlangeZ + 10));
TGeoVolume* vobellows1A = makeBellow("bellows1A", 10, 3.0, 3.97, plieradiusA, 0.03);
Float_t dU = (static_cast<TGeoTube*>(vobellows1A->GetShape()))->GetDZ();
beamPipeAsideSection->AddNode(vobellows1A, 1, new TGeoTranslation(0., 0., kFlangeZ + 2. * dU));
// Comments: removing 1/2 plie (see makeBellow): 0.31= 2*0.17-0.03 and 0.08: free space
Float_t pipeSmallDz = (dU - kFlangeWidth / 2.) / 2.;
TGeoTube* pipeSmall = new TGeoTube("pipeSmall", kConicalBePipeEndOuterRadius - kConicalBerilliumMaxThickness, kConicalBePipeEndOuterRadius, pipeSmallDz);
TGeoVolume* vopipeSmall = new TGeoVolume("voPipeSmall", pipeSmall, kMedAlu2219);
vopipeSmall->SetLineWidth(2);

// ------------------ Outer pipe after flange --------------------
TGeoPcon* pipeOut = new TGeoPcon(0., 360., 2);
pipeOut->DefineSection(0, kFlangeZ + 13.6 - 0.08, kConicalBePipeEndOuterRadius - kConicalBerilliumMaxThickness, kConicalBePipeEndOuterRadius);
pipeOut->DefineSection(1, 714.6, kConicalBePipeEndOuterRadius - kConicalBerilliumMaxThickness, kConicalBePipeEndOuterRadius);
TGeoTube* pipeSmallMo = new TGeoTube(0., kConicalBePipeEndOuterRadius, pipeSmallDz);
TGeoVolume* vopipeSmallMo = new TGeoVolume("voPipeSmallMo", pipeSmallMo, kMedVac);
vopipeSmallMo->SetLineColor(kGreen);
vopipeSmallMo->AddNode(vopipeSmall, 1, gGeoIdentity);

TGeoVolume* OuterPIPE = new TGeoVolume("pipeOut", pipeOut, kMedAlu2219);
barrel->AddNode(OuterPIPE, 1, new TGeoTranslation(0., 30., 0.));
beamPipeAsideSection->AddNode(vopipeSmallMo, 1, new TGeoTranslation(0., 0., kFlangeZ + kFlangeWidth / 2. + pipeSmallDz));

// The end of the barrel volume is at 714.6 cm, after that we start with RB24 volume
TGeoPcon* pipeOutVac = new TGeoPcon(0., 360., 2);
pipeOutVac->DefineSection(0, kFlangeZ + 13.6 - 0.08, 0, kConicalBePipeEndOuterRadius - kConicalBerilliumMaxThickness);
pipeOutVac->DefineSection(1, 714.6, 0., kConicalBePipeEndOuterRadius - kConicalBerilliumMaxThickness);
// ------------------ Outer pipe after flange --------------------
// The end of the barrel volume is at 864.6 cm, after that we start with RB24 volume
Float_t pipeEndZ = 864.6;
Float_t pipeOutDz = (pipeEndZ - (kFlangeZ + 3. * dU)) / 2.;
TGeoTube* pipeOut = new TGeoTube(kConicalBePipeEndOuterRadius - kConicalBerilliumMaxThickness, kConicalBePipeEndOuterRadius, pipeOutDz);
TGeoVolume* OuterPIPE = new TGeoVolume("pipeOut", pipeOut, kMedAlu2219);

TGeoVolume* OuterPIPEVac = new TGeoVolume("pipeOutVac", pipeOutVac, kMedAlu2219);
barrel->AddNode(OuterPIPEVac, 1, new TGeoTranslation(0., 30., 0.));
TGeoTube* pipeOutMo = new TGeoTube(0., kConicalBePipeEndOuterRadius, pipeOutDz);
TGeoVolume* OuterPIPEMo = new TGeoVolume("pipeOutMo", pipeOutMo, kMedVac);
OuterPIPEMo->AddNode(OuterPIPE, 1, gGeoIdentity);
beamPipeAsideSection->AddNode(OuterPIPEMo, 1, new TGeoTranslation(0., 0., pipeEndZ - pipeOutDz));
barrel->AddNode(beamPipeAsideSection, 1, new TGeoTranslation(0., 30., 0.));

//-------------------------------------------------

Expand All @@ -822,19 +810,19 @@ void PipeRun4::ConstructGeometry()

// Copper Tube RB24/1
const float kRB24CuTubeL = 381.5;
const float kRB24cCuTubeL = 155.775 + (28.375 - 18.135);
const float kRB24cCuTubeL = 155.775 - 150.;
const float kRB24bCuTubeL = kRB24CuTubeL - kRB24cCuTubeL;
const float kRB24CuTubeRi = 5.8 / 2.;
const float kRB24CuTubeRo = 6.0 / 2.;
const float kRB24CuTubeFRo = 7.6;
const float kRB24CuTubeFL = 1.86;
const float kRB24CL = 2. * 597.9;
const float kRB24CL = 2. * 597.9 - 150.;

//
// introduce cut at end of barrel 714.6m
//
// outside barrel
TGeoVolume* voRB24cCuTubeM = new TGeoVolume("voRB24cCuTubeM", new TGeoTube(0., kRB24CuTubeRi, kRB24cCuTubeL / 2.), kMedVacNFHC);
TGeoVolume* voRB24cCuTubeM = new TGeoVolume("voRB24cCuTubeM", new TGeoTube(0., kRB24CuTubeRo, kRB24cCuTubeL / 2.), kMedVacNFHC);
TGeoVolume* voRB24cCuTube = new TGeoVolume("voRB24cCuTube", new TGeoTube(kRB24CuTubeRi, kRB24CuTubeRo, kRB24cCuTubeL / 2.), kMedAlu2219);
voRB24cCuTubeM->AddNode(voRB24cCuTube, 1, gGeoIdentity);

Expand Down Expand Up @@ -876,7 +864,7 @@ void PipeRun4::ConstructGeometry()
const float kRB24B1PlieThickness = 0.015; // Plie thickness

const float kRB24B1PlieRadius =
(kRB24B1BellowUndL + (2. * kRB24B1NumberOfPlies - 2.) * kRB24B1PlieThickness) / (4. * kRB24B1NumberOfPlies);
(kRB24B1BellowUndL + 2. * kRB24B1NumberOfPlies * kRB24B1PlieThickness) / (4. * kRB24B1NumberOfPlies + 2.);

const float kRB24B1ProtTubeThickness = 0.02; // Thickness of the protection tube
const float kRB24B1ProtTubeLength = 4.2; // Length of the protection tube
Expand All @@ -892,7 +880,7 @@ void PipeRun4::ConstructGeometry()
//
// Bellow Section
TGeoVolume* voRB24B1Bellow = makeBellow("RB24B1", kRB24B1NumberOfPlies, kRB24B1BellowRi, kRB24B1BellowRo,
kRB24B1BellowUndL, kRB24B1PlieRadius, kRB24B1PlieThickness);
kRB24B1PlieRadius, kRB24B1PlieThickness);
voRB24B1Bellow->SetVisibility(0);
float newRB24B1BellowUndL = 2 * (static_cast<TGeoTube*>(voRB24B1Bellow->GetShape()))->GetDz();

Expand Down Expand Up @@ -2840,13 +2828,11 @@ TGeoPcon* PipeRun4::makeInsulationFromTemplate(TGeoPcon* shape)
return insu;
}

TGeoVolume* PipeRun4::makeBellow(const char* ext, int nc, float rMin, float rMax, float dU, float rPlie,
float dPlie)
TGeoVolume* PipeRun4::makeBellow(const char* ext, int nc, float rMin, float rMax, float rPlie, float dPlie)
{
// nc Number of convolution
// rMin Inner radius of the bellow
// rMax Outer radius of the bellow
// dU Undulation length
// rPlie Plie radius
// dPlie Plie thickness
auto& matmgr = o2::base::MaterialManager::Instance();
Expand Down Expand Up @@ -2896,17 +2882,18 @@ TGeoVolume* PipeRun4::makeBellow(const char* ext, int nc, float rMin, float rMax
asWiggle->AddNode(voWiggleL, 1, new TGeoTranslation(0., 0., z0));
asWiggle->GetShape()->ComputeBBox(); // enforce recomputing of BBox
//
float zBellowTot = nc * (static_cast<TGeoBBox*>(asWiggle->GetShape()))->GetDZ();
TGeoVolume* voBellow = new TGeoVolume(fmt::format("{:s}BellowUS", ext).c_str(), new TGeoTube(rMin, rMax, zBellowTot), kMedVac);
float zBellowTot = nc * (2. * (static_cast<TGeoBBox*>(asWiggle->GetShape()))->GetDZ() - dPlie) + 2. * rPlie;
TGeoVolume* voBellow = new TGeoVolume(fmt::format("{:s}BellowUS", ext).c_str(), new TGeoTube(rMin, rMax, zBellowTot / 2.), kMedVac);
// Positioning of the volumes
z0 = -dU / 2. + rPlie;
z0 = -zBellowTot / 2. + rPlie;
voBellow->AddNode(voWiggleL, 2, new TGeoTranslation(0., 0., z0));
z0 += rPlie;
float zsh = 4. * rPlie - 2. * dPlie;
for (int iw = 0; iw < nc; iw++) {
float zpos = z0 + iw * zsh;
voBellow->AddNode(asWiggle, iw + 1, new TGeoTranslation(0., 0., zpos - dPlie));
}

return voBellow;
}

Expand Down