muonHit Class Reference

#include <muonHit.h>

Inheritance diagram for muonHit:

Collaboration diagram for muonHit:

Public Member Functions
	muonHit ()
	muonHit (Int_t detID, Float_t digi, Bool_t isValid)
	muonHit (muonPoint *p, Double_t t0)
Int_t	DetIDfromXYZ (TVector3 p)
TVector3	XYZfromDetID (Int_t detID)
virtual	~muonHit ()
virtual void	Print () const
TVector3	getPos ()
Bool_t	isValid () const
Double_t	SetMuonTimeRes (Double_t mcTime)
void	setValidity (Bool_t isValid)
Public Member Functions inherited from ShipHit
	ShipHit ()
	ShipHit (Int_t detID, Float_t digi)
virtual	~ShipHit ()
Double_t	GetDigi () const
Int_t	GetDetectorID () const
void	SetDigi (Float_t d)
void	SetDetectorID (Int_t detID)
virtual void	Print (const Option_t *opt="") const
template<class Archive >
void	serialize (Archive &ar, const unsigned int version)

Private Member Functions
	muonHit (const muonHit &point)
muonHit	operator= (const muonHit &point)
void	stInit ()

Private Attributes
Float_t	flag
	flag
Bool_t	hisV

Static Private Attributes
static bool	onlyOnce =false

Additional Inherited Members
Protected Member Functions inherited from ShipHit
	ClassDef (ShipHit, 1)
Protected Attributes inherited from ShipHit
Float_t	fdigi
	digitized detector hit
Int_t	fDetectorID
	Detector unique identifier.

Detailed Description

Definition at line 11 of file muonHit.h.

Constructor & Destructor Documentation

muonHit() [1/4]

muonHit::muonHit

(

)

Default constructor

Definition at line 29 of file muonHit.cxx.

  : ShipHit()
{
 
}

muonHit() [2/4]

muonHit::muonHit	(	Int_t	detID,
		Float_t	digi,
		Bool_t	isValid
	)

Constructor with arguments

Parameters

detID	Detector ID
digi	digitized/measured TDC
flag	True/False, false if there is another hit with smaller tdc

Definition at line 35 of file muonHit.cxx.

  : ShipHit(detID,digi) 
{
  //
  SetDigi(digi);
  setValidity(isV);
}

muonHit() [3/4]

muonHit::muonHit	(	muonPoint *	p,
		Double_t	t0
	)

Definition at line 43 of file muonHit.cxx.

  : ShipHit()
{
//
  TVector3 truePosition = TVector3( p->GetX(), p->GetY(),p->GetZ());
  fdigi = t0 + p->GetTime(); // + drift time,  propagation inside tile + tdc    
  SetDetectorID(DetIDfromXYZ(truePosition));
  SetDigi(SetMuonTimeRes(fdigi));
}

~muonHit()

muonHit::~muonHit ( )

virtual

Destructor

Definition at line 180 of file muonHit.cxx.

                  { 
}

muonHit() [4/4]

muonHit::muonHit ( const muonHit &  point )

private

Copy constructor

Member Function Documentation

DetIDfromXYZ()

Int_t muonHit::DetIDfromXYZ

(

TVector3

p

)

Definition at line 53 of file muonHit.cxx.

{  
    // needs some code to produce a unique detector ID
//
//                    tiles numbering example with present tile dimensions (X=10,Y=20cm):
//                                                            ---------------------------
//                                            muon station 1  |13540|13541...13598|13599|
//                                                           ---------------------------|
//                                           muon station 0  |03540|03541...03598|03599||
// How tiles numbering works:                                |03480|03481...03538|03539|.
// in each station tiles are numbered by rows starting       |                         ||
// from 0. Numbering begins at bottom of the station         ...........................|
// from left to right (looking the muon station from         |                         ||
// em calorimeter). At each tile is added the number         |00060|00061...00118|00119|
// n * 10ˆ4 where n is the station number (0 <= n <= 3).     |00000|00001...00058|00059|
//                                                           ---------------------------
// negative detID means ERROR.
//
    Int_t detID, nStat; // unique detector ID, station number
//
  if (!onlyOnce) {
    stInit();
    onlyOnce = true;
  }
    nStat = -1;
    for(Int_t i=0; i<muStZpos.size(); i++) {
     if (abs(p.Z() - muStZpos[i]) <= muStzMax[i]) { 
       nStat = i;
       break;
     }
    }
//
    if (nStat != -1) {
      detID=(Int_t)((p.X()+muStxMax[nStat])/tileXdim)+tileXn[nStat]
           *(Int_t)((p.Y()+muStyMax[nStat])/tileYdim)+10000*nStat;
    } else {
      detID = -1;
    }
//
    return detID;
}

