- Aug 26, 2018
-
-
Carlos Yero authored
1) THcDCWire.h a) added sigma to THcDCWire constructor which fills fSigmaWire b) added GetSimga method c) added SetSigma method 2) THcDCHit.h a) Added GetWireSigma method 3) THcDriftChamberPlane.h a) Added flag fUsingSigmaPerWire b) Added array fSigmaWire 4) THcDriftChamberPlane.cxx a) In ReadDatabase, set fUsingSigmaPerWire=0 by default with reading of parameter h_using_sigma_per_wire optional for setting fUsingSigmaPerWire. b) If fUsingSigmaPerWire=0 then sets all wires in array fSigmaWire to fSigma ( which is the simga per plane). c) If fUsingSigmaPerWire=1 then reads in the sigma per wire from a parameter file. d) Used new THcDCWire constructor to set the sigma for wire to fSigmaWire[nwire] e) Also change code so that the tzero offset is set in the THcDCWire constructor f) In ProcessHits , use wire->GetTOffSet instead of the fTzeroWire array when setting the time. 5) THcDC.cxx a) In TrackFit, replace fSigma with hit->GetWireSigma()
-
- Jul 17, 2018
-
-
Stephen A. Wood authored
-
- Apr 28, 2018
-
-
hallc-online authored
THcSpacePoint 1) added fSetStubFlag which is intialized to kFALSE in constructor and set to kTRUE when SetStub is called. 2) added method GetSetStubFlag which returns fSetStubFlag. THcDriftchamber 1) In ProcessHits increased fHits reserve to 40 2) In FindSpacePoints use fSpacePoints->Delete(); 3) In LeftRight ,only SetStub if the Stub is fit. THcDC 1) In LinkStubs a) only link stubs if number of total spacepoints <10 b) only link stubs if both spacepoints have SetStubFlag true THcDCHit.h 1) In constructor initalize fLR to 0 THcDCTrack.cxx 1) AddSpacePoint check if fnSP < 10
-
- Apr 16, 2018
-
-
hallc-online authored
THcDriftChamberPlane 1) Remove subtraction of hodoscope starttime from ProcessHits which is called by THcDC::Decode. The hodoscope start time is determined in hodoscope Decode which is called after THcDC::Decode. So drift chamber was using the starttime from the previous event. 2) Create method SubtractStartTime which is called by THcDC:CoarseTrack THcDriftChamber 1) Move fChambers[ic]->ProcessHits() to CoarseTrack 2) Before looping through fChambers[ic]->ProcessHits() in CoarseTrack loop through all planes and call fPlanes[ip]->SubtractStartTime()
-
- Mar 17, 2018
-
-
hallc-online authored
Modify THcDCTrack::TrackFit to calculate residuals for each plane without the plane involved in the track. Add tree variable residualExclPlane
-
- Mar 15, 2018
-
-
hallc-online authored
-
hallc-online authored
Modified THcSpacePoint.h to have public member that keeps track of the spacepoint ID number for each chamber. Used in THcDC::LinkStubs Modified THcDC.cxx and h 1) Add tree variables sp1_id and sp2_id which are the integer ID in the spacepoints in each chamber used in the golden track. The ID refers to index in array of spacepoints created in THcDriftChamber.cxx 2) Modify LinkStub to fill sp1_id and sp2_id for each track 3) For golden track fill sp1_ID_best and sp2_ID_best
-
- Feb 19, 2018
-
-
Stephen A. Wood authored
So that the counts can be reported in report files.
-
Stephen A. Wood authored
First hit above this cut is taken as the reference time. Cuts are per detector and read as an optional parameter in each detector class. TDCs and Flash ADC times have separate cuts. Implement a second method of selecting best reference time. If InitHitList is passed a possitive reference time cut, use the first reference time above that cut. If none of the hits pass the cut, then no reference time is set If InitHitList is passed a negative cut, the cut is taken to be the absolute value of the cut, but a reference time is garuanteed to be set as long as there is at least one reference time hit. If none of the reference time hits pass the cut, then the last hit is taken to be the reference time. (Should be the largest.) Allow reference time cuts for all the spectrometer detector classes Hodoscopes, drift chambers, Aerogels, Gas Cherenkovs, Calorimeters
-
- Feb 15, 2018
-
-
Stephen A. Wood authored
Add End method to each detector that calls MissReport which prints the number of events with missing TDC and ADC reference times.
-
- Dec 06, 2017
-
-
Stephen A. Wood authored
Each detector object will check the global variable X.present each event. If X.present is false, it indicates that the parent spectrometer object determineted that this spectrometer is not in this event, and therefor reference times are not expected. If X.present is false, then the reference time warning messages will be suppressed.
-
- Dec 02, 2017
-
-
Mark Jones authored
Add leaf dc.chisq which is the chisquared/dof of the golden track.
-
- Nov 03, 2017
-
-
Mark Jones authored
Set fHMSStyleChambers =1 according to dc_version parameter. If no dc_version parameter loaded defaults to old HMS style chambers
-
- Nov 02, 2017
-
-
Holly Szumila authored
-
- Oct 09, 2017
-
-
Holly Szumila authored
1) Add X,Y offsets parameter for each plane, instead of X,Y offset per chamber 2) ModifyTHcDriftChamber::CorrectHitTimes to work with SHMS chambers 3) Add GetReadoutSide method to THcDCHit class.
-
- Aug 31, 2017
-
-
Stephen A. Wood authored
-
- Jul 28, 2017
-
-
Eric Pooser authored
-
- Jul 21, 2017
-
-
Stephen A. Wood authored
-
- May 01, 2017
-
-
hallc-online authored
Main purpose is to calculate the DC efficiency per wire using the track information of the golden track. THcDC.h ----------- 1) add method EfficiciencyPerWire 2) Add variables fWire_hit_did and fWire_hit_should THcDC.cxx ------------ 1) fWire_hit_did and fWire_hit_should are fixed length arrays the size of NPlanes that in set in method Init 2) Add variables wireHitDid and wireHitShould to the tree 3) In method SetFocalPlaneBestTrack add call to EfficiencyPerWire method 4) added EfficiencyPerWire method which loops over the hits and compares the wire number for the hit to the the track position at the plane (convert back to wire number). Then fills fWire_hit_did with the wire number if equal to the wire position of the hit equals the track. Also does loop over all planes and fills fWire_hit_should with the wire number of the track in each plane. THcDriftChamberPlane.h and cxx --------------------- 1) Added method CalcWireFromPos
-
- Apr 26, 2017
-
-
hallc-online authored
The purpose is to add a call in THcHallCSpectrometer:TrackCalc to fill variables in THcDC for the golden track. THcHallCSpectrometer.h ---------------------- 1) Add object THcDC* fDC 2) Eliminate object fShower since it was not being used. THcHallCSpectrometer.cxx ---------------------- 1) in ReadDatabase method cast fDC in detector "dc" is defined 2) in FindVertices method add call to fDC->SetFocalPlaneBestTrack(hit_gold_track) where hit_gold_track is the HaTracks array index of the golden track THcDC.h ------- 1) Add method SetFocalPlaneBestTrack THcDC.cxx --------- 1) Eliminate filling of best focal plane quantities in CoarseTrack which defined "best" by lowest chi2 2) Create method SetFocalPlaneBestTrack a)fills the best focal plane quantities using the golden track b) fills fResiduals using golden track
-
- Mar 13, 2017
-
-
hallc-online authored
Add variables fX_fp_best,fY_fp_best,fXp_fp_best,fYp_fp_best to THcDC.h Modified THcDC::CoarseTrack so that the variables are filled with the quantities from the track with the smallest chi-squared. Modified THcDC::TrackFit so that the fResidual per plane is filled with for the track with the lowest chi-squared. Modified THcDC::DefineVariables so that access the variables as dc.x_fp,dc.y_fp,dc.xp_fp and dc.yp_fp
-
- Feb 09, 2017
-
-
Jure Bericic authored
Applied changes in several files to reflect this. Minor fix in `THcRawHodoHit`.
-
- Feb 08, 2017
-
-
Eric Pooser authored
-
Eric Pooser authored
When calculating the residuals, consider all track candidates instead of only the first track. Moreover, eliminate the instances when there was no hit in a plane then the residual was returned as zero.
-
- Dec 15, 2016
-
-
Jure Bericic authored
-
- Nov 23, 2016
-
-
Stephen A. Wood authored
Otherwise InitHistlist doesn't have the detector map and can't do setup of reference time stuff or sanity checks on the detector mapping
-
- Nov 15, 2016
-
-
Jure Bericic authored
Some detector classes lacked the ability to get the appropriate Apparatus character and were only working for HMS.
-
- May 12, 2016
-
-
Stephen A. Wood authored
-
- Jul 31, 2015
-
-
Stephen A. Wood authored
-
- May 16, 2015
-
-
Zafar Ahmed authored
-
- Mar 10, 2015
-
-
Stephen A. Wood authored
Modify THcFormula and THcParmList to match changes to THaFormula (needed for reports) Small changes in THcInterface and detector classes for compatibility with new OO decoder and podd 1.6
-
- Jan 16, 2015
-
-
Stephen A. Wood authored
When it was UInt_t, variable length global variables didn't get the length correctly and there were run time errors.
-
Stephen A. Wood authored
For the most part this is making sure that the loop variable in a for loop is of the same type as the variable that it is compared against (Often a UInt_t).
-
- Jan 14, 2015
-
-
Stephen A. Wood authored
This reverts commit e9ce1040. Conflicts: src/THcHallCSpectrometer.cxx
-
- Jan 08, 2015
-
-
Stephen A. Wood authored
For the most part this is making sure that the loop variable in a for loop is of the same type as the variable that it is compared against (Often a UInt_t).
-
- Jan 07, 2015
-
-
Mark Jones authored
-
- Sep 09, 2014
-
-
Stephen A. Wood authored
-
- Sep 02, 2014
-
-
Mark Jones authored
Found when doing valgrind --tool=memcheck --track-origins=yes --leak-check=full hcana Return errors like below ==14467== Invalid write of size 1 ==14467== at 0x4A0794F: __GI_strcpy (mc_replace_strmem.c:313) ==14467== by 0x57137A3: THcHodoscope::Setup(char const*, char const*) (THcHodoscope.cxx:112) ==14467== by 0x5713A7E: THcHodoscope::Init(TDatime const&) (THcHodoscope.cxx:142) ==14467== by 0x4D38C29: THaApparatus::Init(TDatime const&) (THaApparatus.C:193) ==14467== by 0x4D67E57: THaAnalyzer::InitModules(TList const*, TDatime&, int, char const*) (THaAnalyzer.C:433) ==14467== by 0x4D68B1F: THaAnalyzer::DoInit(THaRunBase*) (THaAnalyzer.C:667) ==14467== by 0x4D68126: THaAnalyzer::Init(THaRunBase*) (THaAnalyzer.C:467) ==14467== by 0x4D6B31A: THaAnalyzer::Process(THaRunBase*) (THaAnalyzer.C:1289) ==14467== by 0x4E21E97: G__haDict_704_0_6(G__value*, char const*, G__param*, int) (haDict.C:23491) ==14467== by 0x6DB6148: Cint::G__ExceptionWrapper(int (*)(G__value*, char const*, G__param*, int), G__value*, char*, G__param*, int) (in /u/apps/root/5.34.13/root/lib/libCint.so.5.34) ==14467== by 0x6E5A3A0: G__execute_call (in /u/apps/root/5.34.13/root/lib/libCint.so.5.34) ==14467== by 0x6E5B201: G__call_cppfunc (in /u/apps/root/5.34.13/root/lib/libCint.so.5.34) ==14467== Address 0xd7636a2 is 0 bytes after a block of size 2 alloc'd ==14467== at 0x4A0674C: operator new[](unsigned long) (vg_replace_malloc.c:305) ==14467== by 0x5713757: THcHodoscope::Setup(char const*, char const*) (THcHodoscope.cxx:111) ==14467== by 0x5713A7E: THcHodoscope::Init(TDatime const&) (THcHodoscope.cxx:142) ==14467== by 0x4D38C29: THaApparatus::Init(TDatime const&) (THaApparatus.C:193) ==14467== by 0x4D67E57: THaAnalyzer::InitModules(TList const*, TDatime&, int, char const*) (THaAnalyzer.C:433) ==14467== by 0x4D68B1F: THaAnalyzer::DoInit(THaRunBase*) (THaAnalyzer.C:667) ==14467== by 0x4D68126: THaAnalyzer::Init(THaRunBase*) (THaAnalyzer.C:467) ==14467== by 0x4D6B31A: THaAnalyzer::Process(THaRunBase*) (THaAnalyzer.C:1289) ==14467== by 0x4E21E97: G__haDict_704_0_6(G__value*, char const*, G__param*, int) (haDict.C:23491) ==14467== by 0x6DB6148: Cint::G__ExceptionWrapper(int (*)(G__value*, char const*, G__param*, int), G__value*, char*, G__param*, int) (in /u/apps/root/5
-
- Aug 28, 2014
-
-
Stephen A. Wood authored
-
- Jul 30, 2014
-
-
Jonathan Stelzleni authored
In tree they are called tnhit and trawhit. Created two new variables for THcDriftChamber, fNthits and fN_True_RawHits. In tree they are called tnhit and trawhit. THcDC defined subdetector THcDriftChamber to have name Ch1 and Ch2 to do used in the tree. THcDriftChamberPlane defined new method GetNRawhits to return fNRawhits for each plane. Added member fNRawhits to THcDriftChamberPlane
-