From 6acf9a0358983008fe029daceeee39302591cf70 Mon Sep 17 00:00:00 2001 From: Samantha Wong Date: Fri, 27 Mar 2026 12:21:12 -0400 Subject: [PATCH] ECPL cutoff parameter corrections --- src/VSpectralFitter.cpp | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/VSpectralFitter.cpp b/src/VSpectralFitter.cpp index bb6572809..30ba079cc 100644 --- a/src/VSpectralFitter.cpp +++ b/src/VSpectralFitter.cpp @@ -165,12 +165,12 @@ bool VSpectralFitter::defineFitFunction() else if( fSpectralFitFunction == 1 ) { cout << "Fitfunction: power law with exponential cutoff" << endl; - sprintf( hname, "[0] * TMath::Power( TMath::Power( 10, x ) / %f, [1] ) * TMath::Exp( -1. * TMath::Power( 10, x ) / [2] )", fSpectralFitFluxNormalisationEnergy ); + sprintf( hname, "[0] * TMath::Power( TMath::Power( 10, x ) / %f, [1] ) * TMath::Exp( -1. * TMath::Power( 10, x ) * [2] )", fSpectralFitFluxNormalisationEnergy ); fFitFunction = new TF1( fFitName.c_str(), hname, log10( fSpectralFitEnergy_min ), log10( fSpectralFitEnergy_max ) ); fFitFunction->SetParameter( 0, 1.e-7 ); fFitFunction->SetParameter( 1, -2. ); - fFitFunction->SetParameter( 2, 10. ); - sprintf( hname, "[0] * TMath::Power( x / %f, [1] ) * TMath::Exp( -1. * x / [2] )", fSpectralFitFluxNormalisationEnergy ); + fFitFunction->SetParameter( 2, 0.1 ); + sprintf( hname, "[0] * TMath::Power( x / %f, [1] ) * TMath::Exp( -1. * x * [2] )", fSpectralFitFluxNormalisationEnergy ); fFitFunction_lin = new TF1( iFitName_lin.c_str(), hname, fSpectralFitEnergy_min, fSpectralFitEnergy_max ); } // broken power law fit @@ -281,7 +281,7 @@ void VSpectralFitter::print() cout << "Gamma = " << fixed << setprecision( 2 ) << fFitFunction->GetParameter( 1 ); cout << " +- " << fFitFunction->GetParError( 1 ) << endl; cout << "Ecut = " << fixed << setprecision( 2 ) << fFitFunction->GetParameter( 2 ); - cout << " +- " << fFitFunction->GetParError( 2 ) << " TeV" << endl; + cout << " +- " << fFitFunction->GetParError( 2 ) << " TeV^-1" << endl; cout << "Chi2 " << setprecision( 2 ) << fFitFunction->GetChisquare(); cout << ", N = " << fFitFunction->GetNDF(); if( fFitFunction->GetNDF() > 0. )