Compare revisions

e974befb · e974befb · e974befb · e974befb · e974befb · e974befb
--- a/src/THcHodoscope.h
+++ b/src/THcHodoscope.h
@@ -76,11 +76,14 @@ public:
  virtual Int_t      FineProcess( TClonesArray& tracks );
  virtual Int_t      End(THaRunBase* run=0);

+  Double_t DetermineTimePeak(Int_t FillFlag);
  void EstimateFocalPlaneTime(void);
  void OriginalTrackEffTest(void);
  void TrackEffTest(void);
+  void CalcCluster(void);
  virtual Int_t      ApplyCorrections( void );
  Double_t GetStartTime() const { return fStartTime; }
+  Double_t GetOffsetTime() const { return fOffsetTime; }
  Bool_t IsStartTimeGood() const {return fGoodStartTime;};
  Int_t GetNfptimes() const {return fNfptimes;};
  Int_t GetScinIndex(Int_t nPlane, Int_t nPaddle);
@@ -192,12 +195,19 @@ protected:
  Bool_t fSHMS;
  Bool_t fGoodStartTime;
  Double_t fStartTime;
+  Double_t fADCStartTime;
+  Double_t fOffsetTime;
  Double_t fFPTimeAll;
  Int_t fNfptimes;
  Bool_t* fPresentP;
-  Double_t fTimeHist_Peak;
-  Double_t fTimeHist_Sigma;
-  Double_t fTimeHist_Hits;
+  Double_t fTimeHist_StartTime_NumPeaks;
+  Double_t fTimeHist_StartTime_Peak;
+  Double_t fTimeHist_StartTime_Sigma;
+  Double_t fTimeHist_StartTime_Hits;
+  Double_t fTimeHist_FpTime_NumPeaks;
+  Double_t fTimeHist_FpTime_Peak;
+  Double_t fTimeHist_FpTime_Sigma;
+  Double_t fTimeHist_FpTime_Hits;

  Double_t     fBeta;

@@ -207,7 +217,7 @@ protected:

  // Potential Hall C parameters.  Mostly here for demonstration
  Int_t fNPlanes;		// Number of planes
-  UInt_t fMaxScinPerPlane,fMaxHodoScin; // max number of scin/plane; product of the first two
+  UInt_t fMaxScinPerPlane,fMaxHodoScin,fTotHodScin; // max number of scin/plane; product of the first two
  Double_t fStartTimeCenter, fStartTimeSlop, fScinTdcToTime;
  Double_t fTofTolerance;
  Int_t fCosmicFlag; //
@@ -216,6 +226,7 @@ protected:
  Double_t fScinTdcMin, fScinTdcMax; // min and max TDC values
  char** fPlaneNames;
  UInt_t* fNPaddle;		// Number of paddles per plane
+  Int_t fTrackBetaIncludeSinglePmtHits;

  Double_t *fHodoNegAdcTimeWindowMin;    
  Double_t *fHodoNegAdcTimeWindowMax;
@@ -265,6 +276,7 @@ protected:

  Int_t        fCheckEvent;
  Int_t        fEventType;
+  Int_t        fEventNum;

  Int_t        fGoodTrack;
  Double_t     fScin2XZpos;
@@ -407,6 +419,10 @@ protected:
  std::vector<Int_t > fNClust;		                // # scins clusters for the plane
  std::vector<std::vector<Int_t> > fClustSize;		                // # scin cluster size
  std::vector<std::vector<Double_t> > fClustPos;		                // # scin cluster position
+  std::vector<Int_t > fNCluster;		                // # scins clusters for the plane
+  std::vector<std::vector<Int_t> > fClusterSize;		                // # scin cluster size
+  std::vector<std::vector<Double_t> > fClusterXPos;		                // # scin cluster position
+  std::vector<std::vector<Double_t> > fClusterYPos;		                // # scin cluster position
  std::vector<Int_t > fThreeScin;	                // # scins three clusters for the plane
  std::vector<Int_t > fGoodScinHitsX;                   // # hits in fid x range
  // Could combine the above into a structure

