minor rebinning and jet axis semaring histos

diff --git a/PWGJE/UserTasks/AliAnalysisTaskJetAntenna.cxx b/PWGJE/UserTasks/AliAnalysisTaskJetAntenna.cxx
index 323251249e8d477aa291c717b48efbee8eebe218..457ebc6b9f649511cc0cfcf15301d8c185a9c789 100644 (file)
--- a/PWGJE/UserTasks/AliAnalysisTaskJetAntenna.cxx
+++ b/PWGJE/UserTasks/AliAnalysisTaskJetAntenna.cxx
@@ -226,11 +226,16 @@ void AliAnalysisTaskJetAntenna::UserCreateOutputObjects()
   fOutputList->Add(fh2JetEntries);
   fh2Circularity=new TH2F("Circcularity","",10,0,1,150,0,150);
   fOutputList->Add(fh2Circularity);
-  Int_t nbinsJet[9]={10,9,7,9,36,10,2,10,10};
+  fh2JetAxisPhi=new TH2F("JetAxisSmearPhi","",9,0,TMath::Pi(),10,-0.5,9.5);
+  fOutputList->Add(fh2JetAxisPhi);
+  
+
+  Int_t nbinsJet[9]={3,9,7,9,36,10,2,10,20};
   Double_t binlowJet[9]= {0,0, 0, 0,-0.5*TMath::Pi(),0,0,-0.5,0};
-  Double_t binupJet[9]= {100,0.9, 150,150,1.5*TMath::Pi(),1,200,9.5,20};
+  Double_t binupJet[9]= {100,0.9, 150,150,1.5*TMath::Pi(),1,200,9.5,200};
   fhnJetTM = new THnSparseF("fhnJetTM", "fhnJetTM; cent;dr;pt_jet;pt_track;phi;circ;nc;pthard",9,nbinsJet,binlowJet,binupJet);
-  Double_t *xPt3=new Double_t[10];
+ 
+ Double_t *xPt3=new Double_t[10];
   xPt3[0] = 0.;
   for(Int_t i = 1;i<=9;i++){
     if(xPt3[i-1]<2)xPt3[i] = xPt3[i-1] + 0.4; // 1 - 5
@@ -250,9 +255,22 @@ void AliAnalysisTaskJetAntenna::UserCreateOutputObjects()
   xPt2[7]=150; 
   
   fhnJetTM->SetBinEdges(2,xPt2);
-  fOutputList->Add(fhnJetTM);
+ 
+
+  Double_t *xPt4=new Double_t[4];
+  xPt4[0] = 0.;
+  xPt4[1]=10;
+  xPt4[2]=30;
+  xPt4[3]=50;
+
+  
+  fhnJetTM->SetBinEdges(0,xPt4);
+
+
+   fOutputList->Add(fhnJetTM);
    delete [] xPt3;
    delete [] xPt2;
+   delete [] xPt4;
   // =========== Switch on Sumw2 for all histos ===========
   for (Int_t i=0; i<fOutputList->GetEntries(); ++i) {
     TH1 *h1 = dynamic_cast<TH1*>(fOutputList->At(i));
@@ -475,11 +493,15 @@ void AliAnalysisTaskJetAntenna::UserExec(Option_t *)
     Double_t areabig=0;
     Double_t phibig=0.;
     Double_t pxbig,pybig,pzbig;
-
+    Double_t phitrue=0;
+    Double_t etatrue=0;
+    Double_t smearphi=0;
 
     AliAODJet* jetbig = (AliAODJet*)(fListJets[0]->At(i));
     etabig  = jetbig->Eta();
     phibig  = jetbig->Phi();
+    phitrue=phibig;
+    etatrue=etabig;
     ptbig   = jetbig->Pt();
     if(ptbig==0) continue;
     areabig = jetbig->EffectiveAreaCharged();
@@ -508,10 +530,22 @@ void AliAnalysisTaskJetAntenna::UserExec(Option_t *)
       
     }
     if(!jetAccepted) continue;
+    etabig=jetmatched->Eta();
+    phibig=jetmatched->Phi();
     pxbig=jetmatched->Px();
     pybig=jetmatched->Py();
     pzbig=jetmatched->Pz();
-    ptbig=jetmatched->Pt()-rho*jetmatched->EffectiveAreaCharged();}
+    ptbig=jetmatched->Pt()-rho*jetmatched->EffectiveAreaCharged();
+
+
+   
+      
+
+    smearphi=RelativePhi(phitrue,jetmatched->Phi());
+    smearphi=TMath::Abs(smearphi);
+   
+
+}
 
 
 
@@ -585,6 +619,10 @@ void AliAnalysisTaskJetAntenna::UserExec(Option_t *)
     fh2Circularity->Fill(circ,ptbig);
     fh2JetEntries->Fill(ptbig,pthardbin);
 
+
+    if(fDoMatching) fh2JetAxisPhi->Fill(smearphi,pthardbin);
+    
+
     
     for (Int_t ip = 0; ip < ParticleList.GetEntries(); ip++) {
       AliVParticle *track = (AliVParticle*)ParticleList.At(ip);

diff --git a/PWGJE/UserTasks/AliAnalysisTaskJetAntenna.h b/PWGJE/UserTasks/AliAnalysisTaskJetAntenna.h
index 726c8f341cda8766c0f3ee3d6979004e2b5ffa3f..594cf55a66be29d44876a59b6e10cc1df06694ea 100644 (file)
--- a/PWGJE/UserTasks/AliAnalysisTaskJetAntenna.h
+++ b/PWGJE/UserTasks/AliAnalysisTaskJetAntenna.h
@@ -158,7 +158,7 @@ private:
     TH1I  *fHistEvtSelection;                  //! event selection statistic
      TH2F*      fh2JetEntries;             //centrality bias of triggers 
      TH2F*      fh2Circularity;             //jet density
-    
+     TH2F*      fh2JetAxisPhi;              //smearing jet axis
      THnSparse   *fhnJetTM;               //Recoil jet spectrum