Hits2Reco method added

diff --git a/TPC/AliTPC.cxx b/TPC/AliTPC.cxx
index 88b7c6cd88845c48d94a1e5fcc4241b8b8209edc..9e31dc761a94adb88cf806e36363acc4af716a36 100644 (file)
--- a/TPC/AliTPC.cxx
+++ b/TPC/AliTPC.cxx
@@ -15,6 +15,9 @@
 
 /*
 $Log$
+Revision 1.38  2001/07/20 14:32:43  kowal2
+Processing of many events possible now
+
 Revision 1.37  2001/06/12 07:17:18  kowal2
 Hits2SDigits method implemented (summable digits)
 
@@ -1321,7 +1324,7 @@ void AliTPC::Hits2DigitsSector(Int_t isec)
 
     TObjArray **row;
     
-    printf("*** Processing sector number %d ***\n",isec);
+    //printf("*** Processing sector number %d ***\n",isec);
 
       Int_t nrows =fTPCParam->GetNRow(isec);
 
@@ -1351,7 +1354,7 @@ void AliTPC::Hits2DigitsSector(Int_t isec)
                Int_t ndig = dig->GetDigitSize(); 
         
  
-     printf("*** Sector, row, compressed digits %d %d %d ***\n",isec,i,ndig);
+       //printf("*** Sector, row, compressed digits %d %d %d ***\n",isec,i,ndig);
        
         fDigitsArray->ClearRow(isec,i);  
 
@@ -2548,7 +2551,45 @@ void AliTPC::FindTrackHitsIntersection(TClonesArray * arr)
 
 }
 //_______________________________________________________________________________
-void AliTPC::Digits2Reco(Int_t eventnumber)
+void AliTPC::Digits2Reco(Int_t firstevent,Int_t lastevent)
 {
-  // empty for a time being
+  // produces rec points from digits and writes them on the root file
+  // AliTPCclusters.root
+
+  TDirectory *cwd = gDirectory;
+
+
+  AliTPCParam *dig=(AliTPCParam *)gDirectory->Get("75x40_100x60");
+  SetParam(dig);
+  cout<<"AliTPC::Digits2Reco: TPC parameteres have been set"<<endl; 
+  TFile *out;
+  if(!gSystem->Getenv("CONFIG_FILE")){
+    out=TFile::Open("AliTPCclusters.root","recreate");
+  }
+  else {
+    const char *tmp1;
+    const char *tmp2;
+    char tmp3[80];
+    tmp1=gSystem->Getenv("CONFIG_FILE_PREFIX");
+    tmp2=gSystem->Getenv("CONFIG_OUTDIR");
+    sprintf(tmp3,"%s%s/AliTPCclusters.root",tmp1,tmp2);
+    out=TFile::Open(tmp3,"recreate");
+  }
+
+  TStopwatch timer;
+  cout<<"AliTPC::Digits2Reco - determination of rec points begins"<<endl;
+  timer.Start();
+  cwd->cd();
+  for(Int_t iev=firstevent;iev<lastevent+1;iev++){
+
+    printf("Processing event %d\n",iev);
+    Digits2Clusters(out,iev);
+  }
+  cout<<"AliTPC::Digits2Reco - determination of rec points ended"<<endl;
+  timer.Stop();
+  timer.Print();
+  out->Close();
+  cwd->cd(); 
+
+
 }

diff --git a/TPC/AliTPC.h b/TPC/AliTPC.h
index 13a6e35111afed0e570b3898abfd0e4abb801288..d3ffdbb0a81762e808240c4977347f0fe3f61b61 100644 (file)
--- a/TPC/AliTPC.h
+++ b/TPC/AliTPC.h
@@ -56,11 +56,11 @@ public:
   virtual void  BuildGeometry();
   virtual void  CreateGeometry() {}
   virtual void  CreateMaterials();
-  virtual void  Hits2Clusters(TFile *of, Int_t eventn);
+  virtual void  Hits2Clusters(TFile *of, Int_t eventn=0);
   virtual void  Hits2ExactClustersSector(Int_t isec); // MI change calculate "exact" cluster position
   virtual void  SDigits2Digits(Int_t eventnumber=0);
   virtual void  Hits2SDigits(Int_t eventnumber=0);
-  virtual void  Digits2Reco(Int_t eventnumber=0);
+  virtual void  Digits2Reco(Int_t firstevent=0,Int_t lastevent=0);
   virtual void  Hits2Digits(Int_t eventnumber=0);   //MI change
   virtual void  Hits2DigitsSector(Int_t isec);  //MI change
   virtual void  Init();