* provided "as is" without express or implied warranty. *
**************************************************************************/
-
#include <TCanvas.h>
#include <TTree.h>
#include <TFile.h>
extern TROOT *gROOT;
+
ClassImp(AliAnalysisTaskLambdaOverK0sJets)
// Global variables:
noas:
- Double_t pNeg = -100.;
Double_t pPos = -100.;
- Double_t dedxNeg = -1000.;
Double_t dedxPos = -1000.;
Double_t nsigPosPion = 0.;
const AliAODPid *pidPos = ptrack->GetDetPid();
if (pidNeg && pidPos) {
- pNeg = pidNeg->GetTPCmomentum();
pPos = pidPos->GetTPCmomentum();
- dedxNeg = pidNeg->GetTPCsignal()/47.;
dedxPos = pidPos->GetTPCsignal()/47.;
}
isTriggerFromV0Daug = TMath::FloorNint( v0Audit->At(1) );
ptV0LP = v0Audit->At(2);
- if(isV0LP)
- fTriggerIsV0->Fill(ptV0LP);
+ fCheckTriggerFromV0Daug->Fill(isTriggerFromV0Daug);
+
+ if(isV0LP) fTriggerIsV0->Fill(ptV0LP);
+
+ if(isTriggerFromV0Daug && isV0LP)
+ fCheckTriggerFromV0Daug->Fill(2);
if(isTriggerFromV0Daug)
fTriggerComingFromDaug->Fill(ptTrigger);
// -------------------
- fCheckTriggerFromV0Daug->Fill(isTriggerFromV0daug);
-
if(!isTriggerFromV0daug && is2ndLoop)
fCheckTriggerFromV0Daug->Fill(-1);
- else if(isTriggerFromV0daug && isV0LP)
- fCheckTriggerFromV0Daug->Fill(isTriggerFromV0daug+isV0LP);
else if(isTriggerFromV0daug && isV0LP && is2ndLoop)
- fCheckTriggerFromV0Daug->Fill(isTriggerFromV0daug+isV0LP+is2ndLoop);
-
+ fCheckTriggerFromV0Daug->Fill(3);
+
// The highest-pt particle is a V0-candidate
if(isV0LP)
fEvents->Fill(8);
Double_t ptMinTrig = 8.,
Double_t ptMaxTrig = 20.,
Double_t etaMaxTrig = 0.75,
+ Double_t rapMaxV0 = 0.75,
Bool_t sepInjec = kTRUE,
Bool_t isMC = kFALSE,
Bool_t usePID = kTRUE){
task->SetCentrality(minCen,maxCen);
task->SetTriggerPt(ptMinTrig,ptMaxTrig);
task->SetTriggerEta(etaMaxTrig);
+ task->SetMaxY(rapMaxV0);
task->SetSeparateInjectedPart(sepInjec);
task->SetMC(isMC);
task->SetPID(usePID);
if (isMC) name+="_mc";
AliAnalysisDataContainer *coutput1 =
- mgr->CreateContainer("LambdaOverK0sRatio", TList::Class(),
+ mgr->CreateContainer(name, TList::Class(),
AliAnalysisManager::kOutputContainer,
name+".root");
class AliAnalysisAlien;\r
\r
void runGridLambdaOverK0sJets(TString runMode = "full", \r
- TString alirootVer = "v5-03-39-AN",\r
- TString rootVer = "v5-33-02b", \r
+ TString alirootVer = "v5-03-46-AN",\r
+ TString rootVer = "v5-34-01-1",\r
TString dataPath = "ESDs/pass2/AOD049/*/AliAOD.root",\r
TString dataDir = "/alice/data/2010/LHC10h",\r
TString workDir = "work",\r
Double_t ptMinTrig = 8.,\r
Double_t ptMaxTrig = 20.,\r
Double_t etaMaxTrig = 0.75,\r
+ Double_t rapMaxV0 = 0.75,\r
Bool_t sepInjec = kTRUE,\r
Bool_t isMC = kFALSE,\r
Bool_t usePID = kTRUE,\r
AliAnalysisTask *pidTask = AddTaskPIDResponse(isMC);\r
//AliAnalysisTask *pidTask = AddTaskPIDResponse(isMC,kTRUE);\r
if(!pidTask) { printf("no PIDtask\n"); return; }\r
- \r
+ \r
// My task\r
gROOT->LoadMacro("AliAnalysisTaskLambdaOverK0sJets.cxx+g"); \r
gROOT->LoadMacro("AddTaskLambdaOverK0sJets.C");\r
- AliAnalysisTaskLambdaOverK0sJets *task = AddTaskLambdaOverK0sJets(name,minCen,maxCen,ptMinTrig,ptMaxTrig,etaMaxTrig,sepInjec,isMC,usePID);\r
+ AliAnalysisTaskLambdaOverK0sJets *task = AddTaskLambdaOverK0sJets(name,minCen,maxCen,ptMinTrig,ptMaxTrig,etaMaxTrig,rapMaxV0,sepInjec,isMC,usePID);\r
// _____________________________________________________ //\r
\r
if (!mgr->InitAnalysis()) return;\r
gSystem->Load("libRAWDatabase.so");\r
gSystem->Load("libSTEER.so");\r
gSystem->Load("libCORRFW.so");\r
- \r
+\r
+ \r
+ gSystem->AddIncludePath("-I$ALICE_ROOT/include");\r
}\r
\r
// ___________________________________________________________________ //\r
//Based on the file $ALICE_ROOT/PWGLF/STRANGENESS/Cascades/macros/runProofCascadePbPb()\r
void runProofLambdaOverK0sJets(TString proofCluster = "xsanchez@skaf.saske.sk",\r
- TString alirootVer = "VO_ALICE@AliRoot::v5-03-33-AN",\r
- TString rootVer = "VO_ALICE@ROOT::v5-33-02b", \r
+ TString alirootVer = "VO_ALICE@AliRoot::v5-03-46-AN",\r
+ TString rootVer = "VO_ALICE@ROOT::v5-34-01-1", \r
TString path = "/alice/data/LHC10h_000138624_p2_AOD049",\r
TString name = "LambdaOverK0sRatio", \r
Double_t minCen = 0.,\r
Double_t ptMinTrig = 8.,\r
Double_t ptMaxTrig = 20.,\r
Double_t etaMaxTrig = 0.75,\r
+ Double_t rapMaxV0 = 0.75,\r
Bool_t sepInjec = kTRUE,\r
Bool_t isMC = kFALSE,\r
Bool_t usePID = kTRUE){\r
// My task\r
gProof->Load("AliAnalysisTaskLambdaOverK0sJets.cxx+g");\r
gROOT->LoadMacro("AddTaskLambdaOverK0sJets.C");\r
- AliAnalysisTaskLambdaOverK0sJets *task = AddTaskLambdaOverK0sJets(name,minCen,maxCen,ptMinTrig,ptMaxTrig,etaMaxTrig,sepInjec,isMC,usePID);\r
+ AliAnalysisTaskLambdaOverK0sJets *task = AddTaskLambdaOverK0sJets(name,minCen,maxCen,ptMinTrig,ptMaxTrig,etaMaxTrig,rapMaxV0,sepInjec,isMC,usePID);\r
\r
// _____________________________________________________ //\r
\r