//by default this runs locally
//With the argument true this submits jobs to the grid
//As written this requires an xml script tag.xml in the ~/et directory on the grid to submit jobs
-void runHadEt(bool submit = false, bool data = false, Int_t dataset = 20100, Bool_t test = kTRUE, Int_t material = 0) {
+void runHadEt(bool submit = false, bool data = false, Int_t dataset = 20100, Int_t test = 1, Int_t material = 0, Bool_t altV0Scale = kFALSE, bool runCompiledVersion = kFALSE, int simflag = 0) {
TStopwatch timer;
timer.Start();
gSystem->Load("libTree.so");
gSystem->Load("libGeom.so");
gSystem->Load("libVMC.so");
+ gSystem->Load("libPhysics.so");
gSystem->Load("libXMLIO.so");
gSystem->Load("libSTEERBase.so");
gSystem->Load("libAOD.so");
gSystem->Load("libANALYSIS");
+ gSystem->Load("libOADB.so");
gSystem->Load("libANALYSISalice");
+ gSystem->Load("libPWGUDbase.so");
gSystem->AddIncludePath("-I$ALICE_ROOT/include");
- gROOT->ProcessLine(".L AliAnalysisEtCuts.cxx+g");
- gROOT->ProcessLine(".L AliAnalysisHadEtCorrections.cxx+g");
- gROOT->ProcessLine(".L AliAnalysisEtCommon.cxx+g");
- gROOT->ProcessLine(".L AliAnalysisHadEt.cxx+g");
- gROOT->ProcessLine(".L AliAnalysisHadEtMonteCarlo.cxx+g");
- gROOT->ProcessLine(".L AliAnalysisHadEtReconstructed.cxx+g");
- gROOT->ProcessLine(".L AliAnalysisEtSelectionContainer.cxx+g");
- gROOT->ProcessLine(".L AliAnalysisEtSelectionHandler.cxx+g");
- gROOT->ProcessLine(".L AliAnalysisTaskTransverseEnergy.cxx+g");
- gROOT->ProcessLine(".L AliAnalysisTaskHadEt.cxx+g");
-
+ gSystem->AddIncludePath("-I$ALICE_ROOT/PWGUD/base");
+ if(runCompiledVersion){
+ cout<<"Using compiled version"<<endl;
+ gSystem->Load("libPWGLFtotEt.so");
+ }
+ else{
+ cout<<"Not using compiled version"<<endl;
+ gROOT->ProcessLine(".L AliAnalysisEtCuts.cxx+g");
+ gROOT->ProcessLine(".L AliAnalysisHadEtCorrections.cxx+g");
+ gROOT->ProcessLine(".L AliAnalysisEtCommon.cxx+g");
+ gROOT->ProcessLine(".L AliAnalysisHadEt.cxx+g");
+ gROOT->ProcessLine(".L AliAnalysisHadEtMonteCarlo.cxx+g");
+ gROOT->ProcessLine(".L AliAnalysisHadEtReconstructed.cxx+g");
+ gROOT->ProcessLine(".L AliAnalysisTaskTransverseEnergy.cxx+g");
+ gROOT->ProcessLine(".L AliAnalysisTaskHadEt.cxx+g");
+ }
char *kTreeName = "esdTree" ;
TChain * chain = new TChain(kTreeName,"myESDTree") ;
TGrid::Connect("alien://") ;
}
bool PbPb = false;
- if(dataset ==20100){
+ if(dataset ==20100 || dataset==2011){
bool PbPb = true;
if(data){
chain->Add("/data/LHC10h/pass2_rev15/10000137366041.860/AliESDs.root");//Data Pb+Pb
chain->Add("/data/LHC10h/pass2_rev15/10000137366041.900/AliESDs.root");//Data Pb+Pb
}
else{
- chain->Add("/data/LHC10h8/137161/999/AliESDs.root");//Hijing Pb+Pb
- chain->Add("/data/LHC10h8/137161/111/AliESDs.root");//Hijing Pb+Pb
- chain->Add("/data/LHC10h8/137161/222/AliESDs.root");//Hijing Pb+Pb
+ //chain->Add("/data/LHC12d3/168464/201/AliESDs.root");//HIJING with embedded signals - works, full acceptance
+ //chain->Add("/data/LHC10h2/137161/001/AliESDs.root");//DPMJET
+ chain->Add("/data/LHC11a9a/137366/001/AliESDs.root");//AMPT
+// chain->Add("/data/LHC10h8/137161/999/AliESDs.root");//Hijing Pb+Pb
+// chain->Add("/data/LHC10h8/137161/111/AliESDs.root");//Hijing Pb+Pb
+// chain->Add("/data/LHC10h8/137161/222/AliESDs.root");//Hijing Pb+Pb
//chain->Add("/data/LHC11a4_bis/137161/999/AliESDs.root");//Hijing Pb+Pb
//chain->Add("/data/LHC10h12/999/AliESDs.root");//Hijing Pb+Pb
}
}
else{
- if(data){
- chain->Add("/data/LHC10dpass2/10000126403050.70/AliESDs.root");//data
- }
- else{
- chain->Add("/data/LHC10d15/1821/AliESDs.root");//simulation p+p
- }
+ if(dataset==2013){//pPb 5 TeV
+ if(data){
+ chain->Add("/data/LHC13b/13000195483082.95/AliESDs.root");
+ }
+ else{
+ //cout<<"Yes I am analyzing the correct file"<<endl;
+ //chain->Add("/data/LHC12c1b/111/AliESDs.root");
+ chain->Add("/data/LHC13b3/9999/AliESDs.root");
+ }
+ }
+ else{
+ if(data){
+ //chain->Add("/data/LHC10dpass2/10000126403050.70/AliESDs.root");//data
+ //chain->Add("/data/LHC11a/11000146856042.90/AliESDs.root");//data pp 2.76 TeV w/SDD pass 2
+ //chain->Add("/data/LHC11a/11000146860043.90/AliESDs.root");//data pp 2.76 TeV w/SDD pass 3
+ chain->Add("/data/LHC10c/10000120683048.90/AliESDs.root");//Data Pb+Pb
+ }
+ else{
+ cerr<<"Hello I am here 75"<<endl;
+ //chain->Add("/data/LHC10d15/1821/AliESDs.root");//simulation p+p
+ chain->Add("/data/LHC11b10a/001/AliESDs.root");
+ }
+ }
}
// Make the analysis manager
if(submit){
gROOT->LoadMacro("CreateAlienHandlerHadEt.C");
- AliAnalysisGrid *alienHandler = CreateAlienHandlerHadEt(dataset,data,test,material);//integer dataset, boolean isData, bool submit-in-test-mode
+ AliAnalysisGrid *alienHandler = CreateAlienHandlerHadEt(dataset,data,test,material,altV0Scale,runCompiledVersion,simflag);//integer dataset, boolean isData, bool submit-in-test-mode, bool use alternatve V0 scaling
if (!alienHandler) return;
mgr->SetGridHandler(alienHandler);
}
+ cerr<<"Hello I am here 91"<<endl;
AliVEventHandler* esdH = new AliESDInputHandler;
mgr->SetInputEventHandler(esdH);
AliMCEventHandler* handler = new AliMCEventHandler;
mgr->SetMCtruthEventHandler(handler);
}
- AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPhysicsSelection.C");
AliPhysicsSelectionTask *physSelTask = AddTaskPhysicsSelection(!data);
+ if(!physSelTask) { Printf("no physSelTask"); return; }
+ AliPhysicsSelection *physSel = physSelTask->GetPhysicsSelection();
+ //physSel->AddCollisionTriggerClass("+CINT1B-ABCE-NOPF-ALL");// #3119 #769");
+ //physSelTask->AddCollisionTriggerClass("kMB");// #3119 #769");
+
+ cerr<<"Hello I am here 108"<<endl;
+
gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskCentrality.C");
AliCentralitySelectionTask *centTask;
cout<<"Setting up centrality for data"<<endl;
}
}
-
+ Int_t pass = 2;
if(dataset==20100){//PbPb 2.76 TeV
- if(data){gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a4_bis.PbPb.ForData.root","corrections.root",kTRUE);}
- else{gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a4_bis.PbPb.ForSimulations.root","corrections.root",kTRUE);}
+ if(simflag==0){
+ if(data){gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a10a_bis.PbPb.ForData.root","corrections.root",kTRUE);}
+ else{gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a10a_bis.PbPb.ForSimulations.root","corrections.root",kTRUE);}
+ if(material==11){//AMPT only available for this run
+ cout<<"Running AMPT"<<endl;
+ gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a10a_bis.PbPb.137366.ForData.root","corrections.root",kTRUE);
+ }
+ }
+ if(simflag==1){
+ if(data){gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a10a_bis.PbPb.Run138442.ForData.root","corrections.root",kTRUE);}
+ else{gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a10a_bis.PbPb.Run138442.ForSimulations.root","corrections.root",kTRUE);}
+ }
+ if(simflag==2){
+ if(data){gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a10a_bis.PbPb.Run138364.ForData.root","corrections.root",kTRUE);}
+ else{gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a10a_bis.PbPb.Run138364.ForSimulations.root","corrections.root",kTRUE);}
+ }
+ if(simflag==3){
+ if(data){gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a10a_bis.PbPb.Run138396.ForData.root","corrections.root",kTRUE);}
+ else{gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a10a_bis.PbPb.Run138396.ForSimulations.root","corrections.root",kTRUE);}
+ }
+ if(simflag==4){
+ if(data){gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a10a_bis.PbPb.Run137722.ForData.root","corrections.root",kTRUE);}
+ else{gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a10a_bis.PbPb.Run137722.ForSimulations.root","corrections.root",kTRUE);}
+ }
+ if(simflag==5){
+ if(data){gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a10a_bis.PbPb.Run137366.ForData.root","corrections.root",kTRUE);}
+ else{gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a10a_bis.PbPb.Run137366.ForSimulations.root","corrections.root",kTRUE);}
+ }
+ if(simflag==6){
+ if(data){gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a10a_bis.PbPb.Run137161.ForData.root","corrections.root",kTRUE);}
+ else{gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a10a_bis.PbPb.Run137161.ForSimulations.root","corrections.root",kTRUE);}
+ }
gSystem->CopyFile("ConfigHadEtMonteCarloPbPb.C","ConfigHadEtMonteCarlo.C",kTRUE);
gSystem->CopyFile("ConfigHadEtReconstructedPbPb.C","ConfigHadEtReconstructed.C",kTRUE);
+ if(data){gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a10a_bis.PbPb.ForData.root","corrections.root",kTRUE);}
+ else{gSystem->CopyFile("rootFiles/corrections/corrections.LHC11a10a_bis.PbPb.ForSimulations.root","corrections.root",kTRUE);}
+ pass = 2;
//centTask->SetPass(1);
}
else{
+ if(dataset==2011){//PbPb 2.76 TeV 2011
+ gSystem->CopyFile("ConfigHadEtMonteCarloPbPb2011.C","ConfigHadEtMonteCarlo.C",kTRUE);
+ gSystem->CopyFile("ConfigHadEtReconstructedPbPb2011.C","ConfigHadEtReconstructed.C",kTRUE);
+ if(data){gSystem->CopyFile("rootFiles/corrections/corrections.LHC13e1abc.PbPb.168464.ForData.root","corrections.root",kTRUE);}
+ else{gSystem->CopyFile("rootFiles/corrections/corrections.LHC13e1abc.PbPb.168464.ForSimulations.root","corrections.root",kTRUE);}
+ }
if(dataset==2009){//pp 900 GeV
gSystem->CopyFile("ConfigHadEtMonteCarlopp900GeV.C","ConfigHadEtMonteCarlo.C",kTRUE);
gSystem->CopyFile("ConfigHadEtReconstructedpp900GeV.C","ConfigHadEtReconstructed.C",kTRUE);
if(data){gSystem->CopyFile("rootFiles/corrections/corrections.LHC11b1a.pp.ForData.root","corrections.root",kTRUE);}
else{gSystem->CopyFile("rootFiles/corrections/corrections.LHC11b1a.pp.ForSimulations.root","corrections.root",kTRUE);}
+ pass = 3;
}
if(dataset==20111){//pp 2.76 TeV
- gSystem->CopyFile("ConfigHadEtMonteCarlopp276TeV.C","ConfigHadEtMonteCarlo.C",kTRUE);
+ if(altV0Scale)gSystem->CopyFile("ConfigHadEtMonteCarlopp276TeVAlt.C","ConfigHadEtMonteCarlo.C",kTRUE);
+ else{gSystem->CopyFile("ConfigHadEtMonteCarlopp276TeV.C","ConfigHadEtMonteCarlo.C",kTRUE);}
gSystem->CopyFile("ConfigHadEtReconstructedpp276TeV.C","ConfigHadEtReconstructed.C",kTRUE);
if(data){gSystem->CopyFile("rootFiles/corrections/corrections.LHC11b10a.pp.ForData.root","corrections.root",kTRUE);}
else{gSystem->CopyFile("rootFiles/corrections/corrections.LHC11b10a.pp.ForSimulations.root","corrections.root",kTRUE);}
+ pass = 4;
}
if(dataset==2010){//pp 7 TeV
gSystem->CopyFile("ConfigHadEtMonteCarlopp7TeV.C","ConfigHadEtMonteCarlo.C",kTRUE);
gSystem->CopyFile("ConfigHadEtReconstructedpp7TeV.C","ConfigHadEtReconstructed.C",kTRUE);
if(data){gSystem->CopyFile("rootFiles/corrections/corrections.LHC10e20.pp.ForData.root","corrections.root",kTRUE);}
else{gSystem->CopyFile("rootFiles/corrections/corrections.LHC10e20.pp.ForSimulations.root","corrections.root",kTRUE);}
+ pass = 2;
+ }
+ if(dataset==2012){//pp 8 TeV
+ gSystem->CopyFile("ConfigHadEtMonteCarlopp8TeV.C","ConfigHadEtMonteCarlo.C",kTRUE);
+ gSystem->CopyFile("ConfigHadEtReconstructedpp8TeV.C","ConfigHadEtReconstructed.C",kTRUE);
+ if(data){gSystem->CopyFile("rootFiles/corrections/corrections.LHC12c1b.pp.ForData.root","corrections.root",kTRUE);}
+ else{gSystem->CopyFile("rootFiles/corrections/corrections.LHC12c1b.pp.ForSimulations.root","corrections.root",kTRUE);}
+ pass = 1;
+ }
+ if(dataset==2013){//pPb 5 TeV
+ gSystem->CopyFile("ConfigHadEtMonteCarlopPb5TeV.C","ConfigHadEtMonteCarlo.C",kTRUE);
+ gSystem->CopyFile("ConfigHadEtReconstructedpPb5TeV.C","ConfigHadEtReconstructed.C",kTRUE);
+ if(data){gSystem->CopyFile("rootFiles/corrections/corrections.LHC13b3.pPb.ForData.root","corrections.root",kTRUE);}
+ else{gSystem->CopyFile("rootFiles/corrections/corrections.LHC13b3.pPb.ForSimulations.root","corrections.root",kTRUE);}
+ pass = 3;
}
}
+
+
+ gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
+ cerr<<"I am adding PID response task 169"<<endl;
+ //AliAnalysisTask *AddTaskPIDResponse(Bool_t isMC=kFALSE, Bool_t autoMCesd=kTRUE,
+// Bool_t tuneOnData=kFALSE, Int_t recoPass=2,
+// Bool_t cachePID=kFALSE, TString detResponse="",
+// Bool_t useTPCEtaCorrection = kFALSE);
+ AliAnalysisTaskPIDResponse *taskPID=AddTaskPIDResponse(!data,kTRUE,kTRUE,pass);
+ //gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDqa.C");
+ //AddTaskPIDqa();
+
AliAnalysisTaskHadEt *task2 = new AliAnalysisTaskHadEt("TaskHadEt",!data);//,recoFile,mcFile);
if(!data) task2->SetMcData();
+ //Add thing here to select collision type!!
+ if(dataset==2013){//pPb 5 TeV
+ //task2->SelectCollisionCandidates(AliVEvent::kAnyINT ) ;
+ task2->SelectCollisionCandidates(AliVEvent::kINT7 ) ;
+ }
+ else{ if(dataset!=20100){task2->SelectCollisionCandidates(AliVEvent::kMB ) ;}}
+ //if(dataset!=20100){task2->SelectCollisionCandidates(AliVEvent::kMB ) ;}
mgr->AddTask(task2);
+ AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("out2", TList::Class(), AliAnalysisManager::kOutputContainer,"Et.ESD.new.sim.root");
+// mgr->ConnectInput(task2,0,cinput1);
+// mgr->ConnectOutput(task2,1,coutput2);
mgr->ConnectInput(task2,0,cinput1);
mgr->ConnectOutput(task2,1,coutput2);