* Author: The ALICE Off-line Project. *
* Contributors are mentioned in the code where appropriate. *
* *
- * Permission to use, copy, modify and distribute this software and its *
+ * Permission to usec, copy, modify and distribute this software and its *
* documentation strictly for non-commercial purposes is hereby granted *
* without fee, provided that the above copyright notice appears in all *
* copies and that both the copyright notice and this permission notice *
// * Go to the next highest PT particle (not already included in a jet) and include in the "jet" all particles
// (not already included in a jet) within the radius R =0.7.
// * Continue until all particles are in a "jet".
-
+ if (fDebug) { printf("AliCDJetfinder::FindJets() %d \n", __LINE__ ); }
AliCdfJetHeader *header = (AliCdfJetHeader*)fHeader;
if (header)
ComputeConesWeight();
-if (fAODwrite) { cout << "Writing AOD" << endl ; WriteJets(); }
+if (fAODwrite) {
+ if(fDebug)cout << "Writing AOD" << endl ;
+ WriteJets();
+ }
if (fAnalyseJets) AnalizeJets();
if (fDebug) jet.Print("");
if (fFromAod && fAODtracksWrite)
- {
+ {
for ( Int_t jetTrack = 0; jetTrack < fNPart; jetTrack++ )
- { if ( fVectParticle[jetTrack]->njet == jetnr ) { jet.AddTrack(fRefArr->At(jetTrack)) ; } }
- }
+ { if ( fVectParticle[jetTrack]->njet == jetnr ) { jet.AddTrack(fRefArr->At(jetTrack)) ; } }
+ }
// tracks REFs written in AOD
AddJet(jet);
}
// If the correction is > jet energy px = py = pz = e = 0
if (corrJets[j].px() == 0 && corrJets[j].py() == 0 && corrJets[j].pz() == 0 && corrJets[j].E() == 0) continue;
- cout<<"********************************** Reconstructed jet(s) (non corrected)"<<endl;
- cout<<"Jet number "<<j+1<<" : "<<"Rapidity : "<<jets[j].rap()<<" Phi : "<<jets[j].phi()<<" pT : "<<jets[j].perp()<<" area : "<<clustSeq.area(jets[j])<<endl;
-// cout<<"px = "<<jets[j].px()<<endl;
-// cout<<"py = "<<jets[j].py()<<endl;
-// cout<<"pz = "<<jets[j].pz()<<endl;
- cout<<"e = "<<jets[j].E()<<endl;
-
- cout<<"********************************** Corrected jet(s)"<<endl;
- cout<<"Jet number "<<j+1<<" : "<<"Rapidity : "<<corrJets[j].rap()<<" Phi : "<<corrJets[j].phi()<<" pT : "<<corrJets[j].perp()<<endl;
-// cout<<"px = "<<corrJets[j].px()<<endl;
-// cout<<"py = "<<corrJets[j].py()<<endl;
-// cout<<"pz = "<<corrJets[j].pz()<<endl;
- cout<<"e = "<<corrJets[j].E()<<endl;
-
+ if(debug){
+ cout<<"********************************** Reconstructed jet(s) (non corrected)"<<endl;
+ cout<<"Jet number "<<j+1<<" : "<<"Rapidity : "<<jets[j].rap()<<" Phi : "<<jets[j].phi()<<" pT : "<<jets[j].perp()<<" area : "<<clustSeq.area(jets[j])<<endl;
+ // cout<<"px = "<<jets[j].px()<<endl;
+ // cout<<"py = "<<jets[j].py()<<endl;
+ // cout<<"pz = "<<jets[j].pz()<<endl;
+ cout<<"e = "<<jets[j].E()<<endl;
+
+ cout<<"********************************** Corrected jet(s)"<<endl;
+ cout<<"Jet number "<<j+1<<" : "<<"Rapidity : "<<corrJets[j].rap()<<" Phi : "<<corrJets[j].phi()<<" pT : "<<corrJets[j].perp()<<endl;
+ // cout<<"px = "<<corrJets[j].px()<<endl;
+ // cout<<"py = "<<corrJets[j].py()<<endl;
+ // cout<<"pz = "<<corrJets[j].pz()<<endl;
+ cout<<"e = "<<corrJets[j].E()<<endl;
+ }
// Go to write AOD info
AliAODJet aodjet (corrJets[j].px(), corrJets[j].py(), corrJets[j].pz(), corrJets[j].E());
if(debug) aodjet.Print("");
for (size_t k = 0; k < jets.size(); k++)
{
- cout<<"********************************** Reconstructed jet(s) (non corrected)"<<endl;
- cout<<"Jet number "<<k+1<<" : "<<"Rapidity : "<<jets[k].rap()<<" Phi : "<<jets[k].phi()<<" pT : "<<jets[k].perp()<<endl;
-// cout<<"px = "<<jets[k].px()<<endl;
-// cout<<"py = "<<jets[k].py()<<endl;
-// cout<<"pz = "<<jets[k].pz()<<endl;
- cout<<"e = "<<jets[k].E()<<endl;
-
+ if(debug){
+ cout<<"********************************** Reconstructed jet(s) (non corrected)"<<endl;
+ cout<<"Jet number "<<k+1<<" : "<<"Rapidity : "<<jets[k].rap()<<" Phi : "<<jets[k].phi()<<" pT : "<<jets[k].perp()<<endl;
+ // cout<<"px = "<<jets[k].px()<<endl;
+ // cout<<"py = "<<jets[k].py()<<endl;
+ // cout<<"pz = "<<jets[k].pz()<<endl;
+ cout<<"e = "<<jets[k].E()<<endl;
+ }
// Go to write AOD info
AliAODJet aodjet (jets[k].px(), jets[k].py(), jets[k].pz(), jets[k].E());
if(debug) aodjet.Print("");