fK0sAssocPtRap = new TH3F("fK0sAssocPtRap","K^{0}_{S} Assoc;p_{T} (GeV/c);y;centrality",nbins,pMin,pMax,30,-1.5,1.5,100,0.,100.);
fOutput->Add(fK0sAssocPtRap);
- fK0sAssocPtMassArm = new TH3F("fK0sAssocPtMassArm","K^{0}_{S} Assoc;Mass (GeV/c^{2})p_{T} (GeV/c);centrality",nbins/2,0.448,0.548,nbins,pMin,pMax,100,0.,100.);
+ fK0sAssocPtMassArm = new TH3F("fK0sAssocPtMassArm","K^{0}_{S} Assoc;Mass (GeV/c^{2})p_{T} (GeV/c);centrality",nbins,0.398,0.598,nbins,pMin,pMax,100,0.,100.);
fOutput->Add(fK0sAssocPtMassArm);
fK0sAssocPtRapEmbeded = new TH3F("fK0sAssocPtRapEmbeded","K^{0}_{S} Assoc - Embeded MC;p_{T} (GeV/c);y;centrality",nbins,pMin,pMax,30,-1.5,1.5,100,0.,100.);
// ****** K0s ******
fK0sMass =
- new TH3F("fK0sMass", "K^{0}_{s}: mass vs p_{T}",nbins/2,0.448,0.548,nbins,pMin,pMax,100,0.,100.);
+ new TH3F("fK0sMass", "K^{0}_{s}: mass vs p_{T}",nbins,0.398,0.598,nbins,pMin,pMax,100,0.,100.);
fK0sMass->GetXaxis()->SetTitle("Mass (GeV/c^2)");
fK0sMass->GetYaxis()->SetTitle("p_{T} (GeV/c)");
fK0sMass->GetZaxis()->SetTitle("centrality");
fK0sMassPtPhi =
new TH3F("fK0sMassPtPhi","K^{0}_{s}: mass vs pt vs #phi;Mass (GeV/c^2);p_{T} (GeV/c);#phi (rad)",
- nbins/2,0.448,0.548,nbins,pMin,pMax,nbinsPhi,0.,2.*TMath::Pi());
+ nbins,0.398,0.598,nbins,pMin,pMax,nbinsPhi,0.,2.*TMath::Pi());
fOutput->Add(fK0sMassPtPhi);
fK0sdPhidEtaPtL[jj*kN1*kNVtxZ + k*kNVtxZ + ll] = new TH3F(hNameHist,"K^{0}_{S}: #Delta#phi vs #Delta#eta vs Inv. Mass",
nbinsdPhi,-TMath::PiOver2(),3*TMath::PiOver2(),
nbinsdEta,-1.5,1.5,
- nbins/2,0.448,0.548);
+ nbins,0.398,0.598);
fK0sdPhidEtaPtL[jj*kN1*kNVtxZ + k*kNVtxZ + ll]->GetXaxis()->SetTitle("#Delta#phi (rad)");
fK0sdPhidEtaPtL[jj*kN1*kNVtxZ + k*kNVtxZ + ll]->GetYaxis()->SetTitle("#Delta#eta");
fK0sdPhidEtaPtL[jj*kN1*kNVtxZ + k*kNVtxZ + ll]->GetZaxis()->SetTitle("Inv. Mass");
// K0s selection
// *******************
if (ctK && (TMath::Abs(v0->RapK0Short())<fYMax) &&
- ( lPtArmV0 > TMath::Abs(0.2*lAlphaV0) ) && ( massK0s > 0.4479 && massK0s < 0.5481 ) ) {
+ ( lPtArmV0 > TMath::Abs(0.2*lAlphaV0) ) && ( massK0s > 0.3979 && massK0s < 0.5981 ) ) {
switch(step) {
case kTriggerCheck:
AliAnalysisTaskLambdaOverK0sJets *AddTaskLambdaOverK0sJets( TString name = "LambdaOverK0sRatio",
TString data = "PbPb2010",
Float_t minCen = 0.,
- Float_t maxCen = 90.,
- Float_t ptMinTrig = 8.,
- Float_t ptMaxTrig = 20.,
- Float_t etaMaxTrig = 0.75,
- Float_t checkIDTrig= kFALSE,
- Float_t rapMaxV0 = 0.75,
+ Float_t maxCen = 40.,
+ Float_t dcaDaug = 1.,
Bool_t sepInjec = kTRUE,
Bool_t isMC = kFALSE,
Bool_t usePID = kTRUE,
Bool_t doQA = kFALSE){
+
+
+ Float_t ptMinTrig = 5.;
+ Float_t ptMaxTrig = 10.;
+ Float_t etaMaxTrig = 0.7;
+ Float_t checkIDTrig = kTRUE;
+ Float_t rapMaxV0 = 0.7;
+
+
AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
if (!mgr) {
::Error("AddTaskCheckCascade", "No analysis manager to connect to.");
task->SetTriggerEta(etaMaxTrig);
task->SetCheckIDTrig(checkIDTrig);
task->SetMaxY(rapMaxV0);
+ task->SetMaxDCADaughter(dcaDaug); // Add to perform systematics
task->SetNSigmaPID(nSigmaPID);
task->SetSeparateInjectedPart(sepInjec);
task->SetMC(isMC);
class AliAnalysisAlien;\r
\r
void runGridLambdaOverK0sJets(TString runMode = "full", \r
- TString alirootVer = "v5-04-59-AN",\r
- TString rootVer = "v5-34-05",\r
+ TString alirootVer = "v5-05-38-AN",\r
+ TString rootVer = "v5-34-11",\r
TString dataPath = "ESDs/pass2/AOD115/*/AliAOD.root",\r
TString dataDir = "/alice/data/2011/LHC11h_2", \r
TString workDir = "test",\r
TString name = "LambdaOverK0sRatio", \r
TString data = "PbPb2011",\r
Float_t minCen = 0.,\r
- Float_t maxCen = 90.,\r
- Float_t ptMinTrig = 5.,\r
- Float_t ptMaxTrig = 10.,\r
- Float_t etaMaxTrig = 0.7,\r
- Float_t rapMaxV0 = 0.7,\r
+ Float_t maxCen = 40.,\r
+ Float_t dcaDaug = 0.5,\r
Bool_t sepInjec = kTRUE,\r
Bool_t isMC = kFALSE,\r
Bool_t usePID = kFALSE,\r
//AliAnalysisTask *pidTask = AddTaskPIDResponse(isMC,kTRUE);\r
if(!pidTask) { printf("no PIDtask\n"); return; }\r
\r
- Float_t checkIDTrig= kTRUE;\r
+ //Float_t checkIDTrig= kTRUE;\r
\r
// My task\r
gROOT->LoadMacro("AliAnalysisTaskLambdaOverK0sJets.cxx++g"); \r
//gSystem->Load("libPWGLFSTRANGENESS");\r
gROOT->LoadMacro("AddTaskLambdaOverK0sJets.C");\r
- AliAnalysisTaskLambdaOverK0sJets *task = AddTaskLambdaOverK0sJets(name,data,minCen,maxCen,ptMinTrig,ptMaxTrig,etaMaxTrig,checkIDTrig,rapMaxV0,sepInjec,isMC,usePID,doQA);\r
+ AliAnalysisTaskLambdaOverK0sJets *task = AddTaskLambdaOverK0sJets(name,data,minCen,maxCen,dcaDaug,sepInjec,isMC,usePID,doQA);\r
// _____________________________________________________ //\r
\r
if (!mgr->InitAnalysis()) return;\r