MaGeGeometryGermaniumCrystalHole.cc

Go to the documentation of this file.

//---------------------------------------------------------------------------//
//bb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nu//
//                                                                           //
//                         MAJORANA Simulation                               //
//                                                                           //
//      This code implementation is the intellectual property of the         //
//      MAJORANA Collaboration. It is based on Geant4, an intellectual       //
//      property of the RD44 GEANT4 collaboration.                           //
//                                                                           //
//                        *********************                              //
//                                                                           //
//    Neither the authors of this software system, nor their employing       //
//    institutes, nor the agencies providing financial support for this      //
//    work  make  any representation or  warranty, express or implied,       //
//    regarding this software system or assume any liability for its use.    //
//    By copying, distributing or modifying the Program (or any work based   //
//    on on the Program) you indicate your acceptance of this statement,     //
//    and all its terms.                                                     //
//                                                                           //
//bb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nu//
//---------------------------------------------------------------------------//
//                                                          
// $Id: MaGeGeometryGermaniumCrystalHole.cc,v 1.1 2004/11/26 09:46:30 pandola Exp $ 
//      
// CLASS IMPLEMENTATION:  MaGeGeometryGermaniumCrystal.cc
//
//---------------------------------------------------------------------------//
// 
//---------------------------------------------------------------------------//
//---------------------------------------------------------------------------//
//
//G4 class headers
#include "globals.hh"
#include "G4LogicalVolume.hh"
#include "G4PVPlacement.hh"
#include "G4VPhysicalVolume.hh"
#include "G4Material.hh"
#include "G4Tubs.hh"
#include "G4VisAttributes.hh"
#include "G4Colour.hh"
#include "G4SDManager.hh"
#include "G4SubtractionSolid.hh"

//---------------------------------------------------------------------------//

#include "gerdageometry/MaGeGeometrySD.hh"
// #include "gerdageometry/MaGeGeometryGermaniumCrystalHoleMessenger.hh"
#include "gerdageometry/MaGeGeometryGermaniumCrystalHole.hh"                  
//Present MJ Class Headers 

//---------------------------------------------------------------------------//

MaGeGeometryGermaniumCrystalHole::MaGeGeometryGermaniumCrystalHole() :
    //GeCrystalInnerRadius(1.0*cm),
GeCrystalInnerRadius(0.0*cm),
GeCrystalOuterRadius(4.0*cm),
GeCrystalHeight(8.0*cm),
GeCrystalDeadlayer(0.01*cm),
GeCrystalHoleDepth(2.5*cm),
GeCrystalHoleOuterRadius(0.5*cm),
theGermaniumCrystalActiveLogical(0),
theGermaniumCrystalDeadlayerLogical(0)
{
  // coded in for the momentum
  //GeCrystalInnerRadius=1.0*cm;
  //GeCrystalOuterRadius=4.0*cm;
  //GeCrystalHeight     =8.0*cm;
  //GeCrystalDeadlayer  =0.01*cm;   
  //  fMessenger = new MaGeGeometryGermaniumCrystalMessenger(this);
//  ConstructCrystal();
}


MaGeGeometryGermaniumCrystalHole::~MaGeGeometryGermaniumCrystalHole()
{;}

void MaGeGeometryGermaniumCrystalHole::ConstructCrystal()
{
  //
  // ---- define material (will put into database later) ---- //
  // 

  G4Material* ge = G4Material::GetMaterial("Ge");
  //
  // ---- define volumes ---- //
  //
    G4double startAngleOfTheTube = 0.*deg;
    G4double spanningAngleOfTheTube = 360.*deg;
    G4double heightofTheTube = 0.5*GeCrystalHeight;
    G4double holeInnerRadius = 0.0*cm;
    G4ThreeVector  vector_hole(0,0,heightofTheTube-GeCrystalHoleDepth+0.5*mm);

    // the additional 0.5mm is to avoid the two volumes (crystal and hole) to 
    // share the same surface, it can create problems with visualization 

    // --- active part of crystal --- //
    G4Tubs*  GeCrystal_ActiveTub = new G4Tubs("gecrystal_activetub",
                                       GeCrystalInnerRadius,
                                       GeCrystalOuterRadius-GeCrystalDeadlayer,
                                       heightofTheTube,
                                       startAngleOfTheTube,
                                       spanningAngleOfTheTube);

    G4Tubs*  GeCrystal_Hole = new G4Tubs("gecrystal_hole",
                                       holeInnerRadius,
                                       GeCrystalHoleOuterRadius,
                                       GeCrystalHoleDepth,
                                       startAngleOfTheTube,
                                       spanningAngleOfTheTube);

    G4SubtractionSolid *GeCrystal_solid = new G4SubtractionSolid("GeCrystal_solid",GeCrystal_ActiveTub,GeCrystal_Hole,0,vector_hole);
 
   theGermaniumCrystalActiveLogical = new G4LogicalVolume(
                                       GeCrystal_solid,ge,
                                       "gecrystal_activelogical",0,0,0);

  // --- inactive part of the crystal --- //
    G4Tubs* GeCrystal_DeadlayerTub = new G4Tubs("gecrystal_deadlayertub",
                                       GeCrystalOuterRadius-GeCrystalDeadlayer,
                                       GeCrystalOuterRadius,
                                       heightofTheTube,
                                       startAngleOfTheTube,
                                       spanningAngleOfTheTube);
    theGermaniumCrystalDeadlayerLogical = new G4LogicalVolume(
                                       GeCrystal_DeadlayerTub,ge,
                                       "gecrystal_deadlayerlogical",0,0,0);

  //--- visualization attributes ---//
    G4VisAttributes* GeCrystalVisAtt=new G4VisAttributes(G4Colour(1.0,0.,0.));
    GeCrystalVisAtt->SetVisibility(true);
    GeCrystalVisAtt->SetForceSolid(true);

    theGermaniumCrystalActiveLogical->SetVisAttributes(GeCrystalVisAtt);
    theGermaniumCrystalDeadlayerLogical->SetVisAttributes(GeCrystalVisAtt);

  //--- sensitive detector definition ---//
  G4SDManager* SDman = G4SDManager::GetSDMpointer();

  G4String geCrystalSDname = "/mydet/gerda/gecrystal";
  MaGeGeometrySD * geCrystalSD = new MaGeGeometrySD(geCrystalSDname);
  SDman->AddNewDetector(geCrystalSD);
  theGermaniumCrystalActiveLogical->SetSensitiveDetector(geCrystalSD);

};

//
                         

---