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 96 of file OpParamAction.h.

Constructor & Destructor Documentation

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

Definition at line 85 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 130 of file OpParamAction.cxx.

  {
  }

Member Function Documentation

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

virtual

Reimplemented from larg4::OpParamAction.

Definition at line 137 of file OpParamAction.cxx.

  {

    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 108 of file OpParamAction.h.

G4ThreeVector larg4::OverlaidWireplanesAction::fPlaneNormal

private

Definition at line 106 of file OpParamAction.h.

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

private

Definition at line 107 of file OpParamAction.h.

---

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

• larsim/v06_53_00/source/larsim/LArG4/OpParamAction.h
• larsim/v06_53_00/source/larsim/LArG4/OpParamAction.cxx