GitLab

In ParallelCyclotronTracker the particles are mapped several times from local to global coordinates and vice versa. The PartBunch::boundp() operation where the particles are redistributed among the cores is always performed in local coordinates except in during the function call ParallelCyclotronTracker::applyPluginElements in case the boolean flag_stripper being true (line 3389 ff.).

if(((*sindex)->first) == ElementBase::STRIPPER)    {
            bool flag_stripper = (static_cast<Stripper *>(((*sindex)->second).second))
                -> checkStripper(itsBunch, turnnumber_m, itsBunch->getT() * 1e9, dt);
            if(flag_stripper) {
                itsBunch->boundp();
                *gmsg << "* Total number of particles after stripping = " << itsBunch->getTotalNum() << endl;
            }
        }

The workflow of ParallelCyclotronTracker::Tracker_Generic()

...
ParallelCyclotronTracker::initDistInGlobalFrame(); // (line 1235) --> particle in global coordinates
ParallelCyclotronTracker::applyPluginElements(dt); // (line 1285) --> PartBunch::boundp() in global coordinates !!!
...
// start tracking

Shouldn't the PartBunch::bounp() operation always be performed in local coordinates?

Best, Matthias