Time scale fix for EMCAL fitters
authoranolivei <antonio.carlos.oliveira.da.silva@cern.ch>
Mon, 2 Feb 2015 16:33:39 +0000 (17:33 +0100)
committergconesab <gustavo.conesa.balbastre@cern.ch>
Tue, 3 Feb 2015 10:46:47 +0000 (11:46 +0100)
EMCAL/EMCALraw/AliCaloRawAnalyzerCrude.cxx
EMCAL/EMCALraw/AliCaloRawAnalyzerFakeALTRO.cxx
EMCAL/EMCALraw/AliCaloRawAnalyzerNN.cxx
EMCAL/EMCALraw/AliCaloRawAnalyzerPeakFinder.cxx

index 27cf481..6278595 100644 (file)
@@ -54,7 +54,7 @@ AliCaloRawAnalyzerCrude::Evaluate(const vector<AliCaloBunchInfo> &bunchvector, c
 
       if(  maxf < fAmpCut  ||  ( maxamp - ped) > fOverflowCut  ) // (maxamp - ped) > fOverflowCut = Close to saturation (use low gain then)
        {
-         return  AliCaloFitResults( maxamp, ped, Ret::kCrude, maxf, timebinOffset);
+         return  AliCaloFitResults( maxamp, ped, Ret::kCrude, maxf, (timebinOffset*TIMEBINWITH)-fL1Phase);
        }
       else if ( maxf >= fAmpCut ) // no if statement needed really; keep for readability
        {
@@ -65,7 +65,7 @@ AliCaloRawAnalyzerCrude::Evaluate(const vector<AliCaloBunchInfo> &bunchvector, c
 
          Float_t chi2 = CalculateChi2(maxf, maxrev, first, last);
          Int_t ndf = last - first - 1; // nsamples - 2
-         return AliCaloFitResults( maxamp, ped, Ret::kCrude, maxf, timebinOffset,
+         return AliCaloFitResults( maxamp, ped, Ret::kCrude, maxf, (timebinOffset*TIMEBINWITH)-fL1Phase,
                                    timebinOffset, chi2, ndf, Ret::kDummy, AliCaloFitSubarray(index, maxrev, first, last) ); 
        } // ampcut
     } // bunch index    
index 3e8107b..1b04cd5 100644 (file)
@@ -67,7 +67,7 @@ AliCaloRawAnalyzerFakeALTRO::Evaluate( const vector<AliCaloBunchInfo>  &bunchvec
     short timebinOffset = maxampindex - (bunchvector.at(index).GetLength()-1);
     if(  maxf < fAmpCut  ||  ( maxamp - ped) > fOverflowCut  ) // (maxamp - ped) > fOverflowCut = Close to saturation (use low gain then)
     {
-      return  AliCaloFitResults( maxamp, ped, Ret::kCrude, maxf, timebinOffset);
+      return  AliCaloFitResults( maxamp, ped, Ret::kCrude, maxf, (timebinOffset*TIMEBINWITH)-fL1Phase);
     }
     else if ( maxf >= fAmpCut )
     {
@@ -100,7 +100,7 @@ AliCaloRawAnalyzerFakeALTRO::Evaluate( const vector<AliCaloBunchInfo>  &bunchvec
              }
              catch (const std::exception & e) {
           AliError( Form("TGraph Fit exception %s, fit status %d", e.what(),tmpStatus) );
-          return AliCaloFitResults( maxamp, ped, Ret::kNoFit, maxf, timebinOffset,
+          return AliCaloFitResults( maxamp, ped, Ret::kNoFit, maxf, (timebinOffset*TIMEBINWITH)-fL1Phase,
                                    timebinOffset, Ret::kDummy, Ret::kDummy,  Ret::kDummy, AliCaloFitSubarray(index, maxrev, first, last) );
              }
         
@@ -116,7 +116,7 @@ AliCaloRawAnalyzerFakeALTRO::Evaluate( const vector<AliCaloBunchInfo>  &bunchvec
         delete graph;
         return AliCaloFitResults( maxamp, ped , Ret::kFitPar,
                                  fTf1->GetParameter(0)/fkEulerSquared,
-                                 tmax,
+                                 (tmax*TIMEBINWITH)-fL1Phase,
                                  timebinOffset,
                                  fTf1->GetChisquare(),
                                  fTf1->GetNDF(),
@@ -129,7 +129,7 @@ AliCaloRawAnalyzerFakeALTRO::Evaluate( const vector<AliCaloBunchInfo>  &bunchvec
            {
              Float_t chi2 = CalculateChi2(maxf, maxrev, first, last);
              Int_t ndf = last - first - 1; // nsamples - 2
-             return AliCaloFitResults( maxamp, ped, Ret::kCrude, maxf, timebinOffset,
+             return AliCaloFitResults( maxamp, ped, Ret::kCrude, maxf, (timebinOffset*TIMEBINWITH)-fL1Phase,
                                  timebinOffset, chi2, ndf, Ret::kDummy, AliCaloFitSubarray(index, maxrev, first, last) ); 
            }
     } // ampcut
index 5e79beb..d53bdd9 100644 (file)
@@ -88,7 +88,7 @@ AliCaloRawAnalyzerNN::Evaluate( const vector<AliCaloBunchInfo> &bunchvector,
   if(  maxf < fAmpCut  ||  ( maxamp - ped) > fOverflowCut  ) // (maxamp - ped) > fOverflowCut = Close to saturation (use low gain then)
     {
       //   cout << __FILE__ << __LINE__<< ":  timebinOffset = " <<  timebinOffset  << "  maxf "<< maxf  << endl; 
-      return  AliCaloFitResults( maxamp, ped, Ret::kCrude, maxf, timebinOffset);
+      return  AliCaloFitResults( maxamp, ped, Ret::kCrude, maxf, (timebinOffset*TIMEBINWITH)-fL1Phase);
     }
 
   int first = 0;
@@ -105,7 +105,7 @@ AliCaloRawAnalyzerNN::Evaluate( const vector<AliCaloBunchInfo> &bunchvector,
          chi2 = CalculateChi2(maxf, maxrev, first, last);
          ndf = last - first - 1; // nsamples - 2
          //      cout << __FILE__ << __LINE__<< ":  timebinOffset = " <<  timebinOffset << "  maxf\t"<< maxf <<endl;
-         return AliCaloFitResults( maxamp, ped, Ret::kCrude, maxf, timebinOffset,
+         return AliCaloFitResults( maxamp, ped, Ret::kCrude, maxf, (timebinOffset*TIMEBINWITH)-fL1Phase,
                                    timebinOffset, chi2, ndf, Ret::kDummy, AliCaloFitSubarray(index, maxrev, first, last) ); 
        }
       else
@@ -124,7 +124,7 @@ AliCaloRawAnalyzerNN::Evaluate( const vector<AliCaloBunchInfo> &bunchvector,
          chi2 = CalculateChi2(amp, tof-timebinOffset+maxrev, first, last);
          ndf = last - first - 1; // nsamples - 2
          //cout << __FILE__ << __LINE__<< ":  tof = " <<  tof << "   amp" << amp <<endl;
-         return AliCaloFitResults( maxamp, ped , Ret::kFitPar, amp , tof, timebinOffset, chi2, ndf,
+         return AliCaloFitResults( maxamp, ped , Ret::kFitPar, amp , (tof*TIMEBINWITH)-fL1Phase, timebinOffset, chi2, ndf,
                                    Ret::kDummy, AliCaloFitSubarray(index, maxrev, first, last) );
 
        }
@@ -133,7 +133,7 @@ AliCaloRawAnalyzerNN::Evaluate( const vector<AliCaloBunchInfo> &bunchvector,
   ndf = last - first - 1; // nsamples - 2
   
   // cout << __FILE__ << __LINE__<< ":  timebinOffset = " << timebinOffset <<  "   maxf ="<< maxf  << endl;
-  return AliCaloFitResults( maxamp, ped, Ret::kCrude, maxf, timebinOffset,
+  return AliCaloFitResults( maxamp, ped, Ret::kCrude, maxf, (timebinOffset*TIMEBINWITH)-fL1Phase,
                            timebinOffset, chi2, ndf, Ret::kDummy, AliCaloFitSubarray(index, maxrev, first, last) ); 
 
 }
index e8bc413..b2c0b7c 100644 (file)
@@ -122,7 +122,7 @@ AliCaloRawAnalyzerPeakFinder::Evaluate( const vector<AliCaloBunchInfo> &bunchvec
  
       if(  maxf < fAmpCut  ||  ( maxamp - ped) > fOverflowCut  ) // (maxamp - ped) > fOverflowCut = Close to saturation (use low gain then)
        {
-         return  AliCaloFitResults( maxamp, ped, Ret::kCrude, maxf, timebinOffset);
+         return  AliCaloFitResults( maxamp, ped, Ret::kCrude, maxf, (timebinOffset*TIMEBINWITH)-fL1Phase);
        }            
       else if ( maxf >= fAmpCut )
        {
@@ -177,7 +177,7 @@ AliCaloRawAnalyzerPeakFinder::Evaluate( const vector<AliCaloBunchInfo> &bunchvec
              tof = timebinOffset - 0.01*tof/fAmp - fL1Phase/TIMEBINWITH; // clock
              Float_t chi2 = CalculateChi2(fAmp, tof-timebinOffset+maxrev, first, last);
              Int_t ndf = last - first - 1; // nsamples - 2
-             return AliCaloFitResults( maxamp, ped , Ret::kFitPar, fAmp, tof, 
+             return AliCaloFitResults( maxamp, ped , Ret::kFitPar, fAmp, (tof*TIMEBINWITH)-fL1Phase,
                                        timebinOffset, chi2, ndf,
                                        Ret::kDummy, AliCaloFitSubarray(index, maxrev, first, last) );  
            }
@@ -185,7 +185,7 @@ AliCaloRawAnalyzerPeakFinder::Evaluate( const vector<AliCaloBunchInfo> &bunchvec
            {
              Float_t chi2 = CalculateChi2(maxf, maxrev, first, last);
              Int_t ndf = last - first - 1; // nsamples - 2
-             return AliCaloFitResults( maxamp, ped , Ret::kCrude, maxf, timebinOffset,
+             return AliCaloFitResults( maxamp, ped , Ret::kCrude, maxf, (timebinOffset*TIMEBINWITH)-fL1Phase,
                                        timebinOffset, chi2, ndf, Ret::kDummy, AliCaloFitSubarray(index, maxrev, first, last) ); 
            }
        } // ampcut