OpParamSD.h

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
/// \file OpParamSD.h
//
/// \author  bjpjones@mit.edu
////////////////////////////////////////////////////////////////////////
//
// 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