Comments, grayspace.

[u/mrichter/AliRoot.git] / T0 / AliT0SetCDB.C

diff --git a/T0/AliT0SetCDB.C b/T0/AliT0SetCDB.C
index adc4bb81e6e942c68e8da4c59aeda93e013b4126..803d6d8ec18f6db70e3b817442d30f2b8c492ab1 100644 (file)
--- a/T0/AliT0SetCDB.C
+++ b/T0/AliT0SetCDB.C
@@ -15,7 +15,8 @@
 
 #include "AliRun.h"
 #include "AliT0CalibData.h"
-#include "AliT0AlignData.h"
+//#include "AliT0AlignData.h"
+#include "AliT0Align.h"
 #include "AliCDBMetaData.h"
 #include "AliCDBId.h"
 #include "AliCDBEntry.h"
@@ -62,9 +63,12 @@ void SetAC()
   DBFolder  ="local://Align";
   firstRun  =  0;
   lastRun   =  10;
-  objFormat = "T0 array Z positions";
+  objFormat = "T0 array  positions";
 
-  
+  AliT0Align *al = new AliT0Align(1,835615);
+  al->Run();
+
+   /* 
   AliT0AlignData *alignda=new AliT0AlignData("T0");
   alignda-> SetZposition (67.9,373);
   alignda->Print();
@@ -87,6 +91,7 @@ void SetAC()
 
    storage->Put(alignda, id, &md);
  }
+   */
 }
 //------------------------------------------------------------------------
 void SetTimeDelay()
@@ -108,13 +113,25 @@ void SetTimeDelay()
   AliT0CalibData *calibda=new AliT0CalibData("T0");
   
    Float_t fTimeDelay  = 1000;
-  
-  
+
+   //for timeDelayDA now we are using Int mean time (+-3RMS fit) 
+   //in number of channel for time in each chanel for 
+   //1000 Hijingpara nparticles=500 & zvertex=0 produced by readDigits.C
+   // this is way to simulate time positon for vertex=0 
+   //with unknown time delay in channel
+   /*
+   Int_t timedelayDA[24] = {501, 509, 511,510 ,510, 510, 509, 509, 510,
+                           510, 509, 508, 511, 510, 509, 508,510,
+                           510, 511, 510, 510, 510, 509, 510};
+ 
+   */
   for(Int_t ipmt=0; ipmt<24; ipmt++) {
     calibda->SetTimeDelayCFD(fTimeDelay+ipmt*100,ipmt);
-    //   calibda->SetTimeDelayLED(fTimeDelay,ipmt);
-    calibda->SetTimeDelayLED(0,ipmt);
+       //   calibda->SetTimeDelayLED(fTimeDelay,ipmt);
+//    calibda->SetTimeDelayDA(timedelayDA[ipmt],ipmt);
+    calibda->SetTimeDelayDA(500,ipmt);
   }
+  calibda->SetMeanT0(499);
   calibda->Print();
   //Store calibration data into database
   AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT");
@@ -159,12 +176,8 @@ void SetWalk()
   TRandom rn;
   
   for(Int_t ipmt=0; ipmt<24; ipmt++) {
-    //    calibda->SetWalk(ipmt,"data/CFD-Amp.txt");
-    //    calibda->SetSlewingLED(ipmt,"data/CFD-LED.txt");
-    //  calibda->SetSlewingRec(ipmt,"data/CFD-LED.txt");
     calibda->SetWalk(ipmt);
-    calibda->SetSlewingLED(ipmt);
-    calibda->SetSlewingRec(ipmt);
+    calibda->SetAmpLEDRec(ipmt);
   }
   //Store calibration data into database
   AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT");
@@ -202,6 +215,12 @@ void GetTimeDelay()
   
   AliT0CalibData *clb = (AliT0CalibData*)entry->GetObject();
   clb->Print();
+  for (Int_t i=0; i<24; i++) {
+   cout<<clb->GetTimeDelayCFD(i)<<" "<<clb->GetTimeDelayDA(i)<<endl;
+   cout<<" equalizing CFD "<<(clb->GetTimeDelayCFD(i)-clb->GetTimeDelayCFD(0))<<endl;
+  cout<<" equalizing DA "<<(clb->GetTimeDelayDA(i)-clb->GetTimeDelayDA(0))<<endl;
+  }
+  
 }
 //------------------------------------------------------------------------
 void GetWalk()
@@ -220,7 +239,7 @@ void GetWalk()
    AliT0CalibData *clb = (AliT0CalibData*)entry->GetObject();
    Int_t ipmt=0;
    //  cin>>" enter channel number">>ipmt;
-   TGraph *gr = clb->GetWalk(ipmt); 
+   TGraph *gr = clb->GetAmpLEDRec(ipmt); 
    gr->Draw("AP");
 }
 //------------------------------------------------------------------------
@@ -261,7 +280,8 @@ void SetLookUp()
 
   AliT0CalibData *calibda=new AliT0CalibData("T0");
 
-  calibda->ReadAsciiLookup("lookUpTable.txt");
+//  calibda->ReadAsciiLookup("lookUpTable.txt");
+  calibda->ReadAsciiLookup("lookUpTable_tanay.txt");
 
   //Store calibration data into database
   AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT");
@@ -304,7 +324,7 @@ void GetLookUp()
   cout<<" AliT0CalibData *clb "<<clb <<endl;
   //cout<<"clb->a="<<clb->GetA()<<endl;
   //  clb->Dump();
-  for (Int_t i=0; i<6; i++) 
-    clb->PrintLookup("all",0,i,4);
+  for (Int_t i=0; i<20; i++) 
+    clb->PrintLookupNames("all",i);
 
 }