--- /dev/null
+// this is example for
+
+
+void first()
+{
+ // base::ProcMgr::instance()->SetRawAnalysis(true);
+ base::ProcMgr::instance()->SetTriggeredAnalysis(true);
+
+ // all new instances get this value
+ base::ProcMgr::instance()->SetHistFilling(4);
+
+ // configure bubbles
+ //hadaq::TdcProcessor::SetBubbleMode(0, 18);
+
+ // this limits used for liner calibrations when nothing else is available
+ hadaq::TdcMessage::SetFineLimits(20, 491);
+
+ // default channel numbers and edges mask
+ hadaq::TrbProcessor::SetDefaults(33, 2);
+
+ // [min..max] range for TDC ids
+ hadaq::TrbProcessor::SetTDCRange(0x0800, 0x12FF);
+
+ // [min..max] range for HUB ids
+ hadaq::TrbProcessor::SetHUBRange(0x8000, 0xc002);
+
+ // when first argument true - TRB/TDC will be created on-the-fly
+ // second parameter is function name, called after elements are created
+ hadaq::HldProcessor* hld = new hadaq::HldProcessor(true, "after_create");
+
+ const char* calname = getenv("CALNAME");
+ if ((calname==0) || (*calname==0)) calname = "calib_";
+ const char* calmode = getenv("CALMODE");
+ int cnt = (calmode && *calmode) ? atoi(calmode) : 10000;
+ //#cnt=10000;
+ const char* caltrig = getenv("CALTRIG");
+ unsigned trig = (caltrig && *caltrig) ? atoi(caltrig) : 0xd;
+ const char* uset = getenv("USETEMP");
+ unsigned use_temp = 0; // 0x80000000;
+ if ((uset!=0) && (*uset!=0) && (strcmp(uset,"1")==0)) use_temp = 0x80000000;
+
+ printf("TDC CALIBRATION MODE: %d, cal trigger: %d\n", cnt, trig);
+
+ //printf("HLD configure calibration calfile:%s cnt:%d trig:%X temp:%X\n", calname, cnt, trig, use_temp);
+
+ // first parameter if filename prefix for calibration files
+ // and calibration mode (empty string - no file I/O)
+ // second parameter is hits count for autocalibration
+ // 0 - only load calibration
+ // -1 - accumulate data and store calibrations only at the end
+ // >0 - automatic calibration after N hits in each active channel
+ // third parameter is trigger type mask used for calibration
+ // (1 << 0xD) - special 0XD trigger with internal pulser, used also for TOT calibration
+ // 0x3FFF - all kinds of trigger types will be used for calibration (excluding 0xE and 0xF)
+ // 0x80000000 in mask enables usage of temperature correction
+ hld->ConfigureCalibration(calname, cnt, (1 << trig) );
+
+ // only accept trigger type 0x1 when storing file
+ // new hadaq::HldFilter(0x1);
+
+ // create ROOT file store
+ //base::ProcMgr::instance()->CreateStore("td.root");
+
+ // 0 - disable store
+ // 1 - std::vector<hadaq::TdcMessageExt> - includes original TDC message
+ // 2 - std::vector<hadaq::MessageFloat> - compact form, without channel 0, stamp as float (relative to ch0)
+ // 3 - std::vector<hadaq::MessageDouble> - compact form, with channel 0, absolute time stamp as double
+ base::ProcMgr::instance()->SetStoreKind(2);
+
+
+ // when configured as output in DABC, one specifies:
+ // <OutputPort name="Output2" url="stream://file.root?maxsize=5000&kind=3"/>
+ printf("test\n");
+
+}
+
+// extern "C" required by DABC to find function from compiled code
+
+extern "C" void after_create(hadaq::HldProcessor* hld)
+{
+ printf("Called after all sub-components are created\n");
+
+ if (hld==0) return;
+
+ for (unsigned k=0;k<hld->NumberOfTRB();k++) {
+ hadaq::TrbProcessor* trb = hld->GetTRB(k);
+ if (trb==0) continue;
+ printf("Configure %s!\n", trb->GetName());
+ trb->SetPrintErrors(100);
+ }
+
+ for (unsigned k=0;k<hld->NumberOfTDC();k++) {
+ hadaq::TdcProcessor* tdc = hld->GetTDC(k);
+ if (tdc==0) continue;
+
+ printf("Configure %s!\n", tdc->GetName());
+
+ // tdc->SetUseLastHit(true);
+
+ // tdc->SetStoreEnabled();
+ // for (unsigned nch=10;nch<tdc->NumChannels();nch++)
+ // tdc->SetRefChannel(nch, nch-1, 0xffff, 20000, -100., 100.);
+
+
+ tdc->SetRefChannel(8, 6, 0xffff, 20000, -100., 100.);
+ tdc->SetRefChannel(6, 1, 0x0840, 20000, -100., 100.);
+ //tdc->SetRefChannel(10, 8, 0xffff, 20000, -1000., 1000.);
+
+ //tdc->SetRefChannel(1, tdc->NumChannels() -1 , 0xffff, 20000, -100., 100.);
+
+
+ }
+}
+
+
--- /dev/null
+#include <stdio.h>
+
+// #include "TTree.h"
+
+#include "base/EventProc.h"
+#include "base/Event.h"
+#include "hadaq/TdcSubEvent.h"
+
+class SecondProc : public base::EventProc {
+protected:
+
+ std::string fTdcId; //!< tdc id where channels will be selected
+ double fHits[33][2]; //!< 33 channel, abstract hits, two dges
+
+ unsigned eventnumber;
+
+ base::H1handle hNumHits; //!< histogram with hits number
+ base::H1handle hDif1; //!< histogram with hits number
+ base::H1handle hDif2; //!< histogram with hits number
+ base::H1handle hDif3; //!< histogram with hits number
+ base::H1handle hDif4; //!< histogram with hits number
+ base::H1handle hDif5; //!< histogram with hits number
+ base::H1handle hDif6; //!< histogram with hits number
+ base::H1handle hDif7; //!< histogram with hits number
+ base::H1handle hTot1; //!< histogram with hits number
+ base::H1handle hTot2; //!< histogram with hits number
+ base::H1handle hTot3; //!< histogram with hits number
+ base::H1handle hTot4; //!< histogram with hits number
+ base::H1handle hTot5; //!< histogram with hits number
+ base::H1handle hTot6; //!< histogram with hits number
+ base::H2handle hUser; //!< user-defined 2D histogram
+
+public:
+ SecondProc(const char* procname, const char* _tdcid) :
+ base::EventProc(procname),
+ fTdcId(_tdcid),
+ hUser(0)
+ {
+ printf("Create %s for %s\n", GetName(), fTdcId.c_str());
+
+ hNumHits = MakeH1("NumHits","Num hits", 100, 0, 100, "number");
+
+ eventnumber = 0;
+
+ /*
+ hDif1 = MakeH1("ToT1","ToT of channel 17 (18)", 1000, 0.5, 3.0, "ns");
+ hDif2 = MakeH1("LED_diff1","LED diff channel 17 - 21", 1000, -2., -1., "ns");
+ hUser = MakeH2("ToT_vs_LED","ToT versus LED difference", 500, -4, 0, 500, -2, -1, "ToT/ns;LED diff/ns");
+ */
+
+ hDif1 = MakeH1("LE1","8 vs. 6", 2800, -10, 10, "ns");
+ hDif2 = MakeH1("LE2","4 vs. 6", 1800, -10, 10, "ns");
+ hDif3 = MakeH1("LE3","6 vs. 2", 1800, -90, 90, "ns");
+ hDif4 = MakeH1("LETR1","2 vs. 4 trip", 1800, -90, 90, "ns");
+ hDif5 = MakeH1("LETR2","4 vs. 6 trip", 1800, -90, 90, "ns");
+ hDif6 = MakeH1("LETR3","6 vs. 2 trip", 1800, -90, 90, "ns");
+
+ hDif7 = MakeH1("LE_CUT", "8 vs. 6 with tot cut", 2800, -10, 10, "ns");
+
+ hTot1 = MakeH1("ToT1","ToT 6", 1800, -90, 90, "ns");
+ hTot2 = MakeH1("ToT2","ToT 8", 1800, -90, 90, "ns");
+ hTot3 = MakeH1("ToT3","ToT 6", 1800, -90, 90, "ns");
+
+ hTot4 = MakeH1("ToT4","ToT 2, all 3", 1800, -90, 90, "ns");
+ hTot5 = MakeH1("ToT5","ToT 4, all 3", 1800, -90, 90, "ns");
+ hTot6 = MakeH1("ToT6","ToT 6, all 3", 1800, -90, 90, "ns");
+
+
+ /*
+ hDif1 = MakeH1("ToT1","ToT of channel 9", 1000, 32, 40, "ns");
+ hDif2 = MakeH1("LED_diff1","LED diff channel 9 - 11", 1000, -3, -1, "ns");
+ hUser = MakeH2("ToT_vs_LED","ToT versus LED difference", 500, 32, 40, 500, -4, -0, "ToT/ns;LED diff/ns");
+ */
+
+ // enable storing already in constructor
+ // SetStoreEnabled();
+ }
+
+ virtual void CreateBranch(TTree* t)
+ {
+ // only called when tree is created in first.C
+ // one can ignore
+ // t->Branch(GetName(), fHits, "hits[8]/D");
+ }
+
+ virtual bool Process(base::Event* ev)
+ {
+ for (unsigned n=0;n<33;n++) {
+ fHits[n][0] = 0.;
+ fHits[n][1] = 0.;
+ }
+
+ eventnumber++;
+ hadaq::TdcSubEventFloat* sub =
+ dynamic_cast<hadaq::TdcSubEventFloat*> (ev->GetSubEvent(fTdcId));
+
+ //printf("%s process sub %p %s\n", GetName(), sub, fTdcId.c_str());
+
+ // when return false, event processing is cancelled
+ if (sub==0) return true;
+
+ double num(0);
+
+ for (unsigned cnt=0;cnt<sub->Size();cnt++) {
+ const hadaq::MessageFloat& msg = sub->msg(cnt);
+
+ unsigned chid = msg.getCh();
+ unsigned edge = msg.getEdge(); // 0 - rising, 1 - falling
+ // if (chid==0) { ch0tm = ext.GetGlobalTime(); continue; }
+
+ // time relative ch0
+ double tm = msg.stamp;
+
+ // use only first hit in each channel
+ if ((chid<33) && (fHits[chid][edge]==0.)) fHits[chid][edge] = tm;
+
+ //printf(" ch:%3d tm:%f, edge:%d\n", chid, tm, edge);
+ num+=1;
+ }
+
+ //printf("tot(%d): %f %f\n", eventnumber, (fHits[9][0] - fHits[9][1]), (fHits[11][0] - fHits[11][1]));
+ //printf("led(%d): %f %f\n", eventnumber, (fHits[9][0] - fHits[11][0]), (fHits[9][1] - fHits[11][1]));
+
+ FillH1(hNumHits, num);
+
+ /*
+ FillH1(hDif1, (fHits[9][1] - fHits[9][0]) );
+ FillH1(hDif2, (fHits[9][0] - fHits[11][0]) );
+ FillH2(hUser, (fHits[9][1] - fHits[9][0]) , (fHits[9][0] - fHits[11][0]) );
+ */
+
+ //printf("tot: %f\n", (fHits[6][1] - fHits[6][0]));
+ double tot2;
+ double tot4;
+ double tot6;
+ double tot8;
+
+ tot2 = fHits[2][1] - fHits[2][0];
+ tot4 = fHits[4][1] - fHits[4][0];
+ tot6 = fHits[6][1] - fHits[6][0];
+ tot8 = fHits[8][1] - fHits[8][0];
+
+
+ if (tot6 > 2 && tot6 < 200)
+ FillH1(hTot1, tot6);
+
+ if (tot8 > 1 && tot8 < 200)
+ FillH1(hTot2, tot8);
+
+ if (tot6 > 35 && tot6 < 55)
+ FillH1(hTot3, tot6);
+
+
+ if( fHits[6][0]!=0 && fHits[2][0]!=0 && fHits[4][0]!=0) {
+ FillH1(hTot4, tot2);
+ FillH1(hTot5, tot4);
+ FillH1(hTot6, tot6);
+ }
+
+
+
+ if (fHits[8][0]!=0 && fHits[6][0]!=0 )
+ FillH1(hDif1, (fHits[8][0] - fHits[6][0]));
+
+ if (fHits[6][0]!=0 && fHits[4][0]!=0 )
+ FillH1(hDif2, (fHits[6][0] - fHits[4][0]));
+
+ if (fHits[6][0]!=0 && fHits[2][0]!=0 )
+ FillH1(hDif3, (fHits[6][0] - fHits[2][0]));
+
+ if (fHits[4][0]!=0 && fHits[2][0]!=0 && fHits[6][0]!=0 )
+ FillH1(hDif4, (fHits[4][0] - fHits[2][0]));
+
+ if (fHits[6][0]!=0 && fHits[4][0]!=0 && fHits[2][0]!=0 )
+ FillH1(hDif5, (fHits[6][0] - fHits[4][0]));
+
+ if (fHits[6][0]!=0 && fHits[2][0]!=0 && fHits[4][0]!=0)
+ FillH1(hDif6, (fHits[6][0] - fHits[2][0]));
+
+
+ if (fHits[8][0]!=0 && fHits[6][0]!=0 &&
+
+ tot6 > 26.1 && tot6 < 26.15 &&
+ tot8 > 25.5 && tot8 < 25.55
+
+ /*
+ tot2 > 7. && tot2 < 12. &&
+ tot4 > 8. && tot4 < 13.0 &&
+ tot6 > 46.0 && tot6 < 50.0
+ */
+
+
+
+
+ )
+ FillH1(hDif7, (fHits[8][0] - fHits[6][0]) );
+
+
+
+
+ // cuts no longer needed - one see only normal triggers here
+ /* if(
+ ((fHits[6][1] - fHits[6][0]) > 38) && ((fHits[6][1] - fHits[6][0])<48) &&
+ ((fHits[4][1] - fHits[4][0]) > 35) && ((fHits[4][1] - fHits[4][0])<42) &&
+ ((fHits[2][1] - fHits[2][0]) > 35) && ((fHits[2][1] - fHits[2][0])<42)
+ ) {
+ //FillH1(hDif1, (fHits[6][0] + fHits[4][0])/2 - fHits[2][0] );
+ FillH1(hDif1, (fHits[6][0] - fHits[2][0] ));
+ //FillH1(hDif2, (fHits[1][1] - fHits[2][1]) );
+ }
+ */
+
+
+ return true;
+ }
+};
+
+
+void second()
+{
+ //new SecondProc("A", "TDC_1133");
+ //new SecondProc("A", "TDC_1580");
+ new SecondProc("A", "TDC_1225");
+}
jspc29 / daqtools.git / commitdiff