Commit 6904b66a authored by Simon Spannagel's avatar Simon Spannagel
Browse files

TransientPropagation: remove special treatment of carriers at surface

parent f5c0e4db

src/modules/TransientPropagation/TransientPropagationModule.cpp

+3 −10
Original line number Diff line number Diff line
@@ -348,16 +348,9 @@ TransientPropagationModule::propagate(Event* event, 
            check_position.z() = last_position.z();

            // Correct for position in z by interpolation to increase precision:

            if(detector_->getModel()->isWithinSensor(static_cast<ROOT::Math::XYZPoint>(check_position))) {

                // FIXME this currently depends in the direction of the drift

                if(position.z() > 0 && type == CarrierType::HOLE) {

                    LOG(DEBUG) << "Not stopping carrier " << type << " at "

                LOG(DEBUG) << "Stopping carrier " << type << " at "

                           << Units::display(static_cast<ROOT::Math::XYZPoint>(position), {"um"});

                } else if(position.z() < 0 && type == CarrierType::ELECTRON) {

                    LOG(DEBUG) << "Not stopping carrier " << type << " at "

                               << Units::display(static_cast<ROOT::Math::XYZPoint>(position), {"um"});

                } else {

                state = CarrierState::HALTED;

                }



                // Carrier left sensor on top or bottom surface, interpolate

                auto z_cur_border = std::fabs(position.z() - model_->getSensorSize().z() / 2.0);