Changeset 214 in Idarraga

Timestamp:

Jul 20, 2011, 1:28:03 AM (13 years ago)

Author:

idarraga

Message:

 

Location:

EUTelraw2mdm

Files:

• USBpixConverterPlugin.cc (added)
• USBpixI4ConverterPlugin.cc (modified) (13 diffs)
• write_mdm.cc (added)

EUTelraw2mdm/USBpixI4ConverterPlugin.cc

@@ -16 +16 @@
 #  include "IMPL/LCGenericObjectImpl.h"
 #  include "UTIL/CellIDEncoder.h"
+
+#include "eudaq/write_mdm.h"
+
 #endif
 
@@ -30 +33 @@
 
 namespace eudaq {
-  // The event type for which this converter plugin will be registered
-  // Modify this to match your actual event type (from the Producer)
-  static const std::string EVENT_TYPE = "USBPIXI4";
-
-  static const unsigned int CHIP_MIN_COL = 1;
-  static const unsigned int CHIP_MAX_COL = 80;
-  static const unsigned int CHIP_MIN_ROW = 1;
-  static const unsigned int CHIP_MAX_ROW = 336;
-  static const unsigned int CHIP_MAX_ROW_NORM = CHIP_MAX_ROW - CHIP_MIN_ROW;    // Maximum ROW normalized (starting with 0)
-  static const unsigned int CHIP_MAX_COL_NORM = CHIP_MAX_COL - CHIP_MIN_COL;
-
-        static int chip_id_offset = 20;
-
-  class USBpixI4ConverterBase {
-  private:
-        unsigned int count_boards;
+// The event type for which this converter plugin will be registered
+// Modify this to match your actual event type (from the Producer)
+static const std::string EVENT_TYPE = "USBPIXI4";
+
+static const unsigned int CHIP_MIN_COL = 1;
+static const unsigned int CHIP_MAX_COL = 80;
+static const unsigned int CHIP_MIN_ROW = 1;
+static const unsigned int CHIP_MAX_ROW = 336;
+static const unsigned int CHIP_MAX_ROW_NORM = CHIP_MAX_ROW - CHIP_MIN_ROW;      // Maximum ROW normalized (starting with 0)
+static const unsigned int CHIP_MAX_COL_NORM = CHIP_MAX_COL - CHIP_MIN_COL;
+
+static int chip_id_offset = 20;
+
+class USBpixI4ConverterBase {
+private:
+        unsigned int count_boards;
         std::vector<unsigned int> board_ids;
 
-  public:
-        unsigned int consecutive_lvl1;
+public:
+
+        mdm * m_mdmhandler;
+
+        USBpixI4ConverterBase () {
+
+                // The constructor will be called on the first event since I need to know the number of DUT's
+                m_mdmhandler = new mdm();
+
+        }
+
+        ~USBpixI4ConverterBase() {
+                int nDet = m_mdmhandler->getndet();
+                for(int detItr = 0 ; detItr < nDet ; detItr++) {
+                        WriteToNtuple::GetInstance("", "", nDet, m_mdmhandler->getDetIdToIndexArray()[detItr], m_mdmhandler->getDetIdToIndexArray())->closeNtuple();
+                }
+        }
+
+        unsigned int consecutive_lvl1;
         unsigned int tot_mode;
         int first_sensor_id;
 
-  int getCountBoards() {
+        int getCountBoards() {
                 return count_boards;
         }
@@ -66 +86 @@
         }
 
-  void getBOREparameters (const Event & ev) {
+        void getBOREparameters (const Event & ev) {
                 count_boards = ev.GetTag ("boards", -1);
                 board_ids.clear();
@@ -201 +221 @@
                 return plane;
         }
-  };
-
-  // Declare a new class that inherits from DataConverterPlugin
-  class USBPixI4ConverterPlugin : public DataConverterPlugin , public USBpixI4ConverterBase {
-
-  public:
-
-    // This is called once at the beginning of each run.
-    // You may extract information from the BORE and/or configuration
-    // and store it in member variables to use during the decoding later.
-    virtual void Initialize(const Event & bore,
-                            const Configuration & cnf) {
-      (void)cnf; // just to suppress a warning about unused parameter cnf
-      getBOREparameters (bore);
-    }
-
-    // This should return the trigger ID (as provided by the TLU)
-    // if it was read out, otherwise it can either return (unsigned)-1,
-    // or be left undefined as there is already a default version.
-    virtual unsigned GetTriggerID(const Event & ev) const {
+};
+
+// Declare a new class that inherits from DataConverterPlugin
+class USBPixI4ConverterPlugin : public DataConverterPlugin , public USBpixI4ConverterBase {
+
+public:
+
+        // This is called once at the beginning of each run.
+        // You may extract information from the BORE and/or configuration
+        // and store it in member variables to use during the decoding later.
+        virtual void Initialize(const Event & bore,
+                        const Configuration & cnf) {
+                (void)cnf; // just to suppress a warning about unused parameter cnf
+                getBOREparameters (bore);
+        }
+
+        // This should return the trigger ID (as provided by the TLU)
+        // if it was read out, otherwise it can either return (unsigned)-1,
+        // or be left undefined as there is already a default version.
+        virtual unsigned GetTriggerID(const Event & ev) const {
 
                 // Make sure the event is of class RawDataEvent
@@ -236 +256 @@
         }
 
-    // Here, the data from the RawDataEvent is extracted into a StandardEvent.
-    // The return value indicates whether the conversion was successful.
-    // Again, this is just an example, adapted it for the actual data layout.
-    virtual bool GetStandardSubEvent(StandardEvent & sev, const Event & ev) const {
+        // Here, the data from the RawDataEvent is extracted into a StandardEvent.
+        // The return value indicates whether the conversion was successful.
+        // Again, this is just an example, adapted it for the actual data layout.
+        virtual bool GetStandardSubEvent(StandardEvent & sev, const Event & ev) const {
                 if (ev.IsBORE()) {
                         return true;
@@ -253 +273 @@
                 }
 
-        return true;
-    }
-
-#if USE_LCIO && USE_EUTELESCOPE
-    // This is where the conversion to LCIO is done
-    virtual lcio::LCEvent * GetLCIOEvent(const Event * /*ev*/) const {
+                return true;
+        }
+
+        //i- #if USE_LCIO && USE_EUTELESCOPE
+        // This is where the conversion to LCIO is done
+        virtual lcio::LCEvent * GetLCIOEvent(const Event * /*ev*/) const {
                 return 0;
-    }
+        }
 
         virtual bool GetLCIOSubEvent(lcio::LCEvent & lcioEvent, const Event & eudaqEvent) const {
-                //std::cout << "getlciosubevent (I4) event " << eudaqEvent.GetEventNumber() << " | " << GetTriggerID(eudaqEvent) << std::endl;
+
+                std::cout << "getlciosubevent (I4) event " << eudaqEvent.GetEventNumber() << " | " << GetTriggerID(eudaqEvent) << std::endl;
+
                 if (eudaqEvent.IsBORE()) {
                         // shouldn't happen
@@ -293 +315 @@
                 const RawDataEvent & ev_raw = dynamic_cast <const RawDataEvent &> (eudaqEvent);
 
+                // First time. Initialize the mdm handlers, number of DUTs won't change in a run.
+                if (!m_mdmhandler->isInitialized()) {
+
+                        // find sensor ID's
+                        std::vector<int> sIDs;
+                        for (size_t chip = 0; chip < ev_raw.NumBlocks(); ++chip) {
+                                // as calculated and stored later in zsDataEncoder["sensorID"]
+                                sIDs.push_back(ev_raw.GetID(chip) + chip_id_offset + first_sensor_id);
+                        }
+
+                        // Initialize the Handler !
+                        m_mdmhandler->init(sIDs, EVENT_TYPE, CHIP_MAX_COL, CHIP_MAX_ROW);
+
+                }
+
                 std::vector< eutelescope::EUTelSetupDescription * >  setupDescription;
 
@@ -320 +357 @@
                         // a single TrackerData object (zsFrame) that will correspond to a single sensor in one event
                         std::auto_ptr< eutelescope::EUTelSparseDataImpl< eutelescope::EUTelAPIXSparsePixel > >
-                                sparseFrame( new eutelescope::EUTelSparseDataImpl< eutelescope::EUTelAPIXSparsePixel > ( zsFrame.get() ) );
+                        sparseFrame( new eutelescope::EUTelSparseDataImpl< eutelescope::EUTelAPIXSparsePixel > ( zsFrame.get() ) );
 
                         unsigned int ToT = 0;
@@ -330 +367 @@
                                 unsigned int Word = (((unsigned int)buffer[i + 3]) << 24) | (((unsigned int)buffer[i + 2]) << 16) | (((unsigned int)buffer[i + 1]) << 8) | (unsigned int)buffer[i];
 
+
                                 if (DATA_HEADER_MACRO(Word)) {
                                         lvl1++;
                                 } else {
+
                                         // First Hit
                                         if (getHitData(Word, false, Col, Row, ToT)) {
@@ -338 +377 @@
                                                 sparseFrame->addSparsePixel( thisHit );
                                                 tmphits.push_back( thisHit );
+
+                                                //std::cout << ".... f -> " << i << " : " << Col << "," << Row << "," << ToT << "," << lvl1-1 << std::endl;
+                                                // hit to mdm
+                                                m_mdmhandler->dofill(zsDataEncoder["sensorID"], Col, Row, ToT, lvl1-1);
                                         }
                                         // Second Hit
@@ -344 +387 @@
                                                 sparseFrame->addSparsePixel( thisHit );
                                                 tmphits.push_back( thisHit );
+
+                                                //std::cout << ".... s -> " << i << " : " << Col << "," << Row << "," << ToT << "," << lvl1-1 << std::endl;
+                                                // hit to mdm
+                                                m_mdmhandler->dofill(zsDataEncoder["sensorID"], Col, Row, ToT, lvl1-1);
                                         }
                                 }
                         }
 
+                        // Write to mdm !  this member clears the matrix after filling
+                        //std::cout << "chip : " << chip << std::endl;
+                        m_mdmhandler->writetontuple(zsDataEncoder["sensorID"], lcioEvent.getEventNumber());
+
                         // write TrackerData object that contains info from one sensor to LCIO collection
                         zsDataCollection->push_back( zsFrame.release() );
@@ -353 +404 @@
                         // clean up
                         for( std::list<eutelescope::EUTelAPIXSparsePixel*>::iterator it = tmphits.begin(); it != tmphits.end(); it++ ){
-                          delete (*it);
+                                delete (*it);
                         }
                 }
@@ -380 +431 @@
                 return true;
 
-    }
-#endif
-
-  private:
-
-    // The constructor can be private, only one static instance is created
-    // The DataConverterPlugin constructor must be passed the event type
-    // in order to register this converter for the corresponding conversions
-    // Member variables should also be initialized to default values here.
-    USBPixI4ConverterPlugin()
-      : DataConverterPlugin(EVENT_TYPE)
-      {}
-
-    // The single instance of this converter plugin
-    static USBPixI4ConverterPlugin m_instance;
-  };
-
-  // Instantiate the converter plugin instance
-  USBPixI4ConverterPlugin USBPixI4ConverterPlugin::m_instance;
+        }
+        //#endif
+
+private:
+
+        // The constructor can be private, only one static instance is created
+        // The DataConverterPlugin constructor must be passed the event type
+        // in order to register this converter for the corresponding conversions
+        // Member variables should also be initialized to default values here.
+        USBPixI4ConverterPlugin()
+        : DataConverterPlugin(EVENT_TYPE)
+        {}
+
+        // The single instance of this converter plugin
+        static USBPixI4ConverterPlugin m_instance;
+};
+
+// Instantiate the converter plugin instance
+USBPixI4ConverterPlugin USBPixI4ConverterPlugin::m_instance;
+
 
 } // namespace eudaq