MUON + CheckCompiler
[u/mrichter/AliRoot.git] / EMCAL / AliCaloRawAnalyzerFitter.cxx
CommitLineData
9e9f5ff6 1// -*- mode: c++ -*-
2
3/**************************************************************************
4 * This file is property of and copyright by the Experimental Nuclear *
5 * Physics Group, Dep. of Physics *
6 * University of Oslo, Norway, 2007 *
7 * *
8 * Author: Per Thomas Hille <perthi@fys.uio.no> for the ALICE HLT Project.*
9 * Contributors are mentioned in the code where appropriate. *
10 * Please report bugs to perthi@fys.uio.no *
11 * *
12 * Permission to use, copy, modify and distribute this software and its *
13 * documentation strictly for non-commercial purposes is hereby granted *
14 * without fee, provided that the above copyright notice appears in all *
15 * copies and that both the copyright notice and this permission notice *
16 * appear in the supporting documentation. The authors make no claims *
17 * about the suitability of this software for any purpose. It is *
18 * provided "as is" without express or implied warranty. *
19 **************************************************************************/
7afee7f0 20
9e9f5ff6 21#include "AliCaloRawAnalyzerFitter.h"
22#include "TF1.h"
9e9f5ff6 23#include <iostream>
7afee7f0 24
9e9f5ff6 25using std::cout;
26using std::endl;
27
9e9f5ff6 28
29AliCaloRawAnalyzerFitter::AliCaloRawAnalyzerFitter(const char *name, const char *nameshort ) :AliCaloRawAnalyzer( name, nameshort),
30 fkEulerSquared(7.389056098930650227),
52fa3af1 31 fTf1(0)
9e9f5ff6 32{
52fa3af1 33 // ctor
34
9e9f5ff6 35 for(int i=0; i < ALTROMAXSAMPLES; i++)
36 {
37 fXaxis[i] = i;
38 }
39
9e9f5ff6 40 fTf1 = new TF1( "myformula", "[0]*((x - [1])/[2])^2*exp(-2*(x -[1])/[2])", 0, 30 );
41
52fa3af1 42 if (fFixTau) // Declared in AliCaloRawAnalyzer
9e9f5ff6 43 {
44 fTf1->FixParameter(2, fTau);
45 }
46 else
47 {
48 fTf1->ReleaseParameter(2); // allow par. to vary
49 fTf1->SetParameter(2, fTau);
50 }
51}
52
53
54AliCaloRawAnalyzerFitter::~AliCaloRawAnalyzerFitter()
55{
52fa3af1 56 // Dtor
57
48b634ff 58 delete fTf1;
9e9f5ff6 59}
60
61
62void
63AliCaloRawAnalyzerFitter::PrintFitResult(const TF1 *f) const
64{
52fa3af1 65 // print fit results
66
9e9f5ff6 67 cout << endl;
68 cout << __FILE__ << __LINE__ << "Using this samplerange we get" << endl;
69 cout << __FILE__ << __LINE__ << "AMPLITUDE = " << f->GetParameter(0)/fkEulerSquared << ",.. !!!!" << endl;
70 cout << __FILE__ << __LINE__ << "TOF = " << f->GetParameter(1) << ",.. !!!!" << endl;
71 cout << __FILE__ << __LINE__ << "NDF = " << f->GetNDF() << ",.. !!!!" << endl;
9e9f5ff6 72 cout << endl << endl;
73}
74
4074cc41 75