#include <splitcalHit.h>

Inheritance diagram for splitcalHit:

[legend]

Collaboration diagram for splitcalHit:

[legend]

Public Member Functions
	splitcalHit ()

	splitcalHit (Int_t detID, Float_t tdc)

	splitcalHit (splitcalPoint *p, Double_t t0)

virtual	~splitcalHit ()

virtual void	Print () const

Float_t	GetTDC () const

void	setInvalid ()

bool	isValid () const

std::string	GetPaddedString (int &id)

std::string	GetDetectorElementName (int &id)

void	Decoder (int &id, int &isPrecision, int &nLayer, int &nModuleX, int &nMdouleY, int &nStrip)

void	Decoder (std::string &encodedID, int &isPrecision, int &nLayer, int &nModuleX, int &nMdouleY, int &nStrip)

void	SetXYZ (double &x, double &y, double &z)

void	SetIDs (int &isPrecision, int &nLayer, int &nModuleX, int &nModuleY, int &nStrip)

void	SetEnergy (double &e)

void	UpdateEnergy (double e)

void	SetIsX (bool x)

void	SetIsY (bool y)

void	SetIsUsed (int u)

void	SetXYZErrors (double xError, double yError, double zError)

void	AddClusterIndex (int i)

void	AddEnergyWeight (double w)

double	GetX ()

double	GetY ()

double	GetZ ()

double	GetEnergy ()

double	GetEnergyForCluster (int i)

int	GetIsPrecisionLayer ()

int	GetLayerNumber ()

int	GetModuleXNumber ()

int	GetModuleYNumber ()

int	GetStripNumber ()

bool	IsX ()

bool	IsY ()

int	IsUsed ()

double	GetXError ()

double	GetYError ()

double	GetZError ()

std::vector< int >	GetClusterIndices ()

std::vector< double >	GetEnergyWeights ()

bool	IsShared ()

double	GetEnergyWeightForIndex (int index)

	splitcalHit (const splitcalHit &point)

splitcalHit	operator= (const splitcalHit &point)

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
	ClassDef (splitcalHit, 3)

Private Attributes
Float_t	flag

double	_x

double	_y

double	_z

double	_xError

double	_yError

double	_zError

double	_energy

int	_isPrecisionLayer

int	_nLayer

int	_nModuleX

int	_nModuleY

int	_nStrip

int	_isUsed

bool	_isX

bool	_isY

std::vector< double >	_vecEnergyWeights

std::vector< int >	_vecClusterIndices

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 10 of file splitcalHit.h.

Constructor & Destructor Documentation

◆ splitcalHit() [1/4]

splitcalHit::splitcalHit ( )

Default constructor

Definition at line 24 of file splitcalHit.cxx.

  : ShipHit()
{
 flag = true;
}

◆ splitcalHit() [2/4]

splitcalHit::splitcalHit	(	Int_t	detID,
		Float_t	tdc
	)

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 30 of file splitcalHit.cxx.

  : ShipHit(detID,tdc)
{
 flag = true;
}

◆ splitcalHit() [3/4]

splitcalHit::splitcalHit	(	splitcalPoint *	p,
		Double_t	t0
	)

