1 void AliTRDdigits2cluster()
4 /////////////////////////////////////////////////////////////////////////
6 // Creates cluster from the digit information.
8 /////////////////////////////////////////////////////////////////////////
10 // Dynamically link some shared libs
11 if (gClassTable->GetID("AliRun") < 0) {
12 gROOT->LoadMacro("loadlibs.C");
14 cout << "Loaded shared libraries" << endl;
17 // Input and output file names
18 Char_t *infile = "galice.root";
19 Char_t *outfile = "AliTRDclusters.root";
21 // Create the clusterizer
22 AliTRDclusterizerV1 *clusterizer =
23 new AliTRDclusterizerV1("clusterizer","Clusterizer class");
26 TFile *parfile = TFile::Open(infile);
27 AliTRDparameter *par = (AliTRDparameter *) parfile->Get("TRDparameter");
29 clusterizer->SetParameter(par);
32 clusterizer->SetVerbose(1);
35 TTree * te = (TTree*)parfile->Get("TE");
36 Int_t nev = (Int_t)te->GetEntries();
38 for(Int_t iev=0;iev<nev;iev++) {
40 // Open the AliRoot file
41 clusterizer->Open(infile,iev);
42 //clusterizer->Open(infile,outfile,0);
46 clusterizer->ReadDigits();
50 clusterizer->MakeClusters();
52 // Write the cluster tree into file AliTRDclusters.root
53 clusterizer->WriteClusters(-1);
56 // Save the clusterizer class in the AliROOT file
57 // clusterizer->Write();