muon Class Reference

#include <muon.h>

Inheritance diagram for muon:

Collaboration diagram for muon:

Public Member Functions
	muon (const char *Name, Bool_t Active)
	muon ()
virtual	~muon ()
virtual void	Initialize ()
virtual Bool_t	ProcessHits (FairVolume *v=0)
virtual void	Register ()
virtual TClonesArray *	GetCollection (Int_t iColl) const
virtual void	Reset ()
void	SetZStationPositions (Double_t z0, Double_t z1, Double_t z2, Double_t z3)
void	SetZFilterPositions (Double_t z0, Double_t z1, Double_t z2)
void	SetActiveThickness (Double_t activeThickness)
void	SetFilterThickness (Double_t filterThickness)
void	SetXMax (Double_t xMax)
void	SetYMax (Double_t yMax)
void	ConstructGeometry ()
muonPoint *	AddHit (Int_t trackID, Int_t detID, TVector3 pos, TVector3 mom, Double_t time, Double_t length, Double_t eLoss, Int_t pdgcode)
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 ()

Private Member Functions
	muon (const muon &)
muon &	operator= (const muon &)
Int_t	InitMedium (const char *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
Double_t	fELoss
	length
Double_t	fM0z
	energy loss
Double_t	fM1z
	z-position of veto station
Double_t	fM2z
	z-position of tracking station 1
Double_t	fM3z
	z-position of tracking station 2
Double_t	fF0z
	z-position of tracking station 3
Double_t	fF1z
	z-position of veto station
Double_t	fF2z
	z-position of tracking station 1
Double_t	fFilterThickness
	z-position of tracking station 2
Double_t	fActiveThickness
Double_t	fXMax
Double_t	fYMax
TClonesArray *	fmuonPointCollection

Detailed Description

Definition at line 13 of file muon.h.

Constructor & Destructor Documentation

muon() [1/3]

muon::muon	(	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 49 of file muon.cxx.

  : FairDetector(name, active, kMuon),
    fTrackID(-1),
    fVolumeID(-1),
    fPos(),
    fMom(),
    fTime(-1.),
    fLength(-1.),
    fELoss(-1),
    fmuonPointCollection(new TClonesArray("muonPoint"))
{
}

muon() [2/3]

muon::muon

(

)

default constructor

Definition at line 36 of file muon.cxx.

  : FairDetector("muon", kTRUE, kMuon),
    fTrackID(-1),
    fVolumeID(-1),
    fPos(),
    fMom(),
    fTime(-1.),
    fLength(-1.),
    fELoss(-1),
    fmuonPointCollection(new TClonesArray("muonPoint"))
{
}

~muon()

muon::~muon ( )

virtual

destructor

Definition at line 62 of file muon.cxx.

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

muon() [3/3]

muon::muon ( const muon &  )

private

Member Function Documentation

AddHit()

muonPoint * muon::AddHit	(	Int_t	trackID,
		Int_t	detID,
		TVector3	pos,
		TVector3	mom,
		Double_t	time,
		Double_t	length,
		Double_t	eLoss,
		Int_t	pdgcode
	)

This method is an example of how to add your own point of type muonPoint to the clones array

Definition at line 254 of file muon.cxx.

{
  TClonesArray& clref = *fmuonPointCollection;
  Int_t size = clref.GetEntriesFast();
  // cout << "muon hit called"<< pos.z()<<endl;
  return new(clref[size]) muonPoint(trackID, detID, pos, mom,
         time, length, eLoss, pdgCode);
}

BeginEvent()

virtual void muon::BeginEvent ( )

inlinevirtual

Definition at line 86 of file muon.h.

{;}

BeginPrimary()

virtual void muon::BeginPrimary ( )

inlinevirtual

Definition at line 83 of file muon.h.

{;}

ConstructGeometry()

void muon::ConstructGeometry

(

)

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.

Definition at line 200 of file muon.cxx.

{
    TGeoVolume *top=gGeoManager->GetTopVolume();
    TGeoVolume *tMuon = new TGeoVolumeAssembly("MuonDetector");
 
    InitMedium("iron");
    InitMedium("Scintillator");
    
    TGeoMedium *Al =gGeoManager->GetMedium("Scintillator");
    TGeoMedium *A2 =gGeoManager->GetMedium("iron");
    
//    TGeoBBox *detbox1 = new TGeoBBox("detbox1", 250, 250, 10);
//    TGeoBBox *detbox2 = new TGeoBBox("detbox2", 245, 245, 10);
    
//    TGeoCompositeShape *detcomp1 = new TGeoCompositeShape("detcomp1", "detbox1-detbox2");    
 
//    TGeoVolume *detmu1 = new TGeoVolume("MuX", detcomp1, Al);
    TGeoVolume *muondet0 = gGeoManager->MakeBox("muondet0", Al, fXMax, fYMax, fActiveThickness);
    TGeoVolume *muondet1 = gGeoManager->MakeBox("muondet1", Al, fXMax, fYMax, fActiveThickness);
    TGeoVolume *muondet2 = gGeoManager->MakeBox("muondet2", Al, fXMax, fYMax, fActiveThickness); 
    TGeoVolume *muondet3 = gGeoManager->MakeBox("muondet3", Al, fXMax, fYMax, fActiveThickness); 
    
    TGeoVolume *muonfilter = gGeoManager->MakeBox("muonfilter", A2, fXMax, fYMax, fFilterThickness);
// 10cm iron to shield against backsplash from cavern
    TGeoVolume *muonshield = gGeoManager->MakeBox("muonshield", A2, fXMax, fYMax, 5.);
 
    AddSensitiveVolume(muondet0);
    AddSensitiveVolume(muondet1);
    AddSensitiveVolume(muondet2);
    AddSensitiveVolume(muondet3);
    muondet0->SetLineColor(kGreen);
    muondet1->SetLineColor(kGreen);
    muondet2->SetLineColor(kGreen);
    muondet3->SetLineColor(kGreen);
    muonfilter->SetLineColor(kBlue);
    Double_t zStartMuon = fM0z;
    Double_t totLength = fM3z-fM0z+2*fActiveThickness+15.;
    Double_t relPos = zStartMuon-fActiveThickness+totLength/2.;
    tMuon->AddNode(muondet0, 1, new TGeoTranslation(0, 0, fM0z-relPos));
    tMuon->AddNode(muonfilter, 0, new TGeoTranslation(0, 0, fF0z-relPos));
    tMuon->AddNode(muondet1, 1, new TGeoTranslation(0, 0, fM1z-relPos));
    tMuon->AddNode(muonfilter, 1, new TGeoTranslation(0, 0, fF1z-relPos));
    tMuon->AddNode(muondet2, 1, new TGeoTranslation(0, 0, fM2z-relPos));
    tMuon->AddNode(muonfilter, 2, new TGeoTranslation(0, 0, fF2z-relPos));
    tMuon->AddNode(muondet3, 1, new TGeoTranslation(0, 0, fM3z-relPos));
    tMuon->AddNode(muonshield, 1, new TGeoTranslation(0, 0, fM3z+fActiveThickness+10.-relPos));
          //finish assembly and position
    top->AddNode(tMuon, 1, new TGeoTranslation(0, 0,relPos));
}

CopyClones()

virtual void muon::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 77 of file muon.h.

                                             {;}

EndOfEvent()

void muon::EndOfEvent ( )

virtual

Definition at line 136 of file muon.cxx.

{
 
  fmuonPointCollection->Clear();
 
}

FinishPrimary()

virtual void muon::FinishPrimary ( )

inlinevirtual

Definition at line 81 of file muon.h.

{;}

FinishRun()

virtual void muon::FinishRun ( )

inlinevirtual

Definition at line 82 of file muon.h.

{;}

GetCollection()

TClonesArray * muon::GetCollection ( Int_t  iColl ) const

virtual

Gets the produced collections

Definition at line 160 of file muon.cxx.

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

Initialize()

void muon::Initialize ( )

virtual

Initialization of the detector is done here

Definition at line 70 of file muon.cxx.

{
  FairDetector::Initialize();
//  FairRuntimeDb* rtdb= FairRun::Instance()->GetRuntimeDb();
//  muonGeoPar* par=(muonGeoPar*)(rtdb->getContainer("muonGeoPar"));
}

InitMedium()

Int_t muon::InitMedium ( const char *  name )

private

Definition at line 77 of file muon.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);
     return -1111;
   }
   TGeoMedium* medium=gGeoManager->GetMedium(name);
  if (medium!=NULL)
    return ShipMedium->getMediumIndex();
 
  return geoBuild->createMedium(ShipMedium);
}

