adding parameters to AddTask (Xitzel)
authorbelikov <belikov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 17 Sep 2012 17:04:15 +0000 (17:04 +0000)
committerbelikov <belikov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 17 Sep 2012 17:04:15 +0000 (17:04 +0000)
PWGLF/STRANGENESS/Correlations/AliAnalysisTaskLambdaOverK0sJets.cxx
PWGLF/STRANGENESS/Correlations/macros/AddTaskLambdaOverK0sJets.C
PWGLF/STRANGENESS/Correlations/macros/runGridLambdaOverK0sJets.C
PWGLF/STRANGENESS/Correlations/macros/runProofLambdaOverK0sJets.C

index 8bec480..dd77cff 100644 (file)
@@ -13,7 +13,6 @@
  * provided "as is" without express or implied warranty.                  * 
  **************************************************************************/
 
-
 #include <TCanvas.h>
 #include <TTree.h>
 #include <TFile.h>
@@ -49,6 +48,7 @@
 
 extern TROOT *gROOT;
 
+
 ClassImp(AliAnalysisTaskLambdaOverK0sJets)
 
 // Global variables:
@@ -1583,9 +1583,7 @@ TArrayD* AliAnalysisTaskLambdaOverK0sJets::V0Loop(AliAODTrack *trkTrig, V0LoopSt
     
   noas:
 
-    Double_t pNeg = -100.;
     Double_t pPos = -100.;
-    Double_t dedxNeg = -1000.;
     Double_t dedxPos = -1000.;
 
     Double_t nsigPosPion   = 0.;
@@ -1598,9 +1596,7 @@ TArrayD* AliAnalysisTaskLambdaOverK0sJets::V0Loop(AliAODTrack *trkTrig, V0LoopSt
       const AliAODPid *pidPos = ptrack->GetDetPid();
       
       if (pidNeg && pidPos) {
-       pNeg = pidNeg->GetTPCmomentum();
        pPos = pidPos->GetTPCmomentum();
-       dedxNeg = pidNeg->GetTPCsignal()/47.; 
        dedxPos = pidPos->GetTPCsignal()/47.; 
       }
      
@@ -2129,8 +2125,12 @@ TArrayD* AliAnalysisTaskLambdaOverK0sJets::TriggerParticle()
   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);
@@ -2326,15 +2326,11 @@ void AliAnalysisTaskLambdaOverK0sJets::UserExec(Option_t *)
  
   // -------------------
    
-  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);
index d289df9..76c9613 100644 (file)
@@ -4,6 +4,7 @@ AliAnalysisTaskLambdaOverK0sJets *AddTaskLambdaOverK0sJets( TString  name      =
                                                            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){
@@ -21,6 +22,7 @@ AliAnalysisTaskLambdaOverK0sJets *AddTaskLambdaOverK0sJets( TString  name      =
   task->SetCentrality(minCen,maxCen);
   task->SetTriggerPt(ptMinTrig,ptMaxTrig);
   task->SetTriggerEta(etaMaxTrig);
+  task->SetMaxY(rapMaxV0);
   task->SetSeparateInjectedPart(sepInjec);
   task->SetMC(isMC);
   task->SetPID(usePID);
@@ -38,7 +40,7 @@ AliAnalysisTaskLambdaOverK0sJets *AddTaskLambdaOverK0sJets( TString  name      =
   if (isMC) name+="_mc";
   
   AliAnalysisDataContainer *coutput1 =  
-    mgr->CreateContainer("LambdaOverK0sRatio", TList::Class(), 
+    mgr->CreateContainer(name, TList::Class(), 
                         AliAnalysisManager::kOutputContainer, 
                         name+".root");
   
index 9f266e6..2b551dd 100644 (file)
@@ -1,8 +1,8 @@
 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
@@ -12,6 +12,7 @@ void runGridLambdaOverK0sJets(TString  runMode    = "full",
                              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
@@ -41,11 +42,11 @@ void runGridLambdaOverK0sJets(TString  runMode    = "full",
   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
@@ -81,7 +82,9 @@ void InitAndLoadLibs() {
     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
index 8ca4ace..f91e7bf 100644 (file)
@@ -1,7 +1,7 @@
 //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
@@ -9,6 +9,7 @@ void runProofLambdaOverK0sJets(TString  proofCluster  = "xsanchez@skaf.saske.sk"
                               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
@@ -56,7 +57,7 @@ void runProofLambdaOverK0sJets(TString  proofCluster  = "xsanchez@skaf.saske.sk"
   // 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