Dicom2RunAction Class Reference

#include <Dicom2RunAction.hh>

Inheritance diagram for Dicom2RunAction:

Public Types
typedef Dicom2Run::Dicom2RunVector	Dicom2RunVector

Public Member Functions
	Dicom2RunAction ()
	Constructor.
virtual	~Dicom2RunAction ()
	Destructor.
virtual G4Run *	GenerateRun () override
virtual void	EndOfRunAction (const G4Run *aRun) override
Public Member Functions inherited from DicomRunAction
	DicomRunAction ()
	Constructor.
	~DicomRunAction ()
	Destructor.
G4Run *	GenerateRun ()
void	BeginOfRunAction (const G4Run *)
void	EndOfRunAction (const G4Run *)
DicomRun *	GetDicomRun () const
void	PrintHeader (std::ostream *out)
std::string	FillString (const std::string &name, char c, G4int n, G4bool back=true)

Additional Inherited Members
Protected Attributes inherited from DicomRunAction
DicomRun *	fDcmrun
std::vector< G4String >	fSDName

Detailed Description

Definition at line 42 of file Dicom2RunAction.hh.

Member Typedef Documentation

Dicom2RunVector

typedef Dicom2Run::Dicom2RunVector Dicom2RunAction::Dicom2RunVector

Definition at line 45 of file Dicom2RunAction.hh.

Constructor & Destructor Documentation

Dicom2RunAction()

Dicom2RunAction::Dicom2RunAction

(

)

Constructor.

Definition at line 47 of file Dicom2RunAction.cc.

: DicomRunAction()
{ }

~Dicom2RunAction()

Dicom2RunAction::~Dicom2RunAction ( )

virtual

Destructor.

Definition at line 53 of file Dicom2RunAction.cc.

{ }

Member Function Documentation

GenerateRun()

G4Run * Dicom2RunAction::GenerateRun ( )

overridevirtual

Definition at line 57 of file Dicom2RunAction.cc.

{
    // Generate new RUN object, which is specially
    // dedicated for MultiFunctionalDetector scheme.
    //  Detail description can be found in Dicom2Run.hh/cc.
    return fDcmrun = new Dicom2Run(fSDName);
}

EndOfRunAction()

void Dicom2RunAction::EndOfRunAction ( const G4Run *  aRun )

overridevirtual

Definition at line 66 of file Dicom2RunAction.cc.

{
    // Lock the output because of the external calls to DicomRunAction
    // otherwise, the output gets too confusing
    G4AutoLock l(G4TypeMutex<Dicom2RunAction>());
 
    G4cout << G4endl;
    G4cout << "[==========================================================="
           << " DICOM "
           << "===========================================================]"
           << G4endl;
    G4cout << G4endl;
 
    DicomRunAction::EndOfRunAction(aRun);
 
    G4cout << G4endl;
    G4cout << "[==========================================================="
           << " DICOM2 "
           << "==========================================================]"
           << G4endl;
    G4cout << G4endl;
 
    G4int nofEvents = aRun->GetNumberOfEvent();
 
    G4StatAnalysis local_total_dose;
 
    const Dicom2Run* dcm2Run = static_cast<const Dicom2Run*>(aRun);
    //--- Dump all scored quantities involved in Dicom2Run.
    for(uintmax_t i = 0; i < fSDName.size(); i++)
    {
        //
        //---------------------------------------------
        // Dump accumulated quantities for this RUN.
        //  (Display only central region of x-y plane)
        //      0       ConcreteSD/DoseDeposit
        //---------------------------------------------
        Dicom2RunVector* DoseDeposit =
                dcm2Run->GetHitsVector(fSDName[i]+"/DoseDeposit");
 
        if(DoseDeposit && DoseDeposit->size() != 0 )
        {
            for(auto itr = DoseDeposit->begin(); itr != DoseDeposit->end(); ++itr)
            {
                // this will sometimes return null pointers
                if(!DoseDeposit->GetObject(itr))
                    continue;
                local_total_dose += (*DoseDeposit->GetObject(itr));
            }
        }
    }
 
    if(IsMaster())
    {
        G4cout << " ###### EndOfRunAction ###### " << G4endl;
        //- Dicom2Run object.
        const Dicom2Run* re02Run = static_cast<const Dicom2Run*>(aRun);
        //--- Dump all scored quantities involved in Dicom2Run.
 
        for(uintmax_t i = 0; i < fSDName.size(); i++)
        {
            //
            //---------------------------------------------
            // Dump accumulated quantities for this RUN.
            //  (Display only central region of x-y plane)
            //      0       ConcreteSD/DoseDeposit
            //---------------------------------------------
            Dicom2RunVector* DoseDeposit =
                    re02Run->GetHitsVector(fSDName[i]+"/DoseDeposit");
 
            G4cout << "============================================================="
                   <<G4endl;
            G4cout << " Number of event processed : "
                   << aRun->GetNumberOfEvent() << G4endl;
            G4cout << "============================================================="
                   <<G4endl;
 
            std::ofstream fileout;
            G4String fname = "dicom2-vector.out";
            fileout.open(fname);
            G4cout << " opened file " << fname << " for dose output" << G4endl;
 
            if(DoseDeposit && DoseDeposit->size() != 0)
            {
                std::ostream *myout = &G4cout;
                PrintHeader(myout);
                for(auto itr = DoseDeposit->begin(); itr != DoseDeposit->end();
                    ++itr)
                {
                    auto _idx = DoseDeposit->GetIndex(itr);
                    G4StatAnalysis* _stat = DoseDeposit->GetObject(itr);
                    if(_stat && _stat->GetHits() > 0)
                    {
                        G4StatAnalysis _tmp_stat = *_stat;
                        _tmp_stat /= CLHEP::gray;
                        fileout << _idx << "     "  << (*_stat) << G4endl;
                    }
                }
                G4cout << "============================================="<<G4endl;
            }
            else
            {
                G4Exception("Dicom2RunAction", "000", JustWarning,
                            "DoseDeposit HitsMap is either a null pointer "
                            "of the HitsMap was empty");
            }
            fileout.close();
            G4cout << " closed file " << fname << " for dose output" << G4endl;
 
        }
    }
 
 
    if (IsMaster())
    {
        // convert to units of Gy
        local_total_dose /= gray;
        G4cout << "--------------------End of Global Run-----------------------"
               << G4endl;
        G4cout << " The run was " << nofEvents << " events " << G4endl;
        G4cout << "      TOTAL DOSE : \t" << local_total_dose << " Gy" << G4endl;
        if(nofEvents > 0)
        {
            local_total_dose /= nofEvents;
            G4cout << " TOTAL DOSE/Bq-s : \t" << local_total_dose << " Gy/Bq-s"
                   << G4endl;
        }
    }
    else
    {
        // convert to units of Gy
        local_total_dose /= gray;
        G4cout << "--------------------End of Local Run------------------------"
               << G4endl;
        G4cout << " The run was " << nofEvents << " events" << G4endl;
        G4cout << "LOCAL TOTAL DOSE : \t" << local_total_dose << " Gy" << G4endl;
        if(nofEvents > 0)
        {
            local_total_dose /= nofEvents;
            G4cout << " LOCAL DOSE/Bq-s : \t" << local_total_dose << " Gy/Bq-s"
                   << G4endl;
        }
    }
 
    G4cout << G4endl;
    G4cout << "Finished : End of Run Action " << aRun->GetRunID() << "\n" << G4endl;
 
}

---

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

• extended/medical/DICOM2/include/Dicom2RunAction.hh
• extended/medical/DICOM2/src/Dicom2RunAction.cc