]> git.uio.no Git - u/mrichter/AliRoot.git/commitdiff
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 8bec480828361bd42098791b8fa2db27b2a0688c..dd77cff6b56dc0241cd67b68a63d4e5bdab272df 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 d289df9a63ccdad679567b86741855952f6b6092..76c9613fc1a8bcb5560f05b5f7a053b4c6ba7e54 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 9f266e65e1f84c83f4917b67a18862bddc494c0b..2b551dde71ff8e900e1b0af6c181f3300f38048d 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 8ca4acedf92f73b45b97439bc6d79eddedc871ec..f91e7bf8da63f218559542686581e25cf83a38c6 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