#include <boxTarget.h>

Inheritance diagram for boxTarget:

Collaboration diagram for boxTarget:

[legend]

Public Member Functions
	boxTarget (const char *Name, Bool_t Active)

	boxTarget ()

virtual	~boxTarget ()

virtual void	Initialize ()

virtual Bool_t	ProcessHits (FairVolume *v=0)

virtual void	Register ()

virtual TClonesArray *	GetCollection (Int_t iColl) const

virtual void	Reset ()

void	ConstructGeometry ()

virtual void	CopyClones (TClonesArray cl1, TClonesArray cl2, Int_t offset)

virtual void	SetSpecialPhysicsCuts ()

virtual void	EndOfEvent ()

virtual void	FinishPrimary ()

virtual void	FinishRun ()

virtual void	BeginPrimary ()

virtual void	PostTrack ()

virtual void	PreTrack ()

virtual void	BeginEvent ()

vetoPoint *	AddHit (Int_t trackID, Int_t detID, TVector3 pos, TVector3 mom, Double_t time, Double_t length, Double_t eLoss, Int_t pdgcode, TVector3 Lpos, TVector3 Lmom)

void	SetTarget (TString material, Float_t L, Bool_t choice)

Private Member Functions
Int_t	InitMedium (TString name)

Private Attributes
Int_t	fTrackID

Int_t	fVolumeID
	track index

TLorentzVector	fPos
	volume id

TLorentzVector	fMom
	position at entrance

Double_t	fTime
	momentum at entrance

Double_t	fLength
	time

TString	fTargetMaterial
	track length

Float_t	fTargetL
	target material

Int_t	index
	target length

Bool_t	fBox

TClonesArray *	fboxTargetPointCollection

Detailed Description

Definition at line 15 of file boxTarget.h.

Constructor & Destructor Documentation

◆ boxTarget() [1/2]

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

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

◆ boxTarget() [2/2]

boxTarget::boxTarget ( )

default constructor

Definition at line 37 of file boxTarget.cxx.

  : FairDetector("boxTarget", kTRUE, kVETO),
    fTrackID(-1),
    fVolumeID(-1),
    fPos(),
    fMom(),
    fTime(-1.),
    fLength(-1.),
    fTargetMaterial("tungsten"),
    fTargetL(1.),
    fboxTargetPointCollection(new TClonesArray("vetoPoint"))
{}

◆ ~boxTarget()

boxTarget::~boxTarget ( )

virtual

destructor

Definition at line 50 of file boxTarget.cxx.

{
  if (fboxTargetPointCollection) {
    fboxTargetPointCollection->Delete();
    delete fboxTargetPointCollection;
  }
}

Member Function Documentation

◆ AddHit()

vetoPoint * boxTarget::AddHit	(	Int_t	trackID,
		Int_t	detID,
		TVector3	pos,
		TVector3	mom,
		Double_t	time,
		Double_t	length,
		Double_t	eLoss,
		Int_t	pdgcode,
		TVector3	Lpos,
		TVector3	Lmom
	)

Definition at line 206 of file boxTarget.cxx.

{
  TClonesArray& clref = *fboxTargetPointCollection;
  Int_t size = clref.GetEntriesFast();
  return new(clref[size]) vetoPoint(trackID, detID, pos, mom,
         time, length, eLoss, pdgCode,Lpos,Lmom);
}

◆ BeginEvent()

virtual void boxTarget::BeginEvent ( )

inlinevirtual

Definition at line 65 of file boxTarget.h.

65{;}

◆ BeginPrimary()

virtual void boxTarget::BeginPrimary ( )

inlinevirtual

Definition at line 62 of file boxTarget.h.

62{;}

◆ ConstructGeometry()

void boxTarget::ConstructGeometry ( )

Create the detector geometry

Definition at line 132 of file boxTarget.cxx.

