rtang/ANASEN_analysis
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

resolution updates

Browse Source
This commit is contained in:
James Szalkie 2026-08-28 12:25:50 -04:00
parent edc7776493
commit af1ea4144a
2 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
.vscode/settings.json vendored
View File

@ -127,5 +127,6 @@
"inspect.C": "cpp" "inspect.C": "cpp"
}, },
"github-enterprise.uri": "https://fsunuc.physics.fsu.edu", "github-enterprise.uri": "https://fsunuc.physics.fsu.edu",
"C_Cpp.default.compilerPath": "/usr/bin/gcc" "C_Cpp.default.compilerPath": "/usr/bin/gcc",
"cmake.sourceDirectory": "/home/jamesszalkie/anasen/Armory/AnasenG4"
} }

6
Armory/anasenMS.cpp
View File

@ -495,7 +495,7 @@ int main(int argc, char **argv){
vertexZ = beamEntranceZ + beamPath_cm * 10.0; // cm -> mm vertexZ = beamEntranceZ + beamPath_cm * 10.0; // cm -> mm
// transverse sampling range from the beam's distance straggle at this energy // transverse sampling range from the beam's distance straggle at this energy
const double sigmaX_mm = std::max(0.0, sigmaXBeam->Eval(beamEnergy)) * 100.0; // cm -> mm const double sigmaX_mm = std::max(0.0, sigmaXBeam->Eval(beamEnergy)) * 10.0; // cm -> mm
//vertexX = 2.0 * sigmaX_mm * gRandom->Rndm() - sigmaX_mm; //vertexX = 2.0 * sigmaX_mm * gRandom->Rndm() - sigmaX_mm;
//vertexY = 2.0 * sigmaX_mm * gRandom->Rndm() - sigmaX_mm; //vertexY = 2.0 * sigmaX_mm * gRandom->Rndm() - sigmaX_mm;
double vertexRangeX = std::sqrt((vertexXRange[1] * vertexXRange[1]) + (sigmaX_mm * sigmaX_mm)); double vertexRangeX = std::sqrt((vertexXRange[1] * vertexXRange[1]) + (sigmaX_mm * sigmaX_mm));
@ -655,7 +655,7 @@ int main(int argc, char **argv){
double hitY = hitSigma.Y(); double hitY = hitSigma.Y();
double hitZ = hitSigma.Z(); double hitZ = hitSigma.Z();
originalEnergy = CalculateOriginalEnergy(hitX, hitY, hitZ, originalEnergy = CalculateOriginalEnergy(std::lround(hitX), std::lround(hitY), std::lround(hitZ),
0, 0, std::lround(vertexZ / 10) * 10, 0, 0, std::lround(vertexZ / 10) * 10,
b, "He", Esx3, b, "He", Esx3,
distance_sx3); distance_sx3);
@ -687,7 +687,7 @@ int main(int argc, char **argv){
//tree1->Fill(); //tree1->Fill();
Kinematics apkin_27Al(26.981538408,4.00260325413,1.00782503224,29.973770136,beamEnergy/26.981538408); //m3 is proton Kinematics apkin_27Al(26.981538408,4.00260325413,1.00782503224,29.973770136,beamEnergy/26.981538408); //m3 is proton
Ex_recon = apkin_27Al.getExc(originalEnergy, std::lround(thetab / 10) * 10); Ex_recon = apkin_27Al.getExc(originalEnergy, std::lround(thetab));
//EBeam_Kin_gs = apkin_27Al.getEbeam_givenQ(Esx3, 0.0, thetab); //EBeam_Kin_gs = apkin_27Al.getEbeam_givenQ(Esx3, 0.0, thetab);
//EBeam_Kin_2_2 = apkin_27Al.getEbeam_givenQ(Esx3, 2.2, thetab); //EBeam_Kin_2_2 = apkin_27Al.getEbeam_givenQ(Esx3, 2.2, thetab);
//EBeam_Kin_3_4 = apkin_27Al.getEbeam_givenQ(Esx3, 3.4, thetab); //EBeam_Kin_3_4 = apkin_27Al.getEbeam_givenQ(Esx3, 3.4, thetab);