Skip to content

Bugfix n roo tracker truth variables #2

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

Closed
nickwp wants to merge 5 commits into from
Closed

Bugfix n roo tracker truth variables #2

Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
00dc276
Initial commit of bugfix to correctly copy neutrino and target info when
mscott201 Jul 22, 2026
ff7ee63
Merge remote-tracking branch 'origin/develop' into bugfix_nRooTracker…
mscott201 Jul 22, 2026
c5c3d01
Commit changes to allow compilation
mscott201 Jul 22, 2026
a6e87d4
Fix bug in ooTracker target energy storage - it should be total energ…
mscott201 Jul 22, 2026
d3e7677
Remove whitespace changes - someone should fix all the incorrect whit…
mscott201 Jul 24, 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
62 changes: 40 additions & 22 deletions src/WCSimPrimaryGeneratorAction.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -616,19 +616,14 @@ void WCSimPrimaryGeneratorAction::GeneratePrimaries(G4Event* anEvent)
G4ThreeVector dir;
double kin_energy;

// Add interaction mode details
TObjString* code = fTmpRootrackerVtx->EvtCode;
mode[0] = code->String().Atoi();

// Now simulate the outgoing particles
for (int i = 0; i < fTmpRootrackerVtx->StdHepN; i++){

// Skip the initial neutrino, target nucleus and target nucleon
int abspdg = abs(fTmpRootrackerVtx->StdHepPdgTemp[i]);
if( i < 3){
if(abspdg > 100000 || abspdg == 2112 || abspdg == 2212){
continue;
}
}
if(abspdg == 12 || abspdg == 14 || abspdg == 16) {
continue;
}
int pdg = fTmpRootrackerVtx->StdHepPdgTemp[i];

xDir=fTmpRootrackerVtx->StdHepP4[i][0]*GeV;
yDir=fTmpRootrackerVtx->StdHepP4[i][1]*GeV;
Expand All @@ -646,19 +641,42 @@ void WCSimPrimaryGeneratorAction::GeneratePrimaries(G4Event* anEvent)

dir = dir*(1./momentum);

particleGun->SetParticleDefinition(particleTable->FindParticle(fTmpRootrackerVtx->StdHepPdgTemp[i]));

kin_energy = fabs(fTmpRootrackerVtx->StdHepP4[i][3])*GeV - particleGun->GetParticleDefinition()->GetPDGMass();
kin_energy = 0;

particleGun->SetParticleEnergy(kin_energy);
particleGun->SetParticlePosition(vtx);
particleGun->SetParticleMomentumDirection(dir);
// Will want to include some beam time structure at some point, but not needed at the moment since we only simulate 1 interaction per events
// particleGun->SetParticleTime(time);
particleGun->GeneratePrimaryVertex(anEvent); //Place vertex in stack
G4cout << "ROOTRACKER: created particle gun with KinEnergy/PDG/Position/Direction: "
<< kin_energy << "/" << fTmpRootrackerVtx->StdHepPdgTemp[i]
<< "/" << vtx << "/" << dir << " i = " << i << G4endl;
// Tag the neutrino and target particles. These aren't simulated in the detector but we need their kinematic info
if(abs(pdg) == 12 || abs(pdg) == 14 || abs(pdg) == 16) {
kin_energy = fabs(fTmpRootrackerVtx->StdHepP4[i][3])*GeV;
SetVtx(vtx);
SetBeamEnergy(kin_energy);
SetBeamDir(dir);
SetBeamPDG(pdg);
}
else if( i < 3){
// These are target nuclei / nucleons. Nuclei have zero kinetic energy
if(abs(pdg) > 100000){
kin_energy = 0;
}
// Now store target nucleon info
else if(abs(pdg) == 2112 || abs(pdg) == 2212){
targetpdgs[0] = pdg;
targetenergies[0] = fabs(fTmpRootrackerVtx->StdHepP4[i][3])*GeV;
targetdirs[0] = dir;
}
}
else{
particleGun->SetParticleDefinition(particleTable->FindParticle(pdg));
kin_energy = fabs(fTmpRootrackerVtx->StdHepP4[i][3])*GeV - particleGun->GetParticleDefinition()->GetPDGMass();

particleGun->SetParticleEnergy(kin_energy);
Comment on lines +667 to +670
particleGun->SetParticlePosition(vtx);
particleGun->SetParticleMomentumDirection(dir);
// Will want to include some beam time structure at some point, but not needed at the moment since we only simulate 1 interaction per events
// particleGun->SetParticleTime(time);
particleGun->GeneratePrimaryVertex(anEvent); //Place vertex in stack
G4cout << "ROOTRACKER: created particle gun with KinEnergy/PDG/Position/Direction: "
<< kin_energy << "/" << fTmpRootrackerVtx->StdHepPdgTemp[i]
<< "/" << vtx << "/" << dir << " i = " << i << G4endl;
}
}//loop over RooTracker particles
}//useRooTracker

Expand Down
Loading