Limits of THnSparse multiplcity are now set on the basis of the system pp, PbPb ...
authorprino <prino@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 25 Jun 2012 09:54:30 +0000 (09:54 +0000)
committerprino <prino@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 25 Jun 2012 09:54:30 +0000 (09:54 +0000)
PWGHF/vertexingHF/AliAnalysisTaskSEDplus.cxx
PWGHF/vertexingHF/AliAnalysisTaskSEDplus.h
PWGHF/vertexingHF/macros/AddTaskDplus.C

index 23516f3..44500b5 100644 (file)
@@ -75,7 +75,8 @@ AliAnalysisTaskSE(),
   fNImpParBins(400),
   fLowerImpPar(-2000.),
   fHigherImpPar(2000.),
-  fDoLS(0)
+  fDoLS(0),
+  fSystem(0)
 {
   // Default constructor
 
@@ -152,7 +153,8 @@ AliAnalysisTaskSEDplus::AliAnalysisTaskSEDplus(const char *name,AliRDHFCutsDplus
   fNImpParBins(400),
   fLowerImpPar(-2000.),
   fHigherImpPar(2000.),
-  fDoLS(0)
+  fDoLS(0),
+  fSystem(0)
 {
   // 
   // Standrd constructor
@@ -1157,9 +1159,12 @@ void AliAnalysisTaskSEDplus::CreateImpactParameterHistos(){
   // Histos for impact paramter study
 
   Int_t nmassbins=GetNBinsHistos();
+  Double_t maxmult;
+  if(fSystem==1) maxmult=5000;
+  else maxmult=200;
   Int_t nbins[6]={nmassbins,200,fNImpParBins,5,50,100};
   Double_t xmin[6]={fLowmasslimit,0.,fLowerImpPar,0.95,0.,-0.5};
-  Double_t xmax[6]={fUpmasslimit,20.,fHigherImpPar,1.,1.,99.5};
+  Double_t xmax[6]={fUpmasslimit,40.,fHigherImpPar,1.,1.,maxmult};
 
   fHistMassPtImpParTC[0]=new THnSparseF("hMassPtImpParAll",
                                        "Mass vs. pt vs.imppar - All",
index 191e995..f63d8de 100644 (file)
@@ -41,6 +41,7 @@ class AliAnalysisTaskSEDplus : public AliAnalysisTaskSE
 
   void SetReadMC(Bool_t readMC=kTRUE){fReadMC=readMC;}
   void SetDoLikeSign(Int_t dols=0){fDoLS=dols;}
+  void SetSystem(Int_t system=0){fSystem=system;}
   void SetCutsDistr(Bool_t cutsDistr=kTRUE){fCutsDistr=cutsDistr;}
   void SetDoImpactParameterHistos(Bool_t doImp=kTRUE){fDoImpPar=doImp;}
   void SetImpactParameterBinning(Int_t nbins, Float_t dmin, Float_t dmax){
@@ -143,8 +144,9 @@ class AliAnalysisTaskSEDplus : public AliAnalysisTaskSE
   Float_t fLowerImpPar;  // lower limit in impact parameter (um)
   Float_t fHigherImpPar; // higher limit in impact parameter (um)
   Int_t  fDoLS;        // flag to do LS analysis
+  Int_t fSystem;   //0=pp,1=PbPb
   
-  ClassDef(AliAnalysisTaskSEDplus,17); // AliAnalysisTaskSE for the MC association of heavy-flavour decay candidates
+  ClassDef(AliAnalysisTaskSEDplus,18); // AliAnalysisTaskSE for the MC association of heavy-flavour decay candidates
 };
 
 #endif
index 6256fc7..533330d 100644 (file)
@@ -56,7 +56,7 @@ AliAnalysisTaskSEDplus *AddTaskDplus(Int_t system=0/*0=pp,1=PbPb*/,
   dplusTask->SetDebugLevel(0);
   dplusTask->SetMassLimits(0.2);
   dplusTask->SetUseBit(kTRUE);
-
+  dplusTask->SetSystem(system);
   if (doSparse) dplusTask->SetDoImpactParameterHistos(kTRUE);
 
   mgr->AddTask(dplusTask);