Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // >> 26 // >> 27 // $Id: G4PolarizedCompton.cc,v 1.3 2006/11/17 14:14:20 vnivanch Exp $ >> 28 // GEANT4 tag $Name: geant4-08-03-patch-01 $ >> 29 // >> 30 // 26 // File name: G4PolarizedCompton 31 // File name: G4PolarizedCompton 27 // 32 // 28 // Author: Andreas Schaelicke 33 // Author: Andreas Schaelicke 29 // based on code by Michel Mair 34 // based on code by Michel Maire / Vladimir IVANTCHENKO 30 // << 31 // Class description 35 // Class description 32 // modified version respecting media and bea << 36 // 33 // using the stokes formalism << 37 // modified version respecting media and beam polarization >> 38 // using the stokes formalism >> 39 // >> 40 // Creation date: 01.05.2005 >> 41 // >> 42 // Modifications: >> 43 // >> 44 // 01-01-05, include polarization description (A.Stahl) >> 45 // 01-01-05, create asymmetry table and determine interactionlength (A.Stahl) >> 46 // 01-05-05, update handling of media polarization (A.Schalicke) >> 47 // 01-05-05, update polarized differential cross section (A.Schalicke) >> 48 // 20-05-05, added polarization transfer (A.Schalicke) >> 49 // 10-06-05, transformation between different reference frames (A.Schalicke) >> 50 // 17-10-05, correct reference frame dependence in GetMeanFreePath (A.Schalicke) >> 51 // 26-07-06, cross section recalculated (P.Starovoitov) >> 52 // 09-08-06, make it work under current geant4 release (A.Schalicke) >> 53 // ----------------------------------------------------------------------------- 34 54 35 #include "G4PolarizedCompton.hh" << 36 55 >> 56 #include "G4PolarizedCompton.hh" 37 #include "G4Electron.hh" 57 #include "G4Electron.hh" 38 #include "G4EmParameters.hh" << 58 39 #include "G4KleinNishinaCompton.hh" << 59 #include "G4StokesVector.hh" 40 #include "G4PhysicsTableHelper.hh" << 41 #include "G4PolarizationManager.hh" 60 #include "G4PolarizationManager.hh" 42 #include "G4PolarizedComptonModel.hh" 61 #include "G4PolarizedComptonModel.hh" 43 #include "G4ProductionCutsTable.hh" 62 #include "G4ProductionCutsTable.hh" 44 #include "G4StokesVector.hh" << 63 #include "G4PhysicsTableHelper.hh" 45 #include "G4SystemOfUnits.hh" << 64 #include "G4KleinNishinaCompton.hh" >> 65 #include "G4PolarizedComptonModel.hh" 46 66 47 //....oooOO0OOooo........oooOO0OOooo........oo 67 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 48 G4PhysicsTable* G4PolarizedCompton::theAsymmet << 49 68 50 G4PolarizedCompton::G4PolarizedCompton(const G 69 G4PolarizedCompton::G4PolarizedCompton(const G4String& processName, 51 G4Proce << 70 G4ProcessType type): 52 : G4VEmProcess(processName, type) << 71 G4VEmProcess (processName, type), 53 , fType(10) << 72 buildAsymmetryTable(true), 54 , fBuildAsymmetryTable(true) << 73 useAsymmetryTable(true), 55 , fUseAsymmetryTable(true) << 74 isInitialised(false), 56 , fIsInitialised(false) << 75 selectedModel(0), 57 { << 76 mType(10), 58 SetStartFromNullFlag(true); << 77 theAsymmetryTable(NULL) 59 SetBuildTableFlag(true); << 78 { 60 SetSecondaryParticle(G4Electron::Electron()) << 79 SetLambdaBinning(90); 61 SetProcessSubType(fComptonScattering); << 80 SetMinKinEnergy(0.1*keV); 62 SetMinKinEnergyPrim(1. * MeV); << 81 SetMaxKinEnergy(100.0*GeV); 63 SetSplineFlag(true); << 64 fEmModel = nullptr; << 65 } 82 } 66 83 67 //....oooOO0OOooo........oooOO0OOooo........oo 84 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 68 G4PolarizedCompton::~G4PolarizedCompton() { Cl << 85 69 << 86 G4PolarizedCompton::~G4PolarizedCompton() 70 //....oooOO0OOooo........oooOO0OOooo........oo << 71 void G4PolarizedCompton::ProcessDescription(st << 72 { 87 { 73 out << "Polarized model for Compton scatteri << 88 if (theAsymmetryTable) { 74 << 75 G4VEmProcess::ProcessDescription(out); << 76 } << 77 << 78 //....oooOO0OOooo........oooOO0OOooo........oo << 79 void G4PolarizedCompton::CleanTable() << 80 { << 81 if(theAsymmetryTable) << 82 { << 83 theAsymmetryTable->clearAndDestroy(); 89 theAsymmetryTable->clearAndDestroy(); 84 delete theAsymmetryTable; 90 delete theAsymmetryTable; 85 theAsymmetryTable = nullptr; << 86 } 91 } 87 } 92 } 88 93 89 //....oooOO0OOooo........oooOO0OOooo........oo 94 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 90 G4bool G4PolarizedCompton::IsApplicable(const << 91 { << 92 return (&p == G4Gamma::Gamma()); << 93 } << 94 95 95 //....oooOO0OOooo........oooOO0OOooo........oo << 96 void G4PolarizedCompton::InitialiseProcess(con 96 void G4PolarizedCompton::InitialiseProcess(const G4ParticleDefinition*) 97 { 97 { 98 if(!fIsInitialised) << 98 if(!isInitialised) { 99 { << 99 isInitialised = true; 100 fIsInitialised = true; << 100 SetBuildTableFlag(true); 101 if(0 == fType) << 101 SetSecondaryParticle(G4Electron::Electron()); 102 { << 102 G4double emin = MinKinEnergy(); 103 if(nullptr == EmModel(0)) << 103 G4double emax = MaxKinEnergy(); 104 { << 104 emModel = new G4PolarizedComptonModel(); 105 SetEmModel(new G4KleinNishinaCompton() << 105 if(0 == mType) selectedModel = new G4KleinNishinaCompton(); 106 } << 106 else if(10 == mType) selectedModel = emModel; 107 } << 107 selectedModel->SetLowEnergyLimit(emin); 108 else << 108 selectedModel->SetHighEnergyLimit(emax); 109 { << 109 AddEmModel(1, selectedModel); 110 fEmModel = new G4PolarizedComptonModel() << 110 } 111 SetEmModel(fEmModel); << 112 } << 113 G4EmParameters* param = G4EmParameters::In << 114 EmModel(0)->SetLowEnergyLimit(param->MinKi << 115 EmModel(0)->SetHighEnergyLimit(param->MaxK << 116 AddEmModel(1, EmModel(0)); << 117 } << 118 } 111 } 119 112 120 //....oooOO0OOooo........oooOO0OOooo........oo 113 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 121 void G4PolarizedCompton::SetModel(const G4Stri << 114 >> 115 void G4PolarizedCompton::PrintInfo() 122 { 116 { 123 if(ss == "Klein-Nishina") << 117 G4cout << " Total cross sections has a good parametrisation" 124 { << 118 << " from 10 KeV to (100/Z) GeV" 125 fType = 0; << 119 << "\n Sampling according " << selectedModel->GetName() << " model" 126 } << 120 << G4endl; 127 if(ss == "Polarized-Compton") << 121 } 128 { << 122 129 fType = 10; << 123 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 130 } << 124 >> 125 void G4PolarizedCompton::SetModel(const G4String& s) >> 126 { >> 127 if(s == "Klein-Nishina") mType = 0; >> 128 if(s == "Polarized-Compton") mType = 10; 131 } 129 } 132 130 >> 131 >> 132 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 133 //....oooOO0OOooo........oooOO0OOooo........oo 133 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 134 G4double G4PolarizedCompton::GetMeanFreePath(c << 134 135 G << 135 136 G << 136 >> 137 G4double G4PolarizedCompton::GetMeanFreePath( >> 138 const G4Track& aTrack, >> 139 G4double previousStepSize, >> 140 G4ForceCondition* condition) 137 { 141 { 138 // *** get unploarised mean free path from l 142 // *** get unploarised mean free path from lambda table *** 139 G4double mfp = << 143 G4double mfp = G4VEmProcess::GetMeanFreePath(aTrack, previousStepSize, condition); 140 G4VEmProcess::GetMeanFreePath(aTrack, prev << 141 144 142 if(theAsymmetryTable && fUseAsymmetryTable & << 145 143 { << 146 if (theAsymmetryTable && useAsymmetryTable) { 144 mfp *= ComputeSaturationFactor(aTrack); << 147 // *** get asymmetry, if target is polarized *** 145 } << 148 const G4DynamicParticle* aDynamicGamma = aTrack.GetDynamicParticle(); 146 if(verboseLevel >= 2) << 149 const G4double GammaEnergy = aDynamicGamma->GetKineticEnergy(); 147 { << 150 const G4StokesVector GammaPolarization = aTrack.GetPolarization(); 148 G4cout << "G4PolarizedCompton::MeanFreePat << 151 const G4ParticleMomentum GammaDirection0 = aDynamicGamma->GetMomentumDirection(); 149 << G4endl; << 152 150 } << 153 G4Material* aMaterial = aTrack.GetMaterial(); 151 return mfp; << 154 G4VPhysicalVolume* aPVolume = aTrack.GetVolume(); >> 155 G4LogicalVolume* aLVolume = aPVolume->GetLogicalVolume(); >> 156 >> 157 // G4Material* bMaterial = aLVolume->GetMaterial(); >> 158 G4PolarizationManager * polarizationManger = G4PolarizationManager::GetInstance(); >> 159 >> 160 const G4bool VolumeIsPolarized = polarizationManger->IsPolarized(aLVolume); >> 161 G4StokesVector ElectronPolarization = polarizationManger->GetVolumePolarization(aLVolume); >> 162 >> 163 if (!VolumeIsPolarized || mfp == DBL_MAX) return mfp; >> 164 >> 165 if (verboseLevel>=2) { >> 166 >> 167 G4cout << " Mom " << GammaDirection0 << G4endl; >> 168 G4cout << " Polarization " << GammaPolarization << G4endl; >> 169 G4cout << " MaterialPol. " << ElectronPolarization << G4endl; >> 170 G4cout << " Phys. Volume " << aPVolume->GetName() << G4endl; >> 171 G4cout << " Log. Volume " << aLVolume->GetName() << G4endl; >> 172 G4cout << " Material " << aMaterial << G4endl; >> 173 } >> 174 >> 175 G4bool isOutRange; >> 176 >> 177 >> 178 G4int midx= CurrentMaterialCutsCoupleIndex(); >> 179 G4PhysicsVector * aVector=(*theAsymmetryTable)(midx); >> 180 >> 181 G4double asymmetry=0; >> 182 if (aVector) >> 183 asymmetry = aVector->GetValue(GammaEnergy, isOutRange); >> 184 else { >> 185 G4cout << " MaterialIndex " << midx << " is out of range \n"; >> 186 } >> 187 if (isOutRange) { >> 188 G4cout << " MaterialIndex " << midx >> 189 << " GammaEnergy=" << GammaEnergy << " is out of range \n"; >> 190 asymmetry=0; >> 191 } >> 192 >> 193 // we have to determine angle between particle motion >> 194 // and target polarisation here >> 195 // circ pol * Vec(ElectronPol)*Vec(PhotonMomentum) >> 196 // both vectors in global reference frame >> 197 >> 198 G4double pol=ElectronPolarization*GammaDirection0; >> 199 >> 200 G4double polProduct = GammaPolarization.p3() * pol; >> 201 mfp *= 1. / ( 1. + polProduct * asymmetry ); >> 202 >> 203 if (verboseLevel>=2) { >> 204 G4cout << " MeanFreePath: " << mfp / mm << " mm " << G4endl; >> 205 G4cout << " Asymmetry: " << asymmetry << G4endl; >> 206 G4cout << " PolProduct: " << polProduct << G4endl; >> 207 } >> 208 } >> 209 >> 210 return mfp; 152 } 211 } 153 212 154 //....oooOO0OOooo........oooOO0OOooo........oo 213 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 155 G4double G4PolarizedCompton::PostStepGetPhysic << 214 156 const G4Track& aTrack, G4double previousStep << 215 void G4PolarizedCompton::PreparePhysicsTable(const G4ParticleDefinition& part) 157 { 216 { 158 // save previous values << 217 G4VEmProcess::PreparePhysicsTable(part); 159 G4double nLength = theNumberOfInteractionLen << 218 if(buildAsymmetryTable) 160 G4double iLength = currentInteractionLength; << 219 theAsymmetryTable = G4PhysicsTableHelper::PreparePhysicsTable(theAsymmetryTable); 161 << 162 // *** compute unpolarized step limit *** << 163 // this changes theNumberOfInteractionLength << 164 G4double x = G4VEmProcess::PostStepGetPhysic << 165 aTrack, previousStepSize, condition); << 166 G4double x0 = x; << 167 G4double satFact = 1.0; << 168 << 169 // *** add corrections on polarisation *** << 170 if(theAsymmetryTable && fUseAsymmetryTable & << 171 { << 172 satFact = ComputeSaturationFact << 173 G4double curLength = currentInteractionLen << 174 G4double prvLength = iLength * satFact; << 175 if(nLength > 0.0) << 176 { << 177 theNumberOfInteractionLengthLeft = << 178 std::max(nLength - previousStepSize / << 179 } << 180 x = theNumberOfInteractionLengthLeft * cur << 181 } << 182 if(verboseLevel >= 2) << 183 { << 184 G4cout << "G4PolarizedCompton::PostStepGPI << 185 << x / mm << " mm;" << G4endl << 186 << " unpolarized valu << 187 << x0 / mm << " mm." << G4endl; << 188 } << 189 return x; << 190 } 220 } 191 221 192 //....oooOO0OOooo........oooOO0OOooo........oo 222 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 193 G4double G4PolarizedCompton::ComputeSaturation << 194 { << 195 G4double factor = 1.0; << 196 << 197 // *** get asymmetry, if target is polarized << 198 const G4DynamicParticle* aDynamicGamma = aTr << 199 const G4double GammaEnergy = aDy << 200 const G4StokesVector GammaPolarization = << 201 G4StokesVector(aTrack.GetPolarization()); << 202 const G4ParticleMomentum GammaDirection0 = << 203 aDynamicGamma->GetMomentumDirection(); << 204 << 205 const G4Material* aMaterial = aTrack.GetMate << 206 G4VPhysicalVolume* aPVolume = aTrack.GetVolu << 207 G4LogicalVolume* aLVolume = aPVolume->GetL << 208 << 209 G4PolarizationManager* polarizationManager = << 210 G4PolarizationManager::GetInstance(); << 211 << 212 const G4bool VolumeIsPolarized = polarizatio << 213 G4StokesVector ElectronPolarization = << 214 polarizationManager->GetVolumePolarization << 215 << 216 if(VolumeIsPolarized) << 217 { << 218 if(verboseLevel >= 2) << 219 { << 220 G4cout << "G4PolarizedCompton::ComputeSa << 221 G4cout << " Mom " << GammaDirection0 << << 222 G4cout << " Polarization " << GammaPolar << 223 G4cout << " MaterialPol. " << ElectronPo << 224 G4cout << " Phys. Volume " << aPVolume-> << 225 G4cout << " Log. Volume " << aLVolume-> << 226 G4cout << " Material " << aMaterial << 227 } << 228 223 229 std::size_t midx = CurrentMa << 230 const G4PhysicsVector* aVector = nullptr; << 231 if(midx < theAsymmetryTable->size()) << 232 { << 233 aVector = (*theAsymmetryTable)(midx); << 234 } << 235 if(aVector) << 236 { << 237 G4double asymmetry = aVector->Value(Gamm << 238 << 239 // we have to determine angle between p << 240 // and target polarisation here << 241 // circ pol * Vec(ElectronPol)*Vec( << 242 // both vectors in global reference fra << 243 << 244 G4double pol = ElectronPolarizati << 245 G4double polProduct = GammaPolarization. << 246 factor /= (1. + polProduct * asymmetry); << 247 if(verboseLevel >= 2) << 248 { << 249 G4cout << " Asymmetry: " << asymme << 250 G4cout << " PolProduct: " << polPro << 251 G4cout << " Factor: " << factor << 252 } << 253 } << 254 else << 255 { << 256 G4ExceptionDescription ed; << 257 ed << "Problem with asymmetry table: mat << 258 << " is out of range or the table is << 259 G4Exception("G4PolarizedComptonModel::Co << 260 JustWarning, ed, ""); << 261 } << 262 } << 263 return factor; << 264 } << 265 224 266 //....oooOO0OOooo........oooOO0OOooo........oo << 267 void G4PolarizedCompton::BuildPhysicsTable(con 225 void G4PolarizedCompton::BuildPhysicsTable(const G4ParticleDefinition& part) 268 { 226 { 269 // *** build (unpolarized) cross section tab 227 // *** build (unpolarized) cross section tables (Lambda) 270 G4VEmProcess::BuildPhysicsTable(part); 228 G4VEmProcess::BuildPhysicsTable(part); 271 if(fBuildAsymmetryTable && fEmModel) << 229 if(buildAsymmetryTable) 272 { << 230 BuildAsymmetryTable(part); 273 G4bool isMaster = true; << 274 const G4PolarizedCompton* masterProcess = << 275 static_cast<const G4PolarizedCompton*>(G << 276 if(masterProcess && masterProcess != this) << 277 { << 278 isMaster = false; << 279 } << 280 if(isMaster) << 281 { << 282 BuildAsymmetryTable(part); << 283 } << 284 } << 285 } 231 } 286 232 287 //....oooOO0OOooo........oooOO0OOooo........oo 233 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 234 >> 235 288 void G4PolarizedCompton::BuildAsymmetryTable(c 236 void G4PolarizedCompton::BuildAsymmetryTable(const G4ParticleDefinition& part) 289 { 237 { 290 // cleanup old, initialise new table << 291 CleanTable(); << 292 theAsymmetryTable = << 293 G4PhysicsTableHelper::PreparePhysicsTable( << 294 << 295 // Access to materials 238 // Access to materials 296 const G4ProductionCutsTable* theCoupleTable << 239 const G4ProductionCutsTable* theCoupleTable= 297 G4ProductionCutsTable::GetProductionCutsTa << 240 G4ProductionCutsTable::GetProductionCutsTable(); 298 G4int numOfCouples = (G4int)theCoupleTable-> << 241 size_t numOfCouples = theCoupleTable->GetTableSize(); 299 if(!theAsymmetryTable) << 242 for(size_t i=0; i<numOfCouples; ++i) { 300 { << 243 if (!theAsymmetryTable) break; 301 return; << 244 if (theAsymmetryTable->GetFlag(i)) { 302 } << 245 303 G4int nbins = LambdaBinning( << 304 G4double emin = MinKinEnergy() << 305 G4double emax = MaxKinEnergy() << 306 G4PhysicsLogVector* aVector = nullptr; << 307 G4PhysicsLogVector* bVector = nullptr; << 308 << 309 for(G4int i = 0; i < numOfCouples; ++i) << 310 { << 311 if(theAsymmetryTable->GetFlag(i)) << 312 { << 313 // create physics vector and fill it 246 // create physics vector and fill it 314 const G4MaterialCutsCouple* couple = << 247 const G4MaterialCutsCouple* couple = theCoupleTable->GetMaterialCutsCouple(i); 315 theCoupleTable->GetMaterialCutsCouple( << 316 // use same parameters as for lambda 248 // use same parameters as for lambda 317 if(!aVector) << 249 G4PhysicsVector* aVector = LambdaPhysicsVector(couple); 318 { << 250 // modelManager->FillLambdaVector(aVector, couple, startFromNull); 319 aVector = new G4PhysicsLogVector(emin, << 320 bVector = aVector; << 321 } << 322 else << 323 { << 324 bVector = new G4PhysicsLogVector(*aVec << 325 } << 326 251 327 for(G4int j = 0; j <= nbins; ++j) << 252 for (G4int j = 0 ; j < LambdaBinning() ; ++j ) { 328 { << 253 G4double lowEdgeEnergy = aVector->GetLowEdgeEnergy(j); 329 G4double energy = bVector->Energy(j); << 254 G4double tasm=0.; 330 G4double tasm = 0.; << 255 G4double asym = ComputeAsymmetry(lowEdgeEnergy, couple, part, 0., tasm); 331 G4double asym = ComputeAsymmetry(ene << 256 aVector->PutValue(j,asym); 332 bVector->PutValue(j, asym); << 333 } 257 } 334 bVector->FillSecondDerivatives(); << 258 335 G4PhysicsTableHelper::SetPhysicsVector(t << 259 G4PhysicsTableHelper::SetPhysicsVector(theAsymmetryTable, i, aVector); 336 } 260 } 337 } 261 } >> 262 338 } 263 } 339 264 340 //....oooOO0OOooo........oooOO0OOooo........oo 265 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 341 G4double G4PolarizedCompton::ComputeAsymmetry( << 266 342 G4double energy, const G4MaterialCutsCouple* << 267 343 const G4ParticleDefinition& aParticle, G4dou << 268 G4double G4PolarizedCompton::ComputeAsymmetry(G4double energy, >> 269 const G4MaterialCutsCouple* couple, >> 270 const G4ParticleDefinition& particle, >> 271 G4double cut, >> 272 G4double & tAsymmetry) 344 { 273 { 345 G4double lAsymmetry = 0.0; 274 G4double lAsymmetry = 0.0; 346 tAsymmetry = 0; << 275 tAsymmetry=0; 347 276 >> 277 // 348 // calculate polarized cross section 278 // calculate polarized cross section 349 G4ThreeVector thePolarization = G4ThreeVecto << 279 // 350 fEmModel->SetTargetPolarization(thePolarizat << 280 G4ThreeVector thePolarization=G4ThreeVector(0.,0.,1.); 351 fEmModel->SetBeamPolarization(thePolarizatio << 281 emModel->SetTargetPolarization(thePolarization); 352 G4double sigma2 = << 282 emModel->SetBeamPolarization(thePolarization); 353 fEmModel->CrossSection(couple, &aParticle, << 283 G4double sigma2=emModel->CrossSection(couple,&particle,energy,cut,energy); 354 284 >> 285 // 355 // calculate unpolarized cross section 286 // calculate unpolarized cross section 356 thePolarization = G4ThreeVector(); << 287 // 357 fEmModel->SetTargetPolarization(thePolarizat << 288 thePolarization=G4ThreeVector(); 358 fEmModel->SetBeamPolarization(thePolarizatio << 289 emModel->SetTargetPolarization(thePolarization); 359 G4double sigma0 = << 290 emModel->SetBeamPolarization(thePolarization); 360 fEmModel->CrossSection(couple, &aParticle, << 291 G4double sigma0=emModel->CrossSection(couple,&particle,energy,cut,energy); 361 << 292 362 // determine asymmetries << 293 // determine assymmetries 363 if(sigma0 > 0.) << 294 if (sigma0>0.) { 364 { << 295 lAsymmetry=sigma2/sigma0-1.; 365 lAsymmetry = sigma2 / sigma0 - 1.; << 366 } 296 } 367 return lAsymmetry; 297 return lAsymmetry; 368 } 298 } >> 299 >> 300 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 301 >> 302 >> 303 std::vector<G4DynamicParticle*>* G4PolarizedCompton::SecondariesPostStep( >> 304 G4VEmModel* model, >> 305 const G4MaterialCutsCouple* couple, >> 306 const G4DynamicParticle* dp) >> 307 { >> 308 std::vector<G4DynamicParticle*>* secondaries = model->SampleSecondaries(couple, dp); >> 309 return secondaries; >> 310 } >> 311 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 369 312