Changeset 281 in svn for trunk/Utilities/Hector/src/H_RecRPObject.cc

Timestamp:

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

Author:

Xavier Rouby

Message:

new Hector version

File:

• trunk/Utilities/Hector/src/H_RecRPObject.cc (modified) (4 diffs)

trunk/Utilities/Hector/src/H_RecRPObject.cc

@@ -1 @@
-/*
----- 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_RecRPObject.cc
@@ -41 +47 @@
                                 l_1(new TF1("l_1","[0] + [1]*x + [2]*x*x ",emin,emax)), 
                                 l_2(new TF1("l_2","[0] + [1]*x + [2]*x*x ",emin,emax)) 
-                                {
-        // emax = -1; /* GeV */
-}
-
-H_RecRPObject::H_RecRPObject(float ss1, float ss2, H_AbstractBeamLine* beam) : emin(0), emax(-1), x1(0), x2(0), y1(0), y2(0), s1(ss1), s2(ss2), 
+{}
+
+H_RecRPObject::H_RecRPObject(const float ss1, const float ss2, const H_AbstractBeamLine* beam) : emin(0), emax(-1), x1(0), x2(0), y1(0), y2(0), s1(ss1), s2(ss2), 
                                 txip(NOT_YET_COMPUTED), tyip(NOT_YET_COMPUTED), energy(NOT_YET_COMPUTED), q2(NOT_YET_COMPUTED), pt(NOT_YET_COMPUTED), 
-                                thebeam(beam),
+                                thebeam(beam->clone()),
                                 f_1(new TF1("f_1","[0] + [1]*x + [2]*x*x ",emin,emax)), 
                                 f_2(new TF1("f_2","[0] + [1]*x + [2]*x*x ",emin,emax)), 
@@ -58 +62 @@
                                 l_1(new TF1("l_1","[0] + [1]*x + [2]*x*x ",emin,emax)), 
                                 l_2(new TF1("l_2","[0] + [1]*x + [2]*x*x ",emin,emax))
-                                {
-        if(ss1==ss2) cout<<"<H_RecRPObject> WARNING : detectors are on same position"<<endl;
-        // emax = -1; /* GeV */
+        {if(ss1==ss2) cout<<"<H_RecRPObject> WARNING : detectors are on same position"<<endl;
+}
+
+H_RecRPObject::H_RecRPObject(const H_RecRPObject& r):
+   emin(r.emin), emax(r.emax), x1(r.x1), x2(r.x2), y1(r.y1), y2(r.y2), s1(r.s1), s2(r.s2),
+   txip(r.txip), tyip(r.tyip), energy(r.energy), q2(r.q2),  pt(r.pt), 
+   //thebeam(r.thebeam->clone()),
+   thebeam(new H_AbstractBeamLine(*(r.thebeam))),
+   f_1(new TF1(*(r.f_1))),  f_2(new TF1(*(r.f_2))),  g_1(new TF1(*(r.g_1))),  g_2(new TF1(*(r.g_2))),
+   d_1(new TF1(*(r.d_1))),  d_2(new TF1(*(r.d_2))),  k_1(new TF1(*(r.k_1))),  k_2(new TF1(*(r.k_2))),
+   l_1(new TF1(*(r.l_1))),  l_2(new TF1(*(r.l_2)))
+ {}
+
+H_RecRPObject& H_RecRPObject::operator=(const H_RecRPObject& r) {
+  if (this == &r) return *this;
+   emin = r.emin, emax = r.emax;
+   x1 = r.x1;   x2 = r.x2; 
+   y1 = r.y1;   y2 = r.y2; 
+   s1 = r.s1;   s2 = r.s2;
+   txip= r.txip; tyip=r.tyip;
+   energy= r.energy; q2= r.q2;  pt= r.pt;
+   //thebeam = r.thebeam->clone();
+   thebeam = new H_AbstractBeamLine(*(r.thebeam));
+   f_1 = new TF1(*(r.f_1));
+   f_2 = new TF1(*(r.f_2));
+   g_1 = new TF1(*(r.g_1));
+   g_2 = new TF1(*(r.g_2));
+   d_1 = new TF1(*(r.d_1));
+   d_2 = new TF1(*(r.d_2));
+   k_1 = new TF1(*(r.k_1));
+   k_2 = new TF1(*(r.k_2));
+   l_1 = new TF1(*(r.l_1));
+   l_2 = new TF1(*(r.l_2));
+  return *this;
 }
 
@@ -308 +343 @@
         return pt;
 }
+
+std::ostream& operator<< (std::ostream& os, const H_RecRPObject& rp) {
+        os << "e_min=" <<  rp.emin << "\t e_max= " << rp.emax << endl;
+        os << "x1="    << rp.x1 << "\t x2= " << rp.x2 << "\t y1=" << rp.y1 << "\t y2=" << rp.y2 
+           << "\t s1=" << rp.s1 << "\t s2=" << rp.s2 << endl;
+        os << "txip=" << rp.txip << "\t tyip=" << rp.tyip << "\t energy=" << rp.energy << "\t q2=" << rp.q2 << "\t pt=" << rp.pt << endl;
+   return os;
+}