X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PWGJE%2FAliAnalysisTaskAj.cxx;h=165c6447ef3544b0f9e1e11fb94b7424a32622d9;hb=489951cd215c791ea1c54717475c98af0303b124;hp=8a456e66f6372d03fe8e3271cf36b6b41d520c15;hpb=5bd732d46e1bfa7c4c8bf6840234c8d627ace49f;p=u%2Fmrichter%2FAliRoot.git diff --git a/PWGJE/AliAnalysisTaskAj.cxx b/PWGJE/AliAnalysisTaskAj.cxx index 8a456e66f63..165c6447ef3 100644 --- a/PWGJE/AliAnalysisTaskAj.cxx +++ b/PWGJE/AliAnalysisTaskAj.cxx @@ -53,6 +53,9 @@ #include "AliAnalysisTaskAj.h" +using std::cout; +using std::endl; + ClassImp(AliAnalysisTaskAj) AliAnalysisTaskAj::AliAnalysisTaskAj() : @@ -479,11 +482,11 @@ void AliAnalysisTaskAj::UserExec(Option_t *) else centValue=aod->GetHeader()->GetCentrality(); if(fDebug) printf("centrality: %f\n", centValue); - if (centValue < fCentMin || centValue > fCentMax){ - fHistEvtSelection->Fill(4); - PostData(1, fOutputList); - return; - } + // if (centValue < fCentMin || centValue > fCentMax){ + // fHistEvtSelection->Fill(4); + // PostData(1, fOutputList); + // return; + // } fHistEvtSelection->Fill(0); @@ -566,8 +569,8 @@ void AliAnalysisTaskAj::UserExec(Option_t *) if(ptcorrj>ptmax){ptmax=ptcorrj; selec=i;}} ///hardest jet selected - if(selec<0){PostData(1, fOutputList); - return;} + if(selec<0){PostData(1, fOutputList); + return;} AliAODJet* jet1 = (AliAODJet*)(fListJets[0]->At(selec)); //What is the hardest constituent track? AliAODTrack* leadtrack1; @@ -777,14 +780,13 @@ void AliAnalysisTaskAj::UserExec(Option_t *) Float_t ppjY = pPerp.Dot(ppJ3); TVector2 vr(ppjX, ppjY) ; //and this is the angle between the particle and the TM axis. - // Float_t phistr = evec.DeltaPhi(vr); + Float_t phistr = evec.DeltaPhi(vr); - Double_t phistr=vr.Phi()-evec.Phi(); + //Double_t phistr=vr.Phi()-evec.Phi(); if(centValue<10.) fh3LocalCoordinates->Fill(ppjX,ppjY,ptcorr2); Double_t deltaEta = eta2-track->Eta(); - if(phistr<-0.5*TMath::Pi()) phistr+=2.*TMath::Pi(); - if(phistr>3./2.*TMath::Pi()) phistr-=2.*TMath::Pi(); + if(deltaPhi<-0.5*TMath::Pi()) deltaPhi+=2.*TMath::Pi(); if(deltaPhi>3./2.*TMath::Pi()) deltaPhi-=2.*TMath::Pi(); @@ -1104,8 +1106,8 @@ void AliAnalysisTaskAj::GetDimParams(Int_t iEntry, TString &label, Int_t &nbins, case 7: label = "deltaPhiTM"; nbins = 60; - xmin = -0.5*pi; - xmax = 1.5*pi; + xmin = 0.; + xmax = 1.3*pi; break;