source: svn/trunk/Utilities/Hector/src/H_Drift.cc@ 235

/*
---- 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_Drift.cc
/// \brief Class aiming at simulating LHC beam drift.

// local includes
#include "H_Drift.h"
#include "H_TransportMatrices.h"


void H_Drift::init() {
        // must be in public section
        element_mat.ResizeTo(MDIM,MDIM);
        setTypeString();
        setMatrix(0,MP,QP);
        return;
}

void H_Drift::printProperties() const {
        cout << typestring << name;
        cout << "\t\t at s = " << fs;
        cout << "\t length = " << element_length;
        cout<<endl;
        if(element_aperture->getType()!=NONE) {
                cout <<"\t aperture type = " << element_aperture->getTypeString();
                element_aperture->printProperties();
        }
        if(element_length<0)  { if(VERBOSE) cout<<"<H_Drift> ERROR : Interpenetration of elements !"<<endl; }
        if(element_length==0) { if(VERBOSE) cout<<"<H_Drift> WARNING : 0-length "<< name << " !" << endl; }
}

void H_Drift::setMatrix(const float eloss, const float p_mass, const float p_charge) {
        element_mat = driftmat(element_length);
        return ;
}

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