// Destructor
}
-//________________________________________________________________________
-Bool_t AliAnalysisTaskEmcalJet::RetrieveEventObjects()
-{
- // Retrieve objects from event.
-
- if (!AliAnalysisTaskEmcal::RetrieveEventObjects())
- return kFALSE;
-
- if (!fJetsName.IsNull() && !fJets) {
- fJets = dynamic_cast<TClonesArray*>(InputEvent()->FindListObject(fJetsName));
- if (!fJets) {
- AliWarning(Form("%s: Could not retrieve jets %s!", GetName(), fJetsName.Data()));
- }
- }
-
- return kTRUE;
-}
-
-//________________________________________________________________________
-Bool_t AliAnalysisTaskEmcalJet::IsJetTrack(AliEmcalJet* jet, Int_t itrack, Bool_t sorted) const
-{
- // Return true if track is in jet.
-
- for (Int_t i = 0; i < jet->GetNumberOfTracks(); ++i) {
- Int_t ijettrack = jet->TrackAt(i);
- if (sorted && ijettrack > itrack)
- return kFALSE;
- if (ijettrack == itrack)
- return kTRUE;
- }
- return kFALSE;
-}
-
-//________________________________________________________________________
-Bool_t AliAnalysisTaskEmcalJet::IsJetCluster(AliEmcalJet* jet, Int_t iclus, Bool_t sorted) const
-{
- // Return true if cluster is in jet.
-
- for (Int_t i = 0; i < jet->GetNumberOfClusters(); ++i) {
- Int_t ijetclus = jet->ClusterAt(i);
- if (sorted && ijetclus > iclus)
- return kFALSE;
- if (ijetclus == iclus)
- return kTRUE;
- }
- return kFALSE;
-}
-
//________________________________________________________________________
Bool_t AliAnalysisTaskEmcalJet::AcceptBiasJet(AliEmcalJet *jet) const
{
+ // Accept jet with a bias.
+
if (jet->MaxTrackPt() < fPtBiasJetTrack && (fAnaType == kTPC || jet->MaxClusterPt() < fPtBiasJetClus))
return kFALSE;
else
{
// Init the analysis.
+ AliAnalysisTaskEmcal::Init();
+
if (fAnaType == kTPC) {
SetEtaLimits(-0.9 + fJetRadius, 0.9 - fJetRadius);
SetPhiLimits(-10, 10);
} else {
AliWarning(Form("%s: Analysis type not recognized! Assuming kTPC!", GetName()));
SetAnaType(kTPC);
- Init();
}
SetInitialized();
}
//________________________________________________________________________
-void AliAnalysisTaskEmcalJet::UserExec(Option_t *)
+Bool_t AliAnalysisTaskEmcalJet::IsJetCluster(AliEmcalJet* jet, Int_t iclus, Bool_t sorted) const
+{
+ // Return true if cluster is in jet.
+
+ for (Int_t i = 0; i < jet->GetNumberOfClusters(); ++i) {
+ Int_t ijetclus = jet->ClusterAt(i);
+ if (sorted && ijetclus > iclus)
+ return kFALSE;
+ if (ijetclus == iclus)
+ return kTRUE;
+ }
+ return kFALSE;
+}
+
+//________________________________________________________________________
+Bool_t AliAnalysisTaskEmcalJet::IsJetTrack(AliEmcalJet* jet, Int_t itrack, Bool_t sorted) const
+{
+ // Return true if track is in jet.
+
+ for (Int_t i = 0; i < jet->GetNumberOfTracks(); ++i) {
+ Int_t ijettrack = jet->TrackAt(i);
+ if (sorted && ijettrack > itrack)
+ return kFALSE;
+ if (ijettrack == itrack)
+ return kTRUE;
+ }
+ return kFALSE;
+}
+
+//________________________________________________________________________
+Bool_t AliAnalysisTaskEmcalJet::RetrieveEventObjects()
{
- if (!fInitialized)
- Init();
+ // Retrieve objects from event.
- AliAnalysisTaskEmcal::UserExec("");
+ if (!AliAnalysisTaskEmcal::RetrieveEventObjects())
+ return kFALSE;
+
+ if (!fJetsName.IsNull() && !fJets) {
+ fJets = dynamic_cast<TClonesArray*>(InputEvent()->FindListObject(fJetsName));
+ if (!fJets) {
+ AliWarning(Form("%s: Could not retrieve jets %s!", GetName(), fJetsName.Data()));
+ }
+ }
+
+ return kTRUE;
}
AliAnalysisTaskEmcalJet(const char *name, Bool_t histo=kFALSE);
virtual ~AliAnalysisTaskEmcalJet();
- virtual void UserExec(Option_t *option);
- virtual void Init();
-
void SetEtaLimits(Float_t min, Float_t max) { fMinEta = min, fMaxEta = max ; }
- void SetInitialized(Bool_t ini = kTRUE) { fInitialized = ini ; }
void SetJetAreaCut(Float_t cut) { fJetAreaCut = cut ; }
void SetJetPtCut(Float_t cut) { fJetPtCut = cut ; }
void SetJetRadius(Float_t r) { fJetRadius = r ; }
void SetPtBiasJetTrack(Float_t b) { fPtBiasJetTrack = b ; }
protected:
-
Bool_t AcceptJet(AliEmcalJet* jet, Bool_t bias = kTRUE, Bool_t upCut = kTRUE) const;
Bool_t AcceptBiasJet(AliEmcalJet* jet) const;
+ void Init();
Bool_t IsJetTrack(AliEmcalJet* jet, Int_t itrack, Bool_t sorted = kTRUE) const;
Bool_t IsJetCluster(AliEmcalJet* jet, Int_t iclus, Bool_t sorted = kTRUE) const;
-
- virtual Bool_t RetrieveEventObjects();
+ Bool_t RetrieveEventObjects();
Float_t fJetRadius; // jet radius
TString fJetsName; // name of jet collection