]> git.uio.no Git - u/mrichter/AliRoot.git/commitdiff
git://git.uio.no / u / mrichter / AliRoot.git / commitdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | patch | inline | side by side (parent: 4826571)
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 patch | blob | blame | history
EMCAL/EMCALraw/AliCaloRawAnalyzerFakeALTRO.cxx patch | blob | blame | history
EMCAL/EMCALraw/AliCaloRawAnalyzerNN.cxx patch | blob | blame | history
EMCAL/EMCALraw/AliCaloRawAnalyzerPeakFinder.cxx patch | blob | blame | history

diff --git a/EMCAL/EMCALraw/AliCaloRawAnalyzerCrude.cxx b/EMCAL/EMCALraw/AliCaloRawAnalyzerCrude.cxx
index 27cf481550837af2b26b5233c0794012b55268ef..62785956e90438ef4bffa8d46f35e0677a344b33 100644 (file)
--- a/EMCAL/EMCALraw/AliCaloRawAnalyzerCrude.cxx
+++ b/EMCAL/EMCALraw/AliCaloRawAnalyzerCrude.cxx
@@ -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    
diff --git a/EMCAL/EMCALraw/AliCaloRawAnalyzerFakeALTRO.cxx b/EMCAL/EMCALraw/AliCaloRawAnalyzerFakeALTRO.cxx
index 3e8107b6af233f8407c5eb20ecaa29a4ccd465db..1b04cd5127e4fca445f8b1cf8ac714a8048e42c8 100644 (file)
--- a/EMCAL/EMCALraw/AliCaloRawAnalyzerFakeALTRO.cxx
+++ b/EMCAL/EMCALraw/AliCaloRawAnalyzerFakeALTRO.cxx
@@ -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
diff --git a/EMCAL/EMCALraw/AliCaloRawAnalyzerNN.cxx b/EMCAL/EMCALraw/AliCaloRawAnalyzerNN.cxx
index 5e79beb2aa9f45b7ea46298f32a699008cfad678..d53bdd9f153450b0ffb5cf74d1830b7d65f89f4f 100644 (file)
--- a/EMCAL/EMCALraw/AliCaloRawAnalyzerNN.cxx
+++ b/EMCAL/EMCALraw/AliCaloRawAnalyzerNN.cxx
@@ -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) ); 
 
 }
diff --git a/EMCAL/EMCALraw/AliCaloRawAnalyzerPeakFinder.cxx b/EMCAL/EMCALraw/AliCaloRawAnalyzerPeakFinder.cxx
index e8bc413026a5babdcd1e56af0b941c1dce375363..b2c0b7ce3c9fddb4586b6fe0b4ef3d2a149b9bcf 100644 (file)
--- a/EMCAL/EMCALraw/AliCaloRawAnalyzerPeakFinder.cxx
+++ b/EMCAL/EMCALraw/AliCaloRawAnalyzerPeakFinder.cxx
@@ -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