AliAnaPi0: Change name of histograms in memory produced during projections, since...
authorgconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 30 Dec 2009 18:10:34 +0000 (18:10 +0000)
committergconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 30 Dec 2009 18:10:34 +0000 (18:10 +0000)
AliAnaCalorimeterQA: Correct name of histogram fhTimeAmp
                                     Fill Invariant mass histograms with pt and not energy in order to compare results with AliAnaPi0

PWG4/PartCorrDep/AliAnaCalorimeterQA.cxx
PWG4/PartCorrDep/AliAnaPi0.cxx

index 658793a..47fc466 100755 (executable)
@@ -427,17 +427,17 @@ TList *  AliAnaCalorimeterQA::GetCreateOutputObjects()
        outputContainer->Add(fh1pOverER02);     
        
        fhIM  = new TH2F ("hIM","Cluster pairs Invariant mass vs reconstructed pair energy",nptbins,ptmin,ptmax,nmassbins,massmin,massmax); 
-       fhIM->SetXTitle("E_{cluster pairs} (GeV) ");
+       fhIM->SetXTitle("p_{T, cluster pairs} (GeV) ");
        fhIM->SetYTitle("M_{cluster pairs} (GeV/c^{2})");
        outputContainer->Add(fhIM);
        
        fhIMCellCut  = new TH2F ("hIMCellCut","Cluster (n cell > 1) pairs Invariant mass vs reconstructed pair energy",nptbins,ptmin,ptmax,nmassbins,massmin,massmax); 
-       fhIMCellCut->SetXTitle("E_{cluster pairs} (GeV) ");
+       fhIMCellCut->SetXTitle("p_{T, cluster pairs} (GeV) ");
        fhIMCellCut->SetYTitle("M_{cluster pairs} (GeV/c^{2})");
        outputContainer->Add(fhIMCellCut);
        
        fhAsym  = new TH2F ("hAssym","Cluster pairs Asymmetry vs reconstructed pair energy",nptbins,ptmin,ptmax,nasymbins,asymmin,asymmax); 
-       fhAsym->SetXTitle("E_{cluster pairs} (GeV) ");
+       fhAsym->SetXTitle("p_{T, cluster pairs} (GeV) ");
        fhAsym->SetYTitle("Asymmetry");
        outputContainer->Add(fhAsym);   
        
@@ -470,7 +470,7 @@ TList *  AliAnaCalorimeterQA::GetCreateOutputObjects()
                fhTimeId->SetYTitle("Cell Absolute Id");
                outputContainer->Add(fhTimeId);
        
-               fhTimeAmp  = new TH2F ("hTimeId","Cell Time vs Absolute Id",nptbins*2,ptmin,ptmax,ntimebins,timemin,timemax); 
+               fhTimeAmp  = new TH2F ("hTimeAmp","Cell Time vs Absolute Id",nptbins*2,ptmin,ptmax,ntimebins,timemin,timemax); 
                fhTimeAmp->SetYTitle("Cell Time (ns)");
                fhTimeAmp->SetXTitle("Cell Energy (GeV)");
                outputContainer->Add(fhTimeAmp);
@@ -548,14 +548,14 @@ TList *  AliAnaCalorimeterQA::GetCreateOutputObjects()
                fhIMMod[imod]  = new TH2F (Form("hIM_Mod%d",imod),
                                                                   Form("Cluster pairs Invariant mass vs reconstructed pair energy in Module %d",imod),
                                                                   nptbins,ptmin,ptmax,nmassbins,massmin,massmax); 
-               fhIMMod[imod]->SetXTitle("E_{cluster pairs} (GeV) ");
+               fhIMMod[imod]->SetXTitle("p_{T, cluster pairs} (GeV) ");
                fhIMMod[imod]->SetYTitle("M_{cluster pairs} (GeV/c^{2})");
                outputContainer->Add(fhIMMod[imod]);
                                
                fhIMCellCutMod[imod]  = new TH2F (Form("hIMCellCut_Mod%d",imod),
                                                                   Form("Cluster (n cells > 1) pairs Invariant mass vs reconstructed pair energy in Module %d",imod),
                                                                   nptbins,ptmin,ptmax,nmassbins,massmin,massmax); 
