Changeset 1365 in svn for trunk/external/Hector/H_VerticalQuadrupole.cc
- Timestamp:
- Apr 16, 2014, 3:56:14 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/external/Hector/H_VerticalQuadrupole.cc
r1360 r1365 1 /* * * * * * * * * * * * * * * * * * * * * * * * * * * * 2 * * 3 * --<--<-- A fast simulator --<--<-- * 4 * / --<--<-- of particle --<--<-- * 5 * ----HECTOR----< * 6 * \ -->-->-- transport through -->-->-- * 7 * -->-->-- generic beamlines -->-->-- * 8 * * 9 * JINST 2:P09005 (2007) * 10 * X Rouby, J de Favereau, K Piotrzkowski (CP3) * 11 * http://www.fynu.ucl.ac.be/hector.html * 12 * * 13 * Center for Cosmology, Particle Physics and Phenomenology * 14 * Universite catholique de Louvain * 15 * Louvain-la-Neuve, Belgium * 16 * * 17 * * * * * * * * * * * * * * * * * * * * * * * * * * * */ 1 /* 2 ---- Hector the simulator ---- 3 A fast simulator of particles through generic beamlines. 4 J. de Favereau, X. Rouby ~~~ hector_devel@cp3.phys.ucl.ac.be 5 6 http://www.fynu.ucl.ac.be/hector.html 7 8 Centre de Physique des Particules et de Phénoménologie (CP3) 9 Université Catholique de Louvain (UCL) 10 */ 18 11 19 12 /// \file H_VerticalQuadrupole.cc … … 23 16 #include "H_VerticalQuadrupole.h" 24 17 25 void H_VerticalQuadrupole::setMatrix(const float eloss, const float p_mass, const float p_charge) {26 if (fk<0) { if(VERBOSE) cout<<" <H_VerticalQuadrupole>ERROR : k1 should be > 0 for H_VerticalQuadrupole (" << name << ")!"<<endl; }27 if (fk !=0 ) element_mat = vquadmat(element_length,fk,eloss, p_mass,p_charge);18 void H_VerticalQuadrupole::setMatrix(const float eloss, const float p_mass, const float p_charge) const { 19 if (fk<0) { if(VERBOSE) cout<<"\t ERROR : k1 should be > 0 for H_VerticalQuadrupole (" << name << ")!"<<endl; } 20 if (fk !=0 ) *element_mat = vquadmat(element_length,fk,eloss, p_mass,p_charge); 28 21 else { 29 element_mat = driftmat(element_length);30 if(VERBOSE) cout<<" <H_VerticalQuadrupole>WARNING : k1= 0, drift-like quadrupole (" << name << ") !" << endl;22 *element_mat = driftmat(element_length); 23 if(VERBOSE) cout<<"\t WARNING : k1= 0, drift-like quadrupole (" << name << ") !" << endl; 31 24 } 32 25 return ; 33 26 } 34 35 H_VerticalQuadrupole* H_VerticalQuadrupole::clone() const {36 H_VerticalQuadrupole* temp_quad = new H_VerticalQuadrupole(name,fs,fk,element_length);37 temp_quad->setAperture(element_aperture);38 temp_quad->setX(xpos);39 temp_quad->setY(ypos);40 temp_quad->setTX(txpos);41 temp_quad->setTY(typos);42 temp_quad->setBetaX(betax);43 temp_quad->setBetaY(betay);44 return temp_quad;45 }
Note:
See TracChangeset
for help on using the changeset viewer.