In This Package:

HistorianStepAction.h

Go to the documentation of this file.

#ifndef HISTORIANSTEPACTION_H
#define HISTORIANSTEPACTION_H

#include "QueriableStepAction.h"
#include "GiGa/GiGaStepActionBase.h"
#include "RuleParser/Queriable.h"
#include "Event/SimTrackReference.h"
#include "Event/SimVertexReference.h"

namespace DayaBay {
  class SimVertex;
  class SimTrack;
  class SimParticleHistory;
}

namespace RuleParser{
  class Rule;
}

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


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

  virtual void UserSteppingAction(const G4Step*);
  
  // 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:
  DayaBay::SimTrackReference    getAncestorTrack()        const;
  DayaBay::SimVertexReference   getAncestorVertex()       const;
  
  // Users are free to override this class, and provide their own way to do things.
  virtual bool    IsInterestingTrack(const G4Track*);
  virtual bool    IsInterestingVertex(const G4Step*);
  
private:
   G4bool   idflag;

  // Configuration.
  std::string              m_TrackSelection;
  std::string              m_VertexSelection;
  G4bool                   m_UseFastMuEnergyCut;  // Defined a flag for using or not use the energy cut for muon simulation.  // Tools and other persistble information
  IHistoryKeeper*               mHistoryKeeper;
  RuleParser::Rule*             mTrackRule;
  RuleParser::Rule*             mVertexRule;
  
  // This is temporary storage for a single processing call.
  DayaBay::SimParticleHistory*     mCurrentHistory;
  DayaBay::SimTrackReference       mCurrentTrackRef;

  StatusCode CreateTrack(const G4Track*, DayaBay::SimTrack* &track);
  StatusCode CreateVertex(DayaBay::SimVertex* &v, 
                          const DayaBay::SimTrackReference& parent); // Only works if vertex hasn't yet been filled.
  
  
  HistorianStepAction (); 
  HistorianStepAction           ( const HistorianStepAction& ); 
  HistorianStepAction& operator=( const HistorianStepAction& ); 
  
  
  enum EHistorianStepParams {
    kPar_HISTORIAN_START = kPar_END_QUERIABLE
    // Tracky things specific to Histories:
    , kPar_ParentPdg                  // int
    , kPar_ParentIndirection          // int
    , kPar_GrandParentPdg             // int
    , kPar_GrandParentIndirection     // int
  
     // Vertex-y things only applicable if we have current history information.
    , kPar_distanceFromLastVertex     // double
    , kPar_timeSinceLastVertex        // double
    , kPar_energyLossSinceLastVertex  // double
    , kPar_angleFromLastVertex        // double
    
  };
};



#endif /* MYSTEPACTION_H */

---