class THaOpticsHRS: public THaOptics


    THaOpticsHRS

  Derived class of THaOptics, which uses matrix expansion algorithm
  The core code is copied from JLab Hall A VDC vertex reconstruction



 Author: Jin Huang <mailto:jinhuang@jlab.org>    Feb 2008
    Modify History:

Function Members (Methods)

public:
THaOpticsHRS(const char* name, const char* desc, THaSpectrometer* pspec, TString strBeamDetectorName = "")
virtual~THaOpticsHRS()
voidTObject::AbstractMethod(const char* method) const
virtual voidTObject::AppendPad(Option_t* option = "")
virtual Int_tTHaOptics::ApplyCorrection(Double_t& Th, Double_t& Ph, Double_t& P, Double_t& PathLen)
virtual Int_tApplyOptics(Double_t trackX, Double_t trackY, Double_t trackTheta, Double_t trackPhi, TVector3 beamPos, TVector3 beamDir, TVector3& P, TVector3& Vertex, Double_t& TargetX, Double_t& TargetY, Double_t& TargetTheta, Double_t& TargetPhi, Double_t& PathLen)
virtual Int_tTHaAnalysisObject::Begin(THaRunBase* r = 0)
virtual voidTObject::Browse(TBrowser* b)
static TClass*Class()
virtual const char*TObject::ClassName() const
virtual voidTHaOptics::Clear(Option_t* opt = "")
virtual TObject*TNamed::Clone(const char* newname = "") const
virtual Int_tTNamed::Compare(const TObject* obj) const
virtual voidTNamed::Copy(TObject& named) const
TVector2THaOptics::CrossingLineLine2D(TVector2 o1, TVector2 d1, TVector2 o2, TVector2 d2)
TVector3THaOptics::CrossingSurfLine3D(TVector3 OSurf, TVector3 NSurf, TVector3 OLine, TVector3 DLine)
static Int_tTHaAnalysisObject::DefineVarsFromList(const void* list, THaAnalysisObject::EType type, THaAnalysisObject::EMode mode, const char* var_prefix, const TObject* obj, const char* prefix, const char* here)
virtual voidTObject::Delete(Option_t* option = "")
virtual Int_tTObject::DistancetoPrimitive(Int_t px, Int_t py)
virtual voidTObject::Draw(Option_t* option = "")
virtual voidTObject::DrawClass() const
virtual TObject*TObject::DrawClone(Option_t* option = "") const
virtual voidTObject::Dump() const
virtual Int_tTHaAnalysisObject::End(THaRunBase* r = 0)
virtual voidTObject::Error(const char* method, const char* msgfmt) const
virtual voidTObject::Execute(const char* method, const char* params, Int_t* error = 0)
virtual voidTObject::Execute(TMethod* method, TObjArray* params, Int_t* error = 0)
virtual voidTObject::ExecuteEvent(Int_t event, Int_t px, Int_t py)
virtual voidTObject::Fatal(const char* method, const char* msgfmt) const
virtual voidTNamed::FillBuffer(char*& buffer)
virtual TObject*TObject::FindObject(const char* name) const
virtual TObject*TObject::FindObject(const TObject* obj) const
static voidTHaAnalysisObject::GeoToSph(Double_t th_geo, Double_t ph_geo, Double_t& th_sph, Double_t& ph_sph)
Double_tTHaOptics::GetCollDist() const
const char*THaAnalysisObject::GetConfig() const
Double_tTHaOptics::GetCurrent() const
virtual const char*THaAnalysisObject::GetDBFileName() const
Int_tTHaAnalysisObject::GetDebug() const
virtual Option_t*TObject::GetDrawOption() const
static Long_tTObject::GetDtorOnly()
virtual const char*TObject::GetIconName() const
TDatimeTHaAnalysisObject::GetInitDate() const
virtual const char*TNamed::GetName() const
virtual char*TObject::GetObjectInfo(Int_t px, Int_t py) const
static Bool_tTObject::GetObjectStat()
virtual Option_t*TObject::GetOption() const
Double_tTHaOptics::GetPcentral() const
const char*THaAnalysisObject::GetPrefix() const
virtual const char*TNamed::GetTitle() const
virtual UInt_tTObject::GetUniqueID() const
virtual Bool_tTObject::HandleTimer(TTimer* timer)
virtual ULong_tTNamed::Hash() const
virtual voidTObject::Info(const char* method, const char* msgfmt) const
virtual Bool_tTObject::InheritsFrom(const char* classname) const
virtual Bool_tTObject::InheritsFrom(const TClass* cl) const
THaAnalysisObject::EStatusTHaAnalysisObject::Init()
virtual THaAnalysisObject::EStatusTHaAnalysisObject::Init(const TDatime& run_time)
virtual Int_tTHaAnalysisObject::InitOutput(THaOutput*)
virtual voidTObject::Inspect() const
static Bool_tTHaAnalysisObject::IntersectPlaneWithRay(const TVector3& xax, const TVector3& yax, const TVector3& org, const TVector3& ray_start, const TVector3& ray_vect, Double_t& length, TVector3& intersect)
voidTObject::InvertBit(UInt_t f)
virtual TClass*IsA() const
virtual Bool_tTObject::IsEqual(const TObject* obj) const
virtual Bool_tTObject::IsFolder() const
Bool_tTHaAnalysisObject::IsInit() const
Bool_tTHaAnalysisObject::IsOK() const
Bool_tTHaAnalysisObject::IsOKOut()
Bool_tTObject::IsOnHeap() const
virtual Bool_tTNamed::IsSortable() const
Bool_tTObject::IsZombie() const
static Int_tTHaAnalysisObject::LoadDB(FILE* file, const TDatime& date, const DBRequest* request, const char* prefix = "", Int_t search = 0)
static Int_tTHaAnalysisObject::LoadDB(FILE* file, const TDatime& date, const TagDef* tags, const char* prefix = "", Int_t search = 0)
static Int_tTHaAnalysisObject::LoadDBvalue(FILE* file, const TDatime& date, const char* tag, Double_t& value)
static Int_tTHaAnalysisObject::LoadDBvalue(FILE* file, const TDatime& date, const char* tag, Int_t& value)
static Int_tTHaAnalysisObject::LoadDBvalue(FILE* file, const TDatime& date, const char* tag, string& text)
static Int_tTHaAnalysisObject::LoadDBvalue(FILE* file, const TDatime& date, const char* tag, TString& text)
virtual voidTNamed::ls(Option_t* option = "") const
voidTObject::MayNotUse(const char* method) const
virtual Bool_tTObject::Notify()
static FILE*THaAnalysisObject::OpenFile(const char* name, const TDatime& date, const char* here = "OpenFile()", const char* filemode = "r", const int debug_flag = 1)
static voidTObject::operator delete(void* ptr)
static voidTObject::operator delete(void* ptr, void* vp)
static voidTObject::operator delete[](void* ptr)
static voidTObject::operator delete[](void* ptr, void* vp)
void*TObject::operator new(size_t sz)
void*TObject::operator new(size_t sz, void* vp)
void*TObject::operator new[](size_t sz)
void*TObject::operator new[](size_t sz, void* vp)
TNamed&TNamed::operator=(const TNamed& rhs)
virtual voidTObject::Paint(Option_t* option = "")
virtual voidTObject::Pop()
virtual voidTNamed::Print(Option_t* option = "") const
virtual Int_tTObject::Read(const char* name)
virtual voidTObject::RecursiveRemove(TObject* obj)
voidTObject::ResetBit(UInt_t f)
virtual voidTObject::SaveAs(const char* filename = "", Option_t* option = "") const
virtual voidTObject::SavePrimitive(basic_ostream<char,char_traits<char> >& out, Option_t* option = "")
static Int_tTHaAnalysisObject::SeekDBconfig(FILE* file, const char* tag, const char* label = "config", Bool_t end_on_tag = false)
static Int_tTHaAnalysisObject::SeekDBdate(FILE* file, const TDatime& date, Bool_t end_on_tag = false)
Bool_tTHaOptics::SetBeamName(const char*)
voidTObject::SetBit(UInt_t f)
voidTObject::SetBit(UInt_t f, Bool_t set)
voidTHaAnalysisObject::SetConfig(const char* label)
virtual voidTHaAnalysisObject::SetDebug(Int_t level)
virtual voidTObject::SetDrawOption(Option_t* option = "")
static voidTObject::SetDtorOnly(void* obj)
virtual voidTHaAnalysisObject::SetName(const char* name)
virtual voidTHaAnalysisObject::SetNameTitle(const char* name, const char* title)
static voidTObject::SetObjectStat(Bool_t stat)
virtual voidTNamed::SetTitle(const char* title = "")
virtual voidTObject::SetUniqueID(UInt_t uid)
virtual voidShowMembers(TMemberInspector& insp, char* parent)
virtual Int_tTNamed::Sizeof() const
static voidTHaAnalysisObject::SphToGeo(Double_t th_sph, Double_t ph_sph, Double_t& th_geo, Double_t& ph_geo)
THaAnalysisObject::EStatusTHaAnalysisObject::Status() const
virtual voidStreamer(TBuffer& b)
voidStreamerNVirtual(TBuffer& b)
virtual voidTObject::SysError(const char* method, const char* msgfmt) const
Bool_tTObject::TestBit(UInt_t f) const
Int_tTObject::TestBits(UInt_t f) const
virtual voidTObject::UseCurrentStyle()
static vector<std::string>THaAnalysisObject::vsplit(const string& s)
virtual voidTObject::Warning(const char* method, const char* msgfmt) const
virtual Int_tTObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0)
virtual Int_tTObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0) const
virtual Int_tTHaOptics::ZeroOrderCorrection(Double_t& Th, Double_t& Ph, Double_t& P, Double_t& PathLen)
protected:
voidCalcMatrix(const double x, vector<THaMatrixElement>& matrix)
Double_tCalcTarget2FPLen(const vector<THaMatrixElement>& matrix, const Double_t powers[][5])
Double_tCalcTargetVar(const vector<THaMatrixElement>& matrix, const double powers[][5])
virtual Int_tTHaAnalysisObject::DefineVariables(THaAnalysisObject::EMode mode = kDefine)
Int_tTHaAnalysisObject::DefineVarsFromList(const VarDef* list, THaAnalysisObject::EMode mode = kDefine, const char* var_prefix = "") const
Int_tTHaAnalysisObject::DefineVarsFromList(const RVarDef* list, THaAnalysisObject::EMode mode = kDefine, const char* var_prefix = "") const
Int_tTHaAnalysisObject::DefineVarsFromList(const void* list, THaAnalysisObject::EType type, THaAnalysisObject::EMode mode, const char* var_prefix = "") const
virtual voidTObject::DoError(int level, const char* location, const char* fmt, va_list va) const
Double_tDoPoly(const int n, const vector<double,allocator<double> >& a, const double x)
THaAnalysisObject*THaAnalysisObject::FindModule(const char* name, const char* classname, bool do_error = true)
virtual THaAnalysisObject*THaOptics::FindModuleNoCheck(const char* name, const char* classname)
static vector<std::string>THaAnalysisObject::GetDBFileList(const char* name, const TDatime& date, const char* here = "GetDBFileList()")
virtual const char*THaAnalysisObject::Here(const char*) const
virtual voidTHaOptics::MakePrefix()
voidTObject::MakeZombie()
virtual FILE*THaAnalysisObject::OpenFile(const TDatime& date)
virtual FILE*THaAnalysisObject::OpenRunDBFile(const TDatime& date)
Double_tPolyInv(const double x1, const double x2, const double xacc, const double y, const int norder, const vector<double,allocator<double> >& a)
static char*THaAnalysisObject::ReadComment(FILE* fp, char* buf, const int len)
virtual Int_tReadDatabase(const TDatime& date)
virtual Int_tTHaOptics::ReadRunDatabase(const TDatime& date)
virtual Int_tTHaAnalysisObject::RemoveVariables()

