1 /**************************************************************************
2 * Copyright(c) 1998-2003, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
17 // Storing digits in a binary file
18 // according to the DDL mapping
19 // To be used in Alice Data Challenges
20 // This class is used by AliTPCDDL.C macro
23 #include <Riostream.h>
24 #include <TObjArray.h>
25 #include "AliTPCBuffer.h"
26 #include "AliSimDigits.h"
31 ClassImp(AliTPCBuffer)
32 //////////////////////////////////////////////////////////////////////////////////////////////////////////////
33 AliTPCBuffer::AliTPCBuffer(const char* fileName){
36 f.open(fileName,ios::binary|ios::out);
38 f.open(fileName,ios::out);
40 // fout=new TFile(fileName,"recreate");
41 // tree=new TTree("tree","Values");
44 remove("TPCdigits.txt");
47 //////////////////////////////////////////////////////////////////////////////////////////////////////////////
48 AliTPCBuffer::~AliTPCBuffer(){
49 // The destructor closes the IO stream
54 //////////////////////////////////////////////////////////////////////////////////////////////////////////////
55 AliTPCBuffer::AliTPCBuffer(const AliTPCBuffer &source):TObject(source){
57 this->fVerbose=source.fVerbose;
60 //////////////////////////////////////////////////////////////////////////////////////////////////////////////
61 AliTPCBuffer& AliTPCBuffer::operator=(const AliTPCBuffer &source){
63 this->fVerbose=source.fVerbose;
66 //////////////////////////////////////////////////////////////////////////////////////////////////////////////
68 void AliTPCBuffer::WriteRow(Int_t eth,AliSimDigits *digrow,Int_t minPad,Int_t maxPad,Int_t flag,Int_t sec,Int_t SubSec,Int_t row){
69 //flag=0 the whole row is written to the root file
70 //flag=1 only value in the range [minPad,MaxPasd] are written to the root file
71 //flag=2 complementary case of 1
75 tree->Branch("sec",&sec,"sec/I");
76 tree->Branch("SubSec",&SubSec,"SubSec/I");
77 tree->Branch("row",&row,"row/I");
78 tree->Branch("Pad",&Pad,"Pad/I");
79 tree->Branch("Dig",&Dig,"Dig/I");
80 tree->Branch("Time",&Time,"Time/I");
83 Dig=digrow->CurrentDigit(); //adc
84 Time=digrow->CurrentRow(); //time
85 Pad =digrow->CurrentColumn(); // pad
86 // cout<<"Sec "<<sec<<" row "<<row<<" Pad "<<Pad<<" Dig "<<Dig<<" Time "<<Time<<endl;
95 if((Pad>=minPad)&&(Pad<=maxPad)){
102 if((Pad<minPad)||(Pad>maxPad)){
110 }while (digrow->Next());
115 //////////////////////////////////////////////////////////////////////////////////////////////////////////////
116 void AliTPCBuffer::WriteRowBinary(Int_t eth,AliSimDigits *digrow,Int_t minPad,Int_t maxPad,Int_t flag,Int_t sec,Int_t SubSec,Int_t row){
117 //It writes TPC digits as par the flag specifications. Being called by AliTPCDDL.C
118 //flag=0 the whole row is written into the file
119 //flag=1 only value in the range [minPad,MaxPasd] are written into the file
120 //flag=2 complementary case of 1
139 ftxt.open("TPCdigits.txt",ios::app);
141 ddlNumber=sec*2+SubSec;
143 ddlNumber=72+(sec-36)*4+SubSec;
146 data.Dig=digrow->CurrentDigit(); //adc
147 data.Time=digrow->CurrentRow(); //time
148 data.Pad =digrow->CurrentColumn(); // pad
150 if (padID!=data.Pad){
151 ftxt<<"S:"<<data.Sec<<" DDL:"<<ddlNumber<<" R:"<<data.Row<<" P:"<<data.Pad<<endl;
158 f.write((char*)(&data),sizeof(data));
160 ftxt<<"A:"<<data.Dig<<" T:"<<data.Time<<endl;
164 if((data.Pad>=minPad)&&(data.Pad<=maxPad)){
165 f.write((char*)(&data),sizeof(data));
167 ftxt<<"A:"<<data.Dig<<" T:"<<data.Time<<endl;
173 if((data.Pad<minPad)||(data.Pad>maxPad)){
174 f.write((char*)(&data),sizeof(data));
176 ftxt<<"A:"<<data.Dig<<" T:"<<data.Time<<endl;
183 }while (digrow->Next());
188 //////////////////////////////////////////////////////////////////////////////////////////////////////////////