larg4::OverlaidWireplanesAction Class Reference

#include "OpParamAction.h"

Inheritance diagram for larg4::OverlaidWireplanesAction:

Public Member Functions
	OverlaidWireplanesAction (std::vector< std::vector< double >>, int)
	~OverlaidWireplanesAction ()
double	GetAttenuationFraction (G4ThreeVector PhotonDirection, G4ThreeVector PhotonPosition)

Private Attributes
G4ThreeVector	fPlaneNormal
std::vector< G4ThreeVector >	fWireDirections
std::vector< double >	fDPRatios

Detailed Description

Definition at line 92 of file OpParamAction.h.

Constructor & Destructor Documentation

larg4::OverlaidWireplanesAction::OverlaidWireplanesAction	(	std::vector< std::vector< double >>	InputVectors,
		int	Orientation
	)

Definition at line 66 of file OpParamAction.cxx.

  {

    G4ThreeVector WireBasis1, WireBasis2;

    if (Orientation == 0) {
      fPlaneNormal = G4ThreeVector(1, 0, 0);
      WireBasis1 = G4ThreeVector(0, 1, 0);
      WireBasis2 = G4ThreeVector(0, 0, 1);
    }
    else if (Orientation == 1) {
      fPlaneNormal = G4ThreeVector(0, 1, 0);
      WireBasis1 = G4ThreeVector(1, 0, 0);
      WireBasis2 = G4ThreeVector(0, 0, 1);
    }
    else if (Orientation == 2) {
      fPlaneNormal = G4ThreeVector(0, 0, 1);
      WireBasis1 = G4ThreeVector(0, 1, 0);
      WireBasis2 = G4ThreeVector(0, 0, 1);
    }
    else {
      throw cet::exception("OpParamAction")
        << "Unrecognized wireplane orientation. Options are 1=Xdrift, 2=Ydrift, 3=Zdrift\n";
    }
    for (size_t i = 0; i != InputVectors.size(); ++i) {
      if (InputVectors.at(i).size() != 3) {
        throw cet::exception("OpParamAction")
          << "Unrecognized wireplane parameter format. Expected vector(3)'s with v[0] = wire "
             "angle, v[1] = wire pitch, v[2] = wire diameter\n";
      }
      else {
        double theta = InputVectors[i][0] * 3.142 / 180.;
        fWireDirections.push_back(cos(theta) * WireBasis1 + sin(theta) * WireBasis2);
        fDPRatios.push_back(InputVectors[i][2] / InputVectors[i][1]);
      }
    }
  }

larg4::OverlaidWireplanesAction::~OverlaidWireplanesAction

(

)

Definition at line 107 of file OpParamAction.cxx.

{}

Member Function Documentation

double larg4::OverlaidWireplanesAction::GetAttenuationFraction ( G4ThreeVector  PhotonDirection, G4ThreeVector  PhotonPosition  )

virtual

Reimplemented from larg4::OpParamAction.

Definition at line 111 of file OpParamAction.cxx.

References util::abs().

  {

    double AttenFraction = 1.;

    for (size_t i = 0; i != fWireDirections.size(); ++i) {
      G4ThreeVector ProjDirection =
        PhotonDirection -
        fWireDirections.at(i) * (fWireDirections.at(i).dot(PhotonDirection.unit()));

      // fWireDirections.at(i).cross(PhotonDirection.cross(fWireDirections.at(i)).unit());
      double CosTheta =
        (ProjDirection.mag() > 0 ? std::abs(fPlaneNormal.dot(ProjDirection)) / ProjDirection.mag() :
                                   1.0);
      if (CosTheta < fDPRatios.at(i))
        return 0;
      else
        AttenFraction *= (1.0 - fDPRatios.at(i) / CosTheta);
    }
    return AttenFraction;
  }

Member Data Documentation

std::vector<double> larg4::OverlaidWireplanesAction::fDPRatios

private

Definition at line 102 of file OpParamAction.h.

G4ThreeVector larg4::OverlaidWireplanesAction::fPlaneNormal

private

Definition at line 100 of file OpParamAction.h.

std::vector<G4ThreeVector> larg4::OverlaidWireplanesAction::fWireDirections

private

Definition at line 101 of file OpParamAction.h.

---

The documentation for this class was generated from the following files:

• larsim/v09_43_00/source/larsim/LegacyLArG4/OpParamAction.h
• larsim/v09_43_00/source/larsim/LegacyLArG4/OpParamAction.cxx