Changeset 281 in svn for trunk/Utilities/Hector/include/H_Beam.h

Timestamp:

Mar 1, 2009, 3:56:11 PM (16 years ago)

Author:

Xavier Rouby

Message:

new Hector version

File:

• trunk/Utilities/Hector/include/H_Beam.h (modified) (6 diffs)

trunk/Utilities/Hector/include/H_Beam.h

@@ -2 +2 @@
 #define _H_Beam_
 
-/*
----- 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)
-*/
+  /* * * * * * * * * * * * * * * * * * * * * * * * * * * *
+ *                                                         *
+*                   --<--<--  A fast simulator --<--<--     *
+*                 / --<--<--     of particle   --<--<--     *
+*  ----HECTOR----<                                          *
+*                 \ -->-->-- transport through -->-->--     *
+*                   -->-->-- generic beamlines -->-->--     *
+*                                                           *
+* JINST 2:P09005 (2007)                                     *
+*      X Rouby, J de Favereau, K Piotrzkowski (CP3)         *
+*       http://www.fynu.ucl.ac.be/hector.html               *
+*                                                           *
+* Center for Cosmology, Particle Physics and Phenomenology  *
+*              Universite catholique de Louvain             *
+*                 Louvain-la-Neuve, Belgium                 *
+ *                                                         *
+   * * * * * * * * * * * * * * * * * * * * * * * * * * * */
 
 /// \file H_Beam.h
@@ -27 +34 @@
 #include "TGraphErrors.h"
 #include "TMultiGraph.h"
+#include "TMath.h"
 
 // local #includes
@@ -56 +64 @@
         //@}
         /// Fills the beam with particles in given position/angle/energy        intervals (flat distribution)
-                //void particleGun(const unsigned int Number_of_particles, const float E_min, const float E_max, const float fs_min, const float fs_max, const float fx_min, const float fx_max, const float fy_min, const float fy_max, const float tx_min, const float tx_max, const float ty_min, const float ty_max, const float p_mass, const double p_charge);
-                void particleGun(const unsigned int Number_of_particles, const float E_min=BE, const float E_max=BE, const float fs_min=0, const float fs_max=0, const float fx_min=0, const float fx_max=0, const float fy_min=0, const float fy_max=0, const float tx_min=-pi/2., const float tx_max=pi/2., const float ty_min=-pi/2., const float ty_max=pi/2., const float p_mass=MP, const double p_charge=QP, const bool flat = true);
+                void particleGun(const unsigned int Number_of_particles, const float E_min=BE, const float E_max=BE, const float fs_min=0, const float fs_max=0, const float fx_min=0, const float fx_max=0, const float fy_min=0, const float fy_max=0, const float tx_min=-TMath::Pi()/2., const float tx_max=TMath::Pi()/2., const float ty_min=-TMath::Pi()/2., const float ty_max=TMath::Pi()/2., const float p_mass=MP, const double p_charge=QP, const bool flat = true, TRandom* r=gRandom);
 
         /// Fills the beam with particles
-        //@{
-                void createBeamParticles(const unsigned int , const double , const double );
-                void createBeamParticles(const unsigned int);
-        //@}
+                void createBeamParticles(const unsigned int Number_of_particles, const double p_mass=MP, const double p_charge=QP, TRandom* r=gRandom);
         /// Fills the beam with particles with incremental offset and no initial transverse momentum
         //@{
@@ -82 +86 @@
                 void add(const H_BeamParticle&);
         /// Compute the position of each particle in the beamline
-        //@{    
-                void computePath(const H_AbstractBeamLine *, const bool);
-                void computePath(const H_AbstractBeamLine *);
-        //@}
+                void computePath(const H_AbstractBeamLine * beamline, const bool NonLinear=false);
         // Photon emission by the particle
-        // @{
-                void emitGamma(const double, const double, const double, const double);
-                void emitGamma(const double, const double);
-        //@}
+                void emitGamma(const double gee, const double gq2, const double phimin=0, const double phimax=2*TMath::Pi());
         // Propagates the beam until a given s
                 void propagate(const float );
@@ -146 +144 @@
                 void printInitialState() const;
         /// Prints the properties for each particle
-                void printProperties() const;
+                void printProperties() const {cout << *this; return;}
         /// Prints the list of the stopping elements in the beamline
                 void printStoppingElements(const vector<H_OpticalElement>&, const vector<int>&) const;
@@ -186 +184 @@
                 unsigned int Nparticles;
 
-        static const float pi = 3.14159265358979312;
+        friend std::ostream& operator<< (std::ostream& os, const H_Beam& be);
 };