18 #include <TClonesArray.h>
19 #include <TObjString.h>
21 #include <TDatabasePDG.h>
29 #include "TParticle.h"
30 #include "AliVCluster.h"
31 #include "AliAODEvent.h"
32 #include "AliAODMCParticle.h"
33 #include "AliEMCALGeoParams.h"
36 #include "AliEMCALGeometry.h"
48 fMinNCells(0), fMinBadDist(0),
49 fHistoECut(0), fCheckSplitDistToBad(0), fFillAngleHisto(kFALSE),
50 fFillTMHisto(kFALSE), fFillTMResidualHisto(kFALSE), fFillSSExtraHisto(kFALSE),
51 fFillMCHisto(kFALSE), fFillSSWeightHisto(kFALSE),
52 fFillNLMDiffCutHisto(kFALSE), fFillEbinHisto(0),
53 fFillMCOverlapHisto(0), fFillNCellHisto(0), fFillIdConvHisto(0),
54 fFillIdEtaHisto(0), fFillHighMultHisto(0),
55 fFillArmenterosHisto(0), fFillThetaStarHisto(0),
56 fSSWeightN(0), fSSECellCutN(0),
57 fNLMSettingN(0), fWSimu(),
58 fClusterMomentum(), fSubClusterMom1(), fSubClusterMom2(),
59 fSubClusterMomSum(), fSubClusterMomBoost(),
60 fPrimaryMom(), fGrandMotherMom(),
61 fMCDaughMom1(), fMCDaughMom2(), fProdVertex(),
63 fhMassAsyCutNLocMax1(0), fhMassAsyCutNLocMax2(0), fhMassAsyCutNLocMaxN(0),
64 fhM02AsyCutNLocMax1(0), fhM02AsyCutNLocMax2(0), fhM02AsyCutNLocMaxN(0),
65 fhMassM02CutNLocMax1(0), fhMassM02CutNLocMax2(0), fhMassM02CutNLocMaxN(0),
66 fhAsymM02CutNLocMax1(0), fhAsymM02CutNLocMax2(0), fhAsymM02CutNLocMaxN(0),
67 fhMassEnCutNLocMax1(0), fhMassEnCutNLocMax2(0), fhMassEnCutNLocMaxN(0),
68 fhM02EnCutNLocMax1(0), fhM02EnCutNLocMax2(0), fhM02EnCutNLocMaxN(0),
69 fhAsymEnCutNLocMax1(0), fhAsymEnCutNLocMax2(0), fhAsymEnCutNLocMaxN(0),
70 fhSplitEFracEnCutNLocMax1(0), fhSplitEFracEnCutNLocMax2(0), fhSplitEFracEnCutNLocMaxN(0),
71 fhMassSplitECutNLocMax1(0), fhMassSplitECutNLocMax2(0), fhMassSplitECutNLocMaxN(0),
72 fhMCGenFracAfterCutsNLocMax1MCPi0(0), fhMCGenFracAfterCutsNLocMax2MCPi0(0), fhMCGenFracAfterCutsNLocMaxNMCPi0(0),
73 fhMCGenSplitEFracAfterCutsNLocMax1MCPi0(0),fhMCGenSplitEFracAfterCutsNLocMax2MCPi0(0),fhMCGenSplitEFracAfterCutsNLocMaxNMCPi0(0),
74 fhNCellMassEHighNLocMax1MCPi0(0), fhNCellM02EHighNLocMax1MCPi0(0),
75 fhNCellMassELowNLocMax1MCPi0(0), fhNCellM02ELowNLocMax1MCPi0(0),
76 fhNCellMassEHighNLocMax2MCPi0(0), fhNCellM02EHighNLocMax2MCPi0(0),
77 fhNCellMassELowNLocMax2MCPi0(0), fhNCellM02ELowNLocMax2MCPi0(0),
78 fhNCellMassEHighNLocMaxNMCPi0(0), fhNCellM02EHighNLocMaxNMCPi0(0),
79 fhNCellMassELowNLocMaxNMCPi0(0), fhNCellM02ELowNLocMaxNMCPi0(0),
80 fhAnglePairPrimPi0RecoNLocMax1(0), fhAnglePairPrimPi0RecoNLocMax2(0), fhAnglePairPrimPi0RecoNLocMaxN(0),
81 fhAnglePairPrimPi0vsRecoNLocMax1(0), fhAnglePairPrimPi0vsRecoNLocMax2(0), fhAnglePairPrimPi0vsRecoNLocMaxN(0),
82 fhAnglePairPrimPi0OverM02NLocMax1(0), fhAnglePairPrimPi0OverM02NLocMax2(0), fhAnglePairPrimPi0OverM02NLocMaxN(0),
83 fhCentralityPi0NLocMax1(0), fhCentralityEtaNLocMax1(0),
84 fhCentralityPi0NLocMax2(0), fhCentralityEtaNLocMax2(0),
85 fhCentralityPi0NLocMaxN(0), fhCentralityEtaNLocMaxN(0),
86 fhEventPlanePi0NLocMax1(0), fhEventPlaneEtaNLocMax1(0),
87 fhEventPlanePi0NLocMax2(0), fhEventPlaneEtaNLocMax2(0),
88 fhEventPlanePi0NLocMaxN(0), fhEventPlaneEtaNLocMaxN(0),
89 fhClusterEtaPhiNLocMax1(0), fhClusterEtaPhiNLocMax2(0), fhClusterEtaPhiNLocMaxN(0),
90 fhPi0EtaPhiNLocMax1(0), fhPi0EtaPhiNLocMax2(0), fhPi0EtaPhiNLocMaxN(0),
91 fhEtaEtaPhiNLocMax1(0), fhEtaEtaPhiNLocMax2(0), fhEtaEtaPhiNLocMaxN(0),
92 fhPi0EPairDiffTimeNLM1(0), fhPi0EPairDiffTimeNLM2(0), fhPi0EPairDiffTimeNLMN(0),
93 fhEtaEPairDiffTimeNLM1(0), fhEtaEPairDiffTimeNLM2(0), fhEtaEPairDiffTimeNLMN(0),
94 fhMCPi0HighNLMPair(0), fhMCPi0LowNLMPair(0),
95 fhMCPi0AnyNLMPair(0), fhMCPi0NoneNLMPair(0),
96 fhMCPi0HighNLMPairNoMCMatch(0), fhMCPi0LowNLMPairNoMCMatch(0),
97 fhMCPi0AnyNLMPairNoMCMatch(0), fhMCPi0NoneNLMPairNoMCMatch(0),
98 fhMCPi0HighNLMPairOverlap(0), fhMCPi0LowNLMPairOverlap(0),
99 fhMCPi0AnyNLMPairOverlap(0), fhMCPi0NoneNLMPairOverlap(0),
100 fhMCPi0HighNLMPairNoMCMatchOverlap(0), fhMCPi0LowNLMPairNoMCMatchOverlap(0),
101 fhMCPi0AnyNLMPairNoMCMatchOverlap(0), fhMCPi0NoneNLMPairNoMCMatchOverlap(0),
102 fhMCPi0DecayPhotonHitHighLM(0), fhMCPi0DecayPhotonAdjHighLM(0),
103 fhMCPi0DecayPhotonHitOtherLM(0), fhMCPi0DecayPhotonAdjOtherLM(0),
104 fhMCPi0DecayPhotonAdjacent(0), fhMCPi0DecayPhotonHitNoLM(0),
105 fhMCPi0DecayPhotonHitHighLMOverlap(0), fhMCPi0DecayPhotonAdjHighLMOverlap(0),
106 fhMCPi0DecayPhotonHitOtherLMOverlap(0), fhMCPi0DecayPhotonAdjOtherLMOverlap(0),
107 fhMCPi0DecayPhotonAdjacentOverlap(0), fhMCPi0DecayPhotonHitNoLMOverlap(0),
108 fhMCEOverlapType(0), fhMCEOverlapTypeMatch(0)
110 for(Int_t i = 0; i < 7; i++)
112 for(Int_t j = 0; j < 2; j++)
271 for(Int_t jj = 0; jj < 4; jj++)
310 for(Int_t nlm = 0; nlm < 3; nlm++)
372 for(Int_t i = 0; i < 2; i++)
379 for(Int_t i = 0; i < 4; i++)
410 for(Int_t nlm = 0; nlm < 3; nlm++)
478 for(Int_t i = 0; i < 20; i++)
494 for(Int_t i = 0; i < 7; i++)
496 for(Int_t j = 0; j < 4; j++)
512 for(Int_t i = 0; i < 5; i++)
514 for(Int_t j = 0; j < 5; j++)
520 for(Int_t k = 0; k < 3; k++)
542 Float_t e1, Float_t e2, Float_t
mass)
548 const UInt_t nc = cluster->GetNCells();
611 for (UInt_t icell = 0; icell < nc; icell++ )
613 list [icell] = cluster->GetCellAbsId(icell);
614 elist[icell] =
GetEMCALCells()->GetCellAmplitude(list[icell]);
618 Int_t nmaxima = nMax;
619 if(nMax==1) nmaxima = nc ;
626 for(Int_t i = 0; i < nmaxima; i++)
637 for(Int_t i = 0; i < nmaxima; i++)
639 if(i==imax)
continue;
660 Int_t mcLabelMax = -1 ;
661 if(imax >=0 && imax < 999)
667 Int_t mcLabelMax2 = -1 ;
668 if(imax2 >=0 && imax2 < 999)
674 Int_t mcLabelclusterMax = cluster->GetLabels()[0];
675 Bool_t matchHighLMAndHighMC = kFALSE;
679 if(mcLabelclusterMax == mcLabelMax && mcLabelclusterMax >= 0)
681 matchHighLMAndHighMC = kTRUE;
687 if(mcLabelclusterMax == mcLabelMax2 && mcLabelclusterMax >= 0)
690 matchHighLMAndHighMC = kTRUE;
695 matchHighLMAndHighMC = kFALSE;
700 Int_t ancPDG = 0, ancStatus = -1;
702 Bool_t high = kFALSE;
723 for(Int_t i = 0; i < nmaxima-1; i++)
728 for(Int_t j = i+1; j < nmaxima; j++)
733 if(mcLabel1 < 0 || mcLabel2 < 0 )
743 if((i==imax && j==imax2) || (j==imax && i==imax2))
748 else if(ancPDG==22 || TMath::Abs(ancPDG)==11)
760 Int_t
pdg = -22222, status = -1;
766 Float_t en = cluster->E();
773 if(matchHighLMAndHighMC)
790 if(matchHighLMAndHighMC)
815 Int_t
pdg = -22222, status = -1;
818 Int_t label = cluster->GetLabel();
820 while( pdg!=111 && label >=0 )
825 if(pdg!=111 || label < 0)
827 AliWarning(
"Mother Pi0 not found!");
835 AliWarning(Form(
"N daughters %d !=2!",nDaugthers));
840 Int_t pdg0 = -22222, status0 = -1; Int_t label0 = -1;
842 Int_t pdg1 = -22222, status1 = -1; Int_t label1 = -1;
845 if(pdg1!=22 || pdg0 != 22)
847 AliWarning(Form(
"Wrong daughters PDG: photon0 %d - photon1 %d",pdg0,pdg1));
870 if(phi0 < 0 ) phi0+=TMath::TwoPi();
871 if(phi1 < 0 ) phi1+=TMath::TwoPi();
873 Int_t absId0=-1, absId1=-1;
877 if(absId0 < 0 || absId1 < 0)
911 if(inlm > 2) inlm = 2;
913 Bool_t match0 = kFALSE;
914 Bool_t match1 = kFALSE;
917 if(imax >= 0 && imax2 >=0 && absId0 > 0 && absId1 > 0 )
919 if (absId0 == list[imax] ) { match0 = kTRUE ; imatch0 = imax ; }
920 else if(absId0 == list[imax2]) { match0 = kTRUE ; imatch0 = imax2 ; }
922 if (absId1 == list[imax] ) { match1 = kTRUE ; imatch1 = imax ; }
923 else if(absId1 == list[imax2]) { match1 = kTRUE ; imatch1 = imax2 ; }
931 Int_t absId0second = -1;
932 Int_t absId1second = -1;
933 Int_t secLabel0 = -1;
934 Int_t secLabel1 = -1;
938 Int_t secpdg = -999999;
939 Int_t secstatus = -1;
940 Int_t secgrandLabel = -1;
942 if(match0) { secLabel0 = label0 ; mcLabel0 = label0 ; }
943 if(match1) { secLabel1 = label1 ; mcLabel1 = label1 ; }
945 if((!match0 || !match1) && mcindex ==
kmcPi0Conv)
947 for (UInt_t ilab = 0; ilab < cluster->GetNLabels(); ilab++ )
949 Int_t mclabel = cluster->GetLabels()[ilab];
953 if(mclabel == label0 || mclabel == label1)
956 if(mclabel == label0 && secLabel0 < 0) { secLabel0 = label0 ; mcLabel0 = label0 ; }
957 if(mclabel == label1 && secLabel1 < 0) { secLabel1 = label1 ; mcLabel1 = label1 ; }
964 Int_t tmplabel = mclabel;
965 while((secLabel0 < 0 || secLabel1 < 0) && tmplabel > 0 )
971 if((secgrandLabel == label0) || (secgrandLabel == label1 ))
974 if(!match0 && mcLabel1 != secgrandLabel) { secLabel0 = mclabel; mcLabel0 = secgrandLabel; }
975 if(!match1 && mcLabel0 != secgrandLabel) { secLabel1 = mclabel; mcLabel1 = secgrandLabel; }
980 tmplabel = secgrandLabel;
985 if(!match0 && secLabel0 > 0)
996 if(absId0second == list[imax] ) { match0 = kTRUE ; imatch0 = imax ; }
997 if(absId0second == list[imax2]) { match0 = kTRUE ; imatch0 = imax2 ; }
1000 if(!match1 && secLabel1 > 0)
1011 if(absId1second == list[imax] ) { match1 = kTRUE ; imatch1 = imax ; }
1012 if(absId1second == list[imax2]) { match1 = kTRUE ; imatch1 = imax2 ; }
1019 if( match0 && match1 )
1028 if(match0 && imatch0 == imax)
1055 if(match0 && imatch0 == imax )
1108 if(match0 && match1)
1118 if(match0 && imatch0 == imax)
1145 if(match0 && imatch0 == imax)
1196 Bool_t matchMCHitOtherLM = kFALSE;
1199 for(Int_t i = 0; i < nmaxima; i++)
1201 if(imax!=i && imax2!=i && absId1 == list[i]) { match1 = kTRUE; matchMCHitOtherLM = kTRUE; }
1207 for(Int_t i = 0; i < nmaxima; i++)
1209 if(imax!=i && imax2!=i && absId0 == list[i]) { match0 = kTRUE; matchMCHitOtherLM = kTRUE; }
1213 if(matchMCHitOtherLM)
1222 if(match0 && imatch0 == imax)
1241 if(match0 && imatch0 == imax)
1262 Bool_t adjacentOther1 = kFALSE;
1265 for(Int_t i = 0; i < nmaxima; i++)
1267 Int_t ieta=-1; Int_t iphi = 0; Int_t rcu = 0;
1278 Bool_t adjacentOther0 = kFALSE;
1281 for(Int_t i = 0; i < nmaxima; i++)
1283 Int_t ieta=-1; Int_t iphi = 0; Int_t rcu = 0;
1294 if((match0 && adjacentOther1) || (match1 && adjacentOther0))
1303 if(match0 && imatch0 == imax)
1322 if(match0 && imatch0 == imax)
1359 Float_t en, Float_t e1, Float_t e2,
1360 Float_t angle, Float_t
mass,
1361 Float_t anglePrim, Float_t m02,
1362 Float_t asym, Int_t pid, Int_t noverlaps)
1369 Bool_t eCutOK= kFALSE;
1370 Int_t inlm = nMax-1;
1371 if(inlm > 2 ) inlm = 2;
1373 if (ensubcut > 0.1 && ensubcut < e1 && ensubcut < e2 ) eCutOK = kTRUE;
1374 else if(ensubcut < 0.1) eCutOK = kTRUE;
1380 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1401 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1422 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1440 if(
IsDataMC() && mcIndex > 0 && mcIndex < 7)
1450 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1461 if((mcIndex ==
kmcPi0 || mcIndex ==
kmcPi0Conv) && !matched && anglePrim > 0)
1477 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1488 if((mcIndex ==
kmcPi0 || mcIndex ==
kmcPi0Conv) && !matched && anglePrim > 0)
1503 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1514 if((mcIndex ==
kmcPi0 || mcIndex ==
kmcPi0Conv) && !matched && anglePrim > 0)
1528 Float_t en, Float_t m02, Int_t pid)
1533 Float_t momentumDaughter1AlongMother = 0.;
1534 Float_t momentumDaughter2AlongMother = 0.;
1536 if (momentumSquaredMother > 0.)
1543 Float_t ptArmSquared = momentumSquaredDaughter1 - momentumDaughter1AlongMother*momentumDaughter1AlongMother;
1546 if (ptArmSquared > 0.)
1547 pTArm = sqrt(ptArmSquared);
1549 Float_t alphaArm = 0.;
1550 if(momentumDaughter1AlongMother +momentumDaughter2AlongMother > 0)
1551 alphaArm = (momentumDaughter1AlongMother -momentumDaughter2AlongMother) / (momentumDaughter1AlongMother + momentumDaughter2AlongMother);
1555 AliDebug(2,Form(
"E %f, alphaArm %f, pTArm %f",en,alphaArm,pTArm));
1562 Bool_t eCutOK= kFALSE;
1563 Int_t inlm = nMax-1;
1564 if(inlm > 2 ) inlm = 2;
1567 if (ensubcut > 0.1 && ensubcut <
fSubClusterMom1.E() && ensubcut < fSubClusterMom2.E() )
1569 else if(ensubcut < 0.1)
1575 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1583 if((m02OK && asyOK) && (asyOn || m02On))
1591 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1597 if(
IsDataMC() && mcIndex > 0 && mcIndex < 7)
1602 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1610 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1618 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1630 Float_t en, Float_t m02, Int_t pid)
1645 Bool_t eCutOK= kFALSE;
1646 Int_t inlm = nMax-1;
1647 if(inlm > 2 ) inlm = 2;
1649 if (ensubcut > 0.1 && ensubcut <
fSubClusterMom1.E() && ensubcut < fSubClusterMom2.E() ) eCutOK = kTRUE;
1650 else if(ensubcut < 0.1) eCutOK = kTRUE;
1658 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1667 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1676 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1682 if(
IsDataMC() && mcIndex > 0 && mcIndex < 7)
1688 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1697 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1706 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1718 Float_t splitFrac, Float_t
mass, Float_t asym, Float_t l0)
1750 Int_t nMax, Float_t
mass, Float_t l0,
1751 Float_t eta, Float_t phi,
1752 Bool_t matched, Int_t mcindex)
1754 Float_t splitFrac = (e1+e2)/en;
1757 if(e1+e2>0) asym = (e1-e2)/(e1+e2);
1766 if(
IsDataMC() && mcindex > 0 && mcindex < 7)
1781 if(
IsDataMC() && mcindex > 0 && mcindex < 7)
1796 else if( nMax == 2 )
1801 if(
IsDataMC() && mcindex > 0 && mcindex < 7)
1816 else if( nMax >= 3 )
1821 if(
IsDataMC() && mcindex > 0 && mcindex < 7)
1843 Float_t e1, Float_t e2,
1844 Int_t nMax, Float_t
mass, Float_t l0,
1845 Bool_t matched, Int_t mcindex)
1847 Float_t efrac = eprim/en;
1848 Float_t efracSplit = 0;
1849 if(e1+e2 > 0) efracSplit = eprim/(e1+e2);
1851 Float_t splitFrac = (e1+e2)/en;
1854 if(e1+e2>0) asym = (e1-e2)/(e1+e2);
1856 Int_t inlm = nMax-1;
1857 if(inlm > 2) inlm = 2;
1865 Bool_t eCutOK = kFALSE;
1867 if (ensubcut > 0.1 && ensubcut < e1 && ensubcut < e2 ) eCutOK = kTRUE;
1868 else if(ensubcut < 0.1) eCutOK = kTRUE;
1878 if(
IsDataMC() && mcindex > 0 && mcindex < 7)
1894 if( !matched && asyOK && asyOn )
1900 if( !matched && m02OK && m02On )
1907 if(!matched && eCutOK && ensubcut > 0.1)
1915 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1918 if(splitFrac > splitFracMin)
1930 else if( nMax == 2 )
1938 if( !matched && asyOK && asyOn )
1944 if( !matched && m02OK && m02On )
1951 if( !matched && eCutOK && ensubcut > 0.1 )
1959 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
1962 if(splitFrac > splitFracMin)
1983 if( !matched && asyOK && asyOn )
1989 if(!matched && m02OK && m02On )
1996 if( !matched && eCutOK && ensubcut > 0.1 )
2004 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
2007 if(splitFrac > splitFracMin)
2021 if(
IsDataMC() && mcindex > 0 && mcindex < 7)
2029 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
2032 if(splitFrac > splitFracMin)
2045 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
2048 if(splitFrac > splitFracMin)
2061 if((m02OK && asyOK) && (asyOn || m02On) && eCutOK)
2064 if(splitFrac > splitFracMin)
2078 Int_t nc, Int_t nMax, Float_t t12diff,
2079 Float_t
mass, Float_t l0,
2080 Float_t eta, Float_t phi,
2081 Bool_t matched, Int_t mcindex)
2084 if(e1+e2>0) asym = (e1-e2)/(e1+e2);
2093 if(
IsDataMC() && mcindex > 0 && mcindex < 7)
2159 if(
IsDataMC() && mcindex > 0 && mcindex < 7)
2196 Int_t nc, Int_t nMax, Float_t t12diff,
2197 Float_t
mass, Float_t l0,
2198 Float_t eta, Float_t phi,
2199 Bool_t matched, Int_t mcindex)
2202 if(e1+e2>0) asym = (e1-e2)/(e1+e2);
2259 if(
IsDataMC() && mcindex > 0 && mcindex < 7)
2290 Float_t
mass, Float_t l0,
2291 Bool_t matched, Int_t mcindex)
2299 else if ( nMax == 2 )
2305 else if ( nMax > 2 )
2312 if(
IsDataMC() && mcindex > 0 && mcindex < 7)
2320 else if ( nMax == 2 )
2326 else if ( nMax > 2 )
2339 Int_t ebin, Int_t mcindex,Int_t noverlaps,
2340 Float_t l0, Float_t
mass,
2341 Int_t nMax, Bool_t matched,
2342 Float_t splitFrac, Float_t asym,
2343 Float_t eprim, Float_t asymGen)
2345 Float_t efrac = eprim/en;
2346 Float_t efracSplit = 0;
2347 if(e1+e2 > 0) efracSplit = eprim/(e1+e2);
2348 Float_t asymDiff = TMath::Abs(asym) - TMath::Abs(asymGen);
2364 if(asym > 0 && !matched)
2393 else if( nMax == 2 )
2399 if(asym > 0 && !matched)
2433 if(asym > 0 && !matched)
2468 Int_t nc, Float_t
mass, Float_t l0,
2469 Float_t asym, Float_t splitFrac,
2470 Int_t inlm, Int_t ebin, Bool_t matched,
2471 Int_t mcindex, Int_t noverlaps)
2491 else if(noverlaps == 1)
2501 else if(noverlaps > 1)
2512 AliWarning(Form(
"n overlaps = %d!!", noverlaps));
2528 else if(noverlaps == 1)
2538 else if(noverlaps > 1)
2549 AliWarning(Form(
"n overlaps in matched = %d!!", noverlaps));
2558 Bool_t matched, Int_t mcindex,
2559 Float_t
mass , Float_t l0)
2567 if (mcindex==
kmcPi0 && !matched)
2581 else if( nMax == 2 )
2587 if (mcindex==
kmcPi0 && !matched)
2601 else if( nMax >= 3 )
2606 if (mcindex==
kmcPi0 && !matched)
2627 Float_t
energy = clus->E();
2628 Float_t m02 = clus->GetM02();
2634 Double_t
mass = 0., angle = 0.;
2635 Int_t absId1 =-1; Int_t absId2 =-1;
2636 Float_t distbad1 =-1; Float_t distbad2 =-1;
2637 Bool_t fidcut1 = 0; Bool_t fidcut2 = 0;
2655 distbad1,distbad2,fidcut1,fidcut2);
2658 AliWarning(
"No local maximum found! It did not pass CaloPID selection criteria");
2663 if(inlm>2) inlm = 2;
2687 Bool_t matched, Int_t mcindex,
2688 Float_t
mass , Int_t ebin)
2690 Float_t en = cluster->E();
2693 Float_t ll0 = 0., ll1 = 0.;
2694 Float_t disp= 0., dispEta = 0., dispPhi = 0.;
2695 Float_t sEta = 0., sPhi = 0., sEtaPhi = 0.;
2698 ll0, ll1, disp, dispEta, dispPhi, sEta, sPhi, sEtaPhi);
2700 Float_t dispAsy = -1;
2701 if(dispEta+dispPhi >0 ) dispAsy = (dispPhi-dispEta) / (dispPhi+dispEta);
2711 if(
IsDataMC() && mcindex > 0 && mcindex < 7)
2726 else if( nMax == 2 )
2734 if(
IsDataMC() && mcindex > 0 && mcindex < 7)
2749 else if( nMax >= 3 )
2757 if(
IsDataMC() && mcindex > 0 && mcindex < 7)
2778 Int_t absId1, Int_t absId2)
2780 AliVCaloCells* cells = 0;
2787 Float_t simuTotWeight = 0;
2796 AliWarning(Form(
"Wrong calculated energy %f",energy));
2801 Float_t amp1 = cells->GetCellAmplitude(absId1);
2803 Float_t amp2 = cells->GetCellAmplitude(absId2);
2806 if(amp1 < amp2) AliWarning(Form(
"Bad local maxima E ordering : id1 E %f, id2 E %f",amp1,amp2));
2807 if(amp1==0 || amp2==0) AliWarning(Form(
"Null E local maxima : id1 E %f, id2 E %f " ,amp1,amp2));
2820 for (Int_t ipos = 0; ipos < clus->GetNCells(); ipos++)
2822 Int_t
id = clus->GetCellsAbsId()[ipos];
2825 Float_t amp = cells->GetCellAmplitude(
id);
2838 if (
id!=absId1 &&
id!=absId2)
2848 Float_t l0org = clus->GetM02();
2849 Float_t l1org = clus->GetM20();
2850 Float_t dorg = clus->GetDispersion();
2861 Float_t l0 = 0., l1 = 0.;
2862 Float_t disp = 0., dEta = 0., dPhi = 0.;
2863 Float_t sEta = 0., sPhi = 0., sEtaPhi = 0.;
2876 clus->SetM02(l0org);
2877 clus->SetM20(l1org);
2878 clus->SetDispersion(dorg);
2883 Float_t l0 = 0., l1 = 0.;
2884 Float_t disp = 0., dEta = 0., dPhi = 0.;
2885 Float_t sEta = 0., sPhi = 0., sEtaPhi = 0.;
2888 dEta, dPhi, sEta, sPhi, sEtaPhi,
fSSECellCut[iec]);
2901 Int_t nMax, Int_t mcindex)
2903 Float_t dZ = cluster->GetTrackDz();
2904 Float_t dR = cluster->GetTrackDx();
2905 Float_t en = cluster->E();
2915 if(TMath::Abs(dR) < 999)
2922 else if( nMax == 2 )
2927 else if( nMax >= 3 )
2933 if(
IsDataMC() && mcindex > 0 && mcindex < 7)
2940 else if( nMax == 2 )
2945 else if( nMax >= 3 )
2954 Bool_t positive = kFALSE;
2955 if(track) positive = (track->Charge()>0);
2966 else if( nMax == 2 )
2971 else if( nMax >= 3 )
2977 if(
IsDataMC() && mcindex > 0 && mcindex < 7)
2984 else if( nMax == 2 )
2989 else if( nMax >= 3 )
3003 else if( nMax == 2 )
3008 else if( nMax >= 3 )
3014 if(
IsDataMC() && mcindex > 0 && mcindex < 7)
3021 else if( nMax == 2 )
3026 else if( nMax >= 3 )
3043 Int_t buffersize = 255;
3044 char onePar[buffersize] ;
3046 snprintf(onePar,buffersize,
"--- AliAnaInsideClusterInvariantMass ---:") ;
3051 snprintf(onePar,buffersize,
"fNLocMaxCutE =%2.2f;",
GetCaloUtils()->GetLocalMaximaCutE()) ;
3053 snprintf(onePar,buffersize,
"fNLocMaxCutEDiff =%2.2f;",
GetCaloUtils()->GetLocalMaximaCutEDiff()) ;
3055 snprintf(onePar,buffersize,
"fMinNCells =%d;",
fMinNCells) ;
3057 snprintf(onePar,buffersize,
"fMinBadDist =%1.1f;",
fMinBadDist) ;
3065 return new TObjString(parList) ;
3074 TList * outputContainer =
new TList() ;
3075 outputContainer->SetName(
"InsideClusterHistos") ;
3093 Bool_t splitOn = kFALSE;
3094 if(
GetCaloPID()->GetSplitEnergyFractionMinimum(0) > 0 ||
3095 GetCaloPID()->GetSplitEnergyFractionMinimum(1) > 0 ||
3096 GetCaloPID()->GetSplitEnergyFractionMinimum(2) > 0) splitOn = kTRUE;
3098 TString ptype[] ={
"",
"#gamma",
"#gamma->e^{#pm}",
"#pi^{0}",
"#pi^{0} (#gamma->e^{#pm})",
"#eta",
"hadron"};
3099 TString pname[] ={
"",
"Photon",
"Conversion",
"Pi0",
"Pi0Conv",
"Eta",
"Hadron"};
3100 TString snlm [] = {
"1",
"2",
"N"};
3102 TString sEBin[] = {
"8 < #it{E} < 12 GeV",
"12 < #it{E} < 16 GeV",
"16 < #it{E} < 20 GeV",
"#it{E} > 20 GeV" };
3108 Int_t nMaxBins = 10;
3110 TString sMatched[] = {
"",
"Matched"};
3118 for(Int_t imatch = 0; imatch < nMatched; imatch++)
3124 fhNLocMaxDiffCut[iE][iDiff][imatch] =
new TH2F(Form(
"hNLocMax_MinE%d_MinDiffE%d%s",iE, iDiff, sMatched[imatch].
Data()),
3125 Form(
"NLM for #it{E}_{LM}>%1.2f, #Delta E=%1.2F %s",
fNLMMinE[iE],
fNLMMinDiff[iDiff],sMatched[imatch].
Data()),
3126 nptbins,ptmin,ptmax, nMaxBins,0,nMaxBins);
3131 fhNLocMaxDiffCutPi0[iE][iDiff][imatch] =
new TH2F(Form(
"hNLocMaxPi0_MinE%d_MinDiffE%d%s",iE, iDiff, sMatched[imatch].
Data()),
3132 Form(
"#pi^{0} NLM for #it{E}_{LM}>%1.2f, #Delta E=%1.2F %s",
3134 nptbins,ptmin,ptmax, nMaxBins,0,nMaxBins);
3139 for(Int_t inlm = 0; inlm < 3; inlm++)
3142 fhM02NLocMaxDiffCut[iE][iDiff][inlm][imatch] =
new TH2F(Form(
"hNLocMaxM02_MinE%d_MinDiffE%d_NLM%s%s",
3143 iE, iDiff, snlm[inlm].
Data(),sMatched[imatch].
Data()),
3144 Form(
"#lambda^{2}_{0} for #it{E}_{LM}>%1.2f, #Delta E=%1.2F NLM %s %s",
3146 nptbins,ptmin,ptmax, ssbins,ssmin,ssmax);
3151 fhMassNLocMaxDiffCut[iE][iDiff][inlm][imatch] =
new TH2F(Form(
"hNLocMaxMass_MinE%d_MinDiffE%d_NLM%s%s",
3152 iE, iDiff, snlm[inlm].
Data(),sMatched[imatch].
Data()),
3153 Form(
"#it{M}_{split} for #it{E}_{LM}>%1.2f, #Delta E=%1.2F NLM %s %s",
3155 nptbins,ptmin,ptmax, mbins,mmin,mmax);
3160 fhM02NLocMaxDiffCutPi0[iE][iDiff][inlm][imatch] =
new TH2F(Form(
"hNLocMaxPi0M02_MinE%d_MinDiffE%d_NLM%s%s",
3161 iE, iDiff, snlm[inlm].
Data(),sMatched[imatch].
Data()),
3162 Form(
"#pi^{0} #lambda^{2}_{0} for #it{E}_{LM}>%1.2f, #Delta E=%1.2F NLM %s %s",
3164 nptbins,ptmin,ptmax, ssbins,ssmin,ssmax);
3170 iE, iDiff, snlm[inlm].
Data(),sMatched[imatch].
Data()),
3171 Form(
"#pi^{0} #it{M}_{split} for #it{E}_{LM}>%1.2f, #Delta E=%1.2F NLM %s %s",
3173 nptbins,ptmin,ptmax, mbins,mmin,mmax);
3183 return outputContainer;
3188 for(Int_t inlm = 0; inlm < 3; inlm++)
3191 Form(
"Invariant mass of splitted cluster with #it{NLM}=%d vs E, 2nd LM close to bad channel",inlm),
3192 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3198 Form(
"#lambda_{0}^{2} for cluster with #it{NLM}=%d vs E, 2nd LM close to bad channel",inlm),
3199 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3205 Form(
"Invariant mass of splitted cluster with #it{NLM}=%d vs E, 2nd LM close to border",inlm),
3206 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3212 Form(
"#lambda_{0}^{2} for cluster with #it{NLM}=%d vs E, 2nd LM close to border",inlm),
3213 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3220 for(Int_t i = 0; i < n; i++)
3222 for(Int_t j = 0; j < nMatched; j++)
3225 fhNLocMax[i][j] =
new TH2F(Form(
"hNLocMax%s%s",pname[i].
Data(),sMatched[j].
Data()),
3226 Form(
"Number of local maxima in cluster %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3227 nptbins,ptmin,ptmax,nMaxBins,0,nMaxBins);
3228 fhNLocMax[i][j] ->SetYTitle(
"#it{N} maxima");
3229 fhNLocMax[i][j] ->SetXTitle(
"#it{E} (GeV)");
3233 Form(
"Number of local maxima in cluster for split cluster 1 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3234 nptbins,ptmin,ptmax,nMaxBins,0,nMaxBins);
3240 Form(
"Number of local maxima in cluster for split cluster 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3241 nptbins,ptmin,ptmax,nMaxBins,0,nMaxBins);
3249 Form(
"Number of local maxima in cluster %s %s, M02 cut",ptype[i].
Data(),sMatched[j].
Data()),
3250 nptbins,ptmin,ptmax,nMaxBins,0,nMaxBins);
3256 Form(
"Number of local maxima in cluster for split cluster 1 %s %s, M02 cut",ptype[i].
Data(),sMatched[j].
Data()),
3257 nptbins,ptmin,ptmax,nMaxBins,0,nMaxBins);
3263 Form(
"Number of local maxima in cluster for split cluster 2 %s %s, M02 cut",ptype[i].
Data(),sMatched[j].
Data()),
3264 nptbins,ptmin,ptmax,nMaxBins,0,nMaxBins);
3271 Form(
"Number of local maxima in pi0 ID cluster %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3272 nptbins,ptmin,ptmax,nMaxBins,0,nMaxBins);
3279 Form(
"Number of local maxima in cluster for split cluster 1 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3280 nptbins,ptmin,ptmax,nMaxBins,0,nMaxBins);
3286 Form(
"Number of local maxima in cluster for split cluster 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3287 nptbins,ptmin,ptmax,nMaxBins,0,nMaxBins);
3294 Form(
"Number of local maxima vs E of split clusters %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3295 nptbins,ptmin,ptmax,nMaxBins,0,nMaxBins);
3302 Form(
"Number of local maxima vs E of split clusters, id as pi0, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3303 nptbins,ptmin,ptmax,nMaxBins,0,nMaxBins);
3311 Form(
"n cells vs E for N max = 1 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3312 nptbins,ptmin,ptmax,ncbins,ncmin,ncmax);
3318 Form(
"n cells vs E for N max = 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3319 nptbins,ptmin,ptmax,ncbins,ncmin,ncmax);
3326 Form(
"n cells vs E for N max > 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3327 nptbins,ptmin,ptmax,ncbins,ncmin,ncmax);
3334 Form(
"Invariant mass of splitted cluster with #it{NLM}=1 vs E, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3335 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3341 Form(
"Invariant mass of splitted cluster with #it{NLM}=2 vs E, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3342 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3348 Form(
"Invariant mass of splitted cluster with NLM>2 vs E, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3349 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3355 Form(
"Invariant mass of splitted cluster with #it{NLM}=1 vs #it{E}_{1}+#it{E}_{2}, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3356 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3362 Form(
"Invariant mass of splitted cluster with #it{NLM}=2 vs #it{E}_{1}+#it{E}_{2}, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3363 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3369 Form(
"Invariant mass of splitted cluster with NLM>2 vs #it{E}_{1}+#it{E}_{2}, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3370 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3376 Form(
"#lambda_{0}^{2} vs E for N max = 1 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3377 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3383 Form(
"#lambda_{0}^{2} vs E for N max = 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3384 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3390 Form(
"#lambda_{0}^{2} vs E for N max > 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3391 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3397 Form(
"Asymmetry of #it{NLM}=1 vs cluster Energy, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3398 nptbins,ptmin,ptmax,200,-1,1);
3399 fhAsymNLocMax1[i][j]->SetYTitle(
"(#it{E}_{1}-#it{E}_{2})/(#it{E}_{1}+#it{E}_{2})");
3404 Form(
"Asymmetry of #it{NLM}=2 vs cluster Energy, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3405 nptbins,ptmin,ptmax,200,-1,1);
3406 fhAsymNLocMax2[i][j]->SetYTitle(
"(#it{E}_{1}-#it{E}_{2})/(#it{E}_{1}+#it{E}_{2})");
3411 Form(
"Asymmetry of NLM>2 vs cluster Energy, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3412 nptbins,ptmin,ptmax,200,-1,1);
3413 fhAsymNLocMaxN[i][j]->SetYTitle(
"(#it{E}_{1}-#it{E}_{2})/(#it{E}_{1}+#it{E}_{2})");
3418 Form(
"(#it{E}_{1}+#it{E}_{2})/#it{E}_{cluster} vs #it{E}_{cluster} for N max = 1 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3419 nptbins,ptmin,ptmax,120,0,1.2);
3425 Form(
"(#it{E}_{1}+#it{E}_{2})/#it{E}_{cluster} vs #it{E}_{cluster} for N max = 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3426 nptbins,ptmin,ptmax,120,0,1.2);
3432 Form(
"(#it{E}_{1}+#it{E}_{2})/#it{E}_{cluster} vs #it{E}_{cluster} for N max > 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3433 nptbins,ptmin,ptmax,120,0,1.2);
3442 fhMassM02CutNLocMax1 =
new TH2F(
"hMassM02CutNLocMax1",
"Invariant mass of splitted cluster with #it{NLM}=1 vs E, M02 cut, no TM",
3443 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3448 fhMassM02CutNLocMax2 =
new TH2F(
"hMassM02CutNLocMax2",
"Invariant mass of splitted cluster with #it{NLM}=2 vs E, M02 cut, no TM",
3449 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3454 fhMassM02CutNLocMaxN =
new TH2F(
"hMassM02CutNLocMaxN",
"Invariant mass of splitted cluster with NLM>2 vs E, M02 cut, no TM",
3455 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3460 fhAsymM02CutNLocMax1 =
new TH2F(
"hAsymM02CutNLocMax1",
"Asymmetry of #it{NLM}=1 vs cluster Energy, M02Cut, no TM", nptbins,ptmin,ptmax,200,-1,1);
3465 fhAsymM02CutNLocMax2 =
new TH2F(
"hAsymM02CutNLocMax2",
"Asymmetry of #it{NLM}=2 vs cluster Energy, M02Cut, no TM", nptbins,ptmin,ptmax,200,-1,1);
3470 fhAsymM02CutNLocMaxN =
new TH2F(
"hAsymM02CutNLocMaxN",
"Asymmetry of NLM>2 vs cluster Energy, M02Cut, no TM", nptbins,ptmin,ptmax,200,-1,1);
3477 fhMassSplitECutNLocMax1 =
new TH2F(
"hMassSplitECutNLocMax1",
"Invariant mass of splitted cluster with #it{NLM}=1 vs E, (#it{E}_{1}+#it{E}_{2})/E cut, M02 cut, no TM",
3478 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3483 fhMassSplitECutNLocMax2 =
new TH2F(
"hMassSplitECutNLocMax2",
"Invariant mass of splitted cluster with #it{NLM}=2 vs E, (#it{E}_{1}+#it{E}_{2})/E cut, M02 cut, no TM",
3484 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3489 fhMassSplitECutNLocMaxN =
new TH2F(
"hMassSplitECutNLocMaxN",
"Invariant mass of splitted cluster with NLM>2 vs E, (#it{E}_{1}+#it{E}_{2})/E cut, M02 cut, no TM",
3490 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3499 fhMassAsyCutNLocMax1 =
new TH2F(
"hMassAsyCutNLocMax1",
"Invariant mass of splitted cluster with #it{NLM}=1 vs E, Asy cut, no TM",
3500 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3505 fhMassAsyCutNLocMax2 =
new TH2F(
"hMassAsyCutNLocMax2",
"Invariant mass of splitted cluster with #it{NLM}=2 vs E, Asy cut, no TM",
3506 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3511 fhMassAsyCutNLocMaxN =
new TH2F(
"hMassAsyCutNLocMaxN",
"Invariant mass of splitted cluster with NLM>2 vs E, Asy cut, no TM",
3512 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3517 fhM02AsyCutNLocMax1 =
new TH2F(
"hM02AsyCutNLocMax1",
"#lambda_{0}^{2} of #it{NLM}=1 vs cluster Energy, AsyCut, no TM",
3518 nptbins,ptmin,ptmax, ssbins,ssmin,ssmax);
3523 fhM02AsyCutNLocMax2 =
new TH2F(
"hM02AsyCutNLocMax2",
"#lambda_{0}^{2} of #it{NLM}=2 vs cluster Energy, AsyCut, no TM",
3524 nptbins,ptmin,ptmax, ssbins,ssmin,ssmax);
3529 fhM02AsyCutNLocMaxN =
new TH2F(
"hM02AsyCutNLocMaxN",
"#lambda_{0}^{2} of NLM>2 vs cluster Energy, AsyCut, no TM",
3530 nptbins,ptmin,ptmax, ssbins,ssmin,ssmax);
3536 if(
GetCaloPID()->GetSubClusterEnergyMinimum(0) > 0.1)
3538 fhMassEnCutNLocMax1 =
new TH2F(
"hMassEnCutNLocMax1",Form(
"Invariant mass of splitted cluster with #it{NLM}=1 vs E, E > %1.1f GeV, no TM",
GetCaloPID()->GetSubClusterEnergyMinimum(0)),
3539 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3544 fhMassEnCutNLocMax2 =
new TH2F(
"hMassEnCutNLocMax2",Form(
"Invariant mass of splitted cluster with #it{NLM}=2 vs E, E > %1.1f GeV, no TM",
GetCaloPID()->GetSubClusterEnergyMinimum(1)),
3545 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3550 fhMassEnCutNLocMaxN =
new TH2F(
"hMassEnCutNLocMaxN",Form(
"Invariant mass of splitted cluster with NLM>2 vs E, E > %1.1f GeV, no TM",
GetCaloPID()->GetSubClusterEnergyMinimum(2)),
3551 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3556 fhM02EnCutNLocMax1 =
new TH2F(
"hM02EnCutNLocMax1",Form(
"#lambda_{0}^{2} of #it{NLM}=1 vs cluster Energy, E > %1.1f GeV, no TM",
GetCaloPID()->GetSubClusterEnergyMinimum(0)),
3557 nptbins,ptmin,ptmax, ssbins,ssmin,ssmax);
3562 fhM02EnCutNLocMax2 =
new TH2F(
"hM02EnCutNLocMax2",Form(
"#lambda_{0}^{2} of #it{NLM}=2 vs cluster Energy, E > %1.1f GeV, no TM",
GetCaloPID()->GetSubClusterEnergyMinimum(1)),
3563 nptbins,ptmin,ptmax, ssbins,ssmin,ssmax);
3568 fhM02EnCutNLocMaxN =
new TH2F(
"hM02EnCutNLocMaxN",Form(
"#lambda_{0}^{2} of NLM>2 vs cluster Energy, E > %1.1f GeV, no TM",
GetCaloPID()->GetSubClusterEnergyMinimum(2)),
3569 nptbins,ptmin,ptmax, ssbins,ssmin,ssmax);
3574 fhAsymEnCutNLocMax1 =
new TH2F(
"hAsymEnCutNLocMax1",Form(
"Asymmetry of #it{NLM}=1 vs cluster Energy, E > %1.1f GeV, no TM",
GetCaloPID()->GetSubClusterEnergyMinimum(0))
3575 , nptbins,ptmin,ptmax,200,-1,1);
3580 fhAsymEnCutNLocMax2 =
new TH2F(
"hAsymEnCutNLocMax2",Form(
"Asymmetry of #it{NLM}=2 vs cluster Energy, E > %1.1f GeV, no TM",
GetCaloPID()->GetSubClusterEnergyMinimum(1))
3581 , nptbins,ptmin,ptmax,200,-1,1);
3586 fhAsymEnCutNLocMaxN =
new TH2F(
"hAsymEnCutNLocMaxN",Form(
"Asymmetry of NLM>2 vs cluster Energy, E > %1.1f GeV, no TM",
GetCaloPID()->GetSubClusterEnergyMinimum(2))
3587 , nptbins,ptmin,ptmax,200,-1,1);
3592 fhSplitEFracEnCutNLocMax1 =
new TH2F(
"hSplitEFracEnCutNLocMax1",Form(
"SplitEFracmetry of #it{NLM}=1 vs cluster Energy, E > %1.1f GeV, no TM",
GetCaloPID()->GetSubClusterEnergyMinimum(0))
3593 , nptbins,ptmin,ptmax,120,0,1.2);
3598 fhSplitEFracEnCutNLocMax2 =
new TH2F(
"hSplitEFracEnCutNLocMax2",Form(
"SplitEFracmetry of #it{NLM}=2 vs cluster Energy, E > %1.1f GeV, no TM",
GetCaloPID()->GetSubClusterEnergyMinimum(1))
3599 , nptbins,ptmin,ptmax,120,0,1.2);
3604 fhSplitEFracEnCutNLocMaxN =
new TH2F(
"hSplitEFracEnCutNLocMaxN",Form(
"SplitEFracmetry of NLM>2 vs cluster Energy, E > %1.1f GeV, no TM",
GetCaloPID()->GetSubClusterEnergyMinimum(2))
3605 , nptbins,ptmin,ptmax,120,0,1.2);
3615 Form(
"Mass vs E, %s %s, for NLM = 1, M02 and asy cut",ptype[i].
Data(),sMatched[j].
Data()),
3616 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3622 Form(
"Mass vs E, %s %s, for NLM = 2, M02 and asy cut",ptype[i].
Data(),sMatched[j].
Data()),
3623 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3629 Form(
"Mass vs E, %s %s, for NLM > 2, M02 and asy cut",ptype[i].
Data(),sMatched[j].
Data()),
3630 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3636 Form(
"Mass vs #it{E}_{1}+#it{E}_{2}, %s %s, for NLM = 1, M02 and asy cut",ptype[i].
Data(),sMatched[j].
Data()),
3637 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3643 Form(
"Mass vs #it{E}_{1}+#it{E}_{2}, %s %s, for NLM = 2, M02 and asy cut",ptype[i].
Data(),sMatched[j].
Data()),
3644 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3650 Form(
"Mass vs #it{E}_{1}+#it{E}_{2}, %s %s, for NLM > 2, M02 and asy cut",ptype[i].
Data(),sMatched[j].
Data()),
3651 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3658 Form(
"(#it{E}_{1}+#it{E}_{2})/#it{E}_{cluster} vs #it{E}_{cluster} for N max = 1, M02 and Asy cut on, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3659 nptbins,ptmin,ptmax,120,0,1.2);
3665 Form(
"(#it{E}_{1}+#it{E}_{2})/#it{E}_{cluster} vs #it{E}_{cluster} for N max = 2, M02 and Asy cut on, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3666 nptbins,ptmin,ptmax,120,0,1.2);
3672 Form(
"(#it{E}_{1}+#it{E}_{2})/#it{E}_{cluster} vs #it{E}_{cluster} for N max > 2, M02 and Asy cut on, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3673 nptbins,ptmin,ptmax,120,0,1.2);
3680 Form(
"Invariant mass of splitted cluster with #it{NLM}=1, #lambda_{0}^{2}, E > 12 GeV,%s %s",ptype[i].
Data(),sMatched[j].
Data()),
3681 ssbins,ssmin,ssmax,mbins,mmin,mmax);
3687 Form(
"Invariant mass of splitted cluster with #it{NLM}=2, #lambda_{0}^{2}, E > 12 GeV, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3688 ssbins,ssmin,ssmax,mbins,mmin,mmax);
3694 Form(
"Invariant mass of splitted cluster with NLM>2, vs #lambda_{0}^{2}, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3695 ssbins,ssmin,ssmax,mbins,mmin,mmax);
3703 Form(
"Invariant mass of splitted cluster with #it{NLM}=1, #sigma_{#eta #eta}^{2}, E > 12 GeV,%s %s",ptype[i].
Data(),sMatched[j].
Data()),
3704 ssbins,ssmin,ssmax,mbins,mmin,mmax);
3710 Form(
"Invariant mass of splitted cluster with #it{NLM}=2 #sigma_{#eta #eta}^{2}, E > 12 GeV, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3711 ssbins,ssmin,ssmax,mbins,mmin,mmax);
3717 Form(
"Invariant mass of splitted cluster with NLM>2, #sigma_{#eta #eta}^{2}, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3718 ssbins,ssmin,ssmax,mbins,mmin,mmax);
3724 Form(
"Invariant mass of 2 highest energy cells #sigma_{#phi #phi}^{2}, E > 12 GeV,%s %s",ptype[i].
Data(),sMatched[j].
Data()),
3725 ssbins,ssmin,ssmax,mbins,mmin,mmax);
3731 Form(
"Invariant mass of 2 local maxima cells #sigma_{#phi #phi}^{2}, E > 12 GeV, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3732 ssbins,ssmin,ssmax,mbins,mmin,mmax);
3738 Form(
"Invariant mass of N>2 local maxima cells vs #sigma_{#phi #phi}^{2}, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3739 ssbins,ssmin,ssmax,mbins,mmin,mmax);
3745 Form(
"Invariant mass of 2 highest energy cells A = (#sigma_{#phi #phi}^{2} - #sigma_{#eta #eta}^{2}) / (#sigma_{#phi #phi}^{2} + #sigma_{#eta #eta}^{2}), E > 12 GeV,%s %s",ptype[i].
Data(),sMatched[j].
Data()),
3746 200,-1,1,mbins,mmin,mmax);
3748 fhMassDispAsyNLocMax1[i][j]->SetXTitle(
"A = (#sigma_{#phi #phi}^{2} - #sigma_{#eta #eta}^{2}) / (#sigma_{#phi #phi}^{2} + #sigma_{#eta #eta}^{2})");
3752 Form(
"Invariant mass of 2 local maxima cells A = (#sigma_{#phi #phi}^{2} - #sigma_{#eta #eta}^{2}) / (#sigma_{#phi #phi}^{2} + #sigma_{#eta #eta}^{2}), E > 12 GeV, %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3753 200,-1,1,mbins,mmin,mmax);
3755 fhMassDispAsyNLocMax2[i][j]->SetXTitle(
"A = (#sigma_{#phi #phi}^{2} - #sigma_{#eta #eta}^{2}) / (#sigma_{#phi #phi}^{2} + #sigma_{#eta #eta}^{2})");
3759 Form(
"Invariant mass of N>2 local maxima cells vsA = (#sigma_{#phi #phi}^{2} - #sigma_{#eta #eta}^{2}) / (#sigma_{#phi #phi}^{2} + #sigma_{#eta #eta}^{2}), %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3760 200,-1,1,mbins,mmin,mmax);
3762 fhMassDispAsyNLocMaxN[i][j]->SetXTitle(
"#it{A} = (#sigma_{#phi #phi}^{2} - #sigma_{#eta #eta}^{2}) / (#sigma_{#phi #phi}^{2} + #sigma_{#eta #eta}^{2})");
3769 Form(
"#lambda_{0}^{2} vs E for N max = 1 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3770 nptbins,ptmin,ptmax,200,0,2);
3776 Form(
"#lambda_{0}^{2} vs E for N max = 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3777 nptbins,ptmin,ptmax,200,0,2);
3783 Form(
"#lambda_{0}^{2} vs E for N max > 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3784 nptbins,ptmin,ptmax,200,0,2);
3790 Form(
"#lambda_{0}^{2} vs E for N max = 1 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3791 nptbins,ptmin,ptmax,200,0,2);
3797 Form(
"#lambda_{0}^{2} vs E for N max = 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3798 nptbins,ptmin,ptmax,200,0,2);
3804 Form(
"#lambda_{0}^{2} vs E for N max > 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3805 nptbins,ptmin,ptmax,200,0,2);
3812 Form(
"#it{E}_{gen} / (#it{E}_{1 split}+#it{E}_{2 split}) vs E for N max = 1 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3813 nptbins,ptmin,ptmax,200,0,2);
3819 Form(
"#it{E}_{gen} / (#it{E}_{1 split}+#it{E}_{2 split}) vs E for N max = 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3820 nptbins,ptmin,ptmax,200,0,2);
3826 Form(
"#it{E}_{gen} / (#it{E}_{1 split}+#it{E}_{2 split}) vs E for N max > 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3827 nptbins,ptmin,ptmax,200,0,2);
3833 Form(
"#it{E}_{gen} / (#it{E}_{1 split}+#it{E}_{2 split}) vs E for N max = 1 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3834 nptbins,ptmin,ptmax,200,0,2);
3840 Form(
"#it{E}_{gen} / (#it{E}_{1 split}+#it{E}_{2 split}) vs E for N max = 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3841 nptbins,ptmin,ptmax,200,0,2);
3847 Form(
"#it{E}_{gen} / (#it{E}_{1 split}+#it{E}_{2 split}) vs E for N max > 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3848 nptbins,ptmin,ptmax,200,0,2);
3854 Form(
"(#it{E}_{1 split}+#it{E}_{2 split})/#it{E}_{reco} vs #it{E}_{gen} / #it{E}_{reco} for N max = 1 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3861 Form(
"(#it{E}_{1 split}+#it{E}_{2 split})/#it{E}_{reco} vs #it{E}_{gen} / #it{E}_{reco} for N max = 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3869 Form(
"(#it{E}_{1 split}+#it{E}_{2 split})/#it{E}_{reco} vs #it{E}_{gen} / #it{E}_{reco} for N max > 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3877 Form(
"#it{E}_{1 split}+#it{E}_{2 split} vs #it{E}_{gen} for N max = 1 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3878 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3884 Form(
"#it{E}_{1 split}+#it{E}_{2 split} vs #it{E}_{gen} for N max = 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3885 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3892 Form(
"#it{E}_{1 split}+#it{E}_{2 split} vs #it{E}_{gen} for N max > 2 %s %s",ptype[i].
Data(),sMatched[j].
Data()),
3893 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3905 Form(
"#lambda_{0}^{2} vs #it{E}, %s, for NLM = 1",ptype[i].
Data()),
3906 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3912 Form(
"#lambda_{0}^{2} vs #it{E}, %s, for NLM = 2",ptype[i].
Data()),
3913 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3919 Form(
"#lambda_{0}^{2} vs #it{E}, %s, for NLM > 2",ptype[i].
Data()),
3920 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3926 Form(
"Mass vs #it{E}, %s, for NLM = 1",ptype[i].
Data()),
3927 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3933 Form(
"Mass vs #it{E} , %s, for NLM = 2",ptype[i].
Data()),
3934 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3940 Form(
"Mass vs #it{E}, %s, for NLM > 2",ptype[i].
Data()),
3941 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3947 Form(
"Mass vs #it{E}_{1}+#it{E}_{2}, %s, for NLM = 1",ptype[i].
Data()),
3948 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3954 Form(
"Mass vs #it{E}_{1}+#it{E}_{2} , %s, for NLM = 2",ptype[i].
Data()),
3955 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3961 Form(
"Mass vs #it{E}_{1}+#it{E}_{2}, %s, for NLM > 2",ptype[i].
Data()),
3962 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3968 Form(
"Asymmetry vs E, %s, for NLM = 1",ptype[i].
Data()),
3969 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3975 Form(
"Asymmetry vs E, %s, for NLM = 2",ptype[i].
Data()),
3976 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3982 Form(
"Asymmetry vs E, %s, for NLM > 2",ptype[i].
Data()),
3983 nptbins,ptmin,ptmax,mbins,mmin,mmax);
3991 Form(
"n cells vs E, %s, for NLM = 1",ptype[i].
Data()),
3992 nptbins,ptmin,ptmax,ncbins,ncmin,ncmax);
3998 Form(
"n cells vs E, %s, for NLM = 2",ptype[i].
Data()),
3999 nptbins,ptmin,ptmax,ncbins,ncmin,ncmax);
4005 Form(
"n cells vs E, %s, for NLM > 2",ptype[i].
Data()),
4006 nptbins,ptmin,ptmax,ncbins,ncmin,ncmax);
4017 Form(
"#lambda_{0}^{2} vs E, %s, for NLM = 1",ptype[i].
Data()),
4018 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4025 Form(
"#lambda_{0}^{2} vs E, %s, for NLM = 2",ptype[i].
Data()),
4026 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4032 Form(
"#lambda_{0}^{2} vs E, %s, for NLM > 2",ptype[i].
Data()),
4033 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4039 Form(
"Mass vs E, %s, for NLM = 1",ptype[i].
Data()),
4040 nptbins,ptmin,ptmax,mbins,mmin,mmax);
4046 Form(
"Mass vs E, %s, for NLM = 2",ptype[i].
Data()),
4047 nptbins,ptmin,ptmax,mbins,mmin,mmax);
4053 Form(
"Mass vs E, %s, for NLM > 2",ptype[i].
Data()),
4054 nptbins,ptmin,ptmax,mbins,mmin,mmax);
4060 Form(
"Asymmetry vs E, %s, for NLM = 1",ptype[i].
Data()),
4061 nptbins,ptmin,ptmax,mbins,mmin,mmax);
4067 Form(
"Asymmetry vs E, %s, for NLM = 2",ptype[i].
Data()),
4068 nptbins,ptmin,ptmax,mbins,mmin,mmax);
4074 Form(
"Asymmetry vs E, %s, for NLM > 2",ptype[i].
Data()),
4075 nptbins,ptmin,ptmax,mbins,mmin,mmax);
4083 Form(
"n cells vs E, %s, for NLM = 1",ptype[i].
Data()),
4084 nptbins,ptmin,ptmax,ncbins,ncmin,ncmax);
4090 Form(
"n cells vs E, %s, for NLM = 2",ptype[i].
Data()),
4091 nptbins,ptmin,ptmax,ncbins,ncmin,ncmax);
4097 Form(
"n cells vs E, %s, for NLM > 2",ptype[i].
Data()),
4098 nptbins,ptmin,ptmax,ncbins,ncmin,ncmax);
4108 Form(
"#lambda_{0}^{2} vs E, %s, for NLM = 1",ptype[i].
Data()),
4109 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4115 Form(
"#lambda_{0}^{2} vs E, %s, for NLM = 2",ptype[i].
Data()),
4116 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4122 Form(
"#lambda_{0}^{2} vs E, %s, for NLM > 2",ptype[i].
Data()),
4123 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4130 Form(
"Mass vs E, %s, for NLM = 1",ptype[i].
Data()),
4131 nptbins,ptmin,ptmax,mbins,mmin,mmax);
4137 Form(
"Mass vs E, %s, for NLM = 2",ptype[i].
Data()),
4138 nptbins,ptmin,ptmax,mbins,mmin,mmax);
4144 Form(
"Mass vs E, %s, for NLM > 2",ptype[i].
Data()),
4145 nptbins,ptmin,ptmax,mbins,mmin,mmax);
4151 Form(
"Asymmetry vs E, %s, for NLM = 1",ptype[i].
Data()),
4152 nptbins,ptmin,ptmax,mbins,mmin,mmax);
4158 Form(
"Asymmetry vs E, %s, for NLM = 2",ptype[i].
Data()),
4159 nptbins,ptmin,ptmax,mbins,mmin,mmax);
4165 Form(
"Asymmetry vs E, %s, for NLM > 2",ptype[i].
Data()),
4166 nptbins,ptmin,ptmax,mbins,mmin,mmax);
4175 for(Int_t j = 0; j < 4; j++)
4178 Form(
"Invariant mass of 2 highest energy cells vs (#it{E}_{1}+#it{E}_{2})/Ecluster, %s, %s",ptype[i].
Data(),sEBin[j].
Data()),
4179 120,0,1.2,mbins,mmin,mmax);
4185 Form(
"Invariant mass of 2 local maxima cells vs (#it{E}_{1}+#it{E}_{2})/Ecluster, %s, %s",ptype[i].
Data(),sEBin[j].
Data()),
4186 120,0,1.2,mbins,mmin,mmax);
4192 Form(
"Invariant mass of N>2 local maxima cells vs (#it{E}_{1}+#it{E}_{2})/Ecluster, %s, %s",ptype[i].
Data(),sEBin[j].
Data()),
4193 120,0,1.2,mbins,mmin,mmax);
4201 Form(
"NLM vs E, %s, %s",ptype[i].
Data(),sEBin[j].
Data()),
4202 200,0,2,nMaxBins,0,nMaxBins);
4208 Form(
"NLM vs E, %s, %s, matched to a track",ptype[i].
Data(),sEBin[j].
Data()),
4209 200,0,2,nMaxBins,0,nMaxBins);
4215 Form(
"Invariant mass of 2 highest energy cells vs E, %s, %s",ptype[i].
Data(),sEBin[j].
Data()),
4216 200,0,2,mbins,mmin,mmax);
4222 Form(
"Invariant mass of 2 local maxima cells vs E, %s, %s",ptype[i].
Data(),sEBin[j].
Data()),
4223 200,0,2,mbins,mmin,mmax);
4229 Form(
"Invariant mass of N>2 local maxima cells vs E, %s, %s",ptype[i].
Data(),sEBin[j].
Data()),
4230 200,0,2,mbins,mmin,mmax);
4236 Form(
"#lambda_{0}^{2} vs E for N max = 1 %s, %s",ptype[i].
Data(),sEBin[j].
Data()),
4237 200,0,2,ssbins,ssmin,ssmax);
4243 Form(
"#lambda_{0}^{2} vs E for N max = 2 %s, %s",ptype[i].
Data(),sEBin[j].
Data()),
4244 200,0,2,ssbins,ssmin,ssmax);
4250 Form(
"#lambda_{0}^{2} vs E for N max > 2 %s, %s",ptype[i].
Data(),sEBin[j].
Data()),
4251 200,0,2,ssbins,ssmin,ssmax);
4265 "E vs Centrality, selected pi0 cluster with #it{NLM}=1",
4266 nptbins,ptmin,ptmax,100,0,100);
4272 "E vs Centrality, selected pi0 cluster with #it{NLM}=2",
4273 nptbins,ptmin,ptmax,100,0,100);
4279 "E vs Centrality, selected pi0 cluster with NLM>1",
4280 nptbins,ptmin,ptmax,100,0,100);
4288 "E vs Centrality, selected pi0 cluster with #it{NLM}=1",
4289 nptbins,ptmin,ptmax,100,0,100);
4295 "E vs Centrality, selected pi0 cluster with #it{NLM}=2",
4296 nptbins,ptmin,ptmax,100,0,100);
4302 "E vs Centrality, selected pi0 cluster with NLM>1",
4303 nptbins,ptmin,ptmax,100,0,100);
4311 fhEventPlanePi0NLocMax1 =
new TH2F(
"hEventPlanePi0NLocMax1",
"E vs Event Plane Angle, selected pi0 cluster with #it{NLM}=1",
4312 nptbins,ptmin,ptmax,100,0,TMath::Pi());
4317 fhEventPlanePi0NLocMax2 =
new TH2F(
"hEventPlanePi0NLocMax2",
"E vs Event Plane Angle, selected pi0 cluster with #it{NLM}=2",
4318 nptbins,ptmin,ptmax,100,0,TMath::Pi());
4323 fhEventPlanePi0NLocMaxN =
new TH2F(
"hEventPlanePi0NLocMaxN",
"E vs Event Plane Angle, selected pi0 cluster with NLM>1",
4324 nptbins,ptmin,ptmax,100,0,TMath::Pi());
4331 fhEventPlaneEtaNLocMax1 =
new TH2F(
"hEventPlaneEtaNLocMax1",
"E vs Event Plane Angle, selected pi0 cluster with #it{NLM}=1",
4332 nptbins,ptmin,ptmax,100,0,TMath::Pi());
4337 fhEventPlaneEtaNLocMax2 =
new TH2F(
"hEventPlaneEtaNLocMax2",
"E vs Event Plane Angle, selected pi0 cluster with #it{NLM}=2",
4338 nptbins,ptmin,ptmax,100,0,TMath::Pi());
4343 fhEventPlaneEtaNLocMaxN =
new TH2F(
"hEventPlaneEtaNLocMaxN",
"E vs Event Plane Angle, selected pi0 cluster with NLM>1",
4344 nptbins,ptmin,ptmax,100,0,TMath::Pi());
4353 for(Int_t i = 0; i < 4; i++)
4356 Form(
"Invariant mass of split clusters vs #lambda_{0}^{2}, #it{NLM}=1, %s",sEBin[i].
Data()),
4357 ssbins,ssmin,ssmax,mbins,mmin,mmax);
4363 Form(
"Invariant mass of split clusters vs #lambda_{0}^{2}, #it{NLM}=2, %s",sEBin[i].
Data()),
4364 ssbins,ssmin,ssmax,mbins,mmin,mmax);
4370 Form(
"Invariant mass of split clusters vs vs #lambda_{0}^{2}, NLM>2, %s",sEBin[i].
Data()),
4371 ssbins,ssmin,ssmax,mbins,mmin,mmax);
4378 Form(
"Invariant mass of split clusters vs split asymmetry, #it{NLM}=1, %s",sEBin[i].
Data()),
4379 200,-1,1,mbins,mmin,mmax);
4385 Form(
"Invariant mass of split clusters vs split asymmetry, #it{NLM}=2, %s",sEBin[i].
Data()),
4386 200,-1,1,mbins,mmin,mmax);
4392 Form(
"Invariant mass of split clusters vs split asymmetry, NLM>2, %s",sEBin[i].
Data()),
4393 200,-1,1,mbins,mmin,mmax);
4401 Form(
"Asymmetry of MC #pi^{0} vs #lambda_{0}^{2}, #it{NLM}=1, %s",sEBin[i].
Data()),
4402 ssbins,ssmin,ssmax,100,0,1);
4408 Form(
"Asymmetry of MC #pi^{0} vs #lambda_{0}^{2}, #it{NLM}=2, %s",sEBin[i].
Data()),
4409 ssbins,ssmin,ssmax,100,0,1);
4415 Form(
"Asymmetry of MC #pi^{0} vs #lambda_{0}^{2}, NLM>2, %s",sEBin[i].
Data()),
4416 ssbins,ssmin,ssmax,100,0,1);
4423 Form(
"Generated vs reconstructed asymmetry of split clusters from pi0, #it{NLM}=1, %s",sEBin[i].
Data()),
4430 Form(
"Generated vs reconstructed asymmetry of split clusters from pi0, #it{NLM}=2, %s",sEBin[i].
Data()),
4437 Form(
"Generated vs reconstructed asymmetry of split clusters from pi0, NLM>2, %s",sEBin[i].
Data()),
4447 Form(
"Invariant mass of 2 highest energy cells #sigma_{#eta #eta}^{2}, %s",sEBin[i].
Data()),
4448 ssbins,ssmin,ssmax,mbins,mmin,mmax);
4454 Form(
"Invariant mass of 2 local maxima cells #sigma_{#eta #eta}^{2}, %s",sEBin[i].
Data()),
4455 ssbins,ssmin,ssmax,mbins,mmin,mmax);
4461 Form(
"Invariant mass of N>2 local maxima cells vs #sigma_{#eta #eta}^{2}, %s",sEBin[i].
Data()),
4462 ssbins,ssmin,ssmax,mbins,mmin,mmax);