- Added method for centrality flattening in 2011 PbPb data analysis (Subhash)
authorfbellini <fbellini@cern.ch>
Thu, 17 Jul 2014 12:58:45 +0000 (14:58 +0200)
committerfbellini <fbellini@cern.ch>
Thu, 17 Jul 2014 13:01:09 +0000 (15:01 +0200)
- Modified axis range for K* output (Subhash)
- Improved macros for D0 analysis (Massimo)

PWGLF/RESONANCES/AliRsnMiniAnalysisTask.cxx
PWGLF/RESONANCES/AliRsnMiniAnalysisTask.h
PWGLF/RESONANCES/macros/mini/AddAnalysisTaskD0.C
PWGLF/RESONANCES/macros/mini/AddAnalysisTaskTPCKStar.C
PWGLF/RESONANCES/macros/mini/ConfigD0.C
PWGLF/RESONANCES/macros/mini/ConfigTPCanalysisKStar.C
PWGLF/RESONANCES/macros/mini/ConfigTPCanalysisKStarMC.C

index 471deb3..cc3137f 100644 (file)
@@ -16,6 +16,7 @@
 #include <TList.h>
 #include <TTree.h>
 #include <TStopwatch.h>
+#include "TRandom.h"
 
 #include "AliLog.h"
 #include "AliEventplane.h"
@@ -49,6 +50,7 @@ AliRsnMiniAnalysisTask::AliRsnMiniAnalysisTask() :
    fUseCentrality(kFALSE),
    fCentralityType("QUALITY"),
    fUseAOD049CentralityPatch(kFALSE),
+   fUseCentralityPatchPbPb2011(0),
    fContinuousMix(kTRUE),
    fNMix(0),
    fMaxDiffMult(10),
@@ -93,6 +95,7 @@ AliRsnMiniAnalysisTask::AliRsnMiniAnalysisTask(const char *name, Bool_t useMC) :
    fUseCentrality(kFALSE),
    fCentralityType("QUALITY"),
    fUseAOD049CentralityPatch(kFALSE),
+  fUseCentralityPatchPbPb2011(0),
    fContinuousMix(kTRUE),
    fNMix(0),
    fMaxDiffMult(10),
