In This Package:

QueriableStepAction.h

Go to the documentation of this file.

#ifndef QUERIABLESTEPACTION_H
#define QUERIABLESTEPACTION_H

#include "GiGa/GiGaStepActionBase.h"
#include "G4DataHelpers/INeutronCaptureInfo.h"
#include "RuleParser/Queriable.h"
#include "Event/SimProcess.h"
#include <map>

namespace RuleParser{
  class Rule;
  class ParameterList;
}

class ITouchableToDetectorElement;
class IHistoryKeeper;
class IDetectorElement;
class ILVolume;
class IPVolume;
class G4Step;
class G4Track;
class G4StepPoint;
class G4HistoryUserTrackInfo;

// Inherit from GiGaBase to pick up Services.

class QueriableStepAction: public RuleParser::Queriable,  public GiGaStepActionBase
{
public:
  QueriableStepAction( const std::string& type   , 
                     const std::string& name   , 
                     const IInterface*  parent ) ;
  virtual ~QueriableStepAction() {};
  
  virtual StatusCode         initialize ();
  virtual StatusCode         finalize   ();

  virtual void UserSteppingAction(const G4Step*) =0 ;  // NB This is a mix-in object!

  virtual StatusCode  GetTrackParameterList(RuleParser::ParameterList&);
  virtual StatusCode  GetVertexParameterList(RuleParser::ParameterList&);
  
  void Reset(const G4Step*, const G4Track*, const G4StepPoint*);
  
  // Things for evaluating selection rules:
  virtual void      queryParam(int id,double&      output) const;
  virtual void      queryParam(int id,std::string& output) const;

  // Things for above, or for filling the track:
  const IDetectorElement*       getDetectorElement()      const;
  int                           getDetectorElementMatch() const;
  const ILVolume*               getLogicalVolume()        const;
  const IPVolume*               getPhysicalVolume()       const;
  const DayaBay::SimProcess&    getProcess()              const;
  double                        getQuenchedEnergy()       const;
  const G4HistoryUserTrackInfo* getUserTrackInfo()        const;
  unsigned int                  getDetectorId()            const;
  unsigned int                  getDetectorId(const IDetectorElement*) const;
  
protected:
  // Configuration:
  std::vector<std::string> m_DetectorElementSearchPath;
  std::string              m_TouchableToDetelem_name;
  std::vector<std::string> m_IdParameterNames;
    
  // This is temporary storage for a single processing call.
  const G4Step*                    mCurrentStep;
  const G4Track*                   mCurrentTrack;
  const G4StepPoint*               mCurrentStepPoint; // Pre or post, depends on what we're doing.
  ITouchableToDetectorElement*     mTouchToDetElem;  

  // Things that are difficult to calculate, and are therefore cached in case
  // the user references them twice.
  mutable const IDetectorElement*  mDetElement;
  mutable int                      mDetElementMatch;
  mutable const ILVolume*          mLogicVolume;
  mutable const IPVolume*          mPhysVolume;
  mutable DayaBay::SimProcess      mProcess;
  mutable double                   mQuenchedEnergy;
  mutable G4HistoryUserTrackInfo*  mUserTrackInfo;
  mutable unsigned int             mDetectorId;
  typedef std::map<const IDetectorElement*,unsigned int> DetectorIdCache_t;
  mutable DetectorIdCache_t        mDetectorIdCache;


 protected:
   enum EParamIds {
     kPar_Unknown = 0
     , kPar_START_VERTEX
     // Applies to end-of-step (Vertex) only
     , kPar_Step_dE                    // double
     , kPar_Step_dE_Ion                // double
     , kPar_Step_qdE                   // double
     , kPar_Step_dx                    // double
     , kPar_Step_dt                    // double
     , kPar_Step_dAngle                // double
     , kPar_capTargetZ                 // int
     , kPar_capTargetA                 // int | No correct A info. 8/2008 Wei
     
     , kPar_Step_E_weighted_x          // double
     , kPar_Step_E_weighted_y          // double
     , kPar_Step_E_weighted_z          // double
     , kPar_Step_E_weighted_t          // double
     , kPar_Step_qE_weighted_x         // double
     , kPar_Step_qE_weighted_y         // double
     , kPar_Step_qE_weighted_z         // double
     , kPar_Step_qE_weighted_t         // double

     , kPar_IsStopping                 // int
     , kPar_IsStarting                 // int
     , kPar_StepNumber                 // int

     , kPar_VolumeChanged              // int as bool
     , kPar_MaterialChanged            // int as bool

     // Both Vertex and Track:
     , kPar_START_TRACK
     , kPar_t                          // double   |
     , kPar_x                          // double   |
     , kPar_y                          // double   |
     , kPar_z                          // double   |
     , kPar_r                          // double   |
     , kPar_local_x                    // double   |
     , kPar_local_y                    // double   |
     , kPar_local_z                    // double   | all of these guys refer to prepoint for tracks
     , kPar_local_r                    // double   | or postpoint for vertices.
     , kPar_LogicalVolumeName          // string
     , kPar_MaterialName               // string
     , kPar_DetectorElement            // custom
     , kPar_DetectorElementName        // string           
     , kPar_DetectorElementMatch       // int
     , kPar_NicheId                    // int
     , kPar_DetectorId                 // int
     , kPar_SiteId                     // int
     , kPar_Site                       // int
     , kPar_AD                         // int
     , kPar_momentum                   // double   |
     , kPar_TotEnergy                  // double   |
     , kPar_KineticEnergy              // double   |
     , kPar_vx                         // double   |
     , kPar_vy                         // double   |
     , kPar_vz                         // double   |
     , kPar_local_vx                   // double   |
     , kPar_local_vy                   // double   |
     , kPar_local_vz                   // double   |
     , kPar_ProcessType                // int      |
     , kPar_ProcessName                // string   |

     , kPar_Pdg                        // int 
     , kPar_Charge                    // int
     , kPar_TrackId                    // int
     , kPar_CreatorPdg                 // int
     , kPar_AncestorPdg                // int
     , kPar_mass                       // double
     , kPar_ParticleName               // string
     , kPar_Prescale                   // custom
     , kPar_CreatorProcessName         // string
     , kPar_END_VERTEX
     , kPar_END_QUERIABLE
   }; 
   
  private:
    QueriableStepAction (); 
    QueriableStepAction           ( const QueriableStepAction& ); 
    QueriableStepAction& operator=( const QueriableStepAction& ); 
    INeutronCaptureInfo* m_capinfo; // A Gaudi Tool passing n-cap info around
    
    // Configuration for Birks constant
    double  m_BirksConstant1;
    double  m_BirksConstant2;
};



#endif /* QUERIABLESTEPACTION_HH */

---