}
else {
Double_t centWidth = (fMaxCent-fMinCent)/(Double_t)fNcentBins;
- fCentBin = TMath::Nint(fCent/centWidth);
+ fCentBin = TMath::FloorNint(fCent/centWidth);
+ if(fCentBin>=fNcentBins) {
+ AliWarning(Form("%s: fCentBin too large: cent = %f fCentBin = %d. Assuming 99", GetName(),fCent,fCentBin));
+ fCentBin = fNcentBins-1;
+ }
if(fCentBin>=fNcentBins) {
AliWarning(Form("%s: fCentBin too large: cent = %f fCentBin = %d. Assuming 99", GetName(),fCent,fCentBin));
fCentBin = fNcentBins-1;
}
+//________________________________________________________________________
+AliEmcalJet* AliAnalysisTaskEmcalJetDev::GetJetFromArray(Int_t j, Int_t c) const {
+ // Get jet j if accepted from container c
+ // If jet not accepted return 0
+
+ AliJetContainer *cont = GetJetContainer(c);
+ if(!cont) {
+ AliError(Form("%s:Container %d not found",GetName(),c));
+ return 0;
+ }
+ AliEmcalJet *jet = cont->GetJet(j);
+
+ return jet;
+
+}
+
//________________________________________________________________________
AliEmcalJet* AliAnalysisTaskEmcalJetDev::GetAcceptJetFromArray(Int_t j, Int_t c) const {
// Get jet j if accepted from container c
AliJetContainer *GetJetContainer(Int_t i) const;
TClonesArray *GetJetArray(Int_t i) const;
+ AliEmcalJet *GetJetFromArray(Int_t j, Int_t c) const;
AliEmcalJet *GetAcceptJetFromArray(Int_t j, Int_t c) const;
Int_t GetNJets(Int_t i) const;
Double_t GetRhoVal(Int_t i) const;
enum JetAcceptanceType {
kTPC = 0, // TPC acceptance
kEMCAL = 1, // EMCal acceptance
- kUser = 2, // User defined acceptance
+ kUser = 2 // User defined acceptance
};
const TString& GetRhoName() const {return fRhoName;}
Double_t GetJetPtCorr(Int_t i) const;
Float_t GetJetRadius() const {return fJetRadius;}
+ Float_t GetJetEtaMin() const {return fJetMinEta;}
+ Float_t GetJetEtaMax() const {return fJetMaxEta;}
+ Float_t GetJetPhiMin() const {return fJetMinPhi;}
+ Float_t GetJetPhiMax() const {return fJetMaxPhi;}
protected:
JetAcceptanceType fJetAcceptanceType; // acceptance type