PndGeoHypGeSingleCluster.cxx

Go to the documentation of this file.

// ******** Implementation of class to build HypGe Double Cluster Geometry *************
//                                      by steinen@kph.uni-mainz.de
#include "PndGeoHypGeSingleCluster.h"
#include "PndGeoHypGeCrystal.h"
#include <iostream>


ClassImp(PndGeoHypGeSingleCluster)

using namespace std;
//double Pi = TMath::Pi();
// -----   Default constructor   -------------------------------------------
PndGeoHypGeSingleCluster::PndGeoHypGeSingleCluster()
{
  // Constructor
}
// -----    constructor to use   -------------------------------------------
PndGeoHypGeSingleCluster::PndGeoHypGeSingleCluster(TGeoMedium *ExtGe, TGeoMedium *ExtAl,TGeoMedium *ExtCu ,Int_t ExtClusterNumber)
{
        MakeTop = true;

        Ge = ExtGe;
        Al = ExtAl;
        Cu = ExtCu;
        ClusterNumber= ExtClusterNumber;
        StartClusterNumber = ExtClusterNumber;
                                        NumberOfCryostatEdges =10;
                                        NumberOfDifferentZValues = 5;
                                        NumberOfXYInput = (NumberOfDifferentZValues*2-1)*2;
                                        NumberOfFrontArb8 = 4;
                                        NumberOfTopArb8 = NumberOfFrontArb8;
        char* InputTitleBuffer;
        ZValuesInputArray = new Double_t[NumberOfDifferentZValues];
        CoordinatesInputArray = new Double_t[NumberOfCryostatEdges*NumberOfXYInput];
        

        CryostatFrontThickness=  0.1;   //ZValuesInputArray[1]-ZValuesInputArray[0];
        CryostatSideHeight = 12;
        CryostatBack1Height = 1;
        CryostatBack2Height = 20;
        ColdFinger1Height = 10;
        ColdFinger1Height = 10;
}



// ------- destructor --------------
PndGeoHypGeSingleCluster::~PndGeoHypGeSingleCluster()
{
 cout << "Destrucor called" << endl;
};

// --------------- build the crystals ---------------------------

void PndGeoHypGeSingleCluster::BuildCrystals(Int_t *CrystalNumber)
{
        //create 2 crystals. The copy number of a crystal is composed of the number of the cluster and a running number for every placed crystal like e.g.: cluster 4, crystal 7: 407 ; cluster 11, crystal 21: 1121
        Crystal1Number = *CrystalNumber;
        *CrystalNumber+=1;
        cout << "Cluster:" << ClusterNumber << "\tCrystal1Number: "<< Crystal1Number<<endl;

        Crystal1 = new PndGeoHypGeCrystal(Ge,Al,100*ClusterNumber+Crystal1Number);
        


        //set positions of the 2 crystals

        
        SingleCrystalZOffset = 0.3;                             //measured in CAD
        
        Crystal1Rotation = new TGeoRotation("Crystal1Rotation", 0,0,0);                 //angles measured in CAD drawing
                Crystal1Rotation ->RegisterYourself();
        

        //maybe some signs needs adaption
        Crystal1Position = new TGeoCombiTrans(0,0,-SingleCrystalZOffset,Crystal1Rotation);
                Crystal1Position->RegisterYourself();
        
        

        Crystal1->PlaceCrystal(ClusterAssembly,Crystal1Position);
};

