[u/mrichter/AliRoot.git] / PWGLF / totEt / AliAnalysisHadEt.cxx

//Created by Christine Nattrass, Rebecca Scott, Irakli Martashvili
//University of Tennessee at Knoxville
//
// This class is designed for the analysis of the hadronic component of 
// transverse energy.  It is used by AliAnalysisTaskHadEt.
// This gets information about the hadronic component of the transverse energy 
// from tracks reconstructed in an event
// it has daughters, AliAnalysisHadEtMonteCarlo and 
// AliAnalysisHadEtReconstructed which loop over either Monte Carlo data or 
// real data to get Et

#include "AliAnalysisHadEt.h"
#include "TMath.h"
#include "TList.h"
#include "TH1F.h"
#include "TH2F.h"
#include <iostream>
#include "AliAnalysisEtCuts.h"
#include "AliMCEvent.h"
#include "AliVEvent.h"
#include "AliStack.h"
#include "AliESDtrackCuts.h"
#include "TDatabasePDG.h"
#include "TParticle.h"
#include "Rtypes.h"
#include "AliPDG.h"
#include "AliPIDResponse.h"
#include "AliTPCPIDResponse.h" 
#include "AliInputEventHandler.h"
#include "AliAnalysisManager.h"
#include "AliLog.h"
#include "AliCentrality.h"

using namespace std;

ClassImp(AliAnalysisHadEt);


Int_t AliAnalysisHadEt::fgnumOfEtaBins = 16;
Float_t AliAnalysisHadEt::fgEtaAxis[17]={-0.78, -0.7, -0.58, -0.46, -0.34, -0.22, -0.12, -0.06, -0.0, 0.06, 0.12, 0.22, 0.34, 0.46, 0.58, 0.7, 0.78};
Int_t AliAnalysisHadEt::fgNumOfPtBins = 111;
Float_t AliAnalysisHadEt::fgPtAxis[117]=
  {0.0,0.01,0.02,0.03,0.04, 0.05, 0.06,0.07,0.08,0.09, 0.10,0.11, .12,0.13, .14,0.15, .16,0.17, .18,0.19,
   0.2, .22, .24, .26, .28, 0.30, 0.32, .34, .36, .38, 0.40, .42, .44, .46, .48,
   0.5, .52, .54, .56, .58, 0.60, 0.62, .64, .66, .68, 0.70, .72, .74, .76, .78,
   .80, .82, .84, .86, .88, 0.90, 0.92, .94, .96, .98, 1.00,1.05, 1.1,1.15, 1.2,
  1.25, 1.3,1.35,1.40,1.45, 1.50, 1.55, 1.6,1.65, 1.7, 1.75, 1.8,1.85, 1.9,1.95,
   2.0, 2.2, 2.4, 2.6, 2.8, 3.00, 3.20, 3.4, 3.6, 3.8, 4.00, 4.2, 4.4, 4.6, 4.8,
   5.0, 5.5, 6.0, 6.5, 7.0, 7.50, 8.00, 8.5, 9.0, 9.5, 10.0,12.0,14.0,16.0,18.0,
  20.0,25.0,30.0,35.0,40.0, 45.0, 50.0};
Float_t AliAnalysisHadEt::fgResAxis[81] = {-0.150,-0.140,-0.130,-0.120,-0.110,-0.100,-0.090,-0.080,-0.070,-0.060,
					   -0.050,-0.045,-0.040,-0.035,-0.030,-0.025,-0.024,-0.023,-0.022,-0.021,
					   -0.020,-0.019,-0.018,-0.017,-0.016,-0.015,-0.014,-0.013,-0.012,-0.011,
					   -0.010,-0.009,-0.008,-0.007,-0.006,-0.005,-0.004,-0.003,-0.002,-0.001,
					   -0.000,0.001,0.002,0.003,0.004,0.005,0.006,0.007,0.008,0.009,
					   0.010,0.011,0.012,0.013,0.014,0.015,0.016,0.017,0.018,0.019,
					   0.020,0.021,0.022,0.023,0.024,0.025,0.030,0.035,0.040,0.045,
					   0.050,0.060,0.070,0.080,0.090,0.100,0.110,0.120,0.130,0.140,
					   0.150,};
// Float_t AliAnalysisHadEt::fgResAxis[31] = {-0.15,-0.14,-0.13,-0.12,-0.11,-0.1,-0.09,-0.08,-0.07,-0.06,-0.05,-0.04,-0.03,-0.02,-0.01,0.0,0.01,0.02,0.03,0.04,0.05,0.06,0.07,0.08,0.09,0.1,0.11,0.12,0.13,0.14,0.15};
Int_t AliAnalysisHadEt::fgNumOfResBins = 80;