Definition at line 36 of file splitcalHit.cxx.

  : ShipHit()
{
 
  flag = true;
 
  double pointX =  p->GetX();
  double pointY =  p->GetY();
  double pointZ =  p->GetZ();
  double pointT =  p->GetTime();
  double pointE =  p->GetEnergyLoss();
  int detID =  p->GetDetectorID();
 
  //fdigi = t0 + t; 
  fdigi = t0 ;  
  // SetDigi(SetTimeRes(fdigi));
  SetDetectorID(detID);
 
  TGeoNavigator* navigator = gGeoManager->GetCurrentNavigator();
  navigator->cd("cave/SplitCalDetector_1");
  TGeoVolume* caloVolume = navigator->GetCurrentVolume();
  // caloVolume->PrintNodes();
 
  std::string stripName = GetDetectorElementName(detID); // it also sets if strip gives x or y coordinate
 
  int isPrec, nL, nMx, nMy, nS;
  Decoder(detID, isPrec, nL, nMx, nMy, nS);
 
  SetIDs(isPrec, nL, nMx, nMy, nS);
 
  TGeoNode* strip = caloVolume->GetNode(stripName.c_str()); 
 
  const Double_t* stripCoordinatesLocal = strip->GetMatrix()->GetTranslation();
  Double_t stripCoordinatesMaster[3] = {0.,0.,0.};
  navigator->LocalToMaster(stripCoordinatesLocal, stripCoordinatesMaster);
 
  TGeoBBox* box = (TGeoBBox*)strip->GetVolume()->GetShape();
  double xHalfLength = box->GetDX();
  double yHalfLength = box->GetDY();
  double zHalfLength = box->GetDZ();
  
  TGeoNode* passiveLayer = caloVolume->GetNode("ECALfilter_200000"); // they are all the same
  TGeoBBox* boxPassive = (TGeoBBox*)passiveLayer->GetVolume()->GetShape();
  double zPassiveHalfLength = box->GetDZ();
 
  // std::cout<< "----------------------"<<std::endl;
  // std::cout<< "-- pointX = " << pointX << std::endl; 
  // std::cout<< "-- pointY = " << pointY << std::endl; 
  // std::cout<< "-- pointZ = " << pointZ << std::endl; 
  // std::cout<< "-- detID = " << detID << std::endl;
  // std::cout<< "-- stripName = " << stripName << std::endl;
  // std::cout<< "-- isPrec = " << isPrec << std::endl;
  // std::cout<< "-- nL = " << nL << std::endl;
  // std::cout<< "-- nMx = " << nMx << std::endl;
  // std::cout<< "-- nMy = " << nMy << std::endl;
  // std::cout<< "-- nS = " << nS << std::endl;    
  // std::cout<< "-- stripCoordinatesLocal[0] = " << stripCoordinatesLocal[0] << std::endl;
  // std::cout<< "-- stripCoordinatesLocal[1] = " << stripCoordinatesLocal[1] << std::endl;
  // std::cout<< "-- stripCoordinatesLocal[2] = " << stripCoordinatesLocal[2] << std::endl;
  // std::cout<< "-- stripCoordinatesMaster[0] = " << stripCoordinatesMaster[0] << std::endl;
  // std::cout<< "-- stripCoordinatesMaster[1] = " << stripCoordinatesMaster[1] << std::endl;
  // std::cout<< "-- stripCoordinatesMaster[2] = " << stripCoordinatesMaster[2] << std::endl;
 
 
  // TGeoNode* check = navigator->FindNode(pointX,pointY,pointZ);
 
  SetEnergy(pointE);
  if (isPrec==1) SetXYZ(pointX,pointY,stripCoordinatesMaster[2]);
  else  SetXYZ(stripCoordinatesMaster[0], stripCoordinatesMaster[1], stripCoordinatesMaster[2]);
  SetXYZErrors(xHalfLength,yHalfLength,2*(zHalfLength+zPassiveHalfLength));
  
 
}

◆ ~splitcalHit()

splitcalHit::~splitcalHit ( )

virtual

Destructor

Definition at line 205 of file splitcalHit.cxx.

205{ }

◆ splitcalHit() [4/4]

splitcalHit::splitcalHit ( const splitcalHit & point )

Copy constructor

Member Function Documentation

◆ AddClusterIndex()

void splitcalHit::AddClusterIndex ( int i )

inline

Definition at line 48 of file splitcalHit.h.

48{_vecClusterIndices.push_back(i);}

splitcalHit::_vecClusterIndices

std::vector< int > _vecClusterIndices

Definition splitcalHit.h:87

◆ AddEnergyWeight()

void splitcalHit::AddEnergyWeight ( double w )

inline

Definition at line 49 of file splitcalHit.h.

49{_vecEnergyWeights.push_back(w);}

splitcalHit::_vecEnergyWeights

std::vector< double > _vecEnergyWeights

Definition splitcalHit.h:86

◆ ClassDef()

splitcalHit::ClassDef	(	splitcalHit	,
		3
	)

private

◆ Decoder() [1/2]

void splitcalHit::Decoder	(	int &	id,
		int &	isPrecision,
		int &	nLayer,
		int &	nModuleX,
		int &	nMdouleY,
		int &	nStrip
	)

Definition at line 172 of file splitcalHit.cxx.

                                                                                                           {
 
  std::string encodedID = GetPaddedString(id);
  Decoder(encodedID, isPrecision, nLayer, nModuleX, nModuleY, nStrip);
 
}

◆ Decoder() [2/2]

void splitcalHit::Decoder	(	std::string &	encodedID,
		int &	isPrecision,
		int &	nLayer,
		int &	nModuleX,
		int &	nMdouleY,
		int &	nStrip
	)

Definition at line 151 of file splitcalHit.cxx.

                                                                                                                        {
 
  std::string subtring;
 
  subtring = encodedID.substr(0, 1);
  isPrecision = atoi(subtring.c_str());
 
  subtring = encodedID.substr(1,3);
  nLayer = atoi(subtring.c_str());
 
  subtring = encodedID.substr(4,1); 
  nModuleX = atoi(subtring.c_str());
 
  subtring = encodedID.substr(5,1); 
  nModuleY = atoi(subtring.c_str());
 
  subtring = encodedID.substr(6,3);
  nStrip = atoi(subtring.c_str());
 
}

