CexmcEventAction Class Reference

#include <CexmcEventAction.hh>

Inheritance diagram for CexmcEventAction:

Collaboration diagram for CexmcEventAction:

Public Member Functions
	CexmcEventAction (CexmcPhysicsManager *physicsManager, G4int verbose=0)
virtual	~CexmcEventAction ()
void	BeginOfEventAction (const G4Event *event)
void	EndOfEventAction (const G4Event *event)
void	BeamParticleChangeHook (void)
void	SetVerboseOnCexmcLevel (G4int value)
void	SetVerboseDrawLevel (G4int value)
CexmcChargeExchangeReconstructor *	GetReconstructor (void)
Public Member Functions inherited from G4UserEventAction
	G4UserEventAction ()
virtual	~G4UserEventAction ()
virtual void	SetEventManager (G4EventManager *value)

Static Public Member Functions
static CexmcEnergyDepositStore *	MakeEnergyDepositStore (const CexmcEnergyDepositDigitizer *digitizer)
static CexmcTrackPointsStore *	MakeTrackPointsStore (const CexmcTrackPointsDigitizer *digitizer)
static void	PrintEnergyDeposit (const CexmcEnergyDepositStore *edStore)
static void	PrintTrackPoints (const CexmcTrackPointsStore *tpStore)
static void	PrintProductionModelData (const CexmcAngularRangeList &angularRanges, const CexmcProductionModelData &pmData)

Private Member Functions
void	PrintReconstructedData (const CexmcAngularRangeList &angularRanges, const CexmcAngularRange &angularGap) const
void	DrawTrajectories (const G4Event *event)
void	DrawTrackPoints (const CexmcTrackPointsStore *tpStore) const
void	DrawReconstructionData (void)
void	UpdateRunHits (const CexmcAngularRangeList &aRangesReal, const CexmcAngularRangeList &aRangesRec, G4bool tpDigitizerHasTriggered, G4bool edDigitizerHasTriggered, G4bool edDigitizerMonitorHasTriggered, G4bool reconstructorHasTriggered, const CexmcAngularRange &aGap)

Private Attributes
CexmcPhysicsManager *	physicsManager
CexmcChargeExchangeReconstructor *	reconstructor
G4double	opKinEnergy
G4int	verbose
G4int	verboseDraw
CexmcEventActionMessenger *	messenger

Additional Inherited Members
Protected Attributes inherited from G4UserEventAction
G4EventManager *	fpEventManager

Detailed Description

Definition at line 61 of file CexmcEventAction.hh.

Constructor & Destructor Documentation

CexmcEventAction()

CexmcEventAction::CexmcEventAction ( CexmcPhysicsManager *  physicsManager, G4int  verbose = 0  )

explicit

Definition at line 94 of file CexmcEventAction.cc.

                                                     :
    physicsManager( physicsManager ), reconstructor( NULL ), opKinEnergy( 0. ),
    verbose( verbose ), verboseDraw( 4 ), messenger( NULL )
{
    G4DigiManager *  digiManager( G4DigiManager::GetDMpointer() );
    digiManager->AddNewModule( new CexmcEnergyDepositDigitizer(
                                                    CexmcEDDigitizerName ) );
    digiManager->AddNewModule( new CexmcTrackPointsDigitizer(
                                                    CexmcTPDigitizerName ) );
    reconstructor = new CexmcChargeExchangeReconstructor(
                                        physicsManager->GetProductionModel() );
    messenger = new CexmcEventActionMessenger( this );
}

Here is the call graph for this function:

~CexmcEventAction()

CexmcEventAction::~CexmcEventAction ( )

virtual

Definition at line 110 of file CexmcEventAction.cc.

{
    delete reconstructor;
    delete messenger;
}

Member Function Documentation

BeamParticleChangeHook()

void CexmcEventAction::BeamParticleChangeHook

(

void

)

Definition at line 117 of file CexmcEventAction.cc.

{
    reconstructor->SetupBeamParticle();
}

Here is the call graph for this function:

Here is the caller graph for this function:

BeginOfEventAction()

void CexmcEventAction::BeginOfEventAction ( const G4Event *  event )

virtual

Reimplemented from G4UserEventAction.

Definition at line 123 of file CexmcEventAction.cc.