-               fhIMCellCutMod[imod]->SetXTitle("E_{cluster pairs} (GeV) ");
+               fhIMCellCutMod[imod]->SetXTitle("p_{T, cluster pairs} (GeV) ");
                fhIMCellCutMod[imod]->SetYTitle("M_{cluster pairs} (GeV/c^{2})");
                outputContainer->Add(fhIMCellCutMod[imod]);
                
@@ -1364,22 +1364,22 @@ void  AliAnaCalorimeterQA::MakeAnalysisFillHistograms()
 
                                //Fill invariant mass histograms
                                //All modules
-                               fhIM  ->Fill((mom+mom2).E(),(mom+mom2).M());
+                               fhIM  ->Fill((mom+mom2).Pt(),(mom+mom2).M());
                                //Single module
                                if(nModule == nModule2 && nModule >=0 && nModule < fNModules)
-                                 fhIMMod[nModule]->Fill((mom+mom2).E(),(mom+mom2).M());
+                                 fhIMMod[nModule]->Fill((mom+mom2).Pt(),(mom+mom2).M());
 
                                //Select only clusters with at least 2 cells
                                if(nCaloCellsPerCluster > 1 && nCaloCellsPerCluster2 > 1) {
                                  //All modules
-                                 fhIMCellCut  ->Fill((mom+mom2).E(),(mom+mom2).M());
+                                 fhIMCellCut  ->Fill((mom+mom2).Pt(),(mom+mom2).M());
                                  //Single modules
                                  if(nModule == nModule2 && nModule >=0 && nModule < fNModules)
-                                   fhIMCellCutMod[nModule]->Fill((mom+mom2).E(),(mom+mom2).M());
+                                   fhIMCellCutMod[nModule]->Fill((mom+mom2).Pt(),(mom+mom2).M());
                                }
 
                                //Asymetry histograms
-                               fhAsym->Fill((mom+mom2).E(),TMath::Abs((mom.E()-mom2.E())/(mom.E()+mom2.E())));
+                               fhAsym->Fill((mom+mom2).Pt(),TMath::Abs((mom.E()-mom2.E())/(mom.E()+mom2.E())));
                                        
                        }// 2nd cluster loop
                }////more than 1 cluster in calorimeter         
@@ -2400,7 +2400,7 @@ void  AliAnaCalorimeterQA::Terminate(TList* outputList)
        
                ctime->cd(3) ; 
                //if(fhTimeAmp->GetEntries() > 0) gPad->SetLogy();
-               fhTimeAmp->Draw("colz");
+               fhTimeAmp->Draw("surf");
        
                sprintf(name,"QA_%s_CellsTime.eps",fCalorimeter.Data());
                ctime->Print(name); printf("Plot: %s\n",name);
index 4e02251..38ead0c 100755 (executable)
@@ -679,31 +679,31 @@ void AliAnaPi0::Terminate(TList* outputList)
 
   cIM->cd(1) ; 
   //gPad->SetLogy();
-  TH1D * hIMAllPt = (TH1D*) fhRe1[0]->ProjectionZ();
+  TH1D * hIMAllPt = (TH1D*) fhRe1[0]->ProjectionZ(Form("IMPtAll_%s",fCalorimeter.Data()));
   hIMAllPt->SetLineColor(2);
   hIMAllPt->SetTitle("No cut on  p_{T, #gamma#gamma} ");
   hIMAllPt->Draw();
 
   cIM->cd(2) ; 
-  TH3F * hRe1Pt5 = (TH3F*)fhRe1[0]->Clone("IMPt5");
+  TH3F * hRe1Pt5 = (TH3F*)fhRe1[0]->Clone(Form("IMPt5_%s",fCalorimeter.Data()));
   hRe1Pt5->GetXaxis()->SetRangeUser(0,5);
