DSTUtil::Target Namespace Reference

Functions
std::string	Name (const double &tgta)
int	Index (const double &tgta)
bool	IsInside (const MIPPVertexSummary *vtx, const MIPPRunSummary *rs)
bool	IsInsideLH2 (const MIPPVertexSummary *vtx, const MIPPRunSummary *rs)
bool	IsInside (const float *x, const MIPPRunSummary *rs)
double	LH2TrgtLen (double &rad)
bool	IsInsideLH2 (const float *x, const MIPPRunSummary *rs)

---

Function Documentation

std::string DSTUtil::Target::Name

(

const double &

tgta

)

Definition at line 1489 of file DSTUtil.cxx.

References A, name, RITargetAZ::ntgt, and RITargetAZ::tgt.

Referenced by MXMLdbTableEntryBuilder::Build(), MXMLdbTableBuilder::Build(), DCCalib::Init(), and main().

{
  for (int i=0 ; i != RITargetAZ::ntgt; ++i) {
    if (fabs(tgta - RITargetAZ::tgt[i].A) < 0.005)
      return RITargetAZ::tgt[i].name;
  }
  cout << "DSTUtil::Target::Name: Error: A=" 
       << tgta << " not found." << endl;
  return "Unknown";
}

int DSTUtil::Target::Index

(

const double &

tgta

)

Definition at line 1502 of file DSTUtil.cxx.

References A, RITargetAZ::ntgt, and RITargetAZ::tgt.

Referenced by main().

{
  for (int i=0 ; i!=RITargetAZ::ntgt; ++i) {
    if (fabs(tgta - RITargetAZ::tgt[i].A) < 0.005)
      return i;
  }
  cout << "DSTUtil::Target::Index: Error: A=" 
       << tgta << " not found." << endl;
  return -1;
}

bool DSTUtil::Target::IsInside	(	const MIPPVertexSummary *	vtx,
		const MIPPRunSummary *	rs
	)

Definition at line 1515 of file DSTUtil.cxx.

References MIPPVertexSummary::x.

{
  return DSTUtil::Target::IsInside(vtx->x, rs);
}

bool DSTUtil::Target::IsInsideLH2	(	const MIPPVertexSummary *	vtx,
		const MIPPRunSummary *	rs
	)

Definition at line 1522 of file DSTUtil.cxx.

References MIPPVertexSummary::x.

Referenced by main().

{
  return DSTUtil::Target::IsInsideLH2(vtx->x, rs);
}

bool DSTUtil::Target::IsInside	(	const float *	x,
		const MIPPRunSummary *	rs
	)

Definition at line 1530 of file DSTUtil.cxx.

References MIPPRunSummary::tgtPos, and MIPPRunSummary::tgtSize.

{
  //Naive (broken) implementation, NEEDS TO IMPROVE, TODO
  const float* tp = rs->tgtPos;
  const float* ts = rs->tgtSize;
  for (int dim = 0; dim !=3; ++dim) {
    if (x[dim] < tp[dim]-ts[dim] || x[dim] > tp[dim]+ts[dim]) 
      return false;
  }
  return true;
}

double DSTUtil::Target::LH2TrgtLen

(

double &

rad

)

Definition at line 1544 of file DSTUtil.cxx.

References par.

Referenced by IsInsideLH2().

                                           {
  //returns the half length in z for a given radius of the active lh2 shrunken volume.

  double  par[5] = {7.0980 , 0.48297  , -1.8404  ,1.7047 ,  -0.54161};
  double hlen = par[0] + par[1]*rad + par[2]*(rad*rad) + 
    par[3]*(rad*rad*rad) + par[4]*(rad*rad*rad*rad);
  return hlen;
}

bool DSTUtil::Target::IsInsideLH2	(	const float *	x,
		const MIPPRunSummary *	rs
	)

Definition at line 1553 of file DSTUtil.cxx.

References LH2TrgtLen(), MIPPRunSummary::tgta, and MIPPRunSummary::tgtPos.

{
  const float* tp = rs->tgtPos;
  //const float* ts = rs->tgtSize;
  const float  ta = rs->tgta;
  bool in = false;
  if(ta<2){
    //in lh2
    double dx = x[0] - tp[0];
    double dy = x[1] - tp[1];
    double dz = x[2] - tp[2];
    double rad = sqrt(dx*dx + dy*dy);
    double tzl = LH2TrgtLen(rad);
    //This is the shrunken LH2 target radius
    if(rad<1.8962 && fabs(dz)<tzl)in = true;

  }else{
    cout<<" Target is not LH2 "<<" atomic weight of target = "<<ta<<endl;
  }
  return in;
}

---