GeoStripVolume.cxx

Go to the documentation of this file.

// $Id: GeoStripVolume.cxx,v 1.12 2008/01/09 06:38:11 schubert Exp $
//
// GeoStripVolume
//
// GeoStripVolume is a strip volume
// Basic shape is a box volume of polystyrene medium
// placed as a node inside a coextrusion box shell
//
// Author:  S. Kasahara 07/04

#include <iostream>
using namespace std;
#include <TGeoCompositeShape.h> // used for split strips
#include <TGeoMatrix.h>
#include <TGeoManager.h>

#include "GeoGeometry/Geo.h"
#include "GeoGeometry/GeoStripVolume.h"
#include "GeoGeometry/GeoGeometry.h"
#include "GeoGeometry/GeoMediumMap.h"
#include "GeoGeometry/GeoMedium.h"

ClassImp(GeoStripVolume)

//_____________________________________________________________________________
GeoStripVolume::GeoStripVolume(GeoGeometry* geo, const char* name,
                               Float_t totallen, Float_t leneastpart,
                               Float_t lenwestpart, Float_t wlsleneast,
                               Float_t wlslenwest, Float_t wlslenbypass,
                               Bool_t isvetoshield) 
  : GeoVolume(geo,name,0), 
    fTotalLen(totallen),fLenEastPart(leneastpart),fLenWestPart(lenwestpart),
    fWlsLenEast(wlsleneast), fWlsLenWest(wlslenwest),
    fWlsLenBypass(wlslenbypass),fWidth(0),fThick(0),fMirrorMask(0),
    fIsVetoShield(isvetoshield) {

  UpdateGlobalManager();

  const GeoMediumMap& medMap = fGeoGeometry -> GetMediumMap();
  
  if ( fIsVetoShield ) SetMedium(medMap.GetMedium(Geo::kVSStpTiO2));
  else SetMedium(medMap.GetMedium(Geo::kStpTiO2));
  
  Double_t scale = Geo::GetScale(fGeoGeometry->GetAppType());
  fWidth = scale*Geo::kStripWidth;
  fThick = scale*Geo::kStripThick;
   
  if ( fGeoGeometry->GetVldRange().GetDetectorMask() == Detector::kNear )
    fMirrorMask = StripEnd::kEast;
  
  // Build strip shape according to full pstyrene + coextrusion dimensions
  TGeoShape* shape = BuildStripShape();

  SetShape(shape);

  // Because shape was unknown at time of TGeoVolume construction, volume
  // was added to run-time list of volumes. Remove and add to normal list
  gGeoManager -> GetListOfGVolumes() -> Remove(this);
  Int_t index = gGeoManager->GetListOfVolumes() -> GetEntriesFast();
  gGeoManager -> GetListOfVolumes() -> AddAtAndExpand(this,index);

  // Activate pstyrene insert visiblity only by default
  this -> SetLineColor(kYellow);
  this -> SetVisibility(kFALSE); 
  this -> VisibleDaughters(kTRUE);
  
  // Build pstyrene only volume and insert this as node in full strip volume
  this -> AddPstyreneNodes();
}

//_____________________________________________________________________________
void GeoStripVolume::Print(Option_t* /* option */) const {
   // print to cout

  UpdateGlobalManager();

  cout << "GeoStripVolume::Print\n"
       << GetName() << endl;
  cout << "Length, Total " << fTotalLen << ", East " << fLenEastPart
       << ", West " << fLenWestPart << endl;
  cout << "Width " << fWidth << " Thick " << fThick << endl;
  cout << "Wls Length, East " << fWlsLenEast << ", West " << fWlsLenWest
       << ", Bypass " << fWlsLenBypass << endl;
  
}