AliAnalysisHadEt::AliAnalysisHadEt() : AliAnalysisEtCommon()
				     ,fPIDResponse(0)
				     ,fSumEt(0)
				     ,fSumEtAcc(0)
				     ,fTotEt(0)
				     ,fTotEtAcc(0)
				     ,fTotNeutralEt(0)
				     ,fTotNeutralEtAcc(0)
				     ,fTotChargedEt(0)
				     ,fTotChargedEtAcc(0)
				     ,fMultiplicity(0)
				     ,fChargedMultiplicity(0)
				     ,fNeutralMultiplicity(0)
				     ,fhistoList(0)
				     ,fGoodEvent(0)
{//default constructor

}

AliAnalysisHadEt::~AliAnalysisHadEt()
{//destructor
  delete fEsdtrackCutsITSTPC;
  delete fEsdtrackCutsITS;
  delete fEsdtrackCutsTPC;
  delete fPIDResponse;
}

Int_t AliAnalysisHadEt::AnalyseEvent(AliVEvent *event)
{ //this line is basically here to eliminate a compiler warning that event is not used.  Making it a virtual function did not work with the plugin.
  AliAnalysisEtCommon::AnalyseEvent(event);
  ResetEventValues();
  return 0;
}


void AliAnalysisHadEt::Init()
{// clear variables, set up cuts and PDG info
  AliAnalysisEtCommon::Init();
  
}

void AliAnalysisHadEt::ResetEventValues()
{//Resets event values of et to zero
  AliAnalysisEtCommon::ResetEventValues();
  fTotEt = 0;
  fTotEtAcc = 0;
  fTotNeutralEt = 0;
  fTotNeutralEtAcc = 0;
  fTotChargedEt  = 0;
  fTotChargedEtAcc = 0;
  fMultiplicity = 0;
  fChargedMultiplicity = 0;
  fNeutralMultiplicity = 0;
  
}
void AliAnalysisHadEt::CreateEtaPtHisto2D(TString name, TString title)
{     //creates a 2-d histogram in eta and phi and adds it to the list of histograms to be saved
  TString *histoname   = new TString();
  TString *histotitle   = new TString();

  histoname->Append(name);
  histotitle->Append(title);

  TH2F *histo = new TH2F(histoname->Data(),histotitle->Data(),fgNumOfPtBins, fgPtAxis, fgnumOfEtaBins, fgEtaAxis);
  histo->SetYTitle("#eta");
  histo->SetXTitle("p_{T}");
  histo->SetZTitle("E_{T}");
  histo->Sumw2();
  fhistoList->Add(histo);
  delete histoname;
  delete histotitle;
    
}
void AliAnalysisHadEt::CreateResolutionPtHisto2D(TString name, TString title, TString xtitle, TString ytitle)
{     //creates a 2-d histogram in eta and phi and adds it to the list of histograms to be saved
  TString *histoname   = new TString();
  TString *histotitle   = new TString();

  histoname->Append(name);
  histotitle->Append(title);

  //TH2F *histo = new TH2F(histoname->Data(),histotitle->Data(),fgNumOfPtBins, fgPtAxis, fgNumOfResBins, fgResAxis);
  TH2F *histo = new TH2F(histoname->Data(),histotitle->Data(),fgNumOfPtBins, fgPtAxis, fgNumOfResBins, fgResAxis);
  histo->SetYTitle(ytitle);
  histo->SetXTitle(xtitle);
  histo->SetZTitle("Number of entries");
  histo->Sumw2();
  fhistoList->Add(histo);
  delete histoname;
  delete histotitle;
    
}
void AliAnalysisHadEt::CreatePtHisto1D(TString name, TString title, TString xtitle, TString ytitle)
{     //creates a 2-d histogram in eta and phi and adds it to the list of histograms to be saved
  TString *histoname   = new TString();
  TString *histotitle   = new TString();

  histoname->Append(name);
  histotitle->Append(title);

  TH1F *histo = new TH1F(histoname->Data(),histotitle->Data(),fgNumOfPtBins, fgPtAxis);
  histo->SetYTitle(ytitle);
  histo->SetXTitle(xtitle);
  histo->Sumw2();
  fhistoList->Add(histo);
  delete histoname;
  delete histotitle;
    
}

