48 fpRegisteredSBPoints()
51 std::set<SBPoint*>::iterator itPt;
52 for(itPt = pSBPoints.begin(); itPt != pSBPoints.end(); ++itPt)
67 std::set<SBPoint*>::iterator itPt;
72 (*itPt)->CleanCluster();
96 std::set<SBPoint* >::iterator itSDSPt;
101 x+= (*itSDSPt)->GetPosition().x();
102 y+= (*itSDSPt)->GetPosition().y();
103 z+= (*itSDSPt)->GetPosition().z();
117 std::set<SBPoint* >::iterator itSDSPt;
122 res += (*itSDSPt)->GetEdep();
132 bool firstStrandTouch =
false;
133 bool secondStrandTouch =
false;
135 std::set<SBPoint* >::iterator itSDSPt;
141 if( ((*itSDSPt)->GetTouchedStrand() == 0 ) && !firstStrandTouch )
143 firstStrandTouch =
true;
144 if(secondStrandTouch)
151 if( ((*itSDSPt)->GetTouchedStrand() == 1 ) && !secondStrandTouch )
153 secondStrandTouch =
true;
180 if(((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)+(z1-z2)*(z1-z2))<=
181 (pMinDist/
nm*pMinDist/
nm))
195 assert((
unsigned int)pMinPts > this->
GetSize());
196 std::vector<SBPoint*>::iterator itPt = pPtsToCheck->begin();
197 while(itPt != pPtsToCheck->end() )
203 if( (!(*itPt)->HasCluster())
209 if(this->
GetSize() >= (
unsigned int)pMinPts)
214 itPt = pPtsToCheck->begin();
227 std::set<SBPoint*>::iterator itClusPt;
233 if((*pPtToCheck != *(*itClusPt)))
260 if(((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)+(z1-z2)*(z1-z2))<=
261 (pMinDist/
nm*pMinDist/
nm))
278 std::set<SBPoint*>::iterator itPt;
279 for(itPt = points.begin(); itPt != points.end(); ++itPt)
ClusterSBPoints(std::set< SBPoint * > pPoints)
std::set< SBPoint * > GetRegistredSBPoints() const
CLHEP::Hep3Vector G4ThreeVector
void FindAllPointsPossible(std::vector< SBPoint * > *ptToCheck, G4int minPts, G4double minDist)
void MergeWith(ClusterSBPoints *)
Float_t y1[n_points_granero]
bool AreOnTheSameCluster(const SBPoint *pPt1, const SBPoint *pPt2, G4double pMinDist)
static constexpr double nm
Definition of the ClusterSBPoints class.
Float_t x1[n_points_granero]
void UpdateDoubleStrand()
bool HasIn(const SBPoint *, G4double)
Float_t y2[n_points_geant4]
G4ThreeVector GetPosition() const
G4ThreeVector GetBarycenter() const
define a cluster of SB Points
void SetCluster(ClusterSBPoints *pCluster)
std::set< SBPoint * > fpRegisteredSBPoints
bool HasInBarycenter(const SBPoint *, G4double)
defines a point of energy deposition which defines a damage to the DNA.
Float_t x2[n_points_geant4]
void AddSBPoint(SBPoint *pSBPoint)
unsigned int GetSize() const