-void AliTOFRawDataRead(Int_t iEvent=0)
+// gSystem->SetIncludePath("-I$ROOTSYS/include -I$ALICE_ROOT/include -I$ALICE_ROOT/RAW -I$ALICE_ROOT/TOF")
+// .L AliTOFRawDataRead.C++
+// AliTOFRawDataRead()
+
+#if !defined(__CINT__) || defined(__MAKECINT__)
+
+// Root include files
+#include "TClonesArray.h"
+
+// AliRoot include files
+#include "AliDAQ.h"
+#include "AliHitMap.h"
+
+#include "AliRawReader.h"
+#include "AliRawReaderFile.h"
+
+#include "AliTOFrawData.h"
+#include "AliTOFRawMap.h"
+#include "AliTOFRawStream.h"
+
+#endif
+
+void AliTOFRawDataRead(Int_t iEvent=0);
+
+void AliTOFRawDataRead(Int_t iEvent)
{
//
// To read TOF raw data
//
+ TClonesArray *fTOFrawData = new TClonesArray("AliTOFrawData",1000);
+ Int_t fPackedDigits=0;
+
Int_t indexDDL = 0;
+ Int_t detectorIndex[5] = {-1, -1, -1, -1, -1};
+ Int_t dummy = -1;
+
+ AliRawReaderFile reader(iEvent);
+ reader.RewindEvents();
+
+ ofstream ftxt;
+ ftxt.open("TOFrawDataReading.txt",ios::app);
+
+ while (reader.NextEvent()) {
+
+ AliTOFRawMap *rawMap = new AliTOFRawMap(fTOFrawData);
+
+ Int_t slot[4] = {-1, -1, -1, -1};
+
+ for (indexDDL = 0; indexDDL < AliDAQ::NumberOfDdls("TOF"); indexDDL++) {
+
+ rawMap->Clear();
+ fTOFrawData->Clear();
+ fPackedDigits = 0;
+
+ printf(" \n \n \n DRM number %2i \n \n \n ", indexDDL);
+
+ reader.Reset();
+ AliTOFRawStream stream(&reader);
+
+ reader.Select("TOF", indexDDL, indexDDL);
+ Bool_t signal = kFALSE;
+
+ while(stream.Next()) {
+
+ signal = (stream.GetSector()!=-1 &&
+ stream.GetPlate()!=-1 &&
+ stream.GetStrip()!=-1 &&
+ stream.GetPadZ()!=-1 &&
+ stream.GetPadX()!=-1);
+
+ if (signal) {
+ printf(" %2i %1i %2i %1i %2i \n", stream.GetSector(), stream.GetPlate(), stream.GetStrip(), stream.GetPadZ(), stream.GetPadX());
+
+ slot[0] = stream.GetTRM();
+ slot[1] = stream.GetTRMchain();
+ slot[2] = stream.GetTDC();
+ slot[3] = stream.GetTDCchannel();
+
+ if (rawMap->TestHit(slot) != kEmpty) {
+
+ AliTOFrawData *rawDigit = static_cast<AliTOFrawData*>(rawMap->GetHit(slot));
+ rawDigit->Update(stream.GetTofBin(), stream.GetToTbin(), stream.GetLeadingEdge(), stream.GetTrailingEdge(), stream.GetPSbit(), stream.GetACQ(), stream.GetErrorFlag());
+
+ }
+ else {
+
+ TClonesArray &arrayTofRawData = *fTOFrawData;
+ new (arrayTofRawData[fPackedDigits++]) AliTOFrawData(stream.GetTRM(), stream.GetTRMchain(), stream.GetTDC(), stream.GetTDCchannel(),
+ stream.GetTofBin(), stream.GetToTbin(), stream.GetLeadingEdge(), stream.GetTrailingEdge(),
+ stream.GetPSbit(), stream.GetACQ(), stream.GetErrorFlag());
+
+ rawMap->SetHit(slot);
+
+ printf(" %6i \n", rawMap->GetHitIndex(slot));
+
+ }
+
+ } // end if (signal)
+
+ } // closed loop on TOF raw data per current DDL file
+
+ printf("\n \n \n end of reading DRM number %2i\n", indexDDL);
+ printf(" packed data %5i\n \n \n ", fTOFrawData->GetEntriesFast());
+
+ for (Int_t iRawData = 0; iRawData<fTOFrawData->GetEntriesFast(); iRawData++) {
+
+ AliTOFrawData *tofRawDatum = (AliTOFrawData*)fTOFrawData->UncheckedAt(iRawData);
- AliRawReaderFile reader(ievt);
+ if (tofRawDatum->GetTOT()==-1 || tofRawDatum->GetTOF()==-1) continue;
- for (indexDDL = 0; indexDDL < 72; indexDDL++) {
+ if (indexDDL<10) ftxt << " " << indexDDL;
+ else ftxt << " " << indexDDL;
+ if (tofRawDatum->GetTRM()<10) ftxt << " " << tofRawDatum->GetTRM();
+ else ftxt << " " << tofRawDatum->GetTRM();
+ ftxt << " " << tofRawDatum->GetTRMchain();
+ if (tofRawDatum->GetTDC()<10) ftxt << " " << tofRawDatum->GetTDC();
+ else ftxt << " " << tofRawDatum->GetTDC();
+ ftxt << " " << tofRawDatum->GetTDCchannel();
- reader.Reset();
- AliTOFRawStream stream(&reader);
- reader.Select(5, indexDDL, indexDDL);
+ stream.EquipmentId2VolumeId(indexDDL, tofRawDatum->GetTRM(), tofRawDatum->GetTRMchain(),
+ tofRawDatum->GetTDC(), tofRawDatum->GetTDCchannel(), detectorIndex);
+ dummy = detectorIndex[3];
+ detectorIndex[3] = detectorIndex[4];
+ detectorIndex[4] = dummy;
- //FILE *fpw = fopen("TOFrawDataRead.txt","w");
+ if (detectorIndex[0]<10) ftxt << " -> " << detectorIndex[0];
+ else ftxt << " -> " << detectorIndex[0];
+ ftxt << " " << detectorIndex[1];
+ if (detectorIndex[2]<10) ftxt << " " << detectorIndex[2];
+ else ftxt << " " << detectorIndex[2];
+ ftxt << " " << detectorIndex[3];
+ if (detectorIndex[4]<10) ftxt << " " << detectorIndex[4];
+ else ftxt << " " << detectorIndex[4];
- while(stream.Next()) {
+ if (tofRawDatum->GetTOT()<10) ftxt << " " << tofRawDatum->GetTOT();
+ else if (tofRawDatum->GetTOT()>=10 && tofRawDatum->GetTOT()<100) ftxt << " " << tofRawDatum->GetTOT();
+ else ftxt << " " << tofRawDatum->GetTOT();
+ if (tofRawDatum->GetTOF()<10) ftxt << " " << tofRawDatum->GetTOF() << endl;
+ else if (tofRawDatum->GetTOF()>=10 && tofRawDatum->GetTOF()<100) ftxt << " " << tofRawDatum->GetTOF() << endl;
+ else if (tofRawDatum->GetTOF()>=100 && tofRawDatum->GetTOF()<1000) ftxt << " " << tofRawDatum->GetTOF() << endl;
+ else ftxt << " " << tofRawDatum->GetTOF() << endl;
- if (stream.GetSector()==0) {
+ } // end loop
- printf("%2i %2i %2i %2i %2i %1i %2i %1i %2i %7i %8i\n",
- stream.GetDDL(),stream.GetTRM(),
- stream.GetTDC(),stream.GetChannel(),
- stream.GetSector(),stream.GetPlate(),
- stream.GetStrip(),stream.GetPadZ(),stream.GetPadX(),
- stream.GetADCbin(),stream.GetTofBin());
- /*
- Int_t iSector = stream.GetSector();
- Int_t iPlate = stream.GetPlate();
- Int_t iStrip = stream.GetStrip();
- Int_t iPadZ = stream.GetPadZ();
- Int_t iPadX = stream.GetPadX();
- Int_t iTof = stream.GetTofBin();
- Int_t iAdc = stream.GetADCbin();
-
- fprintf(fpw,"%2i %1i %2i %1i %2i %8i %7i\n",
- iSector, iPlate, iStrip, iPadZ, iPadX, iTof, iAdc);
- */
+ } // endl loop on DDL files
- }
+ iEvent++;
- }
+ } // end while loop on event
- }
+ ftxt.close();
}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff