source: trunk/source/processes/electromagnetic/muons/src/G4MuMultipleScattering.cc @ 961

Last change on this file since 961 was 961, checked in by garnier, 15 years ago

update processes

File size: 4.7 KB
Line 
1//
2// ********************************************************************
3// * License and Disclaimer                                           *
4// *                                                                  *
5// * The  Geant4 software  is  copyright of the Copyright Holders  of *
6// * the Geant4 Collaboration.  It is provided  under  the terms  and *
7// * conditions of the Geant4 Software License,  included in the file *
8// * LICENSE and available at  http://cern.ch/geant4/license .  These *
9// * include a list of copyright holders.                             *
10// *                                                                  *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work  make  any representation or  warranty, express or implied, *
14// * regarding  this  software system or assume any liability for its *
15// * use.  Please see the license in the file  LICENSE  and URL above *
16// * for the full disclaimer and the limitation of liability.         *
17// *                                                                  *
18// * This  code  implementation is the result of  the  scientific and *
19// * technical work of the GEANT4 collaboration.                      *
20// * By using,  copying,  modifying or  distributing the software (or *
21// * any work based  on the software)  you  agree  to acknowledge its *
22// * use  in  resulting  scientific  publications,  and indicate your *
23// * acceptance of all terms of the Geant4 Software license.          *
24// ********************************************************************
25//
26// $Id: G4MuMultipleScattering.cc,v 1.12 2008/10/16 13:37:04 vnivanch Exp $
27// GEANT4 tag $Name: geant4-09-02-ref-02 $
28//
29// -----------------------------------------------------------------------------
30//
31// GEANT4 Class file
32//
33// File name:     G4MuMultipleScattering
34//
35// Author:        Laszlo Urban
36//
37// Creation date: 24.10.2006 cloned from G4MultipleScattering
38//
39// Modified:
40// 12-02-07 skin can be changed via UI command (VI)
41// 20.03.07 Remove local parameter skin, set facgeom=0.1(V.Ivanchenko)
42//
43// -----------------------------------------------------------------------------
44//
45//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
46//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
47
48#include "G4MuMultipleScattering.hh"
49#include "G4WentzelVIModel.hh"
50#include "G4MscStepLimitType.hh"
51
52//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
53
54using namespace std;
55
56G4MuMultipleScattering::G4MuMultipleScattering(G4double tet,
57                                               const G4String& processName)
58  : G4VMultipleScattering(processName), thetaLimit(tet)
59{
60  dtrl              = 0.05;
61  samplez           = false ; 
62  isInitialized     = false; 
63  SetRangeFactor(0.2);
64  SetLateralDisplasmentFlag(true);
65}
66
67//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
68
69G4MuMultipleScattering::~G4MuMultipleScattering()
70{}
71
72//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
73
74G4bool G4MuMultipleScattering::IsApplicable (const G4ParticleDefinition& p)
75{
76  return (p.GetPDGCharge() != 0.0 && !p.IsShortLived());
77}
78
79//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
80
81void G4MuMultipleScattering::InitialiseProcess(const G4ParticleDefinition* p)
82{
83  // Modification of parameters between runs
84  if(isInitialized) {
85
86    if (p->GetParticleType() != "nucleus" && p->GetPDGMass() < GeV) {
87      mscModel->SetStepLimitType(StepLimitType());
88      mscModel->SetLateralDisplasmentFlag(LateralDisplasmentFlag());
89      mscModel->SetRangeFactor(RangeFactor());
90    }
91    mscModel->SetPolarAngleLimit(PolarAngleLimit());
92    return;
93  }
94
95  if (p->GetParticleType() == "nucleus" || p->GetPDGMass() > GeV) {
96    SetLateralDisplasmentFlag(false);
97    SetBuildLambdaTable(false);
98  }
99
100  // initialisation of the model
101
102  mscModel = new G4WentzelVIModel();
103  mscModel->SetStepLimitType(StepLimitType());
104  mscModel->SetLateralDisplasmentFlag(LateralDisplasmentFlag());
105  mscModel->SetRangeFactor(RangeFactor());
106  mscModel->SetPolarAngleLimit(PolarAngleLimit());
107  mscModel->SetLowEnergyLimit(MinKinEnergy());
108  mscModel->SetHighEnergyLimit(MaxKinEnergy());
109
110  AddEmModel(1,mscModel);
111  isInitialized = true;
112}
113
114//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
115
116void G4MuMultipleScattering::PrintInfo()
117{
118  G4cout << "      RangeFactor= " << RangeFactor()
119         << ", step limit type: " << StepLimitType()
120         << ", lateralDisplacement: " << LateralDisplasmentFlag()
121         << G4endl;
122}
123
124//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
125
Note: See TracBrowser for help on using the repository browser.