* provided "as is" without express or implied warranty. *
**************************************************************************/
-// *******************************************
-// * particle level cuts for azimuthal isotropic *
+// ***************************************************
+// * particle level cuts for azimuthal isotropic *
// * expansion in highly central collisions analysis *
-// * author: Cristian Andrei *
-// * acristian@niham.nipne.ro *
-// * *****************************************
+// * author: Cristian Andrei *
+// * acristian@niham.nipne.ro *
+// ***************************************************
-#include <TF1.h>
#include <TFile.h>
+#include <TF1.h>
-#include "AliAnalysisCentralCutESD.h"
#include "AliESDtrack.h"
-
-class TObject;
-class TParticle;
-
+#include "AliAnalysisCentralCutESD.h"
//____________________________________________________________________
ClassImp(AliAnalysisCentralCutESD)
TFile *f = TFile::Open("$ALICE_ROOT/PWG2/data/PriorProbabilities.root ");
if(!f){
printf("Can't open PWG2 prior probabilities file!\n Exiting ...\n");
- exit(1);
+ return;
}
fElectronFunction = (TF1 *)f->Get("fitElectrons");
fMuonFunction = (TF1 *)f->Get("fitMuons");
if(!track){
printf("AliAnalysisCentralCutESD:IsSelected ->Can't get track!\n");
- exit(1);
+ return kFALSE;
}
if(fReqCharge){
}
if(fReqPID){
- if(!IsA(track, fPartType)) return kFALSE;
+ if(!IsA(track, fPartType)) return kFALSE;
}
return kTRUE;
else{
printf("Unknown PID method!\n");
- exit(1);
+ return kFALSE;
}
if((AliPID::ParticleCode(partType)) != fPartType){
+#ifndef ALIANALYSISCENTRALCUTESD_H
+#define ALIANALYSISCENTRALCUTESD_H
/*
* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved.
* See cxx source for full Copyright notice
* $Id$
*/
-// *********************************************
-// * ESD particle level cuts for azimuthal isotropic *
-// * expansion in highly central collisions analysis *
-// * author: Cristian Andrei *
-// * acristian@niham.nipne.ro *
-// * ********************************************
+// ***************************************************
+// * ESD particle level cuts for azimuthal isotropic *
+// * expansion in highly central collisions analysis *
+// * author: Cristian Andrei *
+// * acristian@niham.nipne.ro *
+// ***************************************************
-
-#ifndef ALI_ANALYSIS_CENTRAL_CUT_ESD_H
-#define ALI_ANALYSIS_CENTRAL_CUT_ESD_H
-
-#include <TString.h>
#include <TPDGCode.h>
#include "AliAnalysisCuts.h"
-
class TObject;
class TList;
class TF1;
+class TString;
class AliESDtrack;
* provided "as is" without express or implied warranty. *
**************************************************************************/
-// *******************************************
-// * ESD event level cuts for azimuthal isotropic *
+// ***************************************************
+// * ESD event level cuts for azimuthal isotropic *
// * expansion in highly central collisions analysis *
-// * author: Cristian Andrei *
-// * acristian@niham.nipne.ro *
-// * *****************************************
+// * author: Cristian Andrei *
+// * acristian@niham.nipne.ro *
+// ***************************************************
#include "TMath.h"
#include <TObjArray.h>
-#include <TList.h>
#include "AliESDEvent.h"
#include "AliESDtrack.h"
#include "AliESDtrackCuts.h"
- #include "AliMultiplicity.h"
-
+#include "AliMultiplicity.h"
#include "AliAnalysisCentralCutESD.h"
-#include "AliAnalysisCentralCutEvtESD.h"
-
-class TObject;
+#include "AliAnalysisCentralCutEvtESD.h"
//____________________________________________________________________
ClassImp(AliAnalysisCentralCutEvtESD)
AliESDEvent *esdEvent = dynamic_cast<AliESDEvent *>(obj);
if(!esdEvent){
- printf("AliAnalysisCentralCutEvtESD:IsSelected ->Can't get ESD event!\n");
- exit(1);
+ printf("AliAnalysisCentralCutEvtESD:IsSelected ->Can't get ESD event!\n");
+ return kFALSE;
}
-
if(fReqMult){
Int_t mult = CalcMult(esdEvent);
if((mult<fMultMin)||(mult>fMultMax)){
}
}
-
-
if(fReqDir){
Double_t dir = CalcDir(esdEvent);
if((dir<fDirMin)||(dir>fDirMax)){
+#ifndef ALIANALYSISCENTRALCUTEVTESD_H
+#define ALIANALYSISCENTRALCUTEVTESD_H
/*
* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved.
* See cxx source for full Copyright notice
* $Id$
*/
-// ********************************************
+// *****************************************************
// * ESD event level cuts for azimuthal isotropic *
// * expansion in highly central collisions analysis *
-// * author: Cristian Andrei *
-// * acristian@niham.nipne.ro *
-// * ******************************************
-
-#ifndef ALIANALYSISCENTRALCUTEVTESD_H
-#define ALIANALYSISCENTRALCUTEVTESD_H
-
+// * author: Cristian Andrei *
+// * acristian@niham.nipne.ro *
+// *****************************************************
#include "AliAnalysisCuts.h"
-
class TObject;
class TObjArray;
class TList;
class AliESDEvent;
-
class AliAnalysisCentralCutEvtESD: public AliAnalysisCuts {
public:
AliAnalysisCentralCutEvtESD(const char *name="AliAnalysisCentralCutEvtESD", const char *title="ESD_Event_cuts");
* provided "as is" without express or implied warranty. *
**************************************************************************/
-// *******************************************
+// ***************************************************
// * MC event level cuts for azimuthal isotropic *
-// * expansion in highly central collisions analysis *
-// * author: Cristian ANDREI *
-// * acristian@niham.nipne.ro *
-// * ******************************************
-
+// * expansion in highly central collisions analysis *
+// * author: Cristian ANDREI *
+// * acristian@niham.nipne.ro *
+// ***************************************************
#include <TParticle.h>
#include "AliStack.h"
#include "AliAnalysisCentralCutEvtMC.h"
-
-
-class TObject;
-
//____________________________________________________________________
ClassImp(AliAnalysisCentralCutEvtMC)
+#ifndef ALIANALYSISCENTRALCUTEVTMC_H
+#define ALIANALYSISCENTRALCUTEVTMC_H
/*
* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved.
* See cxx source for full Copyright notice
* $Id$
*/
-// ********************************************
+// ***************************************************
// * MC event level cuts for azimuthal isotropic *
-// * expansion in highly central collisions analysis *
-// * author: Cristian Andrei *
-// * acristian@niham.nipne.ro *
-// * ******************************************
-
-#ifndef ALI_ANALYSIS_CENTRAL_CUT_EVT_MC_H
-#define ALI_ANALYSIS_CENTRAL_CUT_EVT_MC_H
-
+// * expansion in highly central collisions analysis *
+// * author: Cristian Andrei *
+// * acristian@niham.nipne.ro *
+// ***************************************************
#include "AliAnalysisCuts.h"
-
class TObject;
class AliMCEvent;
class AliStack;
-
class AliAnalysisCentralCutEvtMC: public AliAnalysisCuts {
public:
AliAnalysisCentralCutEvtMC(const char *name="AliAnalysisCentralCutEvtMC", const char *title="MC_cuts");
* provided "as is" without express or implied warranty. *
**************************************************************************/
-// *********************************************
-// * MC particle level cuts for azimuthal isotropic *
-// * expansion in highly central collisions analysis *
-// * author: Cristian Andrei *
-// * acristian@niham.nipne.ro *
-// * ********************************************
+// ***************************************************
+// * MC particle level cuts for azimuthal isotropic *
+// * expansion in highly central collisions analysis *
+// * author: Cristian Andrei *
+// * acristian@niham.nipne.ro *
+// * *************************************************
#include <TParticle.h>
#include "AliMCParticle.h"
#include "AliAnalysisCentralCutMC.h"
-
-class TObject;
-
//____________________________________________________________________
ClassImp(AliAnalysisCentralCutMC)
+#ifndef ALIANALYSISCENTRALCUTMC_H
+#define ALIANALYSISCENTRALCUTMC_H
/*
* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved.
* See cxx source for full Copyright notice
* $Id$
*/
-// *********************************************
-// * MC particle level cuts for azimuthal isotropic *
-// * expansion in highly central collisions analysis *
-// * author: Cristian Andrei *
-// * acristian@niham.nipne.ro *
-// * *******************************************
-
-#ifndef ALI_ANALYSIS_CENTRAL_CUT_MC_H
-#define ALI_ANALYSIS_CENTRAL_CUT_MC_H
-
+// ***************************************************
+// * MC particle level cuts for azimuthal isotropic *
+// * expansion in highly central collisions analysis *
+// * author: Cristian Andrei *
+// * acristian@niham.nipne.ro *
+// * *************************************************
#include "AliAnalysisCuts.h"
-
class TObject;
class TList;
class AliMCEvent;
class AliMCParticle;
-
class AliAnalysisCentralCutMC: public AliAnalysisCuts{
public:
AliAnalysisCentralCutMC(const char *name="AliAnalysisCentralCutMC", const char *title="MC_cuts");
* provided "as is" without express or implied warranty. *
**************************************************************************/
-// **************************************************
-// * pt spectra extrapolation in the 0. - 0.2 region using *
+// *********************************************************
+// * pt spectra extrapolation in the 0. - 0.2 region using *
// * Boltzmann-Gibbs Blast Wave model or Tsallis Blast *
// * Wave model for azimuthal isotropic expansion in *
-// * highly central collisions analysis *
-// * author: Cristian Andrei *
-// * acristian@niham.nipne.ro *
-// * *************************************************
-
+// * highly central collisions analysis *
+// * author: Cristian Andrei *
+// * acristian@niham.nipne.ro *
+// * *******************************************************
#include "TH1D.h"
#include "TFile.h"
if(!data){
printf("Unable to get CFContainer! \n");
- exit(1);
+ return;
}
if(!eff){
printf("No Eff Grid found! \n");
- exit(1);
+ return;
}
// ccorrdata->Divide(1,2);
// ccorrdata->cd(1);
-// ccorrdata->cd(1)->SetLogy();
+// ccorrdata->cd(1)->SetLogy();
AliCFDataGrid *corrdata = new AliCFDataGrid("corrdata","corrected data",*data);
TH1F *hNoEvt = dynamic_cast<TH1F*>(fInputList->FindObject("AliAnalysisCentral_NoEvt"));
if(!hNoEvt){
printf("Unable to get the number of events! \n");
- exit(1);
+ return;
}
Int_t noEvt = (Int_t)(hNoEvt->GetEntries());
printf("\n** No of processed events = %i **\n",noEvt);
-
+
Double_t scale = 1.0/noEvt;
-
+
TH1D *hPtESDCorrNorm = (TH1D*)hPtESDCorr->Clone("hPtESDCorrNorm");
hPtESDCorrNorm->Scale(scale);
timer.Start();
TH1D *ptSpectra = dynamic_cast<TH1D*>(fResultsList->FindObject("hPtESDCorrNorm"));
+ if(!ptSpectra){
+ printf("BoltzmannFit: Can't get the normalized spectrum\n");
+ return;
+ }
+
+ if(!ptSpectra->GetEntries()){
+ printf("BoltzmannFit: The fit data is empty!\n");
+ return;
+ }
TVirtualFitter::SetDefaultFitter("Minuit2");//options Minuit, Minuit2, Fumili, Fumili2
TH1::SetDefaultSumw2();
TH1::AddDirectory(0);
-
+
timer.Start();
TH1D *ptSpectra = dynamic_cast<TH1D*>(fResultsList->FindObject("hPtESDCorrNorm"));
+ if(!ptSpectra){
+ printf("BoltzmannFit: Can't get the normalized spectrum\n");
+ return;
+ }
+
+ if(!ptSpectra->GetEntries()){
+ printf("BoltzmannFit: The fit data is empty!\n");
+ return;
+ }
gStyle->SetOptStat("neRM");
+#ifndef ALIANALYSISCENTRALEXTRAPOLATE_H
+#define ALIANALYSISCENTRALEXTRAPOLATE_H
/*
* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved.
* See cxx source for full Copyright notice
* $Id$
*/
-// **************************************************
-// * pt spectra extrapolation in the 0. - 0.2 region using *
+// *********************************************************
+// * pt spectra extrapolation in the 0. - 0.2 region using *
// * Boltzmann-Gibbs Blast Wave model or Tsallis Blast *
// * Wave model for azimuthal isotropic expansion in *
-// * highly central collisions analysis *
-// * author: Cristian Andrei *
-// * acristian@niham.nipne.ro *
-// * *************************************************
+// * highly central collisions analysis *
+// * author: Cristian Andrei *
+// * acristian@niham.nipne.ro *
+// * *******************************************************
-#ifndef ALI_ANALYSIS_CENTRAL_EXTRAPOLATE_H
-#define ALI_ANALYSIS_CENTRAL_EXTRAPOLATE_H
#include "TObject.h"
* provided "as is" without express or implied warranty. *
**************************************************************************/
-// *************************************
-// * analysis task for azimuthal isotropic *
-// * expansion in highly central collisions *
-// * author: Cristian Andrei *
-// * acristian@niham.nipne.ro *
-// * ***********************************
+// ******************************************
+// * analysis task for azimuthal isotropic *
+// * expansion in highly central collisions *
+// * author: Cristian Andrei *
+// * acristian@niham.nipne.ro *
+// ******************************************
#include "TChain.h"
#include "TTree.h"
// Some cuts (ie MC IsPrimary) need the MC Event info
if (!fCutsList) {
- printf("No particle cut list found!\n\n");
- exit(1);
+ printf("No particle cut list found!\n\n");
+ return;
}
else {
for(Int_t isel=0;isel< 10; isel++){
- if(!fCutsList[isel])continue;
+ if(!fCutsList[isel]) continue;
TObjArrayIter iter(fCutsList[isel]);
AliAnalysisCuts *cut = 0;
TString cutName=cut->GetName();
if(!cutName){
printf("No cutname!\n");
- exit(1);
+ return;
}
Bool_t checkCut=cutName.Contains("AliAnalysisCentralCutMC");
TList *outList = dynamic_cast<TList*>(GetOutputData(0));
if(!outList){
printf("Unable to get output list! \n");
- exit(1);
+ return;
}
AliAnalysisCentralExtrapolate *extPi = new AliAnalysisCentralExtrapolate("extrapolationpi");
+#ifndef ALI_ANALYSIS_TASK_CENTRAL_H
+#define ALI_ANALYSIS_TASK_CENTRAL_H
/*
* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved.
* See cxx source for full Copyright notice
* $Id$
*/
-// *************************************
-// * analysis task for azimuthal isotropic *
-// * expansion in highly central collisions *
-// * author: Cristian Andrei *
-// * acristian@niham.nipne.ro *
-// * ************************************
-
-
-#ifndef ALI_ANALYSIS_TASK_CENTRAL_H
-#define ALI_ANALYSIS_TASK_CENTRAL_H
+// ******************************************
+// * analysis task for azimuthal isotropic *
+// * expansion in highly central collisions *
+// * author: Cristian Andrei *
+// * acristian@niham.nipne.ro *
+// ******************************************
#include "AliAnalysisTask.h"