include/MyElementField.hh

Go to the documentation of this file.

//
// ********************************************************************
// * License and Disclaimer                                           *
// *                                                                  *
// * The  Geant4 software  is  copyright of the Copyright Holders  of *
// * the Geant4 Collaboration.  It is provided  under  the terms  and *
// * conditions of the Geant4 Software License,  included in the file *
// * LICENSE and available at  http://cern.ch/geant4/license .  These *
// * include a list of copyright holders.                             *
// *                                                                  *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work  make  any representation or  warranty, express or implied, *
// * regarding  this  software system or assume any liability for its *
// * use.  Please see the license in the file  LICENSE  and URL above *
// * for the full disclaimer and the limitation of liability.         *
// *                                                                  *
// * This  code  implementation is the result of  the  scientific and *
// * technical work of the GEANT4 collaboration.                      *
// * By using,  copying,  modifying or  distributing the software (or *
// * any work based  on the software)  you  agree  to acknowledge its *
// * use  in  resulting  scientific  publications,  and indicate your *
// * acceptance of all terms of the Geant4 Software license.          *
// ********************************************************************
//
//
//
//

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....

#ifndef MyElementField_h
#define MyElementField_h 1

#include "myglobals.hh"

#include "G4Navigator.hh"
#include "G4TransportationManager.hh"

#include "G4UserLimits.hh"
#include "G4VisAttributes.hh"

//  class MyElementField - interface for the EM field of one element

//  This is the interface class used by GlobalField to compute the field
//  value at a given point[].

//  An element that represents an element with an EM field will
//  derive a class from this one and implement the computation for the
//  element. The construct() function will add the derived object into
//  GlobalField.

class MyElementField {

    private:

        MyElementField& operator=(const MyElementField&);

    public:

        MyElementField(const G4ThreeVector, G4LogicalVolume*);

        void construct();

        virtual ~MyElementField() {
            if (aNavigator) delete aNavigator;
            if (boundingBox) delete boundingBox;
        }

        // Get center of physical field volume
        G4ThreeVector getCenter(){return center;}

        // Set center of physical field volume
        void setCenter(G4ThreeVector);

        // Print Bounding Box
        G4double const*  getBoundingBox();

        void setMaxStep(G4double s) {
            maxStep = s;
            userLimits->SetMaxAllowedStep(maxStep);
            lvolume->SetUserLimits(userLimits);
        }

        G4double getMaxStep() { return maxStep; }

        void setColor(G4String c) {
            color = c;
            lvolume->SetVisAttributes(getVisAttribute(color));
        }

        G4String getColor() { return color; }

        static G4VisAttributes* getVisAttribute(G4String color);

        void setGlobalPoint(const G4double point[4]) {
            if(minX == -DBL_MAX || minX > point[0]) minX = point[0];
            if(minY == -DBL_MAX || minY > point[1]) minY = point[1];
            if(minZ == -DBL_MAX || minZ > point[2]) minZ = point[2];
            if(maxX ==  DBL_MAX || maxX < point[0]) maxX = point[0];
            if(maxY ==  DBL_MAX || maxY < point[1]) maxY = point[1];
            if(maxZ ==  DBL_MAX || maxZ < point[2]) maxZ = point[2];
        }

        bool isInBoundingBox(const G4double point[4]) const {
            if(point[2] < minZ || point[2] > maxZ) return false;
            if(point[0] < minX || point[0] > maxX) return false;
            if(point[1] < minY || point[1] > maxY) return false;
            return true;
        }

        virtual void
            addFieldValue(const G4double point[4], G4double field[6]) const = 0;

        // Get dimensions for the bounding box
        virtual G4double getMaxLength() = 0;
        virtual G4double getMaxWidth()  = 0;
        virtual G4double getMaxHeight() = 0;

    protected:

        G4LogicalVolume* lvolume;

        G4AffineTransform global2local;

        //    MyElementField(const MyElementField&);

    private:

        static G4Navigator* aNavigator;

        G4String color;

        G4ThreeVector center;
        G4double minX, minY, minZ, maxX, maxY,maxZ;

        G4double maxStep;
        G4UserLimits* userLimits;

        G4double* boundingBox;
};
#endif

---