diff --git a/Projects/ComptonTelescope/online/src/CMakeLists.txt b/Projects/ComptonTelescope/online/src/CMakeLists.txt
index 6c94053122fe46778b0e78d3ce94d37bb70778df..7dc613678da0daf19736d960f437c19b10d7111f 100644
--- a/Projects/ComptonTelescope/online/src/CMakeLists.txt
+++ b/Projects/ComptonTelescope/online/src/CMakeLists.txt
@@ -8,7 +8,9 @@ add_library(DecodeD SHARED DecodeD.cpp)
add_library(DecodeT SHARED DecodeT.cpp)
# add the executable
-add_executable(online online.cpp)
+add_executable(online_coinc online_coinc.cpp)
+add_executable(online_dsssd online_dsssd.cpp)
# link library rules
-target_link_libraries(online DecodeR DecodeD DecodeT ${ROOT_LIBRARIES} NPCore NPPhysics NPComptonTelescope)
+target_link_libraries(online_coinc DecodeR DecodeD DecodeT ${ROOT_LIBRARIES} NPCore NPPhysics NPComptonTelescope)
+target_link_libraries(online_dsssd DecodeR DecodeD DecodeT ${ROOT_LIBRARIES} NPCore NPPhysics NPComptonTelescope)
diff --git a/Projects/ComptonTelescope/online/src/online.cpp b/Projects/ComptonTelescope/online/src/online_coinc.cpp
similarity index 100%
rename from Projects/ComptonTelescope/online/src/online.cpp
rename to Projects/ComptonTelescope/online/src/online_coinc.cpp
diff --git a/Projects/ComptonTelescope/online/src/online_dsssd.cpp b/Projects/ComptonTelescope/online/src/online_dsssd.cpp
new file mode 100644
index 0000000000000000000000000000000000000000..c7eb167af505a38df94a41f38d31b0373dd9950c
--- /dev/null
+++ b/Projects/ComptonTelescope/online/src/online_dsssd.cpp
@@ -0,0 +1,520 @@
+// nptool headers
+#include "NPOptionManager.h"
+#include "RootOutput.h"
+#include "NPDetectorManager.h"
+#include "TComptonTelescopeData.h"
+#include "TComptonTelescopePhysics.h"
+
+// root headers
+#include "TCutG.h"
+//#include "TH2.h"
+#include "TFile.h"
+#include "TCanvas.h"
+
+// custom headers
+#include "DecodeR.h"
+#include "DecodeD.h"
+#include "DecodeT.h"
+
+#define __TEST_ZONE__
+#undef __TEST_ZONE__
+#define __USE_CUTG__
+#undef __USE_CUTG__
+
+// C++ headers
+#include <iostream>
+#include <fstream>
+#include <string>
+using namespace std;
+
+////--//--// One-line setter for calorimeter //--//--//
+//
+//void setCTCalorimeter(TComptonTelescopeData* ccamData, DecodeR* D, const int pixelNumber)
+//{
+// ccamData -> SetCTCalorimeter(1, 1, D->getPixelNumber(), D->getTime(), D->getData(), pixelNumber);
+///* ccamData -> SetCTCalorimeterTTowerNbr( 1 );
+// ccamData -> SetCTCalorimeterTDetectorNbr( 1 );//Triggered ASIC number
+// ccamData -> SetCTCalorimeterTChannelNbr( D -> getPixelNumber() );//Pixel that triggered
+// ccamData -> SetCTCalorimeterTTime( D -> getTime() );
+// for (int i = 0; i < pixelNumber; ++i) {//Loop on pixels
+// ccamData -> SetCTCalorimeterETowerNbr(1);
+// ccamData -> SetCTCalorimeterEDetectorNbr( 1 );
+// ccamData -> SetCTCalorimeterEChannelNbr( i );//PMTÂ pixel number
+// ccamData -> SetCTCalorimeterEEnergy( D -> getData()[i] );
+// }//End of loop on pixels*/
+//}
+
+//--//--// One-line setter for DSSSD(s) //--//--//
+
+// One-line setter for the Front of one DSSD
+void setCTTrackerFront(TComptonTelescopeData* ccamData, newframe_t* event, int detNbr, int faceNbr, const int stripNumber)
+{
+ //detNbr and faceNbr are > 0
+ ccamData -> SetCTTrackerFrontTTowerNbr(1);
+ ccamData -> SetCTTrackerFrontTDetectorNbr(detNbr);
+ ccamData -> SetCTTrackerFrontTStripNbr(33);
+ ccamData -> SetCTTrackerFrontTTime(event->timestamp);
+ for (int k = 0; k < stripNumber; k++) {//Loop on strips
+ ccamData -> SetCTTrackerFrontETowerNbr(1);
+ ccamData -> SetCTTrackerFrontEDetectorNbr(detNbr);
+ ccamData -> SetCTTrackerFrontEStripNbr(k);
+ ccamData -> SetCTTrackerFrontEEnergy(event->sample[faceNbr-1][detNbr-1][k]);
+ }//End of loop on strips
+}
+
+// One-line setter for the Back of one DSSD
+void setCTTrackerBack(TComptonTelescopeData* ccamData, newframe_t* event, int detNbr, int faceNbr, const int stripNumber)
+{
+ ccamData -> SetCTTrackerBackTTowerNbr(1);
+ ccamData -> SetCTTrackerBackTDetectorNbr(detNbr);
+ ccamData -> SetCTTrackerBackTStripNbr(33);
+ ccamData -> SetCTTrackerBackTTime(event->timestamp);
+ for (int k = 0; k < stripNumber; k++) {//Loop on strips
+ ccamData -> SetCTTrackerBackETowerNbr(1);
+ ccamData -> SetCTTrackerBackEDetectorNbr(detNbr);
+ ccamData -> SetCTTrackerBackEStripNbr(k);
+ ccamData -> SetCTTrackerBackEEnergy(event->sample[faceNbr-1][detNbr-1][k]);
+ }//End of loop on strips
+}
+
+// One-line setter for DSSSD(s)
+void setCTTracker(TComptonTelescopeData* ccamData, newframe_t* event, vector<int>* nb_asic, vector<int>* chain, const int stripNumber)
+{
+ for (vector<int>::iterator itchain = chain->begin(); itchain != chain->end(); ++itchain) {//Iterates on 2 faces
+ for (vector<int>::iterator itasic = nb_asic->begin(); itasic != nb_asic->end(); ++itasic) {//Iterates on 1 DSSSD
+ if (event->chip_data[*itchain][*itasic]) { // Test if data is present
+ switch (*itchain) {
+ case 0://Assuming 0 is front - to be checked
+ setCTTrackerFront(ccamData, event, *itasic+1, *itchain+1, stripNumber);
+ break;
+ case 1://Assuming 1 is back - to be checked
+ setCTTrackerBack(ccamData, event, *itasic+1, *itchain+1, stripNumber);
+ }//End switch
+ }//End if
+ }//End for
+ }//End for
+}
+
+//--//--//--//--//--//--//--//--//--//--//--//--//--//--//--//
+int main()
+{
+// auto fout = new TFile("pipo.root", "recreate");
+// auto bidim = new TH2F("bidim", "bidim", 2001, -1000, 1000, 2001, -1000, 1000);
+#ifdef __USE_CUTG__
+/* TFile* fcut = new TFile("/disk/proto-data/data/coinc-si/CUT_Compton.root");
+ TCutG* mcut = (TCutG*) fcut -> Get("CUT_Compton");
+ fcut -> Close();
+ cout << fcut << endl;
+ cout << mcut << endl;
+ TCanvas* can = new TCanvas();
+ can->cd();
+ mcut->Draw();*/
+ TCutG *mcut = new TCutG("CUT_Compton",7);
+ mcut->SetVarX("Calor_E");
+ mcut->SetVarY("Half_Energy");
+ mcut->SetTitle("Graph");
+ mcut->SetFillStyle(1000);
+ mcut->SetPoint(0,526905.4,0.6145632);
+ mcut->SetPoint(1,526905.4,0.4996379);
+ mcut->SetPoint(2,592808,0.1947341);
+ mcut->SetPoint(3,629341,0.02586437);
+ mcut->SetPoint(4,627908.3,0.211152);
+ mcut->SetPoint(5,540157.6,0.5981453);
+ mcut->SetPoint(6,526905.4,0.6145632);
+#endif
+
+ ///////////////////////////////////////////////////////////////////////////
+ // configure option manager
+// NPOptionManager::getInstance()->Destroy();
+#ifdef __TEST_ZONE__
+// string arg = "-D ./ComptonCAM.detector -C calibrations.txt -GH -E ./10He.reaction --circular";
+#else
+ string arg = "-D ./ComptonCAM.detector -C calibrations.txt -GH -E ./10He.reaction";
+#endif
+ NPOptionManager::getInstance(arg);
+
+ // open ROOT output file
+ RootOutput::getInstance("OnlineTree.root", "OnlineTree");
+ // get tree pointer
+ auto m_OutputTree = RootOutput::getInstance()->GetTree();
+
+ // configure detector manager
+ string detectorfileName = NPOptionManager::getInstance()->GetDetectorFile();
+ NPL::DetectorManager* m_NPDetectorManager = new NPL::DetectorManager();
+ m_NPDetectorManager->ReadConfigurationFile(detectorfileName);
+ m_NPDetectorManager->InitializeRootOutput();
+
+ // configure spectra server
+ m_NPDetectorManager->SetSpectraServer();
+
+ // instantiate raw ComptonCAM data pointer
+ auto ccamData = new TComptonTelescopeData();
+ // connect raw CCAM data pointer to physics class
+ auto ccamPhys = (TComptonTelescopePhysics*) m_NPDetectorManager->GetDetector("ComptonTelescope");
+ ccamPhys->SetRawDataPointer(ccamData);
+
+ // read data file/flux and fill ccamData object
+ std::cout << "Reading data\n";
+ DecodeR* DR = new DecodeR(false); // Instantiates DecodeR object reading calorimeter data flux
+ DecodeT* DT = new DecodeT(false); // Instantiates DecodeT object reading trigger data flux
+ DecodeD* DD = new DecodeD(false); // Instantiates DecodeD object reading DSSSD(s) data flux
+ newframe_t* event;
+ DD -> setTree("/disk/proto-data/data/coinc-si/bb7_3309-7_cs137-20210205_11h41_coinc_run0_conv.root");
+ int dlen = DD -> getLength();
+
+ //Sets where to look for data in DSSSD root frames
+ vector <int> chain ({0, 1});//Two faces
+ vector <int> nb_asic ({0});//First DSSSD
+
+ int i = 0;// ROSMAP files loop counter
+ int c = 0;// Event counter
+ int cc = 0;
+ // Set some constants
+ const int pixelNumber = 64;
+ const int stripNumber = 32;
+
+#ifdef __TEST_ZONE__
+
+ while (DD -> getCursor() < dlen)
+ {
+ DD -> decodeEvent();
+
+ // Clear raw and physics data
+ m_NPDetectorManager->ClearEventPhysics();
+ m_NPDetectorManager->ClearEventData();
+
+ // Fill data
+ setCTTracker(ccamData, DD -> getEvent(), &nb_asic, &chain, stripNumber);
+
+ // Build physical event
+ m_NPDetectorManager->BuildPhysicalEvent();
+
+ // Fill object in output ROOT file
+ m_OutputTree->Fill();
+
+ // check spectra
+ m_NPDetectorManager->CheckSpectraServer();
+ }
+
+#else
+ // Open data files
+ ifstream iros, itrig;
+ cout << "Loading data files ";
+
+ itrig.open("/disk/proto-data/data/coinc-si/mfm_trigger_202102051141.raw", ios::binary);
+ itrig.seekg(0, ios::end);
+ int tlen = itrig.tellg();
+ itrig.seekg(0, ios::beg);
+ char* tbuff = new char[tlen];
+ itrig.read(tbuff, tlen);
+ itrig.close();
+ cout << "... ";
+
+ iros.open("/disk/proto-data/data/coinc-si/mfm_rdd_rosmap_04_mfm_rosmap_04_2021-02-05_10_41_46.raw.0001", ios::binary);
+ iros.seekg(0, ios::end);
+ int rlen = iros.tellg();
+ iros.seekg(0, ios::beg);
+ char* rbuff = new char[rlen];
+ iros.read(rbuff, rlen);
+ iros.close();
+ cout << "Done" << endl;
+
+// for (int reset=-1000; reset<1001; reset++)
+// {
+// cout << reset << endl;
+ DT -> setRaw(tbuff);
+ DR -> setRaw(rbuff);
+ //DD -> rewind();
+
+ DR -> decodeBlobMFM();
+ DD -> decodeEvent();
+
+ int cr = -17;
+ int cd = 0;
+ c = 0;
+ i = 0;
+ int tr = DR -> getTime();
+ int td = DD -> getTime();
+// int dt = 100;
+ while(DR -> getCursor() < rlen and DD -> getCursor() < dlen)
+ {
+// cout << DR -> getTime() << " " << DD -> getTime() << endl;
+ if (cr == cd) {
+ //if (abs(tr-td+55) < dt) {
+// if (td-tr < 150 and tr-td < 10) { // That one is the real one
+ if (tr-td > 50 and tr-td < 1000) { // Same width, bad position
+ c++;
+ cout << cc << " " << c << "(" << cr << ", " << cd << ") : " << tr << " " << td << endl;
+ //DR -> Dump();
+ //DD -> Dump();
+ //bidim->Fill(reset, td-tr);
+
+ // Clear raw and physics data
+ m_NPDetectorManager->ClearEventPhysics();
+ m_NPDetectorManager->ClearEventData();
+
+ // Fill data
+ ccamData -> SetCTCalorimeter(1, 4, DR->getPixelNumber(), DR->getTime(), DR->getData(), pixelNumber);
+ setCTTracker(ccamData, DD -> getEvent(), &nb_asic, &chain, stripNumber);
+
+ // Build physical event
+ m_NPDetectorManager->BuildPhysicalEvent();
+
+ // Fill object in output ROOT file
+ if (ccamPhys->EventMultiplicity > 0) {
+#ifdef __USE_CUTG__
+ if (mcut->IsInside(ccamPhys->Calor_E[0], ccamPhys->Half_Energy[0])) {
+ //cout << "c" << endl;
+ cc++;
+ m_OutputTree->Fill();
+ }
+#else
+ m_OutputTree->Fill();
+ cc++;
+#endif
+ }
+ //cout << "c" << endl;
+
+ // check spectra
+ m_NPDetectorManager->CheckSpectraServer();
+ //cout << "d" << endl;
+ }
+ if (td < tr) {
+ DD -> decodeEvent();
+ } else {
+ DR -> decodeBlobMFM();
+ }
+ } else if (cr < cd) {
+ DR -> decodeBlobMFM();
+ } else {
+ DD -> decodeEvent();
+ }
+ if (DR -> getTime() < tr) {
+ cr++;
+ }
+ tr = DR -> getTime();
+ if (DD -> getTime() < td) {
+ cd ++;
+ }
+ td = DD -> getTime();
+
+ i++;
+
+ } // End of main loop
+// }
+
+ delete DR;
+ delete DT;
+ delete DD;
+ delete [] rbuff;
+ delete [] tbuff;
+#endif
+
+ // fill spectra
+ m_NPDetectorManager -> WriteSpectra();
+
+// fout->cd();
+// bidim->Write();
+// fout->Close();
+
+ // Essential
+#if __cplusplus > 199711L && NPMULTITHREADING
+ m_NPDetectorManager->StopThread();
+#endif
+ RootOutput::Destroy();
+
+ return 0;
+}
+
+
+/*
+ const bool loopForever = false;
+ //while (loopForever or i<12) // for Am data
+ while (loopForever or i<3) // for Bi data quick analysis
+ //while (loopForever or i<489) // for Bi data all events
+ {
+ // Load a root file and setup DecodeD
+ //DD -> setTree("../data/20200128_11h58_am241_conv.root");
+ DD -> setTree("../data/bb7_3309-7_bi207_20210203_16h25_run8_conv.root");
+ int dlength = DD -> getLength();
+
+ //while (DD -> getCursor() < dlength and (loopForever or i<12))
+ while (DD -> getCursor() < dlength and (loopForever or i<3))
+ {
+ // Load a ROSMAP file
+ std::ifstream is;
+// i = 1;
+ switch (i % 3) {
+ case 3: is.open("./mfm.bin", std::ios::binary); break;
+ case 4: is.open("./133Ba.bin", std::ios::binary); break;
+ case 5: is.open("./241Am.bin", std::ios::binary); break;
+ case 0: is.open("./207Bi.bin", std::ios::binary); break;
+ case 1: is.open("./241Am-1.bin", std::ios::binary); break;
+ case 2: is.open("./241Am-2.bin", std::ios::binary); break;
+ }
+ is.seekg (0, std::ios::end);
+ int length = is.tellg();
+ is.seekg (0, ios::beg);
+ char* buffer = new char [length];
+ is.read(buffer, length);
+ is.close();
+ i++;
+
+ // Setup DecodeR
+ DR -> setRaw(buffer);
+ DR -> decodeRawMFM(); // get rid of the first two (empty) events
+ DR -> decodeRawMFM();
+
+ // Loop on some events
+ while (DD -> getCursor() < dlength and DR -> getCursor() < length)
+ {
+ // Clear raw and physics data
+ m_NPDetectorManager->ClearEventPhysics();
+ m_NPDetectorManager->ClearEventData();
+
+ // Fill calorimeter data
+ DR -> decodeRawMFM();
+ setCTCalorimeter(ccamData, DR, pixelNumber);
+
+ // Fill DSSD data
+ DD -> decodeEvent();
+ event = DD -> getEvent();
+ setCTTracker(ccamData, event, &nb_asic, &chain, stripNumber);
+
+ // Build physical event
+ m_NPDetectorManager->BuildPhysicalEvent();
+
+ // Fill object in output ROOT file
+ m_OutputTree->Fill();
+
+ // check spectra
+ m_NPDetectorManager->CheckSpectraServer();
+
+ c++;
+ //usleep(100);//Simulated 10kHz count rate
+
+ }// End of loop on ROSMAP events
+
+ delete [] buffer;
+
+ }// End of loop on DSSSD data
+
+ }// End of main loop
+
+ // fill spectra
+ m_NPDetectorManager->WriteSpectra();
+
+ // delete Decoders
+ delete DR;
+ delete DD;
+
+ // Essential
+ #if __cplusplus > 199711L && NPMULTITHREADING
+ m_NPDetectorManager->StopThread();
+ #endif
+ RootOutput::Destroy();
+
+ return 0;
+}
+
+*/
+
+
+
+
+/*
+ while (DD -> getCursor() < dlength)
+ {
+ // Clear raw and physics data
+ m_NPDetectorManager->ClearEventPhysics();
+ m_NPDetectorManager->ClearEventData();
+
+ //Read some data
+ DD -> decodeEvent();
+ event = DD -> getEvent();
+ //Fill TComptonTelescopeData here (if possible)
+ for (vector<int>::iterator itchain = chain.begin(); itchain != chain.end(); ++itchain) {//Iterates on 2 faces
+ for (vector<int>::iterator itasic = nb_asic.begin(); itasic != nb_asic.end(); ++itasic) {//Iterates on 1 DSSSD
+ if (event->chip_data[*itchain][*itasic]) { // Test if data is present
+ switch (*itchain) {
+ case 0://Assuming 0 is front - to be checked
+ ccamData -> SetCTTrackerFrontTTowerNbr(1);
+ ccamData -> SetCTTrackerFrontTDetectorNbr(*itasic+1);
+ ccamData -> SetCTTrackerFrontTStripNbr(33);
+ ccamData -> SetCTTrackerFrontTTime(event->timestamp);
+ for (int k = 0; k < stripNumber; k++) {//Loop on strips
+ ccamData -> SetCTTrackerFrontETowerNbr(1);
+ ccamData -> SetCTTrackerFrontEDetectorNbr(*itasic+1);
+ ccamData -> SetCTTrackerFrontEStripNbr(k+1);
+ ccamData -> SetCTTrackerFrontEEnergy(event->sample[*itchain][*itasic][k]);
+ }//End of loop on strips
+ break;
+ case 1://Assuming 1 is back - to be checked
+ ccamData -> SetCTTrackerBackTTowerNbr(1);
+ ccamData -> SetCTTrackerBackTDetectorNbr(*itasic+1);
+ ccamData -> SetCTTrackerBackTStripNbr(33);
+ ccamData -> SetCTTrackerBackTTime(event->timestamp);
+ for (int k = 0; k < stripNumber; k++) {//Loop on strips
+ ccamData -> SetCTTrackerBackETowerNbr(1);
+ ccamData -> SetCTTrackerBackEDetectorNbr(*itasic+1);
+ ccamData -> SetCTTrackerBackEStripNbr(k+1);
+ ccamData -> SetCTTrackerBackEEnergy(event->sample[*itchain][*itasic][k]);
+ }//End of loop on strips
+ }//End switch
+ }//End if
+ }//End for
+ }//End for
+
+ // Build physical event
+ m_NPDetectorManager->BuildPhysicalEvent();
+ //Fill output tree
+ m_OutputTree->Fill();
+ //check spectra
+ m_NPDetectorManager->CheckSpectraServer();
+
+ c++;
+
+ }//End while
+ cout << "Read " << c << " DSSSD events from file" << endl;
+
+ // Read from file
+ D -> setRaw(buffer);
+
+ D -> decodeRawMFM(); // get rid of the first two (empty) events
+ D -> decodeRawMFM();
+
+ while (D -> getCursor() < length)
+ {
+ // Clear raw data and physics objects
+ m_NPDetectorManager->ClearEventPhysics();
+ m_NPDetectorManager->ClearEventData();
+
+ // Read the actual data
+ D -> decodeRawMFM();
+ //D -> Dump();//Optionnal print
+
+ // Set ccamData (a better way is envisionned)
+ for (int i = 0; i < pixelNumber; ++i) {
+ ccamData -> SetCTCalorimeterETowerNbr(1);
+ ccamData -> SetCTCalorimeterEDetectorNbr( 1 );
+ ccamData -> SetCTCalorimeterEChannelNbr( i );//PMTÂ pixel number
+ ccamData -> SetCTCalorimeterEEnergy( D -> getData()[i] );
+ }
+ ccamData -> SetCTCalorimeterTTowerNbr( 1 );
+ ccamData -> SetCTCalorimeterTDetectorNbr( 1 );//Triggered ASIC number
+ ccamData -> SetCTCalorimeterTChannelNbr( D -> getPixelNumber() );//Pixel that triggered
+ ccamData -> SetCTCalorimeterTTime( D -> getTime() );
+ // ccamData -> Dump();
+
+ // Build physical event
+ m_NPDetectorManager->BuildPhysicalEvent();
+
+ // Fill object in output ROOT file
+ m_OutputTree->Fill();
+
+ // check spectra
+ m_NPDetectorManager->CheckSpectraServer();
+
+ c++;
+ usleep(100);//Simulated 10kHz count rate
+ }*/