Skip to content
Snippets Groups Projects
  1. Feb 19, 2018
    • Stephen A. Wood's avatar
      Create parameters to hold missing reference time counts. · 14523959
      Stephen A. Wood authored
        So that the counts can be reported in report files.
      14523959
    • Stephen A. Wood's avatar
      Put a minimum cut on reference times. · 39317b13
      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
      39317b13
  2. Feb 15, 2018
  3. Dec 06, 2017
    • Stephen A. Wood's avatar
      Suppress missing ref time messages when ref time not expected · e2c0eb6f
      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.
      e2c0eb6f
  4. Dec 02, 2017
  5. Nov 03, 2017
    • Mark Jones's avatar
      Modify THcDC::Setup · decdb99f
      Mark Jones authored
      Set fHMSStyleChambers =1 according to dc_version parameter.
      If no dc_version parameter loaded defaults to old HMS style chambers
      decdb99f
  6. Nov 02, 2017
  7. Oct 09, 2017
  8. Aug 31, 2017
  9. Jul 28, 2017
  10. Jul 21, 2017
  11. May 01, 2017
    • hallc-online's avatar
      Modify THcDC and THcDriftChamberPlane to calculate per wire efficiency · 36c76b27
      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
      36c76b27
  12. Apr 26, 2017
    • hallc-online's avatar
      Modify THcDC and THcHallCSpectrometer · d5134c3d
      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
      d5134c3d
  13. Mar 13, 2017
    • hallc-online's avatar
      Modify THcDC::CoarseTrack and THcDC::TrackFit · b140267a
      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
      b140267a
  14. Feb 09, 2017
  15. Feb 08, 2017
  16. Dec 15, 2016
  17. Nov 23, 2016
  18. Nov 15, 2016
  19. May 12, 2016
  20. Jul 31, 2015
  21. May 16, 2015
  22. Mar 10, 2015
  23. Jan 16, 2015
  24. Jan 14, 2015
  25. Jan 08, 2015
  26. Jan 07, 2015
  27. Sep 09, 2014
  28. Sep 02, 2014
    • Mark Jones's avatar
      Add "+1" to length when using new char[] · 085928d1
      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
      085928d1
  29. Aug 28, 2014
  30. Jul 30, 2014
    • Jonathan Stelzleni's avatar
      Created two new variables for THcDC, fNthits and fN_True_RawHits. · d496e1e3
      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
      d496e1e3
  31. Feb 27, 2014
    • Mark Jones's avatar
      Updated drift chamber debug messages · 30ccc1ce
      Mark Jones authored
      add flags debugflagpr and debugflagstubs
      add method THcDC::PrintStubs
      modify THcDC::PrintSpacePoints
      Set THcDC::LinkStubs() to only set fNDCTracks=0 and retrun when
        fNDCTracks>10 for HMS style chambers to match engine.
      30ccc1ce
  32. Jan 30, 2014
  33. Jan 17, 2014
  34. Jan 14, 2014
  35. Dec 19, 2013
  36. Dec 16, 2013
Loading