{
   Double_t cm  = 1;       // cm
   Double_t m   = 100*cm;  //  m
   Double_t mm  = 0.1*cm;  //  mm
 
   TGeoVolume *top=gGeoManager->FindVolumeFast("Detector");
    if(!top)  LOG(ERROR) << "no Detector volume found " ;
 
   InitMedium(fTargetMaterial);
   TGeoMedium *TargetMaterial = gGeoManager->GetMedium(fTargetMaterial);
   InitMedium("vacuums");
   TGeoMedium *vac  = gGeoManager->GetMedium("vacuums");
   InitMedium("Concrete");
   TGeoMedium *concrete = gGeoManager->GetMedium("Concrete");
    if (fBox){
    // for studying absorption length
      TGeoVolume* target  = gGeoManager->MakeBox("Target",TargetMaterial,199.*cm,199.*cm,(fTargetL/2)*cm);
      top->AddNode(target, 1, new TGeoTranslation(0, 0, (fTargetL/2)*cm));
 
      TGeoVolume *sensPlane = gGeoManager->MakeBox("sensPlane",vac,199.*cm,199.*cm,1.*mm);
      sensPlane->SetLineColor(kGreen);
      top->AddNode(sensPlane, 13, new TGeoTranslation(0, 0, (fTargetL+ 0.11)*cm));
      AddSensitiveVolume(sensPlane);
      AddSensitiveVolume(target);
    }else{
    // more realistic setup of coldBox
    //  x1=-1140 mm, x2= 400 mm
    //  y1=100 mm, y2=1114 mm
    //  z1=50 mm, z2=1500mm
    // 20cm x 2cm hole for SciFi cables 
    // 20cm x 20cm hole for cooling pipes
      float x1=-1140;
      float x2= 400;
      float y1=100;
      float y2=1114;
      float z1=50;
      float z2=1500;
      float dx = (x2-x1)/10./2.;
      float dy = (y2-y1)/10./2.;
      float dz = (z2-z1)/10./2.;
      float d = fTargetL;
      float rAir  = 10.;
      float cablesDx = 10;
      float cablesDy = 1;
 
      TGeoBBox* box_I   = new TGeoBBox("box_I",dx,dy,dz);
      double origin[3] = {0,d/2,0};
      TGeoBBox* box_O = new TGeoBBox("box_O",dx+d,dy+d/2.,dz+d,origin);
      TGeoTube* hole_air = new TGeoTube("hole_air",0.,rAir,dx);
      TGeoRotation* RF = new TGeoRotation("R_airTube");
      RF->SetAngles(0.,0.,90.);
      TGeoCombiTrans* CombiTrans1 = new TGeoCombiTrans("T_AirTube",-dx/2,0,10.,RF);
      CombiTrans1->RegisterYourself();
      TGeoCombiTrans* CombiTrans2 = new TGeoCombiTrans("T_CablesDuct",dx/2,0,10.,RF);
      CombiTrans2->RegisterYourself();
      TGeoBBox* hole_cables = new TGeoBBox("hole_cables",cablesDx,cablesDy,dx);
      double floor_thickness  = 25.;
      double originF[3] = {0,-dy-floor_thickness,0};
      TGeoBBox* floor   = new TGeoBBox("floor",1.5*dx,floor_thickness,1.5*dz,originF);
      TGeoCompositeShape* box = new TGeoCompositeShape("box","box_O-box_I-hole_air:T_AirTube-hole_cables:T_CablesDuct");
      TGeoVolume *volBox = new TGeoVolume("vbox",box,TargetMaterial);
      TGeoVolume *volFloor = new TGeoVolume("vfloor",floor,concrete);
      TGeoVolume *sensBox = new TGeoVolume("sensBox",box_I,vac);
      sensBox->SetLineColor(kGreen);
      top->AddNode(sensBox, 13);
      AddSensitiveVolume(sensBox);
      top->AddNode(volBox, 1);
      AddSensitiveVolume(volBox);
      top->AddNode(volFloor, 2);
    }
 
}

◆ CopyClones()

virtual void boxTarget::CopyClones	(	TClonesArray *	cl1,
		TClonesArray *	cl2,
		Int_t	offset
	)

inlinevirtual

