PHG4EtaPhiParameterization.h

Go to the documentation of this file.

#ifndef __PHG4ETAPHIPARAMETERIZATION_H__
#define __PHG4ETAPHIPARAMETERIZATION_H__
 
#include <vector>
#include <Geant4/globals.hh>
#include <Geant4/G4VPVParameterisation.hh>
 
class G4VPhysicalVolume;
class G4Box;
 
// Parameterization to define rings whose size changes with Z to 
// correspond to fixed width in eta.
 
class PHG4EtaPhiParameterization : public G4VPVParameterisation
{ 
public:
  
  PHG4EtaPhiParameterization(
                          unsigned int neta,       // Binning in eta
                          double minEta,  // "
                          double maxEta,  // "
                          unsigned int nphi,
                          double startPhi,
                          double deltaPhi,
                          double radiusIn,  // Radius of inner face of cylinder
                          double radiusOut,  // Radius of outer face of cylinder
                          double centerZ  // overall Z of center of rings
                          );
  
  virtual ~PHG4EtaPhiParameterization();
  
  virtual void Print(std::ostream& os = std::cout) const;
 
  void ComputeTransformation(const G4int copyNo,
                             G4VPhysicalVolume* physVol) const;
  
  void ComputeDimensions(G4Tubs& ring, const G4int copyNo,
                         const G4VPhysicalVolume* physVol) const;
 
  int GetIEta(int copyNo) const { return _ieta.at(copyNo); }
  int GetIPhi(int copyNo) const { return _iphi.at(copyNo); }
 
private:  // Dummy declarations to get rid of warnings ...
 
//   void ComputeDimensions(G4Trd&,const G4int,const G4VPhysicalVolume*) const {}
//   void ComputeDimensions(G4Trap&,const G4int,const G4VPhysicalVolume*) const {}
//   void ComputeDimensions(G4Cons&,const G4int,const G4VPhysicalVolume*) const {}
//   void ComputeDimensions(G4Sphere&,const G4int,const G4VPhysicalVolume*) const {}
//   void ComputeDimensions(G4Orb&,const G4int,const G4VPhysicalVolume*) const {}
//   void ComputeDimensions(G4Torus&,const G4int,const G4VPhysicalVolume*) const {}
//   void ComputeDimensions(G4Para&,const G4int,const G4VPhysicalVolume*) const {}
//   void ComputeDimensions(G4Hype&,const G4int,const G4VPhysicalVolume*) const {}
//   void ComputeDimensions(G4Box&,const G4int,const G4VPhysicalVolume*) const {}
//   void ComputeDimensions(G4Polycone&,const G4int,const G4VPhysicalVolume*) const {}
//   void ComputeDimensions(G4Polyhedra&,const G4int,const G4VPhysicalVolume*) const {}
 
private:
 
  unsigned int _neta;   
  double _minEta;
  double _maxEta;
  unsigned int _nphi;
  double _startPhi;
  double _deltaPhi;
  double _radiusIn;
  double _radiusOut;
  double _centerZ;
  std::vector<double> _zpos; // Z positions of the rings
  std::vector<double> _zhalf; // Z half-widths of the rings
  std::vector<double> _phi0; // Lower edge of phi bins
  std::vector<double> _phi1; // Upper edge of phi bins
  std::vector<int> _ieta;
  std::vector<int> _iphi;
};
 
 
#endif