operator=()

muon & muon::operator= ( const muon &  )

private

PostTrack()

virtual void muon::PostTrack ( )

inlinevirtual

Definition at line 84 of file muon.h.

{;}

PreTrack()

virtual void muon::PreTrack ( )

inlinevirtual

Definition at line 85 of file muon.h.

{;}

ProcessHits()

Bool_t muon::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 97 of file muon.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 muonPoint at exit of active volume
  if ( gMC->IsTrackExiting()    ||
       gMC->IsTrackStop()       ||
       gMC->IsTrackDisappeared()   ) {
    fTrackID  = gMC->GetStack()->GetCurrentTrackNumber();
    //fVolumeID = vol->getMCid();
    //cout << "muon proc "<< fVolumeID<<" "<<vol->GetName()<<" "<<vol->getVolumeId() <<endl;
    //cout << " "<< gGeoManager->FindVolumeFast(vol->GetName())->GetNumber()<<endl;
    fVolumeID = gGeoManager->FindVolumeFast(vol->GetName())->GetNumber();
    if (fELoss == 0. ) { return kFALSE; }
    TParticle* p=gMC->GetStack()->GetCurrentTrack();
    Int_t pdgCode = p->GetPdgCode();
    AddHit(fTrackID, fVolumeID, TVector3(fPos.X(),  fPos.Y(),  fPos.Z()),
           TVector3(fMom.Px(), fMom.Py(), fMom.Pz()), fTime, fLength,
           fELoss,pdgCode);
 
    // Increment number of muon det points in TParticle
    ShipStack* stack = (ShipStack*) gMC->GetStack();
    stack->AddPoint(kMuon);
  }
 
  return kTRUE;
}

