#include <TString.h>
#include <TChain.h>
#include <TTree.h>
-#include <TH1D.h>
-#include <TH2D.h>
-#include <TH3D.h>
+#include <TH1F.h>
+#include <TH2F.h>
+#include <TH3F.h>
+#include <THnSparse.h>
#include <TCanvas.h>
#include <TList.h>
#include <TLorentzVector.h>
#include <TClonesArray.h>
#include <TObjArray.h>
-#include "AliAnalysisTaskSE.h"
+#include "AliAnalysisTaskEmcal.h"
+#include "AliAnalysisTaskEmcalJet.h"
#include "AliAnalysisManager.h"
#include "AliStack.h"
#include "AliESDtrackCuts.h"
#include "AliEMCALRecoUtils.h"
#include "AliVCaloCells.h"
#include "AliPicoTrack.h"
+#include "AliAnalysisTaskEmcal.h"
+#include "AliAnalysisTaskEmcalJet.h"
#include "Rtypes.h"
ClassImp(AliAnalysisTaskFullpAJets)
//________________________________________________________________________
AliAnalysisTaskFullpAJets::AliAnalysisTaskFullpAJets() :
- AliAnalysisTaskSE(),
+ AliAnalysisTaskEmcalJet(),
fOutput(0),
+ flTrack(0),
+ flCluster(0),
fhTrackPt(0),
fhTrackEta(0),
fhTrackPhi(0),
fhCentrality(0),
fhEMCalCellCounts(0),
+ fhChargeAndNeutralEvents(0),
+ fhChargeOnlyEvents(0),
+ fhNeutralOnlyEvents(0),
+ fhNothingEvents(0),
+ fhEMCalChargeAndNeutralEvents(0),
+ fhEMCalChargeOnlyEvents(0),
+ fhEMCalNeutralOnlyEvents(0),
+ fhEMCalNothingEvents(0),
+
fhTrackEtaPhi(0),
fhTrackPhiPt(0),
fhTrackEtaPt(0),
fhClusterEtaPhi(0),
fhClusterPhiPt(0),
fhClusterEtaPt(0),
- fhRhoScale(0),
-
- fhTrackEtaPhiPt(0),
- fhGlobalTrackEtaPhiPt(0),
- fhComplementaryTrackEtaPhiPt(0),
- fhClusterEtaPhiPt(0),
+ fhEMCalEventMult(0),
+ fhTPCEventMult(0),
+ fhEMCalTrackEventMult(0),
fpEMCalEventMult(0),
fpTPCEventMult(0),
- fpRhoScale(0),
fpTrackPtProfile(0),
fpClusterPtProfile(0),
- //fTPCRawJets(0),
- //fEMCalRawJets(0),
- //fRhoFull0(0),
- //fRhoFull1(0),
- //fRhoFull2(0),
- //fRhoFullN(0),
- //fRhoFullDijet(0),
- //fRhoFullkT(0),
- //fRhoFullCMS(0),
- //fRhoCharged0(0),
- //fRhoCharged1(0),
- //fRhoCharged2(0),
- //fRhoChargedN(0),
- //fRhoChargedkT(0),
- //fRhoChargedkTScale(0),
- //fRhoChargedCMS(0),
- fRhoChargedScale(0),
+ fpFullJetEDProfile(0),
+ fpChargedJetEDProfile(0),
+ fpChargedJetEDProfileScaled(0),
+
+ fTPCRawJets(0),
+ fEMCalRawJets(0),
fRhoChargedCMSScale(0),
+ fRhoChargedScale(0),
+ fRhoFull0(0),
+ fRhoFull1(0),
+ fRhoFull2(0),
+ fRhoFullN(0),
+ fRhoFullDijet(0),
+ fRhoFullkT(0),
+ fRhoFullCMS(0),
+ fRhoCharged0(0),
+ fRhoCharged1(0),
+ fRhoCharged2(0),
+ fRhoChargedN(0),
+ fRhoChargedkT(0),
+ fRhoChargedkTScale(0),
+ fRhoChargedCMS(0),
fTPCJet(0),
fTPCFullJet(0),
fEMCALGeometry(0),
fCells(0),
fDoNEF(0),
+ fDoNEFSignalOnly(1),
+ fSignalTrackBias(0),
+ fTrackQA(0),
+ fClusterQA(0),
+ fCalculateRhoJet(0),
+ fDoVertexRCut(0),
+ fMCPartLevel(0),
+ fDoTHnSparse(0),
+ fDoJetRhoDensity(0),
fEMCalPhiMin(1.39626),
fEMCalPhiMax(3.26377),
fEMCalPhiTotal(1.86750),
fParticlePtUp(200.0),
fParticlePtBins(200),
fJetR(0.4),
+ fJetRAccept(0.4),
+ fFullEDJetR(0.6),
+ fChargedEDJetR(0.8),
fJetRForRho(0.5),
fJetAreaCutFrac(0.6),
fJetAreaThreshold(0.30159),
fnEMCalCells(12288),
- fScaleFactor(1.50),
+ fScaleFactor(1.28),
fNColl(7),
fTrackMinPt(0.15),
fClusterMinPt(0.3),
fRhoFull(0),
fRhoCharged(0),
fnTracks(0),
+ fnEMCalTracks(0),
fnClusters(0),
fnCaloClusters(0),
fnAKTFullJets(0),
//________________________________________________________________________
AliAnalysisTaskFullpAJets::AliAnalysisTaskFullpAJets(const char *name) :
- AliAnalysisTaskSE(name),
+ AliAnalysisTaskEmcalJet(name),
fOutput(0),
+ flTrack(0),
+ flCluster(0),
fhTrackPt(0),
fhTrackEta(0),
fhTrackPhi(0),
fhCentrality(0),
fhEMCalCellCounts(0),
+ fhChargeAndNeutralEvents(0),
+ fhChargeOnlyEvents(0),
+ fhNeutralOnlyEvents(0),
+ fhNothingEvents(0),
+ fhEMCalChargeAndNeutralEvents(0),
+ fhEMCalChargeOnlyEvents(0),
+ fhEMCalNeutralOnlyEvents(0),
+ fhEMCalNothingEvents(0),
+
fhTrackEtaPhi(0),
fhTrackPhiPt(0),
fhTrackEtaPt(0),
fhClusterEtaPhi(0),
fhClusterPhiPt(0),
fhClusterEtaPt(0),
- fhRhoScale(0),
-
- fhTrackEtaPhiPt(0),
- fhGlobalTrackEtaPhiPt(0),
- fhComplementaryTrackEtaPhiPt(0),
- fhClusterEtaPhiPt(0),
+ fhEMCalEventMult(0),
+ fhTPCEventMult(0),
+ fhEMCalTrackEventMult(0),
fpEMCalEventMult(0),
fpTPCEventMult(0),
- fpRhoScale(0),
fpTrackPtProfile(0),
fpClusterPtProfile(0),
+ fpFullJetEDProfile(0),
+ fpChargedJetEDProfile(0),
+ fpChargedJetEDProfileScaled(0),
+
fTPCRawJets(0),
fEMCalRawJets(0),
- //fRhoFull0(0),
- //fRhoFull1(0),
- //fRhoFull2(0),
- //fRhoFullN(0),
- //fRhoFullDijet(0),
- //fRhoFullkT(0),
- //fRhoFullCMS(0),
- //fRhoCharged0(0),
- //fRhoCharged1(0),
- //fRhoCharged2(0),
- //fRhoChargedN(0),
- //fRhoChargedkT(0),
- //fRhoChargedkTScale(0),
- //fRhoChargedCMS(0),
- fRhoChargedScale(0),
fRhoChargedCMSScale(0),
+ fRhoChargedScale(0),
+ fRhoFull0(0),
+ fRhoFull1(0),
+ fRhoFull2(0),
+ fRhoFullN(0),
+ fRhoFullDijet(0),
+ fRhoFullkT(0),
+ fRhoFullCMS(0),
+ fRhoCharged0(0),
+ fRhoCharged1(0),
+ fRhoCharged2(0),
+ fRhoChargedN(0),
+ fRhoChargedkT(0),
+ fRhoChargedkTScale(0),
+ fRhoChargedCMS(0),
fTPCJet(0),
fTPCFullJet(0),
fEMCALGeometry(0),
fCells(0),
fDoNEF(0),
+ fDoNEFSignalOnly(1),
+ fSignalTrackBias(0),
+ fTrackQA(0),
+ fClusterQA(0),
+ fCalculateRhoJet(0),
+ fDoVertexRCut(0),
+ fMCPartLevel(0),
+ fDoTHnSparse(0),
+ fDoJetRhoDensity(0),
fEMCalPhiMin(1.39626),
fEMCalPhiMax(3.26377),
fEMCalPhiTotal(1.86750),
fParticlePtUp(200.0),
fParticlePtBins(2000),
fJetR(0.4),
+ fJetRAccept(0.4),
+ fFullEDJetR(0.6),
+ fChargedEDJetR(0.8),
fJetRForRho(0.5),
fJetAreaCutFrac(0.6),
fJetAreaThreshold(0.30159),
fnEMCalCells(12288),
- fScaleFactor(1.50),
+ fScaleFactor(1.28),
fNColl(7),
fTrackMinPt(0.15),
fClusterMinPt(0.3),
fRhoFull(0),
fRhoCharged(0),
fnTracks(0),
+ fnEMCalTracks(0),
fnClusters(0),
fnCaloClusters(0),
fnAKTFullJets(0),
fnEMCalCells=12288; // sMods 1-10 have 24x48 cells, sMods 11&12 have 8x48 cells...
- // Histograms
Int_t TCBins=100;
-
- // QA Plots
- // Hybrid Tracks
- fhTrackPt = new TH1D("fhTrackPt","p_{T} distribution of tracks in event",10*fParticlePtBins,fParticlePtLow,fParticlePtUp);
- fhTrackPt->GetXaxis()->SetTitle("p_{T} (GeV/c)");
- fhTrackPt->GetYaxis()->SetTitle("1/N_{Events} dN/dp_{T}");
- fhTrackPt->Sumw2();
-
- fhTrackPhi = new TH1D("fhTrackPhi","#varphi distribution of tracks in event",TCBins,fTPCPhiMin,fTPCPhiMax);
- fhTrackPhi->GetXaxis()->SetTitle("#varphi");
- fhTrackPhi->GetYaxis()->SetTitle("1/N_{Events} dN/d#varphi");
- fhTrackPhi->Sumw2();
-
- fhTrackEta = new TH1D("fhTrackEta","#eta distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax);
- fhTrackEta->GetXaxis()->SetTitle("#eta");
- fhTrackEta->GetYaxis()->SetTitle("1/N_{Events} dN/d#eta");
- fhTrackEta->Sumw2();
-
- fhTrackEtaPhi = new TH2D("fhTrackEtaPhi","#eta-#varphi distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax,TCBins,fTPCPhiMin,fTPCPhiMax);
- fhTrackEtaPhi->GetXaxis()->SetTitle("#eta");
- fhTrackEtaPhi->GetYaxis()->SetTitle("#varphi");
- fhTrackEtaPhi->GetZaxis()->SetTitle("1/N_{Events} dN/d#etad#varphi");
- fhTrackEtaPhi->Sumw2();
-
- fhTrackPhiPt = new TH2D("fhTrackPhiPt","#varphi-p_{T} distribution of tracks in event",TCBins,fTPCPhiMin,fTPCPhiMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
- fhTrackPhiPt->GetXaxis()->SetTitle("#varphi");
- fhTrackPhiPt->GetYaxis()->SetTitle("p_{T} (GeV/c)");
- fhTrackPhiPt->GetZaxis()->SetTitle("1/N_{Events} dN/d#varphidp_{T}");
- fhTrackPhiPt->Sumw2();
-
- fhTrackEtaPt = new TH2D("fhTrackEtaPt","#eta-p_{T} distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
- fhTrackEtaPt->GetXaxis()->SetTitle("#varphi");
- fhTrackEtaPt->GetYaxis()->SetTitle("p_{T} (GeV/c)");
- fhTrackEtaPt->GetZaxis()->SetTitle("1/N_{Events} dN/d#etadp_{T}");
- fhTrackEtaPt->Sumw2();
-
- fhTrackEtaPhiPt = new TH3D("fhTrackEtaPhiPt","#eta-#varphi-p_{T} distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax,TCBins,fTPCPhiMin,fTPCPhiMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
- fhTrackEtaPhiPt->GetXaxis()->SetTitle("#eta");
- fhTrackEtaPhiPt->GetYaxis()->SetTitle("#varphi");
- fhTrackEtaPhiPt->GetZaxis()->SetTitle("p_{T} (GeV/c)");
- fhTrackEtaPhiPt->Sumw2();
-
- // Global Tracks
- fhGlobalTrackPt = new TH1D("fhGlobalTrackPt","Global p_{T} distribution of tracks in event",10*fParticlePtBins,fParticlePtLow,fParticlePtUp);
- fhGlobalTrackPt->GetXaxis()->SetTitle("p_{T} (GeV/c)");
- fhGlobalTrackPt->GetYaxis()->SetTitle("1/N_{Events} dN/dp_{T}");
- fhGlobalTrackPt->Sumw2();
-
- fhGlobalTrackPhi = new TH1D("fhGlobalTrackPhi","Global #varphi distribution of tracks in event",TCBins,fTPCPhiMin,fTPCPhiMax);
- fhGlobalTrackPhi->GetXaxis()->SetTitle("#varphi");
- fhGlobalTrackPhi->GetYaxis()->SetTitle("1/N_{Events} dN/d#varphi");
- fhGlobalTrackPhi->Sumw2();
-
- fhGlobalTrackEta = new TH1D("fhGlobalTrackEta","Global #eta distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax);
- fhGlobalTrackEta->GetXaxis()->SetTitle("#eta");
- fhGlobalTrackEta->GetYaxis()->SetTitle("1/N_{Events} dN/d#eta");
- fhGlobalTrackEta->Sumw2();
-
- fhGlobalTrackEtaPhi = new TH2D("fhGlobalTrackEtaPhi","Global #eta-#varphi distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax,TCBins,fTPCPhiMin,fTPCPhiMax);
- fhGlobalTrackEtaPhi->GetXaxis()->SetTitle("#eta");
- fhGlobalTrackEtaPhi->GetYaxis()->SetTitle("#varphi");
- fhGlobalTrackEtaPhi->GetZaxis()->SetTitle("1/N_{Events} dN/d#etad#varphi");
- fhGlobalTrackEtaPhi->Sumw2();
-
- fhGlobalTrackPhiPt = new TH2D("fhGlobalTrackPhiPt","Global #varphi-p_{T} distribution of tracks in event",TCBins,fTPCPhiMin,fTPCPhiMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
- fhGlobalTrackPhiPt->GetXaxis()->SetTitle("#varphi");
- fhGlobalTrackPhiPt->GetYaxis()->SetTitle("p_{T} (GeV/c)");
- fhGlobalTrackPhiPt->GetZaxis()->SetTitle("1/N_{Events} dN/d#varphidp_{T}");
- fhGlobalTrackPhiPt->Sumw2();
-
- fhGlobalTrackEtaPt = new TH2D("fhGlobalTrackEtaPt","Global #eta-p_{T} distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
- fhGlobalTrackEtaPt->GetXaxis()->SetTitle("#varphi");
- fhGlobalTrackEtaPt->GetYaxis()->SetTitle("p_{T} (GeV/c)");
- fhGlobalTrackEtaPt->GetZaxis()->SetTitle("1/N_{Events} dN/d#etadp_{T}");
- fhGlobalTrackEtaPt->Sumw2();
-
- fhGlobalTrackEtaPhiPt = new TH3D("fhGlobalTrackEtaPhiPt","Global #eta-#varphi-p_{T} distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax,TCBins,fTPCPhiMin,fTPCPhiMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
- fhGlobalTrackEtaPhiPt->GetXaxis()->SetTitle("#eta");
- fhGlobalTrackEtaPhiPt->GetYaxis()->SetTitle("#varphi");
- fhGlobalTrackEtaPhiPt->GetZaxis()->SetTitle("p_{T} (GeV/c)");
- fhGlobalTrackEtaPhiPt->Sumw2();
-
- // Complementary Tracks
- fhComplementaryTrackPt = new TH1D("fhComplementaryTrackPt","Complementary p_{T} distribution of tracks in event",10*fParticlePtBins,fParticlePtLow,fParticlePtUp);
- fhComplementaryTrackPt->GetXaxis()->SetTitle("p_{T} (GeV/c)");
- fhComplementaryTrackPt->GetYaxis()->SetTitle("1/N_{Events} dN/dp_{T}");
- fhComplementaryTrackPt->Sumw2();
-
- fhComplementaryTrackPhi = new TH1D("fhComplementaryTrackPhi","Complementary #varphi distribution of tracks in event",TCBins,fTPCPhiMin,fTPCPhiMax);
- fhComplementaryTrackPhi->GetXaxis()->SetTitle("#varphi");
- fhComplementaryTrackPhi->GetYaxis()->SetTitle("1/N_{Events} dN/d#varphi");
- fhComplementaryTrackPhi->Sumw2();
-
- fhComplementaryTrackEta = new TH1D("fhComplementaryTrackEta","Complementary #eta distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax);
- fhComplementaryTrackEta->GetXaxis()->SetTitle("#eta");
- fhComplementaryTrackEta->GetYaxis()->SetTitle("1/N_{Events} dN/d#eta");
- fhComplementaryTrackEta->Sumw2();
-
- fhComplementaryTrackEtaPhi = new TH2D("fhComplementaryTrackEtaPhi","Complementary #eta-#varphi distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax,TCBins,fTPCPhiMin,fTPCPhiMax);
- fhComplementaryTrackEtaPhi->GetXaxis()->SetTitle("#eta");
- fhComplementaryTrackEtaPhi->GetYaxis()->SetTitle("#varphi");
- fhComplementaryTrackEtaPhi->GetZaxis()->SetTitle("1/N_{Events} dN/d#etad#varphi");
- fhComplementaryTrackEtaPhi->Sumw2();
-
- fhComplementaryTrackPhiPt = new TH2D("fhComplementaryTrackPhiPt","Complementary #varphi-p_{T} distribution of tracks in event",TCBins,fTPCPhiMin,fTPCPhiMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
- fhComplementaryTrackPhiPt->GetXaxis()->SetTitle("#varphi");
- fhComplementaryTrackPhiPt->GetYaxis()->SetTitle("p_{T} (GeV/c)");
- fhComplementaryTrackPhiPt->GetZaxis()->SetTitle("1/N_{Events} dN/d#varphidp_{T}");
- fhComplementaryTrackPhiPt->Sumw2();
-
- fhComplementaryTrackEtaPt = new TH2D("fhComplementaryTrackEtaPt","Complementary #eta-p_{T} distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
- fhComplementaryTrackEtaPt->GetXaxis()->SetTitle("#varphi");
- fhComplementaryTrackEtaPt->GetYaxis()->SetTitle("p_{T} (GeV/c)");
- fhComplementaryTrackEtaPt->GetZaxis()->SetTitle("1/N_{Events} dN/d#etadp_{T}");
- fhComplementaryTrackEtaPt->Sumw2();
-
- fhComplementaryTrackEtaPhiPt = new TH3D("fhComplementaryTrackEtaPhiPt","Complementary #eta-#varphi-p_{T} distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax,TCBins,fTPCPhiMin,fTPCPhiMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
- fhComplementaryTrackEtaPhiPt->GetXaxis()->SetTitle("#eta");
- fhComplementaryTrackEtaPhiPt->GetYaxis()->SetTitle("#varphi");
- fhComplementaryTrackEtaPhiPt->GetZaxis()->SetTitle("p_{T} (GeV/c)");
- fhComplementaryTrackEtaPhiPt->Sumw2();
-
- // Corrected Calo Clusters
- fhClusterPt = new TH1D("fhClusterPt","p_{T} distribution of clusters in event",10*fParticlePtBins,fParticlePtLow,fParticlePtUp);
- fhClusterPt->GetXaxis()->SetTitle("p_{T} (GeV/c)");
- fhClusterPt->GetYaxis()->SetTitle("1/N_{Events} dN/dp_{T}");
- fhClusterPt->Sumw2();
-
- fhClusterPhi = new TH1D("fhClusterPhi","#varphi distribution of clusters in event",TCBins,fTPCPhiMin,fTPCPhiMax);
- fhClusterPhi->GetXaxis()->SetTitle("#varphi");
- fhClusterPhi->GetYaxis()->SetTitle("1/N_{Events} dN/d#varphi");
- fhClusterPhi->Sumw2();
-
- fhClusterEta = new TH1D("fhClusterEta","#eta distribution of clusters in event",TCBins,fTPCEtaMin,fTPCEtaMax);
- fhClusterEta->GetXaxis()->SetTitle("#eta");
- fhClusterEta->GetYaxis()->SetTitle("1/N_{Events} dN/d#eta");
- fhClusterEta->Sumw2();
-
- fhClusterEtaPhi = new TH2D("fhClusterEtaPhi","#eta-#varphi distribution of clusters in event",TCBins,fEMCalEtaMin,fEMCalEtaMax,TCBins,fEMCalPhiMin,fEMCalPhiMax);
- fhClusterEtaPhi->GetXaxis()->SetTitle("#eta");
- fhClusterEtaPhi->GetYaxis()->SetTitle("#varphi");
- fhClusterEtaPhi->GetZaxis()->SetTitle("1/N_{Events} dN/d#etad#varphi");
- fhClusterEtaPhi->Sumw2();
-
- fhClusterPhiPt = new TH2D("fhClusterPhiPt","#varphi-p_{T} distribution of clusters in event",TCBins,fEMCalPhiMin,fEMCalPhiMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
- fhClusterPhiPt->GetXaxis()->SetTitle("#varphi");
- fhClusterPhiPt->GetYaxis()->SetTitle("p_{T} (GeV/c)");
- fhClusterPhiPt->GetZaxis()->SetTitle("1/N_{Events} dN/d#varphidp_{T}");
- fhClusterPhiPt->Sumw2();
-
- fhClusterEtaPt = new TH2D("fhClusterEtaPt","#eta-p_{T} distribution of clusters in event",TCBins,fEMCalEtaMin,fEMCalEtaMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
- fhClusterEtaPt->GetXaxis()->SetTitle("#varphi");
- fhClusterEtaPt->GetYaxis()->SetTitle("p_{T} (GeV/c)");
- fhClusterEtaPt->GetZaxis()->SetTitle("1/N_{Events} dN/d#etadp_{T}");
- fhClusterEtaPt->Sumw2();
-
- fhClusterEtaPhiPt = new TH3D("fhClusterEtaPhiPt","#eta-#varphi-p_{T} distribution of clusters in event",TCBins,fEMCalEtaMin,fEMCalEtaMax,TCBins,fEMCalPhiMin,fEMCalPhiMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
- fhClusterEtaPhiPt->GetXaxis()->SetTitle("#eta");
- fhClusterEtaPhiPt->GetYaxis()->SetTitle("#varphi");
- fhClusterEtaPhiPt->GetZaxis()->SetTitle("p_{T} (GeV/c)");
- fhClusterEtaPhiPt->Sumw2();
-
- fhCentrality = new TH1D("fhCentrality","Event Centrality Distribution",fCentralityBins*CentralityBinMult,fCentralityLow,fCentralityUp);
- fhCentrality->GetXaxis()->SetTitle(fCentralityTag);
+ Int_t multBins=200;
+ Double_t multLow=0;
+ Double_t multUp=200;
+
+ // Track QA Plots
+ if (fTrackQA==kTRUE)
+ {
+ flTrack = new TList();
+ flTrack->SetName("TrackQA");
+
+ // Hybrid Tracks
+ fhTrackPt = new TH1F("fhTrackPt","p_{T} distribution of tracks in event",10*fParticlePtBins,fParticlePtLow,fParticlePtUp);
+ fhTrackPt->GetXaxis()->SetTitle("p_{T} (GeV/c)");
+ fhTrackPt->GetYaxis()->SetTitle("1/N_{Events} dN/dp_{T}");
+ fhTrackPt->Sumw2();
+
+ fhTrackPhi = new TH1F("fhTrackPhi","#varphi distribution of tracks in event",TCBins,fTPCPhiMin,fTPCPhiMax);
+ fhTrackPhi->GetXaxis()->SetTitle("#varphi");
+ fhTrackPhi->GetYaxis()->SetTitle("1/N_{Events} dN/d#varphi");
+ fhTrackPhi->Sumw2();
+
+ fhTrackEta = new TH1F("fhTrackEta","#eta distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax);
+ fhTrackEta->GetXaxis()->SetTitle("#eta");
+ fhTrackEta->GetYaxis()->SetTitle("1/N_{Events} dN/d#eta");
+ fhTrackEta->Sumw2();
+
+ fhTrackEtaPhi = new TH2F("fhTrackEtaPhi","#eta-#varphi distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax,TCBins,fTPCPhiMin,fTPCPhiMax);
+ fhTrackEtaPhi->GetXaxis()->SetTitle("#eta");
+ fhTrackEtaPhi->GetYaxis()->SetTitle("#varphi");
+ fhTrackEtaPhi->GetZaxis()->SetTitle("1/N_{Events} dN/d#etad#varphi");
+ fhTrackEtaPhi->Sumw2();
+
+ fhTrackPhiPt = new TH2F("fhTrackPhiPt","#varphi-p_{T} distribution of tracks in event",TCBins,fTPCPhiMin,fTPCPhiMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
+ fhTrackPhiPt->GetXaxis()->SetTitle("#varphi");
+ fhTrackPhiPt->GetYaxis()->SetTitle("p_{T} (GeV/c)");
+ fhTrackPhiPt->GetZaxis()->SetTitle("1/N_{Events} dN/d#varphidp_{T}");
+ fhTrackPhiPt->Sumw2();
+
+ fhTrackEtaPt = new TH2F("fhTrackEtaPt","#eta-p_{T} distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
+ fhTrackEtaPt->GetXaxis()->SetTitle("#varphi");
+ fhTrackEtaPt->GetYaxis()->SetTitle("p_{T} (GeV/c)");
+ fhTrackEtaPt->GetZaxis()->SetTitle("1/N_{Events} dN/d#etadp_{T}");
+ fhTrackEtaPt->Sumw2();
+
+ // Global Tracks
+ fhGlobalTrackPt = new TH1F("fhGlobalTrackPt","Global p_{T} distribution of tracks in event",10*fParticlePtBins,fParticlePtLow,fParticlePtUp);
+ fhGlobalTrackPt->GetXaxis()->SetTitle("p_{T} (GeV/c)");
+ fhGlobalTrackPt->GetYaxis()->SetTitle("1/N_{Events} dN/dp_{T}");
+ fhGlobalTrackPt->Sumw2();
+
+ fhGlobalTrackPhi = new TH1F("fhGlobalTrackPhi","Global #varphi distribution of tracks in event",TCBins,fTPCPhiMin,fTPCPhiMax);
+ fhGlobalTrackPhi->GetXaxis()->SetTitle("#varphi");
+ fhGlobalTrackPhi->GetYaxis()->SetTitle("1/N_{Events} dN/d#varphi");
+ fhGlobalTrackPhi->Sumw2();
+
+ fhGlobalTrackEta = new TH1F("fhGlobalTrackEta","Global #eta distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax);
+ fhGlobalTrackEta->GetXaxis()->SetTitle("#eta");
+ fhGlobalTrackEta->GetYaxis()->SetTitle("1/N_{Events} dN/d#eta");
+ fhGlobalTrackEta->Sumw2();
+
+ fhGlobalTrackEtaPhi = new TH2F("fhGlobalTrackEtaPhi","Global #eta-#varphi distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax,TCBins,fTPCPhiMin,fTPCPhiMax);
+ fhGlobalTrackEtaPhi->GetXaxis()->SetTitle("#eta");
+ fhGlobalTrackEtaPhi->GetYaxis()->SetTitle("#varphi");
+ fhGlobalTrackEtaPhi->GetZaxis()->SetTitle("1/N_{Events} dN/d#etad#varphi");
+ fhGlobalTrackEtaPhi->Sumw2();
+
+ fhGlobalTrackPhiPt = new TH2F("fhGlobalTrackPhiPt","Global #varphi-p_{T} distribution of tracks in event",TCBins,fTPCPhiMin,fTPCPhiMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
+ fhGlobalTrackPhiPt->GetXaxis()->SetTitle("#varphi");
+ fhGlobalTrackPhiPt->GetYaxis()->SetTitle("p_{T} (GeV/c)");
+ fhGlobalTrackPhiPt->GetZaxis()->SetTitle("1/N_{Events} dN/d#varphidp_{T}");
+ fhGlobalTrackPhiPt->Sumw2();
+
+ fhGlobalTrackEtaPt = new TH2F("fhGlobalTrackEtaPt","Global #eta-p_{T} distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
+ fhGlobalTrackEtaPt->GetXaxis()->SetTitle("#varphi");
+ fhGlobalTrackEtaPt->GetYaxis()->SetTitle("p_{T} (GeV/c)");
+ fhGlobalTrackEtaPt->GetZaxis()->SetTitle("1/N_{Events} dN/d#etadp_{T}");
+ fhGlobalTrackEtaPt->Sumw2();
+
+ // Complementary Tracks
+ fhComplementaryTrackPt = new TH1F("fhComplementaryTrackPt","Complementary p_{T} distribution of tracks in event",10*fParticlePtBins,fParticlePtLow,fParticlePtUp);
+ fhComplementaryTrackPt->GetXaxis()->SetTitle("p_{T} (GeV/c)");
+ fhComplementaryTrackPt->GetYaxis()->SetTitle("1/N_{Events} dN/dp_{T}");
+ fhComplementaryTrackPt->Sumw2();
+
+ fhComplementaryTrackPhi = new TH1F("fhComplementaryTrackPhi","Complementary #varphi distribution of tracks in event",TCBins,fTPCPhiMin,fTPCPhiMax);
+ fhComplementaryTrackPhi->GetXaxis()->SetTitle("#varphi");
+ fhComplementaryTrackPhi->GetYaxis()->SetTitle("1/N_{Events} dN/d#varphi");
+ fhComplementaryTrackPhi->Sumw2();
+
+ fhComplementaryTrackEta = new TH1F("fhComplementaryTrackEta","Complementary #eta distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax);
+ fhComplementaryTrackEta->GetXaxis()->SetTitle("#eta");
+ fhComplementaryTrackEta->GetYaxis()->SetTitle("1/N_{Events} dN/d#eta");
+ fhComplementaryTrackEta->Sumw2();
+
+ fhComplementaryTrackEtaPhi = new TH2F("fhComplementaryTrackEtaPhi","Complementary #eta-#varphi distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax,TCBins,fTPCPhiMin,fTPCPhiMax);
+ fhComplementaryTrackEtaPhi->GetXaxis()->SetTitle("#eta");
+ fhComplementaryTrackEtaPhi->GetYaxis()->SetTitle("#varphi");
+ fhComplementaryTrackEtaPhi->GetZaxis()->SetTitle("1/N_{Events} dN/d#etad#varphi");
+ fhComplementaryTrackEtaPhi->Sumw2();
+
+ fhComplementaryTrackPhiPt = new TH2F("fhComplementaryTrackPhiPt","Complementary #varphi-p_{T} distribution of tracks in event",TCBins,fTPCPhiMin,fTPCPhiMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
+ fhComplementaryTrackPhiPt->GetXaxis()->SetTitle("#varphi");
+ fhComplementaryTrackPhiPt->GetYaxis()->SetTitle("p_{T} (GeV/c)");
+ fhComplementaryTrackPhiPt->GetZaxis()->SetTitle("1/N_{Events} dN/d#varphidp_{T}");
+ fhComplementaryTrackPhiPt->Sumw2();
+
+ fhComplementaryTrackEtaPt = new TH2F("fhComplementaryTrackEtaPt","Complementary #eta-p_{T} distribution of tracks in event",TCBins,fTPCEtaMin,fTPCEtaMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
+ fhComplementaryTrackEtaPt->GetXaxis()->SetTitle("#varphi");
+ fhComplementaryTrackEtaPt->GetYaxis()->SetTitle("p_{T} (GeV/c)");
+ fhComplementaryTrackEtaPt->GetZaxis()->SetTitle("1/N_{Events} dN/d#etadp_{T}");
+ fhComplementaryTrackEtaPt->Sumw2();
+
+ fhTPCEventMult = new TH2F("fhTPCEventMult","TPC Event Multiplcity vs Centrality",CentralityBinMult*fCentralityBins,fCentralityLow,fCentralityUp,multBins,multLow,multUp);
+ fhTPCEventMult->GetXaxis()->SetTitle(fCentralityTag.Data());
+ fhTPCEventMult->GetYaxis()->SetTitle("Multiplicity");
+ fhTPCEventMult->GetZaxis()->SetTitle("1/N_{Events} dN/dCentdN_{Charged}");
+ fhTPCEventMult->Sumw2();
+
+ fhEMCalTrackEventMult = new TH2F("fhEMCalTrackEventMult","EMCal Track Event Multiplcity vs Centrality",CentralityBinMult*fCentralityBins,fCentralityLow,fCentralityUp,multBins,multLow,multUp);
+ fhEMCalTrackEventMult->GetXaxis()->SetTitle(fCentralityTag.Data());
+ fhEMCalTrackEventMult->GetYaxis()->SetTitle("Multiplicity");
+ fhEMCalTrackEventMult->GetZaxis()->SetTitle("1/N_{Events} dN/dCentdN_{Neutral}");
+ fhEMCalTrackEventMult->Sumw2();
+
+ fpTPCEventMult = new TProfile("fpTPCEventMult","TPC Event Multiplcity vs Centrality",CentralityBinMult*fCentralityBins,fCentralityLow,fCentralityUp);
+ fpTPCEventMult->GetXaxis()->SetTitle(fCentralityTag.Data());
+ fpTPCEventMult->GetYaxis()->SetTitle("Multiplicity");
+
+ // QA::2D Energy Density Profiles for Tracks and Clusters
+ fpTrackPtProfile = new TProfile2D("fpTrackPtProfile","2D Profile of track pT density throughout the TPC",TCBins,fTPCEtaMin,fTPCEtaMax,TCBins,fTPCPhiMin,fTPCPhiMax);
+ fpTrackPtProfile->GetXaxis()->SetTitle("#eta");
+ fpTrackPtProfile->GetYaxis()->SetTitle("#varphi");
+ fpTrackPtProfile->GetZaxis()->SetTitle("p_{T} density (GeV/Area)");
+
+ flTrack->Add(fhTrackPt);
+ flTrack->Add(fhTrackEta);
+ flTrack->Add(fhTrackPhi);
+ flTrack->Add(fhTrackEtaPhi);
+ flTrack->Add(fhTrackPhiPt);
+ flTrack->Add(fhTrackEtaPt);
+ flTrack->Add(fhGlobalTrackPt);
+ flTrack->Add(fhGlobalTrackEta);
+ flTrack->Add(fhGlobalTrackPhi);
+ flTrack->Add(fhGlobalTrackEtaPhi);
+ flTrack->Add(fhGlobalTrackPhiPt);
+ flTrack->Add(fhGlobalTrackEtaPt);
+ flTrack->Add(fhComplementaryTrackPt);
+ flTrack->Add(fhComplementaryTrackEta);
+ flTrack->Add(fhComplementaryTrackPhi);
+ flTrack->Add(fhComplementaryTrackEtaPhi);
+ flTrack->Add(fhComplementaryTrackPhiPt);
+ flTrack->Add(fhComplementaryTrackEtaPt);
+ flTrack->Add(fhTPCEventMult);
+ flTrack->Add(fhEMCalTrackEventMult);
+ flTrack->Add(fpTPCEventMult);
+ flTrack->Add(fpTrackPtProfile);
+ fOutput->Add(flTrack);
+ }
+
+ if (fClusterQA==kTRUE)
+ {
+ flCluster = new TList();
+ flCluster->SetName("ClusterQA");
+
+ fhClusterPt = new TH1F("fhClusterPt","p_{T} distribution of clusters in event",10*fParticlePtBins,fParticlePtLow,fParticlePtUp);
+ fhClusterPt->GetXaxis()->SetTitle("p_{T} (GeV/c)");
+ fhClusterPt->GetYaxis()->SetTitle("1/N_{Events} dN/dp_{T}");
+ fhClusterPt->Sumw2();
+
+ fhClusterPhi = new TH1F("fhClusterPhi","#varphi distribution of clusters in event",TCBins,fTPCPhiMin,fTPCPhiMax);
+ fhClusterPhi->GetXaxis()->SetTitle("#varphi");
+ fhClusterPhi->GetYaxis()->SetTitle("1/N_{Events} dN/d#varphi");
+ fhClusterPhi->Sumw2();
+
+ fhClusterEta = new TH1F("fhClusterEta","#eta distribution of clusters in event",TCBins,fTPCEtaMin,fTPCEtaMax);
+ fhClusterEta->GetXaxis()->SetTitle("#eta");
+ fhClusterEta->GetYaxis()->SetTitle("1/N_{Events} dN/d#eta");
+ fhClusterEta->Sumw2();
+
+ fhClusterEtaPhi = new TH2F("fhClusterEtaPhi","#eta-#varphi distribution of clusters in event",TCBins,fEMCalEtaMin,fEMCalEtaMax,TCBins,fEMCalPhiMin,fEMCalPhiMax);
+ fhClusterEtaPhi->GetXaxis()->SetTitle("#eta");
+ fhClusterEtaPhi->GetYaxis()->SetTitle("#varphi");
+ fhClusterEtaPhi->GetZaxis()->SetTitle("1/N_{Events} dN/d#etad#varphi");
+ fhClusterEtaPhi->Sumw2();
+
+ fhClusterPhiPt = new TH2F("fhClusterPhiPt","#varphi-p_{T} distribution of clusters in event",TCBins,fEMCalPhiMin,fEMCalPhiMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
+ fhClusterPhiPt->GetXaxis()->SetTitle("#varphi");
+ fhClusterPhiPt->GetYaxis()->SetTitle("p_{T} (GeV/c)");
+ fhClusterPhiPt->GetZaxis()->SetTitle("1/N_{Events} dN/d#varphidp_{T}");
+ fhClusterPhiPt->Sumw2();
+
+ fhClusterEtaPt = new TH2F("fhClusterEtaPt","#eta-p_{T} distribution of clusters in event",TCBins,fEMCalEtaMin,fEMCalEtaMax,fParticlePtBins,fParticlePtLow,fParticlePtUp);
+ fhClusterEtaPt->GetXaxis()->SetTitle("#varphi");
+ fhClusterEtaPt->GetYaxis()->SetTitle("p_{T} (GeV/c)");
+ fhClusterEtaPt->GetZaxis()->SetTitle("1/N_{Events} dN/d#etadp_{T}");
+ fhClusterEtaPt->Sumw2();
+
+ fhEMCalEventMult = new TH2F("fhEMCalEventMult","EMCal Event Multiplcity vs Centrality",CentralityBinMult*fCentralityBins,fCentralityLow,fCentralityUp,multBins,multLow,multUp);
+ fhEMCalEventMult->GetXaxis()->SetTitle(fCentralityTag.Data());
+ fhEMCalEventMult->GetYaxis()->SetTitle("Multiplicity");
+ fhEMCalEventMult->GetZaxis()->SetTitle("1/N_{Events} dN/dCentdN_{Neutral}");
+ fhEMCalEventMult->Sumw2();
+
+ fpEMCalEventMult = new TProfile("fpEMCalEventMult","EMCal Event Multiplcity vs Centrality",CentralityBinMult*fCentralityBins,fCentralityLow,fCentralityUp);
+ fpEMCalEventMult->GetXaxis()->SetTitle(fCentralityTag.Data());
+ fpEMCalEventMult->GetYaxis()->SetTitle("Multiplicity");
+
+ fpClusterPtProfile = new TProfile2D("fpClusterPtProfile","2D Profile of cluster pT density throughout the EMCal",TCBins,fEMCalEtaMin,fEMCalEtaMax,TCBins,fEMCalPhiMin,fEMCalPhiMax);
+ fpClusterPtProfile->GetXaxis()->SetTitle("#eta");
+ fpClusterPtProfile->GetYaxis()->SetTitle("#varphi");
+ fpClusterPtProfile->GetZaxis()->SetTitle("p_{T} density (GeV/Area)");
+
+ fhEMCalCellCounts = new TH1F("fhEMCalCellCounts","Distribtuion of cluster counts across the EMCal",fnEMCalCells,1,fnEMCalCells);
+ fhEMCalCellCounts->GetXaxis()->SetTitle("Absoulute Cell Id");
+ fhEMCalCellCounts->GetYaxis()->SetTitle("Counts per Event");
+ fhEMCalCellCounts->Sumw2();
+
+ flCluster->Add(fhClusterPt);
+ flCluster->Add(fhClusterEta);
+ flCluster->Add(fhClusterPhi);
+ flCluster->Add(fhClusterEtaPhi);
+ flCluster->Add(fhClusterPhiPt);
+ flCluster->Add(fhClusterEtaPt);
+ flCluster->Add(fhEMCalEventMult);
+ flCluster->Add(fpEMCalEventMult);
+ flCluster->Add(fpClusterPtProfile);
+ flCluster->Add(fhEMCalCellCounts);
+ fOutput->Add(flCluster);
+ }
+
+ if (fCalculateRhoJet>=0 && fMCPartLevel==kFALSE) // Default Rho & Raw Jet Spectra
+ {
+ fEMCalRawJets = new AlipAJetHistos("EMCalRawJets",fCentralityTag.Data());
+
+ fRhoChargedCMSScale = new AlipAJetHistos("RhoChargedCMSScale",fCentralityTag.Data(),fDoNEF,fDoNEFSignalOnly,fDoTHnSparse);
+ fRhoChargedCMSScale->SetSignalTrackPtBias(fSignalTrackBias);
+ fOutput->Add(fEMCalRawJets->GetOutputHistos());
+ fOutput->Add(fRhoChargedCMSScale->GetOutputHistos());
+
+ }
+ if (fCalculateRhoJet>=1) // Basic Rho & Raw Jet Spectra
+ {
+ fRhoChargedScale = new AlipAJetHistos("RhoChargedScale",fCentralityTag.Data());
+ fRhoChargedScale->SetSignalTrackPtBias(fSignalTrackBias);
+
+ if (fMCPartLevel==kTRUE)
+ {
+ fTPCRawJets = new AlipAJetHistos("TPCRawJets",fCentralityTag.Data());
+ fRhoChargedCMS = new AlipAJetHistos("RhoChargedCMS",fCentralityTag.Data());
+ fRhoChargedCMS->SetSignalTrackPtBias(fSignalTrackBias);
+
+ fOutput->Add(fTPCRawJets->GetOutputHistos());
+ fOutput->Add(fRhoChargedCMS->GetOutputHistos());
+ }
+
+ if (fDoJetRhoDensity == kTRUE)
+ {
+ Double_t fullEDR = fFullEDJetR *100;
+ Double_t chargedEDR = fChargedEDJetR *100;
+ const Int_t fullEDBins = (Int_t) fullEDR;
+ const Int_t chargedEDBins = (Int_t) chargedEDR;
+
+ fpFullJetEDProfile = new TProfile3D("fpFullJetEDProfile","Jet ED Profile for #varphi_{jet}#in(#varphi_{EMCal min} + R,#varphi_{EMCal max} - R) & #eta_{jet}#in(#eta_{EMCal min} + R,#eta_{EMCal max} - R)",fParticlePtBins,fParticlePtLow,fParticlePtUp,fCentralityBins,fCentralityLow,fCentralityUp,fullEDBins,0,fFullEDJetR);
+ fpFullJetEDProfile->GetXaxis()->SetTitle("p_{T,jet}^{ch+em} (GeV)");
+ fpFullJetEDProfile->GetYaxis()->SetTitle(fCentralityTag.Data());
+ fpFullJetEDProfile->GetZaxis()->SetTitle("R");
+
+ fpChargedJetEDProfile = new TProfile3D("fpChargedJetEDProfile","Charged Jet ED Profile for #eta_{jet}#in(#eta_{TPC min} + R,#eta_{TPC max} - R)",fParticlePtBins,fParticlePtLow,fParticlePtUp,fCentralityBins,fCentralityLow,fCentralityUp,chargedEDBins,0,fChargedEDJetR);
+ fpChargedJetEDProfile->GetXaxis()->SetTitle("p_{T,jet}^{ch} (GeV)");
+ fpChargedJetEDProfile->GetYaxis()->SetTitle(fCentralityTag.Data());
+ fpChargedJetEDProfile->GetZaxis()->SetTitle("R");
+
+ fpChargedJetEDProfileScaled = new TProfile3D("fpChargedJetEDProfileScaled","Charged Jet ED Profile x SF for #eta_{jet}#in(#eta_{TPC min} + R,#eta_{TPC max} - R)",fParticlePtBins,fParticlePtLow,fParticlePtUp,fCentralityBins,fCentralityLow,fCentralityUp,chargedEDBins,0,fChargedEDJetR);
+ fpChargedJetEDProfileScaled->GetXaxis()->SetTitle("p_{T,jet}^{ch} (GeV)");
+ fpChargedJetEDProfileScaled->GetYaxis()->SetTitle(fCentralityTag.Data());
+ fpChargedJetEDProfileScaled->GetZaxis()->SetTitle("R");
+
+ fOutput->Add(fpFullJetEDProfile);
+ fOutput->Add(fpChargedJetEDProfile);
+ fOutput->Add(fpChargedJetEDProfileScaled);
+ }
+ fOutput->Add(fRhoChargedScale->GetOutputHistos());
+ }
+ if (fCalculateRhoJet>=2 && fMCPartLevel==kFALSE) // Charged Rho & Raw Jet Spectra
+ {
+ fTPCRawJets = new AlipAJetHistos("TPCRawJets",fCentralityTag.Data());
+ fRhoChargedCMS = new AlipAJetHistos("RhoChargedCMS",fCentralityTag.Data());
+ fRhoChargedCMS->SetSignalTrackPtBias(fSignalTrackBias);
+
+ fOutput->Add(fTPCRawJets->GetOutputHistos());
+ fOutput->Add(fRhoChargedCMS->GetOutputHistos());
+ }
+ if (fCalculateRhoJet>=3 && fMCPartLevel==kFALSE) // Basic Rho & Raw Jet Spectra
+ {
+ fRhoFull0 = new AlipAJetHistos("RhoFull0",fCentralityTag.Data());
+ fRhoFull0->SetSignalTrackPtBias(fSignalTrackBias);
+ fRhoFull1 = new AlipAJetHistos("RhoFull1",fCentralityTag.Data());
+ fRhoFull1->SetSignalTrackPtBias(fSignalTrackBias);
+ fRhoFull2 = new AlipAJetHistos("RhoFull2",fCentralityTag.Data());
+ fRhoFull2->SetSignalTrackPtBias(fSignalTrackBias);
+ fRhoFullN = new AlipAJetHistos("RhoFullN",fCentralityTag.Data());
+ fRhoFullN->SetSignalTrackPtBias(fSignalTrackBias);
+ fRhoFullDijet = new AlipAJetHistos("RhoFullDijet",fCentralityTag.Data());
+ fRhoFullDijet->SetSignalTrackPtBias(fSignalTrackBias);
+ fRhoFullkT = new AlipAJetHistos("RhoFullkT",fCentralityTag.Data());
+ fRhoFullkT->SetSignalTrackPtBias(fSignalTrackBias);
+ fRhoFullCMS = new AlipAJetHistos("RhoFullCMS",fCentralityTag.Data());
+ fRhoFullCMS->SetSignalTrackPtBias(fSignalTrackBias);
+ fRhoCharged0 = new AlipAJetHistos("RhoCharged0",fCentralityTag.Data());
+ fRhoCharged0->SetSignalTrackPtBias(fSignalTrackBias);
+ fRhoCharged1 = new AlipAJetHistos("RhoCharged1",fCentralityTag.Data());
+ fRhoCharged1->SetSignalTrackPtBias(fSignalTrackBias);
+ fRhoCharged2 = new AlipAJetHistos("RhoCharged2",fCentralityTag.Data());
+ fRhoCharged2->SetSignalTrackPtBias(fSignalTrackBias);
+ fRhoChargedN = new AlipAJetHistos("RhoChargedN",fCentralityTag.Data());
+ fRhoChargedN->SetSignalTrackPtBias(fSignalTrackBias);
+ fRhoChargedkT = new AlipAJetHistos("RhoChargedkT",fCentralityTag.Data());
+ fRhoChargedkT->SetSignalTrackPtBias(fSignalTrackBias);
+ fRhoChargedkTScale = new AlipAJetHistos("RhoChargedkTScale",fCentralityTag.Data());
+ fRhoChargedkTScale->SetSignalTrackPtBias(fSignalTrackBias);
+
+ fOutput->Add(fRhoFull0->GetOutputHistos());
+ fOutput->Add(fRhoFull1->GetOutputHistos());
+ fOutput->Add(fRhoFull2->GetOutputHistos());
+ fOutput->Add(fRhoFullN->GetOutputHistos());
+ fOutput->Add(fRhoFullDijet->GetOutputHistos());
+ fOutput->Add(fRhoFullkT->GetOutputHistos());
+ fOutput->Add(fRhoFullCMS->GetOutputHistos());
+ fOutput->Add(fRhoCharged0->GetOutputHistos());
+ fOutput->Add(fRhoCharged1->GetOutputHistos());
+ fOutput->Add(fRhoCharged2->GetOutputHistos());
+ fOutput->Add(fRhoChargedN->GetOutputHistos());
+ fOutput->Add(fRhoChargedkT->GetOutputHistos());
+ fOutput->Add(fRhoChargedkTScale->GetOutputHistos());
+ }
+
+ fhCentrality = new TH1F("fhCentrality","Event Centrality Distribution",fCentralityBins*CentralityBinMult,fCentralityLow,fCentralityUp);
+ fhCentrality->GetXaxis()->SetTitle(fCentralityTag.Data());
fhCentrality->GetYaxis()->SetTitle("1/N_{Events}");
fhCentrality->Sumw2();
- fhEMCalCellCounts = new TH1D("fhEMCalCellCounts","Distribtuion of cluster counts across the EMCal",fnEMCalCells,1,fnEMCalCells);
- fhEMCalCellCounts->GetXaxis()->SetTitle("Absoulute Cell Id");
- fhEMCalCellCounts->GetYaxis()->SetTitle("Counts per Event");
- fhEMCalCellCounts->Sumw2();
-
- Int_t SFBins =100;
- Double_t SFLow=0.0;
- Double_t SFUp=10.0;
-
- fhRhoScale = new TH2D("fhRhoScale","Scaling Factor",SFBins,SFLow,SFUp,CentralityBinMult*fCentralityBins,fCentralityLow,fCentralityUp);
- fhRhoScale->GetXaxis()->SetTitle("Scale Factor");
- fhRhoScale->GetYaxis()->SetTitle("Centrality");
- fhRhoScale->GetZaxis()->SetTitle("Counts");
- fhRhoScale->Sumw2();
-
- // Profiles
- fpEMCalEventMult = new TProfile("fpEMCalEventMult","EMCal Event Multiplcity vs Centrality",CentralityBinMult*fCentralityBins,fCentralityLow,fCentralityUp);
- fpEMCalEventMult->GetXaxis()->SetTitle(fCentralityTag);
- fpEMCalEventMult->GetYaxis()->SetTitle("Multiplicity");
-
- fpTPCEventMult = new TProfile("fpTPCEventMult","TPC Event Multiplcity vs Centrality",CentralityBinMult*fCentralityBins,fCentralityLow,fCentralityUp);
- fpTPCEventMult->GetXaxis()->SetTitle(fCentralityTag);
- fpTPCEventMult->GetYaxis()->SetTitle("Multiplicity");
-
- fpRhoScale = new TProfile("fpRhoScale","Scaling Factor Profile vs Centrality",CentralityBinMult*fCentralityBins,fCentralityLow,fCentralityUp);
- fpRhoScale->GetXaxis()->SetTitle(fCentralityTag);
- fpRhoScale->GetYaxis()->SetTitle("Scale Factor");
-
- // QA::2D Energy Density Profiles for Tracks and Clusters
- fpTrackPtProfile = new TProfile2D("fpTrackPtProfile","2D Profile of track pT density throughout the TPC",TCBins,fTPCEtaMin,fTPCEtaMax,TCBins,fTPCPhiMin,fTPCPhiMax);
- fpTrackPtProfile->GetXaxis()->SetTitle("#eta");
- fpTrackPtProfile->GetYaxis()->SetTitle("#varphi");
- fpTrackPtProfile->GetZaxis()->SetTitle("p_{T} density (GeV/Area)");
-
- fpClusterPtProfile = new TProfile2D("fpClusterPtProfile","2D Profile of cluster pT density throughout the EMCal",TCBins,fEMCalEtaMin,fEMCalEtaMax,TCBins,fEMCalPhiMin,fEMCalPhiMax);
- fpClusterPtProfile->GetXaxis()->SetTitle("#eta");
- fpClusterPtProfile->GetYaxis()->SetTitle("#varphi");
- fpClusterPtProfile->GetZaxis()->SetTitle("p_{T} density (GeV/Area)");
-
- fTPCRawJets = new AlipAJetHistos("TPCRawJets",fCentralityTag);
- fEMCalRawJets = new AlipAJetHistos("EMCalRawJets",fCentralityTag);
-
- /*
- fRhoFull0 = new AlipAJetHistos("RhoFull0",fCentralityTag);
- fRhoFull1 = new AlipAJetHistos("RhoFull1",fCentralityTag);
- fRhoFull2 = new AlipAJetHistos("RhoFull2",fCentralityTag);
- fRhoFullN = new AlipAJetHistos("RhoFullN",fCentralityTag);
- fRhoFullDijet = new AlipAJetHistos("RhoFullDijet",fCentralityTag);
- fRhoFullkT = new AlipAJetHistos("RhoFullkT",fCentralityTag);
- fRhoFullCMS = new AlipAJetHistos("RhoFullCMS",fCentralityTag);
- fRhoCharged0 = new AlipAJetHistos("RhoCharged0",fCentralityTag);
- fRhoCharged1 = new AlipAJetHistos("RhoCharged1",fCentralityTag);
- fRhoCharged2 = new AlipAJetHistos("RhoCharged2",fCentralityTag);
- fRhoChargedN = new AlipAJetHistos("RhoChargedN",fCentralityTag);
- fRhoChargedkT = new AlipAJetHistos("RhoChargedkT",fCentralityTag);
- fRhoChargedkTScale = new AlipAJetHistos("RhoChargedkTScale",fCentralityTag);
- fRhoChargedCMS = new AlipAJetHistos("RhoChargedCMS",fCentralityTag);
- */
- fRhoChargedScale = new AlipAJetHistos("RhoChargedScale",fCentralityTag);
- fRhoChargedCMSScale = new AlipAJetHistos("RhoChargedCMSScale",fCentralityTag,fDoNEF);
-
- fOutput->Add(fhTrackPt);
- fOutput->Add(fhTrackEta);
- fOutput->Add(fhTrackPhi);
- fOutput->Add(fhTrackEtaPhi);
- fOutput->Add(fhTrackPhiPt);
- fOutput->Add(fhTrackEtaPt);
- fOutput->Add(fhTrackEtaPhiPt);
- fOutput->Add(fhGlobalTrackPt);
- fOutput->Add(fhGlobalTrackEta);
- fOutput->Add(fhGlobalTrackPhi);
- fOutput->Add(fhGlobalTrackEtaPhi);
- fOutput->Add(fhGlobalTrackPhiPt);
- fOutput->Add(fhGlobalTrackEtaPt);
- fOutput->Add(fhGlobalTrackEtaPhiPt);
- fOutput->Add(fhComplementaryTrackPt);
- fOutput->Add(fhComplementaryTrackEta);
- fOutput->Add(fhComplementaryTrackPhi);
- fOutput->Add(fhComplementaryTrackEtaPhi);
- fOutput->Add(fhComplementaryTrackPhiPt);
- fOutput->Add(fhComplementaryTrackEtaPt);
- fOutput->Add(fhComplementaryTrackEtaPhiPt);
- fOutput->Add(fhClusterPt);
- fOutput->Add(fhClusterEta);
- fOutput->Add(fhClusterPhi);
- fOutput->Add(fhClusterEtaPhi);
- fOutput->Add(fhClusterPhiPt);
- fOutput->Add(fhClusterEtaPt);
- fOutput->Add(fhClusterEtaPhiPt);
fOutput->Add(fhCentrality);
- fOutput->Add(fhEMCalCellCounts);
- fOutput->Add(fhRhoScale);
-
- fOutput->Add(fpTPCEventMult);
- fOutput->Add(fpEMCalEventMult);
- fOutput->Add(fpRhoScale);
-
- fOutput->Add(fpTrackPtProfile);
- fOutput->Add(fpClusterPtProfile);
-
- fOutput->Add(fTPCRawJets->GetOutputHistos());
- fOutput->Add(fEMCalRawJets->GetOutputHistos());
-
- /*
- fOutput->Add(fRhoFull0->GetOutputHistos());
- fOutput->Add(fRhoFull1->GetOutputHistos());
- fOutput->Add(fRhoFull2->GetOutputHistos());
- fOutput->Add(fRhoFullN->GetOutputHistos());
- fOutput->Add(fRhoFullDijet->GetOutputHistos());
- fOutput->Add(fRhoFullkT->GetOutputHistos());
- fOutput->Add(fRhoFullCMS->GetOutputHistos());
- fOutput->Add(fRhoCharged0->GetOutputHistos());
- fOutput->Add(fRhoCharged1->GetOutputHistos());
- fOutput->Add(fRhoCharged2->GetOutputHistos());
- fOutput->Add(fRhoChargedN->GetOutputHistos());
- fOutput->Add(fRhoChargedkT->GetOutputHistos());
- fOutput->Add(fRhoChargedkTScale->GetOutputHistos());
- fOutput->Add(fRhoChargedCMS->GetOutputHistos());
- */
- fOutput->Add(fRhoChargedScale->GetOutputHistos());
- fOutput->Add(fRhoChargedCMSScale->GetOutputHistos());
+
+ if (fMCPartLevel == kFALSE)
+ {
+ fhChargeAndNeutralEvents = new TH1F("fhChargeAndNeutralEvents","Events with n_{tracks}>0 & n_{clusters}>0 vs Centrality",fCentralityBins*CentralityBinMult,fCentralityLow,fCentralityUp);
+ fhChargeAndNeutralEvents->GetXaxis()->SetTitle(fCentralityTag.Data());
+ fhChargeAndNeutralEvents->GetYaxis()->SetTitle("1/N_{Events}");
+ fhChargeAndNeutralEvents->Sumw2();
+
+ fhChargeOnlyEvents = new TH1F("fhChargeOnlyEvents","Events with n_{tracks}>0 & n_{clusters}=0 vs Centrality",fCentralityBins*CentralityBinMult,fCentralityLow,fCentralityUp);
+ fhChargeOnlyEvents->GetXaxis()->SetTitle(fCentralityTag.Data());
+ fhChargeOnlyEvents->GetYaxis()->SetTitle("1/N_{Events}");
+ fhChargeOnlyEvents->Sumw2();
+
+ fhNeutralOnlyEvents = new TH1F("fhNeutralOnlyEvents","Events with n_{tracks}=0 & n_{clusters}>0 vs Centrality",fCentralityBins*CentralityBinMult,fCentralityLow,fCentralityUp);
+ fhNeutralOnlyEvents->GetXaxis()->SetTitle(fCentralityTag.Data());
+ fhNeutralOnlyEvents->GetYaxis()->SetTitle("1/N_{Events}");
+ fhNeutralOnlyEvents->Sumw2();
+
+ fhNothingEvents = new TH1F("fhNothingEvents","Events with n_{tracks}=n_{clusters}=0 vs Centrality",fCentralityBins*CentralityBinMult,fCentralityLow,fCentralityUp);
+ fhNothingEvents->GetXaxis()->SetTitle(fCentralityTag.Data());
+ fhNothingEvents->GetYaxis()->SetTitle("1/N_{Events}");
+ fhNothingEvents->Sumw2();
+
+ fhEMCalChargeAndNeutralEvents = new TH1F("fhEMCalChargeAndNeutralEvents","Events with n_{emcal tracks}>0 & n_{clusters}>0 vs Centrality",fCentralityBins*CentralityBinMult,fCentralityLow,fCentralityUp);
+ fhEMCalChargeAndNeutralEvents->GetXaxis()->SetTitle(fCentralityTag.Data());
+ fhEMCalChargeAndNeutralEvents->GetYaxis()->SetTitle("1/N_{Events}");
+ fhEMCalChargeAndNeutralEvents->Sumw2();
+
+ fhEMCalChargeOnlyEvents = new TH1F("fhEMCalChargeOnlyEvents","Events with n_{emcal tracks}>0 & n_{clusters}=0 vs Centrality",fCentralityBins*CentralityBinMult,fCentralityLow,fCentralityUp);
+ fhEMCalChargeOnlyEvents->GetXaxis()->SetTitle(fCentralityTag.Data());
+ fhEMCalChargeOnlyEvents->GetYaxis()->SetTitle("1/N_{Events}");
+ fhEMCalChargeOnlyEvents->Sumw2();
+
+ fhEMCalNeutralOnlyEvents = new TH1F("fhEMCalNeutralOnlyEvents","Events with n_{tracks}>0 & n_{emcal tracks}=0 & n_{clusters}>0 vs Centrality",fCentralityBins*CentralityBinMult,fCentralityLow,fCentralityUp);
+ fhEMCalNeutralOnlyEvents->GetXaxis()->SetTitle(fCentralityTag.Data());
+ fhEMCalNeutralOnlyEvents->GetYaxis()->SetTitle("1/N_{Events}");
+ fhEMCalNeutralOnlyEvents->Sumw2();
+
+ fhEMCalNothingEvents = new TH1F("fhEMCalNothingEvents","Events with n_{emcal tracks}=n_{clusters}=0 vs Centrality",fCentralityBins*CentralityBinMult,fCentralityLow,fCentralityUp);
+ fhEMCalNothingEvents->GetXaxis()->SetTitle(fCentralityTag.Data());
+ fhEMCalNothingEvents->GetYaxis()->SetTitle("1/N_{Events}");
+ fhEMCalNothingEvents->Sumw2();
+
+ fOutput->Add(fhChargeAndNeutralEvents);
+ fOutput->Add(fhChargeOnlyEvents);
+ fOutput->Add(fhNeutralOnlyEvents);
+ fOutput->Add(fhNothingEvents);
+ fOutput->Add(fhEMCalChargeAndNeutralEvents);
+ fOutput->Add(fhEMCalChargeOnlyEvents);
+ fOutput->Add(fhEMCalNeutralOnlyEvents);
+ fOutput->Add(fhEMCalNothingEvents);
+ }
// Post data for ALL output slots >0 here,
// To get at least an empty histogram
if (esd)
{
- fEventCentrality=esd->GetCentrality()->GetCentralityPercentile(fCentralityTag);
+ fEventCentrality=esd->GetCentrality()->GetCentralityPercentile(fCentralityTag.Data());
- if (esd->GetPrimaryVertex()->GetNContributors()<1 || (TMath::Abs(esd->GetPrimaryVertex()->GetZ())>fVertexWindow))
- {
- return;
- }
- if (TMath::Sqrt(TMath::Power(esd->GetPrimaryVertex()->GetX(),2)+TMath::Power(esd->GetPrimaryVertex()->GetY(),2))>fVertexMaxR)
+ if (fDoVertexRCut == kTRUE)
{
- return;
+ // Vertex R cut not done in AliAnalysisTaskEmcal
+ if (TMath::Sqrt(TMath::Power(esd->GetPrimaryVertex()->GetX(),2)+TMath::Power(esd->GetPrimaryVertex()->GetY(),2))>fVertexMaxR)
+ {
+ DeleteJetData(0);
+ return;
+ }
}
esd->GetPrimaryVertex()->GetXYZ(fVertex);
}
else if (aod)
{
- fEventCentrality=aod->GetCentrality()->GetCentralityPercentile(fCentralityTag);
+ fEventCentrality=aod->GetCentrality()->GetCentralityPercentile(fCentralityTag.Data());
- if (aod->GetPrimaryVertex()->GetNContributors()<1 || (TMath::Abs(aod->GetPrimaryVertex()->GetZ())>fVertexWindow))
- {
- return;
- }
- if (TMath::Sqrt(TMath::Power(aod->GetPrimaryVertex()->GetX(),2)+TMath::Power(aod->GetPrimaryVertex()->GetY(),2))>fVertexMaxR)
+ if (fDoVertexRCut == kTRUE)
{
- return;
+ // Vertex R cut not done in AliAnalysisTaskEmcal
+ if (TMath::Sqrt(TMath::Power(aod->GetPrimaryVertex()->GetX(),2)+TMath::Power(aod->GetPrimaryVertex()->GetY(),2))>fVertexMaxR)
+ {
+ DeleteJetData(0);
+ return;
+ }
}
-
+
aod->GetPrimaryVertex()->GetXYZ(fVertex);
fCells = (AliVCaloCells*) aod->GetEMCALCells();
fnCaloClusters = aod->GetNumberOfCaloClusters();
else
{
AliError("Cannot get AOD/ESD event. Rejecting Event");
+ DeleteJetData(0);
return;
}
-
+
// Make sure Centrality isn't exactly 100% (to avoid bin filling errors in profile plots. Make it 99.99
if (fEventCentrality>99.99)
{
fEventCentrality=99.99;
}
- fhCentrality->Fill(fEventCentrality);
+ fhCentrality->Fill(fEventCentrality); // Counts total events
+
+
+ // Count Events
+ EventCounts();
- TrackCuts();
// Reject any event that doesn't have any tracks, i.e. TPC is off
if (fnTracks<1)
{
+ if (fTrackQA==kTRUE)
+ {
+ fhTPCEventMult->Fill(fEventCentrality,0.0);
+ fpTPCEventMult->Fill(fEventCentrality,0.0);
+ fhEMCalTrackEventMult->Fill(fEventCentrality,0.0);
+ }
AliWarning("No PicoTracks, Rejecting Event");
+ DeleteJetData(1);
+ PostData(1,fOutput);
return;
}
- ClusterCuts();
if (fnClusters<1)
{
AliInfo("No Corrected CaloClusters, using only charged jets");
- TrackHisto();
+ if (fTrackQA==kTRUE)
+ {
+ TrackHisto();
+ }
InitChargedJets();
GenerateTPCRandomConesPt();
+ if (fClusterQA==kTRUE)
+ {
+ fhEMCalEventMult->Fill(fEventCentrality,0.0);
+ fpEMCalEventMult->Fill(fEventCentrality,0.0);
+ }
+
// Rho's
- /*
- EstimateChargedRho0();
- EstimateChargedRho1();
- EstimateChargedRho2();
- EstimateChargedRhoN();
- EstimateChargedRhokT();
- EstimateChargedRhoCMS();
- */
-
- DeleteJetData(kFALSE);
+ if (fCalculateRhoJet>=1)
+ {
+ if (fDoJetRhoDensity == kTRUE)
+ {
+ ChargedJetEnergyDensityProfile();
+ }
+ if (fMCPartLevel==kTRUE)
+ {
+ EstimateChargedRhoCMS();
+ }
+ }
+ if (fCalculateRhoJet>=2 && fMCPartLevel==kFALSE)
+ {
+ EstimateChargedRhoCMS();
+ }
+ if (fCalculateRhoJet>=3 && fMCPartLevel==kFALSE)
+ {
+ EstimateChargedRho0();
+ EstimateChargedRho1();
+ EstimateChargedRho2();
+ EstimateChargedRhoN();
+ EstimateChargedRhokT();
+ }
+ DeleteJetData(2);
fnEventsCharged++;
-
- PostData(1, fOutput);
+ PostData(1,fOutput);
return;
}
-
- TrackHisto();
- ClusterHisto();
+
+ if (fTrackQA==kTRUE)
+ {
+ TrackHisto();
+ }
+ if (fClusterQA==kTRUE)
+ {
+ ClusterHisto();
+ }
// Prep the jets
InitChargedJets();
- InitFullJets();
GenerateTPCRandomConesPt();
+
+ InitFullJets();
GenerateEMCalRandomConesPt();
- // Rho's
- /*
- EstimateChargedRho0();
- EstimateChargedRho1();
- EstimateChargedRho2();
- EstimateChargedRhoN();
- EstimateChargedRhokT();
- EstimateChargedRhoCMS();
- */
-
- /*
- EstimateFullRho0();
- EstimateFullRho1();
- EstimateFullRho2();
- EstimateFullRhoN();
- EstimateFullRhokT();
- EstimateFullRhoCMS();
-
- EstimateChargedRhokTScale();
- */
- EstimateChargedRhoScale();
- EstimateChargedRhoCMSScale();
-
- // Dijet
- if (IsDiJetEvent()==kTRUE)
+ if (fCalculateRhoJet>=0)
{
- //EstimateFullRhoDijet();
+ EstimateChargedRhoCMSScale();
}
-
+ if (fCalculateRhoJet>=1)
+ {
+ EstimateChargedRhoScale();
+ if (fDoJetRhoDensity == kTRUE)
+ {
+ ChargedJetEnergyDensityProfile();
+ FullJetEnergyDensityProfile();
+ }
+ }
+ if (fCalculateRhoJet>=2)
+ {
+ EstimateChargedRhoCMS();
+ }
+ if (fCalculateRhoJet>=3)
+ {
+ EstimateChargedRho0();
+ EstimateChargedRho1();
+ EstimateChargedRho2();
+ EstimateChargedRhoN();
+ EstimateChargedRhokT();
+
+ EstimateFullRho0();
+ EstimateFullRho1();
+ EstimateFullRho2();
+ EstimateFullRhoN();
+ EstimateFullRhokT();
+ EstimateFullRhoCMS();
+
+ EstimateChargedRhokTScale();
+
+ // Dijet
+ if (IsDiJetEvent()==kTRUE)
+ {
+ EstimateFullRhoDijet();
+ }
+ }
+
// Delete Dynamic Arrays
- DeleteJetData(kTRUE);
- delete fRecoUtil;
+ DeleteJetData(3);
fnEvents++;
-
- PostData(1, fOutput);
+ PostData(1,fOutput);
}
//________________________________________________________________________
{
// Fill a TObjArray with the tracks from a TClonesArray which grabs the picotracks.
Int_t i;
+ Int_t j=0;
fmyTracks = new TObjArray();
for (i=0;i<fOrgTracks->GetEntries();i++)
if (vtrack->Pt()>=fTrackMinPt)
{
fmyTracks->Add(vtrack);
+ if (IsInEMCal(vtrack->Phi(),vtrack->Eta()) == kTRUE)
+ {
+ j++;
+ }
}
}
fnTracks = fmyTracks->GetEntries();
+ fnEMCalTracks = j;
}
void AliAnalysisTaskFullpAJets::ClusterCuts()
Int_t i;
fmyClusters = new TObjArray();
- if(fOrgClusters)
+ TLorentzVector *cluster_vec = new TLorentzVector;
+ for (i=0;i<fOrgClusters->GetEntries();i++)
{
- for (i=0;i<fOrgClusters->GetEntries();i++)
+ AliVCluster* vcluster = (AliVCluster*) fOrgClusters->At(i);
+ vcluster->GetMomentum(*cluster_vec,fVertex);
+
+ if (cluster_vec->Pt()>=fClusterMinPt && vcluster->IsEMCAL()==kTRUE)
{
- AliVCluster* vcluster = (AliVCluster*) fOrgClusters->At(i);
- TLorentzVector *cluster_vec = new TLorentzVector;
- vcluster->GetMomentum(*cluster_vec,fVertex);
-
- if (cluster_vec->Pt()>=fClusterMinPt && vcluster->IsEMCAL()==kTRUE)
+ fmyClusters->Add(vcluster);
+ }
+ }
+ fnClusters = fmyClusters->GetEntries();
+ delete cluster_vec;
+}
+
+void AliAnalysisTaskFullpAJets::EventCounts()
+{
+ TrackCuts();
+ if (fMCPartLevel == kTRUE)
+ {
+ return;
+ }
+
+ ClusterCuts();
+ if (fnTracks>0)
+ {
+ if (fnClusters>0)
+ {
+ fhChargeAndNeutralEvents->Fill(fEventCentrality);
+ }
+ else
+ {
+ fhChargeOnlyEvents->Fill(fEventCentrality);
+ }
+ if (fnEMCalTracks>0)
+ {
+ if (fnClusters>0)
{
- fmyClusters->Add(vcluster);
+ fhEMCalChargeAndNeutralEvents->Fill(fEventCentrality);
+ }
+ else
+ {
+ fhEMCalChargeOnlyEvents->Fill(fEventCentrality);
+ }
+ }
+ else
+ {
+ if (fnClusters>0)
+ {
+ fhEMCalNeutralOnlyEvents->Fill(fEventCentrality);
+ }
+ else
+ {
+ fhEMCalNothingEvents->Fill(fEventCentrality);
}
- delete cluster_vec;
-
}
}
- fnClusters = fmyClusters->GetEntries();
+ else
+ {
+ if (fnClusters>0)
+ {
+ fhNeutralOnlyEvents->Fill(fEventCentrality);
+ }
+ else
+ {
+ fhNothingEvents->Fill(fEventCentrality);
+ }
+ }
}
void AliAnalysisTaskFullpAJets::TrackHisto()
// Fill track histograms: Phi,Eta,Pt
Int_t i,j;
Int_t TCBins=100;
- TH2D *hdummypT= new TH2D("hdummypT","",TCBins,fTPCEtaMin,fTPCEtaMax,TCBins,fTPCPhiMin,fTPCPhiMax); //!
+ TH2F *hdummypT= new TH2F("hdummypT","",TCBins,fTPCEtaMin,fTPCEtaMax,TCBins,fTPCPhiMin,fTPCPhiMax); //!
for (i=0;i<fnTracks;i++)
{
- AliPicoTrack* vtrack =(AliPicoTrack*) fmyTracks->At(i);
+ AliPicoTrack* vtrack = (AliPicoTrack*) fmyTracks->At(i);
fhTrackPt->Fill(vtrack->Pt());
fhTrackEta->Fill(vtrack->Eta());
fhTrackPhi->Fill(vtrack->Phi());
fhTrackEtaPhi->Fill(vtrack->Eta(),vtrack->Phi());
fhTrackPhiPt->Fill(vtrack->Phi(),vtrack->Pt());
fhTrackEtaPt->Fill(vtrack->Eta(),vtrack->Pt());
- fhTrackEtaPhiPt->Fill(vtrack->Eta(),vtrack->Phi(),vtrack->Pt());
// Fill Associated Track Distributions for AOD QA Productions
// Global Tracks
fhGlobalTrackEtaPhi->Fill(vtrack->Eta(),vtrack->Phi());
fhGlobalTrackPhiPt->Fill(vtrack->Phi(),vtrack->Pt());
fhGlobalTrackEtaPt->Fill(vtrack->Eta(),vtrack->Pt());
- fhGlobalTrackEtaPhiPt->Fill(vtrack->Eta(),vtrack->Phi(),vtrack->Pt());
}
// Complementary Tracks
else if (vtrack->GetTrackType()==1)
fhComplementaryTrackEtaPhi->Fill(vtrack->Eta(),vtrack->Phi());
fhComplementaryTrackPhiPt->Fill(vtrack->Phi(),vtrack->Pt());
fhComplementaryTrackEtaPt->Fill(vtrack->Eta(),vtrack->Pt());
- fhComplementaryTrackEtaPhiPt->Fill(vtrack->Eta(),vtrack->Phi(),vtrack->Pt());
}
-
hdummypT->Fill(vtrack->Eta(),vtrack->Phi(),vtrack->Pt());
}
for (i=1;i<=TCBins;i++)
// Fill cluster histograms: Phi,Eta,Pt
Int_t i,j;
Int_t TCBins=100;
- TH2D *hdummypT= new TH2D("hdummypT","",TCBins,fEMCalEtaMin,fEMCalEtaMax,TCBins,fEMCalPhiMin,fEMCalPhiMax); //!
+ TH2F *hdummypT= new TH2F("hdummypT","",TCBins,fEMCalEtaMin,fEMCalEtaMax,TCBins,fEMCalPhiMin,fEMCalPhiMax); //!
Int_t myCellID=-2;
+ TLorentzVector *cluster_vec = new TLorentzVector;
for (i=0;i<fnClusters;i++)
{
AliVCluster* vcluster = (AliVCluster*) fmyClusters->At(i);
- TLorentzVector *cluster_vec = new TLorentzVector;
vcluster->GetMomentum(*cluster_vec,fVertex);
fhClusterPt->Fill(cluster_vec->Pt());
fhClusterEtaPhi->Fill(cluster_vec->Eta(),cluster_vec->Phi());
fhClusterPhiPt->Fill(cluster_vec->Phi(),cluster_vec->Pt());
fhClusterEtaPt->Fill(cluster_vec->Eta(),cluster_vec->Pt());
- fhClusterEtaPhiPt->Fill(cluster_vec->Eta(),cluster_vec->Phi(),cluster_vec->Pt());
hdummypT->Fill(cluster_vec->Eta(),cluster_vec->Phi(),cluster_vec->Pt());
fEMCALGeometry->GetAbsCellIdFromEtaPhi(cluster_vec->Eta(),cluster_vec->Phi(),myCellID);
fhEMCalCellCounts->Fill(myCellID);
- delete cluster_vec;
}
for (i=1;i<=TCBins;i++)
{
}
}
delete hdummypT;
+ delete cluster_vec;
}
void AliAnalysisTaskFullpAJets::InitChargedJets()
fTPCFullJet = new AlipAJetData("fTPCFullJet",kFALSE,fnAKTChargedJets);
fTPCOnlyJet = new AlipAJetData("fTPCOnlyJet",kFALSE,fnAKTChargedJets);
fTPCJetUnbiased = new AlipAJetData("fTPCJetUnbiased",kFALSE,fnAKTChargedJets);
-
+
fTPCJet->SetSignalCut(fTPCJetThreshold);
fTPCJet->SetAreaCutFraction(fJetAreaCutFrac);
fTPCJet->SetJetR(fJetR);
+ fTPCJet->SetSignalTrackPtBias(fSignalTrackBias);
fTPCFullJet->SetSignalCut(fTPCJetThreshold);
fTPCFullJet->SetAreaCutFraction(fJetAreaCutFrac);
fTPCFullJet->SetJetR(fJetR);
+ fTPCFullJet->SetSignalTrackPtBias(fSignalTrackBias);
fTPCOnlyJet->SetSignalCut(fTPCJetThreshold);
fTPCOnlyJet->SetAreaCutFraction(fJetAreaCutFrac);
fTPCOnlyJet->SetJetR(fJetR);
+ fTPCOnlyJet->SetSignalTrackPtBias(fSignalTrackBias);
fTPCJetUnbiased->SetSignalCut(fTPCJetThreshold);
fTPCJetUnbiased->SetAreaCutFraction(fJetAreaCutFrac);
fTPCJetUnbiased->SetJetR(fJetR);
AliEmcalJet *myJet =(AliEmcalJet*) fmyAKTChargedJets->At(i);
fTPCJet->SetIsJetInArray(IsInTPC(fJetR,myJet->Phi(),myJet->Eta(),kFALSE),i);
- fTPCFullJet->SetIsJetInArray(IsInTPC(fJetR,myJet->Phi(),myJet->Eta(),kTRUE),i);
+ fTPCFullJet->SetIsJetInArray(IsInTPC(fJetRAccept,myJet->Phi(),myJet->Eta(),kTRUE),i);
fTPCOnlyJet->SetIsJetInArray(IsInTPCFull(fJetR,myJet->Phi(),myJet->Eta()),i);
fTPCJetUnbiased->SetIsJetInArray(IsInTPC(fJetR,myJet->Phi(),myJet->Eta(),kFALSE),i);
}
+
fTPCJet->InitializeJetData(fmyAKTChargedJets,fnAKTChargedJets);
fTPCFullJet->InitializeJetData(fmyAKTChargedJets,fnAKTChargedJets);
fTPCOnlyJet->InitializeJetData(fmyAKTChargedJets,fnAKTChargedJets);
fTPCkTFullJet->InitializeJetData(fmyKTChargedJets,fnKTChargedJets);
// Raw Charged Jet Spectra
- fTPCRawJets->FillBSJS(fEventCentrality,0.0,fTPCJetThreshold,fmyAKTChargedJets,fTPCFullJet->GetJets(),fTPCFullJet->GetTotalJets());
+ if (fCalculateRhoJet>=2 || (fCalculateRhoJet>=1 && fMCPartLevel==kTRUE))
+ {
+ fTPCRawJets->FillBSJS(fEventCentrality,0.0,fTPCJetThreshold,fmyAKTChargedJets,fTPCFullJet->GetJets(),fTPCFullJet->GetTotalJets());
+ }
}
void AliAnalysisTaskFullpAJets::InitFullJets()
fEMCalJet->SetAreaCutFraction(fJetAreaCutFrac);
fEMCalJet->SetJetR(fJetR);
fEMCalJet->SetNEF(fNEFSignalJetCut);
- fEMCalJet->SetSignalTrackPtBias(kTRUE);
+ fEMCalJet->SetSignalTrackPtBias(fSignalTrackBias);
fEMCalFullJet->SetSignalCut(fEMCalJetThreshold);
fEMCalFullJet->SetAreaCutFraction(fJetAreaCutFrac);
fEMCalFullJet->SetJetR(fJetR);
fEMCalFullJet->SetNEF(fNEFSignalJetCut);
- fEMCalFullJet->SetSignalTrackPtBias(kTRUE);
+ fEMCalFullJet->SetSignalTrackPtBias(fSignalTrackBias);
fEMCalPartJet->SetSignalCut(fEMCalJetThreshold);
fEMCalPartJet->SetAreaCutFraction(fJetAreaCutFrac);
fEMCalPartJet->SetJetR(fJetR);
fEMCalPartJet->SetNEF(fNEFSignalJetCut);
- fEMCalPartJet->SetSignalTrackPtBias(kTRUE);
+ fEMCalPartJet->SetSignalTrackPtBias(fSignalTrackBias);
fEMCalPartJetUnbiased->SetSignalCut(fEMCalJetThreshold);
fEMCalPartJetUnbiased->SetAreaCutFraction(fJetAreaCutFrac);
fEMCalPartJetUnbiased->SetJetR(fJetR);
AliEmcalJet *myJet =(AliEmcalJet*) fmyAKTFullJets->At(i);
fEMCalJet->SetIsJetInArray(IsInEMCal(myJet->Phi(),myJet->Eta()),i);
- fEMCalFullJet->SetIsJetInArray(IsInEMCalFull(fJetR,myJet->Phi(),myJet->Eta()),i);
+ fEMCalFullJet->SetIsJetInArray(IsInEMCalFull(fJetRAccept,myJet->Phi(),myJet->Eta()),i);
fEMCalPartJet->SetIsJetInArray(IsInEMCalPart(fJetR,myJet->Phi(),myJet->Eta()),i);
fEMCalPartJetUnbiased->SetIsJetInArray(IsInEMCalPart(fJetR,myJet->Phi(),myJet->Eta()),i);
}
fEMCalkTFullJet->SetAreaCutFraction(0.25*fJetAreaCutFrac);
fEMCalkTFullJet->SetJetR(fJetR);
fEMCalkTFullJet->SetNEF(fNEFSignalJetCut);
- fEMCalkTFullJet->SetSignalTrackPtBias(kTRUE);
+ fEMCalkTFullJet->SetSignalTrackPtBias(fSignalTrackBias);
for (i=0;i<fnKTFullJets;i++)
{
fEMCalkTFullJet->InitializeJetData(fmyKTFullJets,fnKTFullJets);
// Raw Full Jet Spectra
- fEMCalRawJets->FillBSJS(fEventCentrality,0.0,fEMCalJetThreshold,fmyAKTFullJets,fEMCalFullJet->GetJets(),fEMCalFullJet->GetTotalJets());
+ if (fMCPartLevel==kFALSE)
+ {
+ fEMCalRawJets->FillBSJS(fEventCentrality,0.0,fEMCalJetThreshold,fmyAKTFullJets,fEMCalFullJet->GetJets(),fEMCalFullJet->GetTotalJets());
+ }
}
void AliAnalysisTaskFullpAJets::GenerateTPCRandomConesPt()
Double_t Eta_Center=0.5*(fTPCEtaMin+fTPCEtaMax);
Double_t Phi_Center=0.5*(fTPCPhiMin+fTPCPhiMax);
Int_t event_mult=0;
+ Int_t event_track_mult=0;
Int_t clus_mult=0;
for (i=0;i<fnBckgClusters;i++)
TLorentzVector *dummy= new TLorentzVector;
TLorentzVector *temp_jet= new TLorentzVector;
-
+ TLorentzVector *track_vec = new TLorentzVector;
+
// First, consider the RC with no spatial restrictions
for (j=0;j<fnBckgClusters;j++)
{
AliVTrack* vtrack =(AliVTrack*) fmyTracks->At(i);
if (IsInTPC(fJetR,vtrack->Phi(),vtrack->Eta(),kFALSE)==kTRUE)
{
- TLorentzVector *track_vec = new TLorentzVector;
track_vec->SetPtEtaPhiE(vtrack->Pt(),vtrack->Eta(),vtrack->Phi(),vtrack->E());
if (dummy->DeltaR(*track_vec)<fJetR)
{
E_tracks_total+=vtrack->Pt();
}
- delete track_vec;
}
}
fTPCRCBckgFlucSignal[j]=E_tracks_total;
for (j=0;j<fnBckgClusters;j++)
{
event_mult=0;
+ event_track_mult=0;
clus_mult=0;
E_tracks_total=0.;
for (i=0;i<fnTracks;i++)
{
AliVTrack* vtrack =(AliVTrack*) fmyTracks->At(i);
- if (IsInTPC(fJetR,vtrack->Phi(),vtrack->Eta(),kFALSE)==kTRUE)
+ if (IsInTPC(fJetR,vtrack->Phi(),vtrack->Eta(),kFALSE) == kTRUE)
{
event_mult++;
- TLorentzVector *track_vec = new TLorentzVector;
+ if (IsInEMCal(vtrack->Phi(),vtrack->Eta()) == kTRUE)
+ {
+ event_track_mult++;
+ }
track_vec->SetPtEtaPhiE(vtrack->Pt(),vtrack->Eta(),vtrack->Phi(),vtrack->E());
if (dummy->DeltaR(*track_vec)<fJetR)
{
clus_mult++;
E_tracks_total+=vtrack->Pt();
}
- delete track_vec;
}
}
fTPCRCBckgFluc[j]=E_tracks_total;
}
- fpTPCEventMult->Fill(fEventCentrality,event_mult);
- fTPCRawJets->FillDeltaPt(fEventCentrality,0.0,fJetR,fTPCRCBckgFluc,1);
+ if (fTrackQA==kTRUE)
+ {
+ fhTPCEventMult->Fill(fEventCentrality,event_mult);
+ fpTPCEventMult->Fill(fEventCentrality,event_mult);
+ fhEMCalTrackEventMult->Fill(fEventCentrality,event_track_mult);
+ }
+ if (fCalculateRhoJet>=2 || (fCalculateRhoJet>=1 && fMCPartLevel==kTRUE))
+ {
+ fTPCRawJets->FillDeltaPt(fEventCentrality,0.0,fJetR,fTPCRCBckgFluc,1);
+ }
// For the case of partial exclusion, merely allow a superposition of full and no exclusion with probability p=1/Ncoll
Double_t exclusion_prob;
delete dummy;
delete temp_jet;
+ delete track_vec;
}
void AliAnalysisTaskFullpAJets::GenerateEMCalRandomConesPt()
TLorentzVector *dummy= new TLorentzVector;
TLorentzVector *temp_jet= new TLorentzVector;
-
+ TLorentzVector *track_vec = new TLorentzVector;
+ TLorentzVector *cluster_vec = new TLorentzVector;
+
// First, consider the RC with no spatial restrictions
for (j=0;j<fnBckgClusters;j++)
{
AliVTrack* vtrack =(AliVTrack*) fmyTracks->At(i);
if (IsInEMCal(vtrack->Phi(),vtrack->Eta())==kTRUE)
{
- TLorentzVector *track_vec = new TLorentzVector;
track_vec->SetPtEtaPhiE(vtrack->Pt(),vtrack->Eta(),vtrack->Phi(),vtrack->E());
if (dummy->DeltaR(*track_vec)<fJetR)
{
E_tracks_total+=vtrack->Pt();
}
- delete track_vec;
}
}
for (i=0;i<fnClusters;i++)
{
AliVCluster* vcluster = (AliVCluster*) fmyClusters->At(i);
- TLorentzVector *cluster_vec = new TLorentzVector;
vcluster->GetMomentum(*cluster_vec,fVertex);
if (dummy->DeltaR(*cluster_vec)<fJetR)
{
clus_mult++;
- E_caloclusters_total+=vcluster->E();
+ E_caloclusters_total+=cluster_vec->Pt();
}
- delete cluster_vec;
}
fEMCalRCBckgFlucSignal[j]=E_tracks_total+E_caloclusters_total;
}
if (IsInEMCal(vtrack->Phi(),vtrack->Eta())==kTRUE)
{
event_mult++;
- TLorentzVector *track_vec = new TLorentzVector;
track_vec->SetPtEtaPhiE(vtrack->Pt(),vtrack->Eta(),vtrack->Phi(),vtrack->E());
if (dummy->DeltaR(*track_vec)<fJetR)
{
clus_mult++;
E_tracks_total+=vtrack->Pt();
}
- delete track_vec;
}
}
for (i=0;i<fnClusters;i++)
{
AliVCluster* vcluster = (AliVCluster*) fmyClusters->At(i);
- TLorentzVector *cluster_vec = new TLorentzVector;
vcluster->GetMomentum(*cluster_vec,fVertex);
event_mult++;
if (dummy->DeltaR(*cluster_vec)<fJetR)
{
clus_mult++;
- E_caloclusters_total+=vcluster->E();
+ E_caloclusters_total+=cluster_vec->Pt();
}
- delete cluster_vec;
}
fEMCalRCBckgFluc[j]=E_tracks_total+E_caloclusters_total;
}
- fpEMCalEventMult->Fill(fEventCentrality,event_mult);
+ if (fClusterQA==kTRUE)
+ {
+ fhEMCalEventMult->Fill(fEventCentrality,event_mult);
+ fpEMCalEventMult->Fill(fEventCentrality,event_mult);
+ }
fEMCalRawJets->FillDeltaPt(fEventCentrality,0.0,fJetR,fEMCalRCBckgFluc,1);
// For the case of partial exclusion, merely allow a superposition of full and no exclusion with probability p=1/Ncoll
delete dummy;
delete temp_jet;
+ delete track_vec;
+ delete cluster_vec;
}
-/*
// Charged Rho's
void AliAnalysisTaskFullpAJets::EstimateChargedRho0()
{
Int_t i;
Double_t E_tracks_total=0.0;
- Double_t TPC_rho=0.;
+ Double_t TPC_rho=0.0;
// Loop over all tracks
for (i=0;i<fnTracks;i++)
// Calculate the mean Background density
TPC_rho=E_tracks_total/fTPCArea;
- fRhoCharged=TPC_rho;
+
+ fRhoCharged = TPC_rho;
// Fill Histograms
fRhoCharged0->FillRho(fEventCentrality,TPC_rho);
Double_t E_tracks_total=0.0;
Double_t TPC_rho=0.;
+ TLorentzVector *temp_jet= new TLorentzVector;
+ TLorentzVector *track_vec = new TLorentzVector;
+
if (fTPCJetUnbiased->GetLeadingPt()>=fTPCJetThreshold)
{
AliEmcalJet *myJet =(AliEmcalJet*) fmyAKTChargedJets->At(fTPCJet->GetLeadingIndex());
- TLorentzVector *temp_jet= new TLorentzVector;
myJet->GetMom(*temp_jet);
// Loop over all tracks
AliVTrack* vtrack =(AliVTrack*) fmyTracks->At(i);
if (IsInTPC(fJetR,vtrack->Phi(),vtrack->Eta(),kFALSE)==kTRUE)
{
- TLorentzVector *track_vec = new TLorentzVector;
track_vec->SetPtEtaPhiE(vtrack->Pt(),vtrack->Eta(),vtrack->Phi(),vtrack->E());
if (temp_jet->DeltaR(*track_vec)>fJetRForRho)
{
E_tracks_total+=vtrack->Pt();
}
- delete track_vec;
}
}
- delete temp_jet;
// Calculate the mean Background density
TPC_rho=E_tracks_total/(fTPCArea-AreaWithinTPC(fJetR,myJet->Eta()));
// Calculate the mean Background density
TPC_rho=E_tracks_total/fTPCArea;
}
-
+ delete track_vec;
+ delete temp_jet;
+
// Fill histograms
fRhoCharged1->FillRho(fEventCentrality,TPC_rho);
fRhoCharged1->FillBSJS(fEventCentrality,TPC_rho,fTPCJetThreshold,fmyAKTChargedJets,fTPCFullJet->GetJets(),fTPCFullJet->GetTotalJets());
Int_t i;
Double_t E_tracks_total=0.0;
Double_t TPC_rho=0.;
-
+
+ TLorentzVector *temp_jet1= new TLorentzVector;
+ TLorentzVector *temp_jet2= new TLorentzVector;
+ TLorentzVector *track_vec = new TLorentzVector;
+
if ((fTPCJetUnbiased->GetLeadingPt()>=fTPCJetThreshold) && (fTPCJetUnbiased->GetSubLeadingPt()>=fTPCJetThreshold))
{
AliEmcalJet *myhJet =(AliEmcalJet*) fmyAKTChargedJets->At(fTPCJet->GetLeadingIndex());
- TLorentzVector *temp_jet1= new TLorentzVector;
myhJet->GetMom(*temp_jet1);
AliEmcalJet *myJet =(AliEmcalJet*) fmyAKTChargedJets->At(fTPCJet->GetSubLeadingIndex());
- TLorentzVector *temp_jet2= new TLorentzVector;
myJet->GetMom(*temp_jet2);
// Loop over all tracks
AliVTrack* vtrack =(AliVTrack*) fmyTracks->At(i);
if (IsInTPC(fJetR,vtrack->Phi(),vtrack->Eta(),kFALSE)==kTRUE)
{
- TLorentzVector *track_vec = new TLorentzVector;
track_vec->SetPtEtaPhiE(vtrack->Pt(),vtrack->Eta(),vtrack->Phi(),vtrack->E());
if ((temp_jet1->DeltaR(*track_vec)>fJetRForRho) && (temp_jet2->DeltaR(*track_vec)>fJetRForRho))
{
E_tracks_total+=vtrack->Pt();
}
- delete track_vec;
}
}
- delete temp_jet1;
- delete temp_jet2;
// Calculate the mean Background density
TPC_rho=E_tracks_total/(fTPCArea-AreaWithinTPC(fJetR,myhJet->Eta())-AreaWithinTPC(fJetR,myJet->Eta()));
else if (fTPCJetUnbiased->GetLeadingPt()>=fTPCJetThreshold)
{
AliEmcalJet *myJet =(AliEmcalJet*) fmyAKTChargedJets->At(fTPCJet->GetLeadingIndex());
- TLorentzVector *temp_jet= new TLorentzVector;
- myJet->GetMom(*temp_jet);
+ myJet->GetMom(*temp_jet1);
// Loop over all tracks
for (i=0;i<fnTracks;i++)
AliVTrack* vtrack =(AliVTrack*) fmyTracks->At(i);
if (IsInTPC(fJetR,vtrack->Phi(),vtrack->Eta(),kFALSE)==kTRUE)
{
- TLorentzVector *track_vec = new TLorentzVector;
track_vec->SetPtEtaPhiE(vtrack->Pt(),vtrack->Eta(),vtrack->Phi(),vtrack->E());
- if (temp_jet->DeltaR(*track_vec)>fJetRForRho)
+ if (temp_jet1->DeltaR(*track_vec)>fJetRForRho)
{
E_tracks_total+=vtrack->Pt();
}
- delete track_vec;
}
}
- delete temp_jet;
// Calculate the mean Background density
TPC_rho=E_tracks_total/(fTPCArea-AreaWithinTPC(fJetR,myJet->Eta()));
// Calculate the mean Background density
TPC_rho=E_tracks_total/fTPCArea;
}
-
+ delete temp_jet1;
+ delete temp_jet2;
+ delete track_vec;
+
// Fill histograms
fRhoCharged2->FillRho(fEventCentrality,TPC_rho);
fRhoCharged2->FillBSJS(fEventCentrality,TPC_rho,fTPCJetThreshold,fmyAKTChargedJets,fTPCFullJet->GetJets(),fTPCFullJet->GetTotalJets());
Double_t TPC_rho=0.0;
Double_t jet_area_total=0.0;
+ TLorentzVector *jet_vec= new TLorentzVector;
+ TLorentzVector *track_vec = new TLorentzVector;
+
// First, sum all tracks within the EMCal that are away from jet(s) above Pt Threshold
for (i=0;i<fnTracks;i++)
{
{
track_away_from_jet=kTRUE;
j=0;
- TLorentzVector *track_vec = new TLorentzVector;
track_vec->SetPtEtaPhiE(vtrack->Pt(),vtrack->Eta(),vtrack->Phi(),vtrack->E());
while (track_away_from_jet==kTRUE && j<fTPCJetUnbiased->GetTotalSignalJets())
{
- TLorentzVector *jet_vec= new TLorentzVector;
AliEmcalJet *myJet =(AliEmcalJet*) fmyAKTChargedJets->At(fTPCJetUnbiased->GetSignalJetIndex(j));
myJet->GetMom(*jet_vec);
if (track_vec->DeltaR(*jet_vec)<=fJetRForRho)
{
track_away_from_jet=kFALSE;
}
- delete jet_vec;
j++;
}
if (track_away_from_jet==kTRUE)
{
E_tracks_total+=vtrack->Pt();
}
- delete track_vec;
}
}
}
jet_area_total+=AreaWithinTPC(fJetR,myJet->Eta());
}
}
-
+ delete jet_vec;
+ delete track_vec;
+
// Calculate Rho
TPC_rho = E_tracks_total/(fTPCArea-jet_area_total);
fRhoChargedN->FillLeadingJetPtRho(fTPCFullJet->GetLeadingPt(),TPC_rho);
}
-*/
void AliAnalysisTaskFullpAJets::EstimateChargedRhoScale()
{
Double_t TPC_rho=0.0;
Double_t jet_area_total=0.0;
+ TLorentzVector *jet_vec= new TLorentzVector;
+ TLorentzVector *track_vec = new TLorentzVector;
+
// First, sum all tracks within the EMCal that are away from jet(s) above Pt Threshold
for (i=0;i<fnTracks;i++)
{
{
track_away_from_jet=kTRUE;
j=0;
- TLorentzVector *track_vec = new TLorentzVector;
track_vec->SetPtEtaPhiE(vtrack->Pt(),vtrack->Eta(),vtrack->Phi(),vtrack->E());
while (track_away_from_jet==kTRUE && j<fTPCJetUnbiased->GetTotalSignalJets())
{
- TLorentzVector *jet_vec= new TLorentzVector;
AliEmcalJet *myJet =(AliEmcalJet*) fmyAKTChargedJets->At(fTPCJetUnbiased->GetSignalJetIndex(j));
myJet->GetMom(*jet_vec);
if (track_vec->DeltaR(*jet_vec)<=fJetRForRho)
{
track_away_from_jet=kFALSE;
}
- delete jet_vec;
j++;
}
if (track_away_from_jet==kTRUE)
{
E_tracks_total+=vtrack->Pt();
}
- delete track_vec;
}
}
}
jet_area_total+=AreaWithinTPC(fJetR,myJet->Eta());
}
}
-
+ delete jet_vec;
+ delete track_vec;
+
// Calculate Rho
TPC_rho = E_tracks_total/(fTPCArea-jet_area_total);
TPC_rho*=fScaleFactor;
fRhoChargedScale->FillDeltaPtNColl(fEventCentrality,TPC_rho,fJetR,fEMCalRCBckgFlucNColl,1);
fRhoChargedScale->FillBackgroundFluctuations(fEventCentrality,TPC_rho,fJetR);
fRhoChargedScale->FillLeadingJetPtRho(fEMCalFullJet->GetLeadingPt(),TPC_rho);
- fRhoChargedScale->FillMiscJetStats(fmyAKTFullJets,fEMCalFullJet->GetJets(),fEMCalFullJet->GetTotalJets(),fOrgTracks,fOrgClusters);
-
+ fRhoChargedScale->FillMiscJetStats(fmyAKTFullJets,fEMCalFullJet->GetJets(),fEMCalFullJet->GetTotalJets(),fOrgTracks,fOrgClusters,fVertex);
}
-/*
+
void AliAnalysisTaskFullpAJets::EstimateChargedRhokT()
{
Int_t i;
delete [] RhoArray;
delete [] pTArray;
}
-*/
+
void AliAnalysisTaskFullpAJets::EstimateChargedRhoCMSScale()
{
Int_t i,k;
}
else
{
- //CMSCorrectionFactor = CMSTrackArea/CMSTotalkTArea;
CMSCorrectionFactor = CMSTrackArea/(fTPCPhiTotal*(fTPCEtaTotal-2*fJetR)); // The total physical area should be reduced by the eta cut due to looping over only fully contained kT jets within the TPC
}
kTRho*=CMSCorrectionFactor;
fRhoChargedCMSScale->FillBackgroundFluctuations(fEventCentrality,kTRho,fJetR);
fRhoChargedCMSScale->FillLeadingJetPtRho(fEMCalFullJet->GetLeadingPt(),kTRho);
fRhoChargedCMSScale->DoNEFAnalysis(fNEFSignalJetCut,fEMCalJetThreshold,fmyAKTFullJets,fEMCalFullJet->GetJets(),fEMCalFullJet->GetTotalJets(),fmyClusters,fOrgClusters,fEvent,fEMCALGeometry,fRecoUtil,fCells);
- fRhoChargedCMSScale->FillMiscJetStats(fmyAKTFullJets,fEMCalFullJet->GetJets(),fEMCalFullJet->GetTotalJets(),fOrgTracks,fOrgClusters);
+ fRhoChargedCMSScale->FillMiscJetStats(fmyAKTFullJets,fEMCalFullJet->GetJets(),fEMCalFullJet->GetTotalJets(),fOrgTracks,fOrgClusters,fVertex);
+ fRhoChargedCMSScale->FillJetEventCentrality(fEMCalFullJet->GetLeadingPt(),fEvent);
delete [] RhoArray;
delete [] pTArray;
}
-/*
+
// Full Rho's
void AliAnalysisTaskFullpAJets::EstimateFullRho0()
{
Double_t E_caloclusters_total=0.0;
Double_t EMCal_rho=0.0;
+ TLorentzVector *cluster_vec = new TLorentzVector;
+
// Loop over all tracks
for (i=0;i<fnTracks;i++)
{
for (i=0;i<fnClusters;i++)
{
AliVCluster* vcluster = (AliVCluster*) fmyClusters->At(i);
- TLorentzVector *cluster_vec = new TLorentzVector;
vcluster->GetMomentum(*cluster_vec,fVertex);
E_caloclusters_total+=cluster_vec->Pt();
- //E_caloclusters_total+=0.5*cluster_vec->Pt();
}
-
+ delete cluster_vec;
+
// Calculate the mean Background density
EMCal_rho=(E_tracks_total+E_caloclusters_total)/fEMCalArea;
fRhoFull=EMCal_rho;
// Fill Histograms
- if (fRhoCharged>0)
- {
- fpRhoScale->Fill(fEventCentrality,fRhoFull/fRhoCharged);
- fhRhoScale->Fill(fRhoFull/fRhoCharged,fEventCentrality);
- }
-
fRhoFull0->FillRho(fEventCentrality,EMCal_rho);
fRhoFull0->FillBSJS(fEventCentrality,EMCal_rho,fEMCalJetThreshold,fmyAKTFullJets,fEMCalFullJet->GetJets(),fEMCalFullJet->GetTotalJets());
fRhoFull0->FillDeltaPt(fEventCentrality,EMCal_rho,fJetR,fEMCalRCBckgFluc,1);
Double_t E_caloclusters_total=0.0;
Double_t EMCal_rho=0.0;
+ TLorentzVector *temp_jet= new TLorentzVector;
+ TLorentzVector *track_vec = new TLorentzVector;
+ TLorentzVector *cluster_vec = new TLorentzVector;
+
if (fEMCalPartJetUnbiased->GetLeadingPt()>=fEMCalJetThreshold)
{
AliEmcalJet *myJet =(AliEmcalJet*) fmyAKTFullJets->At(fEMCalPartJetUnbiased->GetLeadingIndex());
- TLorentzVector *temp_jet= new TLorentzVector;
myJet->GetMom(*temp_jet);
// Loop over all tracks
AliVTrack* vtrack =(AliVTrack*) fmyTracks->At(i);
if (IsInEMCal(vtrack->Phi(),vtrack->Eta())==kTRUE)
{
- TLorentzVector *track_vec = new TLorentzVector;
track_vec->SetPtEtaPhiE(vtrack->Pt(),vtrack->Eta(),vtrack->Phi(),vtrack->E());
if (temp_jet->DeltaR(*track_vec)>fJetRForRho)
{
E_tracks_total+=vtrack->Pt();
}
- delete track_vec;
}
}
for (i=0;i<fnClusters;i++)
{
AliVCluster* vcluster = (AliVCluster*) fmyClusters->At(i);
- TLorentzVector *cluster_vec = new TLorentzVector;
vcluster->GetMomentum(*cluster_vec,fVertex);
if (temp_jet->DeltaR(*cluster_vec)>fJetRForRho)
{
- E_caloclusters_total+=vcluster->E();
+ E_caloclusters_total+=cluster_vec->Pt();
}
- delete cluster_vec;
}
- delete temp_jet;
// Calculate the mean Background density
EMCal_rho=(E_tracks_total+E_caloclusters_total)/(fEMCalArea-AreaWithinEMCal(fJetR,myJet->Phi(),myJet->Eta()));
}
for (i=0;i<fnClusters;i++)
{
AliVCluster* vcluster = (AliVCluster*) fmyClusters->At(i);
- E_caloclusters_total+=vcluster->E();
+ vcluster->GetMomentum(*cluster_vec,fVertex);
+ E_caloclusters_total+=cluster_vec->Pt();
}
// Calculate the mean Background density
EMCal_rho=(E_tracks_total+E_caloclusters_total)/fEMCalArea;
}
-
+ delete temp_jet;
+ delete track_vec;
+ delete cluster_vec;
+
// Fill histograms
fRhoFull1->FillRho(fEventCentrality,EMCal_rho);
fRhoFull1->FillBSJS(fEventCentrality,EMCal_rho,fEMCalJetThreshold,fmyAKTFullJets,fEMCalFullJet->GetJets(),fEMCalFullJet->GetTotalJets());
Double_t E_caloclusters_total=0.0;
Double_t EMCal_rho=0.0;
+ TLorentzVector *temp_jet1 = new TLorentzVector;
+ TLorentzVector *temp_jet2 = new TLorentzVector;
+ TLorentzVector *track_vec = new TLorentzVector;
+ TLorentzVector *cluster_vec = new TLorentzVector;
+
if ((fEMCalPartJetUnbiased->GetLeadingPt()>=fEMCalJetThreshold) && (fEMCalPartJetUnbiased->GetSubLeadingPt()>=fEMCalJetThreshold))
{
AliEmcalJet *myhJet =(AliEmcalJet*) fmyAKTFullJets->At(fEMCalPartJetUnbiased->GetLeadingIndex());
- TLorentzVector *temp_jet1 = new TLorentzVector;
myhJet->GetMom(*temp_jet1);
AliEmcalJet *myJet =(AliEmcalJet*) fmyAKTFullJets->At(fEMCalPartJetUnbiased->GetSubLeadingIndex());
- TLorentzVector *temp_jet2 = new TLorentzVector;
myJet->GetMom(*temp_jet2);
// Loop over all tracks
AliVTrack* vtrack =(AliVTrack*) fmyTracks->At(i);
if (IsInEMCal(vtrack->Phi(),vtrack->Eta())==kTRUE)
{
- TLorentzVector *track_vec = new TLorentzVector;
track_vec->SetPtEtaPhiE(vtrack->Pt(),vtrack->Eta(),vtrack->Phi(),vtrack->E());
if ((temp_jet1->DeltaR(*track_vec)>fJetRForRho) && (temp_jet2->DeltaR(*track_vec)>fJetRForRho))
{
E_tracks_total+=vtrack->Pt();
}
- delete track_vec;
}
}
for (i=0;i<fnClusters;i++)
{
AliVCluster* vcluster = (AliVCluster*) fmyClusters->At(i);
- TLorentzVector *cluster_vec = new TLorentzVector;
vcluster->GetMomentum(*cluster_vec,fVertex);
if ((temp_jet1->DeltaR(*cluster_vec)>fJetRForRho) && (temp_jet2->DeltaR(*cluster_vec)>fJetRForRho))
{
- E_caloclusters_total+=vcluster->E();
+ E_caloclusters_total+=cluster_vec->Pt();
}
- delete cluster_vec;
}
- delete temp_jet1;
- delete temp_jet2;
-
+
// Calculate the mean Background density
EMCal_rho=(E_tracks_total+E_caloclusters_total)/(fEMCalArea-AreaWithinEMCal(fJetR,myhJet->Phi(),myhJet->Eta())-AreaWithinEMCal(fJetR,myJet->Phi(),myJet->Eta()));
}
else if (fEMCalPartJetUnbiased->GetLeadingPt()>=fEMCalJetThreshold)
{
AliEmcalJet *myJet =(AliEmcalJet*) fmyAKTFullJets->At(fEMCalPartJetUnbiased->GetLeadingIndex());
- TLorentzVector *temp_jet= new TLorentzVector;
- myJet->GetMom(*temp_jet);
+ myJet->GetMom(*temp_jet1);
// Loop over all tracks
for (i=0;i<fnTracks;i++)
AliVTrack* vtrack =(AliVTrack*) fmyTracks->At(i);
if (IsInEMCal(vtrack->Phi(),vtrack->Eta())==kTRUE)
{
- TLorentzVector *track_vec = new TLorentzVector;
track_vec->SetPtEtaPhiE(vtrack->Pt(),vtrack->Eta(),vtrack->Phi(),vtrack->E());
- if (temp_jet->DeltaR(*track_vec)>fJetRForRho)
+ if (temp_jet1->DeltaR(*track_vec)>fJetRForRho)
{
E_tracks_total+=vtrack->Pt();
}
- delete track_vec;
}
}
for (i=0;i<fnClusters;i++)
{
AliVCluster* vcluster = (AliVCluster*) fmyClusters->At(i);
- TLorentzVector *cluster_vec = new TLorentzVector;
vcluster->GetMomentum(*cluster_vec,fVertex);
- if (temp_jet->DeltaR(*cluster_vec)>fJetRForRho)
+ if (temp_jet1->DeltaR(*cluster_vec)>fJetRForRho)
{
- E_caloclusters_total+=vcluster->E();
+ E_caloclusters_total+=cluster_vec->Pt();
}
- delete cluster_vec;
}
- delete temp_jet;
// Calculate the mean Background density
EMCal_rho=(E_tracks_total+E_caloclusters_total)/(fEMCalArea-AreaWithinEMCal(fJetR,myJet->Phi(),myJet->Eta()));
}
for (i=0;i<fnClusters;i++)
{
AliVCluster* vcluster = (AliVCluster*) fmyClusters->At(i);
- E_caloclusters_total+=vcluster->E();
+ vcluster->GetMomentum(*cluster_vec,fVertex);
+ E_caloclusters_total+=cluster_vec->Pt();
}
// Calculate the mean Background density
EMCal_rho=(E_tracks_total+E_caloclusters_total)/fEMCalArea;
}
-
+ delete temp_jet1;
+ delete temp_jet2;
+ delete track_vec;
+ delete cluster_vec;
+
// Fill histograms
fRhoFull2->FillRho(fEventCentrality,EMCal_rho);
fRhoFull2->FillBSJS(fEventCentrality,EMCal_rho,fEMCalJetThreshold,fmyAKTFullJets,fEMCalFullJet->GetJets(),fEMCalFullJet->GetTotalJets());
Double_t EMCal_rho=0.0;
Double_t jet_area_total=0.0;
+ TLorentzVector *jet_vec= new TLorentzVector;
+ TLorentzVector *track_vec = new TLorentzVector;
+ TLorentzVector *cluster_vec = new TLorentzVector;
+
// First, sum all tracks within the EMCal that are away from jet(s) above Pt Threshold
for (i=0;i<fnTracks;i++)
{
{
track_away_from_jet=kTRUE;
j=0;
- TLorentzVector *track_vec = new TLorentzVector;
track_vec->SetPtEtaPhiE(vtrack->Pt(),vtrack->Eta(),vtrack->Phi(),vtrack->E());
while (track_away_from_jet==kTRUE && j<fEMCalPartJetUnbiased->GetTotalSignalJets())
{
- TLorentzVector *jet_vec= new TLorentzVector;
AliEmcalJet *myJet =(AliEmcalJet*) fmyAKTFullJets->At(fEMCalPartJetUnbiased->GetSignalJetIndex(j));
myJet->GetMom(*jet_vec);
if (track_vec->DeltaR(*jet_vec)<=fJetRForRho)
{
track_away_from_jet=kFALSE;
}
- delete jet_vec;
j++;
}
if (track_away_from_jet==kTRUE)
{
E_tracks_total+=vtrack->Pt();
}
- delete track_vec;
}
}
}
for (i=0;i<fnClusters;i++)
{
AliVCluster* vcluster = (AliVCluster*) fmyClusters->At(i);
+ vcluster->GetMomentum(*cluster_vec,fVertex);
if (fEMCalPartJet->GetTotalSignalJets()<1)
{
- E_caloclusters_total+=vcluster->E();
+ E_caloclusters_total+=cluster_vec->Pt();
}
else
{
cluster_away_from_jet=kTRUE;
j=0;
- TLorentzVector *cluster_vec = new TLorentzVector;
- vcluster->GetMomentum(*cluster_vec,fVertex);
while (cluster_away_from_jet==kTRUE && j<fEMCalPartJetUnbiased->GetTotalSignalJets())
{
- TLorentzVector *jet_vec= new TLorentzVector;
AliEmcalJet *myJet =(AliEmcalJet*) fmyAKTFullJets->At(fEMCalPartJetUnbiased->GetSignalJetIndex(j));
myJet->GetMom(*jet_vec);
if (cluster_vec->DeltaR(*jet_vec)<=fJetRForRho)
{
cluster_away_from_jet=kFALSE;
}
- delete jet_vec;
j++;
}
if (cluster_away_from_jet==kTRUE)
{
- E_caloclusters_total+=vcluster->E();
+ E_caloclusters_total+=cluster_vec->Pt();
}
- delete cluster_vec;
}
}
jet_area_total+=AreaWithinEMCal(fJetR,myJet->Phi(),myJet->Eta());
}
}
-
+ delete jet_vec;
+ delete track_vec;
+ delete cluster_vec;
+
// Calculate Rho
EMCal_rho=(E_tracks_total+E_caloclusters_total)/(fEMCalArea-jet_area_total);
Double_t E_caloclusters_total=0.0;
Double_t EMCal_rho=0.0;
+ TLorentzVector *cluster_vec = new TLorentzVector;
+
// Loop over all tracks
for (i=0;i<fnTracks;i++)
{
for (i=0;i<fnClusters;i++)
{
AliVCluster* vcluster = (AliVCluster*) fmyClusters->At(i);
- E_caloclusters_total+=vcluster->E();
+ vcluster->GetMomentum(*cluster_vec,fVertex);
+ E_caloclusters_total+=cluster_vec->Pt();
}
+ delete cluster_vec;
+
// Calculate the mean Background density
EMCal_rho=(E_tracks_total+E_caloclusters_total)/fEMCalArea;
delete [] RhoArray;
delete [] pTArray;
}
-*/
-void AliAnalysisTaskFullpAJets::DeleteJetData(Bool_t EMCalOn)
+
+void AliAnalysisTaskFullpAJets::FullJetEnergyDensityProfile()
+{
+ Int_t i,j;
+ Double_t delta_R;
+ Double_t dR=0.0;
+ Double_t fullEDR = fFullEDJetR*100;
+ const Int_t fullEDBins = (Int_t) fullEDR;
+ Double_t ED_pT[fullEDBins];
+ dR = fFullEDJetR/Double_t(fullEDBins); // Should be 0.01 be default
+
+ TLorentzVector *jet_vec= new TLorentzVector;
+ TLorentzVector *track_vec = new TLorentzVector;
+ TLorentzVector *cluster_vec = new TLorentzVector;
+
+ for (i=0;i<fEMCalFullJet->GetTotalSignalJets();i++)
+ {
+ AliEmcalJet *myJet = (AliEmcalJet*) fmyAKTFullJets->At(fEMCalFullJet->GetSignalJetIndex(i));
+
+ if (IsInEMCalFull(fFullEDJetR,myJet->Phi(),myJet->Eta())==kTRUE)
+ {
+ for (j=0;j<fullEDBins;j++)
+ {
+ ED_pT[j]=0;
+ }
+ myJet->GetMom(*jet_vec);
+
+ // Sum all tracks in concentric rings around jet vertex
+ for (j=0;j<fnTracks;j++)
+ {
+ AliVTrack* vtrack = (AliVTrack*) fmyTracks->At(j);
+ track_vec->SetPtEtaPhiE(vtrack->Pt(),vtrack->Eta(),vtrack->Phi(),vtrack->E());
+ delta_R=jet_vec->DeltaR(*track_vec);
+ if (delta_R<fFullEDJetR)
+ {
+ ED_pT[TMath::FloorNint((delta_R/fFullEDJetR)*fullEDBins)]+=vtrack->Pt();
+ }
+ }
+
+ // Sum all clusters in concentric rings around jet vertex
+ for (j=0;j<fnClusters;j++)
+ {
+ AliVCluster* vcluster = (AliVCluster*) fmyClusters->At(j);
+ vcluster->GetMomentum(*cluster_vec,fVertex);
+ delta_R=jet_vec->DeltaR(*cluster_vec);
+ if (delta_R<fFullEDJetR)
+ {
+ ED_pT[TMath::FloorNint((delta_R/fFullEDJetR)*fullEDBins)]+=cluster_vec->Pt();
+ }
+ }
+
+ for (j=0;j<fullEDBins;j++)
+ {
+ ED_pT[j] /= TMath::Pi()*dR*dR*(2*j+1);
+ fpFullJetEDProfile->Fill(myJet->Pt(),fEventCentrality,j*dR,ED_pT[j]);
+ }
+ }
+ }
+ delete cluster_vec;
+ delete track_vec;
+ delete jet_vec;
+}
+
+// Although this calculates the charged pT density in concentric rings around a charged jet contained fudically within the TPC (R=0.8), we actually scale the density by the scale factor (fScaleFactor) in order to compare to 'Full Charge Density' and 'Rho Median Occupancy Approach' (RhoChargedCMSScaled). We still keep the charged value for completeness.
+void AliAnalysisTaskFullpAJets::ChargedJetEnergyDensityProfile()
+{
+ Int_t i,j;
+ Double_t delta_R;
+ Double_t dR=0.0;
+ Double_t chargedEDR = fChargedEDJetR*100;
+ const Int_t chargedEDBins = (Int_t) chargedEDR;
+ Double_t ED_pT[chargedEDBins];
+ dR = fChargedEDJetR/Double_t(chargedEDBins); // Should be 0.01 be default
+
+ TLorentzVector *jet_vec= new TLorentzVector;
+ TLorentzVector *track_vec = new TLorentzVector;
+
+ for (i=0;i<fTPCFullJet->GetTotalSignalJets();i++)
+ {
+ AliEmcalJet *myJet = (AliEmcalJet*) fmyAKTChargedJets->At(fTPCFullJet->GetSignalJetIndex(i));
+ if (IsInTPC(fChargedEDJetR,myJet->Phi(),myJet->Eta(),kTRUE)==kTRUE)
+ {
+ for (j=0;j<chargedEDBins;j++)
+ {
+ ED_pT[j]=0;
+ }
+ myJet->GetMom(*jet_vec);
+
+ // Sum all tracks in concentric rings around jet vertex
+ for (j=0;j<fnTracks;j++)
+ {
+ AliVTrack* vtrack = (AliVTrack*) fmyTracks->At(j);
+ track_vec->SetPtEtaPhiE(vtrack->Pt(),vtrack->Eta(),vtrack->Phi(),vtrack->E());
+ delta_R=jet_vec->DeltaR(*track_vec);
+ if (delta_R<fChargedEDJetR)
+ {
+ ED_pT[TMath::FloorNint((delta_R/fChargedEDJetR)*chargedEDBins)]+=vtrack->Pt();
+ }
+ }
+ for (j=0;j<chargedEDBins;j++)
+ {
+ ED_pT[j] /= TMath::Pi()*dR*dR*(2*j+1);
+ fpChargedJetEDProfile->Fill(myJet->Pt(),fEventCentrality,j*dR,ED_pT[j]);
+ fpChargedJetEDProfileScaled->Fill(myJet->Pt(),fEventCentrality,j*dR,fScaleFactor*ED_pT[j]);
+ }
+ }
+ }
+ delete track_vec;
+ delete jet_vec;
+}
+
+void AliAnalysisTaskFullpAJets::DeleteJetData(Int_t delOption)
{
- delete fmyTracks;
- delete fTPCJet;
- delete fTPCFullJet;
- delete fTPCOnlyJet;
- delete fTPCkTFullJet;
- if (EMCalOn==kTRUE)
+ if (delOption ==0)
{
+ delete fRecoUtil;
+
+ fRecoUtil = NULL;
+ }
+ else if (delOption==1)
+ {
+ delete fRecoUtil;
+ delete fmyTracks;
+ delete fmyClusters;
+
+ fRecoUtil = NULL;
+ fmyTracks = NULL;
+ fmyClusters = NULL;
+ }
+ else if (delOption==2)
+ {
+ delete fRecoUtil;
+ delete fmyTracks;
delete fmyClusters;
+ delete fTPCJet;
+ delete fTPCFullJet;
+ delete fTPCOnlyJet;
+ delete fTPCJetUnbiased;
+ delete fTPCkTFullJet;
+
+ fRecoUtil = NULL;
+ fmyTracks = NULL;
+ fmyClusters = NULL;
+ fTPCJet = NULL;
+ fTPCFullJet = NULL;
+ fTPCOnlyJet = NULL;
+ fTPCJetUnbiased = NULL;
+ fTPCkTFullJet = NULL;
+ }
+ if (delOption==3)
+ {
+ delete fRecoUtil;
+ delete fmyTracks;
+ delete fmyClusters;
+ delete fTPCJet;
+ delete fTPCFullJet;
+ delete fTPCOnlyJet;
+ delete fTPCJetUnbiased;
+ delete fTPCkTFullJet;
delete fEMCalJet;
delete fEMCalFullJet;
delete fEMCalPartJet;
+ delete fEMCalPartJetUnbiased;
delete fEMCalkTFullJet;
+
+ fRecoUtil = NULL;
+ fmyTracks = NULL;
+ fmyClusters = NULL;
+ fTPCJet = NULL;
+ fTPCFullJet = NULL;
+ fTPCOnlyJet = NULL;
+ fTPCJetUnbiased = NULL;
+ fEMCalJet = NULL;
+ fEMCalFullJet = NULL;
+ fEMCalPartJet = NULL;
+ fEMCalPartJetUnbiased = NULL;
+ fEMCalkTFullJet = NULL;
}
}
fJetR(0),
fSignalPt(0),
fAreaCutFrac(0.6),
- fNEF(0.9),
- fSignalTrackBias(1),
+ fNEF(1.0),
+ fSignalTrackBias(0),
fPtMaxIndex(0),
fPtMax(0),
fPtSubLeadingIndex(0),
fJetR(0),
fSignalPt(0),
fAreaCutFrac(0.6),
- fNEF(0.9),
- fSignalTrackBias(1),
+ fNEF(1.0),
+ fSignalTrackBias(0),
fPtMaxIndex(0),
fPtMax(0),
fPtSubLeadingIndex(0),
SetSignalCut(0);
SetAreaCutFraction(0.6);
SetJetR(fJetR);
- SetSignalTrackPtBias(1);
+ SetSignalTrackPtBias(0);
}
// Destructor
AliAnalysisTaskFullpAJets::AlipAJetData::~AlipAJetData()
{
- if (fnTotal!=0)
- {
- SetName("");
- SetIsJetsFull(kFALSE);
- SetTotalEntries(0);
- SetTotalJets(0);
- SetTotalSignalJets(0);
- SetLeading(0,0);
- SetSubLeading(0,0);
- SetSignalCut(0);
- SetAreaCutFraction(0);
- SetJetR(0);
- SetNEF(0);
- SetSignalTrackPtBias(kTRUE);
-
- delete [] fJetsIndex;
- delete [] fJetsSCIndex;
- delete [] fIsJetInArray;
- delete [] fJetMaxChargedPt;
- }
+ SetName("");
+ SetIsJetsFull(kFALSE);
+ SetTotalJets(0);
+ SetTotalSignalJets(0);
+ SetLeading(0,0);
+ SetSubLeading(0,0);
+ SetSignalCut(0);
+ SetAreaCutFraction(0);
+ SetJetR(0);
+ SetNEF(0);
+ SetSignalTrackPtBias(kFALSE);
+
+ delete [] fJetsIndex;
+ delete [] fJetsSCIndex;
+ delete [] fIsJetInArray;
+ delete [] fJetMaxChargedPt;
+
+ fnTotal = 0;
}
// User Defined Sub-Routines
fh80100BSPt(0),
fhBSPt(0),
fhBSPtCen(0),
-/* fhBSPtCenLCT(0),*/
fh020BSPtSignal(0),
fh80100BSPtSignal(0),
fhBSPtSignal(0),
fpRho(0),
fpLJetRho(0),
- fhJetConstituentPt(0),
+
+ fhJetPtEtaPhi(0),
fhJetPtArea(0),
+ fhJetConstituentPt(0),
+ fhJetTracksPt(0),
+ fhJetClustersPt(0),
+ fhJetConstituentCounts(0),
+ fhJetTracksCounts(0),
+ fhJetClustersCounts(0),
+ fhJetPtZConstituent(0),
+ fhJetPtZTrack(0),
+ fhJetPtZCluster(0),
+ fhJetPtZLeadingConstituent(0),
+ fhJetPtZLeadingTrack(0),
+ fhJetPtZLeadingCluster(0),
+
+ fhEventCentralityVsZNA(0),
+ fhEventCentralityVsZNAPt(0),
fNEFOutput(0),
- fhNEF(0),
- fhNEFSignal(0),
- fhNEFJetPt(0),
- fhNEFJetPtSignal(0),
- fhNEFEtaPhi(0),
- fhNEFEtaPhiSignal(0),
- fhEtaPhiNEF(0),
- fhNEFTotalMult(0),
- fhNEFTotalMultSignal(0),
- fhNEFNeutralMult(0),
- fhNEFNeutralMultSignal(0),
+ fhJetPtNEF(0),
+ fhJetNEFInfo(0),
+ fhJetNEFSignalInfo(0),
+ fhClusterNEFInfo(0),
+ fhClusterNEFSignalInfo(0),
fhClusterShapeAll(0),
fhClusterPtCellAll(0),
- fhNEFJetPtFCross(0),
- fhNEFZLeadingFCross(0),
- fhNEFTimeCellCount(0),
- fhNEFTimeDeltaTime(0),
fName(0),
fCentralityTag(0),
fLChargedTrackPtLow(0),
fLChargedTrackPtUp(0),
fDoNEFQAPlots(0),
+ fDoNEFSignalOnly(1),
+ fSignalTrackBias(0),
+ fDoTHnSparse(0),
fNEFBins(0),
fNEFLow(0),
fNEFUp(0),
+ fnDimJet(0),
+ fnDimCluster(0),
fEMCalPhiMin(1.39626),
fEMCalPhiMax(3.26377),
fEMCalEtaMin(-0.7),
fh80100BSPt(0),
fhBSPt(0),
fhBSPtCen(0),
-/* fhBSPtCenLCT(0),*/
fh020BSPtSignal(0),
fh80100BSPtSignal(0),
fhBSPtSignal(0),
fpRho(0),
fpLJetRho(0),
- fhJetConstituentPt(0),
+
+ fhJetPtEtaPhi(0),
fhJetPtArea(0),
+ fhJetConstituentPt(0),
+ fhJetTracksPt(0),
+ fhJetClustersPt(0),
+ fhJetConstituentCounts(0),
+ fhJetTracksCounts(0),
+ fhJetClustersCounts(0),
+ fhJetPtZConstituent(0),
+ fhJetPtZTrack(0),
+ fhJetPtZCluster(0),
+ fhJetPtZLeadingConstituent(0),
+ fhJetPtZLeadingTrack(0),
+ fhJetPtZLeadingCluster(0),
+
+ fhEventCentralityVsZNA(0),
+ fhEventCentralityVsZNAPt(0),
fNEFOutput(0),
- fhNEF(0),
- fhNEFSignal(0),
- fhNEFJetPt(0),
- fhNEFJetPtSignal(0),
- fhNEFEtaPhi(0),
- fhNEFEtaPhiSignal(0),
- fhEtaPhiNEF(0),
- fhNEFTotalMult(0),
- fhNEFTotalMultSignal(0),
- fhNEFNeutralMult(0),
- fhNEFNeutralMultSignal(0),
+ fhJetPtNEF(0),
+ fhJetNEFInfo(0),
+ fhJetNEFSignalInfo(0),
+ fhClusterNEFInfo(0),
+ fhClusterNEFSignalInfo(0),
fhClusterShapeAll(0),
fhClusterPtCellAll(0),
- fhNEFJetPtFCross(0),
- fhNEFZLeadingFCross(0),
- fhNEFTimeCellCount(0),
- fhNEFTimeDeltaTime(0),
fName(0),
fCentralityTag(0),
fLChargedTrackPtLow(0),
fLChargedTrackPtUp(0),
fDoNEFQAPlots(0),
+ fDoNEFSignalOnly(1),
+ fSignalTrackBias(0),
+ fDoTHnSparse(0),
fNEFBins(0),
fNEFLow(0),
fNEFUp(0),
+ fnDimJet(0),
+ fnDimCluster(0),
fEMCalPhiMin(1.39626),
fEMCalPhiMax(3.26377),
fEMCalEtaMin(-0.7),
{
SetName(name);
- SetCentralityTag("V0A");
+ SetCentralityTag("ZNA");
SetCentralityRange(100,0,100);
SetPtRange(250,-50,200);
SetRhoPtRange(500,0,50);
Init();
}
-AliAnalysisTaskFullpAJets::AlipAJetHistos::AlipAJetHistos(const char *name, const char *centag, Bool_t doNEF) :
+AliAnalysisTaskFullpAJets::AlipAJetHistos::AlipAJetHistos(const char *name, TString centag, Bool_t doNEF) :
fOutput(0),
fh80100BSPt(0),
fhBSPt(0),
fhBSPtCen(0),
-/* fhBSPtCenLCT(0),*/
fh020BSPtSignal(0),
fh80100BSPtSignal(0),
fhBSPtSignal(0),
fpRho(0),
fpLJetRho(0),
- fhJetConstituentPt(0),
+
+ fhJetPtEtaPhi(0),
fhJetPtArea(0),
+ fhJetConstituentPt(0),
+ fhJetTracksPt(0),
+ fhJetClustersPt(0),
+ fhJetConstituentCounts(0),
+ fhJetTracksCounts(0),
+ fhJetClustersCounts(0),
+ fhJetPtZConstituent(0),
+ fhJetPtZTrack(0),
+ fhJetPtZCluster(0),
+ fhJetPtZLeadingConstituent(0),
+ fhJetPtZLeadingTrack(0),
+ fhJetPtZLeadingCluster(0),
+
+ fhEventCentralityVsZNA(0),
+ fhEventCentralityVsZNAPt(0),
fNEFOutput(0),
- fhNEF(0),
- fhNEFSignal(0),
- fhNEFJetPt(0),
- fhNEFJetPtSignal(0),
- fhNEFEtaPhi(0),
- fhNEFEtaPhiSignal(0),
- fhEtaPhiNEF(0),
- fhNEFTotalMult(0),
- fhNEFTotalMultSignal(0),
- fhNEFNeutralMult(0),
- fhNEFNeutralMultSignal(0),
+ fhJetPtNEF(0),
+ fhJetNEFInfo(0),
+ fhJetNEFSignalInfo(0),
+ fhClusterNEFInfo(0),
+ fhClusterNEFSignalInfo(0),
fhClusterShapeAll(0),
fhClusterPtCellAll(0),
- fhNEFJetPtFCross(0),
- fhNEFZLeadingFCross(0),
- fhNEFTimeCellCount(0),
- fhNEFTimeDeltaTime(0),
fName(0),
fCentralityTag(0),
fLChargedTrackPtLow(0),
fLChargedTrackPtUp(0),
fDoNEFQAPlots(0),
+ fDoNEFSignalOnly(1),
+ fSignalTrackBias(0),
+ fDoTHnSparse(0),
fNEFBins(0),
fNEFLow(0),
fNEFUp(0),
+ fnDimJet(0),
+ fnDimCluster(0),
fEMCalPhiMin(1.39626),
fEMCalPhiMax(3.26377),
fEMCalEtaMin(-0.7),
{
SetName(name);
- SetCentralityTag(centag);
+ SetCentralityTag(centag.Data());
SetCentralityRange(100,0,100);
SetPtRange(250,-50,200);
SetRhoPtRange(500,0,50);
Init();
}
+AliAnalysisTaskFullpAJets::AlipAJetHistos::AlipAJetHistos(const char *name, TString centag, Bool_t doNEF, Bool_t doNEFSignalOnly, Bool_t doTHnSparse) :
+
+ fOutput(0),
+
+ fh020Rho(0),
+ fh80100Rho(0),
+ fhRho(0),
+ fhRhoCen(0),
+ fh020BSPt(0),
+ fh80100BSPt(0),
+ fhBSPt(0),
+ fhBSPtCen(0),
+ fh020BSPtSignal(0),
+ fh80100BSPtSignal(0),
+ fhBSPtSignal(0),
+ fhBSPtCenSignal(0),
+ fh020DeltaPt(0),
+ fh80100DeltaPt(0),
+ fhDeltaPt(0),
+ fhDeltaPtCen(0),
+ fh020DeltaPtSignal(0),
+ fh80100DeltaPtSignal(0),
+ fhDeltaPtSignal(0),
+ fhDeltaPtCenSignal(0),
+ fh020DeltaPtNColl(0),
+ fh80100DeltaPtNColl(0),
+ fhDeltaPtNColl(0),
+ fhDeltaPtCenNColl(0),
+ fh020BckgFlucPt(0),
+ fh80100BckgFlucPt(0),
+ fhBckgFlucPt(0),
+ fhBckgFlucPtCen(0),
+
+ fpRho(0),
+ fpLJetRho(0),
+
+ fhJetPtEtaPhi(0),
+ fhJetPtArea(0),
+ fhJetConstituentPt(0),
+ fhJetTracksPt(0),
+ fhJetClustersPt(0),
+ fhJetConstituentCounts(0),
+ fhJetTracksCounts(0),
+ fhJetClustersCounts(0),
+ fhJetPtZConstituent(0),
+ fhJetPtZTrack(0),
+ fhJetPtZCluster(0),
+ fhJetPtZLeadingConstituent(0),
+ fhJetPtZLeadingTrack(0),
+ fhJetPtZLeadingCluster(0),
+
+ fhEventCentralityVsZNA(0),
+ fhEventCentralityVsZNAPt(0),
+
+ fNEFOutput(0),
+ fhJetPtNEF(0),
+ fhJetNEFInfo(0),
+ fhJetNEFSignalInfo(0),
+ fhClusterNEFInfo(0),
+ fhClusterNEFSignalInfo(0),
+ fhClusterShapeAll(0),
+ fhClusterPtCellAll(0),
+
+ fName(0),
+ fCentralityTag(0),
+ fCentralityBins(0),
+ fCentralityLow(0),
+ fCentralityUp(0),
+ fPtBins(0),
+ fPtLow(0),
+ fPtUp(0),
+ fRhoPtBins(0),
+ fRhoPtLow(0),
+ fRhoPtUp(0),
+ fDeltaPtBins(0),
+ fDeltaPtLow(0),
+ fDeltaPtUp(0),
+ fBckgFlucPtBins(0),
+ fBckgFlucPtLow(0),
+ fBckgFlucPtUp(0),
+ fLJetPtBins(0),
+ fLJetPtLow(0),
+ fLJetPtUp(0),
+ fRhoValue(0),
+ fLChargedTrackPtBins(0),
+ fLChargedTrackPtLow(0),
+ fLChargedTrackPtUp(0),
+ fDoNEFQAPlots(0),
+ fDoNEFSignalOnly(1),
+ fSignalTrackBias(0),
+ fDoTHnSparse(0),
+ fNEFBins(0),
+ fNEFLow(0),
+ fNEFUp(0),
+ fnDimJet(0),
+ fnDimCluster(0),
+ fEMCalPhiMin(1.39626),
+ fEMCalPhiMax(3.26377),
+ fEMCalEtaMin(-0.7),
+ fEMCalEtaMax(0.7)
+
+{
+ SetName(name);
+ SetCentralityTag(centag.Data());
+ SetCentralityRange(100,0,100);
+ SetPtRange(250,-50,200);
+ SetRhoPtRange(500,0,50);
+ SetDeltaPtRange(200,-100,100);
+ SetBackgroundFluctuationsPtRange(100,0,100);
+ SetLeadingJetPtRange(200,0,200);
+ SetLeadingChargedTrackPtRange(100,0,100);
+ SetNEFRange(100,0,1);
+ DoNEFQAPlots(doNEF);
+ DoNEFSignalOnly(doNEFSignalOnly);
+ DoTHnSparse(doTHnSparse);
+
+ Init();
+}
+
// Destructor
AliAnalysisTaskFullpAJets::AlipAJetHistos::~AlipAJetHistos()
{
void AliAnalysisTaskFullpAJets::AlipAJetHistos::Init()
{
+ Int_t i;
+
// Initialize Private Variables
+ Int_t TCBins = 100;
fEMCalPhiMin=(80/(double)360)*2*TMath::Pi();
fEMCalPhiMax=(187/(double)360)*2*TMath::Pi();
fEMCalEtaMin=-0.7;
TString DeltaPtString="";
TString BckgFlucPtString="";
TString CentralityString;
- CentralityString = Form("Centrality (%s)",fCentralityTag);
+ TString EventCentralityString;
+ CentralityString = Form("Centrality (%s) ",fCentralityTag.Data());
+ EventCentralityString = Form("%s vs ZNA Centrality ",fCentralityTag.Data());
// Rho Spectral Plots
RhoString = Form("%d-%d Centrality, Rho Spectrum",0,20);
- fh020Rho = new TH1D("fh020Rho",RhoString,fRhoPtBins,fRhoPtLow,fRhoPtUp);
+ fh020Rho = new TH1F("fh020Rho",RhoString,fRhoPtBins,fRhoPtLow,fRhoPtUp);
fh020Rho->GetXaxis()->SetTitle("p_{T}/Area (GeV/c)");
fh020Rho->GetYaxis()->SetTitle("1/N_{Events} dN/d#rho");
fh020Rho->Sumw2();
RhoString = Form("%d-%d Centrality, Rho Spectrum",80,100);
- fh80100Rho = new TH1D("fh80100Rho",RhoString,fRhoPtBins,fRhoPtLow,fRhoPtUp);
+ fh80100Rho = new TH1F("fh80100Rho",RhoString,fRhoPtBins,fRhoPtLow,fRhoPtUp);
fh80100Rho->GetXaxis()->SetTitle("p_{T}/Area (GeV/c)");
fh80100Rho->GetYaxis()->SetTitle("1/N_{Events} dN/d#rho");
fh80100Rho->Sumw2();
RhoString = Form("%d-%d Centrality, Rho Spectrum",0,100);
- fhRho = new TH1D("fhRho",RhoString,fRhoPtBins,fRhoPtLow,fRhoPtUp);
+ fhRho = new TH1F("fhRho",RhoString,fRhoPtBins,fRhoPtLow,fRhoPtUp);
fhRho->GetXaxis()->SetTitle("p_{T}/Area (GeV/c)");
fhRho->GetYaxis()->SetTitle("1/N_{Events} dN/d#rho");
fhRho->Sumw2();
RhoString = "Rho Spectrum vs Centrality";
- fhRhoCen = new TH2D("fhRhoCen",RhoString,fRhoPtBins,fRhoPtLow,fRhoPtUp,fCentralityBins,fCentralityLow,fCentralityUp);
+ fhRhoCen = new TH2F("fhRhoCen",RhoString,fRhoPtBins,fRhoPtLow,fRhoPtUp,fCentralityBins,fCentralityLow,fCentralityUp);
fhRhoCen->GetXaxis()->SetTitle("p_{T}/Area (GeV/c)");
fhRhoCen->GetYaxis()->SetTitle(Form("%s",CentralityString.Data()));
fhRhoCen->GetZaxis()->SetTitle("1/N_{Events} dN/d#rho");
// Background Subtracted Plots
PtString = Form("%d-%d Centrality, Background Subtracted Jet Spectrum",0,20);
- fh020BSPt = new TH1D("fh020BSPt",PtString,fPtBins,fPtLow,fPtUp);
+ fh020BSPt = new TH1F("fh020BSPt",PtString,fPtBins,fPtLow,fPtUp);
fh020BSPt->GetXaxis()->SetTitle("p_{T} - #rhoA (GeV/c)");
fh020BSPt->GetYaxis()->SetTitle("1/N_{Events} dN/dp_{T}d#etad#varphi");
fh020BSPt->Sumw2();
PtString = Form("%d-%d Centrality, Background Subtracted Jet Spectrum",80,100);
- fh80100BSPt = new TH1D("fh80100BSPt",PtString,fPtBins,fPtLow,fPtUp);
+ fh80100BSPt = new TH1F("fh80100BSPt",PtString,fPtBins,fPtLow,fPtUp);
fh80100BSPt->GetXaxis()->SetTitle("p_{T} - #rhoA (GeV/c)");
fh80100BSPt->GetYaxis()->SetTitle("1/N_{Events} dN/dp_{T}d#etad#varphi");
fh80100BSPt->Sumw2();
PtString = Form("%d-%d Centrality, Background Subtracted Jet Spectrum",0,100);
- fhBSPt = new TH1D("fhBSPt",PtString,fPtBins,fPtLow,fPtUp);
+ fhBSPt = new TH1F("fhBSPt",PtString,fPtBins,fPtLow,fPtUp);
fhBSPt->GetXaxis()->SetTitle("p_{T} - #rhoA (GeV/c)");
fhBSPt->GetYaxis()->SetTitle("1/N_{Events} dN/dp_{T}d#etad#varphi");
fhBSPt->Sumw2();
PtString = "Background Subtracted Jet Spectrum vs Centrality";
- fhBSPtCen = new TH2D("fhBSPtCen",PtString,fPtBins,fPtLow,fPtUp,fCentralityBins,fCentralityLow,fCentralityUp);
+ fhBSPtCen = new TH2F("fhBSPtCen",PtString,fPtBins,fPtLow,fPtUp,fCentralityBins,fCentralityLow,fCentralityUp);
fhBSPtCen->GetXaxis()->SetTitle("p_{T} - #rhoA (GeV/c)");
fhBSPtCen->GetYaxis()->SetTitle(Form("%s",CentralityString.Data()));
fhBSPtCen->GetZaxis()->SetTitle("1/N_{Events} dN/dp_{T}d#etad#varphi");
fhBSPtCen->Sumw2();
- /*
- PtString = "Background Subtracted Jet Spectrum vs Centrality vs Leading Charge Track p_{T}";
- fhBSPtCenLCT = new TH3D("fhBSPtCenLCT",PtString,fPtBins,fPtLow,fPtUp,fCentralityBins,fCentralityLow,fCentralityUp,fLChargedTrackPtBins,fLChargedTrackPtLow,fLChargedTrackPtUp);
- fhBSPtCenLCT->GetXaxis()->SetTitle("p_{T} - #rhoA (GeV/c)");
- fhBSPtCenLCT->GetYaxis()->SetTitle(Form("%s",CentralityString.Data()));
- fhBSPtCenLCT->GetZaxis()->SetTitle("Leading Charged Track p_{T} (GeV/c)");
- fhBSPtCenLCT->Sumw2();
- */
+
PtString = Form("%d-%d Centrality, Background Subtracted Signal Jet Spectrum",0,20);
- fh020BSPtSignal = new TH1D("fh020BSPtSignal",PtString,fPtBins,fPtLow,fPtUp);
+ fh020BSPtSignal = new TH1F("fh020BSPtSignal",PtString,fPtBins,fPtLow,fPtUp);
fh020BSPtSignal->GetXaxis()->SetTitle("p_{T} - #rhoA (GeV/c)");
fh020BSPtSignal->GetYaxis()->SetTitle("1/N_{Events} dN/dp_{T}d#etad#varphi");
fh020BSPtSignal->Sumw2();
PtString = Form("%d-%d Centrality, Background Subtracted Signal Jet Spectrum",80,100);
- fh80100BSPtSignal = new TH1D("fh80100BSPtSignal",PtString,fPtBins,fPtLow,fPtUp);
+ fh80100BSPtSignal = new TH1F("fh80100BSPtSignal",PtString,fPtBins,fPtLow,fPtUp);
fh80100BSPtSignal->GetXaxis()->SetTitle("p_{T} - #rhoA (GeV/c)");
fh80100BSPtSignal->GetYaxis()->SetTitle("1/N_{Events} dN/dp_{T}d#etad#varphi");
fh80100BSPtSignal->Sumw2();
PtString = Form("%d-%d Centrality, Background Subtracted Signal Jet Spectrum",0,100);
- fhBSPtSignal = new TH1D("fhBSPtSignal",PtString,fPtBins,fPtLow,fPtUp);
+ fhBSPtSignal = new TH1F("fhBSPtSignal",PtString,fPtBins,fPtLow,fPtUp);
fhBSPtSignal->GetXaxis()->SetTitle("p_{T} - #rhoA (GeV/c)");
fhBSPtSignal->GetYaxis()->SetTitle("1/N_{Events} dN/dp_{T}d#etad#varphi");
fhBSPtSignal->Sumw2();
PtString = "Background Subtracted Signal Jet Spectrum vs Centrality";
- fhBSPtCenSignal = new TH2D("fhBSPtCenSignal",PtString,fPtBins,fPtLow,fPtUp,fCentralityBins,fCentralityLow,fCentralityUp);
+ fhBSPtCenSignal = new TH2F("fhBSPtCenSignal",PtString,fPtBins,fPtLow,fPtUp,fCentralityBins,fCentralityLow,fCentralityUp);
fhBSPtCenSignal->GetXaxis()->SetTitle("p_{T} - #rhoA (GeV/c)");
fhBSPtCenSignal->GetYaxis()->SetTitle(Form("%s",CentralityString.Data()));
fhBSPtCenSignal->GetZaxis()->SetTitle("1/N_{Events} dN/dp_{T}d#etad#varphi");
// Delta Pt Plots with RC at least 2R away from Leading Signal
DeltaPtString = Form("%d-%d Centrality, #deltap_{T} Spectrum",0,20);
- fh020DeltaPt = new TH1D("fh020DeltaPt",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp);
+ fh020DeltaPt = new TH1F("fh020DeltaPt",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp);
fh020DeltaPt->GetXaxis()->SetTitle("#deltap_{T} (GeV/c)");
fh020DeltaPt->GetYaxis()->SetTitle("Probability Density");
fh020DeltaPt->Sumw2();
DeltaPtString = Form("%d-%d Centrality, #deltap_{T} Spectrum",80,100);
- fh80100DeltaPt = new TH1D("fh80100DeltaPt",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp);
+ fh80100DeltaPt = new TH1F("fh80100DeltaPt",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp);
fh80100DeltaPt->GetXaxis()->SetTitle("#deltap_{T} (GeV/c)");
fh80100DeltaPt->GetYaxis()->SetTitle("Probability Density");
fh80100DeltaPt->Sumw2();
DeltaPtString = Form("%d-%d Centrality, #deltap_{T} Spectrum",0,100);
- fhDeltaPt = new TH1D("fhDeltaPt",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp);
+ fhDeltaPt = new TH1F("fhDeltaPt",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp);
fhDeltaPt->GetXaxis()->SetTitle("#deltap_{T} (GeV/c)");
fhDeltaPt->GetYaxis()->SetTitle("Probability Density");
fhDeltaPt->Sumw2();
DeltaPtString = "#deltap_{T} Spectrum vs Centrality";
- fhDeltaPtCen = new TH2D("fhDeltaPtCen",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp,fCentralityBins,fCentralityLow,fCentralityUp);
+ fhDeltaPtCen = new TH2F("fhDeltaPtCen",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp,fCentralityBins,fCentralityLow,fCentralityUp);
fhDeltaPtCen->GetXaxis()->SetTitle("#deltap_{T} (GeV/c)");
fhDeltaPtCen->GetYaxis()->SetTitle(Form("%s",CentralityString.Data()));
fhDeltaPtCen->GetZaxis()->SetTitle("Probability Density");
// Delta Pt Plots with no spatial restrictions on RC
DeltaPtString = Form("%d-%d Centrality, #deltap_{T} Spectrum",0,20);
- fh020DeltaPtSignal = new TH1D("fh020DeltaPtSignal",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp);
+ fh020DeltaPtSignal = new TH1F("fh020DeltaPtSignal",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp);
fh020DeltaPtSignal->GetXaxis()->SetTitle("#deltap_{T} (GeV/c)");
fh020DeltaPtSignal->GetYaxis()->SetTitle("Probability Density");
fh020DeltaPtSignal->Sumw2();
DeltaPtString = Form("%d-%d Centrality, #deltap_{T} Spectrum",80,100);
- fh80100DeltaPtSignal = new TH1D("fh80100DeltaPtSignal",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp);
+ fh80100DeltaPtSignal = new TH1F("fh80100DeltaPtSignal",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp);
fh80100DeltaPtSignal->GetXaxis()->SetTitle("#deltap_{T} (GeV/c)");
fh80100DeltaPtSignal->GetYaxis()->SetTitle("Probability Density");
fh80100DeltaPtSignal->Sumw2();
DeltaPtString = Form("%d-%d Centrality, #deltap_{T} Spectrum",0,100);
- fhDeltaPtSignal = new TH1D("fhDeltaPtSignal",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp);
+ fhDeltaPtSignal = new TH1F("fhDeltaPtSignal",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp);
fhDeltaPtSignal->GetXaxis()->SetTitle("#deltap_{T} (GeV/c)");
fhDeltaPtSignal->GetYaxis()->SetTitle("Probability Density");
fhDeltaPtSignal->Sumw2();
DeltaPtString = "#deltap_{T} Spectrum vs Centrality";
- fhDeltaPtCenSignal = new TH2D("fhDeltaPtCenSignal",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp,fCentralityBins,fCentralityLow,fCentralityUp);
+ fhDeltaPtCenSignal = new TH2F("fhDeltaPtCenSignal",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp,fCentralityBins,fCentralityLow,fCentralityUp);
fhDeltaPtCenSignal->GetXaxis()->SetTitle("#deltap_{T} (GeV/c)");
fhDeltaPtCenSignal->GetYaxis()->SetTitle(Form("%s",CentralityString.Data()));
fhDeltaPtCenSignal->GetZaxis()->SetTitle("Probability Density");
// Delta Pt Plots with NColl restrictions on RC
DeltaPtString = Form("%d-%d Centrality, #deltap_{T} Spectrum",0,20);
- fh020DeltaPtNColl = new TH1D("fh020DeltaPtNColl",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp);
+ fh020DeltaPtNColl = new TH1F("fh020DeltaPtNColl",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp);
fh020DeltaPtNColl->GetXaxis()->SetTitle("#deltap_{T} (GeV/c)");
fh020DeltaPtNColl->GetYaxis()->SetTitle("Probability Density");
fh020DeltaPtNColl->Sumw2();
DeltaPtString = Form("%d-%d Centrality, #deltap_{T} Spectrum",80,100);
- fh80100DeltaPtNColl = new TH1D("fh80100DeltaPtNColl",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp);
+ fh80100DeltaPtNColl = new TH1F("fh80100DeltaPtNColl",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp);
fh80100DeltaPtNColl->GetXaxis()->SetTitle("#deltap_{T} (GeV/c)");
fh80100DeltaPtNColl->GetYaxis()->SetTitle("Probability Density");
fh80100DeltaPtNColl->Sumw2();
DeltaPtString = Form("%d-%d Centrality, #deltap_{T} Spectrum",0,100);
- fhDeltaPtNColl = new TH1D("fhDeltaPtNColl",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp);
+ fhDeltaPtNColl = new TH1F("fhDeltaPtNColl",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp);
fhDeltaPtNColl->GetXaxis()->SetTitle("#deltap_{T} (GeV/c)");
fhDeltaPtNColl->GetYaxis()->SetTitle("Probability Density");
fhDeltaPtNColl->Sumw2();
DeltaPtString = "#deltap_{T} Spectrum vs Centrality";
- fhDeltaPtCenNColl = new TH2D("fhDeltaPtCenNColl",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp,fCentralityBins,fCentralityLow,fCentralityUp);
+ fhDeltaPtCenNColl = new TH2F("fhDeltaPtCenNColl",DeltaPtString,fDeltaPtBins,fDeltaPtLow,fDeltaPtUp,fCentralityBins,fCentralityLow,fCentralityUp);
fhDeltaPtCenNColl->GetXaxis()->SetTitle("#deltap_{T} (GeV/c)");
fhDeltaPtCenNColl->GetYaxis()->SetTitle(Form("%s",CentralityString.Data()));
fhDeltaPtCenNColl->GetZaxis()->SetTitle("Probability Density");
// Background Fluctuations Pt Plots
BckgFlucPtString = Form("%d-%d Centrality, Background Fluctuation p_{T} Spectrum",0,20);
- fh020BckgFlucPt = new TH1D("fh020BckgFlucPt",PtString,fPtBins,fPtLow,fPtUp);
+ fh020BckgFlucPt = new TH1F("fh020BckgFlucPt",PtString,fPtBins,fPtLow,fPtUp);
fh020BckgFlucPt->GetXaxis()->SetTitle("p_{T} (GeV/c)");
fh020BckgFlucPt->GetYaxis()->SetTitle("1/N_{Events} dN/dp_{T}d#etad#varphi");
fh020BckgFlucPt->Sumw2();
BckgFlucPtString = Form("%d-%d Centrality, Background Fluctuation p_{T} Spectrum",80,100);
- fh80100BckgFlucPt = new TH1D("fh80100BckgFlucPt",BckgFlucPtString,fBckgFlucPtBins,fBckgFlucPtLow,fBckgFlucPtUp);
+ fh80100BckgFlucPt = new TH1F("fh80100BckgFlucPt",BckgFlucPtString,fBckgFlucPtBins,fBckgFlucPtLow,fBckgFlucPtUp);
fh80100BckgFlucPt->GetXaxis()->SetTitle("p_{T} (GeV/c)");
fh80100BckgFlucPt->GetYaxis()->SetTitle("1/N_{Events} dN/dp_{T}d#etad#varphi");
fh80100BckgFlucPt->Sumw2();
BckgFlucPtString = Form("%d-%d Centrality, Background Fluctuation p_{T} Spectrum",0,100);
- fhBckgFlucPt = new TH1D("fhBckgFlucPt",BckgFlucPtString,fBckgFlucPtBins,fBckgFlucPtLow,fBckgFlucPtUp);
+ fhBckgFlucPt = new TH1F("fhBckgFlucPt",BckgFlucPtString,fBckgFlucPtBins,fBckgFlucPtLow,fBckgFlucPtUp);
fhBckgFlucPt->GetXaxis()->SetTitle("p_{T} (GeV/c)");
fhBckgFlucPt->GetYaxis()->SetTitle("1/N_{Events} dN/dp_{T}d#etad#varphi");
fhBckgFlucPt->Sumw2();
BckgFlucPtString = "Background Fluctuation p_{T} Spectrum vs Centrality";
- fhBckgFlucPtCen = new TH2D("fhBckgFlucPtCen",BckgFlucPtString,fBckgFlucPtBins,fBckgFlucPtLow,fBckgFlucPtUp,fCentralityBins,fCentralityLow,fCentralityUp);
+ fhBckgFlucPtCen = new TH2F("fhBckgFlucPtCen",BckgFlucPtString,fBckgFlucPtBins,fBckgFlucPtLow,fBckgFlucPtUp,fCentralityBins,fCentralityLow,fCentralityUp);
fhBckgFlucPtCen->GetXaxis()->SetTitle("#p_{T} (GeV/c)");
fhBckgFlucPtCen->GetYaxis()->SetTitle(Form("%s",CentralityString.Data()));
fhBckgFlucPtCen->GetZaxis()->SetTitle("1/N_{Events} dN/dp_{T}d#etad#varphi");
fpLJetRho->GetXaxis()->SetTitle("Leading Jet p_{T}");
fpLJetRho->GetYaxis()->SetTitle("p_{T}/Area (GeV/c)");
- // Jet pT vs Constituent pT
- fhJetConstituentPt = new TH2D("fhJetConstituentPt","Jet constituents p_{T} distribution",fPtBins,fPtLow,fPtUp,10*fPtBins,fPtLow,fPtUp);
- fhJetConstituentPt->GetXaxis()->SetTitle("Jet p_{T} (GeV/c)");
- fhJetConstituentPt->GetYaxis()->SetTitle("Constituent p_{T} (GeV/c)");
- fhJetConstituentPt->Sumw2();
-
+ // Jet pT, Eta, Phi
+ fhJetPtEtaPhi = new TH3F("fhJetPtEtaPhi","Jet p_{T} vs #eta-#varphi",fPtBins,fPtLow,fPtUp,TCBins,fEMCalEtaMin,fEMCalEtaMax,TCBins,fEMCalPhiMin,fEMCalPhiMax);
+ fhJetPtEtaPhi->GetXaxis()->SetTitle("p_{T} (GeV/c)");
+ fhJetPtEtaPhi->GetYaxis()->SetTitle("#eta");
+ fhJetPtEtaPhi->GetZaxis()->SetTitle("#varphi");
+ fhJetPtEtaPhi->Sumw2();
+
// Jet pT vs Area
Int_t JetPtAreaBins=200;
Double_t JetPtAreaLow=0.0;
Double_t JetPtAreaUp=2.0;
- fhJetPtArea = new TH2D("fhJetPtArea","Jet Area Distribution",fPtBins,fPtLow,fPtUp,JetPtAreaBins,JetPtAreaLow,JetPtAreaUp);
+ fhJetPtArea = new TH2F("fhJetPtArea","Jet Area Distribution",fPtBins,fPtLow,fPtUp,JetPtAreaBins,JetPtAreaLow,JetPtAreaUp);
fhJetPtArea->GetXaxis()->SetTitle("p_{T} (GeV/c)");
fhJetPtArea->GetYaxis()->SetTitle("A_{jet}");
fhJetPtArea->GetZaxis()->SetTitle("1/N_{Events} dN/dA_{jet}dp_{T}");
fhJetPtArea->Sumw2();
+ // Jet pT vs Constituent pT
+ fhJetConstituentPt = new TH2F("fhJetConstituentPt","Jet constituents p_{T} distribution",fPtBins,fPtLow,fPtUp,10*fPtBins,fPtLow,fPtUp);
+ fhJetConstituentPt->GetXaxis()->SetTitle("Jet p_{T} (GeV/c)");
+ fhJetConstituentPt->GetYaxis()->SetTitle("Constituent p_{T} (GeV/c)");
+ fhJetConstituentPt->GetZaxis()->SetTitle("1/N_{Events} dN/dp_{T,jet}dp_{T,track}");
+ fhJetConstituentPt->Sumw2();
+
+ fhJetTracksPt = new TH2F("fhJetTracksPt","Jet constituents Tracks p_{T} distribution",fPtBins,fPtLow,fPtUp,10*fPtBins,fPtLow,fPtUp);
+ fhJetTracksPt->GetXaxis()->SetTitle("Jet p_{T} (GeV/c)");
+ fhJetTracksPt->GetYaxis()->SetTitle("Constituent Track p_{T} (GeV/c)");
+ fhJetTracksPt->GetZaxis()->SetTitle("1/N_{Events} dN/dp_{T,jet}dp_{T,track}");
+ fhJetTracksPt->Sumw2();
+
+ fhJetClustersPt = new TH2F("fhJetClustersPt","Jet constituents Clusters p_{T} distribution",fPtBins,fPtLow,fPtUp,10*fPtBins,fPtLow,fPtUp);
+ fhJetClustersPt->GetXaxis()->SetTitle("Jet p_{T} (GeV/c)");
+ fhJetClustersPt->GetYaxis()->SetTitle("Constituent Cluster p_{T} (GeV/c)");
+ fhJetClustersPt->GetZaxis()->SetTitle("1/N_{Events} dN/dp_{T,jet}dp_{T,cluster}");
+ fhJetClustersPt->Sumw2();
+
+ // Jet pT vs Constituent Counts
+ fhJetConstituentCounts = new TH2F("fhJetConstituentCounts","Jet constituents distribution",fPtBins,fPtLow,fPtUp,TCBins,0,(Double_t)TCBins);
+ fhJetConstituentCounts->GetXaxis()->SetTitle("Jet p_{T} (GeV/c)");
+ fhJetConstituentCounts->GetYaxis()->SetTitle("Constituent Count");
+ fhJetConstituentCounts->GetZaxis()->SetTitle("1/N_{Events} dN/dp_{T,jet}dN_{constituent}");
+ fhJetConstituentCounts->Sumw2();
+
+ fhJetTracksCounts = new TH2F("fhJetTracksCounts","Jet constituents Tracks distribution",fPtBins,fPtLow,fPtUp,TCBins,0,(Double_t)TCBins);
+ fhJetTracksCounts->GetXaxis()->SetTitle("Jet p_{T} (GeV/c)");
+ fhJetTracksCounts->GetYaxis()->SetTitle("Constituent Track Count");
+ fhJetTracksCounts->GetZaxis()->SetTitle("1/N_{Events} dN/dp_{T,jet}dN_{track}");
+ fhJetTracksCounts->Sumw2();
+
+ fhJetClustersCounts = new TH2F("fhJetClustersCounts","Jet constituents Clusters distribution",fPtBins,fPtLow,fPtUp,TCBins,0,(Double_t)TCBins);
+ fhJetClustersCounts->GetXaxis()->SetTitle("Jet p_{T} (GeV/c)");
+ fhJetClustersCounts->GetYaxis()->SetTitle("Constituent Cluster Count");
+ fhJetClustersCounts->GetZaxis()->SetTitle("1/N_{Events} dN/dp_{T,jet}dN_{cluster}");
+ fhJetClustersCounts->Sumw2();
+
+ // Jet pT vs z_{constituent/track/cluster}
+ fhJetPtZConstituent = new TH2F("fhJetPtZConstituent","Jet z_{constituent} distribution",fPtBins,fPtLow,fPtUp,TCBins,0,1.0);
+ fhJetPtZConstituent->GetXaxis()->SetTitle("Jet p_{T} (GeV/c)");
+ fhJetPtZConstituent->GetYaxis()->SetTitle("z_{constituent}");
+ fhJetPtZConstituent->GetZaxis()->SetTitle("1/N_{Events} dN_{jet}/dp_{T,jet}dz_{constituent}");
+ fhJetPtZConstituent->Sumw2();
+
+ fhJetPtZTrack = new TH2F("fhJetPtZTrack","Jet z_{track} distribution",fPtBins,fPtLow,fPtUp,TCBins,0,1.0);
+ fhJetPtZTrack->GetXaxis()->SetTitle("Jet p_{T} (GeV/c)");
+ fhJetPtZTrack->GetYaxis()->SetTitle("z_{track}");
+ fhJetPtZTrack->GetZaxis()->SetTitle("1/N_{Events} dN_{jet}/dp_{T,jet}dz_{track}");
+ fhJetPtZTrack->Sumw2();
+
+ fhJetPtZCluster = new TH2F("fhJetPtZCluster","Jet z_{cluster} distribution",fPtBins,fPtLow,fPtUp,TCBins,0,1.0);
+ fhJetPtZCluster->GetXaxis()->SetTitle("Jet p_{T} (GeV/c)");
+ fhJetPtZCluster->GetYaxis()->SetTitle("z_{cluster}");
+ fhJetPtZCluster->GetZaxis()->SetTitle("1/N_{Events} dN_{jet}/dp_{T,jet}dz_{cluster}");
+ fhJetPtZCluster->Sumw2();
+
+ // Jet pT vs z_Leading{constituent/track/cluster}
+ fhJetPtZLeadingConstituent = new TH2F("fhJetPtZLeadingConstituent","Jet z_{Leading,constituent} distribution",fPtBins,fPtLow,fPtUp,TCBins,0,1.0);
+ fhJetPtZLeadingConstituent->GetXaxis()->SetTitle("Jet p_{T} (GeV/c)");
+ fhJetPtZLeadingConstituent->GetYaxis()->SetTitle("z_{Leading,constituent}");
+ fhJetPtZLeadingConstituent->GetZaxis()->SetTitle("1/N_{Events} dN_{jet}/dp_{T,jet}dz_{constituent}");
+ fhJetPtZLeadingConstituent->Sumw2();
+
+ fhJetPtZLeadingTrack = new TH2F("fhJetPtZLeadingTrack","Jet z_{Leading,track} distribution",fPtBins,fPtLow,fPtUp,TCBins,0,1.0);
+ fhJetPtZLeadingTrack->GetXaxis()->SetTitle("Jet p_{T} (GeV/c)");
+ fhJetPtZLeadingTrack->GetYaxis()->SetTitle("z_{Leading,track}");
+ fhJetPtZLeadingTrack->GetZaxis()->SetTitle("1/N_{Events} dN_{jet}/dp_{T,jet}dz_{track}");
+ fhJetPtZLeadingTrack->Sumw2();
+
+ fhJetPtZLeadingCluster = new TH2F("fhJetPtZLeadingCluster","Jet z_{Leading,cluster} distribution",fPtBins,fPtLow,fPtUp,TCBins,0,1.0);
+ fhJetPtZLeadingCluster->GetXaxis()->SetTitle("Jet p_{T} (GeV/c)");
+ fhJetPtZLeadingCluster->GetYaxis()->SetTitle("z_{Leading,cluster}");
+ fhJetPtZLeadingCluster->GetZaxis()->SetTitle("1/N_{Events} dN_{jet}/dp_{T,jet}dz_{cluster}");
+ fhJetPtZLeadingCluster->Sumw2();
+
+ // Event Centralities vs Leading Jet Pt
+ fhEventCentralityVsZNA = new TH2F("fhEventCentralityVsZNA",EventCentralityString,fCentralityBins,fCentralityLow,fCentralityUp,fCentralityBins,fCentralityLow,fCentralityUp);
+ fhEventCentralityVsZNA->GetXaxis()->SetTitle(Form("%s",CentralityString.Data()));
+ fhEventCentralityVsZNA->GetYaxis()->SetTitle("Centrality (ZNA)");
+ fhEventCentralityVsZNA->GetZaxis()->SetTitle("Probability Density");
+ fhEventCentralityVsZNA->Sumw2();
+
+ EventCentralityString = Form("%s vs ZNA Centrality vs Leading Jet p_{T} ",fCentralityTag.Data());
+ fhEventCentralityVsZNAPt = new TH3F("fhEventCentralityVsZNAPt",EventCentralityString,fCentralityBins,fCentralityLow,fCentralityUp,fCentralityBins,fCentralityLow,fCentralityUp,fPtBins,fPtLow,fPtUp);
+ fhEventCentralityVsZNAPt->GetXaxis()->SetTitle(Form("%s",CentralityString.Data()));
+ fhEventCentralityVsZNAPt->GetYaxis()->SetTitle("Centrality (ZNA)");
+ fhEventCentralityVsZNAPt->GetZaxis()->SetTitle("Leading Jet p_{T} (GeV/c)");
+ fhEventCentralityVsZNAPt->Sumw2();
+
// Neutral Energy Fraction Histograms & QA
if (fDoNEFQAPlots==kTRUE)
{
fNEFOutput->SetOwner();
fNEFOutput->SetName("ListNEFQAPlots");
- Int_t TCBins = 100;
- fhNEF = new TH1D("fhNEF","Neutral Energy Fraction of All Jets",fNEFBins,fNEFLow,fNEFUp);
- fhNEF->GetXaxis()->SetTitle("NEF");
- fhNEF->GetYaxis()->SetTitle("1/N_{Events} dN/dNEF");
- fhNEF->Sumw2();
+ fhJetPtNEF = new TH2F("fhJetPtNEF","Jet p_{T} vs NEF",fPtBins,fPtLow,fPtUp,fNEFBins,fNEFLow,fNEFUp);
+ fhJetPtNEF->GetXaxis()->SetTitle("Jet p_{T} (GeV/c)");
+ fhJetPtNEF->GetYaxis()->SetTitle("Neutral Energy Fraction");
+ fhJetPtNEF->GetZaxis()->SetTitle("1/N_{Events} dN_{jet}/dp_{T}dNEF");
+ fhJetPtNEF->Sumw2();
+
+ SetNEFJetDimensions(10); // Order: nef,Jet Pt,Eta,Phi,Centrality,Constituent mult,Charged mult, Neutral mult, z_leading
+ SetNEFClusterDimensions(11); // Order: nef,Jet Pt,Eta,Phi,Centrality,F_Cross,z_leading,t_cell,deltat_cell,Cell Count, E_Cluster
+
+ Int_t dimJetBins[fnDimJet];
+ Double_t dimJetLow[fnDimJet];
+ Double_t dimJetUp[fnDimJet];
+
+ Int_t dimClusterBins[fnDimCluster];
+ Double_t dimClusterLow[fnDimCluster];
+ Double_t dimClusterUp[fnDimCluster];
+
+ // Establish dimensinality bin counts and bounds
+ // NEF
+ dimJetBins[0] = dimClusterBins[0] = fNEFBins;
+ dimJetLow[0] = dimClusterLow[0] = fNEFLow;
+ dimJetUp[0] = dimClusterUp[0] = fNEFUp;
+
+ // Jet Pt
+ dimJetBins[1] = dimClusterBins[1] = fPtBins;
+ dimJetLow[1] = dimClusterLow[1] = fPtLow;
+ dimJetUp[1] = dimClusterUp[1] = fPtUp;
+
+ // Eta-Phi
+ dimJetBins[2] = dimJetBins[3] = dimClusterBins[2] = dimClusterBins[3] = TCBins;
+ dimJetLow[2] = dimClusterLow[2] = fEMCalEtaMin;
+ dimJetUp[2] = dimClusterUp[2] = fEMCalEtaMax;
+ dimJetLow[3] = dimClusterLow[3] = fEMCalPhiMin;
+ dimJetUp[3] = dimClusterUp[3] = fEMCalPhiMax;
- fhNEFSignal = new TH1D("fhNEFSignal","Neutral Energy Fraction of Signal Jets",fNEFBins,fNEFLow,fNEFUp);
- fhNEFSignal->GetXaxis()->SetTitle("NEF");
- fhNEFSignal->GetYaxis()->SetTitle("1/N_{Events} dN/dNEF");
- fhNEFSignal->Sumw2();
-
- fhNEFJetPt = new TH2D("fhNEFJetPt","Neutral Energy Fraction vs p_{T}^{jet}",fNEFBins,fNEFLow,fNEFUp,fPtBins,fPtLow,fPtUp);
- fhNEFJetPt->GetXaxis()->SetTitle("NEF");
- fhNEFJetPt->GetYaxis()->SetTitle("p_{T}^{jet} (GeV/c)");
- fhNEFJetPt->GetZaxis()->SetTitle("1/N_{Events} dN/dNEFdp_{T}");
- fhNEFJetPt->Sumw2();
-
- fhNEFJetPtSignal = new TH2D("fhNEFJetPtSignal","Neutral Energy Fraction vs p_{T}^{jet}",fNEFBins,fNEFLow,fNEFUp,fPtBins,fPtLow,fPtUp);
- fhNEFJetPtSignal->GetXaxis()->SetTitle("NEF");
- fhNEFJetPtSignal->GetYaxis()->SetTitle("p_{T}^{jet} (GeV/c)");
- fhNEFJetPtSignal->GetZaxis()->SetTitle("1/N_{Events} dN/dNEFdp_{T}");
- fhNEFJetPtSignal->Sumw2();
-
- // Eta-Phi Dependence
- fhNEFEtaPhi = new TH2D("fhNEFEtaPhi","Neutral Energy Fraction #eta-#varphi",TCBins, fEMCalEtaMin,fEMCalEtaMax,TCBins,fEMCalPhiMin,fEMCalPhiMax);
- fhNEFEtaPhi->GetXaxis()->SetTitle("#eta");
- fhNEFEtaPhi->GetYaxis()->SetTitle("#varphi");
- fhNEFEtaPhi->GetZaxis()->SetTitle("1/N{Events} dN/d#etad#varphi");
- fhNEFEtaPhi->Sumw2();
-
- fhNEFEtaPhiSignal = new TH2D("fhNEFEtaPhiSignal","Neutral Energy Fraction #eta-#varphi of Signal Jets",TCBins,fEMCalEtaMin,fEMCalEtaMax,TCBins,fEMCalPhiMin,fEMCalPhiMax);
- fhNEFEtaPhiSignal->GetXaxis()->SetTitle("#eta");
- fhNEFEtaPhiSignal->GetYaxis()->SetTitle("#varphi");
- fhNEFEtaPhiSignal->GetZaxis()->SetTitle("1/N{Events} dN/d#etad#varphi");
- fhNEFEtaPhiSignal->Sumw2();
-
- fhEtaPhiNEF = new TH3D("fhEtaPhiNEF","Neutral Energy Fraction #eta-#varphi",TCBins, fEMCalEtaMin,fEMCalEtaMax,TCBins,fEMCalPhiMin,fEMCalPhiMax,fNEFBins,fNEFLow,fNEFUp);
- fhEtaPhiNEF->GetXaxis()->SetTitle("#eta");
- fhEtaPhiNEF->GetYaxis()->SetTitle("#varphi");
- fhEtaPhiNEF->GetZaxis()->SetTitle("NEF");
- fhEtaPhiNEF->Sumw2();
-
- // Multiplicity Dependence
- fhNEFTotalMult = new TH2D("fhNEFTotalMult","Total Multiplicity Distribution of Constituents",fNEFBins,fNEFLow,fNEFUp,TCBins,0,(Double_t)TCBins);
- fhNEFTotalMult->GetXaxis()->SetTitle("NEF");
- fhNEFTotalMult->GetYaxis()->SetTitle("Multiplicity");
- fhNEFTotalMult->GetZaxis()->SetTitle("1/N_{Events} dN/dNEFdM");
- fhNEFTotalMult->Sumw2();
+ // Centrality
+ dimJetBins[4] = dimClusterBins[4] = fCentralityBins;
+ dimJetLow[4] = dimClusterLow[4] = fCentralityLow;
+ dimJetUp[4] = dimClusterUp[4] = fCentralityUp;
- fhNEFTotalMultSignal = new TH2D("fhNEFTotalMultSignal","Total Multiplicity Distribution of Constituents of Signal Jets",fNEFBins,fNEFLow,fNEFUp,TCBins,0,(Double_t)TCBins);
- fhNEFTotalMultSignal->GetXaxis()->SetTitle("NEF");
- fhNEFTotalMultSignal->GetYaxis()->SetTitle("Multiplicity");
- fhNEFTotalMultSignal->GetZaxis()->SetTitle("1/N_{Events} dN/dNEFdM");
- fhNEFTotalMultSignal->Sumw2();
+ // z_leading
+ dimJetBins[5] = dimClusterBins[5] = TCBins;
+ dimJetLow[5] = dimClusterLow[5] = 0.0;
+ dimJetUp[5] = dimClusterUp[5] = 1.0;
- fhNEFNeutralMult = new TH2D("fhNEFNeutralMult","Neutral Multiplicity Distribution of Constituents",fNEFBins,fNEFLow,fNEFUp,TCBins,0,(Double_t)TCBins);
- fhNEFNeutralMult->GetXaxis()->SetTitle("NEF");
- fhNEFNeutralMult->GetYaxis()->SetTitle("Multiplicity");
- fhNEFNeutralMult->GetZaxis()->SetTitle("1/N_{Events} dN/dNEFdM");
- fhNEFNeutralMult->Sumw2();
+ // Jets Constituent Multiplicity Info {Total,Charged,Neutral}
+ for (i=6;i<9;i++)
+ {
+ dimJetBins[i] = TCBins;
+ dimJetLow[i] = 0.0;
+ dimJetUp[i] = (Double_t)TCBins;
+ }
+
+ // z_leading^track
+ dimJetBins[9] = TCBins;
+ dimJetLow[9] = 0.0;
+ dimJetUp[9] = 1.0;
+
+ // Cluster E
+ dimClusterBins[6] = fPtBins;
+ dimClusterLow[6] = fPtLow;
+ dimClusterUp[6] = fPtUp;
- fhNEFNeutralMultSignal = new TH2D("fhNEFNeutralMultSignal","Neutral Multiplicity Distribution of Constituents of Signal Jets",fNEFBins,fNEFLow,fNEFUp,TCBins,0,(Double_t)TCBins);
- fhNEFNeutralMultSignal->GetXaxis()->SetTitle("NEF");
- fhNEFNeutralMultSignal->GetYaxis()->SetTitle("Multiplicity");
- fhNEFNeutralMultSignal->GetZaxis()->SetTitle("1/N_{Events} dN/dNEFdM");
- fhNEFNeutralMultSignal->Sumw2();
+ // Cluster F_Cross
+ dimClusterBins[7] = TCBins;
+ dimClusterLow[7] = 0.0;
+ dimClusterUp[7] = 1.0;
- // Cluster Shape QA
- fhClusterShapeAll = new TH1D("fhClusterShapeAll","Cluster Shape of all CaloClustersCorr",10*TCBins,0,10*TCBins);
- fhClusterShapeAll->GetXaxis()->SetTitle("Cells");
- fhClusterShapeAll->GetYaxis()->SetTitle("1/N_{Events} dN/dCells");
- fhClusterShapeAll->Sumw2();
-
- fhClusterPtCellAll = new TH2D("fhClusterPtCellAll","Cluster p_{T} vs Cluster Shape of all CaloClustersCorr",fPtBins,fPtLow,fPtUp,10*TCBins,0,10*TCBins);
- fhClusterPtCellAll->GetXaxis()->SetTitle("p_{T} (GeV/c)");
- fhClusterPtCellAll->GetYaxis()->SetTitle("Cells");
- fhClusterPtCellAll->GetZaxis()->SetTitle("1/N_{Events} dN/dp_{T}dCells");
- fhClusterPtCellAll->Sumw2();
-
- fhNEFJetPtFCross = new TH3D("fhNEFJetPtFCross","NEF vs p_{T}^{jet} vs F_{Cross}",fNEFBins,fNEFLow,fNEFUp,fPtBins,fPtLow,fPtUp,TCBins,0,1.0);
- fhNEFJetPtFCross->GetXaxis()->SetTitle("NEF");
- fhNEFJetPtFCross->GetYaxis()->SetTitle("p_{T}^{jet} (GeV/c)");
- fhNEFJetPtFCross->GetZaxis()->SetTitle("F_{Cross}");
- fhNEFJetPtFCross->Sumw2();
-
- fhNEFZLeadingFCross = new TH3D("fhNEFZLeadingFCross","NEF vs z_{Leading} vs F_{Cross}",fNEFBins,fNEFLow,fNEFUp,TCBins,0,1.0,TCBins,0,1.0);
- fhNEFZLeadingFCross->GetXaxis()->SetTitle("NEF");
- fhNEFZLeadingFCross->GetYaxis()->SetTitle("z_{Leading}");
- fhNEFZLeadingFCross->GetZaxis()->SetTitle("F_{Cross}");
- fhNEFZLeadingFCross->Sumw2();
-
- fhNEFTimeCellCount = new TH3D("fhNEFTimeCellCount","NEF vs t_{cell} vs Cell Count",fNEFBins,fNEFLow,fNEFUp,400,-2e-07,2e-07,TCBins,0,100);
- fhNEFTimeCellCount->GetXaxis()->SetTitle("NEF");
- fhNEFTimeCellCount->GetYaxis()->SetTitle("t_{cell} (s)");
- fhNEFTimeCellCount->GetZaxis()->SetTitle("Cell Count");
- fhNEFTimeCellCount->Sumw2();
-
- fhNEFTimeDeltaTime = new TH3D("fhNEFTimeDeltaTime","NEF vs t_{cell} vs #Deltat",fNEFBins,fNEFLow,fNEFUp,400,-2e-07,2e-06,TCBins,0,1e-07);
- fhNEFTimeDeltaTime->GetXaxis()->SetTitle("NEF");
- fhNEFTimeDeltaTime->GetYaxis()->SetTitle("t_{cell} (s)");
- fhNEFTimeDeltaTime->GetZaxis()->SetTitle("#Deltat");
- fhNEFTimeDeltaTime->Sumw2();
+ // Cluster t_cell
+ dimClusterBins[8] = 400;
+ dimClusterLow[8] = -2e-07;
+ dimClusterUp[8] = 2e-07;
+
+ // Cluster delta t_cell
+ dimClusterBins[9] = 100;
+ dimClusterLow[9] = 0.0;
+ dimClusterUp[9] = 1e-07;
- fNEFOutput->Add(fhNEF);
- fNEFOutput->Add(fhNEFSignal);
- fNEFOutput->Add(fhNEFJetPt);
- fNEFOutput->Add(fhNEFJetPtSignal);
- fNEFOutput->Add(fhNEFEtaPhi);
- fNEFOutput->Add(fhNEFEtaPhiSignal);
- fNEFOutput->Add(fhEtaPhiNEF);
- fNEFOutput->Add(fhNEFTotalMult);
- fNEFOutput->Add(fhNEFTotalMultSignal);
- fNEFOutput->Add(fhNEFNeutralMult);
- fNEFOutput->Add(fhNEFNeutralMultSignal);
- fNEFOutput->Add(fhClusterShapeAll);
- fNEFOutput->Add(fhClusterPtCellAll);
- fNEFOutput->Add(fhNEFJetPtFCross);
- fNEFOutput->Add(fhNEFZLeadingFCross);
- fNEFOutput->Add(fhNEFTimeCellCount);
- fNEFOutput->Add(fhNEFTimeDeltaTime);
+ // Cluster Cell Count
+ dimClusterBins[10] = TCBins;
+ dimClusterLow[10] = 0.0;
+ dimClusterUp[10] = 100.0;
+
+ if (fDoTHnSparse == kTRUE)
+ {
+ fhJetNEFSignalInfo = new THnSparseF("fhJetNEFSignalInfo","Signal Jet NEF Information Histogram",fnDimJet,dimJetBins,dimJetLow,dimJetUp);
+ fhJetNEFSignalInfo->Sumw2();
+
+ fhClusterNEFSignalInfo = new THnSparseF("fhClusterNEFSignalInfo","Signal Jet NEF Cluster Information Histogram",fnDimCluster,dimClusterBins,dimClusterLow,dimClusterUp);
+ fhClusterNEFSignalInfo->Sumw2();
+
+ // Cluster Shape QA
+ fhClusterShapeAll = new TH1F("fhClusterShapeAll","Cluster Shape of all CaloClustersCorr",10*TCBins,0,10*TCBins);
+ fhClusterShapeAll->GetXaxis()->SetTitle("Cells");
+ fhClusterShapeAll->GetYaxis()->SetTitle("1/N_{Events} dN/dCells");
+ fhClusterShapeAll->Sumw2();
+
+ fhClusterPtCellAll = new TH2F("fhClusterPtCellAll","Cluster p_{T} vs Cluster Shape of all CaloClustersCorr",fPtBins,fPtLow,fPtUp,10*TCBins,0,10*TCBins);
+ fhClusterPtCellAll->GetXaxis()->SetTitle("p_{T} (GeV/c)");
+ fhClusterPtCellAll->GetYaxis()->SetTitle("Cells");
+ fhClusterPtCellAll->GetZaxis()->SetTitle("1/N_{Events} dN/dp_{T}dCells");
+ fhClusterPtCellAll->Sumw2();
+
+ if (fDoNEFSignalOnly == kFALSE)
+ {
+ fhJetNEFInfo = new THnSparseF("fhJetNEFInfo","Jet NEF Information Histogram",fnDimJet,dimJetBins,dimJetLow,dimJetUp);
+ fhJetNEFInfo->Sumw2();
+
+ fhClusterNEFInfo = new THnSparseF("fhClusterNEFInfo","Jet NEF Cluster Information Histogram",fnDimCluster,dimClusterBins,dimClusterLow,dimClusterUp);
+ fhClusterNEFInfo->Sumw2();
+
+ fNEFOutput->Add(fhJetNEFInfo);
+ fNEFOutput->Add(fhClusterNEFInfo);
+ }
+ fNEFOutput->Add(fhJetNEFSignalInfo);
+ fNEFOutput->Add(fhClusterNEFSignalInfo);
+ fNEFOutput->Add(fhClusterShapeAll);
+ fNEFOutput->Add(fhClusterPtCellAll);
+ }
+ fNEFOutput->Add(fhJetPtNEF);
fOutput->Add(fNEFOutput);
}
fOutput->Add(fh80100BSPt);
fOutput->Add(fhBSPt);
fOutput->Add(fhBSPtCen);
- //fOutput->Add(fhBSPtCenLCT);
fOutput->Add(fh020BSPtSignal);
fOutput->Add(fh80100BSPtSignal);
fOutput->Add(fhBSPtSignal);
fOutput->Add(fhBckgFlucPtCen);
fOutput->Add(fpRho);
fOutput->Add(fpLJetRho);
- fOutput->Add(fhJetConstituentPt);
+ fOutput->Add(fhJetPtEtaPhi);
fOutput->Add(fhJetPtArea);
+ fOutput->Add(fhJetConstituentPt);
+ fOutput->Add(fhJetTracksPt);
+ fOutput->Add(fhJetClustersPt);
+ fOutput->Add(fhJetConstituentCounts);
+ fOutput->Add(fhJetTracksCounts);
+ fOutput->Add(fhJetClustersCounts);
+ fOutput->Add(fhJetPtZConstituent);
+ fOutput->Add(fhJetPtZTrack);
+ fOutput->Add(fhJetPtZCluster);
+ fOutput->Add(fhJetPtZLeadingConstituent);
+ fOutput->Add(fhJetPtZLeadingTrack);
+ fOutput->Add(fhJetPtZLeadingCluster);
+ fOutput->Add(fhEventCentralityVsZNA);
+ fOutput->Add(fhEventCentralityVsZNAPt);
}
void AliAnalysisTaskFullpAJets::AlipAJetHistos::SetName(const char *name)
fName = name;
}
-void AliAnalysisTaskFullpAJets::AlipAJetHistos::SetCentralityTag(const char *name)
+void AliAnalysisTaskFullpAJets::AlipAJetHistos::SetCentralityTag(TString name)
{
- fCentralityTag = name;
+ fCentralityTag = name.Data();
}
void AliAnalysisTaskFullpAJets::AlipAJetHistos::SetCentralityRange(Int_t bins, Double_t low, Double_t up)
fNEFUp=up;
}
+void AliAnalysisTaskFullpAJets::AlipAJetHistos::SetSignalTrackPtBias(Bool_t chargedBias)
+{
+ fSignalTrackBias = chargedBias;
+}
+
+void AliAnalysisTaskFullpAJets::AlipAJetHistos::SetNEFJetDimensions(Int_t n)
+{
+ fnDimJet = n;
+}
+
+void AliAnalysisTaskFullpAJets::AlipAJetHistos::SetNEFClusterDimensions(Int_t n)
+{
+ fnDimCluster = n;
+}
+
+void AliAnalysisTaskFullpAJets::AlipAJetHistos::SetRhoValue(Double_t value)
+{
+ fRhoValue = value;
+}
+
+void AliAnalysisTaskFullpAJets::AlipAJetHistos::DoTHnSparse(Bool_t doTHnSparse)
+{
+ fDoTHnSparse = doTHnSparse;
+}
+
void AliAnalysisTaskFullpAJets::AlipAJetHistos::FillRho(Double_t eventCentrality, Double_t rho)
{
- fRhoValue = rho;
-
+ SetRhoValue(rho);
fhRho->Fill(rho);
fhRhoCen->Fill(rho,eventCentrality);
fpRho->Fill(eventCentrality,rho);
fhBSPt->Fill(tempPt);
fhBSPtCen->Fill(tempPt,eventCentrality);
- //fhBSPtCenLCT->Fill(tempPt,eventCentrality,tempChargedHighPt);
if (eventCentrality<=20)
{
fh020BSPt->Fill(tempPt);
{
fh80100BSPt->Fill(tempPt);
}
- if (tempChargedHighPt>=signalCut)
+ if (fSignalTrackBias==kTRUE)
{
- fhBSPtSignal->Fill(tempPt);
- fhBSPtCenSignal->Fill(tempPt,eventCentrality);
- if (eventCentrality<=20)
+ if (tempChargedHighPt>=signalCut)
{
- fh020BSPtSignal->Fill(tempPt);
+ fhBSPtSignal->Fill(tempPt);
+ fhBSPtCenSignal->Fill(tempPt,eventCentrality);
+ if (eventCentrality<=20)
+ {
+ fh020BSPtSignal->Fill(tempPt);
+ }
+ else if (eventCentrality>=80)
+ {
+ fh80100BSPtSignal->Fill(tempPt);
+ }
}
- else if (eventCentrality>=80)
+ }
+ else
+ {
+ if (tempPt>=signalCut)
{
- fh80100BSPtSignal->Fill(tempPt);
+ fhBSPtSignal->Fill(tempPt);
+ fhBSPtCenSignal->Fill(tempPt,eventCentrality);
+ if (eventCentrality<=20)
+ {
+ fh020BSPtSignal->Fill(tempPt);
+ }
+ else if (eventCentrality>=80)
+ {
+ fh80100BSPtSignal->Fill(tempPt);
+ }
}
}
tempPt=0.0;
fDoNEFQAPlots = doNEFAna;
}
+void AliAnalysisTaskFullpAJets::AlipAJetHistos::DoNEFSignalOnly(Bool_t doNEFSignalOnly)
+{
+ fDoNEFSignalOnly = doNEFSignalOnly;
+}
+
void AliAnalysisTaskFullpAJets::AlipAJetHistos::DoNEFAnalysis(Double_t nefCut, Double_t signalCut, TClonesArray *jetList, Int_t *indexJetList, Int_t nIndexJetList, TObjArray *clusterList, TClonesArray *orgClusterList, AliVEvent *event, AliEMCALGeometry *geometry, AliEMCALRecoUtils *recoUtils, AliVCaloCells *cells)
{
if (fDoNEFQAPlots==kFALSE)
return;
}
+ if (fDoTHnSparse == kFALSE)
+ {
+ Int_t i;
+
+ Double_t nef=0.0;
+ Double_t jetPt=0.0;
+ Double_t tempChargedHighPt=0.0;
+
+ for (i=0;i<nIndexJetList;i++)
+ {
+ AliEmcalJet *myJet = (AliEmcalJet*) jetList->At(indexJetList[i]);
+ tempChargedHighPt = myJet->MaxTrackPt();
+ nef=myJet->NEF();
+ jetPt=myJet->Pt();
+
+ if (fSignalTrackBias==kTRUE)
+ {
+ if (tempChargedHighPt>=signalCut && nef<=nefCut)
+ {
+ fhJetPtNEF->Fill(jetPt,nef);
+ }
+ }
+ else
+ {
+ if (jetPt>=signalCut && nef<=nefCut)
+ {
+ fhJetPtNEF->Fill(jetPt,nef);
+ }
+ }
+
+ nef=0.0;
+ jetPt=0.0;
+ tempChargedHighPt=0.0;
+ }
+ return;
+ }
+
Int_t i,j,k;
+
+ Double_t valJet[fnDimJet];
+ Double_t valCluster[fnDimJet];
+ for (i=0;i<fnDimJet;i++)
+ {
+ valJet[i]=0.0;
+ valCluster[i]=0.0;
+ }
+
Double_t nef=0.0;
+ Double_t jetPt=0.0;
Double_t tempChargedHighPt=0.0;
Double_t eta=0.0;
Double_t phi=0.0;
Int_t totalMult=0;
+ Int_t chargedMult=0;
Int_t neutralMult=0;
Int_t iSupMod = -1, absId = -1, ieta = -1, iphi = -1;
Bool_t shared = kFALSE;
Double_t zLeading=0.0;
+ Double_t zLeadingTrack=0.0;
+ Double_t ECluster=0.0;
Double_t eSeed=0.0;
Double_t tCell=0.0;
Double_t eCross=0.0;
Int_t tempCellID=0;
Double_t tempCellTime=0.0;
+ Double_t event_centrality = event->GetCentrality()->GetCentralityPercentile(fCentralityTag);
+ valJet[4] = valCluster[4] = event_centrality;
+
// First, do Jet QA
for (i=0;i<nIndexJetList;i++)
{
AliEmcalJet *myJet = (AliEmcalJet*) jetList->At(indexJetList[i]);
tempChargedHighPt = myJet->MaxTrackPt();
nef=myJet->NEF();
+ jetPt=myJet->Pt();
eta=myJet->Eta();
phi=myJet->Phi();
totalMult=myJet->GetNumberOfConstituents();
+ chargedMult = myJet->GetNumberOfTracks();
neutralMult=myJet->GetNumberOfClusters();
- zLeading=myJet->MaxClusterPt()/myJet->Pt();
+ zLeading=myJet->MaxPartPt()/myJet->Pt();
+ zLeadingTrack=myJet->MaxTrackPt()/myJet->Pt();
- fhNEF->Fill(nef);
- fhNEFJetPt->Fill(nef,myJet->Pt());
- fhNEFEtaPhi->Fill(eta,phi);
- fhEtaPhiNEF->Fill(eta,phi,nef);
- fhNEFTotalMult->Fill(nef,totalMult);
- fhNEFNeutralMult->Fill(nef,neutralMult);
+ valJet[0] = valCluster[0] = nef;
+ valJet[1] = valCluster[1] = jetPt;
+ valJet[2] = valCluster[2] = eta;
+ valJet[3] = valCluster[3] = phi;
+ valJet[5] = valCluster[5] = zLeading;
+ valJet[9] = zLeadingTrack;
+ valJet[6] = totalMult;
+ valJet[7] = chargedMult;
+ valJet[8] = neutralMult;
+
+ // Supress filling of this histogram due to memory size of THnSparse when running over large datasets
+ if (fDoNEFSignalOnly == kFALSE)
+ {
+ fhJetNEFInfo->Fill(valJet);
+ }
+
+ if (fSignalTrackBias==kTRUE)
+ {
+ if (tempChargedHighPt>=signalCut && nef<=nefCut)
+ {
+ fhJetNEFSignalInfo->Fill(valJet);
+ fhJetPtNEF->Fill(jetPt,nef);
+ }
+ }
+ else
+ {
+ if (jetPt>=signalCut && nef<=nefCut)
+ {
+ fhJetNEFSignalInfo->Fill(valJet);
+ fhJetPtNEF->Fill(jetPt,nef);
+ }
+ }
+
for (j=0;j<neutralMult;j++)
{
AliVCluster* vcluster = (AliVCluster*) orgClusterList->At(myJet->ClusterAt(j));
+ ECluster = vcluster->E();
recoUtils->GetMaxEnergyCell(geometry,cells,vcluster,absId,iSupMod,ieta,iphi,shared);
eSeed = cells->GetCellAmplitude(absId);
tCell = cells->GetCellTime(absId);
eCross = recoUtils->GetECross(absId,tCell,cells,event->GetBunchCrossNumber());
FCross = 1 - eCross/eSeed;
- fhNEFJetPtFCross->Fill(nef,myJet->Pt(),FCross);
- fhNEFZLeadingFCross->Fill(nef,zLeading,FCross);
- fhNEFTimeCellCount->Fill(nef,tCell,vcluster->GetNCells());
-
// Obtain Delta t of Cluster
lowTime=9.99e99;
upTime=-9.99e99;
lowTime=tempCellTime;
}
}
- fhNEFTimeDeltaTime->Fill(nef,tCell,upTime-lowTime);
+ valCluster[6] = ECluster;
+ valCluster[7] = FCross;
+ valCluster[8] = tCell;
+ valCluster[9] = upTime-lowTime;
+ valCluster[10] = vcluster->GetNCells();
+ if (fDoNEFSignalOnly == kFALSE)
+ {
+ fhClusterNEFInfo->Fill(valCluster);
+ }
+
+ if (fSignalTrackBias==kTRUE)
+ {
+ if (tempChargedHighPt>=signalCut && nef<=nefCut)
+ {
+ fhClusterNEFSignalInfo->Fill(valCluster);
+ }
+ }
+ else
+ {
+ if (myJet->Pt()>=signalCut && nef<=nefCut)
+ {
+ fhClusterNEFSignalInfo->Fill(valCluster);
+ }
+ }
tempCellID=0;
tempCellTime=0.0;
eSeed=0.0;
iSupMod=-1,absId=-1,ieta=-1,iphi=-1;
}
- if (tempChargedHighPt>=signalCut)
- {
- if (nef<=nefCut)
- {
- fhNEFSignal->Fill(nef);
- fhNEFJetPtSignal->Fill(nef,myJet->Pt());
- fhNEFEtaPhiSignal->Fill(eta,phi);
- fhNEFTotalMultSignal->Fill(nef,totalMult);
- fhNEFNeutralMultSignal->Fill(nef,neutralMult);
- }
- }
nef=0.0;
+ jetPt=0.0;
tempChargedHighPt=0.0;
eta=0.0;
phi=0.0;
totalMult=0;
+ chargedMult=0;
neutralMult=0;
zLeading=0.0;
+ zLeadingTrack=0.0;
+ ECluster=0.0;
}
// Now do Cluster QA
- // Finally, Cluster QA
for (i=0;i<clusterList->GetEntries();i++)
{
AliVCluster *vcluster = (AliVCluster*) clusterList->At(i);
}
}
-void AliAnalysisTaskFullpAJets::AlipAJetHistos::FillMiscJetStats(TClonesArray *jetList, Int_t *indexJetList, Int_t nIndexJetList, TClonesArray *trackList, TClonesArray *clusterList)
+void AliAnalysisTaskFullpAJets::AlipAJetHistos::FillMiscJetStats(TClonesArray *jetList, Int_t *indexJetList, Int_t nIndexJetList, TClonesArray *trackList, TClonesArray *clusterList, Double_t *vertex)
{
Int_t i,j;
-
+
+ TLorentzVector *cluster_vec = new TLorentzVector;
for (i=0;i<nIndexJetList;i++)
{
AliEmcalJet *myJet = (AliEmcalJet*) jetList->At(indexJetList[i]);
+ fhJetPtEtaPhi->Fill(myJet->Pt(),myJet->Eta(),myJet->Phi());
fhJetPtArea->Fill(myJet->Pt(),myJet->Area());
+ fhJetConstituentCounts->Fill(myJet->Pt(),myJet->GetNumberOfConstituents());
+ fhJetTracksCounts->Fill(myJet->Pt(),myJet->GetNumberOfTracks());
+ fhJetClustersCounts->Fill(myJet->Pt(),myJet->GetNumberOfClusters());
+ fhJetPtZLeadingConstituent->Fill(myJet->Pt(),myJet->MaxPartPt()/myJet->Pt());
+ fhJetPtZLeadingTrack->Fill(myJet->Pt(),myJet->MaxTrackPt()/myJet->Pt());
+ fhJetPtZLeadingCluster->Fill(myJet->Pt(),myJet->MaxClusterPt()/myJet->Pt());
for (j=0;j<myJet->GetNumberOfTracks();j++)
{
AliVTrack *vtrack = (AliVTrack*) myJet->TrackAt(j,trackList);
fhJetConstituentPt->Fill(myJet->Pt(),vtrack->Pt());
+ fhJetTracksPt->Fill(myJet->Pt(),vtrack->Pt());
+ fhJetPtZTrack->Fill(myJet->Pt(),vtrack->Pt()/myJet->Pt());
+ fhJetPtZConstituent->Fill(myJet->Pt(),vtrack->Pt()/myJet->Pt());
}
for (j=0;j<myJet->GetNumberOfClusters();j++)
{
AliVCluster *vcluster = (AliVCluster*) myJet->ClusterAt(j,clusterList);
- fhJetConstituentPt->Fill(myJet->Pt(),vcluster->E());
+ vcluster->GetMomentum(*cluster_vec,vertex);
+ fhJetConstituentPt->Fill(myJet->Pt(),cluster_vec->Pt());
+ fhJetClustersPt->Fill(myJet->Pt(),vcluster->E());
+ fhJetPtZCluster->Fill(myJet->Pt(),cluster_vec->Pt()/myJet->Pt());
+ fhJetPtZConstituent->Fill(myJet->Pt(),cluster_vec->Pt()/myJet->Pt());
}
}
+ delete cluster_vec;
+}
+
+void AliAnalysisTaskFullpAJets::AlipAJetHistos::FillJetEventCentrality(Double_t leadingJetPt, AliVEvent *event)
+{
+ Double_t event_centrality = event->GetCentrality()->GetCentralityPercentile(fCentralityTag);
+ Double_t event_centrality_ZNA = event->GetCentrality()->GetCentralityPercentile("ZNA");
+
+ fhEventCentralityVsZNA->Fill(event_centrality,event_centrality_ZNA);
+ fhEventCentralityVsZNAPt->Fill(event_centrality,event_centrality_ZNA,leadingJetPt);
}
TList* AliAnalysisTaskFullpAJets::AlipAJetHistos::GetOutputHistos()
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff