TSRunAction Class Reference

#include <TSRunAction.hh>

Inheritance diagram for TSRunAction:

Collaboration diagram for TSRunAction:

Public Types
typedef std::tuple< G4double, G4double, G4double >	Compare_t
typedef std::map< G4int, Compare_t >	IDcompare_t
typedef std::map< G4String, IDcompare_t >	TypeCompare_t

Public Member Functions
	TSRunAction ()
virtual	~TSRunAction ()
virtual void	BeginOfRunAction (const G4Run *)
virtual void	EndOfRunAction (const G4Run *)
virtual G4Run *	GenerateRun ()
Public Member Functions inherited from G4UserRunAction
	G4UserRunAction ()
virtual	~G4UserRunAction ()
void	SetMaster (G4bool val=true)
G4bool	IsMaster () const

Private Attributes
TSDetectorConstruction *	fDetector
G4String	fName
TypeCompare_t	fTypeCompare

Additional Inherited Members
Protected Attributes inherited from G4UserRunAction
G4bool	isMaster

Detailed Description

Definition at line 52 of file TSRunAction.hh.

Member Typedef Documentation

Compare_t

typedef std::tuple<G4double, G4double, G4double> TSRunAction::Compare_t

Definition at line 55 of file TSRunAction.hh.

IDcompare_t

typedef std::map<G4int, Compare_t> TSRunAction::IDcompare_t

Definition at line 56 of file TSRunAction.hh.

TypeCompare_t

typedef std::map<G4String, IDcompare_t> TSRunAction::TypeCompare_t

Definition at line 57 of file TSRunAction.hh.

Constructor & Destructor Documentation

TSRunAction()

TSRunAction::TSRunAction

(

)

Definition at line 51 of file TSRunAction.cc.

: fDetector(TSDetectorConstruction::Instance()),
  fName(fDetector->GetMFDName())
{ }

~TSRunAction()

TSRunAction::~TSRunAction ( )

virtual

Definition at line 58 of file TSRunAction.cc.

{ }

Member Function Documentation

BeginOfRunAction()

void TSRunAction::BeginOfRunAction ( const G4Run *  aRun )

virtual

Reimplemented from G4UserRunAction.

Definition at line 70 of file TSRunAction.cc.

{
    G4int evts_to_process = aRun->GetNumberOfEventToBeProcessed();
    G4RunManager::GetRunManager()->SetPrintProgress((evts_to_process > 100)
                                                    ? evts_to_process/100
                                                    : 1);
}

Here is the call graph for this function:

EndOfRunAction()

void TSRunAction::EndOfRunAction ( const G4Run *  aRun )

virtual

Reimplemented from G4UserRunAction.

Definition at line 80 of file TSRunAction.cc.

