SectorField.h

Go to the documentation of this file.

/* 
 *  Copyright (c) 2012, Chris Rogers
 *  All rights reserved.
 *  Redistribution and use in source and binary forms, with or without 
 *  modification, are permitted provided that the following conditions are met: 
 *  1. Redistributions of source code must retain the above copyright notice,
 *     this list of conditions and the following disclaimer. 
 *  2. Redistributions in binary form must reproduce the above copyright notice, 
 *     this list of conditions and the following disclaimer in the documentation 
 *     and/or other materials provided with the distribution.
 *  3. Neither the name of STFC nor the names of its contributors may be used to 
 *     endorse or promote products derived from this software without specific 
 *     prior written permission.
 *
 *  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" 
 *  AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE 
 *  IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE 
 *  ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE 
 *  LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR 
 *  CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF 
 *  SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS 
 *  INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN 
 *  CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) 
 *  ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE 
 *  POSSIBILITY OF SUCH DAMAGE.
 */
 
#ifndef _CLASSIC_FIELDS_SECTORFIELD_HH_
#define _CLASSIC_FIELDS_SECTORFIELD_HH_
 
#include <vector>
 
#include "Fields/Fieldmap.h"
 
class SectorField {
  public:
    SectorField(const std::string& file_name);
 
    virtual ~SectorField();
 
    /* virtual bool getFieldstrengthPolar */
    /*               (const Vector_t &R_p, Vector_t &E_p, Vector_t &B_p) const = 0; */
 
    virtual bool getFieldstrength
                  (const Vector_t &R_c, Vector_t &E_c, Vector_t &B_c) const = 0;
 
    static void convertToPolar(double* position);
 
    static void convertToPolar(const double* position_polar, double* value);
 
    static void convertToCartesian(double* position);
 
    static void convertToCartesian(const double* position_polar, double* value);
 
    virtual std::vector<double> getPolarBoundingBoxMin() const;
 
    virtual std::vector<double> getPolarBoundingBoxMax() const;
 
    void getFieldDimensions(double &zBegin, double &zEnd,
                            double &rBegin, double &rEnd) const;
 
    void getFieldDimensions(double &xIni, double &xFinal,
                            double &yIni, double &yFinal,
                            double &zIni, double &zFinal) const;
 
    inline bool isInBoundingBox(const double R_p[]) const;
 
  protected:
    void setPolarBoundingBox(double bbMinR, double bbMinY, double bbMinPhi,
                             double bbMaxR, double bbMaxY, double bbMaxPhi,
                             double bbTolR, double bbTolY, double bbTolPhi);
 
    std::vector<double> bbMin_m;
    std::vector<double> bbMax_m;
    std::vector<double> polarBBMin_m;
    std::vector<double> polarBBMax_m;
    std::string Filename_m;
   
  private:
    std::vector< std::vector<double> > getCorners
               (double bbMinR, double bbMinPhi, double bbMaxR, double bbMaxPhi);
};
 
bool SectorField::isInBoundingBox(const double R_p[]) const {
    bool inside = R_p[0] >= polarBBMin_m[0] &&
                  R_p[0] <= polarBBMax_m[0] &&
                  R_p[1] >= polarBBMin_m[1] &&
                  R_p[1] <= polarBBMax_m[1] &&
                  R_p[2] >= polarBBMin_m[2] &&
                  R_p[2] <= polarBBMax_m[2];
    // std::cerr << "BB " << inside << "\n" << R_p[0] << " " << R_p[1] << " " << R_p[2] << "\n"
    //          << polarBBMin_m[0] << " " << polarBBMin_m[1] << " " << polarBBMin_m[2] << "\n"
    //          << polarBBMax_m[0] << " " << polarBBMax_m[1] << " " << polarBBMax_m[2] << std::endl;
    return inside;
}
 
#endif