getPos()

TVector3 muonHit::getPos ( )

inline

Definition at line 34 of file muonHit.h.

{return XYZfromDetID(fDetectorID);}

isValid()

Bool_t muonHit::isValid ( ) const

inline

Definition at line 35 of file muonHit.h.

{return hisV;}

operator=()

muonHit muonHit::operator= ( const muonHit &  point )

private

Print()

void muonHit::Print ( ) const

virtual

Output to screen

Definition at line 171 of file muonHit.cxx.

                          {
//
  cout << "-I- muonHit: muon hit " << " in detector " << fDetectorID << endl;
  cout << "  TDC " << fdigi << " ns" << endl;
 
}

SetMuonTimeRes()

Double_t muonHit::SetMuonTimeRes

(

Double_t

mcTime

)

Definition at line 164 of file muonHit.cxx.

                                                {
//
  TRandom3 *rand = new TRandom3(0);
  Double_t cTime = rand->Gaus(mcTime,muonTimeResSigma);
  delete rand;
  return cTime;
}

setValidity()

void muonHit::setValidity

(

Bool_t

isValid

)

Definition at line 178 of file muonHit.cxx.

{hisV = isV;}

stInit()

void muonHit::stInit ( )

private

Definition at line 133 of file muonHit.cxx.

{
//
//
    TGeoShape* muShape; TGeoBBox* muonBox;
    TGeoNavigator* nav = gGeoManager->GetCurrentNavigator();
    Double_t loc[3]={0,0,0}, global[3]={0,0,0};
//
    TString muDet = "cave/MuonDetector_1";
    nav->cd(muDet);
    TGeoNode* node = nav->GetCurrentNode();
    TObjArray* nodes =  node->GetVolume()->GetNodes();
//
    for (Int_t i = 0; i < nodes->GetEntries(); i++) {
      node = (TGeoNode*)nodes->At(i);
      Int_t muStNs = 0;
      if (TString(node->GetName()).Contains("muondet")) {
        nav->cd(muDet+"/"+node->GetName());
        TGeoVolume* volu = node->GetVolume();
        muShape = node->GetVolume()->GetShape();
        muonBox = (TGeoBBox*) muShape;
        muStxMax.push_back(muonBox->GetDX()); muStyMax.push_back(muonBox->GetDY()); 
        muStzMax.push_back(muonBox->GetDZ());
        nav->LocalToMaster(loc,global); muStZpos.push_back(global[2]);
        tileXn.push_back(2*muStxMax.at(muStNs)/tileXdim); 
        tileYn.push_back(2*muStyMax.at(muStNs)/tileYdim);
    muStNs++; // muon stations increment
      }
    }
}

XYZfromDetID()

TVector3 muonHit::XYZfromDetID

(

Int_t

detID

)

Definition at line 95 of file muonHit.cxx.

{
//
// The center of the tile XYZ coordinates are returned
//
// Negative Z coordinate means ERROR
//
    Int_t nStat; TVector3 p; // station number, center tile coordinate
//
  if (!onlyOnce) {
    stInit();
    onlyOnce = true;
  }
    nStat = -1;
    for(Int_t i=0; i<=muStZpos.size(); i++) {
      if (dID < 10000*i) { 
        nStat = i;
        break;
      }
    }
    if (nStat < 1) {// ERROR! dID too low (nStat=0) or too high (nStat=-1)
      p.SetXYZ(-999999,-999999,-999999);
      return p;
    }
 
//
    p.SetZ(muStZpos[--nStat]); // now must be 0 <= nstat <= (n. of stations - 1)
    dID -= 10000*(nStat);
//
    Int_t muXpos, muYpos;
    muXpos = tileXdim*((dID%tileXn[nStat])+0.5)-muStxMax[nStat];
    muYpos = tileYdim*((Int_t)(dID/tileXn[nStat])+0.5)-muStyMax[nStat];
    p.SetXYZ(muXpos, muYpos, muStZpos[nStat]);
//
    return p;
//
}

Member Data Documentation

flag

Float_t muonHit::flag

private

flag

Definition at line 45 of file muonHit.h.

hisV

Bool_t muonHit::hisV

private

Definition at line 50 of file muonHit.h.

onlyOnce

bool muonHit::onlyOnce =false

staticprivate

Definition at line 47 of file muonHit.h.

---

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

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