DicomRun Class Reference

#include <DicomRun.hh>

Inheritance diagram for DicomRun:

Collaboration diagram for DicomRun:

Public Member Functions
	DicomRun ()
	DicomRun (const std::vector< G4String > mfdName)
virtual	~DicomRun ()
virtual void	RecordEvent (const G4Event *)
G4int	GetNumberOfHitsMap () const
G4THitsMap< G4double > *	GetHitsMap (G4int i) const
G4THitsMap< G4double > *	GetHitsMap (const G4String &detName, const G4String &colName) const
G4THitsMap< G4double > *	GetHitsMap (const G4String &fullName) const
void	ConstructMFD (const std::vector< G4String > &)
virtual void	Merge (const G4Run *)
Public Member Functions inherited from G4Run
	G4Run ()
virtual	~G4Run ()
G4int	GetRunID () const
G4int	GetNumberOfEvent () const
G4int	GetNumberOfEventToBeProcessed () const
const G4HCtable *	GetHCtable () const
const G4DCtable *	GetDCtable () const
const G4String &	GetRandomNumberStatus () const
void	SetRunID (G4int id)
void	SetNumberOfEventToBeProcessed (G4int n_ev)
void	SetHCtable (G4HCtable *HCtbl)
void	SetDCtable (G4DCtable *DCtbl)
void	SetRandomNumberStatus (G4String &st)
void	StoreEvent (G4Event *evt)
const std::vector< const G4Event * > *	GetEventVector () const

Additional Inherited Members
Protected Attributes inherited from G4Run
G4int	runID
G4int	numberOfEvent
G4int	numberOfEventToBeProcessed
G4HCtable *	HCtable
G4DCtable *	DCtable
G4String	randomNumberStatus
std::vector< const G4Event * > *	eventVector

Detailed Description

DicomRun class

Example implementation for multi-functional-detector and primitive scorer. This DicomRun class has collections which accumulate a event information into a run information.

Definition at line 50 of file DicomRun.hh.

Constructor & Destructor Documentation

DicomRun::DicomRun

(

)

(Description) DicomRun Class is for accumulating scored quantities which is scored using G4MutiFunctionalDetector and G4VPrimitiveScorer. Accumulation is done using G4THitsMap object.

The constructor DicomRun(const std::vector<G4String> mfdName) needs a vector filled with MultiFunctionalDetector names which was assigned at instantiation of MultiFunctionalDetector(MFD). Then DicomRun constructor automatically scans primitive scorers in the MFD, and obtains collectionIDs of all collections associated to those primitive scorers. Futhermore, the G4THitsMap objects for accumulating during a RUN are automatically created too. (*) Collection Name is same as primitive scorer name.

The resultant information is kept inside DicomRun objects as data members. std::vector<G4String> fCollName; // Collection Name, std::vector<G4int> fCollID; // Collection ID, std::vector<G4THitsMap<G4double>*> fRunMap; // HitsMap for RUN.

The resualtant HitsMap objects are obtain using access method, GetHitsMap(..).

Definition at line 67 of file DicomRun.cc.

: G4Run()
{ }

DicomRun::DicomRun

(

const std::vector< G4String >

mfdName

)

Definition at line 75 of file DicomRun.cc.

                                                   : G4Run()
{
    ConstructMFD(mfdName);
}

Here is the call graph for this function:

DicomRun::~DicomRun ( )

virtual

Definition at line 84 of file DicomRun.cc.

{
  //--- Clear HitsMap for RUN
  G4int Nmap = fRunMap.size();
  for ( G4int i = 0; i < Nmap; i++){
    if(fRunMap[i] ) fRunMap[i]->clear();
  }
  fCollName.clear();
  fCollID.clear();
  fRunMap.clear();
}

Member Function Documentation

void DicomRun::ConstructMFD

(

const std::vector< G4String > &

mfdName

)

Definition at line 100 of file DicomRun.cc.

