RAWDATA closer to reality

[u/mrichter/AliRoot.git] / START / AliSTARTReconstructor.cxx

diff --git a/START/AliSTARTReconstructor.cxx b/START/AliSTARTReconstructor.cxx
index 62c385b1ff5c1ed015e05f0a106a88be3242ba1e..56bd73a6e6e86f060868dfca98f08cc647c74dbb 100644 (file)
--- a/START/AliSTARTReconstructor.cxx
+++ b/START/AliSTARTReconstructor.cxx
@@ -25,56 +25,124 @@
 #include "AliSTARTdigit.h"
 #include "AliSTARTReconstructor.h"
 #include <AliESD.h>
-
+#include <TClonesArray.h>
+#include "AliSTARTRecPoint.h"
+#include "AliRawReader.h"
+#include "AliSTARTRawReader.h"
+#include "AliLog.h"
 ClassImp(AliSTARTReconstructor)
 
-void AliSTARTReconstructor::Reconstruct(AliRunLoader* /*runLoader*/) const 
+  void  AliSTARTReconstructor::ConvertDigits(AliRawReader* rawReader, TTree* digitsTree) const
 {
-// nothing to be done
+  //START raw data-> digits conversion
+ // reconstruct time information from raw data
+  AliSTARTRawReader myrawreader(rawReader,digitsTree);
+  myrawreader.NextThing();
 }
+  void AliSTARTReconstructor::Reconstruct(TTree*digitsTree, TTree*clustersTree) const
+{
+// START digits reconstruction
+// STARTRecPoint writing 
+
+  AliDebug(1,Form("Start DIGITS reconstruction "));
+   Int_t channelWigth=25; //ps
+  TArrayI* fSumMult = new TArrayI(6);
+  Float_t ph2mV = 150./500.;
+  Float_t mV2channel=200000/(25*25);  //5V -> 200ns
+
+  TBranch *brDigits=digitsTree->GetBranch("START");
+  AliSTARTdigit *fDigits = new AliSTARTdigit();
+  if (brDigits) {
+    brDigits->SetAddress(&fDigits);
+  }else{
+    cerr<<"EXEC Branch START digits not found"<<endl;
+    return;
+  }
+  brDigits->GetEntry(0);
+  Int_t besttimeright=channelWigth * (fDigits->BestTimeRight());
+  Int_t besttimeleft=channelWigth * (fDigits->BestTimeLeft());
+
+  //folding with experimental time distribution
+  //  Float_t c = 29.9792; // cm/ns
+  Float_t c = 0.0299792; // cm/ps
+  Float_t lenr=TMath::Sqrt(350*350 + 6.5*6.5);
+  Float_t lenl=TMath::Sqrt(69.7*69.7 + 6.5*6.5);
+  Float_t timeDiff=channelWigth * (fDigits->TimeDiff());
+  Int_t meanTime=channelWigth * (fDigits->MeanTime());
+  Float_t ds=(c*(timeDiff)-(lenr-lenl))/2;
+  AliDebug(2,Form(" timediff in ns %f  real point%f",timeDiff,ds));
+  
 
-void AliSTARTReconstructor::FillESD(AliRunLoader* rl, AliESD *pESD) const
+  fDigits->Print();
+  fDigits->GetSumMult(*fSumMult);
+  Int_t multipl[6]; 
+ for (Int_t i=0; i<6; i++)
+    {
+      Float_t  mult=Float_t (fSumMult->At(i));
+      Float_t   realMultmV=TMath::Exp(mult/mV2channel);
+      multipl[i]=Int_t ((realMultmV/ph2mV)/500+0.5);
+    }
+  AliDebug(2,Form(" multiplicity Abs side %i  multiplicity non-Abs side %i",multipl[1],multipl[2]));
+
+  AliSTARTRecPoint* frecpoints= new AliSTARTRecPoint ();
+  clustersTree->Branch( "START", "AliSTARTRecPoint" ,&frecpoints, 405,1);
+  frecpoints->SetTimeBestRight(besttimeright);
+  frecpoints->SetTimeBestLeft(besttimeleft);
+  frecpoints->SetVertex(ds);
+  frecpoints->SetMeanTime(meanTime);
+  frecpoints->SetMult(multipl[0]);
+  frecpoints->SetMultA(multipl[2]);
+  frecpoints->SetMultC(multipl[1]);
+
+  clustersTree->Fill();
+}
+
+
+void AliSTARTReconstructor::FillESD(AliRunLoader* runLoader, AliESD *pESD) const
 {
+
   /***************************************************
   Resonstruct digits to vertex position
   ****************************************************/
-
-  Float_t c = 0.3;  //speed of light mm/ps
-  Int_t channelWigth=25; //ps
-  if (!rl) {
-    Error("Reconstruct", "No run loader");
+  
+  //  Float_t c = 0.3;  //speed of light mm/ps
+  // Float_t Zposition=0;
+  
+  if (!runLoader) {
+    AliError("Reconstruct >> No run loader");
     return;
   }
+  
+  AliDebug(1,Form("Start FillESD START"));
 
-  if (rl->GetDebug()>1) Info("Reconstruct","START!!!");
-
-  AliSTARTLoader* pStartLoader = (AliSTARTLoader*) rl->GetLoader("STARTLoader");
+  AliSTARTLoader* pStartLoader = (AliSTARTLoader*) runLoader->GetLoader("STARTLoader");
  
-  pStartLoader->LoadDigits();
-  AliSTARTdigit* pDigits=pStartLoader->Digits();
-  if (!pDigits) {
-    Error("Reconstruct", "no digits found");
+  pStartLoader->LoadRecPoints("READ");
+  
+    TTree *treeR = pStartLoader->TreeR();
+  
+   AliSTARTRecPoint* frecpoints= new AliSTARTRecPoint ();
+    if (!frecpoints) {
+    AliError("Reconstruct Fill ESD >> no recpoints found");
     return;
   }
-
-  if (rl->GetDebug()>1) pDigits->Dump();
-  if(pDigits) {
-    Int_t   besttimeright = pDigits->GetBestTimeRight();
-    Int_t besttimeleft  = pDigits->GetBestTimeLeft();
-    Float_t besttimerightPs = Float_t (besttimeright*channelWigth);
-    Float_t besttimeleftPs  = Float_t (besttimeleft*channelWigth);
-   Float_t Zposit=(c*(besttimerightPs-besttimeleftPs)-(3500.-697))/2;
   
+  AliDebug(1,Form("Start FillESD START"));
+   TBranch *brRec = treeR->GetBranch("START");
+    if (brRec) {
+      brRec->SetAddress(&frecpoints);
+    }else{
+      cerr<<"EXEC Branch START rec not found"<<endl;
+      exit(111);
+    } 
+ 
+    brRec->GetEntry(0);
+    Float_t Zposition=frecpoints->GetVertex();
     
-    pESD->SetT0zVertex(Zposit);
-    
-    if (rl->GetDebug()>1) {
-      cout<<" vertex in ESD "<< pESD->GetT0zVertex()<<endl;
-    }
+    pESD->SetT0zVertex(Zposition);
+    pStartLoader->UnloadRecPoints();
     
-  } // vertex in 3 sigma
-  pStartLoader->UnloadDigits();
- }
+} // vertex in 3 sigma