CrystalShapeTest.C

Go to the documentation of this file.

#include "TGeoManager.h"
double Pi = TMath::Pi();



int CrystalShapeTest()
{

gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C");
  
  // Load this libraries
  gSystem->Load("libGeoBase");
  gSystem->Load("libParBase");
  gSystem->Load("libBase");
  gSystem->Load("libPndData");
  gSystem->Load("libPassive");


 FairGeoLoader* geoLoad = new FairGeoLoader("TGeo","FairGeoLoader");
  FairGeoInterface *geoFace = geoLoad->getGeoInterface();
  geoFace->setMediaFile("../../../../geometry/media_pnd.geo");
  geoFace->readMedia();
  geoFace->print();

  FairGeoMedia *Media =  geoFace->getMedia();
  FairGeoBuilder *geobuild=geoLoad->getGeoBuilder();

  FairGeoMedium *medGe  = Media->getMedium("germanium");
  Int_t nmedGe=geobuild->createMedium(medGe);
    
  // aluminum capsule
  FairGeoMedium *medcap  = Media->getMedium("HYPaluminium");
  Int_t nmedcap=geobuild->createMedium(medcap);

        TGeoMedium *Ge = gGeoManager->GetMedium("germanium");
  
  //capsule
  TGeoMedium *Al_sol = gGeoManager->GetMedium("HYPaluminium");


  TGeoManager *geom = (TGeoManager*)gROOT->FindObject("FAIRGeom");
  TGeoVolume *top = new TGeoVolumeAssembly("hpGe");
        TGeoVolume *strange = new TGeoVolumeAssembly("strange");
           
cout<<" geom "<<geom<<endl;           
  geom->SetTopVolume(top);
        
        Double_t GlobalZOffset = -55;
        int ObjektNr = 0;
        
        Double_t                                                CapsuleThickness;
        Double_t                                                VacuumThickness;
        Double_t                                                CrystalLenght;
        Double_t                                                CrystalFrontPgonInnerRadius;
        
        Double_t                                                CrystalInnerPgonRadius;
        //objects to build the main crystal shape
        TGeoTube                                                *CrystalTube;
        TGeoPgon                                                *CrystalPgon;
        //objects to build the inner crystal hole
        TGeoTube                                                *CrystalInnerTube;
        TGeoTranslation                 *CrystalInnerTubeZTranslation;
        TGeoTorus                                               *CrystalInnerTubeTorus;
        TGeoTranslation                 *CrystalInnerTubeTorusZTranslation;
        TGeoTube                                                *CrystalInnerTubeMiniExtraTube;
        TGeoTranslation                 *CrystalInnerTubeMiniExtraTubeZTranslation;
        //objects to get the top rounding of the crystal
        TGeoTube                                                *CrystalTopRoundingTube;
        TGeoTranslation                 *CrystalTopRoundingTubeZTranslation;
        TGeoTorus                                               *CrystalTopRoundingTorus;
        TGeoCompositeShape      *CrystalTopRoundingTooling;
        TGeoTranslation                 *CrystalTopRoundingToolingZTranslation;

        //objects to build the guard ring of the crystal
        TGeoTube                                                *CrystalGuardRingTube;
        TGeoTranslation                 *CrystalGuardRingTubeZTranslation;
        TGeoTorus                                               *CrystalGuardRingTorus;
        TGeoTorus                                               *CrystalGuardRingTorusInner;
        TGeoTranslation                 *CrystalGuardRingTorusZTranslation;
        TGeoTube                                                *CrystalGuardRingMiniExtraTube;
        TGeoTranslation                 *CrystalGuardRingMiniExtraTubeZTranslation;

        //building the whole crystal
        TGeoCompositeShape *CrystalShape;
        TGeoVolume *Crystal;

        
        


        CapsuleThickness = 0.08;
        VacuumThickness = 0.08;
        CrystalLenght = 7.85;
        CrystalFrontPgonInnerRadius = 2.9485;
        
        //building the main shape of the crystal
        CrystalTube = new TGeoTube("CrystalTube",0,3.5,CrystalLenght/2);
        CrystalPgon = new TGeoPgon("CrystalPgon",0.,360.,6,2);
                CrystalPgon->DefineSection(0,-CrystalLenght/2,CrystalFrontPgonInnerRadius,3.8);
                CrystalPgon->DefineSection(1,8-CrystalLenght/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,CrystalLenght/2-3.05);
                CrystalInnerTubeZTranslation->RegisterYourself();
        CrystalInnerTubeTorus = new TGeoTorus("CrystalInnerTubeTorus",0.2,0,0.4,0,360);
        CrystalInnerTubeTorusZTranslation = new TGeoTranslation("CrystalInnerTubeTorusZTranslation",0,0,CrystalLenght/2-3.05*2);
                CrystalInnerTubeTorusZTranslation->RegisterYourself();
        CrystalInnerTubeMiniExtraTube = new TGeoTube("CrystalInnerTubeMiniExtraTube",0,0.2,0.2);
        CrystalInnerTubeMiniExtraTubeZTranslation = new TGeoTranslation("CrystalInnerTubeMiniExtraTubeZTranslation",0,0,CrystalLenght/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,CrystalLenght/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,CrystalLenght/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,CrystalLenght/2-0.15*2);
                CrystalGuardRingTorusZTranslation->RegisterYourself();
        CrystalGuardRingMiniExtraTube = new TGeoTube("CrystalGuardRingMiniExtraTube",1.1,1.2,0.2);
        CrystalGuardRingMiniExtraTubeZTranslation = new TGeoTranslation("CrystalGuardRingMiniExtraTubeZTranslation",0,0,CrystalLenght/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)
        Crystal = new TGeoVolume("Crystal",CrystalShape,Ge);
        Crystal->SetLineColor(kGreen);
        Crystal->SetFillColor(kGreen);

        
//TGeoTranslation *CryTrans = new TGeoTranslation(0,0,5-1.5+CapsuleThickness+VacuumThickness);
//CryTrans->RegisterYourself();
strange->AddNode(Crystal,6);



        top->AddNode(strange,0);
        geom->CloseGeometry();

        
        //strange->Draw("");    strange->Raytrace();
        strange->Draw("ogl");
        
        
  return 0;
}