OpParamSD.h

Go to the documentation of this file.

//
//
// This class represents a partially opaque, parameterized optical volume.
//
// The photon transmission probability is determined as a function of
// photon position and momentum direction by a derived class of
// OpParamAction. The particlar derived class to use is specified in a string given to
// the constructor, along with any parameters required to define that object.
// This implementation allows for generic extensions to new types of wireplanes
// or opaque surfaces for future liquid argon TPC detectors.
//
// On each step of a photon within this volume, the GetAttenuationFraction
// method of the OpParamAction derivative is called to get a photon
// transmission probability.  A fraction of all photons are killed
// accordingly.
//
// This sensitive detector object is attached to physical volumes by the
// OpDetReadoutGeometry class.
//
//
// Ben Jones, MIT,2013
//

#ifndef OpParamSD_h
#define OpParamSD_h 1

#include "Geant4/G4String.hh"
#include "Geant4/G4Types.hh"
#include "Geant4/G4VSensitiveDetector.hh"
#include "Geant4/Randomize.hh"

#include <map>
#include <string>
#include <vector>

class G4HCofThisEvent;
class G4Step;
class G4TouchableHistory;

namespace larg4 {

  class OpParamAction;

  class OpParamSD : public G4VSensitiveDetector {

  public:
    OpParamSD(G4String name,
              std::string ModelName,
              int Orientation,
              std::vector<std::vector<double>> Parameters);
    virtual ~OpParamSD() {}

    // Beginning and end of event
    virtual void Initialize(G4HCofThisEvent*);
    virtual void EndOfEvent(G4HCofThisEvent*) {}

    // Tidy up event in abort
    virtual void clear() {}

    // Run per step in sensitive volume
    virtual G4bool ProcessHits(G4Step*, G4TouchableHistory*);

    // Required but empty
    virtual void DrawAll() {}
    virtual void PrintAll() {}

  private:
    G4bool G4BooleanRand(const G4double prob) const;

    OpParamAction* fOpa;
    std::map<G4int, bool> fPhotonAlreadyCrossed;
  };

  inline G4bool OpParamSD::G4BooleanRand(const G4double prob) const
  {
    /* Returns a random boolean variable with the specified probability */
    return (G4UniformRand() < prob);
  }

}

#endif