#ifndef HGEOPCON_H
#define HGEOPCON_H

#include "hgeobasicshapes.h"

//
// class for Geant-shape PCON
// derived from baseclass HGeoBasicShapes
//
// The functions readPoints(...) read from database or from file and store the
// input in the float-vector vol of class HGeoParam.
// ( point 0:   number of planes perpendicular to the z axis where the dimension 
//              the section is given;
//   point 1:   azimutal angle phi at which the volume begins,
//              opening angle dphi of the volume,;
//   point 2ff: z coordinate of the section,
//              inner radius at position z,
//              outer radius at position z;
//
// The technical coordinate system of a PCON, which sits in CAVE and is not
// rotated, is the laboratory system.
// It's the same as in Geant and therefore the matrix cgRot, which describes
// the transformation between these coordinate systems is the unity matrix.
//

class HGeoPcon : public HGeoBasicShapes {
public:
	HGeoPcon() {
                cgRot=uniRot;
        }
	~HGeoPcon(){}

	// reads input from database
	int readPointsDb(HGeoInfo & geoInfo, HGeoParam & ele);

	// reads input from file
	int readPoints(HGeoInfo & geoInfo, HGeoParam & ele);

	// calculates the Geant-shapeparameters
	void calcShape(HGeoParam & ele, HGeoShapePar & shapePar);

        // calculates the position the Geant coordinate system in the technical
        // coordinate system and invokes posInMother(...) of the base class
	int posShape(HGeoParam & ele, HGeoParam & mo, HGeoShapePar & shapePar);

        // writes parameters to file
        int writePoints(ofstream & fout, HGeoParam & ele);
};

#endif