{

    if(IsMaster())
    {
        G4cout << " ###### EndOfTSRunAction ###### " << G4endl;

        aRun->GetNumberOfEvent();
        std::ofstream fileout;
        G4String fname = "";
        std::stringstream separator;

        separator
        << "============================================================";

        typedef std::set<G4int> IDSet_t;
        IDSet_t IDs;

        //- TSRun object.
        const TSRun* tsRun = static_cast<const TSRun*>(aRun);
        //--- Dump all scored quantities involved in TSRun.

        //---------------------------------------------
        // Dump accumulated quantities for this RUN.
        //---------------------------------------------
        std::vector<G4String> primScorerNames { "EnergyDeposit",
                                                "NumberOfSteps" };
        std::vector<G4String> fnames { "mfd_tl", "mfd_tg" };
        std::vector<G4double> units { CLHEP::eV, CLHEP::keV, 1, 1 };
        std::vector<G4String> unitstr { "keV", "steps" };

        for(unsigned i = 0; i < primScorerNames.size(); ++i)
        {
          for(unsigned j = 0; j < fnames.size(); ++j)
          {
            //----------------------------------------------------------------//
            auto print = [] (std::ostream& fout,
                G4int first, G4double second,
                G4double unit1, G4double unit2, G4String unit2str)
            {
              if(fout)
                  fout <<  first
                  << "     "  << second/unit1
                  << G4endl;

              G4cout
              << "    " << std::setw(10) << first
              << "    " << std::setw(15) << std::setprecision(6)
              << std::fixed << second/unit2 << " " << unit2str
              << G4endl;
              G4cout.unsetf(std::ios::fixed);
            };
            //----------------------------------------------------------------//

            fname = fnames.at(j) + "_" + primScorerNames.at(i) + ".out";
            fileout.open(fname);
            G4cout << separator.str() << G4endl;
            G4cout << " opened file " << fname << " for output" << G4endl;
            G4cout << separator.str() << G4endl;

            G4bool valid = true;
            if(j == 0)
            {
              G4THitsMap<G4double>* hitmap
                  = tsRun->GetHitsMap(fName + "/" + primScorerNames.at(i));
              if(hitmap && hitmap->GetMap()->size() != 0)
              {
                for(auto itr = hitmap->GetMap()->begin();
                    itr != hitmap->GetMap()->end(); itr++)
                {
                  IDs.insert(itr->first);
                  std::get<0>(fTypeCompare[primScorerNames.at(i)][itr->first])
                      = *itr->second/units.at(i);
                  print(fileout, itr->first, *itr->second,
                        units.at(i), units.at(i+1), unitstr.at(i));
                }
              } else {
                valid = false;
              }
            } else {
              G4TAtomicHitsMap<G4double>* hitmap
                  = tsRun->GetAtomicHitsMap(fName + "/" +
                                              primScorerNames.at(i));
              if(hitmap && hitmap->GetMap()->size() != 0)
              {
                for(auto itr = hitmap->GetMap()->begin();
                    itr != hitmap->GetMap()->end(); itr++)
                {
                  IDs.insert(itr->first);
                  std::get<1>(fTypeCompare[primScorerNames.at(i)][itr->first])
                      = *itr->second/units.at(i);
                  print(fileout, itr->first, *itr->second,
                        units.at(i), units.at(i+1), unitstr.at(i));
                }
              } else {
                valid = false;
              }
              if(!valid)
              {
                G4Exception("TSRunAction", "000", JustWarning,
                            G4String(primScorerNames.at(i) +
                            " HitsMap is either not "
                            "created or the HitsMap was empty").c_str());
              }
            }

            fileout.close();
            G4cout << separator.str() << G4endl;
            G4cout << " closed file " << fname << " for output" << G4endl;
          }
          // add the mutex data
          TSRun::MutexHitsMap_t* hitmap
              = tsRun->GetMutexHitsMap(fName + "/" +
                                         primScorerNames.at(i));
          if(hitmap && hitmap->size() != 0)
          {
            for(auto itr = hitmap->begin();
                itr != hitmap->end(); itr++)
            {
              IDs.insert(itr->first);
              std::get<2>(fTypeCompare[primScorerNames.at(i)][itr->first])
                  = itr->second/units.at(i);
            }
          }

        }

        //--------------------------------------------------------------------//
        // Check that the values are equivalent and there are no
        // IDs in one container that aren't in another
        //--------------------------------------------------------------------//

        fname = "mfd_diff.out";
        fileout.open(fname);

        G4cout << separator.str() << G4endl;
        G4cout << " opened file " << fname << " for difference output" << G4endl;
        G4cout << separator.str() << G4endl;

        fileout
        << "    " << std::setw(10) << "ID"
        << "    "
        << std::setw(30) << std::setprecision(12) << std::fixed
        << "MFD value"
        << "    "
        << std::setw(30) << std::setprecision(12) << std::fixed
        << "Atomic Hits Map value"
        << "    "
        << std::setw(30) << std::setprecision(12) << std::fixed
        << "Difference"
        << "    "
        << std::setw(30) << std::setprecision(12) << std::fixed
        << "Diff (MFD - MUTEXED)"
        << "    "
        << std::setw(30) << std::setprecision(12) << std::fixed
        << "Diff (ATOM_HIT_MAP - MUTEXED)"
        << G4endl << G4endl;

        for(auto itr1 = fTypeCompare.begin();
            itr1 != fTypeCompare.end(); ++itr1)
        {
          fileout << "\n\nType = " << itr1->first << "\n" << G4endl;
          for(auto itr2 = itr1->second.begin();
              itr2 != itr1->second.end(); ++itr2)
          {
            fileout
            << "    " << std::setw(10) << itr2->first
            << "    "
            << std::setw(30) << std::setprecision(12) << std::fixed
            << std::get<0>(itr2->second)
            << "    "
            << std::setw(30) << std::setprecision(12) << std::fixed
            << std::get<1>(itr2->second)
            << "    "
            << std::setw(30) << std::setprecision(12) << std::fixed
            << (std::fabs(std::get<0>(itr2->second) - std::get<1>(itr2->second)))
            << "    "
            << std::setw(30) << std::setprecision(12) << std::fixed
            << (std::fabs(std::get<0>(itr2->second) - std::get<2>(itr2->second)))
            << "    "
            << std::setw(30) << std::setprecision(12) << std::fixed
            << (std::fabs(std::get<1>(itr2->second) - std::get<2>(itr2->second)))
            << G4endl;

          }
        }

        fileout.close();
        G4cout << " closed file " << fname << " for difference output" << G4endl;
        G4cout << separator.str() << G4endl;

    }


}

Here is the call graph for this function:

GenerateRun()

G4Run * TSRunAction::GenerateRun ( void  )

virtual

Reimplemented from G4UserRunAction.

Definition at line 63 of file TSRunAction.cc.

{
    return new TSRun(fName);
}

Member Data Documentation

fDetector

TSDetectorConstruction* TSRunAction::fDetector

private

Definition at line 69 of file TSRunAction.hh.

fName

G4String TSRunAction::fName

private

Definition at line 70 of file TSRunAction.hh.

fTypeCompare

TypeCompare_t TSRunAction::fTypeCompare

private

Definition at line 71 of file TSRunAction.hh.

---

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

• TSRunAction.hh
• TSRunAction.cc