Register()

void muon::Register ( )

virtual

Registers the produced collections in FAIRRootManager.

This will create a branch in the output tree called muonPoint, setting the last parameter to kFALSE means: this collection will not be written to the file, it will exist only during the simulation.

Definition at line 145 of file muon.cxx.

{
 
  FairRootManager::Instance()->Register("muonPoint", "muon",
                                        fmuonPointCollection, kTRUE);
 
}

Reset()

void muon::Reset ( )

virtual

has to be called after each event to reset the containers

Definition at line 166 of file muon.cxx.

{
  fmuonPointCollection->Clear();
}

SetActiveThickness()

void muon::SetActiveThickness

(

Double_t

activeThickness

)

Definition at line 184 of file muon.cxx.

{
  fActiveThickness=activeThickness;
}

SetFilterThickness()

void muon::SetFilterThickness

(

Double_t

filterThickness

)

Definition at line 188 of file muon.cxx.

{
  fFilterThickness=filterThickness;
}

SetSpecialPhysicsCuts()

virtual void muon::SetSpecialPhysicsCuts ( )

inlinevirtual

Definition at line 79 of file muon.h.

{;}

SetXMax()

void muon::SetXMax

(

Double_t

xMax

)

Definition at line 192 of file muon.cxx.

{
  fXMax=xMax;
}

SetYMax()

void muon::SetYMax

(

Double_t

yMax

)

Definition at line 196 of file muon.cxx.

{
  fYMax=yMax;
}

SetZFilterPositions()

void muon::SetZFilterPositions	(	Double_t	z0,
		Double_t	z1,
		Double_t	z2
	)

Definition at line 178 of file muon.cxx.

{
  fF0z=z0;
 fF1z=z1;
 fF2z=z2;
}

SetZStationPositions()

void muon::SetZStationPositions	(	Double_t	z0,
		Double_t	z1,
		Double_t	z2,
		Double_t	z3
	)

Definition at line 171 of file muon.cxx.

{
  fM0z=z0;
 fM1z=z1;
 fM2z=z2;
 fM3z=z3;
}

Member Data Documentation

fActiveThickness

Double_t muon::fActiveThickness

private

Definition at line 110 of file muon.h.

fELoss

Double_t muon::fELoss

private

length

Definition at line 100 of file muon.h.

fF0z

Double_t muon::fF0z

private

z-position of tracking station 3

Definition at line 105 of file muon.h.

fF1z

Double_t muon::fF1z

private

z-position of veto station

Definition at line 106 of file muon.h.

fF2z

Double_t muon::fF2z

private

z-position of tracking station 1

Definition at line 107 of file muon.h.

fFilterThickness

Double_t muon::fFilterThickness

private

z-position of tracking station 2

Definition at line 109 of file muon.h.

fLength

Double_t muon::fLength

private

time

Definition at line 99 of file muon.h.

fM0z

Double_t muon::fM0z

private

energy loss

Definition at line 101 of file muon.h.

fM1z

Double_t muon::fM1z

private

z-position of veto station

Definition at line 102 of file muon.h.

fM2z

Double_t muon::fM2z

private

z-position of tracking station 1

Definition at line 103 of file muon.h.

fM3z

Double_t muon::fM3z

private

z-position of tracking station 2

Definition at line 104 of file muon.h.

fMom

TLorentzVector muon::fMom

private

position at entrance

Definition at line 97 of file muon.h.

fmuonPointCollection

TClonesArray* muon::fmuonPointCollection

private

container for data points

Definition at line 115 of file muon.h.

fPos

TLorentzVector muon::fPos

private

volume id

Definition at line 96 of file muon.h.

fTime

Double_t muon::fTime

private

momentum at entrance

Definition at line 98 of file muon.h.

fTrackID

Int_t muon::fTrackID

private

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

Definition at line 94 of file muon.h.

fVolumeID

Int_t muon::fVolumeID

private

track index

Definition at line 95 of file muon.h.

fXMax

Double_t muon::fXMax

private

Definition at line 111 of file muon.h.

fYMax

Double_t muon::fYMax

private

Definition at line 112 of file muon.h.

---

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

• muon/muon.h
• muon/muon.cxx