CexmcReconstructor Class Reference

#include <CexmcReconstructor.hh>

Inheritance diagram for CexmcReconstructor:

Collaboration diagram for CexmcReconstructor:

Public Member Functions
	CexmcReconstructor ()
virtual	~CexmcReconstructor ()
virtual void	Reconstruct (const CexmcEnergyDepositStore *edStore)
void	SetCalorimeterEntryPointDefinitionAlgorithm (CexmcCalorimeterEntryPointDefinitionAlgorithm algo)
void	SetCalorimeterEntryPointDepthDefinitionAlgorithm (CexmcCalorimeterEntryPointDepthDefinitionAlgorithm algo)
void	SetCrystalSelectionAlgorithm (CexmcCrystalSelectionAlgorithm algo)
void	UseInnerRefCrystal (G4bool on=true)
void	SetCalorimeterEntryPointDepth (G4double depth)
CexmcCalorimeterEntryPointDefinitionAlgorithm	GetCalorimeterEntryPointDefinitionAlgorithm (void) const
CexmcCalorimeterEntryPointDepthDefinitionAlgorithm	GetCalorimeterEntryPointDepthDefinitionAlgorithm (void) const
CexmcCrystalSelectionAlgorithm	GetCrystalSelectionAlgorithm (void) const
G4bool	IsInnerRefCrystalUsed (void) const
G4double	GetCalorimeterEntryPointDepth (void) const
const G4ThreeVector &	GetCalorimeterEPLeftPosition (void) const
const G4ThreeVector &	GetCalorimeterEPRightPosition (void) const
const G4ThreeVector &	GetCalorimeterEPLeftDirection (void) const
const G4ThreeVector &	GetCalorimeterEPRightDirection (void) const
const G4ThreeVector &	GetTargetEPPosition (void) const
const G4ThreeVector &	GetTargetEPDirection (void) const
const G4ThreeVector &	GetCalorimeterEPLeftWorldPosition (void) const
const G4ThreeVector &	GetCalorimeterEPRightWorldPosition (void) const
const G4ThreeVector &	GetCalorimeterEPLeftWorldDirection (void) const
const G4ThreeVector &	GetCalorimeterEPRightWorldDirection (void) const
const G4ThreeVector &	GetTargetEPWorldPosition (void) const
const G4ThreeVector &	GetTargetEPWorldDirection (void) const
G4double	GetTheAngle (void) const
G4bool	HasBasicTrigger (void) const
virtual G4bool	HasFullTrigger (void) const

Protected Member Functions
void	ReconstructEntryPoints (const CexmcEnergyDepositStore *edStore)
void	ReconstructTargetPoint (void)
void	ReconstructAngle (void)

Protected Attributes
G4bool	hasBasicTrigger
CexmcCalorimeterEntryPointDefinitionAlgorithm	epDefinitionAlgorithm
CexmcCalorimeterEntryPointDepthDefinitionAlgorithm	epDepthDefinitionAlgorithm
CexmcCrystalSelectionAlgorithm	csAlgorithm
G4bool	useInnerRefCrystal
G4double	epDepth
G4ThreeVector	calorimeterEPLeftPosition
G4ThreeVector	calorimeterEPRightPosition
G4ThreeVector	calorimeterEPLeftDirection
G4ThreeVector	calorimeterEPRightDirection
G4ThreeVector	targetEPPosition
G4ThreeVector	targetEPDirection
G4ThreeVector	calorimeterEPLeftWorldPosition
G4ThreeVector	calorimeterEPRightWorldPosition
G4ThreeVector	calorimeterEPLeftWorldDirection
G4ThreeVector	calorimeterEPRightWorldDirection
G4ThreeVector	targetEPWorldPosition
G4ThreeVector	targetEPWorldDirection
G4double	theAngle
G4double	calorimeterEDLeftAdjacent
G4double	calorimeterEDRightAdjacent
G4bool	collectEDInAdjacentCrystals

Detailed Description

Definition at line 56 of file CexmcReconstructor.hh.

Constructor & Destructor Documentation

