Kamil Sedlak
1e050e6976
- a copyright information added to all relevant *.cc and *.hh files - the file "COPYING" added to the main directory of the musrSim
110 lines
4.6 KiB
C++
110 lines
4.6 KiB
C++
/***************************************************************************
|
|
* musrSim - the program for the simulation of (mainly) muSR instruments. *
|
|
* More info on http://lmu.web.psi.ch/simulation/index.html . *
|
|
* musrSim is based od Geant4 (http://geant4.web.cern.ch/geant4/) *
|
|
* *
|
|
* Copyright (C) 2009 by Paul Scherrer Institut, 5232 Villigen PSI, *
|
|
* Switzerland *
|
|
* *
|
|
* This program is free software; you can redistribute it and/or modify *
|
|
* it under the terms of the GNU General Public License as published by *
|
|
* the Free Software Foundation; either version 2 of the License, or *
|
|
* (at your option) any later version. *
|
|
* *
|
|
* This program is distributed in the hope that it will be useful, *
|
|
* but WITHOUT ANY WARRANTY; without even the implied warranty of *
|
|
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
|
|
* GNU General Public License for more details. *
|
|
* *
|
|
* You should have received a copy of the GNU General Public License *
|
|
* along with this program; if not, write to the Free Software *
|
|
* Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. *
|
|
***************************************************************************/
|
|
|
|
#include "globals.hh"
|
|
#include "G4GeometryManager.hh"
|
|
#include "musrUniformField.hh"
|
|
|
|
musrUniformField::musrUniformField(G4double EMF[6], G4double half_X, G4double half_Y, G4double half_Z, G4LogicalVolume* lv, G4ThreeVector c)
|
|
: F04ElementField(c, lv) {
|
|
// EMF[6] ... constant vector of the field: 3 components magnetic + 3 components electric field
|
|
// half_X, half_Y, half_Z ... half dimenstions of the box, within which the field is defined
|
|
|
|
for (int i = 0; i < 6; i++){
|
|
EMfield[i] = EMF[i];
|
|
}
|
|
|
|
fieldLength = 2.*half_Z;
|
|
fieldWidth = 2.*half_X;
|
|
fieldHeight = 2.*half_Y;
|
|
|
|
G4cout << "\n-----------------------------------------------------------"
|
|
<< "\n Uniform electromagnetic field"
|
|
<< G4endl;
|
|
|
|
G4String volName = lv->GetName().substr(4);
|
|
G4cout << "\n ---> EM field in volume " << volName << " set to:" << G4endl;
|
|
printf (" B = (%0.3g, %0.3g, %0.3g) T, E = (%0.3g, %0.3g, %0.3g) kV/mm\n",
|
|
EMF[0]/tesla, EMF[1]/tesla, EMF[2]/tesla,
|
|
EMF[3]/(kilovolt/mm), EMF[4]/(kilovolt/mm), EMF[5]/(kilovolt/mm));
|
|
G4cout << "-----------------------------------------------------------" << G4endl;
|
|
}
|
|
|
|
|
|
|
|
void musrUniformField::addFieldValue(const G4double point[4],
|
|
G4double field[6]) const
|
|
{
|
|
G4ThreeVector global(point[0],point[1],point[2]);
|
|
G4ThreeVector local;
|
|
|
|
local = global2local.TransformPoint(global);
|
|
|
|
if (isOutside(local)) return;
|
|
|
|
G4ThreeVector B(EMfield[0],EMfield[1],EMfield[2]);
|
|
G4ThreeVector E(EMfield[3],EMfield[4],EMfield[5]);
|
|
|
|
B = global2local.Inverse().TransformAxis(B);
|
|
E = global2local.Inverse().TransformAxis(E);
|
|
|
|
field[0] += B[0];
|
|
field[1] += B[1];
|
|
field[2] += B[2];
|
|
|
|
field[3] += E[0];
|
|
field[4] += E[1];
|
|
field[5] += E[2];
|
|
|
|
//printf (" EM field components: B = (%0.3g, %0.3g, %0.3g) T, E = (%0.3g, %0.3g, %0.3g) kV/mm\n",
|
|
//field[0]/tesla, field[1]/tesla, field[2]/tesla,
|
|
//field[3]/(kilovolt/mm), field[4]/(kilovolt/mm), field[5]/(kilovolt/mm));
|
|
}
|
|
|
|
|
|
|
|
G4double musrUniformField::GetNominalFieldValue() {
|
|
G4double val = std::sqrt(EMfield[0]*EMfield[0]+EMfield[1]*EMfield[1]+EMfield[2]*EMfield[2]+
|
|
EMfield[3]*EMfield[3]+EMfield[4]*EMfield[4]+EMfield[5]*EMfield[5]);
|
|
return val;
|
|
}
|
|
|
|
void musrUniformField::SetNominalFieldValue(G4double newFieldValue){
|
|
EMfield[0] *= newFieldValue;
|
|
EMfield[1] *= newFieldValue;
|
|
EMfield[2] *= newFieldValue;
|
|
EMfield[3] *= newFieldValue;
|
|
EMfield[4] *= newFieldValue;
|
|
EMfield[5] *= newFieldValue;
|
|
}
|
|
|
|
G4bool musrUniformField::isOutside(G4ThreeVector& local) const
|
|
{
|
|
return (std::fabs(local.z()) > fieldLength/2.0 || std::fabs(local.x()) > fieldWidth/2.0 || std::fabs(local.y()) > fieldHeight/2.0);
|
|
}
|
|
|
|
G4bool musrUniformField::isWithin(G4ThreeVector& local) const
|
|
{
|
|
return (std::fabs(local.z()) < fieldLength/2.0 && std::fabs(local.x()) < fieldWidth/2.0 && std::fabs(local.y()) < fieldHeight/2.0);
|
|
}
|