48 G4cout<<
"Regular X-ray TR radiator EM process is called"<<
G4endl ;
79 sigma = 0.5*(aMa + bMb);
98 if (cofMin > kMin) kMin++;
117 for( k = kMin; k <=
kMax; k++ )
120 result = (k - cof1)*(k - cof1)*(k + cof2)*(k + cof2);
122 if( k == kMin && kMin ==
G4int(cofMin) )
124 sum += 0.5*std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/
result;
128 sum += std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/
result;
130 theta2k = std::sqrt(theta2*std::abs(k-cofMin));
136 G4cout<<k<<
" "<<theta2k<<
" "<<std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/result
140 result = 2*( cof1 + cof2 )*( cof1 + cof2 )*sum/
energy;
193 G4complex F1 = (1.0 - Ha)*(1.0 - Hb )/(1.0 - H)
196 G4complex F2 = (1.0-Ha)*(1.0-Ha)*Hb/(1.0-H)/(1.0-H)
201 result = 2.0*std::real(R);
G4double SpectralXTRdEdx(G4double energy) override
std::vector< ExP01TrackerHit * > a
static constexpr double hbarc
G4double GetStackFactor(G4double energy, G4double gamma, G4double varAngle) override
std::complex< G4double > G4complex
G4double GetGasFormationZone(G4double, G4double, G4double)
G4double GetPlateLinearPhotoAbs(G4double)
G4double GetPlateFormationZone(G4double, G4double, G4double)
G4double G4ParticleHPJENDLHEData::G4double result
G4GLOB_DLL std::ostream G4cout
G4RegularXTRadiator(G4LogicalVolume *anEnvelope, G4Material *, G4Material *, G4double, G4double, G4int, const G4String &processName="XTRegularRadiator")
static constexpr double pi
G4double GetGasLinearPhotoAbs(G4double)
G4complex OneInterfaceXTRdEdx(G4double energy, G4double gamma, G4double varAngle)