for(UInt_t ij=0, jetCount=0; ij<jets_incl.size(); ++ij) {
if (jets_incl[ij].perp()<1)
continue;
-#if 0
- AliAODJet *jet = new ((*fJets)[jetCount])
- AliAODJet(jets_incl[ij].px(), jets_incl[ij].py(), jets_incl[ij].pz(), jets_incl[ij].E());
- jet->SetEffArea(fjw.GetJetArea(ij),0);
+ //cout << jets_incl[ij].perp() << " " << jets_incl[ij].eta() << " " << jets_incl[ij].phi() << " " << jets_incl[ij].m() << endl;
+ AliESDJet *jet = new ((*fJets)[jetCount])
+ AliESDJet(jets_incl[ij].perp(), jets_incl[ij].eta(), jets_incl[ij].phi(), jets_incl[ij].m());
+ jet->SetArea(fjw.GetJetArea(ij));
+#if 0 /*tobedone*/
jet->GetRefTracks()->Clear();
vector<fastjet::PseudoJet> constituents = fjw.GetJetConstituents(ij);
Double_t neutralE = 0;
void Terminate(Option_t *option);
void SetAlgo(Int_t a) { fAlgo = a; }
- void SetCaloName(const char *n) { fCaloName = n; }
+ void SetClusName(const char *n) { fCaloName = n; }
void SetHadCorr(Double_t c) { fHadCorr = c; }
void SetJetsName(const char *n) { fJetsName = n; }
void SetTracksName(const char *n) { fTracksName = n; }
protected:
void FindJets(TObjArray *tracks, TObjArray *clus, Int_t algo, Double_t radius);
- TString fTracksName; // name of track collection (if "" use branch)
- TString fCaloName; // name of calo collection
+ TString fTracksName; // name of track collection
+ TString fCaloName; // name of calo cluster collection
TString fJetsName; // name of jet collection
Int_t fAlgo; // algo (0==kt, 1==antikt)
Double_t fRadius; // jet radius