Changeset 14762 for fact/tools/rootmacros/PulseTemplates/pulse.C

Timestamp:

12/21/12 15:33:51 (12 years ago)

Author:

Jens Buss

Message:

risetime function

File:

• fact/tools/rootmacros/PulseTemplates/pulse.C (modified) (6 diffs)

fact/tools/rootmacros/PulseTemplates/pulse.C

@@ -148 +148 @@
     mAmplitude.second   = 0;
     mPhE.second         = 0;
+    mRisetime.first     = 0;
+    mRisetime.second    = 0;
     mOrder              = 0;
     mType               = 0;
@@ -154 +156 @@
     mFitNdf             = 0;
     mChi2               = 0;
+    mTimeResolution     = 0.5;
 }
 
@@ -275 +278 @@
     mAmplitude.first    = mModel->GetMaximum() - mBsl.first;
 //    mAmplitude.second   = mModel->GetMaximum() - mBsl.first;
+    mRisetime.first     = mTimeResolution*CalculateHistEdgeRisetime();
     mFitProb            = mFitResultPtr->Prob();
     mFitNCalls          = mFitResultPtr->NCalls();
@@ -281 +285 @@
 }
 
+double
+Pulse::CalculateHistEdgeRisetime() // in bins
+{
+    //rise-tim for pulse for pulse height from 10% to 90%
+    double riseTime = 0;
+
+    // Get Pulse Maximum value and bin
+    int maxBin = mHisto->GetMaximumBin();
+    double max = mHisto->GetBinContent(maxBin);
+
+    // values of bin with 10% and 90%
+    int tenPctBin       = 0;
+    int ninetyPctBin    = 0;
+
+    // flag points
+    bool finishedLow    = false;
+    bool finishedHigh   = false;
+
+    //2nd runvariable
+    int j = maxBin;
+
+    //loop over from 20 bins bevore max to maximum
+    for (int i = maxBin - 20; i<maxBin; i++){
+        j--;
+        if (j<0) break;
+
+        if (mHisto->GetBinContent(i) >= 0.1*max && !finishedLow)
+        {
+            tenPctBin = i;
+            finishedLow = true;
+        }
+
+
+        if (mHisto->GetBinContent(j) <= 0.9*max && !finishedHigh)
+        {
+            ninetyPctBin = j;
+            finishedHigh = true;
+        }
+
+        if (finishedLow && finishedHigh)
+        {
+            break;
+        }
+    }
+
+    riseTime = mHisto->GetBinLowEdge(ninetyPctBin) - mHisto->GetBinLowEdge(tenPctBin);
+    return riseTime;
+}
+
+
+
+
 // ===========================================================================
 // ACCESS
@@ -287 +343 @@
 TString Pulse::GetName(){ return mName;}
 double  Pulse::GetBsl(){ return mBsl.first;}
+double  Pulse::GetBslErr(){ return mBsl.second;}
 double  Pulse::GetHeight(){ return mHeight.first;}
+double  Pulse::GetHeightErr(){ return mHeight.second;}
 double  Pulse::GetT0(){ return mT0.first;}
+double  Pulse::GetT0Err(){ return mT0.second;}
 double  Pulse::GetT1(){ return mT1.first;}
+double  Pulse::GetT1Err(){ return mT1.second;}
 double  Pulse::GetTau1(){ return mTau1.first;}
+double  Pulse::GetTau1Err(){ return mTau1.second;}
 double  Pulse::GetTau2(){ return mTau2.first;}
+double  Pulse::GetTau2Err(){ return mTau2.second;}
 double  Pulse::GetIntegral(){ return mIntegral.first;}
+double  Pulse::GetIntegralErr(){ return mIntegral.second;}
 double  Pulse::GetAmplitude(){ return mAmplitude.first;}
+double  Pulse::GetAmplitudeErr(){ return mAmplitude.second;}
 int     Pulse::GetPhE(){ return mPhE.first;}
-double  Pulse::GetBslErr(){ return mBsl.second;}
-double  Pulse::GetHeightErr(){ return mHeight.second;}
-double  Pulse::GetT0Err(){ return mT0.second;}
-double  Pulse::GetT1Err(){ return mT1.second;}
-double  Pulse::GetTau1Err(){ return mTau1.second;}
-double  Pulse::GetTau2Err(){ return mTau2.second;}
-double  Pulse::GetIntegralErr(){ return mIntegral.second;}
-double  Pulse::GetAmplitudeErr(){ return mAmplitude.second;}
 int     Pulse::GetPhEErr(){ return mPhE.second;}
+double  Pulse::GetRiseTime(){ return mRisetime.first; }
+double  Pulse::GetRiseTimeErr(){ return mRisetime.second;}
 int     Pulse::GetOrder(){ return mOrder;}
 int     Pulse::GetType(){return mType;}
@@ -312 +370 @@
 double  Pulse::GetFitNdf(){return mFitNdf;}
 double  Pulse::GetChi2(){return mChi2;}
+double  Pulse::GetTimeResolution(){return mTimeResolution;}
+void    Pulse::SetTimeResolution(double timeResolution){mTimeResolution = timeResolution;}