#include <simpleTarget.h>

Inheritance diagram for simpleTarget:

Collaboration diagram for simpleTarget:

[legend]

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

	simpleTarget ()

virtual	~simpleTarget ()

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	SetEnergyCut (Float_t emax)

void	SetOnlyMuons ()

void	SetFastMuon ()

void	SetParameters (TString m, Float_t l, Float_t z)

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

Double_t	fzPos
	length

Double_t	fThick
	zPos

Double_t	fELoss
	thickness

Double_t	fTotalEloss

TString	fMaterial

Float_t	EMax
	material

Bool_t	fOnlyMuons
	max energy to transport

Bool_t	fFastMuon

TFile *	fout
	for fast processing

TClonesArray *	fsimpleTargetPointCollection

Detailed Description

Definition at line 14 of file simpleTarget.h.

Constructor & Destructor Documentation

◆ simpleTarget() [1/2]

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

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

◆ simpleTarget() [2/2]

simpleTarget::simpleTarget ( )

default constructor

Definition at line 35 of file simpleTarget.cxx.

  : FairDetector("simpleTarget", kTRUE, kVETO),
    fTrackID(-1),
    fVolumeID(-1),
    fPos(),
    fMom(),
    fTime(-1.),
    fLength(-1.),
    fThick(-1.),
    fOnlyMuons(kFALSE),
    fFastMuon(kFALSE),
    fzPos(3E8),
    fTotalEloss(0),
    fsimpleTargetPointCollection(new TClonesArray("vetoPoint"))
{}

◆ ~simpleTarget()

simpleTarget::~simpleTarget ( )

virtual

destructor

Definition at line 51 of file simpleTarget.cxx.

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

Member Function Documentation

◆ AddHit()

vetoPoint * simpleTarget::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 146 of file simpleTarget.cxx.

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

◆ BeginEvent()

virtual void simpleTarget::BeginEvent ( )

inlinevirtual

Definition at line 64 of file simpleTarget.h.

64{;}

◆ BeginPrimary()

virtual void simpleTarget::BeginPrimary ( )

inlinevirtual

Definition at line 61 of file simpleTarget.h.

61{;}

◆ ConstructGeometry()

void simpleTarget::ConstructGeometry ( )

Create the detector geometry

Definition at line 128 of file simpleTarget.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(fMaterial);
   TGeoMedium* mat=gGeoManager->GetMedium(fMaterial);
   if (mat==NULL)
     geoBuild->createMedium(ShipMedium);
   mat =gGeoManager->GetMedium(fMaterial);
   TGeoVolume *top=gGeoManager->GetTopVolume();
   TGeoVolume *target = gGeoManager->MakeBox("target",mat,10.*100.,10.*100.,fThick);
   top->AddNode(target, 1, new TGeoTranslation(0, 0, fzPos));
   AddSensitiveVolume(target);
}

◆ CopyClones()

virtual void simpleTarget::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 55 of file simpleTarget.h.

56 {;}

◆ EndOfEvent()

void simpleTarget::EndOfEvent ( )

virtual

Definition at line 121 of file simpleTarget.cxx.

{
  fsimpleTargetPointCollection->Clear();
  fTotalEloss=0;
}

◆ FinishPrimary()

virtual void simpleTarget::FinishPrimary ( )

inlinevirtual

Definition at line 59 of file simpleTarget.h.

59{;}

◆ FinishRun()

virtual void simpleTarget::FinishRun ( )

inlinevirtual

Definition at line 60 of file simpleTarget.h.

60{;}

◆ GetCollection()

TClonesArray * simpleTarget::GetCollection ( Int_t iColl ) const

virtual

Gets the produced collections

Definition at line 171 of file simpleTarget.cxx.

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

◆ Initialize()

void simpleTarget::Initialize ( )

virtual

Initialization of the detector is done here

Definition at line 114 of file simpleTarget.cxx.

{
  FairDetector::Initialize();
  TSeqCollection* fileList=gROOT->GetListOfFiles();
  fout = ((TFile*)fileList->At(0));
}

◆ PostTrack()

virtual void simpleTarget::PostTrack ( )

inlinevirtual

Definition at line 62 of file simpleTarget.h.

62{;}

◆ PreTrack()

void simpleTarget::PreTrack ( )

virtual

Definition at line 157 of file simpleTarget.cxx.

                           {
    if (!fFastMuon){return;}
    if (TMath::Abs(gMC->TrackPid())!=13){
        gMC->StopTrack();
    }
}

◆ ProcessHits()

