hypGeGeoBuilderDEGASStraight40_offset20.C File Reference

#include "TGeoManager.h"

Go to the source code of this file.

Functions
int	hypGeGeoBuilderDEGASStraight40_offset20 (Double_t ZShift=0)

Variables
double	Pi = TMath::Pi()

Function Documentation

int hypGeGeoBuilderDEGASStraight40_offset20

(

Double_t

ZShift = 0

)

Definition at line 5 of file hypGeGeoBuilderDEGASStraight40_offset20.C.

References Double_t, fi, geobuild, geoFace, geoLoad, geom(), gGeoManager, i, Media, outfile, PndGeoHypGeCluster::PlaceCluster(), PndGeoHypGeDEGAS::PrintNodes(), PndGeoHypGeBase::ReadCoordinatesAndAngles(), top, and TString.

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

        TString outfile= "../../../../geometry/hypGeGeoDEGASstraight40_20.root";

  TFile* fi = new TFile(outfile,"RECREATE");

 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);
    

  FairGeoMedium *medcap  = Media->getMedium("HYPaluminium");
  Int_t nmedcap=geobuild->createMedium(medcap);

        TGeoMedium *Ge = gGeoManager->GetMedium("germanium");
  
  
  TGeoMedium *Al = gGeoManager->GetMedium("HYPaluminium");
        
        

  TGeoManager *geom = (TGeoManager*)gROOT->FindObject("FAIRGeom");
  TGeoVolume *top = new TGeoVolumeAssembly("hpGe");
        TGeoVolume *Realtop = new TGeoVolumeAssembly("Realtop");
           
//cout<<" geom "<<geom<<endl;           
  geom->SetTopVolume(top);

        Int_t CrystalNumber = 1;
        Double_t GlobalZOffset = -55;
        Double_t InnerClusterZPosition = -30;
        Double_t OuterClusterZPosition = InnerClusterZPosition + 15;
        
        PndGeoHypGeModifiedBeamPipe *MBP = new PndGeoHypGeModifiedBeamPipe(Al, geom);
        //MBP -> PlaceBeamPipe(Realtop);

        PndGeoHypGeModifiedCTFrame *MCTF = new PndGeoHypGeModifiedCTFrame(Al, geom);
        //MCTF -> PlaceCTFrame(Realtop,new TGeoTranslation(0,0,GlobalZOffset));

        

        PndGeoHypGeDEGAS *DEGAS = new PndGeoHypGeDEGAS(Ge,Al,1);
        
        //Read coordinates and angles from file

        const Int_t NumberOfCluster=20;
        Double_t CoordinatesAndAngles[NumberOfCluster*6];
        string Filename = "DEGASStraightCoordinatesAndAngles40_20.txt";


        DEGAS->ReadCoordinatesAndAngles(Filename,NumberOfCluster,CoordinatesAndAngles);

        for (Int_t i = 0; i < NumberOfCluster*6; i++)
                cout << CoordinatesAndAngles[i] << endl;
        //      Double_t ClusterRadius = 30;

        for (Int_t iCluster = 0; iCluster < 20; iCluster++)
        {

                        
                DEGAS ->PlaceCluster(Realtop,CoordinatesAndAngles[6*(iCluster)],CoordinatesAndAngles[6*(iCluster)+1],CoordinatesAndAngles[6*(iCluster)+2],  GlobalZOffset,CoordinatesAndAngles[6*(iCluster)+3],CoordinatesAndAngles[6*(iCluster)+4],CoordinatesAndAngles[6*(iCluster)+5],&CrystalNumber);
        
                
        }

        
        top->AddNode(Realtop,0);        
        Realtop->PrintNodes();
        geom->CloseGeometry();  
        cout << "# of Nodes: " << Realtop->CountNodes(10,1) << endl;
        DEGAS->PrintNodes(0);
        top->Write();
        fi->Close();
        
        //Realtop->Browse(new TBrowser);
        //c1 = new TCanvas("c1","hypGe",800,600);Realtop->Draw("");//   Realtop->Raytrace();
        Realtop->Draw("ogl");
        

        //c1->x3d();
        //geom->CheckOverlaps(0.0000000001,option="d"); //1 Overlap von Strahlrohr + Kugel ist normal!!!!!!!!!!!!!!!!!
        //geom->PrintOverlaps();
        
        
  return 0;
}

Variable Documentation

double Pi = TMath::Pi()

Definition at line 2 of file hypGeGeoBuilderDEGASStraight40_offset20.C.