////////////////////////////////////////////////////////////////////////
+AliAnalysisTaskJets::AliAnalysisTaskJets():
+ fDebug(0),
+ fJetFinder(0x0),
+ fChain(0x0),
+ fESD(0x0)
+{
+ // Default constructor
+}
+
AliAnalysisTaskJets::AliAnalysisTaskJets(const char* name):
AliAnalysisTask(name, "AnalysisTaskJets"),
fDebug(0),
{
// Default constructor
DefineInput (0, TChain::Class());
- DefineOutput(0, TH1::Class());
+// DefineOutput(0, TTree::Class());
}
-void AliAnalysisTaskJets::ConnectInputData(Option_t */*option*/)
+void AliAnalysisTaskJets::Init()
{
-// Initialisation
-//
- if (fDebug > 1) printf("AnalysisJets::Init() \n");
-
+ // Initialization
+ if (fDebug > 1) printf("AnalysisTaskJets::Init() \n");
// Call configuration file
gROOT->LoadMacro("ConfigJetAnalysis.C");
fJetFinder = (AliJetFinder*) gInterpreter->ProcessLine("ConfigJetAnalysis()");
// Initialise Jet Analysis
fJetFinder->Init();
+}
+
+void AliAnalysisTaskJets::ConnectInputData(Option_t */*option*/)
+{
+// Connect the input data
+//
+ if (fDebug > 1) printf("AnalysisTaskJets::ConnectInputData() \n");
fChain = (TChain*)GetInputData(0);
fJetFinder->ConnectTree(fChain);
fJetFinder->WriteHeaders();
{
// Terminate analysis
//
- printf("AnalysisJets: Terminate() \n");
+ if (fDebug > 1) printf("AnalysisJets: Terminate() \n");
if (fJetFinder) fJetFinder->FinishRun();
}
class AliAnalysisTaskJets : public AliAnalysisTask
{
- public:
- AliAnalysisTaskJets(const char* name);
- virtual ~AliAnalysisTaskJets() {;}
- // Implementation of interface methods
- virtual void ConnectInputData(Option_t *option = "");
- virtual void Exec(Option_t *option);
- virtual void Terminate(Option_t *option);
- virtual void SetDebugLevel(Int_t level) {fDebug = level;}
-
+ public:
+ AliAnalysisTaskJets();
+ AliAnalysisTaskJets(const char* name);
+ virtual ~AliAnalysisTaskJets() {;}
+ // Implementation of interface methods
+ virtual void ConnectInputData(Option_t *option = "");
+ virtual void Init();
+ virtual void Exec(Option_t *option);
+ virtual void Terminate(Option_t *option);
+ virtual void SetDebugLevel(Int_t level) {fDebug = level;}
+
private:
- Int_t fDebug; // Debug flag
- AliJetFinder* fJetFinder; // Pointer to the jet finder
- TChain* fChain; //! chained files
- AliESD* fESD; //! ESD
-
- ClassDef(AliAnalysisTaskJets, 1); // Analysis task for standard jet analysis
+ Int_t fDebug; // Debug flag
+ AliJetFinder* fJetFinder; // Pointer to the jet finder
+ TChain* fChain; //! chained files
+ AliESD* fESD; //! ESD
+
+ ClassDef(AliAnalysisTaskJets, 1); // Analysis task for standard jet analysis
};
#endif
int nMax = fChain->GetEntries();
- printf("\nTotal number of events in chain= %d \n",nMax);
+ printf("\n AliJetESDReader: Total number of events in chain= %d \n",nMax);
// set number of events in header
if (fReaderHeader->GetLastEvent() == -1)
fChain->SetBranchAddress("ESD", &fESD);
Int_t nMax = fChain->GetEntries();
- printf("\nTotal number of events in chain= %5d \n", nMax);
+ printf("\n AliJetESDReader: Total number of events in chain= %5d \n", nMax);
// set number of events in header
if (fReaderHeader->GetLastEvent() == -1)
fReaderHeader->SetLastEvent(nMax);
void SetEMCALGeometry();
void InitParameters();
protected:
- AliJetDummyGeo *fGeom; //!EMCAL Geometry
- TChain *fChain; // chain for reconstructed tracks
- AliESD *fESD; // pointer to esd
- AliJetHadronCorrectionv1 *fHadCorr; // Pointer to Hadron Correction Object
- AliJetGrid *fTpcGrid; // Pointer to grid object
- AliJetGrid *fEmcalGrid; // Pointer to grid object
+ AliJetDummyGeo *fGeom; //! EMCAL Geometry
+ TChain *fChain; //!chain for reconstructed tracks
+ AliESD *fESD; //!pointer to esd
+ AliJetHadronCorrectionv1 *fHadCorr; //! Pointer to Hadron Correction Object
+ AliJetGrid *fTpcGrid; //! Pointer to grid object
+ AliJetGrid *fEmcalGrid; //! Pointer to grid object
Float_t fPtCut; // Pt cut for tracks to minimise background contribution
Int_t fHCorrection; // Hadron correction flag
Int_t fNumUnits; // Number of units in the unit object array
ClassImp(AliJetFinder)
AliJetFinder::AliJetFinder():
+ fJetT(0),
fPlotMode(kFALSE),
fJets(0),
fGenJets(0),
fLeading(0),
fReader(0x0),
+ fHeader(0x0),
fPlots(0x0),
fOut(0x0)
fOut->cd();
char hname[30];
sprintf(hname,"TreeJ%d",i);
- TTree* jetT = new TTree(hname,"AliJet");
- jetT->Branch("FoundJet",&fJets,1000);
- jetT->Branch("GenJet",&fGenJets,1000);
- jetT->Branch("LeadingPart",&fLeading,1000);
- jetT->Fill();
- jetT->Write(hname);
- delete jetT;
+ fJetT = new TTree(hname,"AliJet");
+ fJetT->Branch("FoundJet",&fJets,1000);
+ fJetT->Branch("GenJet",&fGenJets,1000);
+ fJetT->Branch("LeadingPart",&fLeading,1000);
+ fJetT->Fill();
+ fJetT->Write(hname);
+ delete fJetT;
}
////////////////////////////////////////////////////////////////////////
protected:
AliJetFinder(const AliJetFinder& rJetFinder);
AliJetFinder& operator = (const AliJetFinder& rhsf);
-
- Bool_t fPlotMode; // do you want control plots?
- AliJet* fJets; // pointer to jet class
- AliJet* fGenJets; // pointer to generated jets
- AliLeading* fLeading; // pointer to leading particle data
- AliJetReader* fReader; // pointer to reader
- AliJetHeader* fHeader; // pointer to header
- AliJetControlPlots* fPlots; // pointer to control plots
- TFile* fOut; // output file
+ TTree* fJetT; //! pointer to jet tree
+ Bool_t fPlotMode; //! do you want control plots?
+ AliJet* fJets; //! pointer to jet class
+ AliJet* fGenJets; //! pointer to generated jets
+ AliLeading* fLeading; //! pointer to leading particle data
+ AliJetReader* fReader; // pointer to reader
+ AliJetHeader* fHeader; // pointer to header
+ AliJetControlPlots* fPlots; //! pointer to control plots
+ TFile* fOut; //! output file
ClassDef(AliJetFinder,2)
};
AliUA1JetFinder& operator = (const AliUA1JetFinder& ruaf);
AliUA1JetHeader* fHeader; // pointer to jet header
- TH2F * fLego; //! Lego Histo
+ TH2F * fLego; // Lego Histo
ClassDef(AliUA1JetFinder,1)
};