Set histogram names based on macro input
authorSvein Lindal <svein.lindal@gmail.com>
Thu, 23 Jan 2014 15:17:37 +0000 (10:17 -0500)
committerSvein Lindal <svein.lindal@gmail.com>
Thu, 23 Jan 2014 15:18:35 +0000 (10:18 -0500)
only get histograms from filters if they are created

PWGGA/GammaConv/AliAnalysisTaskdPhi.cxx
PWGGA/GammaConv/AliConversionTrackCuts.cxx

index 17aa4a3..b88cad3 100644 (file)
@@ -243,32 +243,32 @@ void AliAnalysisTaskdPhi::UserCreateOutputObjects() {
   for(Int_t igf = 0; igf < fV0Filters[0].GetEntriesFast(); igf ++){
     AliConversionCuts * f = dynamic_cast<AliConversionCuts*>(fV0Filters[0].At(igf));
     if(f) {
-      f->InitCutHistograms(Form("V0Filter_%d", -(igf+1)), kFALSE);
-      fHistograms->Add(f->GetCutHistograms());
+      TList * histograms = f->GetCutHistograms();
+      if(histograms) fHistograms->Add(f->GetCutHistograms());
     }
   }
 
   for(Int_t igf = 0; igf < fV0Filters[1].GetEntriesFast(); igf ++){
     AliConversionCuts * f = dynamic_cast<AliConversionCuts*>(fV0Filters[1].At(igf));
     if(f) {
-      f->InitCutHistograms(Form("V0Filter_%d", igf+1), kFALSE);
-      fHistograms->Add(f->GetCutHistograms());
+      TList * histograms = f->GetCutHistograms();
+      if(histograms) fHistograms->Add(f->GetCutHistograms());
     }
   }
 
   for(Int_t igf = 0; igf < fMesonFilters[0].GetEntriesFast(); igf ++){
     AliConversionMesonCuts * f = dynamic_cast<AliConversionMesonCuts*>(fMesonFilters[0].At(igf));
     if(f) {
-      f->InitCutHistograms(Form("PionFilter_%d", -(igf+1)), kFALSE);
-      fHistograms->Add(f->GetCutHistograms());
+      TList * histograms = f->GetCutHistograms();
+      if(histograms) fHistograms->Add(f->GetCutHistograms());
     }
   }
 
   for(Int_t igf = 0; igf < fMesonFilters[1].GetEntriesFast(); igf ++){
     AliConversionMesonCuts * f = dynamic_cast<AliConversionMesonCuts*>(fMesonFilters[1].At(igf));
     if(f) {
-      f->InitCutHistograms(Form("PionFilter_%d", igf+1), kFALSE);
-      fHistograms->Add(f->GetCutHistograms());
+      TList * histograms = f->GetCutHistograms();
+      if(histograms) fHistograms->Add(f->GetCutHistograms());
     }
   }
 
@@ -558,15 +558,12 @@ void AliAnalysisTaskdPhi::UserExec(Option_t *) {
 
 
   Double_t centrality = 0.0;
-  Double_t eventPlane = 0.0;
   Double_t vertexz = fInputEvent->GetPrimaryVertex()->GetZ();
   if(isAOD) {
     AliAODHeader * header = static_cast<AliAODHeader*>(fInputEvent->GetHeader());
     centrality = header->GetCentrality();
-    eventPlane = header->GetEventplane();
   } else {
     centrality = static_cast<AliESDEvent*>(fInputEvent)->GetCentrality()->GetCentralityPercentile("V0M");
-    eventPlane = fInputEvent->GetEventplane()->GetEventplane("Q");
   }
   
   
@@ -577,7 +574,6 @@ void AliAnalysisTaskdPhi::UserExec(Option_t *) {
   if(DebugLevel () > 4) {
     cout << "centrality: " << centrality <<  " " << GetBin(fAxisCent, centrality) << endl;
     cout << "vertexz: " << vertexz <<  " " << GetBin(fAxisZ, vertexz) << endl;
-    cout << "eventPlane: " << eventPlane <<  " " << endl;
   }
   
   
index 47a6ef5..4568e58 100644 (file)
@@ -359,7 +359,7 @@ TList * AliConversionTrackCuts::CreateHistograms() {
   fHistograms->SetOwner(kTRUE);
   fHistograms->SetName("trackCuts");
 
-  fhPhi = new TH2F("phi", "phi", 5, -0.5, 4.5, 32, 0, TMath::TwoPi());
+  fhPhi = new TH2F(Form("phi_%s", GetName()), Form("phi_%s", GetTitle()), 5, -0.5, 4.5, 32, 0, TMath::TwoPi());
   // TAxis * xax = fhPhi->GetXaxis();
   // for(Int_t i = 0; i < kNCuts; i++){
   //   xax->SetBinLabel(xax->FindFixBin(i), fgkCutNames[i]);
@@ -367,7 +367,7 @@ TList * AliConversionTrackCuts::CreateHistograms() {
   fHistograms->Add(fhPhi);
   
 
-  fhEtaPhi = new TH2F("etahpi", "etaphi", 36, -0.9, 0.9, 32, 0, TMath::TwoPi());
+  fhEtaPhi = new TH2F(Form("etaphi_%s",GetName()), Form("etaphi_%s", GetTitle()), 36, -0.9, 0.9, 32, 0, TMath::TwoPi());
   fHistograms->Add(fhEtaPhi);
 
   // fhPt = new TH2F("pt", "pt", kNCuts+2, kPreCut -0.5, kNCuts + 0.5, 
@@ -381,13 +381,13 @@ TList * AliConversionTrackCuts::CreateHistograms() {
   //  fhPhiPt = new TH2F("phipt", "phipt", 100, 0, 100, 64, 0, TMath::TwoPi());
   //fHistograms->Add(fhPhiPt);
 
-  fhdcaxyPt = new TH2F("dcaxypt", "dcaxypt", 20, 0, 20, 50, -2.5, 2.5);
+  fhdcaxyPt = new TH2F(Form("dcaxypt_%s", GetName()),  Form("dcaxypt_%s", GetTitle()), 20, 0, 20, 50, -2.5, 2.5);
   fHistograms->Add(fhdcaxyPt);
 
-  fhdcazPt = new TH2F("dcazpt", "dcazpt", 20, 0, 20, 70, -3.5, 3.5);
+  fhdcazPt = new TH2F(Form("dcazpt_%s", GetName()),  Form("dcazpt_%s", GetTitle()), 20, 0, 20, 70, -3.5, 3.5);
   fHistograms->Add(fhdcazPt);
 
-  fhdca = new TH2F("dca", "dca", 70, -3.5, 3.5, 50, -2.5, 2.5);
+  fhdca = new TH2F(Form("dca_%s", GetName()),  Form("dca_%s", GetTitle()), 70, -3.5, 3.5, 50, -2.5, 2.5);
   fhdca->SetXTitle("dca z");
   fhdca->SetYTitle("dca xy");