Line data    Source code

       1              : #include "CRPropa.h"
       2              : #include "gtest/gtest.h"
       3              : 
       4              : namespace crpropa {
       5              : 
       6              : /*
       7              :  * Functional test which calculates the particle's gyroradius in a uniform field
       8              :  * r_g = R / (B*c) = 1 EV / (1 nG * c) \approx 1.08*Mpc
       9              :  */
      10            2 : TEST(testFunctionalGroups, gyroradius) {
      11              :         double energy = 1*EeV;
      12              :         double field = 1*nG;
      13              : 
      14            1 :         ParticleState p;
      15            1 :         p.setId(nucleusId(1, 1));
      16            1 :         p.setEnergy(energy);
      17            1 :         p.setPosition(Vector3d(0, 0, 0));
      18            1 :         p.setDirection(Vector3d(0, 0, 1));
      19              : 
      20            1 :         ref_ptr<Candidate> c = new Candidate(p);
      21            1 :         ref_ptr<PropagationCK> propa = new PropagationCK(new UniformMagneticField(Vector3d(field, 0, 0)));
      22            1 :         ref_ptr<ParticleCollector> collector = new ParticleCollector();
      23            1 :         collector->setClone(true);
      24            1 :         ref_ptr<ModuleList> sim = new ModuleList();
      25              : 
      26              :         Vector3d pos;
      27              :         double max_y = 0;
      28              : 
      29            1 :         sim->add(propa);
      30            1 :         sim->add(new MaximumTrajectoryLength(10*Mpc));
      31            1 :         sim->add(collector);
      32              : 
      33            1 :         sim->run(c);
      34              : 
      35           16 :         for (ParticleCollector::iterator itr = collector->begin(); itr != collector->end(); ++itr){
      36           15 :                 pos = (*(itr->get())).current.getPosition();
      37           15 :                 if (max_y < pos.getY())
      38              :                         max_y = pos.getY();
      39              :         }
      40              : 
      41            1 :         EXPECT_NEAR(max_y/2.0, energy/(field * c_light * eplus), 0.01*Mpc);
      42              : 
      43            1 : }
      44              : 
      45            0 : int main(int argc, char **argv) {
      46            0 :         ::testing::InitGoogleTest(&argc, argv);
      47            0 :         return RUN_ALL_TESTS();
      48              : }
      49              : 
      50              : } // namespace crpropa