Add Hodoscope features: Time differences and beta

      Compute time differences between planes
      Calculate beta and Chisquard for each track in the hodoscope.
      New variables: fBeta, fBetaChisq
      New histograms for time differences

examples/output.def

@@ -170,6 +170,12 @@ TH1F hdc2y2_dd 'HDC 2Y2 Drift Distance' H.dc.2y2.dist 300 -0.1 0.6
 TH1F hdc2x2_dd 'HDC 2X2 Drift Distance' H.dc.2x2.dist 300 -0.1 0.6
 
 # Focal Plane times
+TH1F htimedif1 ' Time dif of plane 1 and 2' H.hod.fpTimeDif1 80 -40. 40.
+TH1F htimedif2 ' Time dif of plane 1 and 3' H.hod.fpTimeDif2 80 -40. 40.
+TH1F htimedif3 ' Time dif of plane 1 and 4' H.hod.fpTimeDif3 80 -40. 40.
+TH1F htimedif4 ' Time dif of plane 2 and 3' H.hod.fpTimeDif4 80 -40. 40.
+TH1F htimedif5 ' Time dif of plane 2 and 4' H.hod.fpTimeDif5 80 -40. 40.
+TH1F htimedif6 ' Time dif of plane 3 and 4' H.hod.fpTimeDif6 80 -40. 40.
 TH1F hs1xfptime 'HODO s1x fptime' H.hod.1x.fptime 80 0 80 H.hod.hgoodstarttime
 TH1F hs1yfptime 'HODO s1y fptime' H.hod.1y.fptime 80 0 80 H.hod.hgoodstarttime
 TH1F hs2xfptime 'HODO s2x fptime' H.hod.2x.fptime 80 0 80 H.hod.hgoodstarttime

src/THcHodoscope.h

@@ -7,6 +7,8 @@
 //                                                                           //
 ///////////////////////////////////////////////////////////////////////////////
 
+#include <vector>
+
 #include "TClonesArray.h"
 #include "THaNonTrackingDetector.h"
 #include "THcHitList.h"
@@ -51,17 +53,71 @@ public:
   Double_t GetStartTimeCenter() const {return fStartTimeCenter;}
   Double_t GetStartTimeSlop() const {return fStartTimeSlop;}
   Double_t GetBetaNotrk() const {return fBetaNotrk;}
-  Double_t GetBeta() const {return fBeta;}
+
+  //  Double_t GetBeta() const {return fBeta[];}
+
+  Double_t GetBeta(Int_t iii) const {return fBeta[iii];} // Ahmed
+
   Double_t GetHodoPosSigma(Int_t iii) const {return fHodoPosSigma[iii];}
   Double_t GetHodoNegSigma(Int_t iii) const {return fHodoNegSigma[iii];}
 
   const TClonesArray* GetTrackHits() const { return fTrackProj; }
-  
+
   friend class THaScCalib;
 
   THcHodoscope();  // for ROOT I/O
 protected:
 
+  //  static int MAXHITS = 53;
+
+  Double_t*    FPTime;     // !time at fp from all hits in 1 scintillator plane
+  Double_t     FPTimeDif1; // !time difference at fp between scintillator planes 1 & 2
+  Double_t     FPTimeDif2; // !time difference at fp between scintillator planes 1 & 3
+  Double_t     FPTimeDif3; // !time difference at fp between scintillator planes 1 & 4
+  Double_t     FPTimeDif4; // !time difference at fp between scintillator planes 2 & 3
+  Double_t     FPTimeDif5; // !time difference at fp between scintillator planes 2 & 4
+  Double_t     FPTimeDif6; // !time difference at fp between scintillator planes 3 & 4
+
+  Double_t* fScinSigma;               // Ahmed
+  Double_t* fGoodScinTime;            // Ahmed
+  Double_t* fScinTime;                // Ahmed
+  Double_t* fTime;                    // Ahmed
+  Double_t* adcPh;                    // Ahmed
+  Double_t* fBeta;                    // Ahmed
+  Double_t* fBetaChisq;               // Ahmed
+  Double_t* fTimeAtFP;                // Ahmed
+  Double_t* fPath;                    // Ahmed
+  Double_t* fTimePos;                 // Ahmed
+  Double_t* fTimeNeg;                 // Ahmed
+  Double_t* fScinTimefp;              // Ahmed
+  Double_t* fScinPosTime;             // Ahmed
+  Double_t* fScinNegTime;             // Ahmed
+  Double_t* fSumPlaneTime;            // Ahmed
+
+  Int_t* fHitPaddle;                  // Ahmed
+  Int_t* fNScinHit;                   // Ahmed
+  Int_t* fNPmtHit;                    // Ahmed
+  Int_t* fTimeHist;                   // Ahmed
+  Int_t* fNPlaneTime;                 // Ahmed
+
+  Bool_t* fScinGoodTime;              // Ahmed
+  Bool_t* fKeepPos;                   // Ahmed
+  Bool_t* fKeepNeg;                   // Ahmed
+  Bool_t* fGoodPlaneTime;             // Ahmed
+  Bool_t* fGoodTDCPos;                // Ahmed
+  Bool_t* fGoodTDCNeg;                // Ahmed
+
+  Int_t fGoodTimeIndex;               // Ahmed
+
+  TClonesArray* scinPosADC; // Ahmed
+  TClonesArray* scinNegADC; // Ahmed
+  TClonesArray* scinPosTDC; // Ahmed
+  TClonesArray* scinNegTDC; // Ahmed
+
+  TClonesArray* scinTofPosTDC; // Ahmed
+
+  //  std::vector<bool> myScinGoodTime;  // Ahmed
+
   Int_t fAnalyzePedestals;
 
   // Calibration
@@ -71,8 +127,10 @@ protected:
   Double_t fStartTime; 
   Int_t fNfptimes;
 
+  Double_t fdEdX;
+
   Double_t fBetaNotrk;
-  Double_t fBeta;  
+  //  Double_t fBeta;  
   // Per-event data
 
   // Potential Hall C parameters.  Mostly here for demonstration
@@ -138,13 +196,12 @@ protected:
   
   virtual  Double_t TimeWalkCorrection(const Int_t& paddle,
 					   const ESide side);
-
   void Setup(const char* name, const char* description);
 
-  TClonesArray* scinPosADC; // Ahmed
-  TClonesArray* scinNegADC; // Ahmed
-  TClonesArray* scinPosTDC; // Ahmed
-  TClonesArray* scinNegTDC; // Ahmed
+  //  TClonesArray* scinPosADC; // Ahmed
+  //  TClonesArray* scinNegADC; // Ahmed
+  //  TClonesArray* scinPosTDC; // Ahmed
+  //  TClonesArray* scinNegTDC; // Ahmed
 
   ClassDef(THcHodoscope,0)   // Hodoscope detector
 };