class THaScintillator: public THaNonTrackingDetector
THaScintillator Class for a generic scintillator (hodoscope) consisting of multiple paddles with phototubes on both ends.
Function Members (Methods)
public:
| THaScintillator(const char* name, const char* description = "", THaApparatus* a = NULL) | |
| virtual | ~THaScintillator() |
| void | TObject::AbstractMethod(const char* method) const |
| virtual void | TObject::AppendPad(Option_t* option = "") |
| virtual Int_t | ApplyCorrections() |
| virtual Int_t | THaAnalysisObject::Begin(THaRunBase* r = 0) |
| virtual void | TObject::Browse(TBrowser* b) |
| bool | THaSpectrometerDetector::CalcInterceptCoords(THaTrack* track, Double_t& x, Double_t& y) |
| bool | THaSpectrometerDetector::CalcPathLen(THaTrack* track, Double_t& t) |
| bool | THaSpectrometerDetector::CheckIntercept(THaTrack* track) |
| static TClass* | Class() |
| virtual const char* | TObject::ClassName() const |
| virtual void | Clear(Option_t* = "") |
| virtual TObject* | TNamed::Clone(const char* newname = "") const |
| virtual Int_t | CoarseProcess(TClonesArray& tracks) |
| virtual Int_t | TNamed::Compare(const TObject* obj) const |
| virtual void | TNamed::Copy(TObject& named) const |
| virtual Int_t | Decode(const THaEvData&) |
| static Int_t | THaAnalysisObject::DefineVarsFromList(const void* list, THaAnalysisObject::EType type, THaAnalysisObject::EMode mode, const char* var_prefix, const TObject* obj, const char* prefix, const char* here) |
| virtual void | TObject::Delete(Option_t* option = "")MENU |
| virtual Int_t | TObject::DistancetoPrimitive(Int_t px, Int_t py) |
| virtual void | TObject::Draw(Option_t* option = "") |
| virtual void | TObject::DrawClass() constMENU |
| virtual TObject* | TObject::DrawClone(Option_t* option = "") constMENU |
| virtual void | TObject::Dump() constMENU |
| virtual Int_t | THaAnalysisObject::End(THaRunBase* r = 0) |
| virtual void | TObject::Error(const char* method, const char* msgfmt) const |
| virtual void | TObject::Execute(const char* method, const char* params, Int_t* error = 0) |
| virtual void | TObject::Execute(TMethod* method, TObjArray* params, Int_t* error = 0) |
| virtual void | TObject::ExecuteEvent(Int_t event, Int_t px, Int_t py) |
| virtual void | TObject::Fatal(const char* method, const char* msgfmt) const |
| virtual void | TNamed::FillBuffer(char*& buffer) |
| Int_t | THaDetectorBase::FillDetMap(const vector<Int_t>& values, UInt_t flags = 0, const char* here = "FillDetMap") |
| virtual TObject* | TObject::FindObject(const char* name) const |
| virtual TObject* | TObject::FindObject(const TObject* obj) const |
| virtual Int_t | FineProcess(TClonesArray& tracks) |
| static void | THaAnalysisObject::GeoToSph(Double_t th_geo, Double_t ph_geo, Double_t& th_sph, Double_t& ph_sph) |
| const Double_t* | GetAmplitudes() const |
| THaApparatus* | THaDetector::GetApparatus() const |
| const char* | THaAnalysisObject::GetClassName() const |
| const char* | THaAnalysisObject::GetConfig() const |
| virtual const char* | THaAnalysisObject::GetDBFileName() const |
| Int_t | THaAnalysisObject::GetDebug() const |
| THaDetMap* | THaDetectorBase::GetDetMap() const |
| virtual Option_t* | TObject::GetDrawOption() const |
| static Long_t | TObject::GetDtorOnly() |
| virtual const char* | TObject::GetIconName() const |
| TDatime | THaAnalysisObject::GetInitDate() const |
| virtual const char* | TNamed::GetName() const |
| Int_t | THaDetectorBase::GetNelem() const |
| Int_t | GetNHits() const |
| Int_t | THaNonTrackingDetector::GetNTracks() const |
| static TString& | THaAnalysisObject::GetObjArrayString(const TObjArray* params, Int_t pos) |
| virtual char* | TObject::GetObjectInfo(Int_t px, Int_t py) const |
| static Bool_t | TObject::GetObjectStat() |
| virtual Option_t* | TObject::GetOption() const |
| const TVector3& | THaDetectorBase::GetOrigin() const |
| const char* | THaAnalysisObject::GetPrefix() const |
| const Double_t* | THaDetectorBase::GetSize() const |
| const Double_t* | GetTimes() const |
| virtual const char* | TNamed::GetTitle() const |
| const TClonesArray* | THaNonTrackingDetector::GetTrackHits() const |
| const Double_t* | GetTuncer() const |
| virtual UInt_t | TObject::GetUniqueID() const |
| const TVector3& | THaDetectorBase::GetXax() const |
| Double_t | THaDetectorBase::GetXSize() const |
| const Double_t* | GetYampl() const |
| const TVector3& | THaDetectorBase::GetYax() const |
| Double_t | THaDetectorBase::GetYSize() const |
| const Double_t* | GetYtime() const |
| const TVector3& | THaDetectorBase::GetZax() const |
| Double_t | THaDetectorBase::GetZSize() const |
| virtual Bool_t | TObject::HandleTimer(TTimer* timer) |
| virtual ULong_t | TNamed::Hash() const |
| virtual void | TObject::Info(const char* method, const char* msgfmt) const |
| virtual Bool_t | TObject::InheritsFrom(const char* classname) const |
| virtual Bool_t | TObject::InheritsFrom(const TClass* cl) const |
| virtual THaAnalysisObject::EStatus | Init(const TDatime& run_time) |
| virtual Int_t | THaAnalysisObject::InitOutput(THaOutput*) |
| virtual void | TObject::Inspect() constMENU |
| static Bool_t | THaAnalysisObject::IntersectPlaneWithRay(const TVector3& xax, const TVector3& yax, const TVector3& org, const TVector3& ray_start, const TVector3& ray_vect, Double_t& length, TVector3& intersect) |
| void | TObject::InvertBit(UInt_t f) |
| virtual TClass* | IsA() const |
| virtual Bool_t | TObject::IsEqual(const TObject* obj) const |
| virtual Bool_t | TObject::IsFolder() const |
| Bool_t | THaDetectorBase::IsInActiveArea(Double_t x, Double_t y) const |
| Bool_t | THaAnalysisObject::IsInit() const |
| Bool_t | THaAnalysisObject::IsOK() const |
| Bool_t | THaAnalysisObject::IsOKOut() |
| Bool_t | TObject::IsOnHeap() const |
| virtual Bool_t | THaNonTrackingDetector::IsPid() |
| virtual Bool_t | TNamed::IsSortable() const |
| static bool | THaAnalysisObject::IsTag(const char* buf) |
| virtual Bool_t | THaNonTrackingDetector::IsTracking() |
| Bool_t | TObject::IsZombie() const |
| static Int_t | THaAnalysisObject::LoadDB(FILE* file, const TDatime& date, const DBRequest* request, const char* prefix, Int_t search = 0, const char* here = "THaAnalysisObject::LoadDB") |
| static Int_t | THaAnalysisObject::LoadDBvalue(FILE* file, const TDatime& date, const char* tag, Double_t& value) |
| static Int_t | THaAnalysisObject::LoadDBvalue(FILE* file, const TDatime& date, const char* tag, Int_t& value) |
| static Int_t | THaAnalysisObject::LoadDBvalue(FILE* file, const TDatime& date, const char* tag, string& text) |
| static Int_t | THaAnalysisObject::LoadDBvalue(FILE* file, const TDatime& date, const char* tag, TString& text) |
| virtual void | TNamed::ls(Option_t* option = "") const |
| void | TObject::MayNotUse(const char* method) const |
| virtual Bool_t | TObject::Notify() |
| void | TObject::Obsolete(const char* method, const char* asOfVers, const char* removedFromVers) const |
| static FILE* | THaAnalysisObject::OpenFile(const char* name, const TDatime& date, const char* here = "OpenFile()", const char* filemode = "r", const int debug_flag = 1) |
| static void | TObject::operator delete(void* ptr) |
| static void | TObject::operator delete(void* ptr, void* vp) |
| static void | TObject::operator delete[](void* ptr) |
| static void | TObject::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 void | TObject::Paint(Option_t* option = "") |
| virtual void | TObject::Pop() |
| virtual void | THaAnalysisObject::Print(Option_t* opt = "") const |
| void | THaDetectorBase::PrintDetMap(Option_t* opt = "") const |
| static void | THaAnalysisObject::PrintObjects(Option_t* opt = "") |
| virtual Int_t | TObject::Read(const char* name) |
| static Int_t | THaAnalysisObject::ReadDBline(FILE* fp, char* buf, size_t bufsiz, string& line) |
| virtual void | TObject::RecursiveRemove(TObject* obj) |
| void | TObject::ResetBit(UInt_t f) |
| virtual void | TObject::SaveAs(const char* filename = "", Option_t* option = "") constMENU |
| virtual void | TObject::SavePrimitive(ostream& out, Option_t* option = "") |
| static Int_t | THaAnalysisObject::SeekDBconfig(FILE* file, const char* tag, const char* label = "config", Bool_t end_on_tag = false) |
| static Int_t | THaAnalysisObject::SeekDBdate(FILE* file, const TDatime& date, Bool_t end_on_tag = false) |
| virtual void | THaDetector::SetApparatus(THaApparatus*) |
| void | TObject::SetBit(UInt_t f) |
| void | TObject::SetBit(UInt_t f, Bool_t set) |
| void | THaAnalysisObject::SetConfig(const char* label) |
| virtual void | THaAnalysisObject::SetDebug(Int_t level) |
| virtual void | TObject::SetDrawOption(Option_t* option = "")MENU |
| static void | TObject::SetDtorOnly(void* obj) |
| virtual void | THaAnalysisObject::SetName(const char* name) |
| virtual void | THaAnalysisObject::SetNameTitle(const char* name, const char* title) |
| static void | TObject::SetObjectStat(Bool_t stat) |
| virtual void | TNamed::SetTitle(const char* title = "")MENU |
| virtual void | TObject::SetUniqueID(UInt_t uid) |
| virtual void | ShowMembers(TMemberInspector&) |
| virtual Int_t | TNamed::Sizeof() const |
| static void | THaAnalysisObject::SphToGeo(Double_t th_sph, Double_t ph_sph, Double_t& th_geo, Double_t& ph_geo) |
| THaAnalysisObject::EStatus | THaAnalysisObject::Status() const |
| virtual void | Streamer(TBuffer&) |
| void | StreamerNVirtual(TBuffer& ClassDef_StreamerNVirtual_b) |
| virtual void | TObject::SysError(const char* method, const char* msgfmt) const |
| Bool_t | TObject::TestBit(UInt_t f) const |
| Int_t | TObject::TestBits(UInt_t f) const |
| virtual void | TObject::UseCurrentStyle() |
| static vector<std::string> | THaAnalysisObject::vsplit(const string& s) |
| virtual void | TObject::Warning(const char* method, const char* msgfmt) const |
| virtual Int_t | TObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0) |
| virtual Int_t | TObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0) const |
protected:
| bool | THaSpectrometerDetector::CalcTrackIntercept(THaTrack* track, Double_t& t, Double_t& ycross, Double_t& xcross) |
| Int_t | THaNonTrackingDetector::CalcTrackProj(TClonesArray& tracks) |
| virtual const char* | THaAnalysisObject::ClassNameHere(const char*) const |
| virtual void | THaDetectorBase::DefineAxes(Double_t rotation_angle) |
| virtual Int_t | DefineVariables(THaAnalysisObject::EMode mode = kDefine) |
| Int_t | THaAnalysisObject::DefineVarsFromList(const VarDef* list, THaAnalysisObject::EMode mode = kDefine, const char* var_prefix = "") const |
| Int_t | THaAnalysisObject::DefineVarsFromList(const RVarDef* list, THaAnalysisObject::EMode mode = kDefine, const char* var_prefix = "") const |
| Int_t | THaAnalysisObject::DefineVarsFromList(const void* list, THaAnalysisObject::EType type, THaAnalysisObject::EMode mode, const char* var_prefix = "") const |
| void | DeleteArrays() |
| virtual void | THaAnalysisObject::DoError(int level, const char* location, const char* fmt, va_list va) const |
| THaAnalysisObject* | THaAnalysisObject::FindModule(const char* name, const char* classname, bool do_error = true) |
| static vector<std::string> | THaAnalysisObject::GetDBFileList(const char* name, const TDatime& date, const char* here = "GetDBFileList()") |
| virtual const char* | THaAnalysisObject::Here(const char*) const |
| Int_t | THaAnalysisObject::LoadDB(FILE* f, const TDatime& date, const DBRequest* req, Int_t search = 0) |
| virtual void | THaDetector::MakePrefix() |
| void | TObject::MakeZombie() |
| virtual FILE* | THaAnalysisObject::OpenFile(const TDatime& date) |
| virtual FILE* | THaAnalysisObject::OpenRunDBFile(const TDatime& date) |
| static char* | THaAnalysisObject::ReadComment(FILE* fp, char* buf, const int len) |
| virtual Int_t | ReadDatabase(const TDatime& date) |
| virtual Int_t | THaDetectorBase::ReadGeometry(FILE* file, const TDatime& date, Bool_t required = kFALSE) |
| virtual Int_t | THaAnalysisObject::ReadRunDatabase(const TDatime& date) |
| virtual Int_t | THaAnalysisObject::RemoveVariables() |
| virtual Double_t | TimeWalkCorrection(const Int_t& paddle, const THaScintillator::ESide side) |
Data Members
public:
| enum ESide { | kLeft | |
| kRight | ||
| }; | ||
| enum THaAnalysisObject::EStatus { | kOK | |
| kNotinit | ||
| kInitError | ||
| kFileError | ||
| }; | ||
| enum THaAnalysisObject::EType { | kVarDef | |
| kRVarDef | ||
| }; | ||
| enum THaAnalysisObject::EMode { | kDefine | |
| kDelete | ||
| }; | ||
| enum THaAnalysisObject::EProperties { | kNeedsRunDB | |
| kConfigOverride | ||
| }; | ||
| 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_t | THaAnalysisObject::kBig | = 1.e38; // default junk value |
protected:
| static const Int_t | NDEST | |
| Double_t | fAdcMIP | nominal ADC above pedestal for MIP |
| Double_t* | fAmpl | [fNelem] overall amplitude for the paddle |
| Double_t | fAttenuation | in m^-1: attenuation length of material |
| Double_t | fCn | speed of light in material (meters/second) |
| TString | THaAnalysisObject::fConfig | Configuration to use from database |
| THaScintillator::DataDest | fDataDest[2] | Lookup table for decoder |
| Int_t | THaAnalysisObject::fDebug | Debug level |
| THaDetMap* | THaDetectorBase::fDetMap | Hardware channel map for this detector |
| Int_t* | fHitPad | [fNhit] list of paddles with complete TDC hits |
| TDatime | THaAnalysisObject::fInitDate | Date passed to Init |
| Bool_t | THaAnalysisObject::fIsInit | Flag indicating that ReadDatabase done |
| Bool_t | THaAnalysisObject::fIsSetup | Flag indicating that DefineVariables done. |
| Double_t* | fLA | [fNelem] Array of Left paddles ADC amplitudes |
| Int_t | fLANhit | Number of Left paddles ADC amplitudes |
| Double_t* | fLA_c | [fNelem] Array of Left paddles corrected ADC ampl-s |
| Double_t* | fLA_p | [fNelem] Array of Left paddles ADC minus ped values |
| Double_t* | fLGain | [fNelem] ADC gains for left paddles |
| Double_t* | fLOff | [fNelem] TDC offsets for left paddles |
| Double_t* | fLPed | [fNelem] ADC pedestals for left paddles |
| Double_t* | fLT | [fNelem] Array of Left paddles TDC times (channels) |
| Int_t | fLTNhit | Number of Left paddles TDC times |
| Double_t* | fLT_c | [fNelem] Array of Left PMT corrected TDC times (s) |
| Int_t | fNTWalkPar | number of timewalk correction parameters |
| TString | TNamed::fName | object identifier |
| Int_t | THaDetectorBase::fNelem | Number of detector elements (paddles, mirrors) |
| Int_t | fNhit | Number of paddles with complete TDC hits (l&r) |
| Bool_t | THaAnalysisObject::fOKOut | Flag indicating object-output prepared |
| TVector3 | THaDetectorBase::fOrigin | Center position of detector (m) |
| char* | THaAnalysisObject::fPrefix | Name prefix for global variables |
| UInt_t | THaAnalysisObject::fProperties | Properties of this object (see EProperties) |
| Double_t* | fRA | [fNelem] Array of Right paddles ADC amplitudes |
| Int_t | fRANhit | Number of Right paddles ADC amplitudes |
| Double_t* | fRA_c | [fNelem] Array of Right paddles corrected ADC ampl-s |
| Double_t* | fRA_p | [fNelem] Array of Right paddles ADC minus ped values |
| Double_t* | fRGain | [fNelem] ADC gains for right paddles |
| Double_t* | fROff | [fNelem] TDC offsets for right paddles |
| Double_t* | fRPed | [fNelem] ADC pedestals for right paddles |
| Double_t* | fRT | [fNelem] Array of Right paddles TDC times (channels) |
| Int_t | fRTNhit | Number of Right paddles TDC times |
| Double_t* | fRT_c | [fNelem] Array of Right PMT corrected TDC times (s) |
| Double_t | fResolution | average time resolution per PMT (s) |
| Double_t | THaDetectorBase::fSize[3] | Detector size in x,y,z (m) - x,y are half-widths |
| THaAnalysisObject::EStatus | THaAnalysisObject::fStatus | Initialization status flag |
| Double_t* | fTWalkPar | [fNTWalkPar] time walk correction parameters |
| Double_t | fTdc2T | linear conversion between TDC and time (s/ch) |
| Double_t* | fTime | [fNelem] corrected time for the paddle (s) |
| TString | TNamed::fTitle | object title |
| TClonesArray* | THaNonTrackingDetector::fTrackProj | projection of track(s) onto detector plane |
| Double_t* | fTrigOff | [fNelem] Induced offset of trigger time from |
| TVector3 | THaDetectorBase::fXax | X axis of the detector plane |
| Double_t* | fYa | [fNelem] y-position of hit in paddle from ADC (m) |
| TVector3 | THaDetectorBase::fYax | Y axis of the detector plane |
| Double_t* | fYt | [fNelem] y-position of hit in paddle from TDC (m) |
| TVector3 | THaDetectorBase::fZax | Normal to the detector plane |
| Double_t* | fdTime | [fNelem] uncertainty in time (s) |
Class Charts
Function documentation
THaAnalysisObject::EStatus Init(const TDatime& run_time)
Extra initialization for scintillators: set up DataDest map
Int_t ReadDatabase(const TDatime& date)
Read this detector's parameters from the database file 'fi'. This function is called by THaDetectorBase::Init() once at the beginning of the analysis. 'date' contains the date/time of the run being analyzed.
Int_t DefineVariables(THaAnalysisObject::EMode mode = kDefine)
Initialize global variables and lookup table for decoder
Int_t Decode(const THaEvData& )
Decode scintillator data, correct TDC times and ADC amplitudes, and copy the data to the local data members. This implementation makes the following assumptions about the detector map: - The first half of the map entries corresponds to ADCs, the second half, to TDCs. - The first fNelem detector channels correspond to the PMTs on the right hand side, the next fNelem channels, to the left hand side. (Thus channel numbering for each module must be consecutive.)
Int_t ApplyCorrections()
Apply the ADC/TDC corrections to get the 'REAL' relevant TDC and ADC values. No tracking needs to have been done yet. Permits the dividing up of the decoding step (events could come from a different source) to the applying of corrections. For ease when trying to optimize calibrations
Double_t TimeWalkCorrection(const Int_t& paddle, const THaScintillator::ESide side)
Calculate the time-walk correction. The timewalk might be dependent upon the specific PMT, so information about exactly which PMT fired is required.
Int_t CoarseProcess(TClonesArray& tracks)
Scintillator coarse processing: - Apply timewalk corrections - Calculate rough transverse (y) position and energy deposition for hits for which PMTs on both ends of the paddle fired - Calculate rough track crossing points
Int_t FineProcess(TClonesArray& tracks)
Scintillator fine processing: - Reconstruct coordinates of track cross point with scintillator plane - For each crossing track, determine position residual and paddle number The position residuals are calculated along the x-coordinate (dispersive direction) only. They are useful to determine whether a track has a matching hit ( if abs(dx) <= 0.5*paddle x-width ) and, if so, how close to the edge of the paddle the track crossed. This assumes scintillator paddles oriented along the transverse (non-dispersive, y) direction.