CexmcRunAction Class Reference

#include <CexmcRunAction.hh>

Inheritance diagram for CexmcRunAction:

Public Member Functions
	CexmcRunAction (CexmcPhysicsManager *physicsManager)
G4Run *	GenerateRun (void)
void	EndOfRunAction (const G4Run *run)
Public Member Functions inherited from G4UserRunAction
	G4UserRunAction ()
virtual	~G4UserRunAction ()
virtual void	BeginOfRunAction (const G4Run *aRun)
void	SetMaster (G4bool val=true)
G4bool	IsMaster () const

Static Public Member Functions
static void	PrintResults (const CexmcNmbOfHitsInRanges &nmbOfHitsSampled, const CexmcNmbOfHitsInRanges &nmbOfHitsSampledFull, const CexmcNmbOfHitsInRanges &nmbOfHitsTriggeredRealRange, const CexmcNmbOfHitsInRanges &nmbOfHitsTriggeredRecRange, const CexmcNmbOfHitsInRanges &nmbOfOrphanHits, const CexmcAngularRangeList &angularRanges, G4int nmbOfFalseHitsTriggeredEDT, G4int nmbOfFalseHitsTriggeredRec)

Additional Inherited Members
Protected Attributes inherited from G4UserRunAction
G4bool	isMaster

Detailed Description

Definition at line 54 of file CexmcRunAction.hh.

Constructor & Destructor Documentation

CexmcRunAction::CexmcRunAction ( CexmcPhysicsManager *  physicsManager )

explicit

Definition at line 56 of file CexmcRunAction.cc.

                                                                      :
    physicsManager( physicsManager )
{
}

Member Function Documentation

void CexmcRunAction::EndOfRunAction ( const G4Run *  run )

virtual

Reimplemented from G4UserRunAction.

Definition at line 230 of file CexmcRunAction.cc.

References CexmcWeirdException, G4cout, G4endl, CexmcRun::GetNmbOfFalseHitsTriggeredEDT(), CexmcRun::GetNmbOfFalseHitsTriggeredRec(), CexmcRun::GetNmbOfHitsSampled(), CexmcRun::GetNmbOfHitsSampledFull(), CexmcRun::GetNmbOfHitsTriggeredRealRange(), CexmcRun::GetNmbOfHitsTriggeredRecRange(), CexmcRun::GetNmbOfOrphanHits(), CexmcPhysicsManager::GetProductionModel(), and PrintResults().

{
    const CexmcRun *  theRun( static_cast< const CexmcRun * >( run ) );

    const CexmcNmbOfHitsInRanges &  nmbOfHitsSampled(
                                        theRun->GetNmbOfHitsSampled() );
    const CexmcNmbOfHitsInRanges &  nmbOfHitsSampledFull(
                                        theRun->GetNmbOfHitsSampledFull() );
    const CexmcNmbOfHitsInRanges &  nmbOfHitsTriggeredRealRange(
                                    theRun->GetNmbOfHitsTriggeredRealRange() );
    const CexmcNmbOfHitsInRanges &  nmbOfHitsTriggeredRecRange(
                                    theRun->GetNmbOfHitsTriggeredRecRange() );
    const CexmcNmbOfHitsInRanges &  nmbOfOrphanHits(
                                        theRun->GetNmbOfOrphanHits() );

    CexmcProductionModel *          productionModel(
                                        physicsManager->GetProductionModel() );
    if ( ! productionModel )
            throw CexmcException( CexmcWeirdException );

    const CexmcAngularRangeList &  angularRanges(
                                        productionModel->GetAngularRanges() );

    G4cout << G4endl;
    PrintResults( nmbOfHitsSampled, nmbOfHitsSampledFull,
                  nmbOfHitsTriggeredRealRange, nmbOfHitsTriggeredRecRange,
                  nmbOfOrphanHits, angularRanges,
                  theRun->GetNmbOfFalseHitsTriggeredEDT(),
                  theRun->GetNmbOfFalseHitsTriggeredRec() );
    G4cout << G4endl;
}

G4Run * CexmcRunAction::GenerateRun ( void  )

virtual

Reimplemented from G4UserRunAction.

Definition at line 62 of file CexmcRunAction.cc.

{
    return new CexmcRun;
}

void CexmcRunAction::PrintResults ( const CexmcNmbOfHitsInRanges &  nmbOfHitsSampled, const CexmcNmbOfHitsInRanges &  nmbOfHitsSampledFull, const CexmcNmbOfHitsInRanges &  nmbOfHitsTriggeredRealRange, const CexmcNmbOfHitsInRanges &  nmbOfHitsTriggeredRecRange, const CexmcNmbOfHitsInRanges &  nmbOfOrphanHits, const CexmcAngularRangeList &  angularRanges, G4int  nmbOfFalseHitsTriggeredEDT, G4int  nmbOfFalseHitsTriggeredRec  )

static

Definition at line 68 of file CexmcRunAction.cc.

References G4cout, G4endl, GetAngularGaps(), and G4INCL::CrossSections::total().

Referenced by EndOfRunAction().