void AliAnalysisHadEt::CreateHisto1D(TString name, TString title, TString xtitle, TString ytitle,Int_t xbins, Float_t xlow,Float_t xhigh)
{     //creates a 1d histogram of the given dimensions and adds it to the list of histograms to be saved
  TString *histoname   = new TString();
  TString *histotitle   = new TString();
  histoname->Append(name);
  histotitle->Append(title);
  TH1F *histo = new TH1F(histoname->Data(),histotitle->Data(),xbins,xlow,xhigh);
  histo->SetYTitle(ytitle);
  histo->SetXTitle(xtitle);
  histo->Sumw2();
  fhistoList->Add(histo);
  delete histoname;
  delete histotitle;
}
void AliAnalysisHadEt::CreateIntHisto1D(TString name, TString title, TString xtitle, TString ytitle,Int_t xbins, Int_t xlow,Int_t xhigh)
{     //creates a 1d integer histogram and adds it to the list of histograms to be saved
  TString *histoname   = new TString();
  TString *histotitle   = new TString();
  histoname->Append(name);
  histotitle->Append(title);
  TH1I *histo = new TH1I(histoname->Data(),histotitle->Data(),xbins,xlow,xhigh);
  histo->SetYTitle(ytitle);
  histo->SetXTitle(xtitle);
  histo->Sumw2();
  fhistoList->Add(histo);
  delete histoname;
  delete histotitle;
    
}
void AliAnalysisHadEt::CreateHisto2D(TString name, TString title, TString xtitle, TString ytitle,Int_t xbins, Float_t xlow,Float_t xhigh,Int_t ybins,Float_t ylow,Float_t yhigh)
{     //creates a 2d histogram and adds it to the list of histograms to be saved
  TString *histoname   = new TString();
  TString *histotitle   = new TString();
  histoname->Append(name);
  histotitle->Append(title);
  TH2F *histo = new TH2F(histoname->Data(),histotitle->Data(),xbins,xlow,xhigh,ybins,ylow,yhigh);
  histo->SetYTitle(ytitle);
  histo->SetXTitle(xtitle);
  histo->Sumw2();
  fhistoList->Add(histo);
  delete histoname;
  delete histotitle;
    
}
void AliAnalysisHadEt::CreateIntHisto2D(TString name, TString title, TString xtitle, TString ytitle,Int_t xbins, Int_t xlow,Int_t xhigh,Int_t ybins,Int_t ylow,Int_t yhigh)
{     //creates a 2-d integer histogram and adds it to the list of histograms to be saved
  TString *histoname   = new TString();
  TString *histotitle   = new TString();
  histoname->Append(name);
  histotitle->Append(title);
  TH2I *histo = new TH2I(histoname->Data(),histotitle->Data(),xbins,xlow,xhigh,ybins,ylow,yhigh);
  histo->SetYTitle(ytitle);
  histo->SetXTitle(xtitle);
  histo->Sumw2();
  fhistoList->Add(histo);
  delete histoname;
  delete histotitle;
    
}

void AliAnalysisHadEt::CreateEtaHisto1D(TString name, TString title)
{     //creates 1d histogram in eta and adds it to the list of histograms to be saved
  TString *histoname   = new TString();
  TString *histotitle   = new TString();
  histoname->Append(name);
  histotitle->Append(title);
  TH1F *histo = new TH1F(histoname->Data(),histotitle->Data(),fgnumOfEtaBins, fgEtaAxis);
  histo->SetYTitle("E_{T}");
  histo->SetXTitle("#eta");
  histo->Sumw2();
  fhistoList->Add(histo);
  delete histoname;
  delete histotitle;
}
void AliAnalysisHadEt::FillHisto1D(TString histname, Float_t x, Float_t weight)
{//fills a 1d histogram with the name histoname with the value x and the weight "weight"
  TH1F     *histo; 
  TString  *name   = new TString();
  name->Append(histname);       
  histo = (TH1F *)fhistoList->FindObject(name->Data()); 
  if(histo){
    histo->Fill((Double_t)x, weight);
  }
  else{cerr<<"CorrelationMaker::FillHisto1D: no histogram "<< name->Data()<<endl;}
  delete name;
}
void AliAnalysisHadEt::FillHisto2D(TString histname, Float_t x, Float_t y, Float_t weight)
{//fills a 2d histogram with the name histoname with the value x and the weight "weight"
  TH2F     *histo; 
  TString  *name   = new TString();
  name->Append(histname);       
  histo = (TH2F *)fhistoList->FindObject(name->Data()); 
  if(histo){
    histo->Fill((Double_t)x,(Double_t)y, weight);
  }
  else{cerr<<"CorrelationMaker::FillHisto2D: no histogram "<< name->Data()<<endl;}
  delete name;
}


