]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWG2/SPECTRA/AliAnalysisCentralExtrapolate.cxx
Deleting AcoReco.C macro
[u/mrichter/AliRoot.git] / PWG2 / SPECTRA / AliAnalysisCentralExtrapolate.cxx
index 6907ac7ec5c6570e3efa21805b7ffdd02b4e87e6..1e06267ffd7ca76e292c45e5d4450afb4a389c25 100644 (file)
@@ -59,14 +59,15 @@ AliAnalysisCentralExtrapolate::AliAnalysisCentralExtrapolate(const char *name)
 {
 // Constructor
 
-//  printf("AliAnalysisCentralExtrapolate::AliAnalysisCentralExtrapolate(const char *name)\n");
+       Printf("Running %s!\n", name);
+
        fInputList = new TList();
        fResultsList = new TList();
-       
+
 }
-  
+
 AliAnalysisCentralExtrapolate::~AliAnalysisCentralExtrapolate() {
-       
+
        if(fInputList) delete fInputList;
        if(fResultsList) delete fResultsList;
 
@@ -178,7 +179,7 @@ void AliAnalysisCentralExtrapolate::ApplyEff(){
 //     efficiency->Draw();
 
 
-       TH1F *hNoEvt = dynamic_cast<TH1F*>(fInputList->FindObject("TaskCentral_NoEvt"));
+       TH1D *hNoEvt = dynamic_cast<TH1D*>(fInputList->FindObject("TaskCentral_NoEvt"));
        if(!hNoEvt){
                printf("Unable to get the number of events! \n");
                return;
@@ -320,7 +321,7 @@ void AliAnalysisCentralExtrapolate::TsallisFit(){
 
 
 // create a new, "extended" histogram
-    TH1F *hPtExtTsallis = new TH1F("PtExtTsallis","Pt Corr Norm Ext",25,0.0,5.0);
+    TH1D *hPtExtTsallis = new TH1D("PtExtTsallis","Pt Corr Norm Ext",25,0.0,5.0);
     
     Double_t bin, binerr, test;
     
@@ -365,7 +366,7 @@ Double_t func(Double_t r){
     Double_t rho = TMath::ATanH(betaR);
 
 
-       Double_t mt = sqrt(mass*mass + pt*pt);// !!! par[0]= pt !!!!!
+       mt = sqrt(mass*mass + pt*pt);// !!! par[0]= pt !!!!!
 
     Double_t argI0 = (pt*TMath::SinH(rho))/T; //T = par[1]
        if(argI0>700.0) return 0.0; // !! floating point exception protection
@@ -467,7 +468,7 @@ void AliAnalysisCentralExtrapolate::BoltzmannFit(){
     printf("p2 = %f\tparam2Err = %f\n", p2, param2Err);
 
 
-    TH1F *hPtExtBoltzmann = new TH1F("PtExtBoltzmann","Pt Corr Norm Ext",25,0.0,5.0);
+    TH1D *hPtExtBoltzmann = new TH1D("PtExtBoltzmann","Pt Corr Norm Ext",25,0.0,5.0);
 
     Double_t bin, binerr, test;