]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGLF/SPECTRA/PiKaPr/TestAOD/runGrid.C
Merge branch 'feature-movesplit'
[u/mrichter/AliRoot.git] / PWGLF / SPECTRA / PiKaPr / TestAOD / runGrid.C
index aa7e57005ce5797e4ee0434f4641f387bb3a1f17..33ced85d1f15b204fef7a02aa5ab42b255f8f6c6 100644 (file)
@@ -1,7 +1,7 @@
 class  AliAnalysisManager;
 class  AliAnalysisAlien;
 
-void runGrid(TString mode="test",Bool_t mc=0,Int_t day=15,Int_t month=6, Int_t year=2012) 
+void runGrid(TString mode="test",Bool_t mc=0,Int_t day=19,Int_t month=7, Int_t year=2012) 
 {
   //to be used with Aliroot > v5-03-32-AN
   AliLog::SetGlobalDebugLevel(100);
@@ -9,28 +9,28 @@ void runGrid(TString mode="test",Bool_t mc=0,Int_t day=15,Int_t month=6, Int_t y
   gEnv->SetValue("XSec.GSI.DelegProxy", "2");
   
   // Load common libraries
-  gSystem->Load("libCore.so");
-  gSystem->Load("libTree.so");
-  gSystem->Load("libGeom.so");
-  gSystem->Load("libVMC.so");
-  gSystem->Load("libPhysics.so");
-  gSystem->Load("libMinuit.so"); 
-  gSystem->Load("libGui.so");
-  gSystem->Load("libXMLParser.so");
-  gSystem->Load("libSTEERBase.so");
-  gSystem->Load("libESD.so");
-  gSystem->Load("libCDB.so");
+  gSystem->Load("libCore");
+  gSystem->Load("libTree");
+  gSystem->Load("libGeom");
+  gSystem->Load("libVMC");
+  gSystem->Load("libPhysics");
+  gSystem->Load("libMinuit");
+  gSystem->Load("libGui");
+  gSystem->Load("libXMLParser");
+  gSystem->Load("libSTEERBase");
+  gSystem->Load("libESD");
+  gSystem->Load("libCDB");
   gSystem->Load("libAOD");
   gSystem->Load("libANALYSIS");
   gSystem->Load("libANALYSISalice");
   gSystem->Load("libCORRFW");
-  gSystem->Load("libProof.so");
-  gSystem->Load("libRAWDatabase.so");
-  gSystem->Load("libSTEER.so");
+  gSystem->Load("libProof");
+  gSystem->Load("libRAWDatabase");
+  gSystem->Load("libSTEER");
   //__________________________________________________________________________
   // Use AliRoot includes to compile our task
   gROOT->ProcessLine(".include $ALICE_ROOT/include");
-  gSystem->Load("libPWGLFspectra.so");
+  gSystem->Load("libPWGLFspectra");
   
   // Create and configure the alien handler plugin
   AliAnalysisGrid *alienHandler = CreateAlienHandler(mode,mc,day,month,year);  
@@ -54,16 +54,10 @@ void runGrid(TString mode="test",Bool_t mc=0,Int_t day=15,Int_t month=6, Int_t y
   
   gROOT->LoadMacro("$ALICE_ROOT/PWGLF/SPECTRA/PiKaPr/TestAOD/AddTaskSpectraAOD.C");
   //LOOP OVER SELECTION
-  //                            0    1    2    3    4    5
-  Double_t CentCutMin[4]= {     0,  20,  20,  20};
-  Double_t CentCutMax[4]= {   100,  50,  50,  50};
-  Double_t QvecCutMin[4]=    {  0,   0,   0, 1.5};
-  Double_t QvecCutMax[4]=   { 100, 100, 0.4, 100};
-  Double_t EtaMin[4]={       -0.8,-0.8,-0.8,-0.8};
-  Double_t EtaMax[4]={        0.8, 0.8, 0.8, 0.8};
+  //                            0    1    2    3    4    5    6    7    8    9
   Double_t Nsigmapid=3.;
-  Double_t pt=10.;
-  Double_t p=10.;
+  Double_t pt=10000.;
+  Double_t p=10000.;
   Double_t y=.5;
   Double_t ptTofMatch=.6;
   UInt_t trkbit=1024;
@@ -71,23 +65,40 @@ void runGrid(TString mode="test",Bool_t mc=0,Int_t day=15,Int_t month=6, Int_t y
   Bool_t UseCentPatchAOD049=kFALSE;
   Double_t DCA=100000;
   UInt_t minNclsTPC=70;
-  Int_t nrebin=0;
+  Int_t nrebin=30;
   TString opt="";
+  Int_t centcut=30;
   
-  for(Int_t icut=0;icut<4;icut++){
-    if(icut!=0)continue;
-    AliAnalysisTaskSpectraAOD *taskAOD =AddTaskSpectraAOD(mc,CentCutMin[icut],CentCutMax[icut],QvecCutMin[icut],QvecCutMax[icut],EtaMin[icut],EtaMax[icut],Nsigmapid,pt,p,y,ptTofMatch,trkbit,trkbitQVector,UseCentPatchAOD049,DCA,minNclsTPC,nrebin,opt.Data());
-    
+  for(Int_t icut=0;icut<1;icut++){
+    //if(icut!=0)continue;
+    //MB
+    AliAnalysisTaskSpectraAOD *taskAOD =AddTaskSpectraAOD(mc,centcut,centcut+1,0,100,-0.8,0.8,Nsigmapid,pt,p,y,ptTofMatch,trkbit,trkbitQVector,UseCentPatchAOD049,DCA,minNclsTPC,nrebin,opt);
     taskAOD->GetOutputSlot(1)->GetContainer()->SetName(Form("%s_%s",taskAOD->GetOutputSlot(1)->GetContainer()->GetName(),taskAOD->GetName()));
     taskAOD->GetOutputSlot(2)->GetContainer()->SetName(Form("%s_%s",taskAOD->GetOutputSlot(2)->GetContainer()->GetName(),taskAOD->GetName()));
     taskAOD->GetOutputSlot(3)->GetContainer()->SetName(Form("%s_%s",taskAOD->GetOutputSlot(3)->GetContainer()->GetName(),taskAOD->GetName()));
     taskAOD->GetOutputSlot(4)->GetContainer()->SetName(Form("%s_%s",taskAOD->GetOutputSlot(4)->GetContainer()->GetName(),taskAOD->GetName()));
-    
-    // taskAOD->GetEventCuts()->SetMultiplicityCut(800,2000);
-    // taskAOD->GetEventCuts()->SetVertexCut(-8,10);
-    // taskAOD->GetEventCuts()->PrintCuts();
-    
-  }
+    //taskAOD->GetEventCuts()->SetMultiplicityCut(250,450);
+    taskAOD->GetEventCuts()->PrintCuts();
+    //smallq
+    AliAnalysisTaskSpectraAOD *taskAOD =AddTaskSpectraAOD(mc,centcut,centcut+1,0,0.4,-0.8,0.8,Nsigmapid,pt,p,y,ptTofMatch,trkbit,trkbitQVector,UseCentPatchAOD049,DCA,minNclsTPC,nrebin,opt);
+    taskAOD->GetOutputSlot(1)->GetContainer()->SetName(Form("%s_%s",taskAOD->GetOutputSlot(1)->GetContainer()->GetName(),taskAOD->GetName()));
+    taskAOD->GetOutputSlot(2)->GetContainer()->SetName(Form("%s_%s",taskAOD->GetOutputSlot(2)->GetContainer()->GetName(),taskAOD->GetName()));
+    taskAOD->GetOutputSlot(3)->GetContainer()->SetName(Form("%s_%s",taskAOD->GetOutputSlot(3)->GetContainer()->GetName(),taskAOD->GetName()));
+    taskAOD->GetOutputSlot(4)->GetContainer()->SetName(Form("%s_%s",taskAOD->GetOutputSlot(4)->GetContainer()->GetName(),taskAOD->GetName()));
+    //taskAOD->GetEventCuts()->SetMultiplicityCut(250,450);
+    taskAOD->GetEventCuts()->PrintCuts();
+    taskAOD->GetEventCuts()->SetCentFromV0(kTRUE);
+    //largeq
+    AliAnalysisTaskSpectraAOD *taskAOD =AddTaskSpectraAOD(mc,centcut,centcut+1,1.7,100,-0.8,0.8,Nsigmapid,pt,p,y,ptTofMatch,trkbit,trkbitQVector,UseCentPatchAOD049,DCA,minNclsTPC,nrebin,opt);
+    taskAOD->GetOutputSlot(1)->GetContainer()->SetName(Form("%s_%s",taskAOD->GetOutputSlot(1)->GetContainer()->GetName(),taskAOD->GetName()));
+    taskAOD->GetOutputSlot(2)->GetContainer()->SetName(Form("%s_%s",taskAOD->GetOutputSlot(2)->GetContainer()->GetName(),taskAOD->GetName()));
+    taskAOD->GetOutputSlot(3)->GetContainer()->SetName(Form("%s_%s",taskAOD->GetOutputSlot(3)->GetContainer()->GetName(),taskAOD->GetName()));
+    taskAOD->GetOutputSlot(4)->GetContainer()->SetName(Form("%s_%s",taskAOD->GetOutputSlot(4)->GetContainer()->GetName(),taskAOD->GetName()));
+    //taskAOD->GetEventCuts()->SetMultiplicityCut(250,450);
+    taskAOD->GetEventCuts()->PrintCuts();
+    //increase centrality
+    centcut++;
+    }
   
   mgr->SetDebugLevel(2);
   
@@ -102,8 +113,17 @@ void runGrid(TString mode="test",Bool_t mc=0,Int_t day=15,Int_t month=6, Int_t y
 AliAnalysisGrid* CreateAlienHandler(TString mode="test",Bool_t mc=1,Int_t day=0,Int_t month=0, Int_t year=2012)
 {
   AliAnalysisAlien *plugin = new AliAnalysisAlien();
-  plugin->AddIncludePath("-I. -I$ROOTSYS/include -I$ALICE_ROOT/include -I$ALICE_ROOT/TOF -I$ALICE_ROOT/PWGLF");
+  plugin->AddIncludePath("-I. -I$ROOTSYS/include -I$ALICE_ROOT/include -I$ALICE_ROOT/TOF");
   plugin->SetAdditionalLibs("libPWGLFspectra.so");
+  // gROOT->LoadMacro("AliSpectraAODTrackCuts.cxx+g"); 
+  // gROOT->LoadMacro("AliSpectraAODEventCuts.cxx+g"); 
+  // gROOT->LoadMacro("AliSpectraAODHistoManager.cxx+g"); 
+  // gROOT->LoadMacro("AliSpectraAODPID.cxx+g"); 
+  // gROOT->LoadMacro("AliAnalysisTaskSpectraAOD.cxx+g"); 
+  
+  //plugin->SetAnalysisSource("AliSpectraAODTrackCuts.cxx AliSpectraAODEventCuts.cxx AliSpectraAODHistoManager.cxx AliSpectraAODPID.cxx AliAnalysisTaskSpectraAOD.cxx");
+  //plugin->SetAdditionalLibs("Histograms.h HistogramNames.h HistogramNames.cxx AliSpectraAODTrackCuts.cxx AliSpectraAODEventCuts.cxx AliSpectraAODHistoManager.cxx AliSpectraAODPID.cxx AliAnalysisTaskSpectraAOD.cxx AliSpectraAODTrackCuts.h AliSpectraAODEventCuts.h AliSpectraAODHistoManager.h AliSpectraAODPID.h AliAnalysisTaskSpectraAOD.h");
+
   plugin->SetOverwriteMode();
   plugin->SetExecutableCommand("aliroot -q -b");  
   plugin->SetRunMode(mode.Data());
@@ -111,18 +131,27 @@ AliAnalysisGrid* CreateAlienHandler(TString mode="test",Bool_t mc=1,Int_t day=0,
   //Set versions of used packages
   plugin->SetAPIVersion("V1.1x");
   plugin->SetROOTVersion("v5-33-02b");
-  plugin->SetAliROOTVersion("v5-03-32-AN");
+  plugin->SetAliROOTVersion("v5-03-35-AN");
   // Declare input data to be processed.
   if(mc)
     {
-      plugin->SetGridDataDir("/alice/sim/LHC11a10a_bis");
-      plugin->SetDataPattern("AOD090/*AliAOD.root");
+      // plugin->SetGridDataDir("/alice/sim/LHC11a10a_bis");
+      // plugin->SetDataPattern("AOD090/*AliAOD.root");
+      // plugin->SetRunPrefix(""); 
+      // plugin->SetAnalysisMacro("TaskAODPbPbMC.C");
+      // plugin->SetExecutable("TaskAODPbPbMC.sh");
+      // plugin->SetJDLName("TaskAODPbPbMC.jdl");
+      // //plugin->SetSplitMaxInputFileNumber(500);
+      // plugin->SetGridWorkingDir(Form("/AODPbPb%d%d%d/mc/",day,month,year));
+      
+      plugin->SetGridDataDir(" /alice/sim/2012/LHC12a11e/");
+      plugin->SetDataPattern("AOD081/*AliAOD.root");
       plugin->SetRunPrefix(""); 
       plugin->SetAnalysisMacro("TaskAODPbPbMC.C");
       plugin->SetExecutable("TaskAODPbPbMC.sh");
       plugin->SetJDLName("TaskAODPbPbMC.jdl");
-      //plugin->SetSplitMaxInputFileNumber(500);
-      plugin->SetGridWorkingDir(Form("/AODPbPb%d%d%d/mc/",day,month,year));
+      //plugin->SetSplitMaxInputFileNumber(200);
+      plugin->SetGridWorkingDir(Form("/AODPbPb%d%d%d/mcAMPTNoPhysSel_RangeXLarger/",day,month,year));
     }  
   else
     {
@@ -133,7 +162,7 @@ AliAnalysisGrid* CreateAlienHandler(TString mode="test",Bool_t mc=1,Int_t day=0,
       plugin->SetExecutable("TaskAODPbPbdata.sh");
       plugin->SetJDLName("TaskAODPbPbdata.jdl");
       //plugin->SetSplitMaxInputFileNumber(500);
-      plugin->SetGridWorkingDir(Form("/AODPbPb%d%d%d/data/",day,month,year));
+      plugin->SetGridWorkingDir(Form("/AODPbPb%d%d%d/dataFinalPlot3040/",day,month,year));
     }   
   FILE* listruns=fopen("RunListGrid-AOD086.txt","r");
   Int_t irun;
@@ -142,18 +171,25 @@ AliAnalysisGrid* CreateAlienHandler(TString mode="test",Bool_t mc=1,Int_t day=0,
       fscanf(listruns,"%d\n",&irun);
       plugin->AddRunNumber(irun);
     }
+  // plugin->AddRunNumber(139437);
+  // plugin->AddRunNumber(139038);
+  // plugin->AddRunNumber(138653);
+  // plugin->AddRunNumber(138534);
+  // plugin->AddRunNumber(137686);
   // Declare alien output directory. Relative to working directory.
-  plugin->SetGridOutputDir("output"); // In this case will be $HOME/work/output
-  //plugin->SetNrunsPerMaster(3); // 
+    plugin->SetGridOutputDir("output"); // In this case will be $HOME/work/output
+  //plugin->SetGridOutputDir("/alice/cern.ch/user/l/lmilano/AODPbPb472012/data/output/000/");
+  plugin->SetNrunsPerMaster(24); // 
   
   // Declare the output file names separated by blancs.
   // (can be like: file.root or file.root@ALICE::Niham::File)
   //plugin->SetDefaultOutputs(0);
   //plugin->SetOutputFiles("AnalysisResults.root.root");
-  plugin->SetMergeViaJDL(kTRUE);
+  plugin->SetMergeViaJDL(1);
   plugin->SetOneStageMerging(kFALSE);
   plugin->SetMaxMergeStages(2);
-  
+  //plugin->SetMaxMergeFiles(100);
   plugin->SetTTL(100000);
   // Optionally set input format (default xml-single)
   plugin->SetInputFormat("xml-single");