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

Timestamp:

02/21/13 14:56:09 (12 years ago)

Author:

Jens Buss

Message:

added dynamic startvalue calculation

File:

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

fact/tools/rootmacros/PulseTemplates/pulse.C

@@ -75 +75 @@
 // PUBLIC OPERATIONS
 // ===========================================================================
+
+void
+Pulse::Fit()
+{
+    Fit(mName, mOptions, mType);
+}
+
+void
+Pulse::Fit(
+        TString fitOptions,
+        int     type
+        )
+{
+    Fit(mName, fitOptions, type);
+}
 
 void
@@ -140 +155 @@
     mHisto              = NULL;
     mModel              = NULL;
-    mOptions            = "S";
+    mOptions            = "ERS";
     mBsl.first          = 0;
     mHeight.first       = 0;
@@ -168 +183 @@
     mChi2               = 0;
     mTimeResolution     = 0.5;
+    mVerbosityLvl       = 6;
 }
 
@@ -181 +197 @@
         )
 {
+    int numPara = 7;
     mFitMin = fitMin;
     mFitMax = fitMax;
-    mModel  = new TF1(fitName, shapeFunc, fitMin, fitMax, 7 );
-
+    mModel  = new TF1(fitName, shapeFunc, fitMin, fitMax, numPara );
+    mModel->SetNpx( mHisto->GetNbinsX());
     // ======================================================================
     // Calculate startvaues for fit
@@ -214 +231 @@
     double para[] = {bsl, amplitude, start, stop, tau, tau, phe};
 
+
+    // cout startvalues
+    if (mVerbosityLvl > 3)
+    {
+        for (unsigned int i = 0; i < numPara; i++)
+        {
+            cout << "para[ " << i << "]: " << para[i] << endl;
+        }
+    }
 
 
@@ -246 +272 @@
         )
 {
+    int numPara = 6;
     mFitMin = fitMin;
     mFitMax = fitMax;
-    mModel  = new TF1(fitName, shapeFunc2, fitMin, fitMax, 6 );
-
+    mModel  = new TF1(fitName, shapeFunc2, fitMin, fitMax, numPara );
+    mModel->SetNpx( mHisto->GetNbinsX());
     // ======================================================================
     // Calculate startvaues for fit
-    double tau      = 30.0;
+    double tau  = 30.0;
     double bsl  = 0; //MW der ersten zehn slices
+//    double bsl2  = 0; //MW der ersten zehn slices
     double gain = 10; //startvalue for gain
-    int first_bin   = mHisto->GetXaxis()->GetFirst();
-
-    for (int i = 0; i < 30; i++){
-        bsl += mHisto->GetBinContent(first_bin+0);
-    }
-    bsl = bsl/30;
+    int counter = 0;
 
     double stop     = mHisto->GetMaximumBin(); //pos of max
-    double amplitude= mHisto->GetBinContent(stop);
-    double start    = stop-10; //pos 10 slices before maximum
+    double amplitude=0;
+
+    for (int i = (stop - 3) ; i < (stop + 3) && i <= fitMax && i >= fitMin ; i++)
+    {
+        counter++;
+        amplitude= mHisto->GetBinContent(i);
+    }
+    amplitude=amplitude/counter;
+
+    counter = 0;
+    for (int i = (stop - 15) ; i > fitMin+2 && i <= fitMax ; i--)
+    {
+        if (mHisto->GetBinContent(i) > 0)
+        {
+            break;
+        }
+        counter++;
+        bsl += mHisto->GetBinContent(i);
+    }
+    if (counter <= 0)
+    {
+        counter = 1;
+    }
+    bsl = bsl/counter;
+
+//    counter = 0;
+//    for (int i = stop + 100; i < stop + 200 && i <= fitMax && i >= fitMin ; i++)
+//    {
+//        if (mHisto->GetBinContent(i) > 0.5*amplitude )
+//        {
+//            break;
+//        }
+//        counter++;
+//        bsl2 += mHisto->GetBinContent(i);
+//    }
+//    if (counter <= 0)
+//    {
+//        counter = 1;
+//    }
+//    bsl2 = bsl2/counter;
+
+    double start    = mHisto->GetBinCenter(stop-10); //pos 10 slices before maximum
     int    phe = 0;
     if (mOrder < 0)
@@ -277 +340 @@
 
 
-    double para[] = {bsl, amplitude, start, tau, tau, phe};
+    double para[] = {bsl, 2*amplitude+gain, start, tau*0.1, tau, phe};
     mModel->SetParameters(para);
     mModel->SetParNames("BSL", "A0", "t0", "Tau1", "Tau2", "PhE");
     mModel->SetLineColor(kBlue);
 
+    // cout startvalues
+    if (mVerbosityLvl > 3)
+    {
+        cout << endl;
+        for (unsigned int i = 0; i < numPara; i++)
+        {
+            cout << "para[ " << i << "]: " << para[i] << endl;
+        }
+    }
     mFitResultPtr   = mHisto->Fit(mModel, fitOptions);
     mBsl.first      = mModel->GetParameter(0);
@@ -401 +473 @@
 double  Pulse::GetTimeResolution(){return mTimeResolution;}
 void    Pulse::SetTimeResolution(double timeResolution){mTimeResolution = timeResolution;}
+void    Pulse::SetVerbosity(int verbLvl){mVerbosityLvl = verbLvl;}
 
 
@@ -498 +571 @@
     double tau2     = par[4];
     int    phe      = par[5];
+//    double bsl2      = par[6];
 
     // helper variables