◆ GetClusterIndices()

std::vector< int > splitcalHit::GetClusterIndices ( )

inline

Definition at line 69 of file splitcalHit.h.

69{return _vecClusterIndices;}

◆ GetDetectorElementName()

std::string splitcalHit::GetDetectorElementName ( int & id )

Definition at line 123 of file splitcalHit.cxx.

                                                    {
 
  std::string encodedID = GetPaddedString(id);
  //std::cout << "-- encodedID = " << encodedID <<std::endl;
  int isPrec, nL, nMx, nMy, nS;
  Decoder(encodedID, isPrec, nL, nMx, nMy, nS);
 
  std::string name;
  if (isPrec==1) {
    name = "ECALdet_gas_";
    SetIsX(true);
    SetIsY(true);
  } else if (nL%2==0) {
    name = "stripGivingY_"; 
    SetIsX(false);
    SetIsY(true);
  } else {
    name = "stripGivingX_";
    SetIsX(true);
    SetIsY(false);
  }
  name = name + std::to_string(id);
  // std::cout << "--GetDetectorElementName - name  = " << name <<std::endl;
 
  return name;
 
}

◆ GetEnergy()

double splitcalHit::GetEnergy ( )

inline

Definition at line 54 of file splitcalHit.h.

54{return _energy;}

splitcalHit::_energy

double _energy

Definition splitcalHit.h:82

◆ GetEnergyForCluster()

double splitcalHit::GetEnergyForCluster ( int i )

Definition at line 193 of file splitcalHit.cxx.

                                            {
 
  double unweightedEnergy = GetEnergy();
  double weight = GetEnergyWeightForIndex(i);
  double energy = unweightedEnergy*weight;
  return energy;
 
}

◆ GetEnergyWeightForIndex()

double splitcalHit::GetEnergyWeightForIndex ( int index )

Definition at line 180 of file splitcalHit.cxx.

                                                    {
 
  int iw = 0;
  for(size_t i=0; i<_vecClusterIndices.size(); i++) {
    if (_vecClusterIndices.at(i) == index) {
      iw = i;
      break; 
    }
  }
  return _vecEnergyWeights.at(iw); 
}

◆ GetEnergyWeights()

std::vector< double > splitcalHit::GetEnergyWeights ( )

inline

Definition at line 70 of file splitcalHit.h.

70{return _vecEnergyWeights;}

◆ GetIsPrecisionLayer()

int splitcalHit::GetIsPrecisionLayer ( )

inline

Definition at line 56 of file splitcalHit.h.

56{return _isPrecisionLayer;}

splitcalHit::_isPrecisionLayer

int _isPrecisionLayer

Definition splitcalHit.h:84

◆ GetLayerNumber()

int splitcalHit::GetLayerNumber ( )

inline

Definition at line 57 of file splitcalHit.h.

57{return _nLayer;}

splitcalHit::_nLayer

int _nLayer

Definition splitcalHit.h:84

◆ GetModuleXNumber()

int splitcalHit::GetModuleXNumber ( )

inline

Definition at line 58 of file splitcalHit.h.

58{return _nModuleX;}

splitcalHit::_nModuleX

int _nModuleX

Definition splitcalHit.h:84

◆ GetModuleYNumber()

int splitcalHit::GetModuleYNumber ( )

inline

Definition at line 59 of file splitcalHit.h.

59{return _nModuleY;}

splitcalHit::_nModuleY

int _nModuleY

Definition splitcalHit.h:84

◆ GetPaddedString()

std::string splitcalHit::GetPaddedString ( int & id )

Definition at line 112 of file splitcalHit.cxx.

                                             { 
 
  //zero padded string 
  int totalLength = 9; 
  std::string stringID = std::to_string(id);
  std::string encodedID = std::string(totalLength - stringID.length(), '0') + stringID;
 
 return encodedID;
 
}

◆ GetStripNumber()

int splitcalHit::GetStripNumber ( )

inline

Definition at line 60 of file splitcalHit.h.

60{return _nStrip;}

splitcalHit::_nStrip

int _nStrip

Definition splitcalHit.h:84

◆ GetTDC()

Float_t splitcalHit::GetTDC ( ) const

inline

Definition at line 29 of file splitcalHit.h.

29{return fdigi;}

◆ GetX()

double splitcalHit::GetX ( )

inline

Definition at line 51 of file splitcalHit.h.

