PndGeoHypGeDoubleCluster.cxx

Go to the documentation of this file.

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


ClassImp(PndGeoHypGeDoubleCluster)

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

        Ge = ExtGe;
        Al = ExtAl;
        ClusterNumber= ExtClusterNumber;
        StartClusterNumber = ExtClusterNumber;
        NumberOfCryostatEdges =10;
        NumberOfDifferentZValues = 5;
        NumberOfXYInput = (NumberOfDifferentZValues*2-1)*2;
        NumberOfFrontArb8 = 4;
        NumberOfTopArb8 = NumberOfFrontArb8;
        //ClusterAssembly = new TGeoVolumeAssembly("ClusterAssembly");          //z = 0 only once, every other z value twice (outer and inner); always x and y coordinate

        //read coordinates from file
        char* InputTitleBuffer;
        ZValuesInputArray = new Double_t[NumberOfDifferentZValues];
        CoordinatesInputArray = new Double_t[NumberOfCryostatEdges*NumberOfXYInput];
        string Filename = getenv("VMCWORKDIR");
        Filename +="/hypGe/geometry/PndGeoHypGeDoubleClusterParfile.txt";

        ParametersInputFile.open(Filename.data());
        FillZValues(NumberOfDifferentZValues);
        FillInputArray(NumberOfCryostatEdges,NumberOfXYInput);
        ParametersInputFile.close();

        CryostatFrontThickness=  ZValuesInputArray[1]-ZValuesInputArray[0];
        CryostatLowerArb8Heigt = ZValuesInputArray[2]-ZValuesInputArray[1];
        CryostatMiddleArb8Heigt= ZValuesInputArray[3]-ZValuesInputArray[2];
        CryostatUpperArb8Heigt= ZValuesInputArray[4]-ZValuesInputArray[3];
        CryostatTopThickness=0.5;
}



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

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

void PndGeoHypGeDoubleCluster::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;
        Crystal2Number = *CrystalNumber;
        *CrystalNumber+=1;
        cout << "Cluster:" << ClusterNumber << "\tCrystal1Number: "<< Crystal1Number<<endl;
        cout << "Cluster:" << ClusterNumber << "\tCrystal2Number: "<< Crystal2Number<<endl;

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


        //set positions of the 2 crystals

        DoubleCrystalYOffset = 6.1141/2;                        //measured in CAD
        DoubleCrystalZOffset = 0.7704;                          //measured in CAD
        
        Crystal1Rotation = new TGeoRotation("Crystal1Rotation", 0,-4.125,0);                    //angles measured in CAD drawing
                Crystal1Rotation ->RegisterYourself();
        Crystal2Rotation = new TGeoRotation("Crystal2Rotation", 0,4.125,0);
                Crystal2Rotation ->RegisterYourself();

        //maybe some signs needs adaption
        Crystal1Position = new TGeoCombiTrans(0,-DoubleCrystalYOffset,-DoubleCrystalZOffset,Crystal1Rotation);
                Crystal1Position->RegisterYourself();
        Crystal2Position = new TGeoCombiTrans(0,+DoubleCrystalYOffset,-DoubleCrystalZOffset,Crystal2Rotation);
                Crystal2Position->RegisterYourself();
        

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

