56 std::vector<G4Material*>::const_iterator matite;
57 for( matite = materialTable->begin(); matite != materialTable->end(); matite++ ) {
115 if(
verboseLevel >= 2 )
G4cout <<
" G4LEPTSElasticModel::SampleSecondaries( P1Dir " << P1Dir <<
" P0Dir " << P0Dir <<
" ang " << ang <<
G4endl;
123 G4double P1KinEn = P0KinEn - Energylost;
126 if(
verboseLevel >= 2 )
G4cout <<
" G4LEPTSElasticModel::SampleSecondaries( SetProposedKineticEnergy " << P1KinEn <<
" " << P0KinEn <<
" - " << Energylost <<
G4endl;
144 G4double W = ( (E+MP)*si*si + MT - co*std::sqrt(MT*MT-MP*MP*si*si) ) * E*(E+2*MP)
145 / ( std::pow((E+MP+MT),2) - E*co*co*(E+2*MP) );