// ---------------      build the cryostat --------------------------
void PndGeoHypGeSingleCluster::BuildCryostat()
{
        // the cryostat is a combination front, side and top parts. The single crystal detector is easy and consists of only 2 +x parts

        string CompositeShapeExpression;

        
        //cout << CompositeShapeExpression << endl;
        

        CryostatFront = new TGeoTube("CryostatFront",0,6,CryostatFrontThickness/2);
        CryostatSide = new TGeoTube ("CryostatSide", 6,6.2,CryostatSideHeight/2);
        CryostatBack1 = new TGeoTube ("CryostatBack1", 4.2,6,CryostatBack1Height/2);
        CryostatBack2 = new TGeoTube ("CryostatBack2", 4,4.2,CryostatBack2Height/2);

        
        CryostatFrontTranslation = new TGeoTranslation("CryostatFrontTranslation",0,0,CryostatFrontThickness/2);
        CryostatFrontTranslation->RegisterYourself();
        CryostatSideTranslation = new TGeoTranslation("CryostatSideTranslation",0,0,CryostatSideHeight/2);
        CryostatSideTranslation->RegisterYourself();
        
        CryostatBack1Translation = new TGeoTranslation("CryostatBack1Translation",0,0,CryostatSideHeight+CryostatBack1Height/2);
        CryostatBack1Translation->RegisterYourself();
        CryostatBack2Translation = new TGeoTranslation("CryostatBack2Translation",0,0,CryostatSideHeight+CryostatBack2Height/2);
        CryostatBack2Translation->RegisterYourself();
        

        
        CompositeShapeExpression = "CryostatFront : CryostatFrontTranslation + CryostatSide : CryostatSideTranslation + CryostatBack1 : CryostatBack1Translation + CryostatBack2 : CryostatBack2Translation";
        
        //cout << CompositeShapeExpression << endl;
        //combine the cryostat
        CryostatShape = new TGeoCompositeShape("CryostatShape",CompositeShapeExpression.data() ); 
        Cryostat = new TGeoVolume("Cryostat",CryostatShape,Al);
        Cryostat->SetLineColor(kGreen);
        Cryostat->SetFillColor(kGreen);
        
        // add the coldfinger
        ColdFinger1 = new TGeoTube ("ColdFinger1", 0,1.5,ColdFinger1Height/2);
        ColdFinger2 = new TGeoTube ("ColdFinger2", 0,1.5,ColdFinger2Height/2);

        ColdFinger1Translation= new TGeoTranslation("ColdFinger1Translation",0,0,-SingleCrystalZOffset+10+ColdFinger1Height/2);
                ColdFinger1Translation->RegisterYourself();
        ColdFinger2Translation= new TGeoTranslation("ColdFinger2Translation",0,0,-SingleCrystalZOffset+10+ColdFinger1Height+ColdFinger2Height/2);
                ColdFinger2Translation->RegisterYourself();

        string ColdFingerShapeExpression = "ColdFinger1 : ColdFinger1Translation + ColdFinger2 : ColdFinger2Translation";                       // problem here? no trans if only one part?
        ColdFingerShape = new TGeoCompositeShape("ColdFingerShape",ColdFingerShapeExpression.data() ); 
        ColdFinger = new TGeoVolume("ColdFinger", ColdFingerShape,Cu);
                        //bring the Cryostat in the right position
        CryostatCombiTrans = new TGeoCombiTrans("CryostatCombiTrans",0,0,0,new TGeoRotation("CryostatRotation",0,180,0));
        CryostatCombiTrans->RegisterYourself();
        
        ClusterAssembly->AddNode(Cryostat,ClusterNumber,CryostatCombiTrans);
        ClusterAssembly->AddNode(ColdFinger,ClusterNumber,CryostatCombiTrans);

};

// ---------------  Prints the nodes of a cluster -----------------
void PndGeoHypGeSingleCluster::PrintNodes(Int_t nLevels )
{
        if (nLevels == 0)
        {
                cout << "Nothing to Print, use higher number to go deeper in the tree." << endl;
        }
        if (nLevels > 0)
        {
                cout << "Printing the content of the cluster assembly." << endl;
                cout <<endl << "Number of Nodes in cluster assembly: " << ClusterAssembly->GetNdaughters() << endl;
                ClusterAssembly->PrintNodes();          
        }
        if (nLevels > 1)
        {
                cout << "Printing the content of the crystal assemblies." << endl;
                cout <<endl << "Number of Nodes in crystal assembly 1: " << Crystal1->GetNdaughters() << endl;
                Crystal1->PrintNodes();
                cout <<endl << "Number of Nodes in crystal assembly 2: " << Crystal2->GetNdaughters() << endl;
                Crystal2->PrintNodes();
        }
}