// ---------------      build the cryostat --------------------------
void PndGeoHypGeDoubleCluster::BuildCryostat()
{
        // the cryostat is a combination front, side and top parts. The front and the top is build of 4 TGeoArb8. The side parts are split into 10 lower, 10 middle and 10 upper parts. All are TGeoArb8. This parts are all combined to a TGeoCompositeShape

        string CompositeShapeExpression;
//the front of the cryostat

        CryostatFrontArb8XY_1 = new Double_t[16];
                BuildTArb8FrontArray(CryostatFrontArb8XY_1, 1,2,3,10, NumberOfXYInput);
        CryostatFrontArb8XY_2 = new Double_t[16];
                BuildTArb8FrontArray(CryostatFrontArb8XY_2, 3,4,5,6, NumberOfXYInput);
        CryostatFrontArb8XY_3 = new Double_t[16];
                BuildTArb8FrontArray(CryostatFrontArb8XY_3, 3,6,7,8, NumberOfXYInput);
        CryostatFrontArb8XY_4 = new Double_t[16];
                BuildTArb8FrontArray(CryostatFrontArb8XY_4, 3,8,9,10, NumberOfXYInput);

        CryostatFrontArb8_1 = new TGeoArb8 ("CryostatDoubleFrontArb8_1",CryostatFrontThickness/2,CryostatFrontArb8XY_1);
        CryostatFrontArb8_2 = new TGeoArb8 ("CryostatDoubleFrontArb8_2",CryostatFrontThickness/2,CryostatFrontArb8XY_2);
        CryostatFrontArb8_3 = new TGeoArb8 ("CryostatDoubleFrontArb8_3",CryostatFrontThickness/2,CryostatFrontArb8XY_3);
        CryostatFrontArb8_4 = new TGeoArb8 ("CryostatDoubleFrontArb8_4",CryostatFrontThickness/2,CryostatFrontArb8XY_4);

        CryostatFrontTranslation = new TGeoTranslation("CryostatDoubleFrontTranslation",0,0,CryostatFrontThickness/2);
        CryostatFrontTranslation->RegisterYourself();

        for (Int_t i = 1; i < NumberOfFrontArb8+1; i++)
        {
                char buffer[100];
                sprintf(buffer, "CryostatDoubleFrontArb8_%d : CryostatDoubleFrontTranslation",i);
                CompositeShapeExpression += buffer;
                //if (i<NumberOfFrontArb8)      //uncomment if last part of composite shape
                {
                        CompositeShapeExpression += "+";
                }
        }
                
//the lower part of the cryostat

        CryostatLowerArb8XY_1 = new Double_t[16];
                        BuildTArb8Array(CryostatLowerArb8XY_1, 1,2,1,2, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatLowerArb8XY_2 = new Double_t[16];
                BuildTArb8Array(CryostatLowerArb8XY_2, 2,3,1,2, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatLowerArb8XY_3 = new Double_t[16];
                BuildTArb8Array(CryostatLowerArb8XY_3, 3,4,1,2, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatLowerArb8XY_4 = new Double_t[16];
                BuildTArb8Array(CryostatLowerArb8XY_4, 4,5,1,2, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatLowerArb8XY_5 = new Double_t[16];
                BuildTArb8Array(CryostatLowerArb8XY_5, 5,6,1,2, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatLowerArb8XY_6 = new Double_t[16];
                BuildTArb8Array(CryostatLowerArb8XY_6, 6,7,1,2, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatLowerArb8XY_7 = new Double_t[16];
                BuildTArb8Array(CryostatLowerArb8XY_7, 7,8,1,2, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatLowerArb8XY_8 = new Double_t[16];
                BuildTArb8Array(CryostatLowerArb8XY_8, 8,9,1,2, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatLowerArb8XY_9 = new Double_t[16];
                BuildTArb8Array(CryostatLowerArb8XY_9, 9,10,1,2, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatLowerArb8XY_10 = new Double_t[16];
                BuildTArb8Array(CryostatLowerArb8XY_10, 10,1,1,2, NumberOfDifferentZValues-1, NumberOfXYInput);

                

        CryostatLowerArb8_1 = new TGeoArb8 ("CryostatDoubleLowerArb8_1",CryostatLowerArb8Heigt/2,CryostatLowerArb8XY_1);
        CryostatLowerArb8_2 = new TGeoArb8 ("CryostatDoubleLowerArb8_2",CryostatLowerArb8Heigt/2,CryostatLowerArb8XY_2);
        CryostatLowerArb8_3 = new TGeoArb8 ("CryostatDoubleLowerArb8_3",CryostatLowerArb8Heigt/2,CryostatLowerArb8XY_3);
        CryostatLowerArb8_4 = new TGeoArb8 ("CryostatDoubleLowerArb8_4",CryostatLowerArb8Heigt/2,CryostatLowerArb8XY_4);
        CryostatLowerArb8_5 = new TGeoArb8 ("CryostatDoubleLowerArb8_5",CryostatLowerArb8Heigt/2,CryostatLowerArb8XY_5);
        CryostatLowerArb8_6 = new TGeoArb8 ("CryostatDoubleLowerArb8_6",CryostatLowerArb8Heigt/2,CryostatLowerArb8XY_6);
        CryostatLowerArb8_7 = new TGeoArb8 ("CryostatDoubleLowerArb8_7",CryostatLowerArb8Heigt/2,CryostatLowerArb8XY_7);
        CryostatLowerArb8_8 = new TGeoArb8 ("CryostatDoubleLowerArb8_8",CryostatLowerArb8Heigt/2,CryostatLowerArb8XY_8);
        CryostatLowerArb8_9 = new TGeoArb8 ("CryostatDoubleLowerArb8_9",CryostatLowerArb8Heigt/2,CryostatLowerArb8XY_9);
        CryostatLowerArb8_10 = new TGeoArb8 ("CryostatDoubleLowerArb8_10",CryostatLowerArb8Heigt/2,CryostatLowerArb8XY_10);

        CryostatLowerTranslation = new TGeoTranslation("CryostatDoubleLowerTranslation",0,0,CryostatFrontThickness+CryostatLowerArb8Heigt/2);
        CryostatLowerTranslation->RegisterYourself();

        for (Int_t i = 1; i < NumberOfCryostatEdges+1; i++)
        {
                char buffer[100];
                sprintf(buffer, "CryostatDoubleLowerArb8_%d : CryostatDoubleLowerTranslation",i);
                CompositeShapeExpression += buffer;
                //if (i<NumberOfCryostatEdges)          //uncomment if last part of composite shape
                {
                        CompositeShapeExpression += "+";
                }
        }

//the middle part of the cryostat

        CryostatMiddleArb8XY_1 = new Double_t[16];
                BuildTArb8Array(CryostatMiddleArb8XY_1, 1,2,2,3, NumberOfDifferentZValues-1, NumberOfXYInput);          // from PndGeoHypGeBase
        CryostatMiddleArb8XY_2 = new Double_t[16];
                BuildTArb8Array(CryostatMiddleArb8XY_2, 2,3,2,3, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatMiddleArb8XY_3 = new Double_t[16];
                BuildTArb8Array(CryostatMiddleArb8XY_3, 3,4,2,3, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatMiddleArb8XY_4 = new Double_t[16];
                BuildTArb8Array(CryostatMiddleArb8XY_4, 4,5,2,3, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatMiddleArb8XY_5 = new Double_t[16];
                BuildTArb8Array(CryostatMiddleArb8XY_5, 5,6,2,3, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatMiddleArb8XY_6 = new Double_t[16];
                BuildTArb8Array(CryostatMiddleArb8XY_6, 6,7,2,3, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatMiddleArb8XY_7 = new Double_t[16];
                BuildTArb8Array(CryostatMiddleArb8XY_7, 7,8,2,3, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatMiddleArb8XY_8 = new Double_t[16];
                BuildTArb8Array(CryostatMiddleArb8XY_8, 8,9,2,3, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatMiddleArb8XY_9 = new Double_t[16];
                BuildTArb8Array(CryostatMiddleArb8XY_9, 9,10,2,3, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatMiddleArb8XY_10 = new Double_t[16];
                BuildTArb8Array(CryostatMiddleArb8XY_10, 10,1,2,3, NumberOfDifferentZValues-1, NumberOfXYInput);

                

        CryostatMiddleArb8_1 = new TGeoArb8 ("CryostatDoubleMiddleArb8_1",CryostatMiddleArb8Heigt/2,CryostatMiddleArb8XY_1);
        CryostatMiddleArb8_2 = new TGeoArb8 ("CryostatDoubleMiddleArb8_2",CryostatMiddleArb8Heigt/2,CryostatMiddleArb8XY_2);
        CryostatMiddleArb8_3 = new TGeoArb8 ("CryostatDoubleMiddleArb8_3",CryostatMiddleArb8Heigt/2,CryostatMiddleArb8XY_3);
        CryostatMiddleArb8_4 = new TGeoArb8 ("CryostatDoubleMiddleArb8_4",CryostatMiddleArb8Heigt/2,CryostatMiddleArb8XY_4);
        CryostatMiddleArb8_5 = new TGeoArb8 ("CryostatDoubleMiddleArb8_5",CryostatMiddleArb8Heigt/2,CryostatMiddleArb8XY_5);
        CryostatMiddleArb8_6 = new TGeoArb8 ("CryostatDoubleMiddleArb8_6",CryostatMiddleArb8Heigt/2,CryostatMiddleArb8XY_6);
        CryostatMiddleArb8_7 = new TGeoArb8 ("CryostatDoubleMiddleArb8_7",CryostatMiddleArb8Heigt/2,CryostatMiddleArb8XY_7);
        CryostatMiddleArb8_8 = new TGeoArb8 ("CryostatDoubleMiddleArb8_8",CryostatMiddleArb8Heigt/2,CryostatMiddleArb8XY_8);
        CryostatMiddleArb8_9 = new TGeoArb8 ("CryostatDoubleMiddleArb8_9",CryostatMiddleArb8Heigt/2,CryostatMiddleArb8XY_9);
        CryostatMiddleArb8_10 = new TGeoArb8 ("CryostatDoubleMiddleArb8_10",CryostatMiddleArb8Heigt/2,CryostatMiddleArb8XY_10);

        CryostatMiddleTranslation = new TGeoTranslation("CryostatDoubleMiddleTranslation",0,0,CryostatFrontThickness +CryostatLowerArb8Heigt+CryostatMiddleArb8Heigt/2);
        CryostatMiddleTranslation->RegisterYourself();
        
        for (Int_t i = 1; i < NumberOfCryostatEdges+1; i++)
        {
                char buffer[100];
                sprintf(buffer, "CryostatDoubleMiddleArb8_%d : CryostatDoubleMiddleTranslation",i);
                CompositeShapeExpression += buffer;
                //if (i< NumberOfCryostatEdges)         //uncomment if last part of composite shape
                {
                        CompositeShapeExpression += "+";
                }
        }


        //the upper part of the cryostat

        CryostatUpperArb8XY_1 = new Double_t[16];
                BuildTArb8Array(CryostatUpperArb8XY_1, 1,2,3,4, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatUpperArb8XY_2 = new Double_t[16];
                BuildTArb8Array(CryostatUpperArb8XY_2, 2,3,3,4, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatUpperArb8XY_3 = new Double_t[16];
                BuildTArb8Array(CryostatUpperArb8XY_3, 3,4,3,4, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatUpperArb8XY_4 = new Double_t[16];
                BuildTArb8Array(CryostatUpperArb8XY_4, 4,5,3,4, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatUpperArb8XY_5 = new Double_t[16];
                BuildTArb8Array(CryostatUpperArb8XY_5, 5,6,3,4, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatUpperArb8XY_6 = new Double_t[16];
                BuildTArb8Array(CryostatUpperArb8XY_6, 6,7,3,4, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatUpperArb8XY_7 = new Double_t[16];
                BuildTArb8Array(CryostatUpperArb8XY_7, 7,8,3,4, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatUpperArb8XY_8 = new Double_t[16];
                BuildTArb8Array(CryostatUpperArb8XY_8, 8,9,3,4, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatUpperArb8XY_9 = new Double_t[16];
                BuildTArb8Array(CryostatUpperArb8XY_9, 9,10,3,4, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatUpperArb8XY_10 = new Double_t[16];
                BuildTArb8Array(CryostatUpperArb8XY_10, 10,1,3,4, NumberOfDifferentZValues-1, NumberOfXYInput);

        CryostatUpperArb8_1 = new TGeoArb8 ("CryostatDoubleUpperArb8_1",CryostatUpperArb8Heigt/2,CryostatUpperArb8XY_1);
        CryostatUpperArb8_2 = new TGeoArb8 ("CryostatDoubleUpperArb8_2",CryostatUpperArb8Heigt/2,CryostatUpperArb8XY_2);
        CryostatUpperArb8_3 = new TGeoArb8 ("CryostatDoubleUpperArb8_3",CryostatUpperArb8Heigt/2,CryostatUpperArb8XY_3);
        CryostatUpperArb8_4 = new TGeoArb8 ("CryostatDoubleUpperArb8_4",CryostatUpperArb8Heigt/2,CryostatUpperArb8XY_4);
        CryostatUpperArb8_5 = new TGeoArb8 ("CryostatDoubleUpperArb8_5",CryostatUpperArb8Heigt/2,CryostatUpperArb8XY_5);
        CryostatUpperArb8_6 = new TGeoArb8 ("CryostatDoubleUpperArb8_6",CryostatUpperArb8Heigt/2,CryostatUpperArb8XY_6);
        CryostatUpperArb8_7 = new TGeoArb8 ("CryostatDoubleUpperArb8_7",CryostatUpperArb8Heigt/2,CryostatUpperArb8XY_7);
        CryostatUpperArb8_8 = new TGeoArb8 ("CryostatDoubleUpperArb8_8",CryostatUpperArb8Heigt/2,CryostatUpperArb8XY_8);
        CryostatUpperArb8_9 = new TGeoArb8 ("CryostatDoubleUpperArb8_9",CryostatUpperArb8Heigt/2,CryostatUpperArb8XY_9);
        CryostatUpperArb8_10 = new TGeoArb8 ("CryostatDoubleUpperArb8_10",CryostatUpperArb8Heigt/2,CryostatUpperArb8XY_10);

        CryostatUpperTranslation = new TGeoTranslation("CryostatDoubleUpperTranslation",0,0,CryostatFrontThickness +CryostatLowerArb8Heigt+CryostatMiddleArb8Heigt+CryostatUpperArb8Heigt/2);
        CryostatUpperTranslation->RegisterYourself();
        
        for (Int_t i = 1; i < NumberOfCryostatEdges+1; i++)
        {
                char buffer[100];
                sprintf(buffer, "CryostatDoubleUpperArb8_%d : CryostatDoubleUpperTranslation",i);
                CompositeShapeExpression += buffer;

                if (i<NumberOfCryostatEdges)            
                {
                        CompositeShapeExpression += "+";
                }
                if (i==NumberOfCryostatEdges && MakeTop  )
                {
                        CompositeShapeExpression += "+";
                }
        }

//the top of the cryostat

        CryostatTopArb8XY_1 = new Double_t[16];
                BuildTArb8TopArray(CryostatTopArb8XY_1, 1,2,3,10, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatTopArb8XY_2 = new Double_t[16];
                BuildTArb8TopArray(CryostatTopArb8XY_2, 3,4,5,6, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatTopArb8XY_3 = new Double_t[16];
                BuildTArb8TopArray(CryostatTopArb8XY_3, 3,6,7,8, NumberOfDifferentZValues-1, NumberOfXYInput);
        CryostatTopArb8XY_4 = new Double_t[16];
                BuildTArb8TopArray(CryostatTopArb8XY_4, 3,8,9,10, NumberOfDifferentZValues-1, NumberOfXYInput);

        CryostatTopArb8_1 = new TGeoArb8 ("CryostatDoubleTopArb8_1",CryostatTopThickness/2,CryostatTopArb8XY_1);
        CryostatTopArb8_2 = new TGeoArb8 ("CryostatDoubleTopArb8_2",CryostatTopThickness/2,CryostatTopArb8XY_2);
        CryostatTopArb8_3 = new TGeoArb8 ("CryostatDoubleTopArb8_3",CryostatTopThickness/2,CryostatTopArb8XY_3);
        CryostatTopArb8_4 = new TGeoArb8 ("CryostatDoubleTopArb8_4",CryostatTopThickness/2,CryostatTopArb8XY_4);

        CryostatTopTranslation = new TGeoTranslation("CryostatDoubleTopTranslation",0,0,CryostatFrontThickness +CryostatLowerArb8Heigt+CryostatMiddleArb8Heigt+CryostatUpperArb8Heigt+CryostatTopThickness/2);
        CryostatTopTranslation->RegisterYourself();


        if (MakeTop)
        {
                for (Int_t i = 1; i < NumberOfTopArb8+1; i++)
                {
                        char buffer[100];
                        sprintf(buffer, "CryostatDoubleTopArb8_%d : CryostatDoubleTopTranslation",i);
                        CompositeShapeExpression += buffer;
                        if (i<NumberOfTopArb8)          //uncomment if last part of composite shape
                        {
                                CompositeShapeExpression += "+";
                        }
                }
        }
        
        //cout << CompositeShapeExpression << endl;
        

        //combine the cryostat
        CryostatShape = new TGeoCompositeShape("CryostatShape",CompositeShapeExpression.data() ); 
        Cryostat = new TGeoVolume("Cryostat",CryostatShape,Al);
        Cryostat->SetLineColor(kGreen);
        Cryostat->SetFillColor(kGreen);
                        //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);

};

// ---------------  Prints the nodes of a cluster -----------------
void PndGeoHypGeDoubleCluster::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();
        }
}