#include <TFile.h>
#include <TFolder.h>
#include <TROOT.h>
-#include <TTask.h>
#include <TTree.h>
#include <TVirtualMC.h>
#include <TStopwatch.h>
#include "AliLog.h"
#include "AliMC.h"
#include "AliRun.h"
+#include "AliDAQ.h"
+#include "AliRawReader.h"
-#include "AliTOFDDLRawData.h"
+//#include "AliTOFDDLRawData.h"
#include "AliTOFDigitizer.h"
#include "AliTOFdigit.h"
#include "AliTOFhitT0.h"
#include "AliTOFSDigitizer.h"
#include "AliTOFSDigit.h"
#include "AliTOF.h"
+#include "AliTOFrawData.h"
+#include "AliTOFRawStream.h"
class AliTOFcluster;
-extern TFile *gFile;
extern TROOT *gROOT;
-extern TVirtualMC *gMC;
extern AliRun *gAlice;
-
ClassImp(AliTOF)
-
-AliTOFDDLRawData AliTOF::fgTOFRawWriter = AliTOFDDLRawData();
-
-
//_____________________________________________________________________________
AliTOF::AliTOF():
fFGeom(0x0),
- fDTask(0x0),
- fReTask(0x0),
fSDigits(0x0),
fNSDigits(0),
fReconParticles(0x0),
fIdSens(-1),
fTZero(kFALSE),
fTOFHoles(kTRUE),
- fTOFGeometry(0x0)
+ fTOFGeometry(0x0),
+ fTOFRawWriter(AliTOFDDLRawData())
{
//
// Default constructor
:
AliDetector(name,title),
fFGeom(0x0),
- fDTask(0x0),
- fReTask(0x0),
fSDigits(0x0),
fNSDigits(0),
fReconParticles(0x0),
fIdSens(-1),
fTZero(kFALSE),
fTOFHoles(kTRUE),
- fTOFGeometry(0x0)
+ fTOFGeometry(0x0),
+ fTOFRawWriter(AliTOFDDLRawData())
{
//
// AliTOF standard constructor
}
//____________________________________________________________________________
-void AliTOF::SetTOFSectors(Int_t *sectors)
+void AliTOF::SetTOFSectors(Int_t * const sectors)
{
// Setter for partial/full TOF configuration
void AliTOF::CreateTOFFolders()
{
// create the ALICE TFolder
- // create the ALICE TTasks
// create the ALICE main TFolder
// to be done by AliRun
aliceF->SetOwner() ;
// geometry folder
TFolder * geomF = aliceF->AddFolder("Geometry", "Geometry objects") ;
- TFolder * aliceT = alice->AddFolder("tasks", "Alice tasks Folder") ;
- // make it the owner of the objects that it contains
- aliceT->SetOwner() ;
-
- TTask * aliceDi = new TTask("(S)Digitizer", "Alice SDigitizer & Digitizer") ;
- aliceT->Add(aliceDi);
-
- TTask * aliceRe = new TTask("Reconstructioner", "Alice Reconstructioner") ;
- aliceT->Add(aliceRe);
-
- char * tempo = new char[80] ;
-
- // creates the TOF Digitizer and adds it to alice main (S)Digitizer task
- sprintf(tempo, "%sDigitizers container",GetName() ) ;
- fDTask = new TTask(GetName(), tempo);
- aliceDi->Add(fDTask) ;
-
- // creates the TOF reconstructioner and adds it to alice main Reconstructioner task
- sprintf(tempo, "%sReconstructioner container",GetName() ) ;
- fReTask = new TTask(GetName(), tempo);
- aliceRe->Add(fReTask) ;
-
- delete [] tempo ;
// creates the TOF geometry folder
geomF->AddFolder("TOF", "Geometry for TOF") ;
{
// dtor:
// it remove also the alice folder
- // and task that TOF creates instead of AliRun
/* PH Temporarily commented because of problems
TFolder * alice = (TFolder*)gROOT->GetListOfBrowsables()->FindObject("FPAlice") ;
delete alice;
}
-//_____________________________________________________________________________
-void AliTOF::DrawModule() const
-{
- //
- // Draw a shaded view of the common part of the TOF geometry
- //
-
- AliInfo(" Drawing of AliTOF");
- // Set everything unseen
- gMC->Gsatt("*", "seen", -1);
- //
- // Set ALIC mother transparent
- gMC->Gsatt("ALIC","SEEN",0);
- //
- // Set the volumes visible
- gMC->Gsatt("FTOA","SEEN",1);
- gMC->Gsatt("FTOB","SEEN",1);
- gMC->Gsatt("FTOC","SEEN",1);
- gMC->Gsatt("FLTA","SEEN",1);
- gMC->Gsatt("FLTB","SEEN",1);
- gMC->Gsatt("FLTC","SEEN",1);
- gMC->Gsatt("FSTR","SEEN",1);
- //
- gMC->Gdopt("hide", "on");
- gMC->Gdopt("shad", "on");
- gMC->Gsatt("*", "fill", 7);
- gMC->SetClipBox(".");
- gMC->SetClipBox("*", 0, 1000, -1000, 1000, -1000, 1000);
- gMC->DefaultRange();
- gMC->Gdraw("alic", 40, 30, 0, 12, 9.5, .02, .02);
- gMC->Gdhead(1111, "Time Of Flight");
- gMC->Gdman(18, 4, "MAN");
- gMC->Gdopt("hide","off");
-}
//___________________________________________
void AliTOF::ResetHits ()
/*
// Save the geometry
TDirectory* saveDir = gDirectory;
- AliRunLoader::GetRunLoader()->CdGAFile();
+ AliRunLoader::Instance()->CdGAFile();
fTOFGeometry->Write("TOFGeometry");
saveDir->cd();
*/
AliDetector::MakeBranch(option);
Int_t buffersize = 4000;
- Char_t branchname[10];
- sprintf(branchname,"%s",GetName());
+ const Int_t kSize=10;
+ Char_t branchname[kSize];
+ snprintf(branchname,kSize,"%s",GetName());
const char *oD = strstr(option,"D");
const char *oS = strstr(option,"S");
// AliInfo("Start...");
AliRunLoader * rl = fLoader->GetRunLoader();
+ AliDebug(2,"Initialized runLoader");
AliTOFSDigitizer sd((rl->GetFileName()).Data());
- ToAliDebug(1, sd.Print(""));
+ AliDebug(2,"Initialized TOF sdigitizer");
+ //ToAliDebug(1, sd.Print(""));
+ //AliInfo("ToAliDebug");
- sd.Exec("") ;
+ //sd.Digitize("all") ;
+ sd.Digitize("partial") ;
+
+ AliDebug(2,"I am sorting from AliTOF class");
}
AliTOFSDigitizer sd((rl->GetFileName()).Data(),evNumber1,evNumber2) ;
ToAliDebug(1, sd.Print(""));
- sd.Exec("") ;
+ sd.Digitize("") ;
}
//___________________________________________________________________________
-AliDigitizer* AliTOF::CreateDigitizer(AliRunDigitizer* manager) const
+AliDigitizer* AliTOF::CreateDigitizer(AliDigitizationInput* digInput) const
{
- return new AliTOFDigitizer(manager);
+ AliDebug(2,"I am creating the TOF digitizer");
+ return new AliTOFDigitizer(digInput);
}
//___________________________________________________________________________
-Bool_t AliTOF::CheckOverlap(Int_t* vol, Int_t* digit,Int_t Track)
+Bool_t AliTOF::CheckOverlap(const Int_t * const vol,
+ Int_t* digit,Int_t Track)
{
//
// Checks if 2 or more hits belong to the same pad.
return;
}
- fgTOFRawWriter.Clear();
- fgTOFRawWriter.SetVerbose(0);
- if (fgTOFRawWriter.GetPackedAcquisitionMode()) {
- if(fgTOFRawWriter.GetMatchingWindow()>8192)
+ //AliTOFDDLRawData rawWriter;
+ fTOFRawWriter.Clear();
+ fTOFRawWriter.SetVerbose(0);
+ if (fTOFRawWriter.GetPackedAcquisitionMode()) {
+ if(fTOFRawWriter.GetMatchingWindow()>8192)
AliWarning(Form("You are running in packing mode and the matching window is %.2f ns, i.e. greater than 199.8848 ns",
- fgTOFRawWriter.GetMatchingWindow()*AliTOFGeometry::TdcBinWidth()*1.e-03));
+ fTOFRawWriter.GetMatchingWindow()*AliTOFGeometry::TdcBinWidth()*1.e-03));
}
- fgTOFRawWriter.Clear();
AliDebug(1,"Formatting raw data for TOF");
digits->GetEvent(0);
- fgTOFRawWriter.RawDataTOF(digits->GetBranch("TOF"));
+ fTOFRawWriter.RawDataTOF(digits->GetBranch("TOF"));
fLoader->UnloadDigits();
//
fSDigits = new TClonesArray("AliTOFSDigit", 1000);
}
+//____________________________________________________________________________
+Bool_t AliTOF::Raw2SDigits(AliRawReader* rawReader)
+{
+ //
+ // Converts raw data to sdigits for TOF
+ //
+
+ TStopwatch stopwatch;
+ stopwatch.Start();
+
+ if(!GetLoader()->TreeS()) {MakeTree("S"); MakeBranch("S");}
+ //TClonesArray &aSDigits = *fSDigits;
+
+ AliTOFRawStream tofRawStream = AliTOFRawStream();
+ tofRawStream.Raw2SDigits(rawReader, fSDigits);
+
+ GetLoader()->TreeS()->Fill(); GetLoader()->WriteSDigits("OVERWRITE");//write out sdigits
+ Int_t nSDigits = fSDigits->GetEntries();
+
+ ResetSDigits();
+
+ AliDebug(1, Form("Got %d TOF sdigits", nSDigits));
+ AliDebug(1, Form("Execution time to read TOF raw data and fill TOF sdigit tree : R:%.2fs C:%.2fs",
+ stopwatch.RealTime(),stopwatch.CpuTime()));
+
+ return kTRUE;
+
+}
+
+//____________________________________________________________________________
+void AliTOF::Raw2Digits(AliRawReader* rawReader)
+{
+ //
+ // Converts raw data to digits for TOF
+ //
+
+ TStopwatch stopwatch;
+ stopwatch.Start();
+
+ if(!GetLoader()->TreeD()) {MakeTree("D"); MakeBranch("D");}
+ //TClonesArray &aDigits = *fDigits;
+
+ AliTOFRawStream tofRawStream = AliTOFRawStream();
+ tofRawStream.Raw2Digits(rawReader, fDigits);
+
+ GetLoader()->TreeD()->Fill(); GetLoader()->WriteDigits("OVERWRITE");//write out digits
+ Int_t nDigits = fDigits->GetEntries();
+
+ ResetDigits();
+
+ AliDebug(1, Form("Got %d TOF digits", nDigits));
+ AliDebug(1, Form("Execution time to read TOF raw data and fill TOF digit tree : R:%.2fs C:%.2fs",
+ stopwatch.RealTime(),stopwatch.CpuTime()));
+
+}