Skip to content

Commit 526b5c9

Browse files
committed
- put back initialization of the e3 recoil TLepton in the case of atomic
(coherent) pair conversion. [rtj] - adjust the Pcut thresholds to decrease the rate of clipping in the keep/discard Monte Carlo algorithm, at the expense of efficiency. [rtj]
1 parent 277a134 commit 526b5c9

File tree

1 file changed

+3
-2
lines changed

1 file changed

+3
-2
lines changed

src/GlueXBeamConversionProcess.cc

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -86,8 +86,8 @@ GlueXBeamConversionProcess::GlueXBeamConversionProcess(const G4String &name,
8686
fPairsGeneration = new PairConversionGeneration();
8787
#endif
8888

89-
fPaircohPDF.Pcut = 30;
90-
fTripletPDF.Pcut = 10;
89+
fPaircohPDF.Pcut = 60;
90+
fTripletPDF.Pcut = 15;
9191
}
9292

9393
GlueXBeamConversionProcess::GlueXBeamConversionProcess(
@@ -590,6 +590,7 @@ void GlueXBeamConversionProcess::GenerateBeamPairConversion(const G4Step &step)
590590
// returned as d(sigma)/(dE+ dphi+ d^3qR)
591591
p1.SetMom(q1.Rotate(rockaxis, rockangle));
592592
e2.SetMom(q2.Rotate(rockaxis, rockangle));
593+
e3.SetMom(TThreeVectorReal(0,0,0));
593594
LDouble_t diffXS = fPairsGeneration->DiffXS_pair(gIn, p1, e2);
594595

595596
// Use keep/discard algorithm

0 commit comments

Comments
 (0)