source: ZHANGProjects/ICOSIM/CPP/trunk/source/MagneticCollimator.h @ 2

#ifndef MAGNETICCOLLIMATOR_H
#define MAGNETICCOLLIMATOR_H
#include <iostream>
#include <vector>
#include <string>
#include <cmath>
#include "Collimator.h"
using namespace std;


/*
####################################################################################################################

This class describes all magnetic collimators in the accelerator lattice.

####################################################################################################################
*/


class MagneticCollimator : public Collimator
{

public:

    //=================================================Constructors, destructor=======================================

    MagneticCollimator(const double& ALFX, const double& ALFY, const double& APER_1, const double& APER_2, const double& APER_3, const double& APER_4, const string& APERTYPE, const double& BETX, const double& BETY, const double& DPX, const double& DPY, const double& DX, const double& DY, const string& KEYWORD, const double& L, const double& MUX, const double& MUY, const string& NAME, const double& PTC, const double& PXC, const double& PYC, const double& S, const double& TC, const double& XC, const double& YC, const double& K0L, const double& K0SL, const double& K1L, const double& K1SL, const double& K2L, const double& K2SL, const string& PARENT, const string& meth, const long double& hgap, const long double& hgap2, const double& collang, const long double& pdepth, const long double& pdepth2, const double& tcang, const double& nsig, const double& Bmax, const double& thicknessMagneticField, const double& energyPerIon, const double& mass);

    MagneticCollimator(Element elt, const double& tcang, const double& nsig, const string& meth, const string& material);

    MagneticCollimator(const MagneticCollimator& obj);

    virtual ~MagneticCollimator() {};


    //calculation of B

    double Bcalc(const double& x, const double& y, const double& z, const double& w);


    //describes the passage through a collimator, here magnetic

    void collipass(Particle& p1, double& dpopeff, const double& scaleorbit, const double& R11X, const double& R12X, const double& R21X, const double& R22X, const double& R11Y, const double& R12Y, const double& R21Y, const double& R22Y, const double& dx1, const double& dpx1, const double& dy1, const double& dpy1, const double& delta_s, const double& Apr, const double& Zpr, const double& betgam);

    //displays the parameters of the collimator

    void affiche();


private:

    vector <int> first;//tells you whether each particle hits the collimator on the front end (2), hits the collimator along the first half of the jaw (1) or doesn't hit the first half of the collimator (0)
    vector <int> second;//tells you whether the particle hits the collimator along the last half of the jaw (1), doesn't hit the second part of the collimator (0) or gets lost in the first half of the collimator or is inside the collimator in the second half of the collimator (-1)


};

#endif