Float_t AliAnalysisHadEt::Et(TParticle *part, float mass){//function to calculate et in the same way as it would be calculated in a calorimeter
  if(mass+1000<0.01){//if no mass given return default.  The default argument is -1000
    if(TMath::Abs(part->GetPDG(0)->PdgCode())==2212 || TMath::Abs(part->GetPDG(0)->PdgCode())==2112){
      if(part->GetPDG(0)->PdgCode()==-2212 || part->GetPDG(0)->PdgCode()==-2112){//antiproton or antineutron
	//for antinucleons we specifically want to return the kinetic energy plus twice the rest mass
	return (part->Energy()+part->GetMass())*TMath::Sin(part->Theta());
      }
      if(part->GetPDG(0)->PdgCode()==2212 || part->GetPDG(0)->PdgCode()==2112){//proton or neutron
	//for nucleons we specifically want to return the kinetic energy only
	return (part->Energy()-part->GetMass())*TMath::Sin(part->Theta());
      }
    }
    else{//otherwise go to the default
      return part->Energy()*TMath::Sin(part->Theta());
    }
  }
  else{//otherwise use the mass that was given
    return (TMath::Sqrt(TMath::Power(part->P(),2.0)+TMath::Power(mass,2.0)))*TMath::Sin(part->Theta());
  }
  return 0.0;
}
Float_t AliAnalysisHadEt::Et(Float_t p, Float_t theta, Int_t pid, Short_t charge) const {//function to calculate et in the same way as it would be calculated in a calorimeter
  if(pid==AliAnalysisEtCommon::fgPiPlusCode || pid==AliAnalysisEtCommon::fgPiMinusCode){//Nothing special for pions
    return TMath::Sqrt(p*p + fgPionMass*fgPionMass) * TMath::Sin(theta);
  }
  if(pid==AliAnalysisEtCommon::fgKPlusCode || pid==AliAnalysisEtCommon::fgKMinusCode){//Nothing special for kaons
    return TMath::Sqrt(p*p + AliAnalysisEtCommon::fgKaonMass*AliAnalysisEtCommon::fgKaonMass) * TMath::Sin(theta);
  }
  if(pid==AliAnalysisEtCommon::fgEPlusCode || pid==AliAnalysisEtCommon::fgEMinusCode){//Nothing special for electrons
    return TMath::Sqrt(p*p + AliAnalysisEtCommon::fgElectronMass*AliAnalysisEtCommon::fgElectronMass) * TMath::Sin(theta);
  }
  if(pid==AliAnalysisEtCommon::fgProtonCode || pid==AliAnalysisEtCommon::fgAntiProtonCode){//But for protons we must be careful...
    if(charge<0.0){//antiprotns: kinetic energy plus twice the rest mass
      return (TMath::Sqrt(p*p + AliAnalysisEtCommon::fgProtonMass*AliAnalysisEtCommon::fgProtonMass) + AliAnalysisEtCommon::fgProtonMass) * TMath::Sin(theta);
    }
    if(charge>0.0){//antiprotns: kinetic energy only
      return (TMath::Sqrt(p*p + AliAnalysisEtCommon::fgProtonMass*AliAnalysisEtCommon::fgProtonMass) - AliAnalysisEtCommon::fgProtonMass) * TMath::Sin(theta);
    }
  }
  //cerr<<"Uh-oh!  Et not set properly!"<<endl;
  return 0.0;
}

Float_t AliAnalysisHadEt::TrueP(float pTrec) const {
  if(pTrec>1.0) return pTrec;
  return pTrec/(1-599.334*pTrec+7285.15*pTrec*pTrec)+pTrec;
}

Int_t AliAnalysisHadEt::GetCentralityBin(Int_t numberofbins,AliCentrality *centrality){
  Int_t centralitybin = -1;
  if(numberofbins<21) centralitybin= centrality->GetCentralityClass10(fCentralityMethod);
  else{
    if(numberofbins<41) centralitybin= centrality->GetCentralityClass5(fCentralityMethod);
    else{
      Float_t centpercent = centrality->GetCentralityPercentile(fCentralityMethod);
      centralitybin= centrality->GetCentralityClass5(fCentralityMethod);
      Float_t centralitybinfloat = (centpercent/2.5);
      if(centralitybin>=0){
	centralitybin =(Int_t) (centralitybinfloat);
	//cout<<" centbin "<<centralitybin<<" centrality "<<centpercent<<" "<<centpercent/2.5 <<endl;
      }
    }
  }
  //cout<<" centrality bin "<<centralitybin<<endl;
  return centralitybin;
}
Commit	Line	Data
cf6522d1	1	//Created by Christine Nattrass, Rebecca Scott, Irakli Martashvili
641e1e0c	2	//University of Tennessee at Knoxville
cf6522d1	3	//
	4	// This class is designed for the analysis of the hadronic component of
	5	// transverse energy. It is used by AliAnalysisTaskHadEt.
	6	// This gets information about the hadronic component of the transverse energy
	7	// from tracks reconstructed in an event
	8	// it has daughters, AliAnalysisHadEtMonteCarlo and
	9	// AliAnalysisHadEtReconstructed which loop over either Monte Carlo data or
	10	// real data to get Et
3ce6b879	11
641e1e0c	12	#include "AliAnalysisHadEt.h"
	13	#include "TMath.h"
	14	#include "TList.h"
	15	#include "TH1F.h"
	16	#include "TH2F.h"
	17	#include <iostream>
	18	#include "AliAnalysisEtCuts.h"
020abff0	19	#include "AliMCEvent.h"
641e1e0c	20	#include "AliVEvent.h"
020abff0	21	#include "AliStack.h"
	22	#include "AliESDtrackCuts.h"
	23	#include "TDatabasePDG.h"
	24	#include "TParticle.h"
	25	#include "Rtypes.h"
4b40b2b1	26	#include "AliPDG.h"
9a365626	27	#include "AliPIDResponse.h"
	28	#include "AliTPCPIDResponse.h"
	29	#include "AliInputEventHandler.h"
	30	#include "AliAnalysisManager.h"
	31	#include "AliLog.h"
c692b181	32	#include "AliCentrality.h"
641e1e0c	33
	34	using namespace std;
	35
	36	ClassImp(AliAnalysisHadEt);
	37
	38