Data Members

public:
enum EFPMatrixElemTags { T000
Y000
P000
};
enum { kPORDER
};
enum THaAnalysisObject::EStatus { kOK
kNotinit
kInitError
kFileError
};
enum THaAnalysisObject::EType { kVarDef
kRVarDef
};
enum THaAnalysisObject::EMode { kDefine
kDelete
};
enum THaAnalysisObject::EProperties { kNeedsRunDB
};
enum TObject::EStatusBits { kCanDelete
kMustCleanup
kObjInCanvas
kIsReferenced
kHasUUID
kCannotPick
kNoContextMenu
kInvalidObject
};
enum TObject::[unnamed] { kIsOnHeap
kNotDeleted
kZombie
kBitMask
kSingleKey
kOverwrite
kWriteDelete
};
public:
static const Double_tTHaAnalysisObject::kBig= 1.e38; // default junk value
protected:
TRefTHaOptics::fBeam
TStringTHaAnalysisObject::fConfigConfiguration to use from database
Double_tTHaOptics::fCurrentCentral momentum (GeV)
vector<THaMatrixElement>fDMatrixElems
Int_tTHaAnalysisObject::fDebugDebug level
TStringTHaOptics::fDefBeamName
vector<THaMatrixElement>fFPMatrixElemsmatrix elements used in
TRotationTHaOptics::fHCS2TCS
TDatimeTHaAnalysisObject::fInitDateDate passed to Init
Bool_tTHaAnalysisObject::fIsInitFlag indicating that ReadDatabase done
Bool_tTHaAnalysisObject::fIsSetupFlag indicating that DefineVariables done.
vector<THaMatrixElement>fLMatrixElemsPath-length corrections (meters)
TStringTNamed::fNameobject identifier
Bool_tTHaAnalysisObject::fOKOutFlag indicating object-output prepared
vector<THaMatrixElement>fPMatrixElems
Double_tTHaOptics::fPShift
vector<THaMatrixElement>fPTAMatrixElemsinvolves abs(theta_fp)
Double_tTHaOptics::fPathLenShift
Double_tTHaOptics::fPhShift
char*THaAnalysisObject::fPrefixName prefix for global variables
UInt_tTHaAnalysisObject::fPropertiesProperties of this object (see EProperties)
TRefTHaOptics::fSpecpointer to BigBite
THaAnalysisObject::EStatusTHaAnalysisObject::fStatusInitialization status flag
TRotationTHaOptics::fTCS2HCS
TRotationTHaOptics::fTCS2TRCS
TVector3THaOptics::fTCSOffSetInHCS
vector<THaMatrixElement>fTMatrixElems
TRotationTHaOptics::fTRCS2TCS
TVector3THaOptics::fTRCSOffSetInTCS
Double_tTHaOptics::fThShift
TStringTNamed::fTitleobject title
vector<THaMatrixElement>fYMatrixElems
vector<THaMatrixElement>fYTAMatrixElemsinvolves abs(theta_fp)