CexmcReconstructor::CexmcReconstructor ( )

explicit

Definition at line 50 of file CexmcReconstructor.cc.

                                       : hasBasicTrigger( false ),
    epDefinitionAlgorithm( CexmcEntryPointBySqrtEDWeights ),
    epDepthDefinitionAlgorithm( CexmcEntryPointDepthPlain ),
    csAlgorithm( CexmcSelectAllCrystals ), useInnerRefCrystal( false ),
    epDepth( 0 ), theAngle( 0 ), calorimeterEDLeftAdjacent( 0 ),
    calorimeterEDRightAdjacent( 0 ), collectEDInAdjacentCrystals( false ),
    targetEPInitialized( false ), messenger( NULL )
{
    G4RunManager *      runManager( G4RunManager::GetRunManager() );
    const CexmcSetup *  setup( static_cast< const CexmcSetup * >(
                                runManager->GetUserDetectorConstruction() ) );
    calorimeterGeometry = setup->GetCalorimeterGeometry();
    targetTransform = setup->GetTargetTransform();
    calorimeterLeftTransform = setup->GetCalorimeterLeftTransform();
    calorimeterRightTransform = setup->GetCalorimeterRightTransform();

    messenger = new CexmcReconstructorMessenger( this );
}

Here is the call graph for this function:

CexmcReconstructor::~CexmcReconstructor ( )

virtual

Definition at line 70 of file CexmcReconstructor.cc.

{
    delete messenger;
}

Member Function Documentation

CexmcCalorimeterEntryPointDefinitionAlgorithm CexmcReconstructor::GetCalorimeterEntryPointDefinitionAlgorithm ( void  ) const

inline

Definition at line 248 of file CexmcReconstructor.hh.

{
    return epDefinitionAlgorithm;
}

G4double CexmcReconstructor::GetCalorimeterEntryPointDepth ( void  ) const

inline

Definition at line 276 of file CexmcReconstructor.hh.

{
    return epDepth;
}

CexmcCalorimeterEntryPointDepthDefinitionAlgorithm CexmcReconstructor::GetCalorimeterEntryPointDepthDefinitionAlgorithm ( void  ) const

inline

Definition at line 256 of file CexmcReconstructor.hh.

{
    return epDepthDefinitionAlgorithm;
}

const G4ThreeVector & CexmcReconstructor::GetCalorimeterEPLeftDirection ( void  ) const

inline

Definition at line 297 of file CexmcReconstructor.hh.

{
    return calorimeterEPLeftDirection;
}

const G4ThreeVector & CexmcReconstructor::GetCalorimeterEPLeftPosition ( void  ) const

inline

Definition at line 283 of file CexmcReconstructor.hh.

{
    return calorimeterEPLeftPosition;
}

const G4ThreeVector & CexmcReconstructor::GetCalorimeterEPLeftWorldDirection ( void  ) const

inline

Definition at line 339 of file CexmcReconstructor.hh.

{
    return calorimeterEPLeftWorldDirection;
}

const G4ThreeVector & CexmcReconstructor::GetCalorimeterEPLeftWorldPosition ( void  ) const

inline

Definition at line 325 of file CexmcReconstructor.hh.

{
    return calorimeterEPLeftWorldPosition;
}

const G4ThreeVector & CexmcReconstructor::GetCalorimeterEPRightDirection ( void  ) const

inline

Definition at line 304 of file CexmcReconstructor.hh.

{
    return calorimeterEPRightDirection;
}

const G4ThreeVector & CexmcReconstructor::GetCalorimeterEPRightPosition ( void  ) const

inline

Definition at line 290 of file CexmcReconstructor.hh.

{
    return calorimeterEPRightPosition;
}

const G4ThreeVector & CexmcReconstructor::GetCalorimeterEPRightWorldDirection ( void  ) const

inline

Definition at line 346 of file CexmcReconstructor.hh.

{
    return calorimeterEPRightWorldDirection;
}

const G4ThreeVector & CexmcReconstructor::GetCalorimeterEPRightWorldPosition ( void  ) const