4b40b2b1	39	Int_t AliAnalysisHadEt::fgnumOfEtaBins = 16;
4b40b2b1	40	Float_t AliAnalysisHadEt::fgEtaAxis[17]={-0.78, -0.7, -0.58, -0.46, -0.34, -0.22, -0.12, -0.06, -0.0, 0.06, 0.12, 0.22, 0.34, 0.46, 0.58, 0.7, 0.78};
020abff0	41	Int_t AliAnalysisHadEt::fgNumOfPtBins = 111;
020abff0	42	Float_t AliAnalysisHadEt::fgPtAxis[117]=
641e1e0c	43	{0.0,0.01,0.02,0.03,0.04, 0.05, 0.06,0.07,0.08,0.09, 0.10,0.11, .12,0.13, .14,0.15, .16,0.17, .18,0.19,
	44	0.2, .22, .24, .26, .28, 0.30, 0.32, .34, .36, .38, 0.40, .42, .44, .46, .48,
	45	0.5, .52, .54, .56, .58, 0.60, 0.62, .64, .66, .68, 0.70, .72, .74, .76, .78,
	46	.80, .82, .84, .86, .88, 0.90, 0.92, .94, .96, .98, 1.00,1.05, 1.1,1.15, 1.2,
	47	1.25, 1.3,1.35,1.40,1.45, 1.50, 1.55, 1.6,1.65, 1.7, 1.75, 1.8,1.85, 1.9,1.95,
	48	2.0, 2.2, 2.4, 2.6, 2.8, 3.00, 3.20, 3.4, 3.6, 3.8, 4.00, 4.2, 4.4, 4.6, 4.8,
	49	5.0, 5.5, 6.0, 6.5, 7.0, 7.50, 8.00, 8.5, 9.0, 9.5, 10.0,12.0,14.0,16.0,18.0,
66e4ff36	50	20.0,25.0,30.0,35.0,40.0, 45.0, 50.0};
	51	Float_t AliAnalysisHadEt::fgResAxis[81] = {-0.150,-0.140,-0.130,-0.120,-0.110,-0.100,-0.090,-0.080,-0.070,-0.060,
	52	-0.050,-0.045,-0.040,-0.035,-0.030,-0.025,-0.024,-0.023,-0.022,-0.021,
	53	-0.020,-0.019,-0.018,-0.017,-0.016,-0.015,-0.014,-0.013,-0.012,-0.011,
	54	-0.010,-0.009,-0.008,-0.007,-0.006,-0.005,-0.004,-0.003,-0.002,-0.001,
	55	-0.000,0.001,0.002,0.003,0.004,0.005,0.006,0.007,0.008,0.009,
	56	0.010,0.011,0.012,0.013,0.014,0.015,0.016,0.017,0.018,0.019,
	57	0.020,0.021,0.022,0.023,0.024,0.025,0.030,0.035,0.040,0.045,
	58	0.050,0.060,0.070,0.080,0.090,0.100,0.110,0.120,0.130,0.140,
	59	0.150,};
	60	// Float_t AliAnalysisHadEt::fgResAxis[31] = {-0.15,-0.14,-0.13,-0.12,-0.11,-0.1,-0.09,-0.08,-0.07,-0.06,-0.05,-0.04,-0.03,-0.02,-0.01,0.0,0.01,0.02,0.03,0.04,0.05,0.06,0.07,0.08,0.09,0.1,0.11,0.12,0.13,0.14,0.15};
	61	Int_t AliAnalysisHadEt::fgNumOfResBins = 80;
641e1e0c	62
d6214a64	63
7d2d1773	64	AliAnalysisHadEt::AliAnalysisHadEt() : AliAnalysisEtCommon()
9a365626	65	,fPIDResponse(0)
7d2d1773	66	,fSumEt(0)
	67	,fSumEtAcc(0)
	68	,fTotEt(0)
	69	,fTotEtAcc(0)
	70	,fTotNeutralEt(0)
	71	,fTotNeutralEtAcc(0)
	72	,fTotChargedEt(0)
	73	,fTotChargedEtAcc(0)
	74	,fMultiplicity(0)
	75	,fChargedMultiplicity(0)
	76	,fNeutralMultiplicity(0)
	77	,fhistoList(0)
66e4ff36	78	,fGoodEvent(0)
020abff0	79	{//default constructor
641e1e0c	80
	81	}
	82
	83	AliAnalysisHadEt::~AliAnalysisHadEt()
020abff0	84	{//destructor
951efd81	85	delete fEsdtrackCutsITSTPC;
	86	delete fEsdtrackCutsITS;
	87	delete fEsdtrackCutsTPC;
9a365626	88	delete fPIDResponse;
641e1e0c	89	}
	90
	91	Int_t AliAnalysisHadEt::AnalyseEvent(AliVEvent *event)
020abff0	92	{ //this line is basically here to eliminate a compiler warning that event is not used. Making it a virtual function did not work with the plugin.
7d2d1773	93	AliAnalysisEtCommon::AnalyseEvent(event);
0065531a	94	ResetEventValues();
641e1e0c	95	return 0;
	96	}
	97
