1313// Analysis task to produce smeared pt, eta, phi for electrons/muons in dilepton analysis
1414// Please write to: daiki.sekihata@cern.ch
1515
16- #include < array>
17- #include < string>
18- #include < chrono>
16+ #include " PWGEM/Dilepton/DataModel/dileptonTables.h"
17+ #include " PWGEM/Dilepton/Utils/MomentumSmearer.h"
1918
2019#include " CCDB/BasicCCDBManager.h"
21- #include " Framework/runDataProcessing.h"
22- #include " Framework/AnalysisTask.h"
23- #include " Framework/AnalysisDataModel.h"
24- #include " Framework/ASoAHelpers.h"
2520#include " Framework/ASoA.h"
21+ #include " Framework/ASoAHelpers.h"
22+ #include " Framework/AnalysisDataModel.h"
23+ #include " Framework/AnalysisTask.h"
2624#include " Framework/DataTypes.h"
2725#include " Framework/HistogramRegistry.h"
28- #include " PWGEM/Dilepton/Utils/MomentumSmearer.h"
29- #include " PWGEM/Dilepton/DataModel/dileptonTables.h"
26+ #include " Framework/runDataProcessing.h"
27+
28+ #include < array>
29+ #include < chrono>
30+ #include < string>
3031
3132using namespace o2 ;
3233using namespace o2 ::framework;
@@ -299,7 +300,7 @@ struct ApplySmearing {
299300};
300301
301302struct CheckSmearing {
302- using EMMCParticlesWithSmearing = soa::Join<aod::EMMCParticles, aod::SmearedElectrons, aod::SmearedMuons>;
303+ using EMMCParticlesWithSmearing = soa::Join<aod::EMMCParticles, aod::SmearedElectrons, aod::SmearedMuons>;
303304 using MyCocktailTracks = soa::Join<aod::McParticles, aod::SmearedElectrons, aod::SmearedMuons>;
304305
305306 // Run for electrons or muons
@@ -341,7 +342,6 @@ struct CheckSmearing {
341342
342343 registry.addClone (" Electron/" , " GlobalMuon/" );
343344 registry.addClone (" Electron/" , " StandaloneMuon/" );
344-
345345 }
346346
347347 template <o2::aod::pwgem::dilepton::smearing::EMAnaType type, typename TTracksMC, typename TMCCollisions>
@@ -361,7 +361,7 @@ struct CheckSmearing {
361361
362362 if (std::abs (mctrack.pdgCode ()) == 11 ) { // for electrons
363363 float deltaptoverpt = -1000 .f ;
364- if (mctrack.pt () > 0 .f ){
364+ if (mctrack.pt () > 0 .f ) {
365365 deltaptoverpt = (mctrack.pt () - mctrack.ptSmeared ()) / mctrack.pt ();
366366 }
367367 float deltaeta = mctrack.eta () - mctrack.etaSmeared ();
@@ -379,7 +379,7 @@ struct CheckSmearing {
379379 } else if (std::abs (mctrack.pdgCode ()) == 13 ) { // for muons
380380 float deltaptoverpt = -1000 .f ;
381381 // for standalone muons
382- if (mctrack.pt () > 0 .f ){
382+ if (mctrack.pt () > 0 .f ) {
383383 deltaptoverpt = (mctrack.pt () - mctrack.ptSmeared_sa_muon ()) / mctrack.pt ();
384384 }
385385 float deltaeta = mctrack.eta () - mctrack.etaSmeared_sa_muon ();
@@ -396,7 +396,7 @@ struct CheckSmearing {
396396 registry.fill (HIST (" StandaloneMuon/hCorrelation_Phi" ), mctrack.phi (), mctrack.phiSmeared_sa_muon ());
397397
398398 // for global muons
399- if (mctrack.pt () > 0 .f ){
399+ if (mctrack.pt () > 0 .f ) {
400400 deltaptoverpt = (mctrack.pt () - mctrack.ptSmeared_gl_muon ()) / mctrack.pt ();
401401 }
402402 deltaeta = mctrack.eta () - mctrack.etaSmeared_gl_muon ();
0 commit comments