CexmcEventAction Class Reference

#include <CexmcEventAction.hh>

Inheritance 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 ()
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)

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 ( CexmcPhysicsManager *  physicsManager, G4int  verbose = 0  )

explicit

Definition at line 94 of file CexmcEventAction.cc.

References G4DigiManager::AddNewModule(), CexmcEDDigitizerName, CexmcTPDigitizerName, G4DigiManager::GetDMpointer(), and CexmcPhysicsManager::GetProductionModel().

                                                     :
    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 );
}

CexmcEventAction::~CexmcEventAction ( )

virtual

Definition at line 110 of file CexmcEventAction.cc.

{
    delete reconstructor;
    delete messenger;
}

Member Function Documentation

void CexmcEventAction::BeamParticleChangeHook

(

void

)

Definition at line 117 of file CexmcEventAction.cc.

References CexmcChargeExchangeReconstructor::SetupBeamParticle().

Referenced by CexmcRunManager::BeamParticleChangeHook().

{
    reconstructor->SetupBeamParticle();
}

void CexmcEventAction::BeginOfEventAction ( const G4Event *  event )

virtual

Reimplemented from G4UserEventAction.

Definition at line 123 of file CexmcEventAction.cc.

References CexmcTrackingAction::BeginOfEventAction(), G4RunManager::GetRunManager(), G4RunManager::GetUserTrackingAction(), and CexmcPhysicsManager::ResetNumberOfTriggeredStudiedInteractions().

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

    physicsManager->ResetNumberOfTriggeredStudiedInteractions();
}

void CexmcEventAction::EndOfEventAction ( const G4Event *  event )

virtual

Reimplemented from G4UserEventAction.

Definition at line 799 of file CexmcEventAction.cc.

References CexmcEDDigitizerName, CexmcTPDigitizerName, CexmcWeirdException, CLHEP::HepLorentzVector::cosTheta(), CexmcTrackPointsDigitizer::Digitize(), CexmcEnergyDepositDigitizer::Digitize(), CexmcEventInfo::EdTriggerIsOk(), G4cout, G4endl, GetAngularGaps(), G4DigiManager::GetDMpointer(), G4ParticleDefinition::GetPDGMass(), CexmcPhysicsManager::GetProductionModel(), CexmcChargeExchangeReconstructor::GetProductionModelData(), G4Event::GetUserInformation(), CexmcReconstructor::HasBasicTrigger(), CexmcChargeExchangeReconstructor::HasFullTrigger(), CexmcTrackPointsDigitizer::HasTriggered(), CexmcEnergyDepositDigitizer::HasTriggered(), CexmcTrackPointInfo::IsValid(), MakeEnergyDepositStore(), MakeTrackPointsStore(), CexmcTrackPointInfo::momentumAmp, CexmcEnergyDepositDigitizer::MonitorHasTriggered(), CexmcProductionModelData::outputParticleSCM, CexmcTrackPointInfo::particle, PrintEnergyDeposit(), PrintProductionModelData(), PrintTrackPoints(), CexmcChargeExchangeReconstructor::Reconstruct(), G4Event::SetUserInformation(), CexmcTrackPointsStore::targetTPOutputParticle, and CexmcException::what().

{
    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;
}

CexmcChargeExchangeReconstructor * CexmcEventAction::GetReconstructor ( void  )

inline

Definition at line 177 of file CexmcEventAction.hh.

{
    return reconstructor;
}

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

static

Definition at line 136 of file CexmcEventAction.cc.

References CexmcEnergyDepositDigitizer::GetCalorimeterEDLeft(), CexmcEnergyDepositDigitizer::GetCalorimeterEDLeftCollection(), CexmcEnergyDepositDigitizer::GetCalorimeterEDLeftMaxX(), CexmcEnergyDepositDigitizer::GetCalorimeterEDLeftMaxY(), CexmcEnergyDepositDigitizer::GetCalorimeterEDRight(), CexmcEnergyDepositDigitizer::GetCalorimeterEDRightCollection(), CexmcEnergyDepositDigitizer::GetCalorimeterEDRightMaxX(), CexmcEnergyDepositDigitizer::GetCalorimeterEDRightMaxY(), CexmcEnergyDepositDigitizer::GetMonitorED(), CexmcEnergyDepositDigitizer::GetVetoCounterEDLeft(), and CexmcEnergyDepositDigitizer::GetVetoCounterEDRight().

Referenced by EndOfEventAction().

{
    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 );
}

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

static

Definition at line 165 of file CexmcEventAction.cc.

References CexmcTrackPointsDigitizer::GetCalorimeterTPLeft(), CexmcTrackPointsDigitizer::GetCalorimeterTPRight(), CexmcTrackPointsDigitizer::GetMonitorTP(), CexmcTrackPointsDigitizer::GetTargetTPBeamParticle(), CexmcTrackPointsDigitizer::GetTargetTPNucleusParticle(), CexmcTrackPointsDigitizer::GetTargetTPOutputParticle(), CexmcTrackPointsDigitizer::GetVetoCounterTPLeft(), and CexmcTrackPointsDigitizer::GetVetoCounterTPRight().

Referenced by EndOfEventAction().

{
    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 );
}

void CexmcEventAction::PrintEnergyDeposit ( const CexmcEnergyDepositStore *  edStore )

static

Definition at line 210 of file CexmcEventAction.cc.

References CexmcEnergyDepositStore::calorimeterEDLeft, CexmcEnergyDepositStore::calorimeterEDLeftCollection, CexmcEnergyDepositStore::calorimeterEDRight, CexmcEnergyDepositStore::calorimeterEDRightCollection, G4BestUnit, G4cout, G4endl, CexmcEnergyDepositStore::monitorED, CexmcEnergyDepositStore::vetoCounterEDLeft, and CexmcEnergyDepositStore::vetoCounterEDRight.

Referenced by EndOfEventAction().

{
    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;
}

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

static

Definition at line 259 of file CexmcEventAction.cc.

References G4cout.

Referenced by EndOfEventAction().

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

void CexmcEventAction::PrintTrackPoints ( const CexmcTrackPointsStore *  tpStore )

static

Definition at line 229 of file CexmcEventAction.cc.

References CexmcTrackPointsStore::calorimeterTPLeft, CexmcTrackPointsStore::calorimeterTPRight, python.hepunit::deg, CexmcTrackPointInfo::directionWorld, G4cout, G4endl, G4ParticleDefinition::GetParticleName(), CexmcTrackPointsStore::monitorTP, CexmcTrackPointInfo::particle, CexmcTrackPointsStore::targetTPBeamParticle, CexmcTrackPointsStore::targetTPNucleusParticle, CexmcTrackPointsStore::targetTPOutputParticle, CexmcTrackPointsStore::targetTPOutputParticleDecayProductParticle1, CexmcTrackPointsStore::targetTPOutputParticleDecayProductParticle2, CexmcTrackPointsStore::vetoCounterTPLeft, and CexmcTrackPointsStore::vetoCounterTPRight.

Referenced by EndOfEventAction().

{
    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;
}

void CexmcEventAction::SetVerboseDrawLevel ( G4int  value )

inline

Definition at line 170 of file CexmcEventAction.hh.

Referenced by CexmcEventActionMessenger::SetNewValue().

{
    verboseDraw = value;
}

void CexmcEventAction::SetVerboseOnCexmcLevel ( G4int  value )

inline

Definition at line 164 of file CexmcEventAction.hh.

Referenced by CexmcEventActionMessenger::SetNewValue().

{
    verbose = value;
}

---

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

• CexmcEventAction.hh
• CexmcEventAction.cc