Compare MC with the Anchored data period. Only some part of the functionality as defined in the macro:
AliPhysics_SRC/PWGPP/TPC/macros/tpcMCValidation.C+
In [1]:
%jsroot
In [2]:
gSystem->AddIncludePath("-I$ALICE_ROOT/include/"); //couldn't add include path in .rootr
AliDrawStyle::SetDefaults();
AliDrawStyle::ApplyStyle("figTemplate");
TCanvas *canvasDraw = new TCanvas("canvasDraw","canvasDraw",600,600);
In [3]:
.L $AliPhysics_SRC/PWGPP/TPC/macros/tpcMCValidation.C+
In [4]:
InitTPCMCValidation("LHC15k1a1","passMC","LHC15o", "pass3_lowIR_pidfix",0,0); /// function defined in macro
In [5]:
treeMC->GetListOfFriends()->ls();
treeMC->GetFriend("TPC.Anchor")->GetListOfFriends()->ls();
In [6]:
treeMC->GetListOfAliases()->Print("","*resol*Warning");
In [7]:
treeMC->SetAlias("ratio.meanMIP_PhysAcc","abs(ratio.meanMIP-ratio.meanMIP_RobustMean)<0.5"); // phys. acceptable +-0.5
treeMC->SetAlias("ratio.resolutionMIP_PhysAcc","abs(ratio.resolutionMIP-ratio.resolutionMIP_RobustMean)<0.5"); // phys. acceptable +-0.5
In [8]:
MakeStatusPlot("./", "dcarStatusMC.png", "dcar_Warning", "1");
trendingDraw->fWorkingCanvas->Draw();
In [9]:
MakeStatusPlot("./", "dcazStatusMC.png", "dcaz_Warning", "1");
trendingDraw->fWorkingCanvas->Draw();
In [10]:
trendingDraw->MakePlot(outputDir, "interactionRate.png", "Interaction rate", cRange, "",
"Logbook.averageEventsPerSecond;QA.EVS.interactionRate:run", "defaultCut", "figTemplateTRD",
"figTemplateTRD", 1, 1, 4, kTRUE);
trendingDraw->fWorkingCanvas->Draw();
In [11]:
trendingDraw->MakePlot(outputDir, "matchingTPC-ITSEff.png", "Matching efficiency:MC/Anchor", cRange, "",
"QA.TPC.tpcItsMatchA;QA.TPC.tpcItsMatchC;QA.ITS.EffTOTPt02;QA.ITS.EffTOTPt1;QA.ITS.EffTOTPt10;TPC.Anchor.tpcItsMatchA;TPC.Anchor.tpcItsMatchC;ITS.Anchor.EffTOTPt02;ITS.Anchor.EffTOTPt1;ITS.Anchor.EffTOTPt10:run",
"defaultCut", "figTemplateTRDPair", "figTemplateTRDPair", 1, 1.0, 6, kTRUE);
trendingDraw->fWorkingCanvas->Draw();
In [12]:
trendingDraw->MakePlot(outputDir, "dcar_A_0.png", "dcar_A_0", cRange, "",
"dcar_posA_0;TPC.Anchor.dcar_posA_0;dcar_negA_0;TPC.Anchor.dcar_negA_0:run", "defaultCut",
"figTemplateTRDPair", "figTemplateTRDPair", 1, 1.0, 6, kTRUE);
trendingDraw->fWorkingCanvas->Draw();
In [14]:
canvasDraw->cd();
treeMC->Draw("TPC.Anchor.dcar_posA_0:bz","","colz");
canvasDraw->Draw();
In [ ]: