strawtubes Class Reference

#include <strawtubes.h>

Inheritance diagram for strawtubes:

Collaboration diagram for strawtubes:

Public Member Functions
	strawtubes (const char *Name, Bool_t Active)
	strawtubes (std::string medium)
	strawtubes ()
Bool_t	ProcessHits (FairVolume *v=0) override
void	SetzPositions (Double_t z1, Double_t z2, Double_t z3, Double_t z4)
void	SetApertureArea (Double_t width, Double_t height)
void	SetStrawDiameter (Double_t outer_straw_diameter, Double_t wall_thickness)
void	SetStrawPitch (Double_t straw_pitch, Double_t layer_offset)
void	SetDeltazLayer (Double_t delta_z_layer)
void	SetStereoAngle (Double_t stereo_angle)
void	SetWireThickness (Double_t wire_thickness)
void	SetFrameMaterial (TString frame_material)
void	SetDeltazView (Double_t delta_z_view)
void	SetStationEnvelope (Double_t x, Double_t y, Double_t z)
void	SetStrawResolution (Double_t a, Double_t b)
Double_t	StrawVdrift ()
Double_t	StrawSigmaSpatial ()
void	ConstructGeometry () override
Public Member Functions inherited from SHiP::Detector< strawtubesPoint >
	Detector ()=default
	Detector (const char *Name, Bool_t Active, Int_t detID)
	Detector (const char *Name, Bool_t Active)
	~Detector () override=default
strawtubesPoint *	AddHit (Args &&... args)
void	ConstructGeometry () override=0
void	Initialize () override
void	Reset () override
void	EndOfEvent () override
void	Register () override
TClonesArray *	GetCollection (Int_t iColl) const override
void	UpdatePointTrackIndices (const std::map< Int_t, Int_t > &indexMap) override
	Update track indices in point collection after track filtering.
void	SetSpecialPhysicsCuts () override
void	FinishPrimary () override
void	FinishRun () override
void	BeginPrimary () override
void	PostTrack () override
void	PreTrack () override
void	BeginEvent () override
void	CopyClones (TClonesArray *cl1, TClonesArray *cl2, Int_t offset) override
Public Member Functions inherited from ISTLPointContainer
virtual void	UpdatePointTrackIndices (const std::map< Int_t, Int_t > &indexMap)=0
	Update track indices in point collection after track filtering.
virtual	~ISTLPointContainer ()=default

Static Public Member Functions
static std::array< Int_t, 4 >	StrawDecode (Int_t detID)
static void	StrawEndPoints (Int_t detID, TVector3 &top, TVector3 &bot)

Private Member Functions
	strawtubes (const strawtubes &)=delete
	vacuum box medium
strawtubes &	operator= (const strawtubes &)=delete

Private Attributes
Double_t	f_T1_z
Double_t	f_T2_z
	z-position of tracking station 1
Double_t	f_T3_z
	z-position of tracking station 2
Double_t	f_T4_z
	z-position of tracking station 3
Double_t	f_aperture_width
	z-position of tracking station 4
Double_t	f_aperture_height
	Aperture width (x)
Double_t	f_inner_straw_diameter
	Aperture height (y)
Double_t	f_outer_straw_diameter
	Inner Straw diameter.
Double_t	f_straw_pitch
	Outer Straw diameter.
Double_t	f_offset_layer
	Distance (y) between straws in a layer.
Double_t	f_delta_z_layer
	Offset (y) of straws between layers.
Double_t	f_view_angle
	Distance (z) between layers.
Double_t	f_wire_thickness
	Stereo view angle.
TString	f_frame_material
	Sense wire thickness.
Double_t	f_delta_z_view
	Structure frame material.
Double_t	f_station_width
	Distance (z) between stereo views.
Double_t	f_station_height
	Station envelope width (x)
Double_t	f_station_length
	Station envelope height (y)
Double_t	v_drift
	Station envelope length (z)
Double_t	sigma_spatial
	drift velocity
std::string	fMedium
	spatial resolution

Additional Inherited Members
Protected Attributes inherited from SHiP::Detector< strawtubesPoint >
Int_t	fEventID
Int_t	fTrackID
	event index
Int_t	fVolumeID
	track index
TLorentzVector	fPos
	volume id
TLorentzVector	fMom
	position at entrance
Double_t	fTime
	momentum at entrance
Double_t	fLength
	time
Double_t	fELoss
	length
std::vector< strawtubesPoint > *	fDetPoints
	energy loss
TGeoVolume *	fDetector

Detailed Description

Definition at line 19 of file strawtubes.h.

Constructor & Destructor Documentation

strawtubes() [1/4]

strawtubes::strawtubes	(	const char *	Name,
		Bool_t	Active
	)

Name : Detector Name Active: kTRUE for active detectors (ProcessHits() will be called) kFALSE for inactive detectors

Definition at line 50 of file strawtubes.cxx.

    : Detector(name, active, kStraw) {}

strawtubes() [2/4]

strawtubes::strawtubes ( std::string  medium )

explicit

Definition at line 47 of file strawtubes.cxx.

    : Detector("strawtubes", kTRUE, kStraw), fMedium(medium) {}

strawtubes() [3/4]

strawtubes::strawtubes

(

)

default constructor

Definition at line 44 of file strawtubes.cxx.

    : Detector("strawtubes", kTRUE, kStraw), fMedium("air") {}

strawtubes() [4/4]

strawtubes::strawtubes ( const strawtubes &  )

privatedelete

vacuum box medium

container for data points

Member Function Documentation

ConstructGeometry()

void strawtubes::ConstructGeometry ( )

overridevirtual

Create the detector geometry

If you are using the standard ASCII input for the geometry just copy this and use it for your detector, otherwise you can implement here you own way of constructing the geometry.

Implements SHiP::Detector< strawtubesPoint >.