inline

Definition at line 332 of file CexmcReconstructor.hh.

{
    return calorimeterEPRightWorldPosition;
}

CexmcCrystalSelectionAlgorithm CexmcReconstructor::GetCrystalSelectionAlgorithm ( void  ) const

inline

Definition at line 264 of file CexmcReconstructor.hh.

{
    return csAlgorithm;
}

const G4ThreeVector & CexmcReconstructor::GetTargetEPDirection ( void  ) const

inline

Definition at line 318 of file CexmcReconstructor.hh.

{
    return targetEPDirection;
}

const G4ThreeVector & CexmcReconstructor::GetTargetEPPosition ( void  ) const

inline

Definition at line 311 of file CexmcReconstructor.hh.

{
    return targetEPPosition;
}

const G4ThreeVector & CexmcReconstructor::GetTargetEPWorldDirection ( void  ) const

inline

Definition at line 360 of file CexmcReconstructor.hh.

{
    return targetEPWorldDirection;
}

const G4ThreeVector & CexmcReconstructor::GetTargetEPWorldPosition ( void  ) const

inline

Definition at line 353 of file CexmcReconstructor.hh.

{
    return targetEPWorldPosition;
}

G4double CexmcReconstructor::GetTheAngle ( void  ) const

inline

Definition at line 366 of file CexmcReconstructor.hh.

{
    return theAngle;
}

G4bool CexmcReconstructor::HasBasicTrigger ( void  ) const

inline

Definition at line 372 of file CexmcReconstructor.hh.

{
    return hasBasicTrigger;
}

Here is the caller graph for this function:

G4bool CexmcReconstructor::HasFullTrigger ( void  ) const

virtual

Reimplemented in CexmcChargeExchangeReconstructor.

Definition at line 87 of file CexmcReconstructor.cc.

{
    return hasBasicTrigger;
}

G4bool CexmcReconstructor::IsInnerRefCrystalUsed ( void  ) const

inline

Definition at line 270 of file CexmcReconstructor.hh.

{
    return useInnerRefCrystal;
}

void CexmcReconstructor::Reconstruct ( const CexmcEnergyDepositStore *  edStore )

virtual

Reimplemented in CexmcChargeExchangeReconstructor.

Definition at line 76 of file CexmcReconstructor.cc.

{
    ReconstructEntryPoints( edStore );
    if ( hasBasicTrigger )
        ReconstructTargetPoint();
    if ( hasBasicTrigger )
        ReconstructAngle();
}

Here is the call graph for this function:

void CexmcReconstructor::ReconstructAngle ( void  )

protected

Definition at line 268 of file CexmcReconstructor.cc.