@@ -142,6 +145,7 @@ AliRsnMiniAnalysisTask::AliRsnMiniAnalysisTask(const AliRsnMiniAnalysisTask &cop
    fUseCentrality(copy.fUseCentrality),
    fCentralityType(copy.fCentralityType),
    fUseAOD049CentralityPatch(copy.fUseAOD049CentralityPatch),
+   fUseCentralityPatchPbPb2011(copy.fUseCentralityPatchPbPb2011),
    fContinuousMix(copy.fContinuousMix),
    fNMix(copy.fNMix),
    fMaxDiffMult(copy.fMaxDiffMult),
@@ -197,6 +201,7 @@ AliRsnMiniAnalysisTask &AliRsnMiniAnalysisTask::operator=(const AliRsnMiniAnalys
    fUseCentrality = copy.fUseCentrality;
    fCentralityType = copy.fCentralityType;
    fUseAOD049CentralityPatch = copy.fUseAOD049CentralityPatch;
+   fUseCentralityPatchPbPb2011 = copy.fUseCentralityPatchPbPb2011;
    fContinuousMix = copy.fContinuousMix;
    fNMix = copy.fNMix;
    fMaxDiffMult = copy.fMaxDiffMult;
@@ -822,17 +827,19 @@ Double_t AliRsnMiniAnalysisTask::ComputeCentrality(Bool_t isESD)
 //
 
    if (fUseCentrality) {
-
-      if ((!fUseMC) && (!isESD) && (fUseAOD049CentralityPatch)) {
-         return ApplyCentralityPatchAOD049();
-      } else {
-         AliCentrality *centrality = fInputEvent->GetCentrality();
+     if ((!fUseMC) && (fUseCentralityPatchPbPb2011)) {
+       return ApplyCentralityPatchPbPb2011();//
+    }
+     if ((!fUseMC) && (!isESD) && (fUseAOD049CentralityPatch)) {
+       return ApplyCentralityPatchAOD049();
+     } else {
+       AliCentrality *centrality = fInputEvent->GetCentrality();
          if (!centrality) {
-            AliError("Cannot compute centrality!");
-            return -1.0;
+          AliError("Cannot compute centrality!");
+          return -1.0;
          }
          return centrality->GetCentralityPercentile(fCentralityType.Data());
-      }
+     }
    } else {
       if (!fCentralityType.CompareTo("TRACKS"))
          return fInputEvent->GetNumberOfTracks();
@@ -1187,7 +1194,54 @@ Bool_t AliRsnMiniAnalysisTask::EventsMatch(AliRsnMiniEvent *event1, AliRsnMiniEv
    }
 }
 
-
+//---------------------------------------------------------------------
+Double_t AliRsnMiniAnalysisTask::ApplyCentralityPatchPbPb2011(){
+  //This part rejects randomly events such that the centrality gets flat for LHC11h Pb-Pb data
+  //for 0-5% and 10-20% centrality bin
+  
+  if (fCentralityType!="V0M") {
+    AliWarning("Wrong value (not centrality from V0).");
+    return -999.0;
+  }
+  
+  AliCentrality *centrality = fInputEvent->GetCentrality();
+  if (!centrality) {
+    AliWarning("Cannot get centrality from AOD event.");
+    return -999.0;
+  }
+  
+  
+  Double_t cent = (Float_t)(centrality->GetCentralityPercentile("V0M"));               
+  Double_t rnd_hc, testf, ff, N1, N2;
+
+  if(fUseCentralityPatchPbPb2011==5){
+    N1 = 1.9404e+06;
+    N2 = 1.56435e+06;
+    ff = 5.04167e+06 - 1.49885e+06*cent + 2.35998e+05*cent*cent -1.22873e+04*cent*cent*cent;
+  }
+  
+  if(fUseCentralityPatchPbPb2011==1020){
+    N1 = 1.56435e+06;
+    N2 = 4.20e+05;
+    ff = 1.68062e+08 - 5.19673e+07*cent + 6.4068e+06*cent*cent + 6.4068e+06*cent*cent*cent - 392687*cent*cent*cent*cent - 145.07*cent*cent*cent*cent*cent; 
+  }
+
+  testf = ( N2 + (N1-ff) ) / N1;
+  rnd_hc = gRandom->Rndm();
+
+  //AliDebugClass(1, Form("Flat Centrality %d", fUseCentralityPatchPbPb2011));
+
+  if (rnd_hc < 0 || rnd_hc > 1 ) 
+    {
+      AliWarning("Wrong Random number generated");
+      return -999.0;
+    }
+  
+  if (rnd_hc < testf)
+    return cent;
+  else
+    return -999.0;
+}
 //---------------------------------------------------------------------
 Double_t AliRsnMiniAnalysisTask::ApplyCentralityPatchAOD049()
 {
index 69970e0..50afb09 100644 (file)
@@ -39,6 +39,7 @@ public:
    void                UseESDTriggerMask(UInt_t trgMask)     {fTriggerMask = trgMask;}
    void                UseCentrality(const char *type)    {fUseCentrality = kTRUE; fCentralityType = type; fCentralityType.ToUpper();}
    void                SetUseCentralityPatch(Bool_t isAOD049) {fUseAOD049CentralityPatch = isAOD049;}
+   void                SetUseCentralityPatchPbPb2011(Int_t centralityPatchPbPb2011) {fUseCentralityPatchPbPb2011 = centralityPatchPbPb2011;}
    void                UseMultiplicity(const char *type)  {fUseCentrality = kFALSE; fCentralityType = type; fCentralityType.ToUpper();}
    void                UseContinuousMix()                 {fContinuousMix = kTRUE;}
    void                UseBinnedMix()                     {fContinuousMix = kFALSE;}
@@ -78,6 +79,7 @@ private:
    Double_t ComputeCentrality(Bool_t isESD);
    Double_t ComputeMultiplicity(Bool_t isESD,TString type);
    Double_t ApplyCentralityPatchAOD049();
+   Double_t ApplyCentralityPatchPbPb2011();
    void     FillTrueMotherESD(AliRsnMiniEvent *event);
    void     FillTrueMotherAOD(AliRsnMiniEvent *event);
    void     StoreTrueMother(AliRsnMiniPair *pair, AliRsnMiniEvent *event);
@@ -89,6 +91,7 @@ private:
    Bool_t               fUseCentrality;   //  if true, use centrality for event, otherwise use multiplicity
    TString              fCentralityType;  //  definition used to choose what centrality or multiplicity to use
    Bool_t               fUseAOD049CentralityPatch; //flag to enable AOD049 centrality patch
+   Int_t                fUseCentralityPatchPbPb2011; //for PbPb 2011 centrality flattening
 
    Bool_t               fContinuousMix;   //  mixing --> technique chosen (continuous or binned)
    Int_t                fNMix;            //  mixing --> required number of mixes
@@ -123,7 +126,7 @@ private:
    Bool_t              fKeepDfromBOnly;     // flag to keep only the charm particles that comes from beauty decays (specific for D meson analysis)
    Bool_t              fRejectIfNoQuark;    // flag to remove events not generated with PYTHIA
 
-   ClassDef(AliRsnMiniAnalysisTask, 8);   // AliRsnMiniAnalysisTask
+   ClassDef(AliRsnMiniAnalysisTask, 9);   // AliRsnMiniAnalysisTask
 };
 
 
index fa80388..65d01a2 100644 (file)
@@ -20,6 +20,7 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskD0
    Bool_t      isMC,
    Bool_t      isPP,
    Bool_t      ispPb,
+   Bool_t      monitor = kTRUE,
    TString     centrality = "V0M",
    Int_t       aodFilterBit = 5,  
    Float_t     nsigmaTPCPi = 3.0,
@@ -39,9 +40,9 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskD0
    Bool_t      maxDCAcutFixed = kFALSE,
    Bool_t      ptdepPIDcut = kFALSE,
    Bool_t      fixedYcut = kTRUE,
+   Bool_t      checkpileup = kFALSE,
    Bool_t      SPDpileup = kTRUE,
-   Bool_t      checkFeedDown = kTRUE,
-   Bool_t      checkQuark = kTRUE,
+   Bool_t      doCalculationInMC = kTRUE,
    UShort_t    originDselection = 0,
    Int_t       nmix = 5,
    Double_t    minYlab =  -0.5,
@@ -63,8 +64,10 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskD0
   // retrieve analysis manager
   //
   Float_t     cutV = 10.0;
+  Bool_t      checkFeedDown = kTRUE;
+  Bool_t      checkQuark = kTRUE;
   Int_t       aodN = 0;
-  
   
   UInt_t trigger = 0;
   
@@ -130,12 +133,14 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskD0
    // - 3rd argument --> minimum required number of contributors
    // - 4th argument --> tells if TPC stand-alone vertexes must be accepted
    AliRsnCutPrimaryVertex *cutVertex = new AliRsnCutPrimaryVertex("cutVertex", cutV, 0, kFALSE);
-   if(SPDpileup == kTRUE)cutVertex->SetCheckPileUp(kTRUE);
-   AliRsnCutEventUtils *eventUtils = new AliRsnCutEventUtils("cutEventUtils", kFALSE, kFALSE); 
-   if(SPDpileup == kFALSE)eventUtils->SetUseMVPlpSelection(kTRUE); 
+   if(checkpileup == kTRUE){
+       if(SPDpileup == kTRUE)cutVertex->SetCheckPileUp(kTRUE);
+       AliRsnCutEventUtils *eventUtils = new AliRsnCutEventUtils("cutEventUtils", kFALSE, kFALSE); 
+       if(SPDpileup == kFALSE)eventUtils->SetUseMVPlpSelection(kTRUE);
+   } 
                 
-   ::Info("AddAnalysisTaskD0", Form("Checking Pile up? %s", (isPP || ispPb ? "yes" : "no") ));
-   ::Info("AddAnalysisTaskD0", Form("Which Method? %s", (SPDpileup ? "SPD" : "Multi Vertex") ));
+   ::Info("AddAnalysisTaskD0", Form("Checking Pile up? %s", (checkpileup ? "yes" : "no") ));
+   ::Info("AddAnalysisTaskD0", Form("Which Method? %s", (checkpileup? (SPDpileup ? "SPD" : "Multi Vertex") : "none" )));
    
    // define and fill cut set for event cut
    AliRsnCutSet *eventCuts = new AliRsnCutSet("eventCuts", AliRsnTarget::kEvent);
@@ -159,6 +164,12 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskD0
      outMult->AddAxis(multID, 400, 0.0, 400.0);
    else
      outMult->AddAxis(multID, 100, 0.0, 100.0);
+     
+   //tracklets
+   //Int_t trackletID = task->CreateValue(AliRsnMiniValue::kTracklets, kFALSE);
+   //AliRsnMiniOutput *outTracklets = task->CreateOutput("eventTracklets", "HIST", "EVENT");
+   //outTracklets->AddAxis(trackletID, 400, 0.0, 400.0);
+   
    
    //event plane (only for PbPb)
    Int_t planeID = task->CreateValue(AliRsnMiniValue::kPlaneAngle, kFALSE);
@@ -225,7 +236,7 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskD0
        Printf("========================== MC analysis - PID cuts used");
    } else 
      Printf("========================== DATA analysis - PID cuts used");
-   if (!ConfigD0(task, isPP, isMC, nsigmaTPCPi, nsigmaTPCKa, nsigmaTOFPi, nsigmaTOFKa, aodFilterBit, trackDCAcutMax, trackDCAcutMin, trackDCAZcutMax, NTPCcluster, minSPDclt, minpt, maxSisters, checkP,  minDCAcutFixed, maxDCAcutFixed, ptdepPIDcut, checkFeedDown, checkQuark, originDselection, mineta, maxeta, min_inv_mass, max_inv_mass, bins, "", cutsPairY, cutsPair)) return 0x0;
+   if (!ConfigD0(task, isPP, isMC, monitor, nsigmaTPCPi, nsigmaTPCKa, nsigmaTOFPi, nsigmaTOFKa, aodFilterBit, trackDCAcutMax, trackDCAcutMin, trackDCAZcutMax, NTPCcluster, minSPDclt, minpt, maxSisters, checkP,  minDCAcutFixed, maxDCAcutFixed, ptdepPIDcut, checkFeedDown, checkQuark, doCalculationInMC, originDselection, mineta, maxeta, min_inv_mass, max_inv_mass, bins, "", cutsPairY, cutsPair)) return 0x0;
    
    //
    // -- CONTAINERS --------------------------------------------------------------------------------
index b644530..1275b41 100644 (file)
@@ -27,7 +27,7 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskTPCKStar
    Bool_t      enableMonitor = kTRUE,
    Bool_t      IsMcTrueOnly = kFALSE,
    UInt_t      triggerMask = AliVEvent::kMB,
-   //Bool_t      is2011PbPb = kFALSE,
+   Int_t      PbPb2011CentFlat = 0,
    Int_t       nmix = 0,
    Float_t     maxDiffVzMix = 1.0,
    Float_t     maxDiffMultMix = 10.0,
@@ -54,13 +54,12 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskTPCKStar
      Printf(Form("========== SETTING USE CENTRALITY PATCH AOD049 : %s", (aodN==49)? "yes" : "no"));
      task->SetUseCentralityPatch(aodN==49);
    }
