/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
* See cxx source for full Copyright notice */
+
+/* $Id$ */
//---------------------------------------------------------------------
// SISCone (FastJet v2.3.4) finder algorithm interface
// Finder Header Class
// Author: swensy.jangal@ires.in2p3.fr
//---------------------------------------------------------------------
-
-
-#include "fastjet/AreaDefinition.hh"
-#include "fastjet/ClusterSequenceArea.hh"
+#ifndef __CINT__
#include "fastjet/JetDefinition.hh"
+#else
+namespace fastjet {
+ class Strategy;
+ class RecombinationScheme;
+}
+#endif
#include "AliJetHeader.h"
class AliSISConeJetHeader : public AliJetHeader
virtual ~AliSISConeJetHeader() { }
// Getters
- Bool_t GetBGMode() const {return fBGMode;}
- Bool_t GetDebug() const {return fDebug;}
+ Bool_t GetBGMode() const {return fBGMode;}
+ Int_t GetActiveAreaRepeats() const {return fActiveAreaRepeats;}
+ Int_t GetAreaTypeNumber() const {return fAreaTypeNumber;}
+ Int_t GetBGAlgorithm() const {return fBGAlgo;}
+ Int_t GetNPassMax() const {return fNPassMax;}
+ Int_t GetSplitMergeScale() const {return fSplitMergeScaleNumber;}
+ Double_t GetGhostEtaMax() const {return fGhostEtaMax;}
+ Double_t GetGhostArea() const {return fGhostArea;}
+ Double_t GetEffectiveRFact() const {return fEffectiveRFact;}
+ Double_t GetRapMax() const {return fRapMax;}
+ Double_t GetRapMin() const {return fRapMin;}
+ Double_t GetPhiMax() const {return fPhiMax;}
+ Double_t GetPhiMin() const {return fPhiMin;}
+ Double_t GetConeRadius() const {return fConeRadius;}
+ Double_t GetOverlapThreshold() const {return fOverlapThreshold;}
+ Double_t GetPtProtojetMin() const {return fPtProtoJetMin;}
+ Double_t GetRForRho() const {return fRRho;}
+ Double_t GetCaching() const {return fCaching;}
+ Double_t GetSplitMergeStoppingScale() const {return fSplitMergeStoppingScale;}
+ Double_t GetMinJetPt() const {return fMinJetPt;}
+ Double_t GetGridScatter() const {return fGridScatter;}
+ Double_t GetKtScatter() const {return fKtScatter;}
+ Double_t GetMeanGhostKt() const {return fMeanGhostKt;}
- Int_t GetActiveAreaRepeats() const {return fActiveAreaRepeats;}
- Int_t GetAreaTypeNumber() const {return fAreaTypeNumber;}
- Int_t GetBGAlgorithm() const {return fBGAlgo;}
- Int_t GetNPassMax() const {return fNPassMax;}
- Int_t GetSplitMergeScale() const {return fSplitMergeScaleNumber;}
+ // Setters
+ void SetBGAlgorithm(Int_t value) {fBGAlgo = value;}
+ void SetBGMode(Bool_t value) {fBGMode = value;}
+ void SetCaching(Bool_t value) {fCaching = value;}
+ void SetComment(TString com) {fComment=com;}
+ void SetComment(const char* com) {AliJetHeader::SetComment(com);}
+ void SetGhostEtaMax(Double_t f) {fGhostEtaMax = f;}
+ void SetGhostArea(Double_t f) {fGhostArea = f;}
+ void SetActiveAreaRepeats(Int_t f) {fActiveAreaRepeats =f;}
+ void SetAreaTypeNumber(Int_t f) {fAreaTypeNumber = f;}
+ void SetEffectiveRFact(Double_t value) {fEffectiveRFact = value;}
+ void SetConeRadius(Double_t value) {fConeRadius = value;}
+ void SetMinJetPt(Double_t value) {fMinJetPt = value;}
+ void SetNPassMax(Int_t value) {fNPassMax = value;}
+ void SetOverlapThreshold(Double_t value) {fOverlapThreshold = value;}
+ void SetPhiRange(Double_t fmin, Double_t fmax) {fPhiMin = fmin; fPhiMax = fmax;}
+ void SetPtProtojetMin(Double_t value) {fPtProtoJetMin = value;}
+ void SetRapRange(Double_t fmin, Double_t fmax) {fRapMin = fmin; fRapMax = fmax;}
+ void SetRForRho(Double_t value) {fRRho = value;}
+ void SetSplitMergeScale(Int_t value) {fSplitMergeScaleNumber = value;}
+ void SetSplitMergeStoppingScale(Double_t value) {fSplitMergeStoppingScale = value;}
+ void SetGridScatter(Double_t value) {fGridScatter = value;}
+ void SetKtScatter(Double_t value) {fKtScatter = value;}
+ void SetMeanGhostKt(Double_t value) {fMeanGhostKt = value;}
- Double_t GetGhostEtaMax() const {return fGhostEtaMax;}
- Double_t GetGhostArea() const {return fGhostArea;}
- Double_t GetEffectiveRFact() const {return fEffectiveRFact;}
- Double_t GetRapMax() const {return fRapMax;}
- Double_t GetRapMin() const {return fRapMin;}
- Double_t GetPhiMax() const {return fPhiMax;}
- Double_t GetPhiMin() const {return fPhiMin;}
- Double_t GetConeRadius() const {return fConeRadius;}
- Double_t GetOverlapThreshold() const {return fOverlapThreshold;}
- Double_t GetPtProtojetMin() const {return fPtProtoJetMin;}
- Double_t GetCaching() const {return fCaching;}
- Double_t GetSplitMergeStoppingScale() const {return fSplitMergeStoppingScale;}
- Double_t GetMinJetPt() const {return fMinJetPt;}
- Double_t GetGridScatter() const {return fGridScatter;}
- Double_t GetKtScatter() const {return fKtScatter;}
- Double_t GetMeanGhostKt() const {return fMeanGhostKt;}
+ // Added for background stuff
+ fastjet::Strategy GetStrategy() const {return fStrategy;}
+ fastjet::RecombinationScheme GetRecombScheme() const {return fRecombScheme;}
+ Double_t GetRparamBkg() const {return fRparamBkg;}
+ Bool_t Use4VectorArea() const {return fkUse4VectorArea;}
- // Setters
- void SetBGAlgorithm(Int_t value) {fBGAlgo = value;}
- void SetBGMode(Bool_t value) {fBGMode = value;}
- void SetDebug(Bool_t value) {fDebug = value;}
- void SetCaching(Bool_t value) {fCaching = value;}
- void SetComment(TString com) {fComment=com;}
- void SetComment(const char* com) {AliJetHeader::SetComment(com);}
- void SetGhostEtaMax(Double_t f) {fGhostEtaMax = f;}
- void SetGhostArea(Double_t f) {fGhostArea = f;}
- void SetActiveAreaRepeats(Int_t f) {fActiveAreaRepeats =f;}
- void SetAreaTypeNumber(Int_t f) {fAreaTypeNumber = f;}
- void SetEffectiveRFact(Double_t value) {fEffectiveRFact = value;}
- void SetConeRadius(Double_t value) {fConeRadius = value;}
- void SetMinJetPt(Double_t value) {fMinJetPt = value;}
- void SetNPassMax(Int_t value) {fNPassMax = value;}
- void SetOverlapThreshold(Double_t value) {fOverlapThreshold = value;}
- void SetPhiRange(Double_t fmin, Double_t fmax) {fPhiMin = fmin; fPhiMax = fmax;}
- void SetPtProtojetMin(Double_t value) {fPtProtoJetMin = value;}
- void SetRapRange(Double_t fmin, Double_t fmax) {fRapMin = fmin; fRapMax = fmax;}
- void SetSplitMergeScale(Int_t value) {fSplitMergeScaleNumber = value;}
- void SetSplitMergeStoppingScale(Double_t value) {fSplitMergeStoppingScale = value;}
- void SetGridScatter(Double_t value) {fGridScatter = value;}
- void SetKtScatter(Double_t value) {fKtScatter = value;}
- void SetMeanGhostKt(Double_t value) {fMeanGhostKt = value;}
+ void SetUse4VectorArea() {fkUse4VectorArea = kTRUE;}
+ void SetStrategy(fastjet::Strategy f) {fStrategy = f;}
+ void SetRecombScheme(fastjet::RecombinationScheme f) {fRecombScheme = f;}
+ void SetRparamBkg(Double_t f) {fRparamBkg = f;}
// others
- void PrintParameters() const;
+ void PrintParameters() const;
protected:
- Int_t fActiveAreaRepeats; // How many times do you want to caculate active areas?
- Int_t fAreaTypeNumber; // Kind of area
- Int_t fBGAlgo; // Algorithm for rho calculus
- Bool_t fBGMode; // Do we subtract BG or not?
- Bool_t fCaching; // Do we record found cones for this set of data?
- Double_t fConeRadius; // Cone radius
- Bool_t fDebug; // Debug option
- Double_t fEffectiveRFact; // Radius for Voronoi diagram
- Double_t fGhostEtaMax; // Maximum eta in which a ghost can be generated
- Double_t fGhostArea; // Area of one ghost
- Double_t fGridScatter; // fractional random fluctuations of the position of the ghosts on the y-phi grid
- Double_t fKtScatter; // fractional random fluctuations of the tranverse momentum of the ghosts on the y-phi grid
- Double_t fMeanGhostKt; // average transverse momentum of the ghosts.
- Double_t fMinJetPt; // Jet minimum energy
- Int_t fNPassMax; // maximum number of passes
- Double_t fOverlapThreshold; // overlap parameter
- Double_t fPhiMax, fPhiMin; // Phi range
- Double_t fPtProtoJetMin; // pT min of protojets
- Double_t fRapMax, fRapMin; // Eta range
- Int_t fSplitMergeScaleNumber; // Kind of recombination in split/merge procedure, there's only one
- Double_t fSplitMergeStoppingScale; // Stopping scale for split/merge procedure in case of area calculus
+ Int_t fActiveAreaRepeats; // How many times do you want to caculate active areas?
+ Int_t fAreaTypeNumber; // Kind of area
+ Int_t fBGAlgo; // Algorithm for rho calculus
+ Bool_t fCaching; // Do we record found cones for this set of data?
+ Double_t fConeRadius; // Cone radius
+ Double_t fEffectiveRFact; // Radius for Voronoi diagram
+ Double_t fGhostEtaMax; // Maximum eta in which a ghost can be generated
+ Double_t fGhostArea; // Area of one ghost
+ Double_t fGridScatter; // fractional random fluctuations of the position of the ghosts on the y-phi grid
+ Double_t fKtScatter; // fractional random fluctuations of the tranverse momentum of the ghosts on the y-phi grid
+ Double_t fMeanGhostKt; // average transverse momentum of the ghosts.
+ Double_t fMinJetPt; // Jet minimum energy
+ Int_t fNPassMax; // maximum number of passes
+ Double_t fOverlapThreshold; // overlap parameter
+ Double_t fPhiMax, fPhiMin; // Phi range
+ Double_t fPtProtoJetMin; // pT min of protojets
+ Double_t fRapMax, fRapMin; // Eta range
+ Double_t fRRho; // Radius to determine rho
+ Int_t fSplitMergeScaleNumber; // Kind of recombination in split/merge procedure, there's only one
+ Double_t fSplitMergeStoppingScale; // Stopping scale for split/merge procedure in case of area calculus
+
+ // Added for background
+ Double_t fRparamBkg; //R param for bkg calculation
+ fastjet::Strategy fStrategy; // fastjet::Best;
+ fastjet::RecombinationScheme fRecombScheme; // fastjet::BIpt_scheme;
+ Bool_t fkUse4VectorArea; // Toggle use of 4-vector area
- ClassDef(AliSISConeJetHeader,2)
+ ClassDef(AliSISConeJetHeader,5) // SISCONE header class
};
#endif