Added mother histo to Kstar config for TOF analysis + monitor histo binning changed

diff --git a/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskTOFKStar.C b/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskTOFKStar.C
index cf26a62e3e7b4d82c215fd2b44e3da0151347167..6abed6007beff0ccf7d94024d83658d0be428e06 100644 (file)
--- a/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskTOFKStar.C
+++ b/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskTOFKStar.C
@@ -103,8 +103,7 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskTOFKStar
    //event plane (only for PbPb)
    Int_t planeID = task->CreateValue(AliRsnMiniValue::kPlaneAngle, kFALSE);
    AliRsnMiniOutput *outPlane = task->CreateOutput("eventPlane", "HIST", "EVENT");
-   if (!isPP)
-     outPlane->AddAxis(planeID, 180, 0.0, TMath::Pi());
+   outPlane->AddAxis(planeID, 180, 0.0, TMath::Pi());
    
    //
    // -- PAIR CUTS (common to all resonances) ------------------------------------------------------

diff --git a/PWGLF/RESONANCES/macros/mini/AddMonitorOutput.C b/PWGLF/RESONANCES/macros/mini/AddMonitorOutput.C
index 93c5102bbf39e4561b4cac5b7b50b49b83784a31..68505721a905e3a0b87664cc16bbc20de3d3ad7c 100644 (file)
--- a/PWGLF/RESONANCES/macros/mini/AddMonitorOutput.C
+++ b/PWGLF/RESONANCES/macros/mini/AddMonitorOutput.C
@@ -11,7 +11,7 @@ void AddMonitorOutput(Bool_t useMCMon = 0, TObjArray *mon=0,TString opt="",AliRs
   
   // Multiplicity/centrality
   AliRsnValueEvent *multi = new AliRsnValueEvent("multi",AliRsnValueEvent::kMult);
-  multi->SetBins(0.0, 100.0, 1);
+  multi->SetBins(0.0, 400.0, 1);
   // Momentum
   AliRsnValueDaughter *axisMomTPC = new AliRsnValueDaughter("pTPC", AliRsnValueDaughter::kPtpc);
   axisMomTPC->SetBins(0.0, 10.0, 0.01);
@@ -22,33 +22,33 @@ void AddMonitorOutput(Bool_t useMCMon = 0, TObjArray *mon=0,TString opt="",AliRs
   axisMomPt->SetBins(0.0,10.0,0.01);
   // Eta
   AliRsnValueDaughter *axisMomEta = new AliRsnValueDaughter("eta", AliRsnValueDaughter::kEta);
-  axisMomEta->SetBins(-1.0, 1.0, 0.05);
+  axisMomEta->SetBins(-1.0, 1.0, 0.1);
 
   /* dEdx tpc */
   AliRsnValueDaughter *axisSigTPC = new AliRsnValueDaughter("sTPC", AliRsnValueDaughter::kTPCsignal);
   axisSigTPC->SetBins(0.0, 500.0, 2.0);
   // kTPCnsigmaPi
   AliRsnValueDaughter *axisTPCnsigmaPi = new AliRsnValueDaughter("pi", AliRsnValueDaughter::kTPCnsigmaPi);
-  axisTPCnsigmaPi->SetBins(-50.,50., 0.1);
+  axisTPCnsigmaPi->SetBins(-10.,10., 0.1);
   // kTPCnsigmaK
   AliRsnValueDaughter *axisTPCnsigmaK = new AliRsnValueDaughter("K", AliRsnValueDaughter::kTPCnsigmaK);
-  axisTPCnsigmaK->SetBins(-50.,50., 0.1);
+  axisTPCnsigmaK->SetBins(-10.,10., 0.1);
   // kTPCnsigmaP
   AliRsnValueDaughter *axisTPCnsigmaP = new AliRsnValueDaughter("p", AliRsnValueDaughter::kTPCnsigmaP);
-  axisTPCnsigmaP->SetBins(-50.,50., 0.1);
+  axisTPCnsigmaP->SetBins(-10.,10., 0.1);
 
   /* tof signal - time */
   AliRsnValueDaughter *axisSigTOF = new AliRsnValueDaughter("sTOF", AliRsnValueDaughter::kTOFsignal);
   axisSigTOF->SetBins(10000.0, 50000.0, 250);//in ps
   // kTOFnsigmaPi
   AliRsnValueDaughter *axisTOFnsigmaPi = new AliRsnValueDaughter("pi", AliRsnValueDaughter::kTOFnsigmaPi);
-  axisTOFnsigmaPi->SetBins(-50.,50., 0.1);
+  axisTOFnsigmaPi->SetBins(-10.,10., 0.1);
   // kTOFnsigmaK
   AliRsnValueDaughter *axisTOFnsigmaK = new AliRsnValueDaughter("K", AliRsnValueDaughter::kTOFnsigmaK);
-  axisTOFnsigmaK->SetBins(-50.,50., 0.1);
+  axisTOFnsigmaK->SetBins(-10.,10., 0.1);
   // kTOFnsigmaP
   AliRsnValueDaughter *axisTOFnsigmaP = new AliRsnValueDaughter("p", AliRsnValueDaughter::kTOFnsigmaP);
-  axisTOFnsigmaP->SetBins(-50.,50., 0.1);
+  axisTOFnsigmaP->SetBins(-10.,10., 0.1);
    
   /****************************************************************/
   /***************       MONITOR AOB           ********************/
@@ -167,7 +167,7 @@ void AddMonitorOutput(Bool_t useMCMon = 0, TObjArray *mon=0,TString opt="",AliRs
     // Eta
     AliRsnValueDaughter *axisEtaMC = new AliRsnValueDaughter("etaMC", AliRsnValueDaughter::kEta);
     axisEtaMC->SetUseMCInfo(kTRUE);
-    axisEtaMC->SetBins(-1.0,1.0,0.01);
+    axisEtaMC->SetBins(-1.0,1.0,0.1);
 
     // output: 2D histo for kine acceptance
     AliRsnListOutput *outMonitorPtVsEtaMC = new AliRsnListOutput("Pt_VsEtaMC", AliRsnListOutput::kHistoDefault);

diff --git a/PWGLF/RESONANCES/macros/mini/ConfigTOFanalysisKStar.C b/PWGLF/RESONANCES/macros/mini/ConfigTOFanalysisKStar.C
index 81f15b2e05a7eb2542f97808aae94aed532a1949..5bde65af77b027d5cd50df198e8586542bdb343b 100644 (file)
--- a/PWGLF/RESONANCES/macros/mini/ConfigTOFanalysisKStar.C
+++ b/PWGLF/RESONANCES/macros/mini/ConfigTOFanalysisKStar.C
@@ -111,5 +111,25 @@ Bool_t ConfigTOFanalysisKStar
     // out->AddAxis(yID, 10, -0.5, 0.5);
     
   }   
+  
+  if (isMC){   
+   // create output
+    AliRsnMiniOutput *outm = task->CreateOutput(Form("kstar_Mother%s", suffix), "SPARSE", "MOTHER");
+   // selection settings
+   outm->SetDaughter(0, AliRsnDaughter::kKaon);
+   outm->SetDaughter(1, AliRsnDaughter::kPion);
+   outm->SetMotherPDG(313);
+   outm->SetMotherMass(0.89594);
+   // pair cuts
+   outm->SetPairCuts(cutsPair);
+   // binnings
+   outm->AddAxis(imID, 90, 0.6, 1.5);
+   outm->AddAxis(ptID, 100, 0.0, 10.0);
+   
+   if (!isPP)
+     outm->AddAxis(centID, 100, 0.0, 100.0);
+   else 
+     outm->AddAxis(centID, 400, 0.0, 400.0);
+  }
   return kTRUE;
 }