Int_t selec=-1;
Double_t ptmax=-10;
Double_t areaj=0;
- Double_t phij=0;
+ // Double_t phij=0;
Double_t etaj=0;
Double_t ptj=0;
Double_t ptcorrj=0;
for(Int_t i=0; i<fListJets[0]->GetEntries(); ++i){
AliAODJet* jetj = (AliAODJet*)(fListJets[0]->At(i));
etaj = jetj->Eta();
- phij = jetj->Phi();
+ // phij = jetj->Phi();
ptj = jetj->Pt();
if(ptj==0) continue;
areaj = jetj->EffectiveAreaCharged();
if(j==selec) continue;
AliAODJet* jetb = (AliAODJet*)(fListJets[0]->At(j));
etaj = jetb->Eta();
- phij = jetb->Phi();
+ // phij = jetb->Phi();
ptj = jetb->Pt();
if(ptj<=0) continue;
areaj = jetb->EffectiveAreaCharged();
Float_t ptMax = 0.;
Float_t etaMax = 0.;
Float_t phiMax = 0.;
- Int_t iMax = -1;
+ // Int_t iMax = -1;
//1st loop over all tracks
for(int it = 0;it<nT;++it){
AliVParticle *track = (AliVParticle*)ParticleList.At(it);
// max pt
if (pt > ptMax) {
ptMax = pt;
- iMax = it;
+ // iMax = it;
etaMax = deta;
phiMax = dphi;
}
Double_t ptmax=-10;
Double_t dphi=0;
- Double_t dif=0;
+ // Double_t dif=0;
Int_t iCount=0;
for(int it = 0;it < aod->GetNumberOfTracks();++it){
AliAODTrack *tr = aod->GetTrack(it);
if(TMath::Abs(dphi)<TMath::Pi()-0.2) continue;
if(tr->Pt()>ptmax){ ptmax=tr->Pt();
index=iCount-1;
- dif=dphi; }}
+ // dif=dphi;
+ }}
return index;
Double_t phibig=0.;
Double_t etasmall=0;
Double_t ptsmall=0;
- Double_t areasmall=0;
+ // Double_t areasmall=0;
Double_t phismall=0.;
Int_t iCount=0;
Int_t trigJet=-1;
Int_t trigBBTrack=-1;
- Int_t trigInTrack=-1;
+ // Int_t trigInTrack=-1;
fRPAngle = aod->GetHeader()->GetEventplane();
if(fHardest==0 || fHardest==1){
if(iCount==0){
trigJet=i;
trigBBTrack=nT;
- trigInTrack=ippt;
+ // trigInTrack=ippt;
iCount=iCount+1;}
etasmall = jetsmall->Eta();
phismall = jetsmall->Phi();
ptsmall = jetsmall->Pt();
- areasmall = jetsmall->EffectiveAreaCharged();
+ // areasmall = jetsmall->EffectiveAreaCharged();
Double_t tmpDeltaR=(phismall-phibig)*(phismall-phibig)+(etasmall-etabig)*(etasmall-etabig);
tmpDeltaR=TMath::Sqrt(tmpDeltaR);
//Fraction in the jet core
Int_t index=-1;
Double_t ptmax=-10;
Double_t dphi=0;
- Double_t dif=0;
+ // Double_t dif=0;
Int_t iCount=0;
for(int it = 0;it < aod->GetNumberOfTracks();++it){
AliAODTrack *tr = aod->GetTrack(it);
if(TMath::Abs(dphi)<TMath::Pi()-0.6) continue;
if(tr->Pt()>ptmax){ ptmax=tr->Pt();
index=iCount-1;
- dif=dphi; }}
+ // dif=dphi;
+ }}
return index;
if(fDebug > 2) printf("AliAnalysisTaskJetProperties::FillJetProperties() \n");
for(Int_t iJet=0; iJet < jetlist->GetEntries(); iJet++){
- Float_t JetPt;Float_t JetPhi;Float_t JetEta;Float_t JetE;
+ Float_t JetPt;Float_t JetPhi;Float_t JetEta; // Float_t JetE;
AliAODJet *jet = dynamic_cast<AliAODJet*>(jetlist->At(iJet));
if(!jet)continue;
JetEta = jet->Eta();
JetPhi = jet->Phi();
JetPt = jet->Pt();
- JetE = jet->E();
+ // JetE = jet->E();
fh2EtaJet ->Fill(JetPt,JetEta);
fh2PhiJet ->Fill(JetPt,JetPhi);
fh2PtJet ->Fill(JetPt,JetPt);
if(fTrackType==kTrackUndef)continue;
Float_t TrackEta=-99.0; Float_t TrackPt=-99.0; Float_t TrackPhi=-99.0;
Float_t FF=-99.0; Float_t DelEta=-99.0; Float_t DelPhi=-99.0;
- Float_t DelR=-99.0; Float_t AreaJ=-99.0; Float_t Ksi=-99.0;
+ Float_t DelR=-99.0; // Float_t AreaJ=-99.0;
+ Float_t Ksi=-99.0;
if(fTrackType==kTrackAOD){
AliAODTrack *trackaod = dynamic_cast<AliAODTrack*>(fTrackListJet->At(j));
if(!trackaod)continue;
DelPhi = TMath::Abs(JetPhi - TrackPhi);
if(DelPhi>TMath::Pi())DelPhi = TMath::Abs(DelPhi-TMath::TwoPi());
DelR = TMath::Sqrt(DelEta*DelEta + DelPhi*DelPhi);
- AreaJ = TMath::Pi()*DelR*DelR;
+ // AreaJ = TMath::Pi()*DelR*DelR;
fh2EtaTrack ->Fill(JetPt,TrackEta);
fh2PhiTrack ->Fill(JetPt,TrackPhi);
fh2PtTrack ->Fill(JetPt,TrackPt);
if(incrementJetPt) fh1JetPt->Fill(jetPt);
//fh3TrackPt->Fill(jetPt,invM,trackPt);//Fill(x,y,z)
- invM = 0;
- Double_t z = 0.;
- if(jetPt>0) z = trackPt / jetPt;
+ // invM = 0;
+ // Double_t z = 0.;
+ // if(jetPt>0) z = trackPt / jetPt;
// Double_t xi = 0;
//if(z>0) xi = TMath::Log(1/z);