Obsolete medium 609 commented out (I.Hrivnacova). Sum of G10 fractions normalized...

[u/mrichter/AliRoot.git] / TRD / AliTRDdigits2cluster.C

diff --git a/TRD/AliTRDdigits2cluster.C b/TRD/AliTRDdigits2cluster.C
index f295779df279c22c231c95ae3981d0c8eda1671d..d2e763a035a96b549f81ecf1d7b01ce52f014fba 100644 (file)
--- a/TRD/AliTRDdigits2cluster.C
+++ b/TRD/AliTRDdigits2cluster.C
@@ -22,25 +22,37 @@ void AliTRDdigits2cluster()
   AliTRDclusterizerV1 *clusterizer = 
     new AliTRDclusterizerV1("clusterizer","Clusterizer class"); 
 
-  // Define the parameter object
-  // If no external parameter object is defined, 
-  // default parameter will be used
-  AliTRDparameter *parameter = new AliTRDparameter("TRDparameter"
-                                                 ,"TRD parameter class");
-  clusterizer->SetParameter(parameter);
-
-  // Open the AliRoot file 
-  clusterizer->Open(infile,outfile,0);
-
-  // Load the digits
-  clusterizer->ReadDigits();
+  // Read the parameter
+  TFile *parfile = TFile::Open(infile);
+  AliTRDparameter *par = (AliTRDparameter *) parfile->Get("TRDparameter"); 
+  par->ReInit();
+  clusterizer->SetParameter(par);
+
+  // Set the parameter
+  clusterizer->SetVerbose(1);
+
+  //Number of events
+  TTree * te = (TTree*)parfile->Get("TE");
+  Int_t nev = (Int_t)te->GetEntries();
+
+  for(Int_t iev=0;iev<nev;iev++) {
+
+    // Open the AliRoot file 
+    clusterizer->Open(infile,iev);
+    //clusterizer->Open(infile,outfile,0);
+
+
+    // Load the digits
+    clusterizer->ReadDigits();
+    clusterizer->Dump();
  
-  // Find the cluster
-  clusterizer->MakeClusters();
+    // Find the cluster
+    clusterizer->MakeClusters();
 
-  // Write the cluster tree into file AliTRDclusters.root
-  clusterizer->WriteClusters(-1);
+    // Write the cluster tree into file AliTRDclusters.root
+    clusterizer->WriteClusters(-1);
 
+  }
   // Save the clusterizer class in the AliROOT file
   // clusterizer->Write();