CCalRotationMatrixFactory Class Reference

#include <CCalRotationMatrixFactory.hh>

Public Member Functions
	~CCalRotationMatrixFactory ()
G4RotationMatrix *	findMatrix (const G4String &)
G4RotationMatrix *	AddMatrix (const G4String &name, G4double th1, G4double phi1, G4double th2, G4double phi2, G4double th3, G4double phi3)

Static Public Member Functions
static CCalRotationMatrixFactory *	getInstance ()
static CCalRotationMatrixFactory *	getInstance (const G4String &rotfile)
static void	setFileName (const G4String &rotfile)

Detailed Description

Definition at line 45 of file CCalRotationMatrixFactory.hh.

Constructor & Destructor Documentation

CCalRotationMatrixFactory::~CCalRotationMatrixFactory

(

)

Definition at line 84 of file CCalRotationMatrixFactory.cc.

                                                     {
  G4RotationMatrixTableIterator i;
  for(i=theMatrices.begin(); i != theMatrices.end(); ++i) {
    delete (*i).second;
  };
  theMatrices.clear();
}

Member Function Documentation

G4RotationMatrix * CCalRotationMatrixFactory::AddMatrix	(	const G4String &	name,
		G4double	th1,
		G4double	phi1,
		G4double	th2,
		G4double	phi2,
		G4double	th3,
		G4double	phi3
	)

Definition at line 105 of file CCalRotationMatrixFactory.cc.

References python.hepunit::deg, G4cout, G4endl, CLHEP::HepRotation::invert(), and CLHEP::HepRotation::rotateAxes().

Referenced by CCalG4Ecal::constructIn().

                                             {
  G4double sinth1, sinth2,  sinth3, costh1, costh2, costh3;
  G4double sinph1, sinph2, sinph3, cosph1, cosph2, cosph3;
  G4double TH1 = th1/deg, TH2 = th2/deg, TH3 = th3/deg;
  G4double PH1 = phi1/deg, PH2 = phi2/deg, PH3 = phi3/deg;
        
  if (TH1 == 0.0 || TH1 == 360) {
    sinth1 = 0.0; costh1 = 1.0;
  } else if (TH1 == 90.0 || TH1 == -270) {
    sinth1 = 1.0; costh1 = 0.0;
  } else if (TH1 == 180.0 || TH1 == -180.0) {
    sinth1 = 0.0; costh1 = -1.0;
  } else if (TH1 == 270.0 || TH1 == -90.0) {
    sinth1 = -1.0; costh1 = 0.0;
  } else {
    sinth1 = std::sin(th1); costh1 = std::cos(th1);
  }
  
  if (TH2 == 0.0 || TH2 == 360) {
    sinth2 = 0.0; costh2 = 1.0;
  } else if (TH2 == 90.0 || TH2 == -270) {
    sinth2 = 1.0; costh2 = 0.0;
  } else if (TH2 == 180.0 || TH2 == -180.0) {
    sinth2 = 0.0; costh2 = -1.0;
  } else if (TH2 == 270.0 || TH2 == -90.0) {
    sinth2 = -1.0; costh2 = 0.0;
  } else {
    sinth2 = std::sin(th2); costh2 = std::cos(th2);
  }
        
  if (TH3 == 0.0 || TH3 == 360) {
    sinth3 = 0.0; costh3 = 1.0;
  } else if (TH3 == 90.0 || TH2 == -270) {
    sinth3 = 1.0; costh3 = 0.0;
  } else if (TH3 == 180.0 || TH3 == -180.0) {
    sinth3 = 0.0; costh3 = -1.0;
  } else if (TH3 == 270.0 || TH3 == -90.0) {
    sinth3 = -1.0; costh3 = 0.0;
  } else {
    sinth3 = std::sin(th3); costh3 = std::cos(th3);
  }
      
  if (PH1 == 0.0 || PH1 == 360) {
    sinph1 = 0.0; cosph1 = 1.0;
  } else if (PH1 == 90.0 || PH1 == -270) {
    sinph1 = 1.0; cosph1 = 0.0;
  } else if (PH1 == 180.0 || PH1 == -180.0) {
    sinph1 = 0.0; cosph1 = -1.0;
  } else if (PH1 == 270.0 || PH1 == -90.0) {
    sinph1 = -1.0; cosph1 = 0.0;
  } else {
    sinph1 = std::sin(phi1); cosph1 = std::cos(phi1);
  }

  if (PH2 == 0.0 || PH2 == 360) {
    sinph2 = 0.0; cosph2 = 1.0;
  } else if (PH2 == 90.0 || PH2 == -270) {
    sinph2 = 1.0; cosph2 = 0.0;
  } else if (PH2 == 180.0 || PH2 == -180.0) {
    sinph2 = 0.0; cosph2 = -1.0;
  } else if (PH2 == 270.0 || PH2 == -90.0) {
    sinph2 = -1.0; cosph2 = 0.0;
  } else {
    sinph2 = std::sin(phi2); cosph2 = std::cos(phi2);
  }
        
  if (PH3 == 0.0 || PH3 == 360) {
    sinph3 = 0.0; cosph3 = 1.0;
  } else if (PH3 == 90.0 || PH3 == -270) {
    sinph3 = 1.0; cosph3 = 0.0;
  } else if (PH3 == 180.0 || PH3 == -180.0) {
    sinph3 = 0.0; cosph3 = -1.0;
  } else if (PH3 == 270.0 || PH3 == -90.0) {
    sinph3 = -1.0; cosph3 = 0.0;
  } else {
    sinph3 = std::sin(phi3); cosph3 = std::cos(phi3);
  }
                    
  //xprime axis coordinates
  CLHEP::Hep3Vector xprime(sinth1*cosph1,sinth1*sinph1,costh1);
  //yprime axis coordinates
  CLHEP::Hep3Vector yprime(sinth2*cosph2,sinth2*sinph2,costh2);
  //zprime axis coordinates
  CLHEP::Hep3Vector zprime(sinth3*cosph3,sinth3*sinph3,costh3);

#ifdef ddebug
  G4cout << xprime << '\t';    G4cout << yprime << '\t';    G4cout << zprime << G4endl;
#endif
  G4RotationMatrix *rotMat = new G4RotationMatrix();
  rotMat->rotateAxes(xprime, yprime, zprime);
  if (*rotMat == G4RotationMatrix()) {
    // G4cerr << "WARNING: Matrix " << name << " will not be created as a rotation matrix." 
    // G4cerr << "WARNING: Matrix " << name << " is = identity matrix. It will not be created as a rotation matrix." << G4endl;
    delete rotMat;
    rotMat=0;
  } else {
    rotMat->invert();
    theMatrices[name]=rotMat;
#ifdef ddebug
    G4cout << *rotMat << G4endl;
#endif
  }

  return rotMat;
}