641e1e0c	98
641e1e0c	99	void AliAnalysisHadEt::Init()
0065531a	100	{// clear variables, set up cuts and PDG info
7d2d1773	101	AliAnalysisEtCommon::Init();
9a365626	102
641e1e0c	103	}
	104
	105	void AliAnalysisHadEt::ResetEventValues()
020abff0	106	{//Resets event values of et to zero
7d2d1773	107	AliAnalysisEtCommon::ResetEventValues();
0065531a	108	fTotEt = 0;
	109	fTotEtAcc = 0;
	110	fTotNeutralEt = 0;
	111	fTotNeutralEtAcc = 0;
	112	fTotChargedEt = 0;
	113	fTotChargedEtAcc = 0;
	114	fMultiplicity = 0;
	115	fChargedMultiplicity = 0;
	116	fNeutralMultiplicity = 0;
	117
4998becf	118	}
641e1e0c	119	void AliAnalysisHadEt::CreateEtaPtHisto2D(TString name, TString title)
020abff0	120	{ //creates a 2-d histogram in eta and phi and adds it to the list of histograms to be saved
641e1e0c	121	TString *histoname = new TString();
	122	TString *histotitle = new TString();
	123
	124	histoname->Append(name);
	125	histotitle->Append(title);
641e1e0c	126
020abff0	127	TH2F *histo = new TH2F(histoname->Data(),histotitle->Data(),fgNumOfPtBins, fgPtAxis, fgnumOfEtaBins, fgEtaAxis);
641e1e0c	128	histo->SetYTitle("#eta");
	129	histo->SetXTitle("p_{T}");
	130	histo->SetZTitle("E_{T}");
	131	histo->Sumw2();
020abff0	132	fhistoList->Add(histo);
641e1e0c	133	delete histoname;
	134	delete histotitle;
	135
66e4ff36	136	}
	137	void AliAnalysisHadEt::CreateResolutionPtHisto2D(TString name, TString title, TString xtitle, TString ytitle)
	138	{ //creates a 2-d histogram in eta and phi and adds it to the list of histograms to be saved
	139	TString *histoname = new TString();
	140	TString *histotitle = new TString();
	141
	142	histoname->Append(name);
	143	histotitle->Append(title);
	144
	145	//TH2F *histo = new TH2F(histoname->Data(),histotitle->Data(),fgNumOfPtBins, fgPtAxis, fgNumOfResBins, fgResAxis);
	146	TH2F *histo = new TH2F(histoname->Data(),histotitle->Data(),fgNumOfPtBins, fgPtAxis, fgNumOfResBins, fgResAxis);
	147	histo->SetYTitle(ytitle);
	148	histo->SetXTitle(xtitle);
	149	histo->SetZTitle("Number of entries");
	150	histo->Sumw2();
	151	fhistoList->Add(histo);
	152	delete histoname;
	153	delete histotitle;
	154
	155	}
	156	void AliAnalysisHadEt::CreatePtHisto1D(TString name, TString title, TString xtitle, TString ytitle)
	157	{ //creates a 2-d histogram in eta and phi and adds it to the list of histograms to be saved
	158	TString *histoname = new TString();
	159	TString *histotitle = new TString();
	160
	161	histoname->Append(name);
	162	histotitle->Append(title);
	163
	164	TH1F *histo = new TH1F(histoname->Data(),histotitle->Data(),fgNumOfPtBins, fgPtAxis);
	165	histo->SetYTitle(ytitle);
	166	histo->SetXTitle(xtitle);
	167	histo->Sumw2();
	168	fhistoList->Add(histo);
	169	delete histoname;
	170	delete histotitle;
	171
641e1e0c	172	}
	173
	174	void AliAnalysisHadEt::CreateHisto1D(TString name, TString title, TString xtitle, TString ytitle,Int_t xbins, Float_t xlow,Float_t xhigh)
020abff0	175	{ //creates a 1d histogram of the given dimensions and adds it to the list of histograms to be saved
641e1e0c	176	TString *histoname = new TString();
641e1e0c	177	TString *histotitle = new TString();
641e1e0c	178	histoname->Append(name);
641e1e0c	179	histotitle->Append(title);
641e1e0c	180	TH1F *histo = new TH1F(histoname->Data(),histotitle->Data(),xbins,xlow,xhigh);
	181	histo->SetYTitle(ytitle);
	182	histo->SetXTitle(xtitle);
	183	histo->Sumw2();
020abff0	184	fhistoList->Add(histo);
641e1e0c	185	delete histoname;
641e1e0c	186	delete histotitle;
641e1e0c	187	}
641e1e0c	188	void AliAnalysisHadEt::CreateIntHisto1D(TString name, TString title, TString xtitle, TString ytitle,Int_t xbins, Int_t xlow,Int_t xhigh)
020abff0	189	{ //creates a 1d integer histogram and adds it to the list of histograms to be saved
641e1e0c	190	TString *histoname = new TString();
641e1e0c	191	TString *histotitle = new TString();
641e1e0c	192	histoname->Append(name);
641e1e0c	193	histotitle->Append(title);
641e1e0c	194	TH1I *histo = new TH1I(histoname->Data(),histotitle->Data(),xbins,xlow,xhigh);
	195	histo->SetYTitle(ytitle);
	196	histo->SetXTitle(xtitle);
	197	histo->Sumw2();
020abff0	198	fhistoList->Add(histo);
641e1e0c	199	delete histoname;
	200	delete histotitle;
	201
	202	}
	203	void AliAnalysisHadEt::CreateHisto2D(TString name, TString title, TString xtitle, TString ytitle,Int_t xbins, Float_t xlow,Float_t xhigh,Int_t ybins,Float_t ylow,Float_t yhigh)
