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);
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);
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;
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);
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());
//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
{
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;
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");