G4RotationMatrix * CCalRotationMatrixFactory::findMatrix

(

const G4String &

rot

)

!!Maybe a treatment on not-found case needed.

Definition at line 92 of file CCalRotationMatrixFactory.cc.

Referenced by CCalG4Ecal::constructIn().

                                                                            {
  G4RotationMatrix* retrot=0;
  //Rotation :NULL is no rotation so a null pointer is returned
  if (rot != ":NULL") {
    //retrot untouched if rot not found!!!
    G4RotationMatrixTableIterator it = theMatrices.find(rot);
    if (it != theMatrices.end())
      retrot = (*it).second;
  }
  
  return retrot; //!!!Maybe a treatment on not-found case needed.
}

CCalRotationMatrixFactory * CCalRotationMatrixFactory::getInstance ( )

static

Definition at line 60 of file CCalRotationMatrixFactory.cc.

References G4cerr, and G4endl.

Referenced by CCalDetectorConstruction::Construct(), CCalG4Ecal::constructIn(), and getInstance().

                                                                 {
  if (file=="") {
    G4cerr << "ERROR: You haven't defined which file to use for materials in "
     << "CCalRotationMatrixFactory::getInstance(G4String)" << G4endl;
    return 0;
  }

  if (instance==0) {
    instance = new CCalRotationMatrixFactory;
    return instance;
  }
  else
    return instance;
}

CCalRotationMatrixFactory * CCalRotationMatrixFactory::getInstance ( const G4String &  rotfile )

static

Definition at line 46 of file CCalRotationMatrixFactory.cc.

References G4cerr, G4endl, and getInstance().

                                                                                         {
  if (rotfile=="" || rotfile==file)
    return getInstance();
  else if (file="") {
    file=rotfile;
    return getInstance();
  } else {
    G4cerr << "ERROR: Trying to get Rotation Matrices from " << rotfile 
     << " when previously were retrieved from " << file <<"." << G4endl;
    return 0;
  }
}

void CCalRotationMatrixFactory::setFileName ( const G4String &  rotfile )

static

Definition at line 75 of file CCalRotationMatrixFactory.cc.

References G4cerr, and G4endl.

                                                                   {
  if (rotfile!=file && file!="") {
    G4cerr << "ERROR: Trying to change Rotation Matrices file name to " 
     << rotfile << "." << G4endl;
    G4cerr << "       Using previous file: " << file << G4endl;
  }
  file=rotfile;
}

---

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

• CCalRotationMatrixFactory.hh
• CCalRotationMatrixFactory.cc