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 //___________________________________________________________
34 AliTPCBuffer::AliTPCBuffer():TObject(),
43 //____________________________________________________________
44 AliTPCBuffer::AliTPCBuffer(const char* fileName):TObject(),
51 f.open(fileName,ios::binary|ios::out);
53 f.open(fileName,ios::out);
55 // fout=new TFile(fileName,"recreate");
56 // tree=new TTree("tree","Values");
58 remove("TPCdigits.txt");
61 //////////////////////////////////////////////////////////////////////////////////////////////////////////////
62 AliTPCBuffer::~AliTPCBuffer(){
63 // The destructor closes the IO stream
68 //////////////////////////////////////////////////////////////////////////////////////////////////////////////
69 AliTPCBuffer::AliTPCBuffer(const AliTPCBuffer &source):TObject(source),
75 this->fVerbose=source.fVerbose;
78 //////////////////////////////////////////////////////////////////////////////////////////////////////////////
79 AliTPCBuffer& AliTPCBuffer::operator=(const AliTPCBuffer &source){
81 this->fVerbose=source.fVerbose;
84 //////////////////////////////////////////////////////////////////////////////////////////////////////////////
86 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){
87 //flag=0 the whole row is written to the root file
88 //flag=1 only value in the range [minPad,MaxPasd] are written to the root file
89 //flag=2 complementary case of 1
93 tree->Branch("sec",&sec,"sec/I");
94 tree->Branch("SubSec",&SubSec,"SubSec/I");
95 tree->Branch("row",&row,"row/I");
96 tree->Branch("Pad",&Pad,"Pad/I");
97 tree->Branch("Dig",&Dig,"Dig/I");
98 tree->Branch("Time",&Time,"Time/I");
101 Dig=digrow->CurrentDigit(); //adc
102 Time=digrow->CurrentRow(); //time
103 Pad =digrow->CurrentColumn(); // pad
104 // cout<<"Sec "<<sec<<" row "<<row<<" Pad "<<Pad<<" Dig "<<Dig<<" Time "<<Time<<endl;
113 if((Pad>=minPad)&&(Pad<=maxPad)){
120 if((Pad<minPad)||(Pad>maxPad)){
128 }while (digrow->Next());
133 //////////////////////////////////////////////////////////////////////////////////////////////////////////////
134 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){
135 //It writes TPC digits as par the flag specifications. Being called by AliTPCDDL.C
136 //flag=0 the whole row is written into the file
137 //flag=1 only value in the range [minPad,MaxPasd] are written into the file
138 //flag=2 complementary case of 1
152 if (!digrow->First()) return;
157 ftxt.open("TPCdigits.txt",ios::app);
159 ddlNumber=sec*2+SubSec;
161 ddlNumber=72+(sec-36)*4+SubSec;
164 data.Dig=digrow->CurrentDigit(); //adc
165 data.Time=digrow->CurrentRow(); //time
166 data.Pad =digrow->CurrentColumn(); // pad
168 if (padID!=data.Pad){
169 ftxt<<"S:"<<data.Sec<<" DDL:"<<ddlNumber<<" R:"<<data.Row<<" P:"<<data.Pad<<endl;
176 f.write((char*)(&data),sizeof(data));
178 ftxt<<"A:"<<data.Dig<<" T:"<<data.Time<<endl;
182 if((data.Pad>=minPad)&&(data.Pad<=maxPad)){
183 f.write((char*)(&data),sizeof(data));
185 ftxt<<"A:"<<data.Dig<<" T:"<<data.Time<<endl;
191 if((data.Pad<minPad)||(data.Pad>maxPad)){
192 f.write((char*)(&data),sizeof(data));
194 ftxt<<"A:"<<data.Dig<<" T:"<<data.Time<<endl;
201 }while (digrow->Next());
206 //////////////////////////////////////////////////////////////////////////////////////////////////////////////