51{return _x;}

splitcalHit::_x

double _x

Definition splitcalHit.h:81

◆ GetXError()

double splitcalHit::GetXError ( )

inline

Definition at line 64 of file splitcalHit.h.

64{return _xError;}

splitcalHit::_xError

double _xError

Definition splitcalHit.h:81

◆ GetY()

double splitcalHit::GetY ( )

inline

Definition at line 52 of file splitcalHit.h.

52{return _y;}

splitcalHit::_y

double _y

Definition splitcalHit.h:81

◆ GetYError()

double splitcalHit::GetYError ( )

inline

Definition at line 65 of file splitcalHit.h.

65{return _yError;}

splitcalHit::_yError

double _yError

Definition splitcalHit.h:81

◆ GetZ()

double splitcalHit::GetZ ( )

inline

Definition at line 53 of file splitcalHit.h.

53{return _z;}

splitcalHit::_z

double _z

Definition splitcalHit.h:81

◆ GetZError()

double splitcalHit::GetZError ( )

inline

Definition at line 66 of file splitcalHit.h.

66{return _zError;}

splitcalHit::_zError

double _zError

Definition splitcalHit.h:81

◆ IsShared()

bool splitcalHit::IsShared ( )

inline

Definition at line 71 of file splitcalHit.h.

71{return GetClusterIndices().size()>1; }

splitcalHit::GetClusterIndices

std::vector< int > GetClusterIndices()

Definition splitcalHit.h:69

◆ IsUsed()

int splitcalHit::IsUsed ( )

inline

Definition at line 63 of file splitcalHit.h.

63{return _isUsed;}

splitcalHit::_isUsed

int _isUsed

Definition splitcalHit.h:84

◆ isValid()

bool splitcalHit::isValid ( ) const

inline

Definition at line 31 of file splitcalHit.h.

31{return flag;}

◆ IsX()

bool splitcalHit::IsX ( )

inline

Definition at line 61 of file splitcalHit.h.

61{return _isX;}

splitcalHit::_isX

bool _isX

Definition splitcalHit.h:85

◆ IsY()

bool splitcalHit::IsY ( )

inline

Definition at line 62 of file splitcalHit.h.

62{return _isY;}

splitcalHit::_isY

bool _isY

Definition splitcalHit.h:85

◆ operator=()

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

◆ Print()

void splitcalHit::Print ( ) const

virtual

Output to screen

Definition at line 209 of file splitcalHit.cxx.

{
  //  cout << "-I- splitcalHit: splitcal hit " << " in detector " << fDetectorID << endl;
  //  cout << "  TDC " << fdigi << " ns" << endl;
  // std::cout<< "-I- splitcalHit: " <<std::endl;
  std::cout<< "------- " <<std::endl;
  std::cout<< "    (x,y,z) = " 
       << _x << " +- " << _xError << " ,  "
       << _y << " +- " << _yError << " ,  "
       << _z << " +- " << _zError <<std::endl;  
  std::cout<< "    isP, nL, nMx, nMy, nS = " 
       << _isPrecisionLayer << " , "  
       << _nLayer << " , "  
       << _nModuleX << " , " 
       << _nModuleY << " , " 
       << _nStrip << std::endl;
   std::cout<< "------- " <<std::endl;
}

◆ SetEnergy()

void splitcalHit::SetEnergy ( double & e )

inline

Definition at line 40 of file splitcalHit.h.

40{_energy = e;}

◆ SetIDs()

void splitcalHit::SetIDs	(	int &	isPrecision,
		int &	nLayer,
		int &	nModuleX,
		int &	nModuleY,
		int &	nStrip
	)

inline

Definition at line 39 of file splitcalHit.h.

39{_isPrecisionLayer = isPrecision; _nLayer = nLayer; _nModuleX = nModuleX; _nModuleY=nModuleY; _nStrip = nStrip;}

◆ setInvalid()

void splitcalHit::setInvalid ( )

inline

Definition at line 30 of file splitcalHit.h.

30{flag = false;}

◆ SetIsUsed()

void splitcalHit::SetIsUsed ( int u )

inline

Definition at line 44 of file splitcalHit.h.

44{_isUsed = u;}

◆ SetIsX()

void splitcalHit::SetIsX ( bool x )

inline

Definition at line 42 of file splitcalHit.h.

42{_isX = x;}

geometry_config.x

x

Definition geometry_config.py:180

◆ SetIsY()

void splitcalHit::SetIsY ( bool y )

inline

Definition at line 43 of file splitcalHit.h.

43{_isY = y;}

geometry_config.y

y

Definition geometry_config.py:179

◆ SetXYZ()