{

  G4SDManager* SDman = G4SDManager::GetSDMpointer();
  //=================================================
  //  Initalize RunMaps for accumulation.
  //  Get CollectionIDs for HitCollections.
  //=================================================
  G4int Nmfd = mfdName.size();
  for ( G4int idet = 0; idet < Nmfd ; idet++){  // Loop for all MFD.
    G4String detName = mfdName[idet];
    //--- Seek and Obtain MFD objects from SDmanager.
    G4MultiFunctionalDetector* mfd =
      (G4MultiFunctionalDetector*)(SDman->FindSensitiveDetector(detName));
    //
    if ( mfd ){
      //--- Loop over the registered primitive scorers.
      for (G4int icol = 0; icol < mfd->GetNumberOfPrimitives(); icol++){
        // Get Primitive Scorer object.
        G4VPrimitiveScorer* scorer = mfd->GetPrimitive(icol);
        // collection name and collectionID for HitsCollection,
        // where type of HitsCollection is G4THitsMap in case 
        // of primitive scorer.
        // The collection name is given by <MFD name>/<Primitive 
        // Scorer name>.
        G4String collectionName = scorer->GetName();
        G4String fullCollectionName = detName+"/"+collectionName;
        G4int    collectionID = SDman->GetCollectionID(fullCollectionName);
        //
        if ( collectionID >= 0 ){
          G4cout << "++ "<<fullCollectionName<< " id " << collectionID 
                 << G4endl;
          // Store obtained HitsCollection information into data members.
          // And, creates new G4THitsMap for accumulating quantities during RUN.
          fCollName.push_back(fullCollectionName);
          fCollID.push_back(collectionID);
          fRunMap.push_back(new G4THitsMap<G4double>(detName,collectionName));
        } else {
          G4cout << "** collection " << fullCollectionName << " not found. "
                 <<G4endl;
        }
      }
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4THitsMap<G4double>* DicomRun::GetHitsMap ( G4int  i ) const

inline

Definition at line 71 of file DicomRun.hh.

{return fRunMap[i];}

Here is the caller graph for this function:

G4THitsMap< G4double > * DicomRun::GetHitsMap	(	const G4String &	detName,
		const G4String &	colName
	)		const

Definition at line 210 of file DicomRun.cc.

{
  G4String fullName = detName+"/"+colName;
  return GetHitsMap(fullName);
}

Here is the call graph for this function:

G4THitsMap< G4double > * DicomRun::GetHitsMap

(

const G4String &

fullName

)

const

Definition at line 221 of file DicomRun.cc.

{
  
  //G4THitsMap<G4double>* hitsmap = 0;
  
  G4int Ncol = fCollName.size();
  for ( G4int i = 0; i < Ncol; i++){
    if ( fCollName[i] == fullName ){
      return fRunMap[i];
            //if(hitsmap) { *hitsmap += *fRunMap[i]; }
            //if(!hitsmap) { hitsmap = fRunMap[i]; }
    }
  }
  
  //if(hitsmap) { return hitsmap; }
  
  G4Exception("DicomRun", fullName.c_str(), JustWarning,
              "GetHitsMap failed to locate the requested HitsMap");
  return NULL;
}

Here is the call graph for this function:

G4int DicomRun::GetNumberOfHitsMap ( ) const

inline

Definition at line 67 of file DicomRun.hh.

{return fRunMap.size();}

void DicomRun::Merge ( const G4Run *  aRun )

virtual

Reimplemented from G4Run.

Definition at line 188 of file DicomRun.cc.

{
  const DicomRun* localRun = static_cast<const DicomRun*>(aRun);
  Copy(fCollName, localRun->fCollName);
  Copy(fCollID, localRun->fCollID);
  unsigned ncopies = Copy(fRunMap, localRun->fRunMap);
  // copy function returns the fRunMap size if all data is copied
  // so this loop isn't executed the first time around
  G4cout << "DicomRun :: Num copies = " << ncopies << G4endl;
  for(unsigned i = ncopies; i < fRunMap.size(); ++i) {
    *fRunMap[i] += *localRun->fRunMap[i];
  }
  G4Run::Merge(aRun);
}

Here is the call graph for this function:

void DicomRun::RecordEvent ( const G4Event *  aEvent )

virtual

Reimplemented from G4Run.

Definition at line 151 of file DicomRun.cc.

{
  G4Run::RecordEvent(aEvent);
  
  //G4cout << "Dicom Run :: Recording event " << aEvent->GetEventID() 
  //<< "..." << G4endl;
  //=============================
  // HitsCollection of This Event
  //============================
  G4HCofThisEvent* HCE = aEvent->GetHCofThisEvent();
  if (!HCE) return;
  
  //=======================================================
  // Sum up HitsMap of this Event  into HitsMap of this RUN
  //=======================================================
  G4int Ncol = fCollID.size();
  for ( G4int i = 0; i < Ncol ; i++ ){  // Loop over HitsCollection
    G4THitsMap<G4double>* EvtMap = 0;
    if ( fCollID[i] >= 0 ){           // Collection is attached to HCE
      EvtMap = static_cast<G4THitsMap<G4double>*>(HCE->GetHC(fCollID[i]));
    }else{
      G4cout <<" Error EvtMap Not Found "<< i << G4endl;
    }
    if ( EvtMap )  {
      //=== Sum up HitsMap of this event to HitsMap of RUN.===
      *fRunMap[i] += *EvtMap;
      //G4cout << "Summing EvtMap into RunMap at " << i << "..." << G4endl;
      //======================================================
    } //else {
    //G4cout << "Null pointer to EvtMap at " << i << "..." << G4endl;
    //}
  }
    
}

Here is the call graph for this function:

---

The documentation for this class was generated from the following files:

• source/geant4.10.03.p02/examples/extended/medical/DICOM/include/DicomRun.hh
• source/geant4.10.03.p02/examples/extended/medical/DICOM/src/DicomRun.cc