020abff0	204	{ //creates a 2d histogram and adds it to the list of histograms to be saved
641e1e0c	205	TString *histoname = new TString();
641e1e0c	206	TString *histotitle = new TString();
641e1e0c	207	histoname->Append(name);
641e1e0c	208	histotitle->Append(title);
641e1e0c	209	TH2F *histo = new TH2F(histoname->Data(),histotitle->Data(),xbins,xlow,xhigh,ybins,ylow,yhigh);
	210	histo->SetYTitle(ytitle);
	211	histo->SetXTitle(xtitle);
	212	histo->Sumw2();
020abff0	213	fhistoList->Add(histo);
641e1e0c	214	delete histoname;
	215	delete histotitle;
	216
	217	}
	218	void AliAnalysisHadEt::CreateIntHisto2D(TString name, TString title, TString xtitle, TString ytitle,Int_t xbins, Int_t xlow,Int_t xhigh,Int_t ybins,Int_t ylow,Int_t yhigh)
020abff0	219	{ //creates a 2-d integer histogram and adds it to the list of histograms to be saved
641e1e0c	220	TString *histoname = new TString();
641e1e0c	221	TString *histotitle = new TString();
641e1e0c	222	histoname->Append(name);
641e1e0c	223	histotitle->Append(title);
641e1e0c	224	TH2I *histo = new TH2I(histoname->Data(),histotitle->Data(),xbins,xlow,xhigh,ybins,ylow,yhigh);
	225	histo->SetYTitle(ytitle);
	226	histo->SetXTitle(xtitle);
	227	histo->Sumw2();
020abff0	228	fhistoList->Add(histo);
641e1e0c	229	delete histoname;
	230	delete histotitle;
	231
	232	}
	233
	234	void AliAnalysisHadEt::CreateEtaHisto1D(TString name, TString title)
020abff0	235	{ //creates 1d histogram in eta and adds it to the list of histograms to be saved
641e1e0c	236	TString *histoname = new TString();
641e1e0c	237	TString *histotitle = new TString();
641e1e0c	238	histoname->Append(name);
641e1e0c	239	histotitle->Append(title);
020abff0	240	TH1F *histo = new TH1F(histoname->Data(),histotitle->Data(),fgnumOfEtaBins, fgEtaAxis);
641e1e0c	241	histo->SetYTitle("E_{T}");
	242	histo->SetXTitle("#eta");
	243	histo->Sumw2();
020abff0	244	fhistoList->Add(histo);
641e1e0c	245	delete histoname;
641e1e0c	246	delete histotitle;
641e1e0c	247	}
641e1e0c	248	void AliAnalysisHadEt::FillHisto1D(TString histname, Float_t x, Float_t weight)
020abff0	249	{//fills a 1d histogram with the name histoname with the value x and the weight "weight"
641e1e0c	250	TH1F *histo;
641e1e0c	251	TString *name = new TString();
641e1e0c	252	name->Append(histname);
020abff0	253	histo = (TH1F *)fhistoList->FindObject(name->Data());
641e1e0c	254	if(histo){
	255	histo->Fill((Double_t)x, weight);
	256	}
	257	else{cerr<<"CorrelationMaker::FillHisto1D: no histogram "<< name->Data()<<endl;}
	258	delete name;
	259	}
	260	void AliAnalysisHadEt::FillHisto2D(TString histname, Float_t x, Float_t y, Float_t weight)
020abff0	261	{//fills a 2d histogram with the name histoname with the value x and the weight "weight"
641e1e0c	262	TH2F *histo;
641e1e0c	263	TString *name = new TString();
641e1e0c	264	name->Append(histname);
020abff0	265	histo = (TH2F *)fhistoList->FindObject(name->Data());
641e1e0c	266	if(histo){
	267	histo->Fill((Double_t)x,(Double_t)y, weight);
	268	}
	269	else{cerr<<"CorrelationMaker::FillHisto2D: no histogram "<< name->Data()<<endl;}
	270	delete name;
	271	}
	272
	273