--- a/src/THcParmList.cxx
+++ b/src/THcParmList.cxx
--- a/src/THcRawAdcHit.cxx
+++ b/src/THcRawAdcHit.cxx
@@ -258,6 +258,10 @@ void THcRawAdcHit::SetRefTime(Int_t refTime) {
  fHasRefTime = kTRUE;
 }

+void THcRawAdcHit::SetRefDiffTime(Int_t refDiffTime) {
+  fRefDiffTime = refDiffTime;
+}
+
 void THcRawAdcHit::SetSample(Int_t data) {
  if (fNSamples >= fMaxNSamples) {
    // throw std::out_of_range("`THcRawAdcHit::SetSample`: too many samples!");
@@ -481,6 +485,22 @@ Int_t THcRawAdcHit::GetRefTime() const {
  }
 }

-Bool_t THcRawAdcHit::HasRefTime() const { return fHasRefTime; }
+
+Int_t THcRawAdcHit::GetRefDiffTime() const {
+  if (fHasRefTime) {
+    return fRefDiffTime;
+  }
+  else {
+    TString msg = TString::Format(
+      "`THcRawAdcHit::GetRefTime`: Reference time not available!"
+    );
+    throw std::runtime_error(msg.Data());
+  }
+}
+
+Bool_t THcRawAdcHit::HasRefTime() const {
+  return fHasRefTime;
+}
+

 ClassImp(THcRawAdcHit)
--- a/src/THcRawAdcHit.h
+++ b/src/THcRawAdcHit.h
@@ -16,12 +16,14 @@ class THcRawAdcHit : public podd2::HitLogging<TObject> {
    void SetData(Int_t data);
    void SetSample(Int_t data);
    void SetRefTime(Int_t refTime);
+    void SetRefDiffTime(Int_t refDiffTime);
    void SetDataTimePedestalPeak(
      Int_t data, Int_t time, Int_t pedestal, Int_t peak
    );

    Int_t GetRawData(UInt_t iPulse=0) const;
    Double_t GetF250_PeakPedestalRatio() {return fPeakPedestalRatio;};
+    Int_t GetF250_NPedestalSamples() {return fNPedestalSamples;};

    Double_t GetAverage(UInt_t iSampleLow, UInt_t iSampleHigh) const;
    Int_t GetIntegral(UInt_t iSampleLow, UInt_t iSampleHigh) const;
@@ -41,6 +43,7 @@ class THcRawAdcHit : public podd2::HitLogging<TObject> {
    Int_t GetPulseTimeRaw(UInt_t iPulse=0) const;
    Int_t GetSampleRaw(UInt_t iSample=0) const;
    Int_t GetRefTime() const;
+    Int_t GetRefDiffTime() const;

    Double_t GetPed() const;
    Double_t GetPulseInt(UInt_t iPulse=0) const;
@@ -82,6 +85,7 @@ class THcRawAdcHit : public podd2::HitLogging<TObject> {
    Int_t fPulseTime[fMaxNPulses];
    Int_t fSample[fMaxNSamples]; // the big buffer
    Int_t fRefTime;
+    Int_t fRefDiffTime;

    Bool_t fHasMulti;
    Bool_t fHasRefTime;

--- a/src/THcRawDCHit.cxx
+++ b/src/THcRawDCHit.cxx
@@ -59,6 +59,17 @@ void THcRawDCHit::SetReference(Int_t signal, Int_t reference) {
  }
 }

+void THcRawDCHit::SetReferenceDiff(Int_t signal, Int_t reference) {
+  if (signal == 0) {
+    fTdcHit.SetRefDiffTime(reference);
+  }
+  else {
+    throw std::out_of_range(
+      "`THcRawDCHit::SetReference`: only signal `0` available!"
+    );
+  }
+}
+

 Int_t THcRawDCHit::GetData(Int_t signal) {
  if (signal == 0) {
@@ -95,6 +106,17 @@ Int_t THcRawDCHit::GetReference(Int_t signal) {
  }
 }

+Int_t THcRawDCHit::GetReferenceDiff(Int_t signal) {
+  if (signal == 0) {
+    return fTdcHit.GetRefDiffTime();
+  }
+  else {
+    throw std::out_of_range(
+      "`THcRawDCHit::GetReference`: only signal `0` available!"
+    );
+  }
+}
+

 THcRawHit::ESignalType THcRawDCHit::GetSignalType(Int_t signal) {
  if (signal == 0) {

--- a/src/THcRawDCHit.h
+++ b/src/THcRawDCHit.h
@@ -18,10 +18,12 @@ class THcRawDCHit : public THcRawHit {

    virtual void SetData(Int_t signal, Int_t data);
    virtual void SetReference(Int_t signal, Int_t reference);
+    virtual void SetReferenceDiff(Int_t signal, Int_t reference);

    virtual Int_t GetData(Int_t signal);
    virtual Int_t GetRawData(Int_t signal);
    virtual Int_t GetReference(Int_t signal);
+    virtual Int_t GetReferenceDiff(Int_t signal);
    virtual ESignalType GetSignalType(Int_t signal);
    virtual Int_t GetNSignals();


--- a/src/THcRawHit.h
+++ b/src/THcRawHit.h
@@ -38,9 +38,10 @@ public:
  virtual ESignalType GetSignalType(Int_t /* signal */) {return kUndefined;}
  virtual Int_t GetNSignals() { return 1;}

-  virtual void SetReference(Int_t /* signal */, Int_t /* reference */) {};
-  virtual Bool_t HasReference(Int_t /* signal */) {return kFALSE;};
-  virtual Int_t GetReference(Int_t /* signal */) {return 0;};
+  virtual void SetReference(Int_t signal, Int_t reference) {};
+  virtual void SetReferenceDiff(Int_t signal, Int_t reference) {};
+  virtual Bool_t HasReference(Int_t signal) {return kFALSE;};
+  virtual Int_t GetReference(Int_t signal) {return 0;};

  virtual void SetF250Params(Int_t /* NSA */, Int_t /* NSB */, Int_t /* NPED */) {};


--- a/src/THcRawHodoHit.cxx
+++ b/src/THcRawHodoHit.cxx
@@ -105,6 +105,18 @@ void THcRawHodoHit::SetReference(Int_t signal, Int_t reference) {
  }
 }

+void THcRawHodoHit::SetReferenceDiff(Int_t signal, Int_t referenceDiff) {
+  if (signal < fNAdcSignals) {
+    fAdcHits[signal].SetRefDiffTime(referenceDiff);
+  } else if (signal < fNAdcSignals+fNTdcSignals) {
+    fTdcHits[signal-fNAdcSignals].SetRefDiffTime(referenceDiff);
+  } else {
+    throw std::out_of_range(
+      "`THcRawHodoHit::SetReference`: only signals `2` and `3` available!"
+    );
+  }
+}
+

 Int_t THcRawHodoHit::GetData(Int_t signal) {
  if (0 <= signal && signal < fNAdcSignals) {
@@ -147,6 +159,17 @@ Int_t THcRawHodoHit::GetReference(Int_t signal) {
  }
 }

+Int_t THcRawHodoHit::GetReferenceDiff(Int_t signal) {
+  if (fNAdcSignals <= signal && signal < fNAdcSignals+fNTdcSignals) {
+    return fTdcHits[signal-fNAdcSignals].GetRefDiffTime();
+  }
+  else {
+    throw std::out_of_range(
+      "`THcRawHodoHit::GetReference`: only signals `2` and `3` available!"
+    );
+  }
+}
+

 THcRawHit::ESignalType THcRawHodoHit::GetSignalType(Int_t signal) {
  if (0 <= signal && signal < fNAdcSignals) {

--- a/src/THcRawHodoHit.h
+++ b/src/THcRawHodoHit.h
@@ -25,10 +25,12 @@ class THcRawHodoHit : public THcRawHit {
      Int_t signal, Int_t data, Int_t time, Int_t pedestal, Int_t peak
    );
    virtual void SetReference(Int_t signal, Int_t reference);
+    virtual void SetReferenceDiff(Int_t signal, Int_t referenceDiff);

    virtual Int_t GetData(Int_t signal);
    virtual Int_t GetRawData(Int_t signal);
    virtual Int_t GetReference(Int_t signal);
+    virtual Int_t GetReferenceDiff(Int_t signal);
    virtual ESignalType GetSignalType(Int_t signal);
    virtual Int_t GetNSignals();


--- a/src/THcRawShowerHit.cxx
+++ b/src/THcRawShowerHit.cxx
@@ -92,6 +92,16 @@ void THcRawShowerHit::SetReference(Int_t signal, Int_t reference) {
  }
 }

+void THcRawShowerHit::SetReferenceDiff(Int_t signal, Int_t reference) {
+  if (signal < fNAdcSignals) {
+    fAdcHits[signal].SetRefDiffTime(reference);
+  } else {
+    throw std::out_of_range(
+      "`THcRawHodoHit::SetReference`: only signals `2` and `3` available!"
+    );
+  }
+}
+

 Int_t THcRawShowerHit::GetData(Int_t signal) {
  if (0 <= signal && signal < fNAdcSignals) {

--- a/src/THcRawShowerHit.h
+++ b/src/THcRawShowerHit.h
@@ -22,6 +22,7 @@ class THcRawShowerHit : public THcRawHit {
      Int_t signal, Int_t data, Int_t time, Int_t pedestal, Int_t peak
    );
    virtual void SetReference(Int_t signal, Int_t reference);
+    virtual void SetReferenceDiff(Int_t signal, Int_t reference);

    virtual Int_t GetData(Int_t signal);
    virtual Int_t GetRawData(Int_t signal);

--- a/src/THcRawTdcHit.cxx
+++ b/src/THcRawTdcHit.cxx
--- a/src/THcRawTdcHit.h
+++ b/src/THcRawTdcHit.h
@@ -14,10 +14,12 @@ class THcRawTdcHit : public TObject {

    void SetTime(Int_t time);
    void SetRefTime(Int_t refTime);
+    void SetRefDiffTime(Int_t refDiffTime);

    Int_t GetTimeRaw(UInt_t iHit=0) const;
    Int_t GetTime(UInt_t iHit=0) const;
    Int_t GetRefTime() const;
+    Int_t GetRefDiffTime() const;

    Bool_t HasRefTime() const;

@@ -30,6 +32,7 @@ class THcRawTdcHit : public TObject {

    Int_t fTime[fMaxNHits];
    Int_t fRefTime;
+    Int_t fRefDiffTime;

    Bool_t fHasRefTime;
    UInt_t fNHits;

--- a/src/THcScalerEvtHandler.cxx
+++ b/src/THcScalerEvtHandler.cxx
--- a/src/THcScintillatorPlane.cxx
+++ b/src/THcScintillatorPlane.cxx
--- a/src/THcScintillatorPlane.h
+++ b/src/THcScintillatorPlane.h
--- a/src/THcShower.cxx
+++ b/src/THcShower.cxx
--- a/src/THcShower.h
+++ b/src/THcShower.h
--- a/src/THcShowerArray.cxx
+++ b/src/THcShowerArray.cxx
--- a/src/THcShowerArray.h
+++ b/src/THcShowerArray.h
@@ -150,6 +150,7 @@ protected:
  static const Int_t kADCSampIntDynPed=3;
  Double_t *fAdcTimeWindowMin ;
  Double_t *fAdcTimeWindowMax ;
+  Int_t *fPedDefault ;
  Double_t fAdcThreshold ;
  Double_t fAdcTdcOffset;
No results found