{
    G4RunManager *         runManager( G4RunManager::GetRunManager() );
    CexmcTrackingAction *  trackingAction
            ( static_cast< CexmcTrackingAction * >(
                        const_cast< G4UserTrackingAction * >(
                                    runManager->GetUserTrackingAction() ) ) );
    trackingAction->BeginOfEventAction();

    physicsManager->ResetNumberOfTriggeredStudiedInteractions();
}

Here is the call graph for this function:

DrawReconstructionData()

void CexmcEventAction::DrawReconstructionData ( void  )

private

Definition at line 645 of file CexmcEventAction.cc.

{
    G4VisManager *  visManager( static_cast< G4VisManager * >(
                                    G4VVisManager::GetConcreteInstance() ) );
    if ( ! visManager || ! visManager->GetCurrentGraphicsSystem() )
        return;

    G4Circle  circle( reconstructor->GetTargetEPWorldPosition() );
    circle.SetScreenSize( CexmcSmallCircleScreenSize );
    circle.SetFillStyle( G4Circle::filled );
    G4VisAttributes  visAttributes( CexmcRecTrackPointsMarkerColour );
    circle.SetVisAttributes( visAttributes );
    visManager->Draw( circle );

    circle.SetScreenSize( CexmcBigCircleScreenSize );
    circle.SetPosition( reconstructor->GetCalorimeterEPLeftWorldPosition() );
    visManager->Draw( circle );

    circle.SetPosition( reconstructor->GetCalorimeterEPRightWorldPosition() );
    visManager->Draw( circle );
}

Here is the call graph for this function:

Here is the caller graph for this function:

DrawTrackPoints()

void CexmcEventAction::DrawTrackPoints ( const CexmcTrackPointsStore *  tpStore ) const

private

Definition at line 587 of file CexmcEventAction.cc.

