Skip to content

Fix IPA EndRing parameters and add debugging for wire support #1915

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
oksuzian merged 2 commits into from
Aug 4, 2026
Merged

Fix IPA EndRing parameters and add debugging for wire support #1915

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
4d270bb
Made Changes to the IPA so that we get IPA EndRing output. Corrected …
Aug 3, 2026
f857837
Adressed PR changes: accidentally
Aug 3, 2026
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 Mu2eG4/geom/geom_run1.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -63,7 +63,7 @@ double mu2e.detectorSystemZ0 = 10171.; // mm G4BL: (17730-7292=9801 mm)
#include "Offline/Mu2eG4/geom/STM_v09.txt" // (muon) stopping target monitor

// Proton Absorber
#include "Offline/Mu2eG4/geom/protonAbsorber_cylindrical_v04.txt"
#include "Offline/Mu2eG4/geom/protonAbsorber_cylindrical_v05.txt"
#include "Offline/Mu2eG4/geom/degrader_v02.txt" // pion degrader. Off by default

#include "Offline/Mu2eG4/geom/ProductionTarget_Hayman_v2_2.txt"
Expand Down
2 changes: 1 addition & 1 deletion Mu2eG4/geom/geom_run1_a_stickman.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -63,7 +63,7 @@ double mu2e.detectorSystemZ0 = 10171.; // mm G4BL: (17730-7292=9801 mm)
#include "Offline/Mu2eG4/geom/STM_v09.txt" // (muon) stopping target monitor

// Proton Absorber
#include "Offline/Mu2eG4/geom/protonAbsorber_cylindrical_v04.txt"
#include "Offline/Mu2eG4/geom/protonAbsorber_cylindrical_v05.txt"
#include "Offline/Mu2eG4/geom/degrader_v02.txt" // pion degrader. Off by default

#include "Offline/Mu2eG4/geom/ProductionTarget_Stickman_v1_0.txt"
Expand Down
12 changes: 12 additions & 0 deletions Mu2eG4/geom/protonAbsorber_cylindrical_v05.txt
View file
Open in desktop
Comment thread
AndrewEdmonds11 marked this conversation as resolved.
Original file line number Diff line number Diff line change
@@ -0,0 +1,12 @@
//
// Conical proton absorber geometry file.
// Correced Dimensions and endRing Material
// See Doc-57487-v3
//

#include "Offline/Mu2eG4/geom/protonAbsorber_cylindrical_v04.txt"

//Corrected Dimensions - overwriting previous entries
double protonabsorber.ipa.endRingHalfLength = 3.175; // 1/4" long in the z-direction
double protonabsorber.ipa.endRingRadialLength = 19.05; // 3/4" long in the r-direction
string protonabsorber.ipa.endRingMaterial = "IPAPolystyrene"; //Changed Material
11 changes: 10 additions & 1 deletion Mu2eG4/src/ConstructMaterials.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
//
//
// Construct materials requested by the run-time configuration system.
//

Expand Down Expand Up @@ -1334,6 +1334,15 @@ namespace mu2e {
CorrugatedPolypropylene->AddElement( getElementOrThrow("H"), 6);
}

mat = uniqueMaterialOrThrow("IPAPolystyrene");
{
// We used an explicity density here
// Reflects the corrected density for End Ring material
G4double IPAEffectiveDensity = 0.426*CLHEP::g/CLHEP::cm3; //explicit density refer to Doc-57487-v3
G4Material* IPAPolystyrene = new G4Material( mat.name, IPAEffectiveDensity, 2);
IPAPolystyrene->AddElement( getElementOrThrow("C"), 8);
IPAPolystyrene->AddElement( getElementOrThrow("H"), 8);
}

//G10-FR4 used for printed board of the I-Tracker
// G10 http://personalpages.to.infn.it/~tosello/EngMeet/ITSmat/SDD/SDD_G10FR4.html
Expand Down
3 changes: 3 additions & 0 deletions Mu2eG4/src/constructProtonAbsorber.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1614,6 +1614,9 @@ namespace mu2e {
ostringstream wirename ; wirename << "IPAsupport_set" << iS << "_wire" << ++iW ;

const double rStartOfWire = pabs1rOut0+(supportWire.originInMu2e().z()-zstartOfIPA)/pabs1len*(pabs1rOut1-pabs1rOut0);
if ( verbosityLevel > 0){
std::cout << "Set " << iSet << ", Wire " << iWire << ": rStartOfWire = " << rStartOfWire << "\n";
}
const double wirePhi = iWire * 360.*CLHEP::deg / ipaSup->nWiresPerSet() + wire_angle_offset ;
CLHEP::Hep3Vector additionalOffset ( (supportWire.halfLength()+0.005+rStartOfWire) *
std::cos(wirePhi),
Expand Down