-  TH1D * hIMPt5 = (TH1D*) hRe1Pt5->Project3D("z");
+  TH1D * hIMPt5 = (TH1D*) hRe1Pt5->Project3D(Form("IMPt5_%s_pz",fCalorimeter.Data()));
   hIMPt5->SetLineColor(2);  
   hIMPt5->SetTitle("0 < p_{T, #gamma#gamma} < 5 GeV/c");
   hIMPt5->Draw();
   
   cIM->cd(3) ; 
-  TH3F * hRe1Pt10 =  (TH3F*)fhRe1[0]->Clone("IMPt10");
+  TH3F * hRe1Pt10 =  (TH3F*)fhRe1[0]->Clone(Form("IMPt10_%s",fCalorimeter.Data()));
   hRe1Pt10->GetXaxis()->SetRangeUser(5,10);
-  TH1D * hIMPt10 = (TH1D*) hRe1Pt10->Project3D("z");
+  TH1D * hIMPt10 = (TH1D*) hRe1Pt10->Project3D(Form("IMPt10_%s_pz",fCalorimeter.Data()));
   hIMPt10->SetLineColor(2);  
   hIMPt10->SetTitle("5 < p_{T, #gamma#gamma} < 10 GeV/c");
   hIMPt10->Draw();
   
   cIM->cd(4) ; 
-  TH3F * hRe1Pt20 =  (TH3F*)fhRe1[0]->Clone("IMPt20");
+  TH3F * hRe1Pt20 =  (TH3F*)fhRe1[0]->Clone(Form("IMPt20_%s",fCalorimeter.Data()));
   hRe1Pt20->GetXaxis()->SetRangeUser(10,20);
-  TH1D * hIMPt20 = (TH1D*) hRe1Pt20->Project3D("z");
+  TH1D * hIMPt20 = (TH1D*) hRe1Pt20->Project3D(Form("IMPt20_%s_pz",fCalorimeter.Data()));
   hIMPt20->SetLineColor(2);  
   hIMPt20->SetTitle("10 < p_{T, #gamma#gamma} < 20 GeV/c");
   hIMPt20->Draw();
@@ -725,7 +725,7 @@ void AliAnaPi0::Terminate(TList* outputList)
   hPt->Draw();
 
   cPt->cd(2) ; 
-  TH3F * hRe1IM1 = (TH3F*)fhRe1[0]->Clone("PtIM1");
+  TH3F * hRe1IM1 = (TH3F*)fhRe1[0]->Clone(Form("Pt1_%s",fCalorimeter.Data()));
   hRe1IM1->GetZaxis()->SetRangeUser(0.05,0.21);
   TH1D * hPtIM1 = (TH1D*) hRe1IM1->Project3D("x");
   hPtIM1->SetLineColor(2);  
@@ -733,7 +733,7 @@ void AliAnaPi0::Terminate(TList* outputList)
   hPtIM1->Draw();
   
   cPt->cd(3) ; 
-  TH3F * hRe1IM2 = (TH3F*)fhRe1[0]->Clone("PtIM2");
+  TH3F * hRe1IM2 = (TH3F*)fhRe1[0]->Clone(Form("Pt2_%s",fCalorimeter.Data()));
   hRe1IM2->GetZaxis()->SetRangeUser(0.09,0.17);
   TH1D * hPtIM2 = (TH1D*) hRe1IM2->Project3D("x");
   hPtIM2->SetLineColor(2);  
@@ -741,7 +741,7 @@ void AliAnaPi0::Terminate(TList* outputList)
   hPtIM2->Draw();
 
   cPt->cd(4) ; 
-  TH3F * hRe1IM3 = (TH3F*)fhRe1[0]->Clone("PtIM3");
+  TH3F * hRe1IM3 = (TH3F*)fhRe1[0]->Clone(Form("Pt3_%s",fCalorimeter.Data()));
   hRe1IM3->GetZaxis()->SetRangeUser(0.11,0.15);
   TH1D * hPtIM3 = (TH1D*) hRe1IM1->Project3D("x");
   hPtIM3->SetLineColor(2);