AliAnalysisTaskEmcalJetSpectraMECpA::AliAnalysisTaskEmcalJetSpectraMECpA() :
AliAnalysisTaskEmcalJet("spectra",kFALSE),
fHistRhovsCent(0),
AliAnalysisTaskEmcalJetSpectraMECpA::AliAnalysisTaskEmcalJetSpectraMECpA() :
AliAnalysisTaskEmcalJet("spectra",kFALSE),
fHistRhovsCent(0),
AliAnalysisTaskEmcalJetSpectraMECpA::AliAnalysisTaskEmcalJetSpectraMECpA(const char *name) :
AliAnalysisTaskEmcalJet(name,kTRUE),
fHistRhovsCent(0),
AliAnalysisTaskEmcalJetSpectraMECpA::AliAnalysisTaskEmcalJetSpectraMECpA(const char *name) :
AliAnalysisTaskEmcalJet(name,kTRUE),
fHistRhovsCent(0),
AliAnalysisTaskEmcalJet::UserCreateOutputObjects();
fHistRhovsCent = new TH2F("RhovsCent", "RhovsCent", 100, 0.0, 100.0, 500, 0, 500);
AliAnalysisTaskEmcalJet::UserCreateOutputObjects();
fHistRhovsCent = new TH2F("RhovsCent", "RhovsCent", 100, 0.0, 100.0, 500, 0, 500);
- fHistRhoScvsCent = new TH2F("RhoScvsCent", "RhoScvsCent", 100, 0.0, 100.0, 500, 0, 500);
+ fHistRhoScvsCent = new TH2F("RhoScvsCent", "RhoScvsCent", 100, 0.0, 100.0, 500, 0, 500);
fHistNjetvsCent = new TH2F("NjetvsCent", "NjetvsCent", 100, 0.0, 100.0, 100, 0, 100);
TString name;
TString title;
fHistNjetvsCent = new TH2F("NjetvsCent", "NjetvsCent", 100, 0.0, 100.0, 100, 0, 100);
TString name;
TString title;
name = TString(Form("JetPtvsTrackPt_%i",i));
title = TString(Form("Jet pT vs Leading Track pT cent bin %i",i));
fHistJetPtvsTrackPt[i] = new TH2F(name,title,1000,-500,500,100,0,100);
name = TString(Form("JetPtvsTrackPt_%i",i));
title = TString(Form("Jet pT vs Leading Track pT cent bin %i",i));
fHistJetPtvsTrackPt[i] = new TH2F(name,title,1000,-500,500,100,0,100);
- const Int_t nTrack = fTracks->GetEntriesFast();
- for (int i = 0;i<nTrack;i++){
- AliVParticle *track = static_cast<AliVParticle*>(fTracks->At(i));
- if (!track)
- continue;
- fHistTrackPt[centbin]->Fill(track->Pt());
+ if (fTracks){
+ const Int_t nTrack = fTracks->GetEntriesFast();
+ for (int i = 0;i<nTrack;i++){
+ AliVParticle *track = static_cast<AliVParticle*>(fTracks->At(i));
+ if (!track)
+ continue;
+ fHistTrackPt[centbin]->Fill(track->Pt());
+ }
fHistEP0[centbin]->Fill(fEPV0);
fHistEP0A[centbin]->Fill(fEPV0A);
fHistEP0C[centbin]->Fill(fEPV0C);
fHistEP0[centbin]->Fill(fEPV0);
fHistEP0A[centbin]->Fill(fEPV0A);
fHistEP0C[centbin]->Fill(fEPV0C);
fRho = GetRhoFromEvent(fRhoName);
fRhoVal = fRho->GetVal();
fHistRhovsCent->Fill(fCent,fRhoVal);
fRho = GetRhoFromEvent(fRhoName);
fRhoVal = fRho->GetVal();
fHistRhovsCent->Fill(fCent,fRhoVal);
Double_t jetPt = -500;
jetPt = jet->Pt()-jet->Area()*fRhoVal;
fHistJetPtvsTrackPt[centbin]->Fill(jetPt,jet->MaxTrackPt());
Double_t jetPt = -500;
jetPt = jet->Pt()-jet->Area()*fRhoVal;
fHistJetPtvsTrackPt[centbin]->Fill(jetPt,jet->MaxTrackPt());
if(fRhoScVal) {
Double_t jetPtSc = -500;
jetPtSc = jet->Pt()-jet->Area()*fRhoScVal;
fHistJetPtScvsTrackPt[centbin]->Fill(jetPtSc,jet->MaxTrackPt());
if(fRhoScVal) {
Double_t jetPtSc = -500;
jetPtSc = jet->Pt()-jet->Area()*fRhoScVal;
fHistJetPtScvsTrackPt[centbin]->Fill(jetPtSc,jet->MaxTrackPt());