* provided "as is" without express or implied warranty. *
**************************************************************************/
-//---------------------------------------------------------------------
+//------------------------------------------------------------------------
// Jet reader base class
// manages the reading of input for jet algorithms
-// Author: jgcn@mda.cinvestav.mx
-//---------------------------------------------------------------------
+// Authors: jgcn@mda.cinvestav.mx
+// magali.estienne@IReS.in2p3.fr
+//------------------------------------------------------------------------
+// root
#include <TClonesArray.h>
-
+#include <TRefArray.h>
+#include "TTask.h"
+//AliRoot
#include "AliJetReader.h"
#include "AliJetReaderHeader.h"
-#include "AliESD.h"
+#include "AliESDEvent.h"
#include "AliHeader.h"
+#include "AliJetFillUnitArrayTracks.h"
+#include "AliJetFillUnitArrayEMCalDigits.h"
+#include "AliJetUnitArray.h"
+#include "AliJetHadronCorrectionv1.h"
ClassImp(AliJetReader)
////////////////////////////////////////////////////////////////////////
-
-AliJetReader::AliJetReader()
-{
+AliJetReader::AliJetReader():
// Constructor
- fChain = 0;
- fChainMC = 0;
- fESD = 0;
- fMomentumArray = new TClonesArray("TLorentzVector",2000);
- fArrayMC = 0;
- fAliHeader = 0;
+ fChain(0),
+ fMomentumArray(new TClonesArray("TLorentzVector",2000)),
+ fArrayMC(0),
+ fFillUnitArray(new TTask("fillUnitArray","Fill unit array jet finder")),
+ fESD(0),
+ fReaderHeader(0),
+ fSignalFlag(0),
+ fCutFlag(0),
+ fUnitArray(new TClonesArray("AliJetUnitArray",60000)),
+ fRefArray(new TRefArray()),
+ fUnitArrayNoCuts(new TClonesArray("AliJetUnitArray",60000)),
+ fArrayInitialised(0),
+ fFillUAFromTracks(new AliJetFillUnitArrayTracks()),
+ fFillUAFromEMCalDigits(new AliJetFillUnitArrayEMCalDigits()),
+ fNumCandidate(0),
+ fNumCandidateCut(0)
+{
+ // Default constructor
+ fSignalFlag = TArrayI();
+ fCutFlag = TArrayI();
}
////////////////////////////////////////////////////////////////////////
AliJetReader::~AliJetReader()
{
// Destructor
- delete fChain;
- delete fChainMC;
- delete fESD;
- delete fAliHeader;
if (fMomentumArray) {
fMomentumArray->Delete();
delete fMomentumArray;
}
+
+ if (fUnitArray) {
+ fUnitArray->Delete();
+ delete fUnitArray;
+ }
+
+ if (fUnitArrayNoCuts) {
+ fUnitArrayNoCuts->Delete();
+ delete fUnitArrayNoCuts;
+ }
+
+ if (fFillUnitArray) {
+ fFillUnitArray->Delete();
+ delete fFillUnitArray;
+ }
delete fArrayMC;
+
}
{
if (fMomentumArray) fMomentumArray->Clear();
+ if (fFillUnitArray) fFillUnitArray->Clear();
}