Class Charts

Inheritance Chart:
TObject
←
TNamed
←
THaAnalysisObject
←
THaOptics
←
THaOpticsHRS

Function documentation

THaOpticsHRS(const char* name, const char* desc, THaSpectrometer* pspec, TString strBeamDetectorName = "")
 constructor function of THaOpticsHRS
~THaOpticsHRS()



{}






Int_t ReadDatabase(const TDatime& date)




void CalcMatrix(const double x, vector<THaMatrixElement>& matrix)


 calculates the values of the matrix elements for a given location
 by evaluating a polynomial in x of order it->order with
 coefficients given by it->poly





Double_t CalcTargetVar(const vector<THaMatrixElement>& matrix, const double powers[][5])


 calculates the value of a variable at the target
 the x-dependence is already in the matrix, so only 1-3 (or np) used





Double_t CalcTarget2FPLen(const vector<THaMatrixElement>& matrix, const Double_t powers[][5])


 calculates distance from the nominal target position (z=0)
 to the transport plane





Int_t ApplyOptics(Double_t trackX, Double_t trackY, Double_t trackTheta, Double_t trackPhi, TVector3 beamPos, TVector3 beamDir, TVector3& P, TVector3& Vertex, Double_t& TargetX, Double_t& TargetY, Double_t& TargetTheta, Double_t& TargetPhi, Double_t& PathLen)


apply optics : do the job





Double_t DoPoly(const int n, const vector<double,allocator<double> >& a, const double x)




Double_t PolyInv(const double x1, const double x2, const double xacc, const double y, const int norder, const vector<double,allocator<double> >& a)







Author:  Jin Huang     Feb 2008

Last update: Tue Jul  7 19:26:17 2009





This page has been automatically generated. If you have any comments or suggestions about the page layout send a mail to ROOT support, or contact the developers with any questions or problems regarding ROOT.