//_____________________________________________________________________________
TGeoShape* GeoStripVolume::BuildStripShape() const {
  // Private method used to build shape for this strip.  
  // If strip is split, shape will be a composite shape consisting of 
  // a TGeoBBox for the total length - TGeoBBox for the hole region.  If 
  // strip is whole, shape will be a TGeoBBox.

  UpdateGlobalManager();
  
  TGeoShape* shape = 0;

  if ( IsSplit() ) {
    // strip is split
    std::string boxname = std::string(GetName())+"_box";
    new TGeoBBox(boxname.c_str(),0.5*fTotalLen,0.5*fWidth,0.5*fThick);
    Float_t lengap = fTotalLen-fLenEastPart-fLenWestPart;
    std::string gapname = std::string(GetName())+"_gap";
    new TGeoBBox(gapname.c_str(),0.5*lengap,0.5*fWidth,0.5*fThick);
    Float_t xgap = 0.5*(fLenEastPart-fLenWestPart);
    TGeoTranslation* trGap = new TGeoTranslation(gapname.c_str(),xgap,0,0);
    trGap -> RegisterYourself();

    std::string compdef = boxname + "-" + gapname + ":" + gapname;
    shape = new TGeoCompositeShape(GetName(),compdef.c_str());
  }
  else {
    shape = new TGeoBBox(GetName(),0.5*fTotalLen,0.5*fWidth,0.5*fThick);
  }     

  return shape;

} 

//_____________________________________________________________________________
void GeoStripVolume::AddPstyreneNodes() {
  // Private method used to insert pstyrene nodes for this strip.

  UpdateGlobalManager();

  const GeoMediumMap& medMap = fGeoGeometry -> GetMediumMap();

  Geo::EDetComponent detcomp = Geo::kStpScint;
  if ( fIsVetoShield ) detcomp = Geo::kVSStpScint;

  Float_t scale = Geo::GetScale(fGeoGeometry->GetAppType());
  Float_t coexthick = Geo::kCoexThick*scale;
  std::string stripname = std::string(GetName());
  Float_t pstythick = fThick - 2.*coexthick;
  Float_t pstywidth = fWidth - 2.*coexthick;

  if ( IsSplit() ) {
    // strip is split
    std::string pstyeastname = stripname+"_pstyeast";
    std::string pstywestname = stripname+"_pstywest";

    TGeoVolume* eastvol = gGeoManager -> MakeBox(pstyeastname.c_str(),
                  medMap.GetMedium(detcomp),
                  0.5*fLenEastPart-coexthick,0.5*pstywidth,0.5*pstythick);
    TGeoVolume* westvol = gGeoManager -> MakeBox(pstywestname.c_str(),
                  medMap.GetMedium(detcomp),
                  0.5*fLenWestPart-coexthick,0.5*pstywidth,0.5*pstythick);
    eastvol -> SetLineColor(kCyan); 
    eastvol -> SetVisibility(kTRUE);
    westvol -> SetLineColor(kCyan); 
    westvol -> SetVisibility(kTRUE);
    Float_t xeast = 0.5*(-fTotalLen+fLenEastPart); // east(-)
    Float_t xwest = 0.5*(fTotalLen-fLenWestPart); // west(+)
    this->AddNode(eastvol,1,
                  new TGeoTranslation(pstyeastname.c_str(),xeast,0,0));
    this->AddNode(westvol,1,
                  new TGeoTranslation(pstywestname.c_str(),xwest,0,0));
  }
  else {
    std::string pstyname = stripname+"_psty";
    TGeoVolume* vol = gGeoManager -> MakeBox(pstyname.c_str(),
                  medMap.GetMedium(detcomp),
                  0.5*fTotalLen-coexthick,0.5*pstywidth,0.5*pstythick);
    vol -> SetLineColor(kCyan);
    vol -> SetVisibility(kTRUE);
    this -> AddNode(vol,1,gGeoIdentity);
  }     

  return;

} 

//_____________________________________________________________________________
Float_t GeoStripVolume::GetWlsLen(StripEnd::StripEnd_t end) const {
   // Return length of wls (east,west, or total)

  UpdateGlobalManager();

  switch(end) {

  case (StripEnd::kEast): return fWlsLenEast;

  case (StripEnd::kWest): return fWlsLenWest;

  case (StripEnd::kWhole): {      
    if ( !fWlsLenBypass ) {
      return fWlsLenEast + fWlsLenWest + fTotalLen;
    }
    else {
      return fWlsLenEast + fWlsLenWest + fWlsLenBypass 
           + fLenEastPart + fLenWestPart;
    }
  }    
  default: return 0;

  } // end of switch
}

//_____________________________________________________________________________
Float_t GeoStripVolume::GetLenPart(StripEnd::StripEnd_t end) const {
   // Return length of strip (east,west, or total)

  UpdateGlobalManager();

  switch(end) {

  case (StripEnd::kEast): return fLenEastPart;

  case (StripEnd::kWest): return fLenWestPart;

  case (StripEnd::kWhole): return fTotalLen;

  default: return 0;

  } // end of switch
  
}

---