// ******** Implementation of class to build HypGe encapsulated HPGe Crystal needed to build cluster detectors *************
//					by steinen@kph.uni-mainz.de

#include "PndGeoHypGeCrystal.h"




ClassImp(PndGeoHypGeCrystal);
Double_t Pi = TMath::Pi();
// -----   Default constructor   -------------------------------------------
PndGeoHypGeCrystal::PndGeoHypGeCrystal(){};


// ----  constructor to use ----------------------------------
PndGeoHypGeCrystal::PndGeoHypGeCrystal(TGeoMedium *ExtGe, TGeoMedium *ExtAl,Int_t ExtCrystalNumber)
{
	Ge = ExtGe;
	Al = ExtAl;
	CrystalNumber = ExtCrystalNumber;

	
//building a capsule
	OuterPgonInnerRadius = 3.0571;
	CapsuleThickness = 0.08;
	InnerScale = (OuterPgonInnerRadius-CapsuleThickness)/OuterPgonInnerRadius;
	CapsuleLength = 10;
	
	// building parts of capsule shape
	CapsuleOuterTube = new TGeoTube("CapsuleOuterTube",0.,3.75,CapsuleLength/2);
	CapsuleOuterPgon = new TGeoPgon("CapsuleOuterPgon",0.,360.,6,2);;
		CapsuleOuterPgon->DefineSection(0,-CapsuleLength/2,OuterPgonInnerRadius,3.8);
		CapsuleOuterPgon->DefineSection(1,4.9094,3.75,3.8);
	CapsuleInnerTube = new TGeoTube("CapsuleInnerTube",0,3.75*InnerScale,5);;
	CapsuleInnerPgon = new TGeoPgon("CapsuleInnerPgon",0.,360.,6,2);
		CapsuleInnerPgon->DefineSection(0,-CapsuleLength/2,OuterPgonInnerRadius-CapsuleThickness,3.8);
		CapsuleInnerPgon->DefineSection(1,4.9094,3.75-CapsuleThickness,3.8);
	CapsuleTopCover = new TGeoTube("CapsuleTopCover",0,3.75,CapsuleThickness/2);
	CapsuleTopZTranslation = new TGeoTranslation("CapsuleTopZTranslation",0,0,5-CapsuleThickness/2);
		CapsuleTopZTranslation->RegisterYourself();

	//combining simple shapes
	CapsuleOuterShape = new TGeoCompositeShape("CapsuleOuterShape","(CapsuleOuterTube-CapsuleOuterPgon)");
	CapsuleInnerShape = new TGeoCompositeShape("CapsuleInnerShape","(CapsuleInnerTube-CapsuleInnerPgon)");
	CapsuleInnerZTranslation = new TGeoTranslation("CapsuleInnerZTranslation",0,0,CapsuleThickness);
		CapsuleInnerZTranslation->RegisterYourself();
	CapsuleCompleteShape = new TGeoCompositeShape("CapsuleCompleteShape", "CapsuleOuterShape-CapsuleInnerShape:CapsuleInnerZTranslation+CapsuleTopCover:CapsuleTopZTranslation");

	//building capsule volume
	CapsuleName= new char[11];
	sprintf(CapsuleName,"Capsule",(int)CrystalNumber);
	Capsule = new TGeoVolume(CapsuleName,CapsuleCompleteShape,Al);
	Capsule->SetTransparency(50);
	Capsule->SetLineColor(19);

	
//building a crystal

	CapsuleThickness = 0.08;		
	VacuumThickness = 0.08;
	CrystalLength = 7.85;
	CrystalFrontPgonInnerRadius = 2.9485;
	
	//building the main shape of the crystal
	CrystalTube = new TGeoTube("CrystalTube",0,3.5,CrystalLength/2);
	CrystalPgon = new TGeoPgon("CrystalPgon",0.,360.,6,2);
		CrystalPgon->DefineSection(0,-CrystalLength/2,CrystalFrontPgonInnerRadius,3.8);
		CrystalPgon->DefineSection(1,8-CrystalLength/2,CrystalFrontPgonInnerRadius+8*TMath::Tan(Pi/180*4.12),3.8);

	//building the inner hole of the crystal
	CrystalInnerTube = new TGeoTube("CrystalInnerTube",0,0.6,3.05);
	CrystalInnerTubeZTranslation = new TGeoTranslation("CrystalInnerTubeZTranslation",0,0,CrystalLength/2-3.05);
		CrystalInnerTubeZTranslation->RegisterYourself();
	CrystalInnerTubeTorus = new TGeoTorus("CrystalInnerTubeTorus",0.2,0,0.4,0,360);
	CrystalInnerTubeTorusZTranslation = new TGeoTranslation("CrystalInnerTubeTorusZTranslation",0,0,CrystalLength/2-3.05*2);
		CrystalInnerTubeTorusZTranslation->RegisterYourself();
	CrystalInnerTubeMiniExtraTube = new TGeoTube("CrystalInnerTubeMiniExtraTube",0,0.2,0.2);
	CrystalInnerTubeMiniExtraTubeZTranslation = new TGeoTranslation("CrystalInnerTubeMiniExtraTubeZTranslation",0,0,CrystalLength/2-3.05*2-0.2);
		CrystalInnerTubeMiniExtraTubeZTranslation ->RegisterYourself();

	//building the top rounding of the crystal
	CrystalTopRoundingTube = new TGeoTube("CrystalTopRoundingTube",2.9,3.5,0.3);
	CrystalTopRoundingTubeZTranslation = new TGeoTranslation("CrystalTopRoundingTubeZTranslation",0,0,0.3);
		CrystalTopRoundingTubeZTranslation->RegisterYourself();
	CrystalTopRoundingTorus = new TGeoTorus("CrystalTopRoundingTorus",2.9,0,0.6,0,360);
	CrystalTopRoundingTooling = new TGeoCompositeShape ("CrystalTopRoundingTooling"	,"CrystalTopRoundingTube:CrystalTopRoundingTubeZTranslation-CrystalTopRoundingTorus");
	CrystalTopRoundingToolingZTranslation = new TGeoTranslation("CrystalTopRoundingToolingZTranslation",0,0,CrystalLength/2-0.6);
		CrystalTopRoundingToolingZTranslation->RegisterYourself();

	//building the guard ring of the crystal
	CrystalGuardRingTube = new TGeoTube("CrystalGuardRingTube",0.9,1.4,0.15);
	CrystalGuardRingTubeZTranslation = new TGeoTranslation("CrystalGuardRingTubeZTranslation",0,0,CrystalLength/2-0.15);
		CrystalGuardRingTubeZTranslation->RegisterYourself();
	CrystalGuardRingTorus = new TGeoTorus("CrystalGuardRingTorus",1.2,0,0.2,0,360);
	CrystalGuardRingTorusInner = new TGeoTorus("CrystalGuardRingTorusInner",1.1,0,0.2,0,360);
	CrystalGuardRingTorusZTranslation = new TGeoTranslation("CrystalGuardRingTorusZTranslation",0,0,CrystalLength/2-0.15*2);
		CrystalGuardRingTorusZTranslation->RegisterYourself();
	CrystalGuardRingMiniExtraTube = new TGeoTube("CrystalGuardRingMiniExtraTube",1.1,1.2,0.2);
	CrystalGuardRingMiniExtraTubeZTranslation = new TGeoTranslation("CrystalGuardRingMiniExtraTubeZTranslation",0,0,CrystalLength/2-0.15*2);
		CrystalGuardRingMiniExtraTubeZTranslation ->RegisterYourself();

	//combining the final shape of the crystal		
	CrystalShape= new TGeoCompositeShape("CrystalShape","(CrystalTube-CrystalPgon)-CrystalTopRoundingTooling:CrystalTopRoundingToolingZTranslation-(CrystalInnerTube:CrystalInnerTubeZTranslation+CrystalInnerTubeTorus:CrystalInnerTubeTorusZTranslation+CrystalInnerTubeMiniExtraTube:CrystalInnerTubeMiniExtraTubeZTranslation)-(CrystalGuardRingTube:CrystalGuardRingTubeZTranslation+CrystalGuardRingTorus:CrystalGuardRingTorusZTranslation+CrystalGuardRingTorusInner:CrystalGuardRingTorusZTranslation+CrystalGuardRingMiniExtraTube:CrystalGuardRingMiniExtraTubeZTranslation)"); //	(main shape) - (inner tube) - (top rounding) - (guard ring)

	CrystalName= new char[11];
	sprintf(CrystalName,"Crystal",CrystalNumber);
	Crystal = new TGeoVolume(CrystalName,CrystalShape,Ge);
	Crystal->SetLineColor(kBlue);
	Crystal->SetFillColor(kBlue);

	
//combining capsule and crystal

	CrystalPlusCapsule = new TGeoVolumeAssembly("CrystalPlusCapsule");
  CapsuleCombiTrans = new TGeoCombiTrans(0,0,-CapsuleLength/2, new TGeoRotation("CapsuleRotation",0,180,0));
		CapsuleCombiTrans->RegisterYourself();
  CrystalInsideCapsuleCombiTrans = new TGeoCombiTrans(0,0,-(CrystalLength/2+CapsuleThickness+VacuumThickness), new TGeoRotation("CrystalRotation",0,180,0));
		CrystalInsideCapsuleCombiTrans->RegisterYourself();

	CrystalPlusCapsule->AddNode(Capsule, CrystalNumber, CapsuleCombiTrans);
	CrystalPlusCapsule->AddNode(Crystal,CrystalNumber,CrystalInsideCapsuleCombiTrans); 


};

Int_t	PndGeoHypGeCrystal::GetCrystalNumber()
{
	return CrystalNumber;
}
Int_t PndGeoHypGeCrystal::GetNdaughters()
{
	CrystalPlusCapsule->GetNdaughters();
};


void PndGeoHypGeCrystal::PrintNodes()
{
	CrystalPlusCapsule->PrintNodes();
}

void PndGeoHypGeCrystal::PlaceCrystal(TGeoVolume *top, TGeoMatrix* CrystalPlaceAndDirectionTranslation)
{
	top->AddNode(CrystalPlusCapsule,CrystalNumber,CrystalPlaceAndDirectionTranslation);
};