{
    G4ThreeVector  epLeft( calorimeterEPLeftWorldPosition -
                           targetEPWorldPosition );
    G4ThreeVector  epRight( calorimeterEPRightWorldPosition -
                            targetEPWorldPosition );
    theAngle = epLeft.angle( epRight );

    hasBasicTrigger = true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void CexmcReconstructor::ReconstructEntryPoints ( const CexmcEnergyDepositStore *  edStore )

protected

Definition at line 93 of file CexmcReconstructor.cc.

{
    G4int     columnLeft( edStore->calorimeterEDLeftMaxX );
    G4int     rowLeft( edStore->calorimeterEDLeftMaxY );
    G4int     columnRight( edStore->calorimeterEDRightMaxX );
    G4int     rowRight( edStore->calorimeterEDRightMaxY );
    G4double  crystalLength( calorimeterGeometry.crystalLength );

    if ( useInnerRefCrystal )
    {
        TransformToAdjacentInnerCrystal( columnLeft, rowLeft );
        TransformToAdjacentInnerCrystal( columnRight, rowRight );
    }

    calorimeterEPLeftPosition.setX( 0 );
    calorimeterEPLeftPosition.setY( 0 );
    calorimeterEPLeftPosition.setZ( -crystalLength / 2 + epDepth );
    calorimeterEPLeftDirection.setX( 0 );
    calorimeterEPLeftDirection.setY( 0 );
    calorimeterEPLeftDirection.setZ( 0 );
    calorimeterEPRightPosition.setX( 0 );
    calorimeterEPRightPosition.setY( 0 );
    calorimeterEPRightPosition.setZ( -crystalLength / 2 + epDepth );
    calorimeterEPRightDirection.setX( 0 );
    calorimeterEPRightDirection.setY( 0 );
    calorimeterEPRightDirection.setZ( 0 );

    G4bool  edInAdjacentCrystalsCollected( false );

    switch ( epDefinitionAlgorithm )
    {
    case CexmcEntryPointInTheCenter :
        break;
    case CexmcEntryPointInTheCenterOfCrystalWithMaxED :
        {
            G4int     nCrystalsInColumn(
                                    calorimeterGeometry.nCrystalsInColumn );
            G4int     nCrystalsInRow( calorimeterGeometry.nCrystalsInRow );
            G4double  crystalWidth( calorimeterGeometry.crystalWidth );
            G4double  crystalHeight( calorimeterGeometry.crystalHeight );

            calorimeterEPLeftPosition.setX(
                ( G4double( columnLeft ) - G4double( nCrystalsInRow ) / 2 ) *
                                        crystalWidth  + crystalWidth / 2 );
            calorimeterEPLeftPosition.setY(
                ( G4double( rowLeft ) - G4double( nCrystalsInColumn ) / 2 ) *
                                        crystalHeight + crystalHeight / 2 );
            calorimeterEPRightPosition.setX(
                ( G4double( columnRight ) - G4double( nCrystalsInRow ) / 2 ) *
                                        crystalWidth  + crystalWidth / 2 );
            calorimeterEPRightPosition.setY(
                ( G4double( rowRight ) - G4double( nCrystalsInColumn ) / 2 ) *
                                        crystalHeight + crystalHeight / 2 );
        }
        break;
    case CexmcEntryPointByLinearEDWeights :
    case CexmcEntryPointBySqrtEDWeights :
        {
            G4double  x( 0 );
            G4double  y( 0 );

            CalculateWeightedEPPosition( edStore->calorimeterEDLeftCollection,
                                         rowLeft, columnLeft, x, y,
                                         calorimeterEDLeftAdjacent );
            calorimeterEPLeftPosition.setX( x );
            calorimeterEPLeftPosition.setY( y );
            CalculateWeightedEPPosition( edStore->calorimeterEDRightCollection,
                                         rowRight, columnRight, x, y,
                                         calorimeterEDRightAdjacent );
            calorimeterEPRightPosition.setX( x );
            calorimeterEPRightPosition.setY( y );

            if ( csAlgorithm == CexmcSelectAdjacentCrystals )
                edInAdjacentCrystalsCollected = true;
        }
        break;
    default :
        break;
    }

    switch ( epDepthDefinitionAlgorithm )
    {
    case CexmcEntryPointDepthPlain :
        break;
    case CexmcEntryPointDepthSphere :
        {
            G4double  calorimeterEPLeftRadiusOfTheSphere(
                             calorimeterLeftTransform.NetTranslation().mag() +
                             calorimeterEPLeftPosition.z() );
            G4double  calorimeterEPLeftRadiusOfTheSphere2(
                                      calorimeterEPLeftRadiusOfTheSphere *
                                      calorimeterEPLeftRadiusOfTheSphere );
            G4double  calorimeterEPLeftPositionX2(
                                            calorimeterEPLeftPosition.x() *
                                            calorimeterEPLeftPosition.x() );
            G4double  calorimeterEPLeftPositionY2(
                                            calorimeterEPLeftPosition.y() *
                                            calorimeterEPLeftPosition.y() );
            G4double  calorimeterEPLeftPositionZOffset(
                           calorimeterEPLeftRadiusOfTheSphere - std::sqrt(
                                  calorimeterEPLeftRadiusOfTheSphere2 -
                                  calorimeterEPLeftPositionX2 -
                                  calorimeterEPLeftPositionY2 ) );
            G4double  calorimeterEPRightRadiusOfTheSphere(
                              calorimeterRightTransform.NetTranslation().mag() +
                              calorimeterEPRightPosition.z() );
            G4double  calorimeterEPRightRadiusOfTheSphere2(
                                       calorimeterEPRightRadiusOfTheSphere *
                                       calorimeterEPRightRadiusOfTheSphere );
            G4double  calorimeterEPRightPositionX2(
                                            calorimeterEPRightPosition.x() *
                                            calorimeterEPRightPosition.x() );
            G4double  calorimeterEPRightPositionY2(
                                            calorimeterEPRightPosition.y() *
                                            calorimeterEPRightPosition.y() );
            G4double  calorimeterEPRightPositionZOffset(
                            calorimeterEPRightRadiusOfTheSphere - std::sqrt(
                                    calorimeterEPRightRadiusOfTheSphere2 -
                                    calorimeterEPRightPositionX2 -
                                    calorimeterEPRightPositionY2 ) );
            calorimeterEPLeftPosition.setZ( calorimeterEPLeftPosition.z() -
                                            calorimeterEPLeftPositionZOffset );
            calorimeterEPRightPosition.setZ( calorimeterEPRightPosition.z() -
                                         calorimeterEPRightPositionZOffset );
        }
        break;
    default :
        break;
    }

    calorimeterEPLeftWorldPosition = calorimeterLeftTransform.TransformPoint(
                                                 calorimeterEPLeftPosition );
    calorimeterEPLeftWorldDirection = calorimeterLeftTransform.TransformAxis(
                                                 calorimeterEPLeftDirection );
    calorimeterEPRightWorldPosition = calorimeterRightTransform.TransformPoint(
                                                 calorimeterEPRightPosition );
    calorimeterEPRightWorldDirection = calorimeterRightTransform.TransformAxis(
                                                 calorimeterEPRightDirection );

    if ( collectEDInAdjacentCrystals && ! edInAdjacentCrystalsCollected )
    {
        CollectEDInAdjacentCrystals( edStore->calorimeterEDLeftCollection,
                                     rowLeft, columnLeft,
                                     calorimeterEDLeftAdjacent );
        CollectEDInAdjacentCrystals( edStore->calorimeterEDRightCollection,
                                     rowRight, columnRight,
                                     calorimeterEDRightAdjacent );
    }

    hasBasicTrigger = true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void CexmcReconstructor::ReconstructTargetPoint ( void  )

protected

Definition at line 247 of file CexmcReconstructor.cc.

{
    if ( ! targetEPInitialized )
    {
        targetEPWorldPosition = targetTransform.TransformPoint(
                                                    G4ThreeVector( 0, 0, 0 ) );
        targetEPWorldDirection.setX( 0 );
        targetEPWorldDirection.setY( 0 );
        targetEPWorldDirection.setZ( 1 );

        targetEPPosition = targetTransform.Inverse().TransformPoint(
                                                    targetEPWorldPosition );
        targetEPDirection = targetTransform.Inverse().TransformAxis(
                                                    targetEPWorldDirection );
        targetEPInitialized = true;
    }

    hasBasicTrigger = true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void CexmcReconstructor::SetCalorimeterEntryPointDefinitionAlgorithm ( CexmcCalorimeterEntryPointDefinitionAlgorithm  algo )

inline

Definition at line 212 of file CexmcReconstructor.hh.

{
    epDefinitionAlgorithm = algo;
}

Here is the caller graph for this function:

void CexmcReconstructor::SetCalorimeterEntryPointDepth ( G4double  depth )

inline

Definition at line 240 of file CexmcReconstructor.hh.

{
    epDepth = depth;
}

Here is the caller graph for this function:

void CexmcReconstructor::SetCalorimeterEntryPointDepthDefinitionAlgorithm ( CexmcCalorimeterEntryPointDepthDefinitionAlgorithm  algo )

inline

Definition at line 220 of file CexmcReconstructor.hh.

{
    epDepthDefinitionAlgorithm = algo;
}

Here is the caller graph for this function:

void CexmcReconstructor::SetCrystalSelectionAlgorithm ( CexmcCrystalSelectionAlgorithm  algo )

inline

Definition at line 227 of file CexmcReconstructor.hh.

{
    csAlgorithm = algo;
}

Here is the caller graph for this function:

void CexmcReconstructor::UseInnerRefCrystal ( G4bool  on = true )

inline

Definition at line 234 of file CexmcReconstructor.hh.

{
    useInnerRefCrystal = on;
}

Here is the caller graph for this function:

Member Data Documentation

G4double CexmcReconstructor::calorimeterEDLeftAdjacent

protected

Definition at line 190 of file CexmcReconstructor.hh.

G4double CexmcReconstructor::calorimeterEDRightAdjacent

protected

Definition at line 192 of file CexmcReconstructor.hh.

G4ThreeVector CexmcReconstructor::calorimeterEPLeftDirection

protected

Definition at line 167 of file CexmcReconstructor.hh.

G4ThreeVector CexmcReconstructor::calorimeterEPLeftPosition

protected

Definition at line 163 of file CexmcReconstructor.hh.

G4ThreeVector CexmcReconstructor::calorimeterEPLeftWorldDirection

protected

Definition at line 179 of file CexmcReconstructor.hh.

G4ThreeVector CexmcReconstructor::calorimeterEPLeftWorldPosition

protected

Definition at line 175 of file CexmcReconstructor.hh.

G4ThreeVector CexmcReconstructor::calorimeterEPRightDirection

protected

Definition at line 169 of file CexmcReconstructor.hh.

G4ThreeVector CexmcReconstructor::calorimeterEPRightPosition

protected

Definition at line 165 of file CexmcReconstructor.hh.

G4ThreeVector CexmcReconstructor::calorimeterEPRightWorldDirection

protected

Definition at line 181 of file CexmcReconstructor.hh.

G4ThreeVector CexmcReconstructor::calorimeterEPRightWorldPosition

protected

Definition at line 177 of file CexmcReconstructor.hh.

G4bool CexmcReconstructor::collectEDInAdjacentCrystals

protected

Definition at line 194 of file CexmcReconstructor.hh.

CexmcCrystalSelectionAlgorithm CexmcReconstructor::csAlgorithm

protected

Definition at line 156 of file CexmcReconstructor.hh.

CexmcCalorimeterEntryPointDefinitionAlgorithm CexmcReconstructor::epDefinitionAlgorithm

protected

Definition at line 151 of file CexmcReconstructor.hh.

G4double CexmcReconstructor::epDepth

protected

Definition at line 160 of file CexmcReconstructor.hh.

CexmcCalorimeterEntryPointDepthDefinitionAlgorithm CexmcReconstructor::epDepthDefinitionAlgorithm

protected

Definition at line 154 of file CexmcReconstructor.hh.

G4bool CexmcReconstructor::hasBasicTrigger

protected

Definition at line 148 of file CexmcReconstructor.hh.

G4ThreeVector CexmcReconstructor::targetEPDirection

protected

Definition at line 173 of file CexmcReconstructor.hh.

G4ThreeVector CexmcReconstructor::targetEPPosition

protected

Definition at line 171 of file CexmcReconstructor.hh.

G4ThreeVector CexmcReconstructor::targetEPWorldDirection

protected

Definition at line 185 of file CexmcReconstructor.hh.

G4ThreeVector CexmcReconstructor::targetEPWorldPosition

protected

Definition at line 183 of file CexmcReconstructor.hh.

G4double CexmcReconstructor::theAngle

protected

Definition at line 187 of file CexmcReconstructor.hh.

G4bool CexmcReconstructor::useInnerRefCrystal

protected

Definition at line 158 of file CexmcReconstructor.hh.

---

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

• source/geant4.10.03.p03/examples/advanced/ChargeExchangeMC/include/CexmcReconstructor.hh
• source/geant4.10.03.p03/examples/advanced/ChargeExchangeMC/src/CexmcReconstructor.cc