The following methods can be implemented if you need to make any optional action in your detector during the transport.

Definition at line 56 of file boxTarget.h.

57 {;}

◆ EndOfEvent()

void boxTarget::EndOfEvent ( )

virtual

Definition at line 122 of file boxTarget.cxx.

{
 
  fboxTargetPointCollection->Clear();
 
}

◆ FinishPrimary()

virtual void boxTarget::FinishPrimary ( )

inlinevirtual

Definition at line 60 of file boxTarget.h.

60{;}

◆ FinishRun()

void boxTarget::FinishRun ( )

virtual

Definition at line 129 of file boxTarget.cxx.

129 {

130}

◆ GetCollection()

TClonesArray * boxTarget::GetCollection ( Int_t iColl ) const

virtual

Gets the produced collections

Definition at line 224 of file boxTarget.cxx.

{
  if (iColl == 0) { return fboxTargetPointCollection; }
  else { return NULL; }
}

◆ Initialize()

void boxTarget::Initialize ( )

virtual

Initialization of the detector is done here

Definition at line 117 of file boxTarget.cxx.

{
  FairDetector::Initialize();
}

◆ InitMedium()

Int_t boxTarget::InitMedium ( TString name )

private

Definition at line 99 of file boxTarget.cxx.

{
   static FairGeoLoader *geoLoad=FairGeoLoader::Instance();
   static FairGeoInterface *geoFace=geoLoad->getGeoInterface();
   static FairGeoMedia *media=geoFace->getMedia();
   static FairGeoBuilder *geoBuild=geoLoad->getGeoBuilder();
 
   FairGeoMedium *ShipMedium=media->getMedium(name);
 
   if (!ShipMedium)
     Fatal("InitMedium","Material %s not defined in media file.", name.Data());
   TGeoMedium* medium=gGeoManager->GetMedium(name);
   if (medium)
     return ShipMedium->getMediumIndex();
   return geoBuild->createMedium(ShipMedium);
}

◆ PostTrack()

virtual void boxTarget::PostTrack ( )

inlinevirtual

Definition at line 63 of file boxTarget.h.

63{;}

◆ PreTrack()

virtual void boxTarget::PreTrack ( )

inlinevirtual

Definition at line 64 of file boxTarget.h.

64{;}

◆ ProcessHits()

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

virtual

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

This method is called from the MC stepping

Definition at line 58 of file boxTarget.cxx.

{
  if ( gMC->IsTrackEntering() ) {
    fTime   = gMC->TrackTime() * 1.0e09;
    fLength = gMC->TrackLength();
    gMC->TrackPosition(fPos);
    gMC->TrackMomentum(fMom);
    LOG(DEBUG) <<"track enters"<<gMC->GetStack()->GetCurrentTrack()->GetPdgCode();
  }
  if (!gMC->IsTrackEntering() ) { 
        LOG(DEBUG) <<"track is not entering"<<gMC->GetStack()->GetCurrentTrack()->GetPdgCode();
  }
  if ( gMC->IsTrackExiting()    ||
       gMC->IsTrackStop()       ||
       gMC->IsTrackDisappeared()   ) {
       LOG(DEBUG) <<"track stopped"<<gMC->GetStack()->GetCurrentTrack()->GetPdgCode();
 
    fTrackID  = gMC->GetStack()->GetCurrentTrackNumber();
    Int_t veto_uniqueId;
    gMC->CurrentVolID(veto_uniqueId);
 
    TParticle* p=gMC->GetStack()->GetCurrentTrack();
    Int_t pdgCode = p->GetPdgCode();
    TLorentzVector Pos;
    gMC->TrackPosition(Pos);
    TLorentzVector Mom;
    gMC->TrackMomentum(Mom);
    Double_t xmean = (fPos.X()+Pos.X())/2. ;
    Double_t ymean = (fPos.Y()+Pos.Y())/2. ;
    Double_t zmean = (fPos.Z()+Pos.Z())/2. ;
    AddHit(fTrackID, veto_uniqueId, TVector3(xmean, ymean,  zmean),
           TVector3(fMom.Px(), fMom.Py(), fMom.Pz()), fTime, fLength,
           0.,pdgCode,TVector3(Pos.X(), Pos.Y(), Pos.Z()),TVector3(Mom.Px(), Mom.Py(), Mom.Pz()) );
    // Increment number of veto det points in TParticle
    ShipStack* stack = (ShipStack*) gMC->GetStack();
    stack->AddPoint(kVETO);
  }
  return kTRUE;
}