020abff0	274	Float_t AliAnalysisHadEt::Et(TParticle *part, float mass){//function to calculate et in the same way as it would be calculated in a calorimeter
020abff0	275	if(mass+1000<0.01){//if no mass given return default. The default argument is -1000
641e1e0c	276	if(TMath::Abs(part->GetPDG(0)->PdgCode())==2212 \|\| TMath::Abs(part->GetPDG(0)->PdgCode())==2112){
	277	if(part->GetPDG(0)->PdgCode()==-2212 \|\| part->GetPDG(0)->PdgCode()==-2112){//antiproton or antineutron
	278	//for antinucleons we specifically want to return the kinetic energy plus twice the rest mass
	279	return (part->Energy()+part->GetMass())*TMath::Sin(part->Theta());
	280	}
0e866ddc	281	if(part->GetPDG(0)->PdgCode()==2212 \|\| part->GetPDG(0)->PdgCode()==2112){//proton or neutron
641e1e0c	282	//for nucleons we specifically want to return the kinetic energy only
	283	return (part->Energy()-part->GetMass())*TMath::Sin(part->Theta());
	284	}
	285	}
	286	else{//otherwise go to the default
	287	return part->Energy()*TMath::Sin(part->Theta());
	288	}
	289	}
	290	else{//otherwise use the mass that was given
	291	return (TMath::Sqrt(TMath::Power(part->P(),2.0)+TMath::Power(mass,2.0)))*TMath::Sin(part->Theta());
	292	}
	293	return 0.0;
	294	}
464aa50c	295	Float_t AliAnalysisHadEt::Et(Float_t p, Float_t theta, Int_t pid, Short_t charge) const {//function to calculate et in the same way as it would be calculated in a calorimeter
7d2d1773	296	if(pid==AliAnalysisEtCommon::fgPiPlusCode \|\| pid==AliAnalysisEtCommon::fgPiMinusCode){//Nothing special for pions
4b40b2b1	297	return TMath::Sqrt(pp + fgPionMassfgPionMass) * TMath::Sin(theta);
3ce6b879	298	}
7d2d1773	299	if(pid==AliAnalysisEtCommon::fgKPlusCode \|\| pid==AliAnalysisEtCommon::fgKMinusCode){//Nothing special for kaons
7d2d1773	300	return TMath::Sqrt(pp + AliAnalysisEtCommon::fgKaonMassAliAnalysisEtCommon::fgKaonMass) * TMath::Sin(theta);
3ce6b879	301	}
7d2d1773	302	if(pid==AliAnalysisEtCommon::fgEPlusCode \|\| pid==AliAnalysisEtCommon::fgEMinusCode){//Nothing special for electrons
7d2d1773	303	return TMath::Sqrt(pp + AliAnalysisEtCommon::fgElectronMassAliAnalysisEtCommon::fgElectronMass) * TMath::Sin(theta);
3ce6b879	304	}
7d2d1773	305	if(pid==AliAnalysisEtCommon::fgProtonCode \|\| pid==AliAnalysisEtCommon::fgAntiProtonCode){//But for protons we must be careful...
3ce6b879	306	if(charge<0.0){//antiprotns: kinetic energy plus twice the rest mass
7d2d1773	307	return (TMath::Sqrt(pp + AliAnalysisEtCommon::fgProtonMassAliAnalysisEtCommon::fgProtonMass) + AliAnalysisEtCommon::fgProtonMass) * TMath::Sin(theta);
3ce6b879	308	}
3ce6b879	309	if(charge>0.0){//antiprotns: kinetic energy only
7d2d1773	310	return (TMath::Sqrt(pp + AliAnalysisEtCommon::fgProtonMassAliAnalysisEtCommon::fgProtonMass) - AliAnalysisEtCommon::fgProtonMass) * TMath::Sin(theta);
3ce6b879	311	}
3ce6b879	312	}
66e4ff36	313	//cerr<<"Uh-oh! Et not set properly!"<<endl;
3ce6b879	314	return 0.0;
3ce6b879	315	}
7d2d1773	316
d263b18f	317	Float_t AliAnalysisHadEt::TrueP(float pTrec) const {
66e4ff36	318	if(pTrec>1.0) return pTrec;
	319	return pTrec/(1-599.334pTrec+7285.15pTrec*pTrec)+pTrec;
	320	}
c692b181	321
	322	Int_t AliAnalysisHadEt::GetCentralityBin(Int_t numberofbins,AliCentrality *centrality){
	323	Int_t centralitybin = -1;
	324	if(numberofbins<21) centralitybin= centrality->GetCentralityClass10(fCentralityMethod);
	325	else{
	326	if(numberofbins<41) centralitybin= centrality->GetCentralityClass5(fCentralityMethod);
	327	else{
	328	Float_t centpercent = centrality->GetCentralityPercentile(fCentralityMethod);
	329	centralitybin= centrality->GetCentralityClass5(fCentralityMethod);
4d0983d4	330	Float_t centralitybinfloat = (centpercent/2.5);
	331	if(centralitybin>=0){
	332	centralitybin =(Int_t) (centralitybinfloat);
	333	//cout<<" centbin "<<centralitybin<<" centrality "<<centpercent<<" "<<centpercent/2.5 <<endl;
	334	}
c692b181	335	}
	336	}
	337	//cout<<" centrality bin "<<centralitybin<<endl;
	338	return centralitybin;
	339	}