//do generic subtraction if requested
if(fDoGenericSubtraction) {
-#ifdef FASTJET_CONTRIB
+#ifdef FASTJET_VERSION
std::vector<fastjet::contrib::GenericSubtractorInfo> jetMassInfo = fjw.GetGenSubtractorInfoJetMass();
UInt_t n = (UInt_t)jetMassInfo.size();
if(n>ij && n>0) {
//run constituent subtractor if requested
if(fDoConstituentSubtraction) {
-#ifdef FASTJET_CONTRIB
+#ifdef FASTJET_VERSION
if(!fJetsSub) AliWarning(Form("No jet branch to write to for subtracted jets. fJetsSubName: %s",fJetsSubName.Data()));
else {
std::vector<fastjet::PseudoJet> jets_sub;
Double_t GetJetSubtractedPt (UInt_t idx) const;
virtual std::vector<double> GetSubtractedJetsPts(Double_t median_pt = -1, Bool_t sorted = kFALSE);
Bool_t GetLegacyMode() { return fLegacyMode; }
-#ifdef FASTJET_CONTRIB
+#ifdef FASTJET_VERSION
const std::vector<fastjet::contrib::GenericSubtractorInfo> GetGenSubtractorInfoJetMass() const {return fGenSubtractorInfoJetMass;}
const std::vector<fastjet::PseudoJet> GetConstituentSubtrJets() const { return fConstituentSubtrJets; }
#endif
Bool_t fUseArea4Vector; //!
#ifdef FASTJET_VERSION
fastjet::JetMedianBackgroundEstimator *fBkrdEstimator; //!
-#endif
-#ifdef FASTJET_CONTRIB
+ //from contrib package
fastjet::contrib::GenericSubtractor *fGenSubtractor; //!
std::vector<fastjet::contrib::GenericSubtractorInfo> fGenSubtractorInfoJetMass; //!
#endif
, fUseArea4Vector (kFALSE)
#ifdef FASTJET_VERSION
, fBkrdEstimator (0)
-#endif
-#ifdef FASTJET_CONTRIB
, fGenSubtractor (0)
, fGenSubtractorInfoJetMass ( )
#endif
delete fClustSeq;
#ifdef FASTJET_VERSION
if (fBkrdEstimator) delete fBkrdEstimator;
-#endif
-#ifdef FASTJET_CONTRIB
if (fGenSubtractor) delete fGenSubtractor;
#endif
}
delete fClustSeq; fClustSeq = 0;
#ifdef FASTJET_VERSION
if (fBkrdEstimator) delete fBkrdEstimator ; fBkrdEstimator = 0;
-#endif
-#ifdef FASTJET_CONTRIB
if (fGenSubtractor) delete fGenSubtractor ; fGenSubtractor = 0;
#endif
}
//_________________________________________________________________________________________________
Int_t AliFJWrapper::DoGenericSubtractionJetMass() {
//Do generic subtraction for jet mass
-#ifdef FASTJET_CONTRIB
+#ifdef FASTJET_VERSION
if(fUseExternalBkg) fGenSubtractor = new fj::contrib::GenericSubtractor(fRho,fRhom);
else fGenSubtractor = new fj::contrib::GenericSubtractor(fBkrdEstimator);
//_________________________________________________________________________________________________
Int_t AliFJWrapper::DoConstituentSubtraction() {
//Do constituent subtraction
-#ifdef FASTJET_CONTRIB
+#ifdef FASTJET_VERSION
fj::contrib::ConstituentSubtractor *subtractor;
if(fUseExternalBkg)
subtractor = new fj::contrib::ConstituentSubtractor(fRho,fRhom,kFALSE,kTRUE);
#include <fastjet/tools/JetMedianBackgroundEstimator.hh>
#include <fastjet/tools/BackgroundEstimatorBase.hh>
#include <fastjet/tools/Subtractor.hh>
-#ifdef FASTJET_CONTRIB
+//from contrib package
#include <fastjet/contrib/GenericSubtractor.hh>
#include <fastjet/contrib/ShapeWithComponents.hh>
#include <fastjet/contrib/ConstituentSubtractor.hh>
#endif
#endif
-#endif
#endif