Definition at line 184 of file strawtubes.cxx.

                                   {
  TGeoVolume* top = gGeoManager->GetTopVolume();
  ShipGeo::InitMedium("mylar");
  TGeoMedium* mylar = gGeoManager->GetMedium("mylar");
  ShipGeo::InitMedium("STTmix8020_1bar");
  TGeoMedium* sttmix8020_1bar = gGeoManager->GetMedium("STTmix8020_1bar");
  ShipGeo::InitMedium("tungsten");
  TGeoMedium* tungsten = gGeoManager->GetMedium("tungsten");
  ShipGeo::InitMedium(f_frame_material);
  TGeoMedium* FrameMatPtr = gGeoManager->GetMedium(f_frame_material);
  ShipGeo::InitMedium(fMedium.c_str());
  TGeoMedium* med = gGeoManager->GetMedium(fMedium.c_str());
 
  gGeoManager->SetVisLevel(4);
  gGeoManager->SetTopVisible();
 
  // Epsilon to avoid overlapping volumes
  Double_t eps = 0.0001;
  // Straw (half) length
  Double_t straw_length = f_aperture_width + 2. * eps;
  // Width of frame: standard HEA 500 I-beam width
  Double_t frame_width = 49.;
  // Offset due to floor space limitation
  Double_t floor_offset = 14.;
 
  Double_t rmin, rmax, T_station_z;
 
  // Arguments of boxes are half-lengths
  [[maybe_unused]] TGeoBBox* detbox1 = new TGeoBBox(
      "detbox1", f_aperture_width + frame_width,
      f_aperture_height + frame_width - floor_offset / 2., f_station_length);
  [[maybe_unused]] TGeoBBox* detbox2 = new TGeoBBox(
      "detbox2", straw_length + eps,
      f_aperture_height +
          TMath::Tan(f_view_angle * TMath::Pi() / 180.0) * straw_length * 2 +
          f_offset_layer / TMath::Cos(f_view_angle * TMath::Pi() / 180.0) + eps,
      f_station_length + eps);
  TGeoTranslation* move_up =
      new TGeoTranslation("move_up", 0, floor_offset / 2., 0);
  move_up->RegisterYourself();
 
  // Composite shape to create frame
  TGeoCompositeShape* detcomp1 =
      new TGeoCompositeShape("detcomp1", "(detbox1:move_up)-detbox2");
 
  // Volume: straw
  rmin = f_inner_straw_diameter / 2.;
  rmax = f_outer_straw_diameter / 2.;
  // Third argument is half-length of tube
  TGeoTube* straw_tube = new TGeoTube("straw", rmin, rmax, straw_length);
  TGeoVolume* straw = new TGeoVolume("straw", straw_tube, mylar);
  straw->SetLineColor(4);
  straw->SetVisibility(kTRUE);
 
  // Volume: gas
  rmin = f_wire_thickness / 2. + eps;
  rmax = f_inner_straw_diameter / 2. - eps;
  TGeoTube* gas_tube = new TGeoTube("gas", rmin, rmax, straw_length - 2. * eps);
  TGeoVolume* gas = new TGeoVolume("gas", gas_tube, sttmix8020_1bar);
  gas->SetLineColor(5);
  // Only the gas is sensitive
  AddSensitiveVolume(gas);
 
  // Volume: wire
  rmin = 0.;
  rmax = f_wire_thickness / 2.;
  TGeoTube* wire_tube =
      new TGeoTube("wire", rmin, rmax, straw_length - 4. * eps);
  TGeoVolume* wire = new TGeoVolume("wire", wire_tube, tungsten);
  wire->SetLineColor(6);
 
  // Tracking stations
  // statnb = station number; vnb = view number; lnb = layer number; snb = straw
  // number
 
  // Station box to contain all components
  TGeoBBox* statbox =
      new TGeoBBox("statbox", f_station_width,
                   f_station_height - floor_offset / 2., f_station_length);
 
  f_frame_material.ToLower();
 
  for (Int_t statnb = 1; statnb < 5; statnb++) {
    // Tracking station loop
    TString nmstation = "Tr";
    std::stringstream ss;
    ss << statnb;
    nmstation = nmstation + ss.str();
    switch (statnb) {
      case 1:
        T_station_z = f_T1_z;
        break;
      case 2:
        T_station_z = f_T2_z;
        break;
      case 3:
        T_station_z = f_T3_z;
        break;
      case 4:
        T_station_z = f_T4_z;
        break;
      default:
        T_station_z = f_T1_z;
    }
 
    TGeoVolume* vol = new TGeoVolume(nmstation, statbox, med);
    // z-translate the station to its (absolute) position
    top->AddNode(vol, statnb,
                 new TGeoTranslation(0, floor_offset / 2., T_station_z));
 
    TGeoVolume* statframe =
        new TGeoVolume(nmstation + "_frame", detcomp1, FrameMatPtr);
    vol->AddNode(statframe, statnb * 1e6,
                 new TGeoTranslation(0, -floor_offset / 2., 0));
    statframe->SetLineColor(kRed);
 
    for (Int_t vnb = 0; vnb < 4; vnb++) {
      // View loop
      const Double_t angle = [&] {
        switch (vnb) {
          case 1:
            return f_view_angle;
          case 2:
            return -f_view_angle;
          default:
            return 0.;
        }
      }();
      const TString nmview = [&] {
        switch (vnb) {
          case 0:
            return nmstation + "_y1";
          case 1:
            return nmstation + "_u";
          case 2:
            return nmstation + "_v";
          case 3:
            return nmstation + "_y2";
          default:
            return nmstation + "_y1";
        }
      }();
 
      // Adjustments in the stereo views
      // stereo_growth: extension of stereo views beyond aperture
      // stereo_pitch: straw pitch in stereo views
      // offset_layer: layer offset in stereo views
      // straws_per_layer: number of straws in one layer with stereo extension
      // If angle == 0., all numbers return the case of non-stereo views.
      const Double_t stereo_growth =
          TMath::Tan(TMath::Abs(angle) * TMath::Pi() / 180.0) * straw_length;
      const Double_t stereo_pitch =
          f_straw_pitch / TMath::Cos(TMath::Abs(angle) * TMath::Pi() / 180.0);
      const Double_t offset_layer =
          f_offset_layer / TMath::Cos(TMath::Abs(angle) * TMath::Pi() / 180.0);
      const Int_t straws_per_layer =
          std::ceil(2 * (f_aperture_height + stereo_growth) / stereo_pitch);
 
      for (Int_t lnb = 0; lnb < 2; lnb++) {
        // Layer loop
        TString nmlayer = nmview + "_layer_";
        nmlayer += lnb;
        TGeoBBox* layer = new TGeoBBox(
            "layer box", straw_length + eps / 4,
            f_aperture_height + stereo_growth * 2 + offset_layer + eps / 4,
            f_outer_straw_diameter / 2. + eps / 4);
        TGeoVolume* layerbox = new TGeoVolume(nmlayer, layer, med);
 
        // The layer box sits in the viewframe.
        // Hence, z-translate the layer w.r.t. the view
        vol->AddNode(
            layerbox, statnb * 1e6 + vnb * 1e5 + lnb * 1e4,
            new TGeoTranslation(0, -floor_offset / 2.,
                                (vnb - 3. / 2.) * f_delta_z_view +
                                    (lnb - 1. / 2.) * f_delta_z_layer));
 
        TGeoRotation r6s;
        TGeoTranslation t6s;
        for (Int_t snb = 1; snb <= straws_per_layer; snb++) {
          // Straw loop
          // y-translate the straw to its position
          t6s.SetTranslation(0,
                             f_aperture_height + stereo_growth -
                                 (snb - 1. / 2.) * stereo_pitch +
                                 lnb * offset_layer,
                             0);
          // Rotate the straw with stereo angle
          r6s.SetAngles(90 + angle, 90, 0);
          TGeoCombiTrans c6s(t6s, r6s);
          TGeoHMatrix* h6s = new TGeoHMatrix(c6s);
          layerbox->AddNode(
              straw, statnb * 1e6 + vnb * 1e5 + lnb * 1e4 + 1e3 + snb, h6s);
          layerbox->AddNode(
              gas, statnb * 1e6 + vnb * 1e5 + lnb * 1e4 + 2e3 + snb, h6s);
          layerbox->AddNode(
              wire, statnb * 1e6 + vnb * 1e5 + lnb * 1e4 + 3e3 + snb, h6s);
          // End of straw loop
        }
 
        if (lnb == 1) {
          // Add one more straw at the bottom of the second layer to cover
          // aperture entirely
          t6s.SetTranslation(0,
                             f_aperture_height + stereo_growth -
                                 (straws_per_layer - 1. / 2.) * stereo_pitch -
                                 lnb * offset_layer,
                             0);
          r6s.SetAngles(90 + angle, 90, 0);
          TGeoCombiTrans c6s(t6s, r6s);
          TGeoHMatrix* h6s = new TGeoHMatrix(c6s);
          layerbox->AddNode(
              straw,
              statnb * 1e6 + vnb * 1e5 + lnb * 1e4 + 1e3 + straws_per_layer + 1,
              h6s);
          layerbox->AddNode(
              gas,
              statnb * 1e6 + vnb * 1e5 + lnb * 1e4 + 2e3 + straws_per_layer + 1,
              h6s);
          layerbox->AddNode(
              wire,
              statnb * 1e6 + vnb * 1e5 + lnb * 1e4 + 3e3 + straws_per_layer + 1,
              h6s);
        }
        // End of layer loop
      }
      // End of view loop
    }
    // End of tracking station loop
  }
}

