GMT

StRoot/StAnalysisMaker/StAnalysisMaker.cxx

@@ -1141,13 +1141,14 @@ void StAnalysisMaker::summarizeEvent(StEvent *event, Int_t mEventCounter) {
   if (event->numberOfPsds()) {
     LOG_QA << "# PSDs:                " << event->numberOfPsds() << endm;
   }
-#ifdef _ST_GMT_HIT_H_
+
   if (event->gmtCollection() && event->gmtCollection()->getNumHits()) {
     LOG_QA << "# GMT hits:           " << event->gmtCollection()->getNumHits() 
-	   << " points: " << event->gmtCollection()->getNumPoints()
+	   << ", points: " << event->gmtCollection()->getNumPoints()
+     << ", strips: " << event->gmtCollection()->getNumStrips()
 	   << endm;
   }
-#endif /* _ST_GMT_HIT_H_ */
+
   if (event->rpsCollection()) {
     LOG_QA << "# RPS hits:            " << event->rpsCollection()->clusters().size() << endm;
   }

StRoot/StBFChain/BigFullChain.h

@@ -1381,6 +1381,13 @@ Bfc_st BFC[] = { // standard chains
    ,                                                                                "MTD hit maker",kFALSE},
   {"mtdTrkMask","","","db","StMtdTrackingMaskMaker"      ,"StMtdEvtFilterMaker","MTD track masking",kFALSE},
 
+  // GMT
+  {"gmt"        ,"","","gmtDat,gmtClu"                                      ,"","","Gmt data Chain",kFALSE},
+  {"gmtDat"     ,"","","event","StGmtRawMaker","StGmtRawMaker",                   "GMT Data reader",kFALSE},
+  {"gmtClu"   ,"","","gmtutil","StGmtClusterMaker","Spectrum,StGmtClusterMaker","GMT cluster maker",kFALSE},
+  {"gmtCosmics" ,"","","Cosmics,gmt","","",  "Save only events with GMT clusters and Cosmic tracks",kFALSE},
+  {"gmtClusTree","","","","","",                               "WARNING *** Option is OBSOLETE ***",kFALSE},
+
   // EPD
   {"epdHit",    "", "", "epdDb,event",            "StEpdHitMaker", "StEpdHitMaker","EPD hit maker", kFALSE},
 
@@ -1782,6 +1789,10 @@ Bfc_st BFC[] = { // standard chains
   {"mtdCalib"    ,"","","db",                 "StMtdCalibMaker","StMtdCalibMaker","MTD calibration",kFALSE},
   {"mtdEvtFilt"  ,"","","db",       "StMtdEvtFilterMaker","StMtdEvtFilterMaker","MTD event filter" ,kFALSE},
 
+  // GMT
+  {"gmtMatch"                                 ,"","","", "","","WARNING *** Option is OBSOLETE ***",kFALSE},
+  {"gmtPlotter"                               ,"","","", "","","WARNING *** Option is OBSOLETE ***",kFALSE},
+
   {"FindVtxSeed"   ,"FindVtxSeed"   ,"","globT,MuDSTDeps,picoDst","StVertexSeedMaker"
    ,                                   "StPass0CalibMaker",          "Performs vertex seed finding",kFALSE},
   {"FindEvtVtxSeed","FindEvtVtxSeed","","MuDSTDeps,picoDst","StEvtVtxSeedMaker"
@@ -1913,6 +1924,8 @@ Bfc_st BFC[] = { // standard chains
    ,                              "StSvtPoolEventT,StSvtPoolSvtMatchedTree","Create SvtMatchedTree",kFALSE},
   {"LAna"        ,"","","in,detDb,StEvent,tpcDb","StLaserAnalysisMaker"
    ,                                                   "StLaserAnalysisMaker","Laser data Analysis",kFALSE},
+  // GMT
+  {"gmtAligner"    ,"","","detDb",             "StGmtAligner","StGmtAligner","GMT cluster plotting",kFALSE},
   {"EandBDir","","","in,StEvent,TpcHitMover,nodefault"
    ,   "StEandBDirMaker","MathMore,Spectrum,StEandBDirMaker",                   "get E&B direction",kFALSE},
   {"SpinTag"     ,"","",""                                               ,"","",STAR_CHAIN_OBSOLETE,kFALSE},

StRoot/StBFChain/StBFChain.cxx

@@ -550,6 +550,11 @@ Int_t StBFChain::Instantiate()
       if (GetOption("Cosmics"))   mk->SetAttr("Cosmics"    ,kTRUE);
       mk->PrintAttr();
     }
+
+    if (maker== "StGmtClusterMaker") {
+      if (GetOption("gmtCosmics"))  mk->SetAttr("gmtCosmics",  kTRUE);
+    }
+
     if (maker=="StKFVertexMaker" && GetOption("laserIT"))   mk->SetAttr("laserIT"    ,kTRUE);
     //		Sti(ITTF) end
     if (maker=="StGenericVertexMaker") {

StRoot/StEvent/StContainers.cxx

@@ -198,6 +198,9 @@
 #include "StTrackPidTraits.h"
 #include "StV0Vertex.h"
 #include "StXiVertex.h"
+#include "StGmtHit.h"
+#include "StGmtPoint.h"
+#include "StGmtStrip.h"
 
 StCollectionImp(BTofHit)
 StCollectionImp(BTofRawHit)
@@ -277,3 +280,6 @@ StCollectionImp(TrackNode)
 StCollectionImp(TrackPidTraits)
 StCollectionImp(V0Vertex)
 StCollectionImp(XiVertex)
+StCollectionImp(GmtHit)
+StCollectionImp(GmtStrip)
+StCollectionImp(GmtPoint)

StRoot/StEvent/StContainers.h

@@ -214,6 +214,9 @@ class StTrackNode;
 class StTrackPidTraits;
 class StV0Vertex;
 class StXiVertex;
+class StGmtHit;
+class StGmtStrip;
+class StGmtPoint;
 
 StCollectionDef(BTofHit)
 StCollectionDef(BTofRawHit)
@@ -293,6 +296,9 @@ StCollectionDef(TrackNode)
 StCollectionDef(TrackPidTraits)
 StCollectionDef(V0Vertex)
 StCollectionDef(XiVertex)
+StCollectionDef(GmtHit)
+StCollectionDef(GmtStrip)
+StCollectionDef(GmtPoint)
 
 
 #endif

StRoot/StEvent/StEvent.cxx

@@ -256,6 +256,7 @@
 #include "StTrackNode.h"
 #include "StTrack.h"
 #include "StFwdTrackCollection.h"
+#include "StGmtCollection.h"
 
 #ifndef ST_NO_NAMESPACES
 using std::swap;
@@ -908,6 +909,23 @@ StEvent::fgtCollection() const
     return fgtCollection;
 }
 
+StGmtCollection*
+StEvent::gmtCollection()
+{
+    StGmtCollection *gmtCollection = 0;
+    _lookup(gmtCollection, mContent);
+    return gmtCollection;
+}
+
+const StGmtCollection*
+StEvent::gmtCollection() const
+{
+    StGmtCollection *gmtCollection = 0;
+    _lookup(gmtCollection, mContent);
+    return gmtCollection;
+}
+
+
 StIstHitCollection*
 StEvent::istHitCollection()
 {
@@ -1456,6 +1474,12 @@ StEvent::setFgtCollection(StFgtCollection* val)
     _lookupAndSet(val, mContent);
 }
 
+void
+StEvent::setGmtCollection(StGmtCollection* val)
+{
+    _lookupAndSet(val, mContent);
+}
+
 void
 StEvent::setIstHitCollection(StIstHitCollection* val)
 {
@@ -1629,6 +1653,7 @@ void StEvent::statistics()
     cout << "\t# of hits in EMC:            " << (emcCollection() ? emcCollection()->barrelPoints().size() : 0) << endl;
     cout << "\t# of hits in EEMC:           " << (emcCollection() ? emcCollection()->endcapPoints().size() : 0) << endl;
     cout << "\t# of hits in FGT:            " << (fgtCollection() ? fgtCollection()->getNumHits() : 0) << endl;
+    cout << "\t# of hits in GMT:            " << (gmtCollection() ? gmtCollection()->getNumHits() : 0) << endl;
     cout << "\t# of hits in RICH:           " << (richCollection() ? richCollection()->getRichHits().size() : 0) << endl;
     cout << "\t# of PSDs:                   " << numberOfPsds() << endl;
 }

StRoot/StEvent/StEvent.h

@@ -228,6 +228,7 @@ class StIstHitCollection;
 class StFstEvtCollection;
 class StFstHitCollection;
 class StFwdTrackCollection;
+class StGmtCollection;
 
 class StEvent : public StXRefMain {
 public:
@@ -318,6 +319,8 @@ class StEvent : public StXRefMain {
     const StTriggerIdCollection*        triggerIdCollection() const;
     StTriggerData*                      triggerData();
     const StTriggerData*                triggerData() const;
+    StGmtCollection*                    gmtCollection();
+    const StGmtCollection*              gmtCollection() const;
 
     StSPtrVecTrackDetectorInfo&         trackDetectorInfo();
     const StSPtrVecTrackDetectorInfo&   trackDetectorInfo() const;
@@ -412,6 +415,7 @@ class StEvent : public StXRefMain {
     void removePsd(StPsd*);
     void addHitCollection(StSPtrVecHit* p, const Char_t *name);
     void removeHitCollection(const Char_t *name);
+    void setGmtCollection(StGmtCollection*);
 
     virtual Bool_t Notify();
 

StRoot/StEvent/StEventClusteringHints.cxx

@@ -202,6 +202,7 @@ StEventClusteringHints::StEventClusteringHints()
     setBranch("StRpsCollection",             "evt_aux",      7);
     setBranch("StFttCollection",             "evt_aux",      7);
     setBranch("StFstEvtCollection",          "evt_aux",      7);
+    setBranch("StGmtCollection",             "evt_aux",      7);
     setBranch("StSsdHitCollection",          "evt_hits",     8);
     setBranch("StSstHitCollection",          "evt_hits",     8);
     setBranch("StSvtHitCollection",          "evt_hits",     8);
@@ -211,6 +212,7 @@ StEventClusteringHints::StEventClusteringHints()
     setBranch("StTpcHitCollection",          "evt_hits",     8);
     setBranch("StFtpcHitCollection",         "evt_hits",     8);
     setBranch("StRnDHitCollection",          "evt_hits",     8);
+    setBranch("StGmtHitCollection",          "evt_hits",     8);
     setBranch("StHltEvent",                  "evt_hlt",      9);
     setBranch("StFgtCollection",             "evt_fgt",      9);
 

StRoot/StEvent/StEventTypes.h

@@ -169,6 +169,8 @@
  **************************************************************************/
 #ifndef StEventTypes_hh
 #define StEventTypes_hh
+
+#include "StContainers.h"
 
 #include "StBbcTriggerDetector.h"
 #include "StCalibrationVertex.h"
@@ -330,6 +332,8 @@
 #include "StFgtHitCollection.h"
 #include "StFgtPoint.h"
 #include "StFgtPointCollection.h"
+#include "StGmtHit.h"
+#include "StGmtCollection.h"
 #include "StIstHit.h"
 #include "StIstSensorHitCollection.h" 
 #include "StIstLadderHitCollection.h"

StRoot/StEvent/StGmtCollection.cxx

@@ -0,0 +1,58 @@
+/**
+ * \class StGmtColection
+ * \brief Holds collections of GMT data
+ * 
+ * GMT data collection for StEvent (based on StFgtCollection)
+ *
+ * \author K.S. Engle, Jan. 2013
+ * \author Richard Witt ([email protected]), Jan. 2013
+ * \author Grigory Nigmatkulov ([email protected]), Dec. 2020
+ */
+
+// Load header of the GMT collection
+#include "StGmtCollection.h"
+
+ClassImp(StGmtCollection)
+
+//________________
+StGmtCollection::StGmtCollection() : StObject() {
+  // Set the module field for some of the collections
+  for( int i=0; i<kGmtNumModules; ++i ) {
+    mStripCollection[i].setModule( i );
+    mHitCollection[i].setModule( i );
+  }
+}
+
+//________________
+StGmtCollection::~StGmtCollection() {
+  /* empty */
+}
+
+//__________
+void StGmtCollection::Clear( Option_t *opt ) {
+  for( int i=0; i<kGmtNumModules; ++i ){
+    mStripCollection[i].Clear( opt );
+    mHitCollection[i].Clear( opt );
+  }
+  mPointCollection.Clear( opt );
+}
+
+//________________
+size_t StGmtCollection::getNumStrips() const { 
+  size_t n = 0;
+  for( const StGmtStripCollection* ptr = &mStripCollection[0]; 
+       ptr != &mStripCollection[kGmtNumModules]; ++ptr ) {
+    n += ptr->getNumStrips();
+  }
+  return n;
+}
+
+//________________
+size_t StGmtCollection::getNumHits() const {
+  size_t n = 0;
+  for( const StGmtHitCollection* ptr = &mHitCollection[0]; 
+       ptr != &mHitCollection[kGmtNumModules]; ++ptr ) {
+    n += ptr->getNumHits();
+  }
+  return n;
+}

StRoot/StEvent/StGmtCollection.h

@@ -0,0 +1,90 @@
+/**
+ * \class StGmtCollection
+ * \brief Holds collections of GMT data
+ * 
+ * GMT data collection for StEvent (based on StFgtCollection)
+ *
+ * \author K.S. Engle, Jan. 2013
+ * \author Richard Witt ([email protected]), Jan. 2013
+ * \author Grigory Nigmatkulov ([email protected]), Dec. 2020
+ */
+
+#ifndef StGmtCollection_hh
+#define StGmtCollection_hh
+
+// StRoot headers
+#include "StObject.h"
+#include "StGmtStripCollection.h"
+#include "StGmtHitCollection.h"
+#include "StGmtPointCollection.h"
+#include "StEnumerations.h"
+
+//________________
+class StGmtCollection : public StObject {
+ public:
+  /// Constructor
+  StGmtCollection();
+  // StGmtCollection( const StGmtCollection& other );            ---> use default
+  // StGmtCollection& operator=( const StGmtCollection& other ); ---> use default 
+
+  /// Destructor
+  ~StGmtCollection();
+
+  /// Number of modules
+  size_t getNumModules() const                             { return kGmtNumModules; }
+  /// Number total number of strips
+  size_t getNumStrips() const;
+  /// Number of strips in the i-th module
+  size_t getNumStrips( unsigned short moduleIdx) const
+  { return (moduleIdx < kGmtNumModules ? mStripCollection[moduleIdx].getNumStrips() : 0 ); }
+  /// Number total number of hits
+  size_t getNumHits() const;
+  /// Number of hits in the i-th module
+  size_t getNumHits( unsigned short moduleIdx ) const
+  { return (moduleIdx < kGmtNumModules ? mHitCollection[moduleIdx].getNumHits() : 0 ); }
+  /// Number of points
+  size_t getNumPoints() const
+  { return mPointCollection.getNumPoints(); }
+
+  // note: ownership of all pointers is retained by the containers.
+  // Do not delete any pointers received from this class.
+
+  /// Pointer to the GMT strip collection in the i-th module
+  StGmtStripCollection* getStripCollection( unsigned short moduleIdx )
+  { return (moduleIdx < kGmtNumModules ? &mStripCollection[moduleIdx] : 0 ); }
+  /// Pointer to the GMT strip collection in the i-th module  
+  const StGmtStripCollection* getStripCollection( unsigned short moduleIdx ) const
+  { return (moduleIdx < kGmtNumModules ? &mStripCollection[moduleIdx] : 0 ); }
+
+  /// Pointer to the GMT hit collection in the i-th module  
+  StGmtHitCollection* getHitCollection( unsigned short moduleIdx )
+  { return (moduleIdx < kGmtNumModules ? &mHitCollection[moduleIdx] : 0 ); }
+  /// Pointer to the GMT hit collection in the i-th module  
+  const StGmtHitCollection* getHitCollection( unsigned short moduleIdx ) const
+  { return (moduleIdx < kGmtNumModules ? &mHitCollection[moduleIdx] : 0 ); }
+
+  /// Pointer to the GMT point collection
+  StGmtPointCollection* getPointCollection()             { return &mPointCollection; }
+  /// Pointer to the GMT point collection
+  const StGmtPointCollection* getPointCollection() const { return &mPointCollection; }
+
+  /// Clear method  
+  void Clear( Option_t *opts = "" );
+
+ protected:
+
+  /// Friend class of the StMuDSTMaker
+  friend class StMuDstMaker;
+
+  /// GMT strip collections for all modules
+  StGmtStripCollection mStripCollection[kGmtNumModules];
+  /// GMT hit collections for all modules
+  StGmtHitCollection mHitCollection[kGmtNumModules];
+  /// GMT point collection
+  StGmtPointCollection mPointCollection;
+
+ private:   
+  ClassDef(StGmtCollection,1)
+}; 
+
+#endif // #define StGmtCollection_hh