fHistJetsCorrPtArea(0),
fHistPtDEtaDPhiTrackClus(0),
fHistPtDEtaDPhiClusTrack(0),
+ fHistClustDx(0),
+ fHistClustDz(0),
fJetsCont(0),
fTracksCont(0),
fCaloClustersCont(0)
fHistJetsCorrPtArea(0),
fHistPtDEtaDPhiTrackClus(0),
fHistPtDEtaDPhiClusTrack(0),
+ fHistClustDx(0),
+ fHistClustDz(0),
fJetsCont(0),
fTracksCont(0),
fCaloClustersCont(0)
fTracksCont = GetParticleContainer(0);
fCaloClustersCont = GetClusterContainer(0);
}
- fTracksCont->SetClassName("AliVTrack");
- fCaloClustersCont->SetClassName("AliAODCaloCluster");
+ if(fTracksCont) fTracksCont->SetClassName("AliVTrack");
+ if(fCaloClustersCont) fCaloClustersCont->SetClassName("AliVCluster");
TString histname;
fHistPtDEtaDPhiClusTrack = new TH3F(histname.Data(),Form("%s;#it{p}_{T}^{clus};#Delta#eta;#Delta#varphi",histname.Data()),100,0.,100.,100,-0.1,0.1,100,-0.1,0.1);
fOutput->Add(fHistPtDEtaDPhiClusTrack);
+ fHistClustDx = new TH1F("fHistClustDx","fHistClustDx;Dx",1000,0.,1.);
+ fOutput->Add(fHistClustDx);
+
+ fHistClustDz = new TH1F("fHistClustDz","fHistClustDz;Dz",1000,0.,1.);
+ fOutput->Add(fHistClustDz);
+
PostData(1, fOutput); // Post data for ALL output slots > 0 here.
}
track = static_cast<AliVTrack*>(fTracksCont->GetNextAcceptParticle());
}
}
-
+
if (fCaloClustersCont) {
AliVCluster *cluster = fCaloClustersCont->GetNextAcceptCluster(0);
while(cluster) {
TLorentzVector nPart;
cluster->GetMomentum(nPart, fVertex);
fHistClustersPt[fCentBin]->Fill(nPart.Pt());
-
+ Double_t dx = cluster->GetTrackDx();
+ Double_t dz = cluster->GetTrackDz();
+ fHistClustDx->Fill(dx);
+ fHistClustDz->Fill(dz);
cluster = fCaloClustersCont->GetNextAcceptCluster();
}
}
Float_t ptLeading = fJetsCont->GetLeadingHadronPt(jet);
fHistJetsPtLeadHad[fCentBin]->Fill(jet->Pt(), ptLeading);
- Float_t corrPt = jet->Pt() - fJetsCont->GetRhoVal() * jet->Area();
- fHistJetsCorrPtArea[fCentBin]->Fill(corrPt, jet->Area());
-
+ if (fHistJetsCorrPtArea[fCentBin]) {
+ Float_t corrPt = jet->Pt() - fJetsCont->GetRhoVal() * jet->Area();
+ fHistJetsCorrPtArea[fCentBin]->Fill(corrPt, jet->Area());
+ }
jet = fJetsCont->GetNextAcceptJet();
}