]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - T0/AliT0Reconstructor.cxx
protection against electronic failture
[u/mrichter/AliRoot.git] / T0 / AliT0Reconstructor.cxx
index 688a4066e646ea75dc984804d5e187f0994e8bc6..2663cb9b68be9afb5f18c63cb9bdd7e783181f22 100644 (file)
 #include "AliESDfriend.h"
 #include "AliESDTZEROfriend.h"
 #include "AliLog.h"
+#include "AliCDBEntry.h" 
+#include "AliCDBManager.h"
+#include "AliCTPTimeParams.h"
+#include "AliLHCClockPhase.h"
 
 #include <TArrayI.h>
 #include <TGraph.h>
@@ -55,10 +59,26 @@ ClassImp(AliT0Reconstructor)
                                              fLatencyL1(0),
                                              fLatencyL1A(0),
                                              fLatencyL1C(0),
+                                            fGRPdelays(0),
                                              fESDTZEROfriend(NULL)
 
 {
   //constructor
+  AliCDBEntry *entry = AliCDBManager::Instance()->Get("GRP/CTP/CTPtiming");
+  if (!entry) AliFatal("CTP timing parameters are not found in OCDB !");
+  AliCTPTimeParams *ctpParams = (AliCTPTimeParams*)entry->GetObject();
+  Float_t l1Delay = (Float_t)ctpParams->GetDelayL1L0()*25.0;
+
+  AliCDBEntry *entry1 = AliCDBManager::Instance()->Get("GRP/CTP/TimeAlign");
+  if (!entry1) AliFatal("CTP time-alignment is not found in OCDB !");
+  AliCTPTimeParams *ctpTimeAlign = (AliCTPTimeParams*)entry1->GetObject();
+  l1Delay += ((Float_t)ctpTimeAlign->GetDelayL1L0()*25.0);
+  AliCDBEntry *entry4 = AliCDBManager::Instance()->Get("GRP/Calib/LHCClockPhase");
+  if (!entry4) AliFatal("LHC clock-phase shift is not found in OCDB !");
+  AliLHCClockPhase *phase = (AliLHCClockPhase*)entry4->GetObject();
+
+  fGRPdelays = l1Delay - phase->GetMeanPhase();
 
   fParam = AliT0Parameters::Instance();
   fParam->Init();
@@ -76,7 +96,7 @@ ClassImp(AliT0Reconstructor)
   fLatencyL1A = fParam->GetLatencyL1A();
   fLatencyL1C = fParam->GetLatencyL1C();
   fLatencyHPTDC = fParam->GetLatencyHPTDC();
-  AliDebug(10,Form(" LatencyL1 %f latencyL1A %f latencyL1C %f latencyHPTDC %f \n",fLatencyL1, fLatencyL1A, fLatencyL1C, fLatencyHPTDC));
+  AliDebug(2,Form(" LatencyL1 %f latencyL1A %f latencyL1C %f latencyHPTDC %f \n",fLatencyL1, fLatencyL1A, fLatencyL1C, fLatencyHPTDC));
   
   // fdZonC = TMath::Abs(fParam->GetZPositionShift("T0/C/PMT1"));
   //fdZonA = TMath::Abs(fParam->GetZPositionShift("T0/A/PMT15"));
@@ -94,7 +114,7 @@ void AliT0Reconstructor::Reconstruct(TTree*digitsTree, TTree*clustersTree) const
   
 {
   // T0 digits reconstruction
-  Int_t refAmp = GetRecoParam()->GetRefAmp();
+  Int_t refAmp = Int_t (GetRecoParam()->GetRefAmp());
   
   TArrayI * timeCFD = new TArrayI(24); 
   TArrayI * timeLED = new TArrayI(24); 
@@ -149,7 +169,7 @@ void AliT0Reconstructor::Reconstruct(TTree*digitsTree, TTree*clustersTree) const
       else
        adc[ipmt] = 0;
       
-      time[ipmt] = fCalib-> WalkCorrection(refAmp, ipmt, adc[ipmt],  timeCFD->At(ipmt)) ;
+     time[ipmt] = fCalib-> WalkCorrection(refAmp, ipmt, Int_t(adc[ipmt]),  timeCFD->At(ipmt)) ;
             
       Double_t sl = Double_t(timeLED->At(ipmt) - timeCFD->At(ipmt));
       //    time[ipmt] = fCalib-> WalkCorrection( refAmp,ipmt, Int_t(sl),  timeCFD->At(ipmt) ) ;
@@ -235,8 +255,11 @@ void AliT0Reconstructor::Reconstruct(AliRawReader* rawReader, TTree*recTree) con
   //
   // reference amplitude and time ref. point from reco param
 
-  Int_t refAmp = GetRecoParam()->GetRefAmp();
-  Int_t refPoint = GetRecoParam()->GetRefPoint();
+  Float_t refAmp = GetRecoParam()->GetRefAmp();
+  Int_t refPoint = 0;
+  //Bad channel
+  Int_t badpmt = GetRecoParam()->GetRefPoint();
+  Int_t low[110], high[110];
 
   Int_t allData[110][5];
   
@@ -249,7 +272,9 @@ void AliT0Reconstructor::Reconstruct(AliRawReader* rawReader, TTree*recTree) con
   Float_t meanVertex = 0;
   for (Int_t i0=0; i0<105; i0++)
     {
-      for (Int_t j0=0; j0<5; j0++) allData[i0][j0]=0;  
+      for (Int_t j0=0; j0<5; j0++) allData[i0][j0]=0; 
+      low[i0] = Int_t (GetRecoParam()->GetLow(i0));    
+      high[i0] = Int_t (GetRecoParam()->GetHigh(i0));
     }
    
   Double32_t besttimeA=9999999;
@@ -270,9 +295,15 @@ void AliT0Reconstructor::Reconstruct(AliRawReader* rawReader, TTree*recTree) con
     AliDebug(1,Form(" no raw data found!!"));
   else
     {  
-      Int_t fBCID=Int_t (rawReader->GetBCID());
+   for (Int_t i=0; i<24; i++)
+    {
+      timeCFD[i]=0; timeLED[i]=0; chargeQT0[i]=0; chargeQT1[i]=0;
+    }
+     Int_t fBCID=Int_t (rawReader->GetBCID());
       Int_t trmbunch= myrawreader.GetTRMBunchID();
-      if(type == 7 &&( trmbunch -fBCID )==37 ) {  //only physics 
+      AliDebug(10,Form(" CDH BC ID %i, TRM BC ID %i \n", fBCID, trmbunch ));
+      if(type == 7  ) {  //only physics 
        for (Int_t i=0; i<105; i++) {
        for (Int_t iHit=0; iHit<5; iHit++) 
          {
@@ -287,39 +318,79 @@ void AliT0Reconstructor::Reconstruct(AliRawReader* rawReader, TTree*recTree) con
        
        //       Int_t meanT0 = fParam->GetMeanT0();
        
-       
        for (Int_t in=0; in<12; in++)  
          {
-           timeCFD[in] = allData[in+1][0] ;
-           timeCFD[in+12] = allData[in+56+1][0] ;
-           timeLED[in] = allData[in+12+1][0] ;
-           timeLED[in+12] = allData[in+68+1][0] ;
+           for (Int_t iHit=0; iHit<5; iHit++) 
+             {
+               if(allData[in+1][iHit] > low[in+1] && 
+                  allData[in+1][iHit] < high[in+1])
+                 {
+                   timeCFD[in] = allData[in+1][iHit] ; 
+                   break;
+                 }
+             }
+           for (Int_t iHit=0; iHit<5; iHit++) 
+             {
+               if(allData[in+12+1][iHit] > low[in+12+1] && 
+                  allData[in+1+12][iHit] < high[in+12+1])
+                 {
+                       timeLED[in] = allData[in+12+1][iHit] ;
+                       break;
+                 }
+             }
+           for (Int_t iHit=0; iHit<5; iHit++) 
+             {
+               if(allData[in+1+56][iHit] > low[in+1+56] && 
+                  allData[in+1+56][iHit] < high[in+1+56])
+                 {
+                   timeCFD[in+12] = allData[in+56+1][iHit] ;
+                   break;
+                 }
+             }
+           
+           for (Int_t iHit=0; iHit<5; iHit++) 
+             {
+               if(allData[in+1+68][iHit] > low[in+1+68] && 
+                  allData[in+1+68][iHit] < high[in+1+68])
+                 {
+                   timeLED[in+12] = allData[in+68+1][iHit] ;
+                   break;
+                 }
+             }
            AliDebug(10, Form(" readed i %i cfdC %i cfdA %i ledC %i ledA%i ",
                              in, timeCFD[in],timeCFD[in+12],timeLED[in], 
                              timeLED[in+12]));   
+           
          }
        
+       
        for (Int_t in=0; in<12;  in++)
          {
            chargeQT0[in]=allData[2*in+25][0];
            chargeQT1[in]=allData[2*in+26][0];
-         }
-       
+           AliDebug(10, Form(" readed Raw %i %i %i",
+                             in, chargeQT0[in],chargeQT1[in]));
+         }     
        for (Int_t in=12; in<24;  in++)
          {
            chargeQT0[in]=allData[2*in+57][0];
            chargeQT1[in]=allData[2*in+58][0];
+           AliDebug(10, Form(" readed Raw %i %i %i",
+                             in, chargeQT0[in],chargeQT1[in]));
+           
          }
        
-
-       for (Int_t in=0; in<24; in++)  
-         AliDebug(10, Form(" readed Raw %i %i %i %i %i",
-                           in, timeLED[in],timeCFD[in],chargeQT0[in],chargeQT1[in]));
-        onlineMean = allData[49][0];       
-       
+       for (Int_t iHit=0; iHit<5; iHit++) 
+         {
+             if(allData[49][iHit] > low[49] && 
+                allData[49][iHit] < high[49]){
+               onlineMean = allData[49][iHit];
+               break;
+             }
+         }
        Double32_t time[24], adc[24],  noncalibtime[24];
        for (Int_t ipmt=0; ipmt<24; ipmt++) {
-        if(timeCFD[ipmt]>0 && timeLED[ipmt]>0){
+         if(timeCFD[ipmt]>0 ){
           //for simulated data
             //for physics  data
           if(( chargeQT0[ipmt] - chargeQT1[ipmt])>0)  {
@@ -327,9 +398,7 @@ void AliT0Reconstructor::Reconstruct(AliRawReader* rawReader, TTree*recTree) con
           }
           else
             adc[ipmt] = 0;
-          
-
-          time[ipmt] = fCalib-> WalkCorrection(refAmp, ipmt, adc[ipmt], timeCFD[ipmt] ) ;
+          time[ipmt] = fCalib-> WalkCorrection(Int_t (refAmp), ipmt, Int_t(adc[ipmt]), timeCFD[ipmt] ) ;
           
           Double_t sl = timeLED[ipmt] - timeCFD[ipmt];
           // time[ipmt] = fCalib-> WalkCorrection( refAmp,ipmt, Int_t(sl), timeCFD[ipmt] ) ;
@@ -339,13 +408,11 @@ void AliT0Reconstructor::Reconstruct(AliRawReader* rawReader, TTree*recTree) con
           Double_t qtMip = ((TGraph*)fQTC.At(ipmt))->Eval(adc[ipmt]);
           AliDebug(10,Form("  Amlitude in MIPS LED %f ; QTC %f;  in channels %f\n ",ampMip,qtMip, adc[ipmt]));
           //bad peak removing
-          if(sl<550) {
             frecpoints->SetTime(ipmt, Float_t(time[ipmt]) );
             // frecpoints->SetTime(ipmt,Double32_t(timeCFD[ipmt]));
             frecpoints->SetAmp(ipmt, Double32_t( qtMip)); //for cosmic &pp beam 
             frecpoints->SetAmpLED(ipmt, Double32_t(ampMip));        
             noncalibtime[ipmt]= Double32_t (timeCFD[ipmt]);
-          }
         }
         else {
           time[ipmt] = 0;
@@ -355,7 +422,7 @@ void AliT0Reconstructor::Reconstruct(AliRawReader* rawReader, TTree*recTree) con
        }
        fESDTZEROfriend->SetT0timeCorr(noncalibtime) ;     
        for (Int_t ipmt=0; ipmt<12; ipmt++){
-        if(time[ipmt] > 1 )
+        if(time[ipmt] > 1 && ipmt != badpmt && adc[ipmt]>0.1 )
           {
             if(time[ipmt]<besttimeC){
               besttimeC=time[ipmt]; //timeC
@@ -365,7 +432,7 @@ void AliT0Reconstructor::Reconstruct(AliRawReader* rawReader, TTree*recTree) con
        }
        for ( Int_t ipmt=12; ipmt<24; ipmt++)
         {
-          if(time[ipmt] > 1 )
+          if(time[ipmt] > 1 && ipmt != badpmt && adc[ipmt]>0.1)
             {
               if(time[ipmt]<besttimeA) {
                 besttimeA=time[ipmt]; //timeA
@@ -374,14 +441,14 @@ void AliT0Reconstructor::Reconstruct(AliRawReader* rawReader, TTree*recTree) con
             }
         }
        if(besttimeA < 999999) 
-        frecpoints->SetTimeBestA(besttimeA * channelWidth - 1000.*fLatencyHPTDC + 1000.*fLatencyL1A);
+        frecpoints->SetTimeBestA(besttimeA * channelWidth - 1000.*fLatencyHPTDC + 1000.*fLatencyL1A - 1000.*fGRPdelays);
        if( besttimeC < 999999 ) 
-        frecpoints->SetTimeBestC(besttimeC * channelWidth - 1000.*fLatencyHPTDC +1000.*fLatencyL1C);
+        frecpoints->SetTimeBestC(besttimeC * channelWidth - 1000.*fLatencyHPTDC +1000.*fLatencyL1C - 1000.*fGRPdelays);
        AliDebug(10,Form(" pmtA %i besttimeA %f ps, pmtC %i besttimeC %f ps",
                       pmtBestA,besttimeA, pmtBestC,  besttimeC));
         if(besttimeA <999999 && besttimeC < 999999 ){
         timeDiff = ( besttimeA - besttimeC)* 0.001* channelWidth + fLatencyL1A - fLatencyL1C;
-        timeclock = channelWidth * Float_t( besttimeA+besttimeC)/2. - 1000.*fLatencyHPTDC + 1000.*fLatencyL1;  
+        timeclock = channelWidth * Float_t( besttimeA+besttimeC)/2. - 1000.*fLatencyHPTDC + 1000.*fLatencyL1 - 1000.*fGRPdelays;  
         meanTime = (besttimeA+besttimeC-2.*Float_t(ref))/2.;
         vertex =  meanVertex - c*(timeDiff)/2. ; //+ (fdZonA - fdZonC)/2; 
        }
@@ -397,10 +464,16 @@ void AliT0Reconstructor::Reconstruct(AliRawReader* rawReader, TTree*recTree) con
       Int_t trchan[5]= {50,51,52,55,56};
       for (Int_t i=0; i<5; i++) tr[i]=false; 
       for (Int_t itr=0; itr<5; itr++) {
-       if(allData[trchan[itr]][0]>0) tr[itr]=true;
-       frecpoints->SetT0Trig(tr);
+       for (Int_t iHit=0; iHit<5; iHit++) 
+         {
+           Int_t trr=trchan[itr];
+           if(allData[trr][iHit] > low[trr] &&
+              allData[trr][iHit] < high[trr]&&
+              allData[trr][iHit] > 0) tr[itr]=true;
+         }
       }
-    } // if (else )raw data
+      frecpoints->SetT0Trig(tr);
+     } // if (else )raw data
   recTree->Fill();
   if(frecpoints) delete frecpoints;
 }
@@ -431,10 +504,9 @@ void AliT0Reconstructor::Reconstruct(AliRawReader* rawReader, TTree*recTree) con
     currentVertex = vertex->GetZ();
     
     ncont = vertex->GetNContributors();
-    // cout<<" spdver "<<spdver<<" ncont "<<ncont<<endl;
+    //  cout<<"@@ spdver "<<spdver<<" ncont "<<ncont<<endl;
     if(ncont>2 ) {
       shift = currentVertex/c;
-      //         cout<<" vertex shif "<<shift<<" vertex "<<spdver<<" IsFromVertexer3D  "<<fverSPD->IsFromVertexer3D()<<endl;
     }
   }
   TTree *treeR = clustersTree;
@@ -471,7 +543,8 @@ void AliT0Reconstructor::Reconstruct(AliRawReader* rawReader, TTree*recTree) con
     if ( time[i] >1) {
       ampQTC[i] = frecpoints -> GetAmp(i);
       amp[i] = frecpoints -> AmpLED(i);
-    }
+      AliDebug(1,Form("T0: time %f  ampQTC %f ampLED %f \n", time[i], ampQTC[i], amp[i]));
+   }
   }
   Int_t trig= frecpoints ->GetT0Trig();
   pESD->SetT0Trig(trig);
@@ -491,18 +564,18 @@ void AliT0Reconstructor::Reconstruct(AliRawReader* rawReader, TTree*recTree) con
     if (fr) {
       AliDebug(1, Form("Writing TZERO friend data to ESD tree"));
 
-      if (ncont>2) {
+      //     if (ncont>2) {
        tcorr = fESDTZEROfriend->GetT0timeCorr();
        for ( Int_t i=0; i<24; i++) {
-          timecorr[i]=tcorr[i];
-         if(i<12 && time[i]>1) timecorr[i] -=  shift/channelWidth;
-         if(i>11 && time[i]>1) timecorr[i] +=  shift/channelWidth;
+         if(i<12 && time[i]>1) timecorr[i] = tcorr[i] -  shift/channelWidth;
+         if(i>11 && time[i]>1) timecorr[i] = tcorr[i] +  shift/channelWidth;
+         if(time[i]>1)   AliDebug(1,Form("T0 friend : time %f  ampQTC %f ampLED %f \n", timecorr[i], ampQTC[i], amp[i]));
        }
        fESDTZEROfriend->SetT0timeCorr( timecorr) ;
        fESDTZEROfriend->SetT0ampLEDminCFD(amp);
        fESDTZEROfriend->SetT0ampQTC(ampQTC);
-       fr->SetTZEROfriend(fESDTZEROfriend);
-      }//
+       fr->SetTZEROfriend(fESDTZEROfriend);
+       //      }//
     }
   }