-
-   //if(is2011PbPb)
+   //centrality flatening patch LHC11h
+   if(PbPb2011CentFlat)
+     task->SetUseCentralityPatchPbPb2011(PbPb2011CentFlat);
    //task->SelectCollisionCandidates(AliVEvent::kMB | AliVEvent::kCentral | AliVEvent::kSemiCentral);
-   //else
    task->SelectCollisionCandidates(triggerMask);
-
-
+   
    if (isPP) 
      task->UseMultiplicity("QUALITY");
    else
index f7862dc..6f0e7c9 100644 (file)
@@ -12,7 +12,8 @@ Bool_t ConfigD0
 (  
    AliRsnMiniAnalysisTask *task, 
    Bool_t                  isPP,
-   Bool_t                  isMC,  
+   Bool_t                  isMC,
+   Bool_t                  monitor = kTRUE,  
    Float_t                nsigmaTPCPi = 3.0,
    Float_t                nsigmaTPCKa = 3.0,
    Float_t                nsigmaTOFPi = 2.0,
@@ -31,6 +32,7 @@ Bool_t ConfigD0
    Bool_t                  ptdepPIDcut = kFALSE,
    Bool_t                         checkFeedDown = kTRUE,
    Bool_t                         checkQuark = kTRUE,
+   Bool_t                         doCalculationInMC = kTRUE,
    UShort_t                       originDselection = 0,
    Float_t                 mineta = -0.8,
    Float_t                 maxeta = 0.8,
@@ -76,7 +78,7 @@ Bool_t ConfigD0
    cutQuality->SetDCAZmax(trackDCAZcutMax);
    cutQuality->SetSPDminNClusters(minSPDclt);
    cutQuality->SetITSminNClusters(0);
-   cutQuality->SetITSmaxChi2(1E+20);
+   cutQuality->SetITSmaxChi2(36);
    cutQuality->SetTPCmaxChi2(4.0);
    cutQuality->SetRejectKinkDaughters();
    cutQuality->Print();
@@ -110,7 +112,7 @@ Bool_t ConfigD0
    cutQuality->SetDCAZmax(trackDCAZcutMax);
    cutQuality->SetSPDminNClusters(minSPDclt);
    cutQuality->SetITSminNClusters(0);
-   cutQuality->SetITSmaxChi2(1E+20);
+   cutQuality->SetITSmaxChi2(36);
    cutQuality->SetTPCmaxChi2(4.0);
    cutQuality->SetRejectKinkDaughters();
    cutQuality->Print();
@@ -148,17 +150,21 @@ Bool_t ConfigD0
    // [1] = mixing
    // [2] = like ++
    // [3] = like --
-   Bool_t   use     [8] = { 1      ,  1       ,  1      ,  1       ,  1        ,  1        ,  1       ,  1       };
-   Bool_t   useIM   [8] = { 1      ,  1       ,  1      ,  1       ,  1        ,  1        ,  1       ,  1       };
-   TString  name    [8] = {"Unlike1", "Unlike2", "Mixing1", "Mixing2", "RotateK1", "RotateK2", "LikePP" , "LikeMM" };
-   TString  comp    [8] = {"PAIR"   , "PAIR"   , "MIX"  , "MIX"    , "ROTATE1" , "ROTATE1" , "PAIR"   , "PAIR"   };
-   TString  output  [8] = {"SPARSE" , "SPARSE" , "SPARSE" , "SPARSE" , "SPARSE"  , "SPARSE"  , "SPARSE" , "SPARSE" };
-   Char_t   charge1 [8] = {'-'     , '+'      , '-'     , '+'      , '-'       , '+'       , '+'      , '-'      };
-   Char_t   charge2 [8] = {'+'     , '-'      , '+'     , '-'      , '+'       , '-'       , '+'      , '-'      };
-   Int_t    cutID1  [8] = { iCutK   ,  iCutK   ,  iCutK   ,  iCutK   ,  iCutK  ,  iCutK    ,  iCutK   ,  iCutK   };
-   Int_t    cutID2  [8] = { iCutPi  ,  iCutPi  ,  iCutPi  ,  iCutPi  ,  iCutPi ,  iCutPi   ,  iCutPi  ,  iCutPi  };
-   Int_t    ipdg    [8] = { 421     , -421     ,  421   , -421     ,  421      , -421      ,  421     , -421     };
-   Double_t mass    [8] = { 1.86486 ,  1.86486 ,  1.86486 ,  1.86486 ,  1.86486  ,  1.86486  ,  1.86486 ,  1.86486 };
+   
+   
+   if(!isMC || doCalculationInMC == kTRUE){
+   
+   Bool_t   use     [8] = { 1      ,  1       ,  1       ,  1       ,  1        ,  1        ,  1       ,  1      };
+   Bool_t   useIM   [8] = { 1      ,  1       ,  1       ,  1       ,  1        ,  1        ,  1       ,  1      };
+   TString  name    [8] = {"Unlike1", "Unlike2", "Mixing1", "Mixing2", "RotateK1", "RotateK2", "LikePP" , "LikeMM"};
+   TString  comp    [8] = {"PAIR"   , "PAIR"   , "MIX"   , "MIX"    , "ROTATE1" , "ROTATE1" , "PAIR"   , "PAIR"  };
+   TString  output  [8] = {"SPARSE" , "SPARSE" , "SPARSE" , "SPARSE" , "SPARSE"  , "SPARSE"  , "SPARSE" , "SPARSE"};
+   Char_t   charge1 [8] = {'-'     , '+'      , '-'      , '+'      , '-'       , '+'       , '+'      , '-'     };
+   Char_t   charge2 [8] = {'+'     , '-'      , '+'      , '-'      , '+'       , '-'       , '+'      , '-'     };
+   Int_t    cutID1  [8] = { iCutK   ,  iCutK   ,  iCutK   ,  iCutK   ,  iCutK   ,  iCutK    ,  iCutK   ,  iCutK  };
+   Int_t    cutID2  [8] = { iCutPi  ,  iCutPi  ,  iCutPi  ,  iCutPi  ,  iCutPi  ,  iCutPi   ,  iCutPi  ,  iCutPi };
+   Int_t    ipdg    [8] = { 421     , -421     ,  421    , -421     ,  421      , -421      ,  421     , -421    };
+   Double_t mass    [8] = { 1.86486 ,  1.86486 ,  1.86486 ,  1.86486 ,  1.86486  ,  1.86486  ,  1.86486 ,  1.86486};
    
    for (Int_t i = 0; i < 8; i++) {
       if (!use[i]) continue;
@@ -199,7 +205,9 @@ Bool_t ConfigD0
       else out->AddAxis(centID, 400, 0.0, 400.0);
    }
    
+   }
    
+   if(monitor == kTRUE){
    AddMonitorOutput_PionEta(cutSetPi->GetMonitorOutput());
    AddMonitorOutput_PionY(cutSetPi->GetMonitorOutput());
    AddMonitorOutput_PionMinPt(cutSetPi->GetMonitorOutput());
@@ -215,6 +223,8 @@ Bool_t ConfigD0
    AddMonitorOutput_KaonTPC_PIDCut(cutSetK->GetMonitorOutput());
    AddMonitorOutput_KaonTOF_PIDCut(cutSetK->GetMonitorOutput());
    AddMonitorOutput_KaonNTPC(cutSetK->GetMonitorOutput());
+   }
+   
    
    if (isMC) {
    
@@ -277,63 +287,6 @@ Bool_t ConfigD0
    else out->AddAxis(centID, 400, 0.0, 400.0);
    
    
-   // INVARIANT RESOLUTION
-   
-   TString mode = "SPARSE";
-   
-   // create output
-   AliRsnMiniOutput *out = task->CreateOutput("D0_Res1", mode.Data(), "TRUE");
-   // selection settings
-   out->SetDaughter(0, AliRsnDaughter::kKaon);
-   out->SetDaughter(1, AliRsnDaughter::kPion);
-   out->SetCharge(0, '-');
-   out->SetCharge(1, '+');
-   out->SetMotherPDG(421);
-   out->SetMotherMass(1.86486);
-   // pair cuts
-   out->SetPairCuts(cutsPair);
-   out->SetMaxNSisters(maxSisters);
-   out->SetCheckMomentumConservation(checkP);
-   out->SetCheckFeedDown(checkFeedDown);
-   out->SetRejectCandidateIfNotFromQuark(checkQuark);
-   out->SetDselection(originDselection);
-   // binnings
-   out->AddAxis(resID, 200, -0.02, 0.02);
-   out->AddAxis(ptID, 200, 0.0, 20.0);
-   //out->AddAxis(yID, 100, -1, 1);
-   //out->AddAxis(dcapID, 100, -0.001, 0.001);
-   //out->AddAxis(nsistID, 10, 0, 5);
-
-   if (!isPP) out->AddAxis(centID, 100, 0.0, 100.0);
-   else out->AddAxis(centID, 400, 0.0, 400.0);
-   
-   // create output
-   AliRsnMiniOutput *out = task->CreateOutput("D0_Res2", mode.Data(), "TRUE");
-   // selection settings
-   out->SetDaughter(0, AliRsnDaughter::kKaon);
-   out->SetDaughter(1, AliRsnDaughter::kPion);
-   out->SetCharge(0, '+');
-   out->SetCharge(1, '-');
-   out->SetMotherPDG(-421);
-   out->SetMotherMass(1.86486);
-   // pair cuts
-   out->SetPairCuts(cutsPair);
-   out->SetMaxNSisters(maxSisters);
-   out->SetCheckMomentumConservation(checkP);
-   out->SetCheckFeedDown(checkFeedDown);
-   out->SetRejectCandidateIfNotFromQuark(checkQuark);
-   out->SetDselection(originDselection);
-   // binnings
-   out->AddAxis(resID, 200, -0.02, 0.02);
-   out->AddAxis(ptID, 200, 0.0, 20.0);
-   //out->AddAxis(yID, 100, -1, 1);
-   //out->AddAxis(dcapID, 100, -0.001, 0.001);
-   //out->AddAxis(nsistID, 10, 0, 5);
-
-   if (!isPP) out->AddAxis(centID, 100, 0.0, 100.0);
-   else out->AddAxis(centID, 400, 0.0, 400.0);
-   
-   
    // GENERATED MOTHERS
    
    TString mode = "SPARSE";
index 6f635d4..4793a8f 100644 (file)
@@ -97,7 +97,7 @@ Bool_t ConfigTPCanalysisKStar
       out->AddAxis(resID, 200, -0.02, 0.02);
     
     // axis Y: transverse momentum
-    out->AddAxis(ptID, 300, 0.0, 30.0);
+    out->AddAxis(ptID, 200, 0.0, 20.0);
     
     // axis Z: centrality-multiplicity
     if (!isPP)
index 9c97846..61a6b7e 100644 (file)
@@ -99,7 +99,7 @@ Bool_t ConfigTPCanalysisKStarMC
     //out->AddAxis(resID, 200, -0.02, 0.02);
     
     // axis Y: transverse momentum
-    out->AddAxis(ptID, 300, 0.0, 30.0);
+    out->AddAxis(ptID, 200, 0.0, 20.0);
     
     // axis Z: centrality-multiplicity
     if (!isPP)
@@ -127,7 +127,7 @@ Bool_t ConfigTPCanalysisKStarMC
     outm->SetPairCuts(cutsPair);
     // binnings
     outm->AddAxis(imID, 90, 0.6, 1.5);
-    outm->AddAxis(ptID, 300, 0.0, 30.0);
+    outm->AddAxis(ptID, 200, 0.0, 20.0);
     if (!isPP){
     outm->AddAxis(centID, 100, 0.0, 100.0);
     }   else {