◆ Register()

void boxTarget::Register ( )

virtual

Registers the produced collections in FAIRRootManager.

Definition at line 217 of file boxTarget.cxx.

{
 
  FairRootManager::Instance()->Register("vetoPoint", "veto",
                                        fboxTargetPointCollection, kTRUE);
}

◆ Reset()

void boxTarget::Reset ( )

virtual

has to be called after each event to reset the containers

Definition at line 229 of file boxTarget.cxx.

{
  fboxTargetPointCollection->Clear();
}

◆ SetSpecialPhysicsCuts()

virtual void boxTarget::SetSpecialPhysicsCuts ( )

inlinevirtual

Definition at line 58 of file boxTarget.h.

58{;}

◆ SetTarget()

void boxTarget::SetTarget	(	TString	material,
		Float_t	L,
		Bool_t	choice
	)

inline

Definition at line 71 of file boxTarget.h.

71{fTargetMaterial = material; fTargetL=L; fBox=choice;}

conditionsDatabase.demo.drifttubeMonitoring.L

int L

Definition drifttubeMonitoring.py:297

study_GammaConv.material

dict material

Definition study_GammaConv.py:24

Member Data Documentation

◆ fBox

Bool_t boxTarget::fBox

private

Definition at line 86 of file boxTarget.h.

◆ fboxTargetPointCollection

TClonesArray* boxTarget::fboxTargetPointCollection

private

container for data points

Definition at line 88 of file boxTarget.h.

◆ fLength

Double_t boxTarget::fLength

private

time

Definition at line 82 of file boxTarget.h.

◆ fMom

TLorentzVector boxTarget::fMom

private

position at entrance

Definition at line 80 of file boxTarget.h.

◆ fPos

TLorentzVector boxTarget::fPos

private

volume id

Definition at line 79 of file boxTarget.h.

◆ fTargetL

Float_t boxTarget::fTargetL

private

target material

Definition at line 84 of file boxTarget.h.

◆ fTargetMaterial

TString boxTarget::fTargetMaterial

private

track length

Definition at line 83 of file boxTarget.h.

◆ fTime

Double_t boxTarget::fTime

private

momentum at entrance

Definition at line 81 of file boxTarget.h.

◆ fTrackID

Int_t boxTarget::fTrackID

private

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

Definition at line 77 of file boxTarget.h.

◆ fVolumeID

Int_t boxTarget::fVolumeID

private

track index

Definition at line 78 of file boxTarget.h.

◆ index

Int_t boxTarget::index

private

target length

Definition at line 85 of file boxTarget.h.

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

shipLHC/boxTarget.h
shipLHC/boxTarget.cxx

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ boxTarget() [1/2]

◆ boxTarget() [2/2]

◆ ~boxTarget()

Member Function Documentation

◆ AddHit()

◆ BeginEvent()

◆ BeginPrimary()

◆ ConstructGeometry()

◆ CopyClones()

◆ EndOfEvent()

◆ FinishPrimary()

◆ FinishRun()

◆ GetCollection()

◆ Initialize()

◆ InitMedium()

◆ PostTrack()

◆ PreTrack()

◆ ProcessHits()

◆ Register()

◆ Reset()

◆ SetSpecialPhysicsCuts()

◆ SetTarget()

Member Data Documentation

◆ fBox

◆ fboxTargetPointCollection

◆ fLength

◆ fMom

◆ fPos

◆ fTargetL

◆ fTargetMaterial

◆ fTime

◆ fTrackID

◆ fVolumeID

◆ index