operator=()

strawtubes & strawtubes::operator= ( const strawtubes &  )

privatedelete

ProcessHits()

Bool_t strawtubes::ProcessHits ( FairVolume *  v = 0 )

override

this method is called for each step during simulation (see FairMCApplication::Stepping())

This method is called from the MC stepping

Definition at line 53 of file strawtubes.cxx.

                                              {
  // Set parameters at entrance of volume. Reset ELoss.
  if (gMC->IsTrackEntering()) {
    fELoss = 0.;
    fTime = gMC->TrackTime() * 1.0e09;
    fLength = gMC->TrackLength();
    gMC->TrackPosition(fPos);
    gMC->TrackMomentum(fMom);
  }
  // Sum energy loss for all steps in the active volume
  fELoss += gMC->Edep();
 
  // Create strawtubesPoint at exit of active volume
  if (gMC->IsTrackExiting() || gMC->IsTrackStop() ||
      gMC->IsTrackDisappeared()) {
    if (fELoss == 0.) {
      return kFALSE;
    }
    TParticle* p = gMC->GetStack()->GetCurrentTrack();
    Int_t pdgCode = p->GetPdgCode();
    fTrackID = gMC->GetStack()->GetCurrentTrackNumber();
    fEventID = gMC->CurrentEvent();
    Int_t straw_uniqueId;
    gMC->CurrentVolID(straw_uniqueId);
    if (fVolumeID == straw_uniqueId) {
      // std::cout << pdgCode<< " same volume again ? "<< straw_uniqueId << "
      // exit:" << gMC->IsTrackExiting() << " stop:" << gMC->IsTrackStop() << "
      // disappeared:" << gMC->IsTrackDisappeared()<< std::endl;
      return kTRUE;
    }
    fVolumeID = straw_uniqueId;
    // # d = |pq . u x v|/|u x v|
    TVector3 bot, top;
    StrawEndPoints(straw_uniqueId, bot, top);
    TLorentzVector Pos;
    gMC->TrackPosition(Pos);
    Double_t xmean = (fPos.X() + Pos.X()) / 2.;
    Double_t ymean = (fPos.Y() + Pos.Y()) / 2.;
    Double_t zmean = (fPos.Z() + Pos.Z()) / 2.;
    TVector3 pq = TVector3(top.x() - xmean, top.y() - ymean, top.z() - zmean);
    TVector3 u =
        TVector3(bot.x() - top.x(), bot.y() - top.y(), bot.z() - top.z());
    TVector3 v =
        TVector3(fPos.X() - Pos.X(), fPos.Y() - Pos.Y(), fPos.Z() - Pos.Z());
    TVector3 uCrossv = u.Cross(v);
    Double_t dist2Wire = fabs(pq.Dot(uCrossv)) / (uCrossv.Mag() + 1E-8);
    Double_t deltaTrackLength = gMC->TrackLength() - fLength;
    AddHit(fEventID, fTrackID, straw_uniqueId, TVector3(xmean, ymean, zmean),
           TVector3(fMom.Px(), fMom.Py(), fMom.Pz()), fTime, deltaTrackLength,
           fELoss, pdgCode, dist2Wire);
    if (dist2Wire > f_inner_straw_diameter / 2) {
      std::cout << "addhit " << dist2Wire << " straw id " << straw_uniqueId
                << " pdgcode " << pdgCode << " dot prod " << pq.Dot(uCrossv)
                << std::endl;
      std::cout << " exit:" << gMC->IsTrackExiting()
                << " stop:" << gMC->IsTrackStop()
                << " disappeared:" << gMC->IsTrackDisappeared() << std::endl;
      std::cout << " entry:" << fPos.X() << " " << fPos.Y() << " " << fPos.Z()
                << std::endl;
      std::cout << " exit:" << Pos.X() << " " << Pos.Y() << " " << Pos.Z()
                << std::endl;
      std::cout << " mean:" << xmean << " " << ymean << " " << zmean
                << std::endl;
      std::cout << " bot:" << bot.x() << " " << bot.y() << " " << bot.z()
                << std::endl;
      std::cout << " top:" << top.x() << " " << top.y() << " " << top.z()
                << std::endl;
      pq.Print();
      u.Print();
      v.Print();
      uCrossv.Print();
    }
    // Increment number of strawtubes det points in TParticle
    ShipStack* stack = dynamic_cast<ShipStack*>(gMC->GetStack());
    stack->AddPoint(kStraw);
  }
  return kTRUE;
}