{
    /* there are 7 auxiliary columns:
     * 1. acc real, [floating point number]
     * 2. triggered real,
     * 3. total hits sampled and monitored,
     * 4. acc reconstructed, [floating point number]
     * 5. triggered reconstructed,
     * 6. total hits sampled and monitored (identical with #3),
     * 7. total hits sampled.
     * As far as #3 and #6 are identical, nmbOfAuxColumns = 6 */
    const size_t                               nmbOfAuxColumns( 6 );
    const std::streamsize                      prec( 8 );
    std::vector< std::vector< std::string > >  auxStrings;
    size_t                                     maxSize[ nmbOfAuxColumns ];

    for ( size_t  i( 0 ); i < nmbOfAuxColumns; ++i )
        maxSize[ i ] = 0;

    /* addition of 2 (for '0.') for acceptances, that are floating point
     * numbers, is correct as far as ios::fixed will be used, and no negative
     * values are expected, and values will be less than 1. */
    maxSize[ 0 ] = prec + 2;
    maxSize[ 3 ] = prec + 2;

    for ( CexmcAngularRangeList::const_iterator  k( angularRanges.begin() );
                                                k != angularRanges.end(); ++k )
    {
        G4int     total( 0 );
        G4int     totalFull( 0 );
        G4int     triggered( 0 );
        G4double  acc( std::numeric_limits< G4double >::quiet_NaN() );

        CexmcNmbOfHitsInRanges::const_iterator  found(
                                        nmbOfHitsSampled.find( k->index ) );
        if ( found != nmbOfHitsSampled.end() )
        {
            total = found->second;
            acc = 0;
        }

        found = nmbOfHitsSampledFull.find( k->index );
        if ( found != nmbOfHitsSampledFull.end() )
        {
            totalFull = found->second;
        }

        G4double  accSave( acc );
        found = nmbOfHitsTriggeredRealRange.find( k->index );
        if ( found != nmbOfHitsTriggeredRealRange.end() )
        {
            triggered = found->second;
            if ( total > 0 )
                acc = G4double( triggered ) / total;
        }

        std::ostringstream  auxStringStream[ nmbOfAuxColumns ];

        for ( size_t  i( 0 ); i < nmbOfAuxColumns; ++i )
        {
            auxStringStream[ i ].precision( prec );
            auxStringStream[ i ].flags( std::ios::fixed );
        }

        G4int  i( 0 );

        auxStringStream[ i ] << acc;
        auxStringStream[ ++i ] << triggered;
        size_t  size( auxStringStream[ i ].str().size() );
        maxSize[ i ] = maxSize[ i ] > size ? maxSize[ i ] : size;
        auxStringStream[ ++i ] << total;
        size = auxStringStream[ i ].str().size();
        maxSize[ i ] = maxSize[ i ] > size ? maxSize[ i ] : size;

        triggered = 0;
        acc = accSave;
        found = nmbOfHitsTriggeredRecRange.find( k->index );
        if ( found != nmbOfHitsTriggeredRecRange.end() )
        {
            triggered = found->second;
            if ( total > 0 )
                acc = G4double( triggered ) / total;
        }

        auxStringStream[ ++i ] << acc;
        auxStringStream[ ++i ] << triggered;
        size = auxStringStream[ i ].str().size();
        maxSize[ i ] = maxSize[ i ] > size ? maxSize[ i ] : size;
        auxStringStream[ ++i ] << totalFull;
        size = auxStringStream[ i ].str().size();
        maxSize[ i ] = maxSize[ i ] > size ? maxSize[ i ] : size;

        std::vector< std::string >  auxString( nmbOfAuxColumns );

        for ( size_t  i( 0 ); i < nmbOfAuxColumns; ++i )
            auxString[ i ] = auxStringStream[ i ].str();

        auxStrings.push_back( auxString );
    }

    G4cout << " --- Setup acceptances (range | real (trg / mon) | "
              "rec (trg / mon / all)):" << G4endl;

    G4int  i( 0 );
    for ( CexmcAngularRangeList::const_iterator  k( angularRanges.begin() );
                                                k != angularRanges.end(); ++k )
    {
        G4cout << "       " << *k;
        G4int  j( 0 );
        G4cout << "  | " << std::setw( maxSize[ j ] );
        G4cout << auxStrings[ i ][ j++ ];
        G4cout << " ( " << std::setw( maxSize[ j ] );
        G4cout << auxStrings[ i ][ j++ ];
        G4cout << " / " << std::setw( maxSize[ j ] );
        G4cout << auxStrings[ i ][ j++ ];
        G4cout << " )  | " << std::setw( maxSize[ j ] );
        G4cout << auxStrings[ i ][ j++ ];
        G4cout << " ( " << std::setw( maxSize[ j ] );
        G4cout << auxStrings[ i ][ j++ ];
        G4cout << " / " << std::setw( maxSize[ 2 ] );
        G4cout << auxStrings[ i ][ 2 ];
        G4cout << " / " << std::setw( maxSize[ j ] );
        G4cout << auxStrings[ i++ ][ j++ ] << " )" << G4endl;
    }

    CexmcAngularRangeList  angularGaps;
    GetAngularGaps( angularRanges, angularGaps );

    if ( ! angularGaps.empty() )
    {
        G4cout << "    orphans detected: " << G4endl;
        for ( CexmcAngularRangeList::const_iterator  k( angularGaps.begin() );
                                                k != angularGaps.end(); ++k )
        {
            G4cout << "       " << *k;
            G4int     total( 0 );

            CexmcNmbOfHitsInRanges::const_iterator  found(
                                            nmbOfOrphanHits.find( k->index ) );
            if ( found != nmbOfHitsSampled.end() )
            {
                total = found->second;
            }
            G4cout << " " << total << G4endl;
        }
    }

    G4cout << "       ---" << G4endl;
    G4cout << "       False hits (edt | rec):  " <<
        nmbOfFalseHitsTriggeredEDT << " | " << nmbOfFalseHitsTriggeredRec <<
        G4endl;
}

---

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

• CexmcRunAction.hh
• CexmcRunAction.cc