source: svn/trunk/Utilities/Hector/src/H_HorizontalQuadrupole.cc@ 96

/*
---- Hector the simulator ----
   A fast simulator of particles through generic beamlines.
   J. de Favereau, X. Rouby ~~~ hector_devel@cp3.phys.ucl.ac.be

        http://www.fynu.ucl.ac.be/hector.html

   Centre de Physique des Particules et de Phénoménologie (CP3)
   Université Catholique de Louvain (UCL)
*/

/// \file H_HorizontalQuadrupole.cc
/// \brief Classes aiming at simulating horizontal quadrupoles in beamline

// local #includes
#include "H_HorizontalQuadrupole.h"

void H_HorizontalQuadrupole::setMatrix(const float eloss, const float p_mass, const float p_charge) {
        if (fk>0)  { if(VERBOSE) cout<<"<H_HorizontalQuadrupole> ERROR : k1 should be < 0 (" << name << ")!"<<endl; }
                if (fk !=0 ) element_mat = hquadmat(element_length,fk,eloss, p_mass, p_charge);
                else  {
                        element_mat = driftmat(element_length);
                        if(VERBOSE) cout<<"<H_HorizontalQuadrupole> WARNING : k1= 0 ; drift-like quadrupole (" << name << ") !" << endl;
                }
        return ;
}

H_HorizontalQuadrupole* H_HorizontalQuadrupole::clone() const {
        H_HorizontalQuadrupole* temp_quad = new H_HorizontalQuadrupole(name,fs,fk,element_length);
        temp_quad->setAperture(element_aperture);
        temp_quad->setX(xpos);
        temp_quad->setY(ypos);
        temp_quad->setTX(txpos);
        temp_quad->setTY(typos);
        temp_quad->setBetaX(betax);
        temp_quad->setBetaY(betay);
        return temp_quad;
}