SetApertureArea()

void strawtubes::SetApertureArea	(	Double_t	width,
		Double_t	height
	)

Aperture width (x)

Aperture height (y)

Definition at line 141 of file strawtubes.cxx.

                                                                {
  f_aperture_width = width;    
  f_aperture_height = height;  
}

SetDeltazLayer()

void strawtubes::SetDeltazLayer

(

Double_t

delta_z_layer

)

Distance (z) between layers

Definition at line 158 of file strawtubes.cxx.

                                                      {
  f_delta_z_layer = delta_z_layer;  
}

SetDeltazView()

void strawtubes::SetDeltazView

(

Double_t

delta_z_view

)

Distance (z) between stereo views

Definition at line 170 of file strawtubes.cxx.

                                                    {
  f_delta_z_view = delta_z_view;  
}

SetFrameMaterial()

void strawtubes::SetFrameMaterial

(

TString

frame_material

)

Structure frame material

Definition at line 174 of file strawtubes.cxx.

                                                        {
  f_frame_material = frame_material;  
}

SetStationEnvelope()

void strawtubes::SetStationEnvelope	(	Double_t	x,
		Double_t	y,
		Double_t	z
	)

Station envelope width (x)

Station envelope height (y)

Station envelope length (z)

Definition at line 178 of file strawtubes.cxx.

                                                                      {
  f_station_width = x;   
  f_station_height = y;  
  f_station_length = z;  
}

