Changeset 268 in Idarraga

Timestamp:

Mar 8, 2012, 3:10:30 PM (12 years ago)

Author:

idarraga

Message:

check hit collection valid

File:

• allpix/src/AllPixTrackerSD.cc (modified) (7 diffs)

allpix/src/AllPixTrackerSD.cc

@@ -75 +75 @@
 
 /*
- * Second constructor for sensitive devices wich are
+ * Second constructor for sensitive devices which are
  * not actual pixel detectors
  */
-AllPixTrackerSD::AllPixTrackerSD(G4String name,
-                G4ThreeVector absPos,
-                G4RotationMatrix * rot)
+AllPixTrackerSD::AllPixTrackerSD(G4String name, G4ThreeVector absPos, G4RotationMatrix * rot)
 :G4VSensitiveDetector(name)
 {
@@ -103 +101 @@
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
+/**
+ * Runs once per event.  The hitsCollections pointers are retrieved
+ *      I have a set of pointers to control that I get them right
+ */
 void AllPixTrackerSD::Initialize(G4HCofThisEvent* HCE)
 {
 
+        // Create the hit collection
         hitsCollection = new AllPixTrackerHitsCollection
                         (SensitiveDetectorName, collectionName[0]);
 
-        //G4cout << "TrackerSD : " << SensitiveDetectorName << " " << collectionName[0] << G4endl;
-
         static G4int HCID = -1;
         HCID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
 
+        // Add to hits collection of this event
         HCE->AddHitsCollection( HCID, hitsCollection );
 
+        // Insert the pointer in a set to check its existence later
+        // Normally there is only one instance of AllPixTrackerSD per sensitive volume
+        m_hitsCollectionSet.insert(hitsCollection);
+
 }
 
@@ -124 +130 @@
 {
 
+        // Check first if we are working in a valid HitCollection
+        try {
+                if( m_hitsCollectionSet.find(hitsCollection) == m_hitsCollectionSet.end()) // not found !
+                        throw hitsCollection;
+        } catch (AllPixTrackerHitsCollection * h) {
+                G4cout << "[WARNING] The following pointer to hitsCollection is invalid : " << h << G4endl;
+                G4cout << "          The available set contains " << m_hitsCollectionSet.size() << " collections." << G4endl;
+                G4cout << "          Trying to recover by ignoring this hit where this hitCollection" << G4endl;
+                G4cout << "          has been given. AllPixTrackerSD::ProcessHits returns false." << G4endl;
+                return false;
+        }
+
+        // Work with the Hit
         G4double edep = aStep->GetTotalEnergyDeposit();
         if(edep==0.) return false;
@@ -139 +158 @@
         G4ThreeVector correctedPos(0,0,0);
 
-        if(m_thisIsAPixelDetector){
+        if (m_thisIsAPixelDetector) {
                 // This positions are global, I will bring them to pixel-centered frame
                 // I can use the physical volumes for that
@@ -174 +193 @@
                 copyIDx_pre  = touchablepre->GetCopyNumber(1);
         }
+
         // Look at the touchablepost only if in the same volume, i.e. in the sensitive Si Box
         // If the hit is in a different pixel, it is still the same phys volume
@@ -222 +242 @@
         newHit->SetParentVolumeName(aTrack->GetLogicalVolumeAtVertex()->GetName());
 
+        //G4cout << "hitsCollection : " << hitsCollection << G4endl;
+        //G4cout << "     entries --> " << hitsCollection->entries() << G4endl;
         hitsCollection->insert(newHit);
-
         //newHit->Print();
         //newHit->Draw();
@@ -238 +259 @@
         if(NbHits > 0)
                 G4cout << "--------> Hits Collection : " << collectionName[0] << " has " << NbHits << " hits " << G4endl;
+
+        // clear the Set of pointers to hitCollection used for verification
+        m_hitsCollectionSet.clear();
 
 }