Bool_t simpleTarget::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 59 of file simpleTarget.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 vetoPoint at exit of active volume
  if ( gMC->IsTrackExiting()    ||
       gMC->IsTrackStop()       ||
       gMC->IsTrackDisappeared()   ) {
    if (fELoss == 0. ) { return kFALSE; }
   // if (fOnlyMuons and fPos.Z()<140){ used for LiquidKrypton study
    if (fOnlyMuons){
       fTotalEloss+=fELoss;
       TClonesArray& clref = *fsimpleTargetPointCollection;
       new (clref[0]) vetoPoint(0, 0, TVector3(0, 0,  0), TVector3(0, 0,  0),
         0, 0, fTotalEloss, 0,TVector3(0, 0,  0), TVector3(0, 0,  0));
       return kTRUE;}
    
    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. ;
    //cout << veto_uniqueId << " :(" << xmean << ", " << ymean << ", " << zmean << "): " << fELoss << endl;
    AddHit(fTrackID, veto_uniqueId, TVector3(xmean, ymean,  zmean),
           TVector3(fMom.Px(), fMom.Py(), fMom.Pz()), fTime, fLength,
           fELoss,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 simpleTarget::Register ( )

virtual

Registers the produced collections in FAIRRootManager.

Definition at line 164 of file simpleTarget.cxx.

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

◆ Reset()

void simpleTarget::Reset ( )

virtual

has to be called after each event to reset the containers

Definition at line 176 of file simpleTarget.cxx.

{
  fsimpleTargetPointCollection->Clear();
}

◆ SetEnergyCut()

void simpleTarget::SetEnergyCut ( Float_t emax )

inline

Definition at line 70 of file simpleTarget.h.

70{EMax=emax;}// min energy to be copied to Geant4

simpleTarget::EMax

Float_t EMax

material

Definition simpleTarget.h:91

makeSNDGenieEvents.emax

emax

Definition makeSNDGenieEvents.py:140

◆ SetFastMuon()

void simpleTarget::SetFastMuon ( )

inline

Definition at line 72 of file simpleTarget.h.

72{fFastMuon=kTRUE;}

◆ SetOnlyMuons()

void simpleTarget::SetOnlyMuons ( )

inline

Definition at line 71 of file simpleTarget.h.

71{fOnlyMuons=kTRUE;}

◆ SetParameters()

void simpleTarget::SetParameters	(	TString	m,
		Float_t	l,
		Float_t	z
	)

inline

Definition at line 73 of file simpleTarget.h.

73{fMaterial=m;fThick=l;fzPos=z;}

m

Double_t m

Definition ShipBellField.cxx:20

ana_ShipMuon.l

dict l

Definition ana_ShipMuon.py:494

charm-geometry_config.z

z

Definition charm-geometry_config.py:7

◆ SetSpecialPhysicsCuts()

virtual void simpleTarget::SetSpecialPhysicsCuts ( )

inlinevirtual

Definition at line 57 of file simpleTarget.h.

57{;}

Member Data Documentation

◆ EMax

Float_t simpleTarget::EMax

private

material

Definition at line 91 of file simpleTarget.h.

◆ fELoss

Double_t simpleTarget::fELoss

private

thickness

Definition at line 88 of file simpleTarget.h.

◆ fFastMuon

Bool_t simpleTarget::fFastMuon

private

Definition at line 93 of file simpleTarget.h.

◆ fLength

Double_t simpleTarget::fLength

private

time

Definition at line 85 of file simpleTarget.h.

◆ fMaterial

TString simpleTarget::fMaterial

private

Definition at line 90 of file simpleTarget.h.

◆ fMom

TLorentzVector simpleTarget::fMom

private

position at entrance

Definition at line 83 of file simpleTarget.h.

◆ fOnlyMuons

Bool_t simpleTarget::fOnlyMuons

private

max energy to transport

Definition at line 92 of file simpleTarget.h.

◆ fout

TFile* simpleTarget::fout

private

for fast processing

Definition at line 94 of file simpleTarget.h.

◆ fPos

TLorentzVector simpleTarget::fPos

private

volume id

Definition at line 82 of file simpleTarget.h.

◆ fsimpleTargetPointCollection

TClonesArray* simpleTarget::fsimpleTargetPointCollection

private

container for data points

Definition at line 96 of file simpleTarget.h.

◆ fThick

Double_t simpleTarget::fThick

private

zPos

Definition at line 87 of file simpleTarget.h.

◆ fTime

Double_t simpleTarget::fTime

private

momentum at entrance

Definition at line 84 of file simpleTarget.h.

◆ fTotalEloss

Double_t simpleTarget::fTotalEloss

private

Definition at line 89 of file simpleTarget.h.

◆ fTrackID

Int_t simpleTarget::fTrackID

private

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

Definition at line 80 of file simpleTarget.h.

◆ fVolumeID

Int_t simpleTarget::fVolumeID

private

track index

Definition at line 81 of file simpleTarget.h.

◆ fzPos

Double_t simpleTarget::fzPos

private

length

Definition at line 86 of file simpleTarget.h.

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

muonShieldOptimization/simpleTarget.h
muonShieldOptimization/simpleTarget.cxx

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ simpleTarget() [1/2]

◆ simpleTarget() [2/2]

◆ ~simpleTarget()

Member Function Documentation

◆ AddHit()

◆ BeginEvent()

◆ BeginPrimary()

◆ ConstructGeometry()

◆ CopyClones()

◆ EndOfEvent()

◆ FinishPrimary()

◆ FinishRun()

◆ GetCollection()

◆ Initialize()

◆ PostTrack()

◆ PreTrack()

◆ ProcessHits()

◆ Register()

◆ Reset()

◆ SetEnergyCut()

◆ SetFastMuon()

◆ SetOnlyMuons()

◆ SetParameters()

◆ SetSpecialPhysicsCuts()

Member Data Documentation

◆ EMax

◆ fELoss

◆ fFastMuon

◆ fLength

◆ fMaterial

◆ fMom

◆ fOnlyMuons

◆ fout

◆ fPos

◆ fsimpleTargetPointCollection

◆ fThick

◆ fTime

◆ fTotalEloss

◆ fTrackID

◆ fVolumeID

◆ fzPos