SetStereoAngle()

void strawtubes::SetStereoAngle

(

Double_t

stereo_angle

)

Stereo view angle

Definition at line 162 of file strawtubes.cxx.

                                                     {
  f_view_angle = stereo_angle;  
}

SetStrawDiameter()

void strawtubes::SetStrawDiameter	(	Double_t	outer_straw_diameter,
		Double_t	wall_thickness
	)

Outer straw diameter

Inner straw diameter

Definition at line 146 of file strawtubes.cxx.

                                                           {
  f_outer_straw_diameter = outer_straw_diameter;  
  f_inner_straw_diameter =
      outer_straw_diameter - 2 * wall_thickness;  
}

SetStrawPitch()

void strawtubes::SetStrawPitch	(	Double_t	straw_pitch,
		Double_t	layer_offset
	)

Distance (y) between straws in a layer

Offset (y) of straws between layers

Definition at line 153 of file strawtubes.cxx.

                                                                          {
  f_straw_pitch = straw_pitch;    
  f_offset_layer = layer_offset;  
}

SetStrawResolution()

void strawtubes::SetStrawResolution ( Double_t  a, Double_t  b  )

inline

Definition at line 50 of file strawtubes.h.

                                                  {
    v_drift = a;
    sigma_spatial = b;
  }

SetWireThickness()

void strawtubes::SetWireThickness

(

Double_t

wire_thickness

)

Sense wire thickness

Definition at line 166 of file strawtubes.cxx.

                                                         {
  f_wire_thickness = wire_thickness;  
}

SetzPositions()

void strawtubes::SetzPositions	(	Double_t	z1,
		Double_t	z2,
		Double_t	z3,
		Double_t	z4
	)

z-position of tracking station 1

z-position of tracking station 2

z-position of tracking station 3

z-position of tracking station 4

Definition at line 133 of file strawtubes.cxx.

                                            {
  f_T1_z = z1;  
  f_T2_z = z2;  
  f_T3_z = z3;  
  f_T4_z = z4;  
}

StrawDecode()

std::array< Int_t, 4 > strawtubes::StrawDecode ( Int_t  detID )

static

Definition at line 421 of file strawtubes.cxx.

                                                      {
  Int_t statnb, vnb, lnb, snb;
  statnb = detID / 1e6;
  vnb = (detID - statnb * 1e6) / 1e5;
  lnb = (detID - statnb * 1e6 - vnb * 1e5) / 1e4;
  snb = detID - statnb * 1e6 - vnb * 1e5 - lnb * 1e4 - 2e3;
 
  if (statnb < 1 || statnb > 4 || vnb < 0 || vnb > 3 || lnb < 0 || lnb > 1 ||
      snb < 1 || snb > 317) {
    LOG(warning) << "Invalid strawtubes detID:";
    LOG(warning) << detID << " -> station: " << statnb << ", view: " << vnb
                 << ", layer: " << lnb << ", straw: " << snb;
    LOG(warning) << "strawtubes detID is 7-digit!";
    return {0, -1, -1, 0};
  } else {
    return {statnb, vnb, lnb, snb};
  }
}

StrawEndPoints()

void strawtubes::StrawEndPoints ( Int_t  detID, TVector3 &  top, TVector3 &  bot  )

static

Definition at line 443 of file strawtubes.cxx.

{
  const auto [statnb, vnb, lnb, snb] = StrawDecode(fDetectorID);
  TString stat = "Tr";
  stat += statnb;
  stat += "_";
  stat += statnb;
  const TString view = [&] {
    switch (vnb) {
      case 1:
        return TString("_u");
      case 2:
        return TString("_v");
      case 3:
        return TString("_y2");
      default:
        return TString("_y1");
    }
  }();
  TGeoNavigator* nav = gGeoManager->GetCurrentNavigator();
  TString prefix = "Tr";
  prefix += statnb;
  prefix += view;
  prefix += "_";
  TString layer = prefix + "layer_";
  layer += lnb;
  layer += "_";
  layer += statnb;
  layer += vnb;
  layer += lnb;
  layer += "0000";
  TString wire = "wire_";
  wire += fDetectorID + 1e3;
  TString path = "/";
  path += stat;
  path += "/";
  path += layer;
  path += "/";
  path += wire;
  Bool_t rc = nav->cd(path);
  if (!rc) {
    LOG(warning) << "strawtubes::StrawDecode, TGeoNavigator failed" << path;
    return;
  }
  TGeoNode* W = nav->GetCurrentNode();
  TGeoTube* S = dynamic_cast<TGeoTube*>(W->GetVolume()->GetShape());
  Double_t top[3] = {0, 0, S->GetDZ()};
  Double_t bot[3] = {0, 0, -S->GetDZ()};
  Double_t Gtop[3], Gbot[3];
  nav->LocalToMaster(top, Gtop);
  nav->LocalToMaster(bot, Gbot);
  vtop.SetXYZ(Gtop[0], Gtop[1], Gtop[2]);
  vbot.SetXYZ(Gbot[0], Gbot[1], Gbot[2]);
}

