Adding histograms for smearing normized to p_T,Gen
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 26 Jul 2010 15:48:02 +0000 (15:48 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 26 Jul 2010 15:48:02 +0000 (15:48 +0000)
PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.cxx
PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.h

index fccdfbb..8770fa8 100644 (file)
@@ -112,6 +112,7 @@ AliAnalysisTaskJetSpectrum2::AliAnalysisTaskJetSpectrum2(): AliAnalysisTaskSE(),
   fh2TracksLeadingJetPhiPt(0x0),
   fh2JetPtJetPhi(0x0),
   fh2TrackPtTrackPhi(0x0),
+  fh2RelPtFGen(0x0),
   fh2DijetDeltaPhiPt(0x0),      
   fh2DijetAsymPt(0x0),          
   fh2DijetAsymPtCut(0x0),       
@@ -193,6 +194,7 @@ AliAnalysisTaskJetSpectrum2::AliAnalysisTaskJetSpectrum2(const char* name):
   fh2TracksLeadingJetPhiPt(0x0),
   fh2JetPtJetPhi(0x0),
   fh2TrackPtTrackPhi(0x0),
+  fh2RelPtFGen(0x0),
   fh2DijetDeltaPhiPt(0x0),      
   fh2DijetAsymPt(0x0),          
   fh2DijetAsymPtCut(0x0),       
@@ -371,7 +373,8 @@ void AliAnalysisTaskJetSpectrum2::UserCreateOutputObjects()
   fh2JetPtJetPhi = new TH2F("fh2JetPtJetPhi","Reconstructed jet phi vs. pt",nBinPt,binLimitsPt,nBinPhi2,binLimitsPhi2);
   fh2TrackPtTrackPhi = new TH2F("fh2TrackPtTrackPhi","Reconstructed track phi vs. pt",nBinPt,binLimitsPt,nBinPhi2,binLimitsPhi2);
 
-
+  fh2RelPtFGen = new TH2F("fh2RelPtFGen",";p_{T,gen};#Delta p_{T}/p_{T,Gen}",nBinPt,binLimitsPt,120,-1.2,1.2);
+  
   for(int ij = 0;ij < kMaxJets;++ij){
     fh1PtRecIn[ij] = new TH1F(Form("fh1PtRecIn_j%d",ij),"rec p_T input ;p_{T,rec}",nBinPt,binLimitsPt);
     fh1PtGenIn[ij] = new TH1F(Form("fh1PtGenIn_j%d",ij),"found p_T input ;p_{T,gen}",nBinPt,binLimitsPt);
@@ -466,6 +469,7 @@ void AliAnalysisTaskJetSpectrum2::UserCreateOutputObjects()
     fHistList->Add(fhnCorrelationPhiZRec);
     fHistList->Add(fh2JetPtJetPhi);
     fHistList->Add(fh2TrackPtTrackPhi);
+    fHistList->Add(fh2RelPtFGen);
 
     fHistList->Add(fh2DijetDeltaPhiPt);       
     fHistList->Add(fh2DijetAsymPt);       
@@ -1068,6 +1072,10 @@ void AliAnalysisTaskJetSpectrum2::UserExec(Option_t */*option*/)
       if(TMath::Abs(etaRec)<fRecEtaWindow){
        fhnJetContainer[kStep3+kMaxStep]->Fill(container,eventW);
        fhnCorrelation->Fill(container);
+       if(ptGen>0){
+         Float_t delta = (ptRec-ptGen)/ptGen;
+         fh2RelPtFGen->Fill(ptGen,delta,eventW);
+       }
        if(fhnCorrelationPhiZRec)fhnCorrelationPhiZRec->Fill(containerPhiZ);
 
       }// if etarec in window
index dac0159..c9f86eb 100644 (file)
@@ -153,6 +153,7 @@ class AliAnalysisTaskJetSpectrum2 : public AliAnalysisTaskSE
     TH2F*         fh2TracksLeadingJetPhiPt; // track correlation with leading track
     TH2F*         fh2JetPtJetPhi;           // Phi distribution of accepted jets 
     TH2F*         fh2TrackPtTrackPhi;       // phi distribution of accepted tracks
+    TH2F*         fh2RelPtFGen;             // relative difference between generated and found 
     TH2F*         fh2PhiPt[kMaxJets];    // delta phi correlation of tracks with the jet      
     TH2F*         fh2PhiEta[kMaxJets];   // eta   phi correlation of tracks with the jet      
     TH2F*         fh2RhoPtRec[kMaxJets];    // jet shape variable rho
@@ -164,6 +165,8 @@ class AliAnalysisTaskJetSpectrum2 : public AliAnalysisTaskSE
     TH2F*         fh2FragGen[kMaxJets];     // fragmentation function
     TH2F*         fh2FragLnGen[kMaxJets];   // fragmetation in xi
 
+
+
     // Dijet histos
     TH2F*   fh2DijetDeltaPhiPt;      // dijet delta phi vs pt
     TH2F*   fh2DijetAsymPt;          // dijet asym vs pt
@@ -178,7 +181,7 @@ class AliAnalysisTaskJetSpectrum2 : public AliAnalysisTaskSE
     TList *fHistList; // Output list
    
 
-    ClassDef(AliAnalysisTaskJetSpectrum2, 5) // Analysis task for standard jet analysis
+    ClassDef(AliAnalysisTaskJetSpectrum2, 6) // Analysis task for standard jet analysis
 };
  
 #endif