In This Package:

DsFastMuonStackAction.h

Go to the documentation of this file.

#ifndef DSOPSTACKACTION_H
#define DSOPSTACKACTION_H 1

#include <iostream>
#include <map>
#include <vector>

#include "globals.hh"
#include "G4UserStackingAction.hh"
#include "GaudiAlg/GaudiTool.h"
#include "GiGa/GiGaStackActionBase.h"
#include "GaudiKernel/RndmGenerators.h"

#ifdef DSFMSA_DEBUGPOSITION
class TTree; 
class TFile; 
#endif

class G4Track;
class G4VTouchable;
class IGeometryInfo;
class ICoordSysSvc;
class IDetectorElement;
using namespace std;

class DsFastMuonStackAction :  public GiGaStackActionBase
{
  public:

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

  virtual G4ClassificationOfNewTrack ClassifyNewTrack( const G4Track* aTrack);
  virtual void NewStage();
  virtual void PrepareNewEvent();
    
  private:
    //    IGeometryInfo* m_geo;
    ICoordSysSvc* m_csvc;
    class DECounter {
        public:
            inline DECounter(G4int maxcount=0, G4double aWeight=1);
            inline void operator++(G4int);
            inline void reset(G4bool all=false); 
            G4bool full()  { return filled; }
            G4int maxCount;
            G4int count;
            G4double nTracks;
            G4bool  filled;
            G4double weight;
    };
    std::map<IDetectorElement*, DECounter> fCounter;
    std::map<const G4VTouchable*, IDetectorElement*> fTouchablesCache;
    std::map<IDetectorElement*, IDetectorElement*> fInsideCache;
    G4bool fTrackRemainingPhotons;

    std::vector<std::string> m_detectors;
    std::vector<G4int>       m_limits;
    std::vector<G4double>    m_weights;
    Rndm::Numbers m_uni;

#ifdef DSFMSA_DEBUGPOSITION
    TTree* tree;
    TFile* file;
    double xx,yy,zz;
    int ii;
    std::map<IDetectorElement*, int> tmpMap;
    int event;
#endif
};


DsFastMuonStackAction::DECounter::DECounter(G4int maxcount, G4double aWeight)
    : maxCount(maxcount)
    , count(0)
    , nTracks(0.)
    , filled(false)
    , weight(aWeight) 
{}

void DsFastMuonStackAction::DECounter::reset(G4bool all){
  if (all) {
      filled=false; 
      count=0;
      nTracks=0.;
  }
  else if ( !filled ) count=0; 
}

void DsFastMuonStackAction::DECounter::operator++(G4int){
    count++;
    if ( count>=maxCount ) filled=true;
}

#endif

---