#include "AliVZEROBuffer.h"
#include "AliRunDigitizer.h"
#include "AliVZEROdigit.h"
+#include "AliDAQ.h"
+#include "AliRawReader.h"
+#include "AliVZERORawStream.h"
ClassImp(AliVZERO)
-
+ //__________________________________________________________________
+AliVZERO::AliVZERO(): AliDetector(),
+ fIdSens1(0),
+ fThickness(0.),
+ fThickness1(0.),
+ fMaxStepQua(0.),
+ fMaxStepAlu(0.),
+ fMaxDestepQua(0.),
+ fMaxDestepAlu(0.)
+{
+/// Default Constructor
+
+ AliDebug(1,Form("default (empty) ctor this=%p",this));
+ fIshunt = 0;
+}
//_____________________________________________________________________________
AliVZERO::AliVZERO(const char *name, const char *title)
- : AliDetector(name,title)
+ : AliDetector(name,title),
+ fIdSens1(0),
+ fThickness(4.4),
+ fThickness1(2.0),
+ fMaxStepQua(0.05),
+ fMaxStepAlu(0.01),
+ fMaxDestepQua(-1.0),
+ fMaxDestepAlu(-1.0)
{
- //
+
// Standard constructor for VZERO Detector
- //
+
+ AliDebug(1,Form("ctor this=%p",this));
// fIshunt = 1; // All hits are associated with primary particles
gAlice->GetMCApp()->AddHitList(fHits);
- fThickness = 4.4; // total thickness of the V0R box in cm
- fThickness1 = 2.0; // thickness of scintillating cells in cm
-
- fMaxStepQua = 0.05;
- fMaxStepAlu = 0.01;
+// fThickness = 4.4; // total thickness of the V0R box in cm
+// fThickness1 = 2.0; // thickness of scintillating cells in cm
+//
+// fMaxStepQua = 0.05;
+// fMaxStepAlu = 0.01;
+//
+// fMaxDestepQua = -1.0;
+// fMaxDestepAlu = -1.0;
- fMaxDestepQua = -1.0;
- fMaxDestepAlu = -1.0;
- SetMarkerColor(kRed);
}
//_____________________________________________________________________________
// Builds VZEROgetter (AliLoader type)
// if detector wants to use customized getter, it must overload this method
//
-
- Info("MakeLoader","Creating AliVZEROLoader. Top folder is %s.",topfoldername);
+// Info("MakeLoader","Creating AliVZEROLoader. Top folder is %s.",topfoldername);
+
+ AliDebug(1,Form("Creating AliVZEROLoader, Top folder is %s ",topfoldername));
fLoader = new AliVZEROLoader(GetName(),topfoldername);
return fLoader;
}
fVZERO->SetTreeAddress();
digits->GetBranch("VZERODigit")->SetAddress(&VZEROdigits);
- const char *fileName = "VZERO_3584.ddl";
+ const char *fileName = AliDAQ::DdlFileName("VZERO",0);
AliVZEROBuffer* buffer = new AliVZEROBuffer(fileName);
// Verbose level
ofstream ftxt;
buffer->SetVerbose(0);
- Int_t fVerbose = buffer->GetVerbose();
+ Int_t verbose = buffer->GetVerbose();
+
+ // Get Trigger information first
+ // Read trigger inputs from trigger-detector object
+ AliDataLoader * dataLoader = fLoader->GetDigitsDataLoader();
+ if( !dataLoader->IsFileOpen() )
+ dataLoader->OpenFile( "READ" );
+ AliTriggerDetector* trgdet = (AliTriggerDetector*)dataLoader->GetDirectory()->Get( "Trigger" );
+ UInt_t triggerInfo = 0;
+ if(trgdet) {
+ triggerInfo = trgdet->GetMask() & 0xffff;
+ }
+ else {
+ AliError(Form("There is no trigger object for %s",fLoader->GetName()));
+ }
+ buffer->WriteTriggerInfo((UInt_t)triggerInfo);
+ // Now write the channel information: charge+time
+ // We assume here an ordered (by PMNumber) array of
+ // digits!!
Int_t nEntries = Int_t(digits->GetEntries());
-
for (Int_t i = 0; i < nEntries; i++) {
fVZERO->ResetDigits();
Int_t ndig = VZEROdigits->GetEntriesFast();
if(ndig == 0) continue;
- if(fVerbose == 2) {ftxt.open("VZEROdigits.txt",ios::app);}
+ if(verbose == 2) {ftxt.open("VZEROdigits.txt",ios::app);}
for(Int_t k=0; k<ndig; k++){
AliVZEROdigit* fVZERODigit = (AliVZEROdigit*) VZEROdigits->At(k);
- Int_t ADC = fVZERODigit->ADC();
- Int_t cell = fVZERODigit->CellNumber();
- Int_t Time = fVZERODigit->Time();
- if(fVerbose == 1) { cout <<"DDL: "<<fileName<< "\tdigit number: "<< k<<"\tcell: "
- <<cell<<"\tADC: "<< ADC << "\tTime: "<< Time << endl;}
- if(fVerbose == 2) {
- ftxt<<"DDL: "<<fileName<< "\tdigit number: "<< k<<"\tcell: "
- <<cell<<"\tADC: "<< ADC << endl;
+ Int_t ADC = fVZERODigit->ADC();
+ Int_t PMNumber = fVZERODigit->PMNumber();
+ Int_t Time = fVZERODigit->Time();
+ if(verbose == 1) { cout <<"DDL: "<<fileName<< "\tdigit number: "<< k<<"\tPM number: "
+ <<PMNumber<<"\tADC: "<< ADC << "\tTime: "<< Time << endl;}
+ if(verbose == 2) {
+ ftxt<<"DDL: "<<fileName<< "\tdigit number: "<< k<<"\tPM number: "
+ <<PMNumber<<"\tADC: "<< ADC << "\tTime: "<< Time << endl;
}
- buffer->WriteBinary(cell, ADC, Time);
+ buffer->WriteChannel(PMNumber, ADC, Time);
}
- if(fVerbose==2) ftxt.close();
+ if(verbose==2) ftxt.close();
}
+ buffer->WriteScalers();
+ buffer->WriteMBInfo();
+
delete buffer;
fLoader->UnloadDigits();
}
+//_____________________________________________________________________________
+Bool_t AliVZERO::Raw2SDigits(AliRawReader* rawReader){
+ // Converts the VZERO raw data into digits
+ // The method is used for merging simulated and
+ // real data events
+ TStopwatch timer;
+ timer.Start();
+
+ if(!fLoader) {
+ AliError("no VZERO loader found");
+ return kFALSE; }
+
+ TTree* treeD = fLoader->TreeD();
+ if(!treeD) {
+ fLoader->MakeTree("D");
+ treeD = fLoader->TreeD(); }
+
+ AliVZEROdigit digit;
+ AliVZEROdigit* pdigit = &digit;
+ const Int_t kBufferSize = 4000;
+
+ treeD->Branch("VZERO", "AliVZEROdigit", &pdigit, kBufferSize);
+
+ rawReader->Reset();
+ AliVZERORawStream* rawStream = new AliVZERORawStream(rawReader);
+
+ if (!rawStream->Next()) return kFALSE; // No VZERO data found
+
+ for(Int_t i=0; i<64; i++) {
+ new(pdigit) AliVZEROdigit(i, (Int_t)rawStream->GetADC(i), (Int_t)rawStream->GetTime(i));
+ treeD->Fill();
+ }
+
+// Checks if everything is OK by printing results
+
+// for(int i=0;i<64;i++) {
+// printf("Channel %d : %d %d \n",i,rawStream->GetADC(i),rawStream->GetTime(i)); }
+// treeD->Print(); printf(" \n");
+
+ fLoader->WriteDigits("OVERWRITE");
+ fLoader->UnloadDigits();
+
+ delete rawStream;
+
+ timer.Stop();
+ timer.Print();
+}
+