StrawSigmaSpatial()

Double_t strawtubes::StrawSigmaSpatial ( )

inline

Definition at line 55 of file strawtubes.h.

{ return sigma_spatial; }

StrawVdrift()

Double_t strawtubes::StrawVdrift ( )

inline

Definition at line 54 of file strawtubes.h.

{ return v_drift; }

Member Data Documentation

f_aperture_height

Double_t strawtubes::f_aperture_height

private

Aperture width (x)

Definition at line 69 of file strawtubes.h.

f_aperture_width

Double_t strawtubes::f_aperture_width

private

z-position of tracking station 4

Definition at line 68 of file strawtubes.h.

f_delta_z_layer

Double_t strawtubes::f_delta_z_layer

private

Offset (y) of straws between layers.

Definition at line 74 of file strawtubes.h.

f_delta_z_view

Double_t strawtubes::f_delta_z_view

private

Structure frame material.

Definition at line 78 of file strawtubes.h.

f_frame_material

TString strawtubes::f_frame_material

private

Sense wire thickness.

Definition at line 77 of file strawtubes.h.

f_inner_straw_diameter

Double_t strawtubes::f_inner_straw_diameter

private

Aperture height (y)

Definition at line 70 of file strawtubes.h.

f_offset_layer

Double_t strawtubes::f_offset_layer

private

Distance (y) between straws in a layer.

Definition at line 73 of file strawtubes.h.

f_outer_straw_diameter

Double_t strawtubes::f_outer_straw_diameter

private

Inner Straw diameter.

Definition at line 71 of file strawtubes.h.

f_station_height

Double_t strawtubes::f_station_height

private

Station envelope width (x)

Definition at line 80 of file strawtubes.h.

f_station_length

Double_t strawtubes::f_station_length

private

Station envelope height (y)

Definition at line 81 of file strawtubes.h.

f_station_width

Double_t strawtubes::f_station_width

private

Distance (z) between stereo views.

Definition at line 79 of file strawtubes.h.

f_straw_pitch

Double_t strawtubes::f_straw_pitch

private

Outer Straw diameter.

Definition at line 72 of file strawtubes.h.

f_T1_z

Double_t strawtubes::f_T1_z

private

Track information to be stored until the track leaves the active volume.

Definition at line 64 of file strawtubes.h.

f_T2_z

Double_t strawtubes::f_T2_z

private

z-position of tracking station 1

Definition at line 65 of file strawtubes.h.

f_T3_z

Double_t strawtubes::f_T3_z

private

z-position of tracking station 2

Definition at line 66 of file strawtubes.h.

f_T4_z

Double_t strawtubes::f_T4_z

private

z-position of tracking station 3

Definition at line 67 of file strawtubes.h.

f_view_angle

Double_t strawtubes::f_view_angle

private

Distance (z) between layers.

Definition at line 75 of file strawtubes.h.

f_wire_thickness

Double_t strawtubes::f_wire_thickness

private

Stereo view angle.

Definition at line 76 of file strawtubes.h.

fMedium

std::string strawtubes::fMedium

private

spatial resolution

Definition at line 84 of file strawtubes.h.

sigma_spatial

Double_t strawtubes::sigma_spatial

private

drift velocity

Definition at line 83 of file strawtubes.h.

v_drift

Double_t strawtubes::v_drift

private

Station envelope length (z)

Definition at line 82 of file strawtubes.h.

---

The documentation for this class was generated from the following files:

• strawtubes/strawtubes.h
• strawtubes/strawtubes.cxx