GMT (#680)

Commit GMT-related classes and macroses. Also, update StarDb/geometry
for GMT with the latest alignment.

---------

Co-authored-by: Grigory Nigmatkulov <[email protected]>
Co-authored-by: Grigory Nigmatkulov <[email protected]>
Co-authored-by: Yuri Fisyak <[email protected]>
Co-authored-by: Dmitri Smirnov <[email protected]>

Author: 5 people

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

@@ -1429,6 +1429,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","","","","","",                                                STAR_CHAIN_OBSOLETE,kFALSE},
+
   // EPD
   {"epdHit",    "", "", "epdDb,event",            "StEpdHitMaker", "StEpdHitMaker","EPD hit maker", kFALSE},
 
@@ -1964,6 +1971,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", "StGmtAlignmentMaker","StGmtAlignmentMaker","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,56 @@
+/**
+ * \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"
+
+//________________
+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;
+}