why No effect change of particle energy or field or De Dh ON NBodyChargeCloud plots ?

[rajxabc]

documentation for NBodyChargeCloud is bit confusing
Why No effect change of particle energy or field or De Dh ON NBodyChargeCloud plots ?
Is this suppose to plot animation "simulate!(evt, sim, diffusion = true, self_repulsion = true)"?
Can you suggest a way to show on plot e and h concentration density as a whole from geant4 events ?
Your support is much appreciated .
Thanks

[fhagemann]

Im not quite sure what you mean..
Are you maybe affected by this: #461
Could be that the different charge cloud are not interacting with each other.

we simulate group effects using point charges, not using concentrations (?)

[rajxabc]

Checking #461 Thanks

[fhagemann]

If that does not help, I might need some more information about what exactly you mean with „NBodyChargeCloud Plots“

[rajxabc]

1. I am unable to plot the animated plot "Combination of Group Effects" from the document.

center = CartesianPoint{T}([0,0,0])
energy = 1460u"keV"
nbcc = NBodyChargeCloud(center, energy, 100)
evt = Event(nbcc)
simulate!(evt, sim, diffusion = true, self_repulsion = true)

2. I figured it out "NBodyChargeCloud plots " adjust its radius automatically , that's why appear same .
3. Allow automatic grouping of hits based on a group_distance when creating an Event #461 seems good , apologies for confusion .
   Thank you

[fhagemann]

1. This could be because the drift is terminated in the absence of an external electric field:
   

   SolidStateDetectors.jl/src/ChargeDrift/ChargeDrift.jl

   Lines 101 to 110 in 3d62146

   	function _add_fieldvector_drift!(step_vectors::Vector{CartesianVector{T}}, current_pos::Vector{CartesianPoint{T}}, done::Vector{Bool},
   	electric_field::Interpolations.Extrapolation{<:StaticVector{3}, 3}, det::SolidStateDetector{T}, ::Type{S})::Nothing where {T, S}
   	for n in eachindex(step_vectors)
   	if !done[n]
   	step_vectors[n] += get_velocity_vector(electric_field, _convert_point(current_pos[n], S))
   	done[n] = (step_vectors[n] == CartesianVector{T}(0,0,0))
   	end
   	end
   	nothing
   	end

   
   You could explicitly comment out that line by adding this to your code

   using Interpolation, StaticArrays
   using SolidStateDetectors: CartesianVector, CartesianPoint, get_velocity_vector, _convert_point
   
   function SolidStateDetectors._add_fieldvector_drift!(step_vectors::Vector{CartesianVector{T}}, current_pos::Vector{CartesianPoint{T}}, done::Vector{Bool}, 
     electric_field::Interpolations.Extrapolation{<:StaticVector{3}, 3}, det::SolidStateDetector{T}, ::Type{S})::Nothing where {T, S}
     for n in eachindex(step_vectors)
        if !done[n]
            step_vectors[n] += get_velocity_vector(electric_field, _convert_point(current_pos[n], S))
            # done[n] = (step_vectors[n] == CartesianVector{T}(0,0,0)) # this line is commented out
        end
     end
     nothing
   end

[fhagemann]

There is also an effort to somewhat refactor the drift, which will introduce a keyword that will let the user decide whether the drift should be ended in regions with no electric field: #476.
Can you let us know if this solved your issue?

[fhagemann]

Another question: how do you define your Event. Sometimes you cannot see effects if the charge cloud only consists of one charge.