{
    G4VisManager *  visManager( static_cast< G4VisManager * >(
                                    G4VVisManager::GetConcreteInstance() ) );
    if ( ! visManager || ! visManager->GetCurrentGraphicsSystem() )
        return;

    G4Circle         circle;
    G4VisAttributes  visAttributes( CexmcTrackPointsMarkerColour );
    circle.SetScreenSize( CexmcSmallCircleScreenSize );
    circle.SetFillStyle( G4Circle::filled );
    circle.SetVisAttributes( visAttributes );

    if ( tpStore->monitorTP.IsValid() )
    {
        circle.SetPosition( tpStore->monitorTP.positionWorld );
        visManager->Draw( circle );
    }

    if ( tpStore->targetTPBeamParticle.IsValid() )
    {
        circle.SetPosition( tpStore->targetTPBeamParticle.positionWorld );
        visManager->Draw( circle );
    }

    if ( tpStore->targetTPOutputParticle.IsValid() )
    {
        circle.SetPosition( tpStore->targetTPOutputParticle.positionWorld );
        visManager->Draw( circle );
    }

    if ( tpStore->vetoCounterTPLeft.IsValid() )
    {
        circle.SetPosition( tpStore->vetoCounterTPLeft.positionWorld );
        visManager->Draw( circle );
    }

    if ( tpStore->vetoCounterTPRight.IsValid() )
    {
        circle.SetPosition( tpStore->vetoCounterTPRight.positionWorld );
        visManager->Draw( circle );
    }

    if ( tpStore->calorimeterTPLeft.IsValid() )
    {
        circle.SetPosition( tpStore->calorimeterTPLeft.positionWorld );
        visManager->Draw( circle );
    }

    if ( tpStore->calorimeterTPRight.IsValid() )
    {
        circle.SetPosition( tpStore->calorimeterTPRight.positionWorld );
        visManager->Draw( circle );
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

DrawTrajectories()

void CexmcEventAction::DrawTrajectories ( const G4Event *  event )

private

Definition at line 564 of file CexmcEventAction.cc.

{
    G4VisManager *  visManager( static_cast< G4VisManager * >(
                                    G4VVisManager::GetConcreteInstance() ) );
    if ( ! visManager || ! visManager->GetCurrentGraphicsSystem() )
        return;

    G4int                    nTraj( 0 );
    G4TrajectoryContainer *  trajContainer( event->GetTrajectoryContainer() );

    if ( ! trajContainer )
        return;

    nTraj = trajContainer->entries();

    for ( int  i( 0 ); i < nTraj; ++i )
    {
        G4VTrajectory *  traj( ( *trajContainer )[ i ] );
        traj->DrawTrajectory();
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

EndOfEventAction()

void CexmcEventAction::EndOfEventAction ( const G4Event *  event )

virtual

Reimplemented from G4UserEventAction.

Definition at line 799 of file CexmcEventAction.cc.

{
    G4DigiManager *                digiManager( G4DigiManager::GetDMpointer() );
    CexmcEnergyDepositDigitizer *  energyDepositDigitizer(
            static_cast< CexmcEnergyDepositDigitizer * >( digiManager->
                                FindDigitizerModule( CexmcEDDigitizerName ) ) );
    CexmcTrackPointsDigitizer *    trackPointsDigitizer(
            static_cast< CexmcTrackPointsDigitizer * >( digiManager->
                                FindDigitizerModule( CexmcTPDigitizerName ) ) );

    energyDepositDigitizer->Digitize();
    trackPointsDigitizer->Digitize();

    G4bool  edDigitizerMonitorHasTriggered(
                                energyDepositDigitizer->MonitorHasTriggered() );
    G4bool  edDigitizerHasTriggered( false );

    CexmcEventInfo *  eventInfo( static_cast< CexmcEventInfo * >(
                                                event->GetUserInformation() ) );
    if ( ! eventInfo || eventInfo->EdTriggerIsOk() )
        edDigitizerHasTriggered = energyDepositDigitizer->HasTriggered();

    G4bool  tpDigitizerHasTriggered( trackPointsDigitizer->HasTriggered() );
    G4bool  reconstructorHasBasicTrigger( false );
    G4bool  reconstructorHasFullTrigger( false );

    CexmcEnergyDepositStore *  edStore( MakeEnergyDepositStore(
                                                    energyDepositDigitizer ) );
    CexmcTrackPointsStore *    tpStore( MakeTrackPointsStore(
                                                    trackPointsDigitizer ) );

    try
    {
        CexmcProductionModel *  productionModel(
                                        physicsManager->GetProductionModel() );

        if ( ! productionModel )
            throw CexmcException( CexmcWeirdException );

        const CexmcAngularRangeList &     angularRanges(
                                productionModel->GetAngularRanges() );
        const CexmcAngularRangeList &     triggeredAngularRanges(
                                productionModel->GetTriggeredAngularRanges() );
        const CexmcProductionModelData &  pmData(
                                productionModel->GetProductionModelData() );

        if ( edDigitizerHasTriggered )
        {
            reconstructor->Reconstruct( edStore );
            reconstructorHasBasicTrigger = reconstructor->HasBasicTrigger();
            reconstructorHasFullTrigger = reconstructor->HasFullTrigger();
        }

        CexmcAngularRangeList  triggeredRecAngularRanges;

        if ( reconstructorHasBasicTrigger )
        {
            for ( CexmcAngularRangeList::const_iterator
                  k( angularRanges.begin() ); k != angularRanges.end(); ++k )
            {
                G4double  cosTheta( reconstructor->GetProductionModelData().
                                    outputParticleSCM.cosTheta() );
                if ( cosTheta <= k->top && cosTheta > k->bottom )
                    triggeredRecAngularRanges.push_back( CexmcAngularRange(
                                                k->top, k->bottom, k->index ) );
            }
        }

        CexmcAngularRange  angularGap( 0.0, 0.0, 0 );
        if ( triggeredRecAngularRanges.empty() )
        {
            CexmcAngularRangeList  angularGaps;
            GetAngularGaps( angularRanges, angularGaps );
            for ( CexmcAngularRangeList::const_iterator
                    k( angularGaps.begin() ); k != angularGaps.end(); ++k )
            {
                G4double  cosTheta( reconstructor->GetProductionModelData().
                                    outputParticleSCM.cosTheta() );
                if ( cosTheta <= k->top && cosTheta > k->bottom )
                {
                    angularGap = *k;
                    break;
                }
            }
        }

        UpdateRunHits( triggeredAngularRanges, triggeredRecAngularRanges,
                       tpDigitizerHasTriggered, edDigitizerHasTriggered,
                       edDigitizerMonitorHasTriggered,
                       reconstructorHasFullTrigger, angularGap );

        if ( verbose > 0 )
        {
            G4bool  printMessages( verbose > 3 ||
                        ( ( verbose == 1 ) && tpDigitizerHasTriggered ) ||
                        ( ( verbose == 2 ) && edDigitizerHasTriggered ) ||
                        ( ( verbose == 3 ) && ( tpDigitizerHasTriggered ||
                                                edDigitizerHasTriggered ) ) );
            if ( printMessages )
            {
                G4cout << "Event " << event->GetEventID() << G4endl;
                if ( tpDigitizerHasTriggered )
                {
                    PrintTrackPoints( tpStore );
                    PrintProductionModelData( triggeredAngularRanges, pmData );
                }
                if ( reconstructorHasBasicTrigger )
                    PrintReconstructedData( triggeredRecAngularRanges,
                                            angularGap );
                if ( edDigitizerHasTriggered )
                    PrintEnergyDeposit( edStore );
            }
        }

        if ( verboseDraw > 0 )
        {
            G4bool  drawTrajectories( verboseDraw > 3 ||
                        ( ( verboseDraw == 1 ) && tpDigitizerHasTriggered ) ||
                        ( ( verboseDraw == 2 ) && edDigitizerHasTriggered ) ||
                        ( ( verboseDraw == 3 ) && ( tpDigitizerHasTriggered ||
                                                edDigitizerHasTriggered ) ) );
            if ( drawTrajectories )
            {
                DrawTrajectories( event );
                if ( tpDigitizerHasTriggered )
                    DrawTrackPoints( tpStore );
                if ( reconstructorHasBasicTrigger )
                    DrawReconstructionData();
            }
        }

#ifdef CEXMC_USE_PERSISTENCY
        if ( edDigitizerHasTriggered || tpDigitizerHasTriggered )
        {
            SaveEventFast( event, tpDigitizerHasTriggered,
                           edDigitizerHasTriggered,
                           edDigitizerMonitorHasTriggered,
                           pmData.outputParticleSCM.cosTheta() );
            SaveEvent( event, edDigitizerHasTriggered, edStore, tpStore,
                       pmData );
        }
#endif

#ifdef CEXMC_USE_ROOT
        /* opKinEnergy will be used in several histos */
        if ( tpStore->targetTPOutputParticle.IsValid() )
        {
            opKinEnergy = CexmcGetKinEnergy(
                    tpStore->targetTPOutputParticle.momentumAmp,
                    tpStore->targetTPOutputParticle.particle->GetPDGMass() );
        }

        if ( edDigitizerHasTriggered )
            FillEDTHistos( edStore, triggeredAngularRanges );

        /* fill TPT histos only when the monitor has triggered because events
         * when it was missed have less value for us */
        if ( tpDigitizerHasTriggered && edDigitizerMonitorHasTriggered )
            FillTPTHistos( tpStore, pmData, triggeredAngularRanges );

        if ( reconstructorHasBasicTrigger )
            FillRTHistos( reconstructorHasFullTrigger, edStore, tpStore,
                          pmData, triggeredAngularRanges );
#endif

        G4Event *  theEvent( const_cast< G4Event * >( event ) );
        if ( eventInfo )
        {
            delete eventInfo;
            theEvent->SetUserInformation( NULL );
        }
        theEvent->SetUserInformation( new CexmcEventInfo(
                                                edDigitizerHasTriggered,
                                                tpDigitizerHasTriggered,
                                                reconstructorHasFullTrigger ) );
    }
    catch ( CexmcException &  e )
    {
        G4cout << e.what() << G4endl;
    }
    catch ( ... )
    {
        G4cout << "Unknown exception caught" << G4endl;
    }

    delete edStore;
    delete tpStore;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetReconstructor()

CexmcChargeExchangeReconstructor * CexmcEventAction::GetReconstructor ( void  )

inline

Definition at line 177 of file CexmcEventAction.hh.

{
    return reconstructor;
}

Here is the caller graph for this function:

MakeEnergyDepositStore()

CexmcEnergyDepositStore * CexmcEventAction::MakeEnergyDepositStore ( const CexmcEnergyDepositDigitizer *  digitizer )

static

Definition at line 136 of file CexmcEventAction.cc.

{
    G4double  monitorED( digitizer->GetMonitorED() );
    G4double  vetoCounterEDLeft( digitizer->GetVetoCounterEDLeft() );
    G4double  vetoCounterEDRight( digitizer->GetVetoCounterEDRight() );
    G4double  calorimeterEDLeft( digitizer->GetCalorimeterEDLeft() );
    G4double  calorimeterEDRight( digitizer->GetCalorimeterEDRight() );
    G4int     calorimeterEDLeftMaxX( digitizer->GetCalorimeterEDLeftMaxX() );
    G4int     calorimeterEDLeftMaxY( digitizer->GetCalorimeterEDLeftMaxY() );
    G4int     calorimeterEDRightMaxX( digitizer->GetCalorimeterEDRightMaxX() );
    G4int     calorimeterEDRightMaxY( digitizer->GetCalorimeterEDRightMaxY() );

    const CexmcEnergyDepositCalorimeterCollection &
              calorimeterEDLeftCollection(
                            digitizer->GetCalorimeterEDLeftCollection() );
    const CexmcEnergyDepositCalorimeterCollection &
              calorimeterEDRightCollection(
                            digitizer->GetCalorimeterEDRightCollection() );

    /* ATTENTION: return object in heap - must be freed by caller! */
    return new CexmcEnergyDepositStore( monitorED, vetoCounterEDLeft,
                    vetoCounterEDRight, calorimeterEDLeft, calorimeterEDRight,
                    calorimeterEDLeftMaxX, calorimeterEDLeftMaxY,
                    calorimeterEDRightMaxX, calorimeterEDRightMaxY,
                    calorimeterEDLeftCollection, calorimeterEDRightCollection );
}

Here is the call graph for this function:

Here is the caller graph for this function:

MakeTrackPointsStore()

CexmcTrackPointsStore * CexmcEventAction::MakeTrackPointsStore ( const CexmcTrackPointsDigitizer *  digitizer )

static

Definition at line 165 of file CexmcEventAction.cc.

{
    const CexmcTrackPointInfo &
                monitorTP( digitizer->GetMonitorTP() );
    const CexmcTrackPointInfo &
                targetTPBeamParticle(
                    digitizer->GetTargetTPBeamParticle() );
    const CexmcTrackPointInfo &
                targetTPOutputParticle(
                    digitizer->GetTargetTPOutputParticle() );
    const CexmcTrackPointInfo &
                targetTPNucleusParticle(
                    digitizer->GetTargetTPNucleusParticle() );
    const CexmcTrackPointInfo &
                targetTPOutputParticleDecayProductParticle1(
                    digitizer->
                    GetTargetTPOutputParticleDecayProductParticle( 0 ) );
    const CexmcTrackPointInfo &
                targetTPOutputParticleDecayProductParticle2(
                    digitizer->
                    GetTargetTPOutputParticleDecayProductParticle( 1 ) );
    const CexmcTrackPointInfo &
                vetoCounterTPLeft(
                    digitizer->GetVetoCounterTPLeft() );
    const CexmcTrackPointInfo &
                vetoCounterTPRight(
                    digitizer->GetVetoCounterTPRight() );
    const CexmcTrackPointInfo &
                calorimeterTPLeft(
                    digitizer->GetCalorimeterTPLeft() );
    const CexmcTrackPointInfo &
                calorimeterTPRight(
                    digitizer->GetCalorimeterTPRight() );

    /* ATTENTION: return object in heap - must be freed by caller! */
    return new CexmcTrackPointsStore( monitorTP, targetTPBeamParticle,
                  targetTPOutputParticle, targetTPNucleusParticle,
                  targetTPOutputParticleDecayProductParticle1,
                  targetTPOutputParticleDecayProductParticle2,
                  vetoCounterTPLeft, vetoCounterTPRight,
                  calorimeterTPLeft, calorimeterTPRight );
}

Here is the call graph for this function:

Here is the caller graph for this function:

PrintEnergyDeposit()

void CexmcEventAction::PrintEnergyDeposit ( const CexmcEnergyDepositStore *  edStore )

static

Definition at line 210 of file CexmcEventAction.cc.

{
    G4cout << " --- Energy Deposit" << G4endl;
    G4cout << "       monitor : " <<
            G4BestUnit( edStore->monitorED, "Energy" ) << G4endl;
    G4cout << "        vc (l) : " <<
            G4BestUnit( edStore->vetoCounterEDLeft, "Energy" ) << G4endl;
    G4cout << "        vc (r) : " <<
            G4BestUnit( edStore->vetoCounterEDRight, "Energy" ) << G4endl;
    G4cout << "       cal (l) : " <<
            G4BestUnit( edStore->calorimeterEDLeft, "Energy" );
    G4cout << edStore->calorimeterEDLeftCollection;
    G4cout << "       cal (r) : " <<
            G4BestUnit( edStore->calorimeterEDRight, "Energy" );
    G4cout << edStore->calorimeterEDRightCollection;
}

Here is the caller graph for this function:

PrintProductionModelData()

void CexmcEventAction::PrintProductionModelData ( const CexmcAngularRangeList &  angularRanges, const CexmcProductionModelData &  pmData  )

static

Definition at line 259 of file CexmcEventAction.cc.

{
    G4cout << " --- Triggered angular ranges: " << angularRanges;
    G4cout << " --- Production model data: " << pmData;
}

Here is the caller graph for this function:

PrintReconstructedData()

void CexmcEventAction::PrintReconstructedData ( const CexmcAngularRangeList &  angularRanges, const CexmcAngularRange &  angularGap  ) const

private

Definition at line 268 of file CexmcEventAction.cc.

{
    G4cout << " --- Reconstructed data: " << G4endl;
    G4cout << "       -- entry points:" << G4endl;
    G4cout << "              left: " << G4BestUnit(
        reconstructor->GetCalorimeterEPLeftPosition(), "Length" ) << G4endl;
    G4cout << "             right: " << G4BestUnit(
        reconstructor->GetCalorimeterEPRightPosition(), "Length" ) << G4endl;
    G4cout << "            target: " << G4BestUnit(
        reconstructor->GetTargetEPPosition(), "Length" ) << G4endl;
    G4cout << "       -- the angle: " << reconstructor->GetTheAngle() / deg <<
        " deg" << G4endl;
    G4cout << "       -- mass of the output particle: " << G4BestUnit(
        reconstructor->GetOutputParticleMass(), "Energy" ) << G4endl;
    G4cout << "       -- mass of the nucleus output particle: " << G4BestUnit(
        reconstructor->GetNucleusOutputParticleMass(), "Energy" ) << G4endl;
    if ( reconstructor->IsMassCutUsed() )
    {
        if ( reconstructor->HasMassCutTriggered() )
            G4cout << "            < mass cut passed >" << G4endl;
        else
            G4cout << "            < mass cut failed >" << G4endl;
    }
    if ( reconstructor->IsAbsorbedEnergyCutUsed() )
    {
        if ( reconstructor->HasAbsorbedEnergyCutTriggered() )
            G4cout << "            < absorbed energy cut passed >" << G4endl;
        else
            G4cout << "            < absorbed energy cut failed >" << G4endl;
    }
    const CexmcProductionModelData &  pmData(
                                    reconstructor->GetProductionModelData() );
    G4cout << "       -- production model data: " << pmData;
    G4cout << "       -- triggered angular ranges: ";
    if ( triggeredRecAngularRanges.empty() )
        G4cout << "< orphan detected, gap " << angularGap << " >" << G4endl;
    else
        G4cout << triggeredRecAngularRanges;
}

Here is the call graph for this function:

Here is the caller graph for this function:

PrintTrackPoints()

void CexmcEventAction::PrintTrackPoints ( const CexmcTrackPointsStore *  tpStore )

static

Definition at line 229 of file CexmcEventAction.cc.

{
    if ( ! tpStore )
        return;

    G4cout << " --- Track Points" << G4endl;
    G4cout << "       monitor : " << tpStore->monitorTP << G4endl;
    G4cout << "        target : " << tpStore->targetTPBeamParticle << G4endl;
    G4cout << "               : " << tpStore->targetTPOutputParticle << G4endl;
    G4cout << "               : " << tpStore->targetTPNucleusParticle << G4endl;
    G4cout << "               : " <<
            tpStore->targetTPOutputParticleDecayProductParticle1 << G4endl;
    G4cout << "               : " <<
            tpStore->targetTPOutputParticleDecayProductParticle2 << G4endl;
    G4cout << "        vc (l) : " << tpStore->vetoCounterTPLeft << G4endl;
    G4cout << "        vc (r) : " << tpStore->vetoCounterTPRight << G4endl;
    G4cout << "       cal (l) : " << tpStore->calorimeterTPLeft << G4endl;
    G4cout << "       cal (r) : " << tpStore->calorimeterTPRight << G4endl;
    G4cout << "      ---" << G4endl;
    G4cout << "      angle between the " <<
        tpStore->targetTPOutputParticle.particle->GetParticleName() <<
        " decay products : " <<
        tpStore->targetTPOutputParticleDecayProductParticle1.directionWorld.
            angle(
        tpStore->targetTPOutputParticleDecayProductParticle2.directionWorld ) /
            deg << " deg" << G4endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetVerboseDrawLevel()

void CexmcEventAction::SetVerboseDrawLevel ( G4int  value )

inline

Definition at line 170 of file CexmcEventAction.hh.

{
    verboseDraw = value;
}

Here is the caller graph for this function:

SetVerboseOnCexmcLevel()

void CexmcEventAction::SetVerboseOnCexmcLevel ( G4int  value )

inline

Definition at line 164 of file CexmcEventAction.hh.

{
    verbose = value;
}

Here is the caller graph for this function:

UpdateRunHits()

void CexmcEventAction::UpdateRunHits ( const CexmcAngularRangeList &  aRangesReal, const CexmcAngularRangeList &  aRangesRec, G4bool  tpDigitizerHasTriggered, G4bool  edDigitizerHasTriggered, G4bool  edDigitizerMonitorHasTriggered, G4bool  reconstructorHasTriggered, const CexmcAngularRange &  aGap  )

private

Definition at line 668 of file CexmcEventAction.cc.

{
    G4RunManager *    runManager( G4RunManager::GetRunManager() );
    const CexmcRun *  run( static_cast< const CexmcRun * >(
                                                runManager->GetCurrentRun() ) );
    CexmcRun *        theRun( const_cast< CexmcRun * >( run ) );

    if ( tpDigitizerHasTriggered )
    {
        for ( CexmcAngularRangeList::const_iterator  k( aRangesReal.begin() );
                                                k != aRangesReal.end(); ++k )
        {
            theRun->IncrementNmbOfHitsSampledFull( k->index );
            if ( edDigitizerMonitorHasTriggered )
                theRun->IncrementNmbOfHitsSampled( k->index );
            if ( reconstructorHasFullTrigger )
                theRun->IncrementNmbOfHitsTriggeredRealRange( k->index );
        }
        if ( reconstructorHasFullTrigger )
        {
            if ( aRangesRec.empty() )
            {
                theRun->IncrementNmbOfOrphanHits( aGap.index );
            }
            else
            {
                for ( CexmcAngularRangeList::const_iterator
                        k( aRangesRec.begin() ); k != aRangesRec.end(); ++k )
                {
                    theRun->IncrementNmbOfHitsTriggeredRecRange( k->index );
                }
            }
        }
    }
    else
    {
        if ( edDigitizerHasTriggered )
            theRun->IncrementNmbOfFalseHitsTriggeredEDT();
        if ( reconstructorHasFullTrigger )
            theRun->IncrementNmbOfFalseHitsTriggeredRec();
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

messenger

CexmcEventActionMessenger* CexmcEventAction::messenger

private

Definition at line 160 of file CexmcEventAction.hh.

opKinEnergy

G4double CexmcEventAction::opKinEnergy

private

Definition at line 153 of file CexmcEventAction.hh.

physicsManager

CexmcPhysicsManager* CexmcEventAction::physicsManager

private

Definition at line 149 of file CexmcEventAction.hh.

reconstructor

CexmcChargeExchangeReconstructor* CexmcEventAction::reconstructor

private

Definition at line 151 of file CexmcEventAction.hh.

verbose

G4int CexmcEventAction::verbose

private

Definition at line 156 of file CexmcEventAction.hh.

verboseDraw

G4int CexmcEventAction::verboseDraw

private

Definition at line 158 of file CexmcEventAction.hh.

---

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

• CexmcEventAction.hh
• CexmcEventAction.cc