void splitcalHit::SetXYZ	(	double &	x,
		double &	y,
		double &	z
	)

inline

Definition at line 38 of file splitcalHit.h.

38{_x = x; _y = y; _z = z;}

charm-geometry_config.z

z

Definition charm-geometry_config.py:7

◆ SetXYZErrors()

void splitcalHit::SetXYZErrors	(	double	xError,
		double	yError,
		double	zError
	)

inline

Definition at line 45 of file splitcalHit.h.

45{_xError = xError; _yError = yError; _zError = zError;}

◆ UpdateEnergy()

void splitcalHit::UpdateEnergy ( double e )

inline

Definition at line 41 of file splitcalHit.h.

41{_energy = _energy+e;}

Member Data Documentation

◆ _energy

double splitcalHit::_energy

private

Definition at line 82 of file splitcalHit.h.

◆ _isPrecisionLayer

int splitcalHit::_isPrecisionLayer

private

Definition at line 84 of file splitcalHit.h.

◆ _isUsed

int splitcalHit::_isUsed

private

Definition at line 84 of file splitcalHit.h.

◆ _isX

bool splitcalHit::_isX

private

Definition at line 85 of file splitcalHit.h.

◆ _isY

bool splitcalHit::_isY

private

Definition at line 85 of file splitcalHit.h.

◆ _nLayer

int splitcalHit::_nLayer

private

Definition at line 84 of file splitcalHit.h.

◆ _nModuleX

int splitcalHit::_nModuleX

private

Definition at line 84 of file splitcalHit.h.

◆ _nModuleY

int splitcalHit::_nModuleY

private

Definition at line 84 of file splitcalHit.h.

◆ _nStrip

int splitcalHit::_nStrip

private

Definition at line 84 of file splitcalHit.h.

◆ _vecClusterIndices

std::vector<int > splitcalHit::_vecClusterIndices

private

Definition at line 87 of file splitcalHit.h.

◆ _vecEnergyWeights

std::vector<double > splitcalHit::_vecEnergyWeights

private

Definition at line 86 of file splitcalHit.h.

◆ _x

double splitcalHit::_x

private

Definition at line 81 of file splitcalHit.h.

◆ _xError

double splitcalHit::_xError

private

Definition at line 81 of file splitcalHit.h.

◆ _y

double splitcalHit::_y

private

Definition at line 81 of file splitcalHit.h.

◆ _yError

double splitcalHit::_yError

private

Definition at line 81 of file splitcalHit.h.

◆ _z

double splitcalHit::_z

private

Definition at line 81 of file splitcalHit.h.

◆ _zError

double splitcalHit::_zError

private

Definition at line 81 of file splitcalHit.h.

◆ flag

Float_t splitcalHit::flag

private

Definition at line 80 of file splitcalHit.h.

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

splitcal/splitcalHit.h
splitcal/splitcalHit.cxx

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ splitcalHit() [1/4]

◆ splitcalHit() [2/4]

◆ splitcalHit() [3/4]

◆ ~splitcalHit()

◆ splitcalHit() [4/4]

Member Function Documentation

◆ AddClusterIndex()

◆ AddEnergyWeight()

◆ ClassDef()

◆ Decoder() [1/2]

◆ Decoder() [2/2]

◆ GetClusterIndices()

◆ GetDetectorElementName()

◆ GetEnergy()

◆ GetEnergyForCluster()

◆ GetEnergyWeightForIndex()

◆ GetEnergyWeights()

◆ GetIsPrecisionLayer()

◆ GetLayerNumber()

◆ GetModuleXNumber()

◆ GetModuleYNumber()

◆ GetPaddedString()

◆ GetStripNumber()

◆ GetTDC()

◆ GetX()

◆ GetXError()

◆ GetY()

◆ GetYError()

◆ GetZ()

◆ GetZError()

◆ IsShared()

◆ IsUsed()

◆ isValid()

◆ IsX()

◆ IsY()

◆ operator=()

◆ Print()

◆ SetEnergy()

◆ SetIDs()

◆ setInvalid()

◆ SetIsUsed()

◆ SetIsX()

◆ SetIsY()

◆ SetXYZ()

◆ SetXYZErrors()

◆ UpdateEnergy()

Member Data Documentation

◆ _energy

◆ _isPrecisionLayer

◆ _isUsed

◆ _isX

◆ _isY

◆ _nLayer

◆ _nModuleX

◆ _nModuleY

◆ _nStrip

◆ _vecClusterIndices

◆ _vecEnergyWeights

◆ _x

◆ _xError

◆ _y

◆ _yError

◆ _z

◆ _zError

◆ flag