17 #include <TObjArray.h>
18 #include <TParticle.h>
19 #include <TDatabasePDG.h>
21 #include <TObjString.h>
27 #include "AliVCaloCells.h"
29 #include "AliVCluster.h"
30 #include "AliVTrack.h"
31 #include "AliVEvent.h"
32 #include "AliVEventHandler.h"
33 #include "AliAODMCParticle.h"
38 #include "AliPHOSGeoUtils.h"
39 #include "AliEMCALGeometry.h"
52 fFillAllCellTimeHisto(kTRUE),
53 fFillAllPosHisto(kFALSE), fFillAllPosHisto2(kFALSE),
55 fFillAllTMHisto(kTRUE),
56 fFillAllPi0Histo(kTRUE), fFillInvMassOpenAngle(kFALSE),
57 fFillPi0PairDiffTime(kFALSE), fFillInvMassInEMCALWithPHOSDCalAcc(kFALSE),
58 fFillEBinAcceptanceHisto(kFALSE),
59 fCorrelate(kTRUE), fStudyBadClusters(kFALSE),
60 fStudyClustersAsymmetry(kFALSE), fStudyExotic(kFALSE),
61 fStudyWeight(kFALSE), fStudyTCardCorrelation(kFALSE),
62 fStudyM02Dependence (kFALSE),
65 fNModules(12), fNRCU(2),
66 fNMaxCols(48), fNMaxRows(24),
67 fTimeCutMin(-10000), fTimeCutMax(10000),
69 fEMCALCellAmpMin(0), fPHOSCellAmpMin(0),
70 fEMCALClusterM02Min(0),
71 fEMCALClusterNCellMin(0), fPHOSClusterNCellMin(0),
75 fInvMassMinECut(0), fInvMassMaxECut(0),
76 fInvMassMinM02Cut(0), fInvMassMaxM02Cut(0),
77 fInvMassMaxOpenAngle(0), fInvMassMaxTimeDifference(0),
80 fExoNECrossCuts(0), fExoECrossCuts(),
81 fExoNDTimeCuts(0), fExoDTimeCuts(),
83 fClusterMomentum(), fClusterMomentum2(),
85 fConstantTimeShift(0),
89 fhEtaPhi(0), fhEtaPhiE(0),
90 fhECharged(0), fhPtCharged(0),
91 fhPhiCharged(0), fhEtaCharged(0),
92 fhEtaPhiCharged(0), fhEtaPhiECharged(0),
95 fhIM(0), fhIMSame(0), fhIMDiff(0),
96 fhIMDCAL(0), fhIMDCALSame(0), fhIMDCALDiff(0),
97 fhIMDCALPHOS(0), fhIMDCALPHOSSame(0),
98 fhIMEMCALPHOS(0), fhIMEMCALPHOSSame(0),
100 fhOpAngle(0), fhIMvsOpAngle(0),
101 fhNCellsPerCluster(0), fhNCellsPerClusterNoCut(0),
105 fhClusterTimeEnergy(0), fhCellTimeSpreadRespectToCellMax(0),
106 fhCellIdCellLargeTimeSpread(0), fhClusterPairDiffTimeE(0), fhClusterPairDiffTimeESameMod(0),
107 fhClusterMaxCellCloseCellRatio(0), fhClusterMaxCellCloseCellDiff(0),
108 fhClusterMaxCellDiff(0), fhClusterMaxCellDiffNoCut(0),
110 fhClusterMaxCellECross(0),
112 fhLambda0(0), fhLambda1(0), fhNLocMax(0),
114 fhEnergyTMEtaResidual1Cell(0), fhEnergyTMPhiResidual1Cell(0),
115 fhColRowExoticHighE1CellPosTime(0), fhColRowExoticHighE1CellNegTime(0), fhColRowExoticHighE1CellNulTime(0),
116 fhEnergyTMEtaResidualExotic(0), fhEnergyTMPhiResidualExotic(0),
117 fhColRowExoticHighEPosTime(0), fhColRowExoticHighENegTime(0), fhColRowExoticHighENulTime(0),
118 fhColRowHighEPosTime(0), fhColRowHighENegTime(0), fhColRowHighENulTime(0),
121 fhBadClusterEnergy(0), fhBadClusterTimeEnergy(0), fhBadClusterEtaPhi(0),
122 fhBadClusterPairDiffTimeE(0), fhBadCellTimeSpreadRespectToCellMax(0),
123 fhBadClusterMaxCellCloseCellRatio(0), fhBadClusterMaxCellCloseCellDiff(0), fhBadClusterMaxCellDiff(0),
124 fhBadClusterMaxCellDiffAverageTime(0), fhBadClusterMaxCellDiffWeightedTime(0),
125 fhBadClusterMaxCellECross(0),
126 fhBadClusterLambda0(0), fhBadClusterLambda1(0),
127 fhBadClusterDeltaIEtaDeltaIPhiE0(0), fhBadClusterDeltaIEtaDeltaIPhiE2(0),
128 fhBadClusterDeltaIEtaDeltaIPhiE6(0), fhBadClusterDeltaIA(0),
131 fhRNCells(0), fhXNCells(0),
132 fhYNCells(0), fhZNCells(0),
136 fhRCellE(0), fhXCellE(0),
137 fhYCellE(0), fhZCellE(0),
139 fhDeltaCellClusterRNCells(0), fhDeltaCellClusterXNCells(0),
140 fhDeltaCellClusterYNCells(0), fhDeltaCellClusterZNCells(0),
141 fhDeltaCellClusterRE(0), fhDeltaCellClusterXE(0),
142 fhDeltaCellClusterYE(0), fhDeltaCellClusterZE(0),
145 fhNCells(0), fhNCellsCutAmpMin(0),
146 fhAmplitude(0), fhAmpId(0),
147 fhEtaPhiAmpCell(0), fhEtaPhiCell(0),
149 fhTimeId(0), fhTimeL1UnCorrId(0), fhTimeAmp(0),
150 fhAmpIdLowGain(0), fhTimeIdLowGain(0), fhTimeAmpLowGain(0),
154 fhEMCALPHOSCorrNClusters(0), fhEMCALPHOSCorrEClusters(0),
155 fhEMCALPHOSCorrNCells(0), fhEMCALPHOSCorrECells(0),
156 fhEMCALDCALCorrNClusters(0), fhEMCALDCALCorrEClusters(0),
157 fhEMCALDCALCorrNCells(0), fhEMCALDCALCorrECells(0),
158 fhDCALPHOSCorrNClusters(0), fhDCALPHOSCorrEClusters(0),
159 fhDCALPHOSCorrNCells(0), fhDCALPHOSCorrECells(0),
160 fhCaloV0SCorrNClusters(0), fhCaloV0SCorrEClusters(0),
161 fhCaloV0SCorrNCells(0), fhCaloV0SCorrECells(0),
162 fhCaloV0MCorrNClusters(0), fhCaloV0MCorrEClusters(0),
163 fhCaloV0MCorrNCells(0), fhCaloV0MCorrECells(0),
164 fhCaloTrackMCorrNClusters(0), fhCaloTrackMCorrEClusters(0),
165 fhCaloTrackMCorrNCells(0), fhCaloTrackMCorrECells(0),
166 fhCaloCenNClusters(0), fhCaloCenEClusters(0),
167 fhCaloCenNCells(0), fhCaloCenECells(0),
168 fhCaloEvPNClusters(0), fhCaloEvPEClusters(0),
169 fhCaloEvPNCells(0), fhCaloEvPECells(0),
172 fhEMod(0), fhAmpMod(0),
173 fhEWeirdMod(0), fhAmpWeirdMod(0),
175 fhNClustersMod(0), fhNCellsMod(0),
176 fhNCellsSumAmpPerMod(0), fhNClustersSumEnergyPerMod(0),
177 fhNCellsPerClusterMod(0), fhNCellsPerClusterModNoCut(0),
178 fhNCellsPerClusterWeirdMod(0), fhNCellsPerClusterWeirdModNoCut(0),
180 fhGridCells(0), fhGridCellsE(0), fhGridCellsTime(0),
181 fhGridCellsLowGain(0), fhGridCellsELowGain(0), fhGridCellsTimeLowGain(0),
182 fhTimeAmpPerRCU(0), fhIMMod(0),
185 fhECellClusterRatio(0), fhECellClusterLogRatio(0),
186 fhEMaxCellClusterRatio(0), fhEMaxCellClusterLogRatio(0),
187 fhECellTotalRatio(0), fhECellTotalLogRatio(0),
188 fhECellTotalRatioMod(0), fhECellTotalLogRatioMod(0),
190 fhExoL0ECross(0), fhExoL1ECross(0),
193 fhRecoMCE(), fhRecoMCPhi(), fhRecoMCEta(),
194 fhRecoMCDeltaE(), fhRecoMCRatioE(),
195 fhRecoMCDeltaPhi(), fhRecoMCDeltaEta(),
198 fhGenMCE(), fhGenMCPt(), fhGenMCEtaPhi(),
199 fhGenMCAccE(), fhGenMCAccPt(), fhGenMCAccEtaPhi(),
202 fhEMVxyz(0), fhEMR(0),
203 fhHaVxyz(0), fhHaR(0),
204 fh1EOverP(0), fh2dR(0),
205 fh2EledEdx(0), fh2MatchdEdx(0),
206 fh1EOverPR02(0), fh1EleEOverP(0),
207 fhMCEle1EOverP(0), fhMCEle1dR(0), fhMCEle2MatchdEdx(0),
208 fhMCChHad1EOverP(0), fhMCChHad1dR(0), fhMCChHad2MatchdEdx(0),
209 fhMCNeutral1EOverP(0), fhMCNeutral1dR(0), fhMCNeutral2MatchdEdx(0),
210 fhMCEle1EOverPR02(0), fhMCChHad1EOverPR02(0), fhMCNeutral1EOverPR02(0),
211 fhMCEle1EleEOverP(0), fhMCChHad1EleEOverP(0), fhMCNeutral1EleEOverP(0),
212 fhTrackMatchedDEtaNeg(0), fhTrackMatchedDPhiNeg(0), fhTrackMatchedDEtaDPhiNeg(0),
213 fhTrackMatchedDEtaPos(0), fhTrackMatchedDPhiPos(0), fhTrackMatchedDEtaDPhiPos(0),
214 fhTrackMatchedDEtaNegMod(0), fhTrackMatchedDPhiNegMod(0),
215 fhTrackMatchedDEtaPosMod(0), fhTrackMatchedDPhiPosMod(0),
216 fhClusterTimeEnergyM02(0), fhCellTimeSpreadRespectToCellMaxM02(0),
217 fhClusterMaxCellCloseCellRatioM02(0), fhClusterMaxCellCloseCellDiffM02(0),
218 fhClusterMaxCellDiffM02(0), fhClusterMaxCellECrossM02(0), fhNCellsPerClusterM02(0)
221 for(
Int_t i =0; i < 12; i++)
223 for(
Int_t j = 0; j < 4; j++)
225 for(
Int_t k = 0; k < 3; k++)
232 for(
Int_t j = 0; j < 5; j++)
249 for (
Int_t ie = 0; ie < 10 ; ie++)
252 for (
Int_t idt = 0; idt < 5 ; idt++)
266 for(
Int_t i = 0; i < 7; i++)
277 for(
Int_t i = 0; i < 4; i++)
287 for(
Int_t i = 0; i < 14; i++)
297 for(
Int_t tm = 0; tm < 2; tm++)
366 for(
Int_t i = 0; i < 6; i++)
368 for(
Int_t j = 0; j < 6; j++)
441 for(
Int_t i = 0; i < 14; i++)
460 for(
Int_t j = 0; j < 6; j++)
462 for(
Int_t k = 0; k < 6; k++)
483 for(
Int_t i = 0; i < 12; i++)
540 if(phi < 0) phi += TMath::TwoPi();
551 for(
Int_t iclus2 = 0; iclus2 < caloClusters->GetEntriesFast(); iclus2++ )
553 AliVCluster* clus2 = (AliVCluster*) caloClusters->At(iclus2);
555 if(clus->GetID() == clus2->GetID())
continue;
560 if(
IsGoodCluster(absIdMax2, clus->GetM02(), clus->GetNCells(), cells) && clus2->GetM02() > 0.1 )
562 Double_t tof2 = clus2->GetTOF()*1.e9;
580 for (
Int_t ipos = 0; ipos < clus->GetNCells(); ipos++)
582 Int_t absId = clus->GetCellsAbsId()[ipos];
583 if(absId!=absIdMax && cells->GetCellAmplitude(absIdMax) > 0.01)
585 Float_t frac = cells->GetCellAmplitude(absId)/cells->GetCellAmplitude(absIdMax);
592 Double_t time = cells->GetCellTime(absId);
610 AliVCaloCells* cells,
618 for (
Int_t ipos = 0; ipos < clus->GetNCells(); ipos++)
620 Int_t id = clus->GetCellsAbsId()[ipos];
623 amp = cells->GetCellAmplitude(
id);
642 Int_t ncells = clus->GetNCells();
644 for (
Int_t ipos = 0; ipos < ncells; ipos++)
646 id = clus ->GetCellsAbsId()[ipos];
647 amp = cells->GetCellAmplitude(
id);
648 time = cells->GetCellTime(
id);
656 wTime += time*1e9 * w;
660 if(ncells > 0) aTime /= ncells;
663 if(wTot > 0) wTime /= wTot;
666 timeAverages[0] = aTime;
667 timeAverages[1] = wTime;
675 Int_t ncells = cells->GetNumberOfCells();
689 nCellsInModule[imod] = 0 ;
690 eCellsInModule[imod] = 0.;
693 Int_t icol = -1, icolAbs = -1;
694 Int_t irow = -1, irowAbs = -1;
704 for (
Int_t iCell = 0; iCell < cells->GetNumberOfCells(); iCell++)
706 if ( cells->GetCellNumber(iCell) < 0 )
continue;
708 AliDebug(2,Form(
"Cell : amp %f, absId %d", cells->GetAmplitude(iCell), cells->GetCellNumber(iCell)));
714 AliDebug(2,Form(
"\t module %d, column %d (%d), row %d (%d)", nModule,icolAbs,icol,irowAbs,irow));
716 if(nModule < fNModules)
723 if(
GetCaloUtils()->GetEMCALChannelStatus(nModule,icol,irow))
continue;
727 if(
GetCaloUtils()->GetPHOSChannelStatus(nModule,icol,irow) )
continue;
731 amp = cells->GetAmplitude(iCell)*recalF;
732 time = cells->GetTime(iCell);
733 id = cells->GetCellNumber(iCell);
734 highG = cells->GetCellHighGain(
id);
748 timeL1UnCorr *= 1.0e9;
753 AliDebug(1,Form(
"Remove cell with Time %f",time));
759 GetCaloUtils()->GetEMCALRecoUtils()->IsExoticCell(
id, cells, bc))
continue;
762 if(amp>=10)binWidthCorrection=1.0/4;
763 if(amp>=20)binWidthCorrection=1.0/10;
792 nCellsInModule[nModule]++ ;
844 if ( cellphi < 0 ) cellphi+=TMath::TwoPi();
858 Float_t rcell = TMath::Sqrt(cellpos[0]*cellpos[0]+cellpos[1]*cellpos[1]);
866 Int_t relId[4], module;
874 Float_t rcell = TMath::Sqrt(xyz.X()*xyz.X()+xyz.Y()*xyz.Y());
894 AliDebug(1,Form(
"Module %d, calo %s, N cells %d, sum Amp %f", imod,
GetCalorimeterString().
Data(), nCellsInModule[imod], eCellsInModule[imod]));
905 for (
Int_t iCell = 0; iCell < cells->GetNumberOfCells(); iCell++)
907 if ( cells->GetCellNumber(iCell) < 0 )
continue;
909 AliDebug(2,Form(
"Cell : amp %f, absId %d", cells->GetAmplitude(iCell), cells->GetCellNumber(iCell)));
913 AliDebug(2,Form(
"\t module %d, column %d, row %d", nModule,icol,irow));
915 if(nModule < fNModules)
922 if(
GetCaloUtils()->GetEMCALChannelStatus(nModule, icol, irow))
continue;
926 if(
GetCaloUtils()->GetPHOSChannelStatus (nModule, icol, irow) )
continue;
930 amp = cells->GetAmplitude(iCell)*recalF;
931 time = cells->GetTime(iCell);
932 id = cells->GetCellNumber(iCell);
939 GetReader()->GetInputEvent()->GetBunchCrossNumber());
947 AliDebug(1,Form(
"Remove cell with Time %f",time));
953 GetCaloUtils()->GetEMCALRecoUtils()->IsExoticCell(
id, cells, bc))
continue;
965 if(eCellsInModule[nModule] > 0)
967 Float_t ratioMod = amp/eCellsInModule[nModule];
976 delete [] nCellsInModule;
977 delete [] eCellsInModule;
985 Int_t nCaloCellsPerCluster = clus->GetNCells();
987 UShort_t * indexList = clus->GetCellsAbsId();
990 clus->GetPosition(pos);
995 for (
Int_t ipos = 0; ipos < nCaloCellsPerCluster; ipos++)
998 Int_t absId = indexList[ipos];
1015 Float_t r = TMath::Sqrt(pos[0] *pos[0] + pos[1] * pos[1] );
1016 Float_t rcell = TMath::Sqrt(cellpos[0]*cellpos[0] + cellpos[1]* cellpos[1]);
1024 Int_t relId[4], module;
1040 Float_t r = TMath::Sqrt(pos[0] * pos[0] + pos[1] * pos[1] );
1041 Float_t rcell = TMath::Sqrt(xyz.X() * xyz.X() + xyz.Y() * xyz.Y());
1065 Int_t icol = -1, irow = -1, iRCU = -1, icolAbs = -1, irowAbs = -1;
1069 Int_t ncells = clus->GetNCells();
1070 Double_t time = clus->GetTOF()*1.e9;
1072 Float_t deta = clus->GetTrackDz();
1073 Float_t dphi = clus->GetTrackDx();
1079 if(energy >= 5 && energy < 8)
1081 else if(energy >= 8)
1094 else if(exoticity > 0.97)
1096 if(energy >= 5 && energy < 8)
1098 else if(energy >= 8)
1111 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1113 Int_t absId = clus->GetCellsAbsId()[ipos];
1115 Int_t rowDiff = -100, colDiff = -100;
1119 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
1122 if ( exoticity > 0.97 )
1126 if(energy >= 5 && energy < 8)
1128 else if(energy >= 8)
1133 if(energy >= 5 && energy < 8)
1135 else if(energy >= 8)
1141 else if ( energy > 8 )
1151 if ( clus->GetDistanceToBadChannel() < 5 )
return ;
1154 Int_t etaRegion = -1, phiRegion = -1;
1157 if ( etaRegion !=0 ) return ;
1162 Int_t absIdList[ncells];
1169 Int_t nCellWithWeight = 1;
1177 Float_t eCellMax = cells->GetCellAmplitude(absIdMax);
1178 Double_t tCellMax = cells->GetCellTime(absIdMax);
1200 Float_t eCellSameRowSameTCardNearCol = 0.;
1201 Float_t eCellSameRowDiffTCardNearCol = 0.;
1202 Double_t tCellSameRowSameTCardNearCol = 0.;
1203 Double_t tCellSameRowDiffTCardNearCol = 0.;
1210 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1212 Int_t absId = clus->GetCellsAbsId()[ipos];
1214 Float_t eCell = cells->GetCellAmplitude(absId);
1215 Double_t tCell = cells->GetCellTime(absId);
1225 if( absId == absIdMax || weight < 0.01 )
continue;
1252 Int_t rowDiff = -100, colDiff = -100;
1255 Int_t indexType = -1;
1260 if(TMath::Abs(rowDiff) == 1) nearRow = kTRUE;
1261 if(TMath::Abs(colDiff) == 1) nearCol = kTRUE;
1263 if ( rowDiff == 0 && colDiff != 0 )
1265 if ( nearCol ) indexType = 6;
1271 else if ( rowDiff != 0 && colDiff == 0 )
1273 if ( nearRow ) indexType = 8;
1281 if ( nearRow && nearCol) indexType = 10;
1282 else indexType = 11;
1292 if ( rowDiff == 0 && colDiff != 0 )
1294 if ( nearCol ) indexType = 0;
1298 else if ( rowDiff != 0 && colDiff == 0 )
1300 if ( nearRow ) indexType = 2;
1305 if ( nearCol && nearRow ) indexType = 4;
1310 if ( rowDiff == 0 && TMath::Abs(colDiff) == 1 )
1314 eCellSameRowSameTCardNearCol = eCell;
1315 tCellSameRowSameTCardNearCol = tCell;
1319 eCellSameRowDiffTCardNearCol = eCell;
1320 tCellSameRowDiffTCardNearCol = tCell;
1326 Float_t eCellDiff = eCellMax - eCell;
1327 Float_t eClusDiff = energy - eCell;
1330 Float_t tCellDiff = tCellMax - tCell;
1338 if(exoticity > 0.97)
1348 if ( exoticity > 0.97 )
1351 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
1356 if(energy >= 5 && energy < 8)
1358 else if(energy >= 8)
1363 if(energy >= 5 && energy < 8)
1365 else if(energy >= 8)
1370 if ( eCell > emax2nd )
1417 Float_t ratioNcells = nCellWithWeight/(ncells*1.);
1422 if ( nCellWithWeight <= 1 )
return;
1430 if(energy >= 5 && energy < 8)
1432 else if(energy >= 8)
1435 if(exoticity > 0.97)
1437 if(energy >= 5 && energy < 8)
1439 else if(energy >= 8)
1483 if(eCellSameRowSameTCardNearCol > 0 && eCellSameRowDiffTCardNearCol > 0)
1485 Float_t eDiff = eCellSameRowSameTCardNearCol - eCellSameRowDiffTCardNearCol ;
1486 Float_t tDiff = tCellSameRowSameTCardNearCol - tCellSameRowDiffTCardNearCol ;
1493 if(exoticity > 0.97)
1503 Int_t nCorrInd = nCorr;
1504 if(nCorr > 4) nCorrInd = 5;
1506 Int_t nCorrNoInd = nCorrNo;
1507 if(nCorrNoInd > 4) nCorrNoInd = 5;
1554 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1556 Int_t absId = clus->GetCellsAbsId()[ipos];
1558 Float_t eCell = cells->GetCellAmplitude(absId);
1559 Double_t tCell = cells->GetCellTime(absId);
1569 if( absId == absIdMax || weight < 0.01 )
continue;
1571 Float_t tDiffMaxSecondary = tCellMax - tCell;
1574 if(exoticity > 0.97)
1607 Double_t tof2 = clus2->GetTOF()*1.e9;
1615 if( absIdMax == absIdMax2
1616 || !
IsGoodCluster(absIdMax2, clus2->GetM02(), clus2->GetNCells(), cells)
1648 if ( energy >= 5 && energy < 8)
1650 else if ( energy >= 8 )
1666 if ( exoticity > 0.97 )
1668 if ( energy >= 5 && energy < 8)
1670 else if ( energy >= 8 )
1733 if(exoticity > 0.97 )
1747 Float_t nCellRat = nCorr*1. / ((nCorr+nCorrNo)*1.);
1749 if(exoticity > 0.97 )
1838 if( clus->GetNCells() <= 3 )
return;
1846 for (
Int_t ipos = 0; ipos < clus->GetNCells(); ipos++)
1848 Int_t absId = clus->GetCellsAbsId()[ipos];
1853 if(dIphi < TMath::Abs(iphi-iphiMax)) dIphi = TMath::Abs(iphi-iphiMax);
1857 if(dIeta < TMath::Abs(ieta-ietaMax)) dIeta = TMath::Abs(ieta-ietaMax);
1861 Int_t ietaShift = ieta;
1862 Int_t ietaMaxShift = ietaMax;
1863 if (ieta > ietaMax) ietaMaxShift+=48;
1864 else ietaShift +=48;
1865 if(dIeta < TMath::Abs(ietaShift-ietaMaxShift)) dIeta = TMath::Abs(ietaShift-ietaMaxShift);
1869 Float_t dIA = 1.*(dIphi-dIeta)/(dIeta+dIphi);
1890 Int_t nLabel = clus->GetNLabels();
1891 Int_t* labels = clus->GetLabels();
1940 Double_t tof = clus->GetTOF()*1.e9;
1963 Int_t nCaloCellsPerCluster = clus->GetNCells();
1966 for(
Int_t iclus2 = 0; iclus2 < caloClusters->GetEntriesFast(); iclus2++ )
1968 AliVCluster* clus2 = (AliVCluster*) caloClusters->At(iclus2);
1970 if( clus->GetID() == clus2->GetID() )
continue;
1972 if( clus->GetM02() > 0.01 && clus2->GetM02() > 0.01 )
1976 Double_t tof2 = clus2->GetTOF()*1.e9;
1981 if ( nModule2 == nModule )
1986 if(nCaloCellsPerCluster > 1)
2000 for (
Int_t ipos = 0; ipos < nCaloCellsPerCluster; ipos++)
2002 Int_t absId = clus->GetCellsAbsId()[ipos];
2003 if( absId == absIdMax || cells->GetCellAmplitude(absIdMax) < 0.01 )
continue;
2005 Float_t frac = cells->GetCellAmplitude(absId)/cells->GetCellAmplitude(absIdMax);
2017 Double_t time = cells->GetCellTime(absId);
2037 AliDebug(1,Form(
"cluster: E %2.3f, pT %2.3f, eta %2.3f, phi %2.3f",e,pt,eta,phi*TMath::RadToDeg()));
2052 Int_t icol = -1, irow = -1, iRCU = -1, icolAbs = -1, irowAbs = -1;
2081 clus->GetPosition(pos);
2093 Float_t rxyz = TMath::Sqrt(pos[0]*pos[0]+pos[1]*pos[1]);
2116 AliVCaloCells* cells)
2119 Int_t *labels = 0x0;
2120 Int_t nCaloClusters = caloClusters->GetEntriesFast() ;
2121 Int_t nCaloClustersAccepted = 0 ;
2122 Int_t nCaloCellsPerCluster = 0 ;
2134 nClustersInModule[imod] = 0;
2135 energyInModule [imod] = 0;
2141 for(
Int_t iclus = 0; iclus < nCaloClusters; iclus++)
2143 AliDebug(1,Form(
"Cluster: %d/%d, data %d",iclus+1,nCaloClusters,
GetReader()->GetDataType()));
2145 AliVCluster* clus = (AliVCluster*) caloClusters->At(iclus);
2152 Double_t tof = clus->GetTOF()*1.e9;
2157 AliDebug(1,Form(
"Remove cluster with TOF %f",tof));
2170 nLabel = clus->GetNLabels();
2171 labels = clus->GetLabels();
2177 nCaloCellsPerCluster = clus->GetNCells();
2183 Double_t tmax = cells->GetCellTime(absIdMax);
2196 if(nModule >=0 && nModule < fNModules)
2205 Float_t ampMax = cells->GetCellAmplitude(absIdMax);
2211 Bool_t goodCluster =
IsGoodCluster(absIdMax, clus->GetM02(), nCaloCellsPerCluster, cells);
2214 if(ampMax > 0.01) eCrossFrac = 1-
GetECross(absIdMax,cells)/ampMax;
2216 AliDebug(1,Form(
"Accept cluster? %d",goodCluster));
2229 maxCellFraction, eCrossFrac, tmax);
2231 nCaloClustersAccepted++;
2240 nClustersInModule[nModule]++;
2242 energyInModule [nModule] += clus->E();
2252 Int_t mcOK = kFALSE;
2254 if(
IsDataMC() && nLabel > 0 && labels)
2266 && nCaloClusters > 1
2267 && nCaloCellsPerCluster > 1
2268 &&
GetCaloUtils()->GetNumberOfLocalMaxima(clus,cells) == 1
2284 AliDebug(1,Form(
"Module %d calo %s clusters %d, sum E %f", imod,
GetCalorimeterString().
Data(), nClustersInModule[imod], energyInModule[imod]));
2292 delete [] nClustersInModule;
2293 delete [] energyInModule;
2303 if(!labels || nLabels<=0)
2305 AliWarning(Form(
"Strange, labels array %p, n labels %d", labels,nLabels));
2309 AliDebug(1,Form(
"Primaries: nlabels %d",nLabels));
2314 if(phi < 0) phi +=TMath::TwoPi();
2316 AliAODMCParticle * aodprimary = 0x0;
2317 TParticle * primary = 0x0;
2320 Int_t label = labels[0];
2324 AliDebug(1,Form(
" *** bad label ***: label %d", label));
2343 AliDebug(1,Form(
"*** large label ***: label %d, n tracks %d", label,
GetMCStack()->GetNtrack()));
2349 pdg0 = TMath::Abs(primary->GetPdgCode());
2351 status = primary->GetStatusCode();
2352 vxMC = primary->Vx();
2353 vyMC = primary->Vy();
2354 iParent = primary->GetFirstMother();
2356 AliDebug(1,
"Cluster most contributing mother:");
2357 AliDebug(1,Form(
"\t Mother label %d, pdg %d, %s, status %d, parent %d",iMother, pdg0, primary->GetName(),status, iParent));
2365 pdg = TMath::Abs(primary->GetPdgCode());
2367 AliDebug(2,
"Converted cluster!. Find before conversion:");
2369 while((pdg == 22 || pdg == 11) && status != 1)
2371 Int_t iMotherOrg = iMother;
2374 status = primary->GetStatusCode();
2375 pdg = TMath::Abs(primary->GetPdgCode());
2376 iParent = primary->GetFirstMother();
2380 if(pdg==111 || pdg==221)
2382 primary =
GetMCStack()->Particle(iMotherOrg);
2392 AliDebug(2,Form(
"\t pdg %d, index %d, %s, status %d",pdg, iMother, primary->GetName(),status));
2395 AliDebug(1,
"Converted Cluster mother before conversion:");
2396 AliDebug(1,Form(
"\t Mother label %d, pdg %d, %s, status %d, parent %d",iMother, pdg, primary->GetName(), status, iParent));
2404 AliDebug(2,
"Overlapped Meson decay!, Find it:");
2406 while(pdg != 111 && pdg != 221)
2411 status = primary->GetStatusCode();
2412 pdg = TMath::Abs(primary->GetPdgCode());
2413 iParent = primary->GetFirstMother();
2415 if( iParent < 0 )
break;
2417 AliDebug(2,Form(
"\t pdg %d, %s, index %d",pdg, primary->GetName(),iMother));
2421 AliWarning(
"Tagged as Overlapped photon but meson not found, why?");
2426 AliDebug(2,Form(
"Overlapped %s decay, label %d",primary->GetName(),iMother));
2429 eMC = primary->Energy();
2431 phiMC = primary->Phi();
2432 etaMC = primary->Eta();
2433 pdg = TMath::Abs(primary->GetPdgCode());
2434 charge = (
Int_t) TDatabasePDG::Instance()->GetParticle(pdg)->Charge();
2437 else if(
GetReader()->ReadAODMCParticles() &&
2443 AliFatal(
"MCParticles not available!");
2447 pdg0 = TMath::Abs(aodprimary->GetPdgCode());
2449 status = aodprimary->IsPrimary();
2450 vxMC = aodprimary->Xv();
2451 vyMC = aodprimary->Yv();
2452 iParent = aodprimary->GetMother();
2454 AliDebug(1,
"Cluster most contributing mother:");
2455 AliDebug(1,Form(
"\t Mother label %d, pdg %d, Primary? %d, Physical Primary? %d, parent %d",
2456 iMother, pdg0, aodprimary->IsPrimary(), aodprimary->IsPhysicalPrimary(), iParent));
2461 AliDebug(2,
"Converted cluster!. Find before conversion:");
2465 pdg = TMath::Abs(aodprimary->GetPdgCode());
2467 while ((pdg == 22 || pdg == 11) && !aodprimary->IsPhysicalPrimary())
2469 Int_t iMotherOrg = iMother;
2472 status = aodprimary->IsPrimary();
2473 iParent = aodprimary->GetMother();
2474 pdg = TMath::Abs(aodprimary->GetPdgCode());
2478 if( pdg == 111 || pdg == 221 )
2490 AliDebug(2,Form(
"\t pdg %d, index %d, Primary? %d, Physical Primary? %d",
2491 pdg, iMother, aodprimary->IsPrimary(), aodprimary->IsPhysicalPrimary()));
2494 AliDebug(1,
"Converted Cluster mother before conversion:");
2495 AliDebug(1,Form(
"\t Mother label %d, pdg %d, parent %d, Primary? %d, Physical Primary? %d",
2496 iMother, pdg, iParent, aodprimary->IsPrimary(), aodprimary->IsPhysicalPrimary()));
2504 AliDebug(2,Form(
"Overlapped Meson decay!, Find it: PDG %d, mom %d",pdg, iMother));
2506 while(pdg != 111 && pdg != 221)
2510 status = aodprimary->IsPrimary();
2511 iParent = aodprimary->GetMother();
2512 pdg = TMath::Abs(aodprimary->GetPdgCode());
2514 if( iParent < 0 )
break;
2516 AliDebug(2,Form(
"\t pdg %d, index %d",pdg, iMother));
2520 AliWarning(
"Tagged as Overlapped photon but meson not found, why?");
2525 AliDebug(2,Form(
"Overlapped %s decay, label %d",aodprimary->GetName(),iMother));
2528 status = aodprimary->IsPrimary();
2529 eMC = aodprimary->E();
2531 phiMC = aodprimary->Phi();
2532 etaMC = aodprimary->Eta();
2533 pdg = TMath::Abs(aodprimary->GetPdgCode());
2534 charge = aodprimary->Charge();
2538 Float_t rVMC = TMath::Sqrt(vxMC*vxMC + vyMC*vyMC);
2539 if( ( pdg == 22 || TMath::Abs(pdg) == 11 ) && status != 1 )
2572 else if(charge == 0)
2587 if( mcIndex >= 0 && mcIndex < 7 && e > 0.5 && eMC > 0.5)
2598 if( primary || aodprimary )
return kTRUE ;
2614 if(phi < 0) phi +=TMath::TwoPi();
2632 Double_t dedx = track->GetTPCsignal();
2633 Int_t nITS = track->GetNcls(0);
2634 Int_t nTPC = track->GetNcls(1);
2635 Bool_t positive = kFALSE;
2636 if(track) positive = (track->Charge()>0);
2639 Float_t deta = clus->GetTrackDz();
2640 Float_t dphi = clus->GetTrackDx();
2641 Double_t dR = TMath::Sqrt(dphi*dphi + deta*deta);
2644 if( TMath::Abs(dphi) < 999 )
2681 if(dedx > 60 && dedx < 100)
2691 if(e > 0.5 && tmom > 0.5)
2697 if(dR < 0.02 && eOverP > 0.6 && eOverP < 1.2
2698 && clus->GetNCells() > 1 && nITS > 3 && nTPC > 20)
2706 Double_t charge = TDatabasePDG::Instance()->GetParticle(pdg)->Charge();
2708 if(TMath::Abs(pdg) == 11)
2732 else if(charge == 0)
2757 if(!caloClustersEMCAL || !caloClustersPHOS)
2759 AliDebug(1,Form(
"PHOS (%p) or EMCAL (%p) clusters array not available, do not correlate",caloClustersPHOS,caloClustersEMCAL));
2767 if(!cellsEMCAL || !cellsPHOS)
2769 AliDebug(1,Form(
"PHOS (%p) or EMCAL (%p) cells array ot available, do not correlate",cellsPHOS,cellsEMCAL));
2778 Float_t sumClusterEnergyEMCAL = 0;
2779 Float_t sumClusterEnergyDCAL = 0;
2780 Float_t sumClusterEnergyPHOS = 0;
2784 AliVCluster* cluster = 0;
2785 for(iclus = 0 ; iclus < caloClustersEMCAL->GetEntriesFast() ; iclus++)
2787 cluster = (AliVCluster*)caloClustersEMCAL->At(iclus);
2788 Float_t energy = cluster->E();
2790 if( energy < 0.5 )
continue;
2792 if(cluster->GetCellsAbsId()[0] < 12288)
2795 sumClusterEnergyEMCAL +=
energy;
2800 sumClusterEnergyDCAL +=
energy;
2804 for(iclus = 0 ; iclus < caloClustersPHOS ->GetEntriesFast(); iclus++)
2806 cluster = (AliVCluster*) caloClustersPHOS->At(iclus);
2808 energy = cluster->E();
2810 if( energy < 0.5 )
continue;
2813 sumClusterEnergyPHOS +=
energy;
2817 Int_t ncellsEMCAL = 0 ;
2818 Int_t ncellsDCAL = 0 ;
2819 Int_t ncellsPHOS = 0;
2821 Float_t sumCellEnergyEMCAL = 0;
2822 Float_t sumCellEnergyDCAL = 0;
2823 Float_t sumCellEnergyPHOS = 0;
2825 for(icell = 0 ; icell < cellsEMCAL->GetNumberOfCells() ; icell++)
2827 Float_t amp = cellsEMCAL->GetAmplitude(icell);
2828 Int_t cellId = cellsEMCAL->GetCellNumber(icell);
2835 sumCellEnergyEMCAL += amp;
2840 sumCellEnergyDCAL += amp;
2844 for(icell = 0 ; icell < cellsPHOS->GetNumberOfCells(); icell++)
2846 Float_t amp = cellsPHOS->GetAmplitude(icell);
2847 Int_t cellId = cellsPHOS->GetCellNumber(icell);
2849 if ( cellId < 0 ) continue ;
2854 sumCellEnergyPHOS += amp;
2880 Int_t ncl = nclPHOS;
2881 Float_t sumClusterEnergy = sumClusterEnergyPHOS;
2882 Int_t ncells = ncellsPHOS;
2883 Float_t sumCellEnergy = sumCellEnergyPHOS;
2887 ncl = nclEMCAL + nclDCAL;
2888 sumClusterEnergy = sumClusterEnergyEMCAL + sumClusterEnergyDCAL;
2889 ncells = ncellsEMCAL + ncellsDCAL;
2890 sumCellEnergy = sumCellEnergyEMCAL + sumCellEnergyDCAL;
2918 AliDebug(1,
"Correlate():");
2919 AliDebug(1,Form(
"\t EMCAL: N cells %d, N clusters %d, summed E cells %f, summed E clusters %f",
2920 ncellsEMCAL,nclEMCAL, sumCellEnergyEMCAL,sumClusterEnergyEMCAL));
2921 AliDebug(1,Form(
"\t DCAL : N cells %d, N clusters %d, summed E cells %f, summed E clusters %f",
2922 ncellsDCAL,nclDCAL, sumCellEnergyDCAL,sumClusterEnergyDCAL));
2923 AliDebug(1,Form(
"\t PHOS : N cells %d, N clusters %d, summed E cells %f, summed E clusters %f",
2924 ncellsPHOS,nclPHOS,sumCellEnergyPHOS,sumClusterEnergyPHOS));
2925 AliDebug(1,Form(
"\t V0 : Signal %d, Multiplicity %d, Track Multiplicity %d", v0S,v0M,trM));
2926 AliDebug(1,Form(
"\t centrality : %f, Event plane angle %f", cen,ep));
2935 const Int_t buffersize = 255;
2936 char onePar[buffersize] ;
2938 snprintf(onePar,buffersize,
"--- AliAnaCalorimeterQA ---:") ;
2965 return new TObjString(parList) ;
2972 AliVCluster *clus, AliVCaloCells* cells)
2976 AliDebug(1,Form(
"Low amplitude energy %f",ampMax));
2983 Int_t nc = clus->GetNCells();
3016 for (
Int_t icell = 0; icell < clus->GetNCells(); icell++)
3018 Int_t absId = clus->GetCellsAbsId()[icell];
3019 Double_t time = cells->GetCellTime(absId);
3045 outputContainer->SetName(
"QAHistos") ;
3104 fhE =
new TH1F (
"hE",
"#it{E} reconstructed clusters ", nptbins*5,ptmin,ptmax*5);
3105 fhE->SetXTitle(
"#it{E} (GeV)");
3106 outputContainer->Add(
fhE);
3108 fhPt =
new TH1F (
"hPt",
"#it{p}_{T} reconstructed clusters", nptbins,ptmin,ptmax);
3109 fhPt->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3110 outputContainer->Add(
fhPt);
3112 fhPhi =
new TH1F (
"hPhi",
"#varphi reconstructed clusters ",nphibins,phimin,phimax);
3113 fhPhi->SetXTitle(
"#varphi (rad)");
3114 outputContainer->Add(
fhPhi);
3116 fhEta =
new TH1F (
"hEta",
"#eta reconstructed clusters ",netabins,etamin,etamax);
3117 fhEta->SetXTitle(
"#eta ");
3118 outputContainer->Add(
fhEta);
3124 fhEtaPhiE =
new TH3F (
"hEtaPhiE",
"#eta vs #varphi vs energy, reconstructed clusters",
3125 netabins,etamin,etamax,nphibins,phimin,phimax,nptbins,ptmin,ptmax);
3133 fhEtaPhi =
new TH2F (
"hEtaPhi",
"#eta vs #varphi for #it{E} > 0.5 GeV, reconstructed clusters",
3134 netabins,etamin,etamax,nphibins,phimin,phimax);
3136 fhEtaPhi->SetYTitle(
"#varphi (rad)");
3142 nptbins,ptmin,ptmax, ntimebins,timemin,timemax);
3148 nptbins,ptmin,ptmax, tdbins,tdmin,tdmax);
3154 nptbins,ptmin,ptmax, tdbins,tdmin,tdmax);
3165 fhLambda0 =
new TH2F (
"hLambda0",
"shower shape, #lambda^{2}_{0} vs E",
3166 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3167 fhLambda0->SetXTitle(
"#it{E}_{cluster} (GeV)");
3168 fhLambda0->SetYTitle(
"#lambda^{2}_{0}");
3171 fhLambda1 =
new TH2F (
"hLambda1",
"shower shape, #lambda^{2}_{1} vs E",
3172 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3173 fhLambda1->SetXTitle(
"#it{E}_{cluster} (GeV)");
3174 fhLambda1->SetYTitle(
"#lambda^{2}_{1}");
3178 nptbins,ptmin,ptmax,10,0,10);
3179 fhNLocMax->SetXTitle(
"#it{E}_{cluster} (GeV)");
3186 (
"hColRowExoticHighE1CellPosTime",
3187 "column vs row, 1 cell, E > 8 GeV, t > 5 ns",
3188 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3194 (
"hColRowExoticHighEPosTime",
3195 "column vs row, exo > 0.97, E > 8 GeV, t > 5 ns",
3196 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3202 (
"hColRowHighEPosTime",
3203 "column vs row, exo < 0.97, E > 8 GeV, t > 5 ns",
3204 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3210 (
"hColRowExoticHighE1CellNegTime",
3211 "column vs row, 1 cell, E > 8 GeV, t < -5 ns",
3212 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3218 (
"hColRowExoticHighENegTime",
3219 "column vs row, exo > 0.97, E > 8 GeV, t < -5 ns",
3220 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3226 (
"hColRowHighENegTime",
3227 "column vs row, exo < 0.97, E > 8 GeV, t < -5 ns",
3228 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3234 (
"hColRowExoticHighE1CellNulTime",
3235 "column vs row, 1 cell, E > 8 GeV, -5 < t < 5 ns",
3236 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3242 (
"hColRowExoticHighENulTime",
3243 "column vs row, exo > 0.97, E > 8 GeV, -5 < t < 5 ns",
3244 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3250 (
"hColRowHighENulTime",
3251 "column vs row, exo < 0.97, E > 8 GeV, -5 < t < 5 ns",
3252 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3257 TString add[] = {
"",
"TrackMatched"};
3258 for(
Int_t tm = 0; tm < 2; tm++)
3261 (Form(
"hEnergyTime1Cell%s",add[tm].
Data()),
3262 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
3263 nptbins,ptmin,ptmax,300,-150,150);
3269 (Form(
"hEnergyTimeExotic%s",add[tm].
Data()),
3270 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
3271 nptbins,ptmin,ptmax,300,-150,150);
3277 (Form(
"hColRowExoticLowE1Cell%s",add[tm].
Data()),
3278 Form(
"column vs row, 1 cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3279 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3285 (Form(
"hColRowExoticHighE1Cell%s",add[tm].
Data()),
3286 Form(
"column vs row, 1 cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3287 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3293 (Form(
"hColRowExoticLowE%s",add[tm].
Data()),
3294 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3295 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3301 (Form(
"hColRowExoticHighE%s",add[tm].
Data()),
3302 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3303 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3309 (Form(
"hColRowExotic2ndCellDiffLowE%s",add[tm].
Data()),
3310 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3311 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3317 (Form(
"hColRowExotic2ndCellDiffHighE%s",add[tm].
Data()),
3318 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3319 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3325 (Form(
"hColRowExotic2ndCellSameLowE%s",add[tm].
Data()),
3326 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3327 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3333 (Form(
"hColRowExotic2ndCellSameHighE%s",add[tm].
Data()),
3334 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3335 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3341 (Form(
"hColRowTCardCorrNoSelectionExoticLowE%s",add[tm].
Data()),
3342 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3343 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3349 (Form(
"hColRowTCardCorrNoSelectionExoticHighE%s",add[tm].
Data()),
3350 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3351 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3357 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffLowE%s",add[tm].
Data()),
3358 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3359 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3365 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffHighE%s",add[tm].
Data()),
3366 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3367 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3373 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameLowE%s",add[tm].
Data()),
3374 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3375 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3381 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameHighE%s",add[tm].
Data()),
3382 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3383 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3390 (Form(
"hColRowTCardCorrNoSelectionLowE%s",add[tm].
Data()),
3391 Form(
"column vs row, max E cell for TCard correlation selected clusters, 5 < E < 8 GeV %s",add[tm].
Data()),
3392 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3398 (Form(
"hColRowTCardCorrNoSelectionHighE%s",add[tm].
Data()),
3399 Form(
"column vs row, max E cell for TCard correlation selected clusters, E > 8 GeV %s",add[tm].
Data()),
3400 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3408 (Form(
"hNCellsTCardCorrNoSelection%s",add[tm].
Data()),
3409 Form(
"# custer # cells vs #it{E} %s",add[tm].
Data()),
3410 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
3416 (Form(
"hNCellsTCardCorrWithWeightNoSelection%s",add[tm].
Data()),
3417 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
3418 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
3424 (Form(
"hNCellsTCardCorrRatioWithWeightNoSelection%s",add[tm].
Data()),
3425 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
3426 nptbins,ptmin,ptmax, 100,0,1);
3432 (Form(
"hLambda0TCardCorrNoSelection%s",add[tm].
Data()),
3433 Form(
"#lambda^{2}_{0} vs #it{E} %s",add[tm].
Data()),
3434 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3440 (Form(
"hLambda1TCardCorrNoSelection%s",add[tm].
Data()),
3441 Form(
"#lambda^{2}_{1} vs #it{E} %s",add[tm].
Data()),
3442 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3448 (Form(
"hLambda0NLM1TCardCorrNoSelection%s",add[tm].
Data()),
3449 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=1 %s",add[tm].
Data()),
3450 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3456 (Form(
"hLambda1NLM1TCardCorrNoSelection%s",add[tm].
Data()),
3457 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=1 %s",add[tm].
Data()),
3458 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3464 (Form(
"hLambda0NLM2TCardCorrNoSelection%s",add[tm].
Data()),
3465 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=2 %s",add[tm].
Data()),
3466 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3472 (Form(
"hLambda1NLM2TCardCorrNoSelection%s",add[tm].
Data()),
3473 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=2 %s",add[tm].
Data()),
3474 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3481 (Form(
"hLambdaRTCardCorrNoSelection%s",add[tm].
Data()),
3482 Form(
"#lambda^{1}_{0}/#lambda^{2}_{0} vs #it{E} %s",add[tm].
Data()),
3483 nptbins,ptmin,ptmax,110,0,1.1);
3489 (Form(
"hNLocMaxTCardCorrNoSelection%s",add[tm].
Data()),
3490 Form(
"#it{n}_{LM} vs E %s",add[tm].
Data()),
3491 nptbins,ptmin,ptmax,10,0,10);
3497 (Form(
"hEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3498 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
3499 nptbins,ptmin,ptmax,100,0,1);
3505 (Form(
"hEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3506 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3507 nptbins,ptmin,ptmax,100,0,1);
3513 (Form(
"hEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3514 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3515 nptbins,ptmin,ptmax,100,0,1);
3521 (Form(
"hE2ndRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3522 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
3523 nptbins,ptmin,ptmax,100,0,0.5);
3529 (Form(
"hE2ndRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3530 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3531 nptbins,ptmin,ptmax,100,0,0.5);
3537 (Form(
"hE2ndRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3538 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3539 nptbins,ptmin,ptmax,100,0,0.5);
3546 (Form(
"hE2ndEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3547 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
3548 nptbins,ptmin,ptmax,100,0,1);
3554 (Form(
"hE2ndEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3555 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3556 nptbins,ptmin,ptmax,100,0,1);
3562 (Form(
"hE2ndEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3563 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3564 nptbins,ptmin,ptmax,100,0,1);
3571 (Form(
"hECellClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3572 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
3573 nptbins,ptmin,ptmax,100,0,0.5);
3579 (Form(
"hECellClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3580 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3581 nptbins,ptmin,ptmax,100,0,0.5);
3587 (Form(
"hECellClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3588 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3589 nptbins,ptmin,ptmax,100,0,0.5);
3595 (Form(
"hLogECellNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3596 Form(
"log(#it{E}_{cell}) vs E_{cluster}, #it{n}_{LM}=1, w > 0.01 %s",add[tm].
Data()),
3597 nptbins,ptmin,ptmax,150,-3,3);
3603 (Form(
"hLogECellNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3604 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}=2, w > 0.01 %s",add[tm].
Data()),
3605 nptbins,ptmin,ptmax,150,-3,3);
3611 (Form(
"hLogECellNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3612 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}>2, w > 0.01 %s",add[tm].
Data()),
3613 nptbins,ptmin,ptmax,150,-3,3);
3620 (Form(
"hECellWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3621 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
3622 nptbins,ptmin,ptmax,90,0,4.5);
3628 (Form(
"hECellWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3629 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3630 nptbins,ptmin,ptmax,90,0,4.5);
3636 (Form(
"hECellWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3637 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster}))vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3638 nptbins,ptmin,ptmax,90,0,4.5);
3644 (Form(
"hExoticTCardCorrNoSelection%s",add[tm].
Data()),
3645 Form(
"exoticity vs #it{E} %s",add[tm].
Data()),
3646 nptbins,ptmin,ptmax,200,-1,1);
3652 (Form(
"hNCellsTCardSameAndDiff%s",add[tm].
Data()),
3653 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, %s",add[tm].
Data()),
3654 nptbins,ptmin,ptmax,100,0,1);
3660 (Form(
"hNCellsTCardSameAndDiffFraction_Exotic%s",add[tm].
Data()),
3661 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, exo > 0.97 %s",add[tm].
Data()),
3662 nptbins,ptmin,ptmax,100,0,1);
3669 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterE%s",add[tm].
Data()),
3670 Form(
"#Delta row = 0, |#Delta col = 1|, with respect to leading cell, #it{E}_{cell}^{same TCard}-#it{E}_{cell}^{diff TCard} vs #it{E}_{cluster} %s",add[tm].
Data()),
3671 nptbins,ptmin,ptmax,200,-10,10);
3677 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterE%s",add[tm].
Data()),
3678 Form(
"#Delta row = 0, |#Delta col = 1|, with respect to leading cell, #it{t}_{cell}^{same TCard}-#it{t}_{cell}^{diff TCard} vs #it{E}_{cluster} %s",add[tm].
Data()),
3679 nptbins,ptmin,ptmax,200,-100,100);
3685 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxE%s",add[tm].
Data()),
3686 Form(
"#Delta row = 0, |#Delta col = 1|, with respect to leading cell, #it{E}_{cell}^{same TCard}-#it{E}_{cell}^{diff TCard} vs #it{E}_{cell max} %s",add[tm].
Data()),
3687 nptbins,ptmin,ptmax,200,-10,10);
3693 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxE%s",add[tm].
Data()),
3694 Form(
"#Delta row = 0, |#Delta col = 1|, with respect to leading cell, #it{t}_{cell}^{same TCard}-#it{t}_{cell}^{diff TCard} vs #it{E}_{cell max} %s",add[tm].
Data()),
3695 nptbins,ptmin,ptmax,200,-100,100);
3701 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterEExo%s",add[tm].
Data()),
3702 Form(
"#Delta row = 0, |#Delta col = 1|, with respect to leading cell, #it{E}_{cell}^{same TCard}-#it{E}_{cell}^{diff TCard} vs #it{E}_{cluster}, exo > 0.97 %s",add[tm].
Data()),
3703 nptbins,ptmin,ptmax,200,-10,10);
3709 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterEExo%s",add[tm].
Data()),
3710 Form(
"#Delta row = 0, |#Delta col = 1|, with respect to leading cell, #it{t}_{cell}^{same TCard}-#it{t}_{cell}^{diff TCard} vs #it{E}_{cluster}, exo > 0.97 %s",add[tm].
Data()),
3711 nptbins,ptmin,ptmax,200,-100,100);
3717 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxEExo%s",add[tm].
Data()),
3718 Form(
"#Delta row = 0, |#Delta col = 1|, with respect to leading cell, #it{E}_{cell}^{same TCard}-#it{E}_{cell}^{diff TCard} vs #it{E}_{cell max}, exo > 0.97 %s",add[tm].
Data()),
3719 nptbins,ptmin,ptmax,200,-10,10);
3725 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxEExo%s",add[tm].
Data()),
3726 Form(
"#Delta row = 0, |#Delta col = 1|, with respect to leading cell, #it{t}_{cell}^{same TCard}-#it{t}_{cell}^{diff TCard} vs #it{E}_{cell max}, exo > 0.97 %s",add[tm].
Data()),
3727 nptbins,ptmin,ptmax,200,-100,100);
3732 for(
Int_t i = 0; i < 6; i++)
3734 for(
Int_t j = 0; j < 6; j++)
3737 (Form(
"hLambda0TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3738 Form(
"#lambda^{2}_{0} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3739 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3745 (Form(
"hLambda1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3746 Form(
"#lambda^{2}_{1} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3747 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3753 (Form(
"hLambda0NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3754 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=1, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3755 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3761 (Form(
"hLambda1NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3762 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=1, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3763 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3769 (Form(
"hLambda0NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3770 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=2, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3771 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3777 (Form(
"hLambda1NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3778 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=2, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3779 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3794 (Form(
"hNLocMaxTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3795 Form(
"#it{n}_{LM} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3796 nptbins,ptmin,ptmax,10,0,10);
3802 (Form(
"hEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3803 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs #it{E}_{cluster}, #it{n}_{LM} = 1, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3804 nptbins,ptmin,ptmax,100,0,1);
3810 (Form(
"hEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3811 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs #it{E}_{cluster}, #it{n}_{LM} = 2, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3812 nptbins,ptmin,ptmax,100,0,1);
3818 (Form(
"hEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3819 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs #it{E}_{cluster}, #it{n}_{LM} > 2, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3820 nptbins,ptmin,ptmax,100,0,1);
3826 (Form(
"hE2ndRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3827 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cluster} vs #it{E}_{cluster}, #it{n}_{LM} = 1, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3828 nptbins,ptmin,ptmax,100,0,0.5);
3834 (Form(
"hE2ndRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3835 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs #it{E}_{cluster}, #it{n}_{LM} = 2, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3836 nptbins,ptmin,ptmax,100,0,0.5);
3842 (Form(
"hE2ndRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3843 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs #it{E}_{cluster}, #it{n}_{LM} > 2, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3844 nptbins,ptmin,ptmax,100,0,0.5);
3850 (Form(
"hE2ndEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3851 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cell}^{max} vs #it{E}_{cluster}, #it{n}_{LM} = 1, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3852 nptbins,ptmin,ptmax,100,0,1);
3858 (Form(
"hE2ndEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3859 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs #it{E}_{cluster}, #it{n}_{LM} = 2, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3860 nptbins,ptmin,ptmax,100,0,1);
3866 (Form(
"hE2ndEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3867 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs #it{E}_{cell}^{max}, #it{n}_{LM} > 2, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3868 nptbins,ptmin,ptmax,100,0,1);
3875 (Form(
"hECellClusRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3876 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs #it{E}_{cluster}, #it{n}_{LM} = 1, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3877 nptbins,ptmin,ptmax,100,0,0.5);
3883 (Form(
"hECellClusRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3884 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs #it{E}_{cluster}, #it{n}_{LM} = 2, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3885 nptbins,ptmin,ptmax,100,0,0.5);
3891 (Form(
"hECellClusRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3892 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs #it{E}_{cluster}, #it{n}_{LM} > 2, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3893 nptbins,ptmin,ptmax,100,0,0.5);
3899 (Form(
"hLogECellNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3900 Form(
"log(#it{E}_{cell}) vs #it{E}_{cluster}, #it{n}_{LM} = 1, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3901 nptbins,ptmin,ptmax,150,-3,3);
3907 (Form(
"hLogECellNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3908 Form(
"log(#it{E}_{cell}) vs #it{E}_{cluster}, #it{n}_{LM} = 2, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3909 nptbins,ptmin,ptmax,150,-3,3);
3915 (Form(
"hLogECellNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3916 Form(
"log(#it{E}_{cell}) vs #it{E}_{cluster}, #it{n}_{LM} > 2, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3917 nptbins,ptmin,ptmax,150,-3,3);
3923 (Form(
"hECellWeightNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3924 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs #it{E}_{cluster}, #it{n}_{LM} = 1, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3925 nptbins,ptmin,ptmax,90,0,4.5);
3931 (Form(
"hECellWeightNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3932 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs #it{E}_{cluster}, #it{n}_{LM} = 2, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3933 nptbins,ptmin,ptmax,90,0,4.5);
3939 (Form(
"hECellWeightNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3940 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs #it{E}_{cluster}, #it{n}_{LM} > 2, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3941 nptbins,ptmin,ptmax,90,0,4.5);
3947 (Form(
"hMassEClusTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3948 Form(
"#it{M}_{#gamma #gamma} vs #it{E}_{cluster}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3949 nptbins,ptmin,ptmax,nmassbins,massmin,massmax);
3963 (Form(
"hExoticTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3964 Form(
"exoticity vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3965 nptbins,ptmin,ptmax,200,-1,1);
3971 (Form(
"hTimeDiffTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3972 Form(
"#it{t}_{cell}^{max}-#it{t}_{cell}^{other} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3973 nptbins,ptmin,ptmax,300,-150,150);
3979 (Form(
"hTimeDiffExoTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3980 Form(
"#it{t}_{cell}^{max}-#it{t}_{cell}^{other} vs #it{E}, N cells with w > 0.01, exoticity > 0.97, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3981 nptbins,ptmin,ptmax,300,-150,150);
3987 (Form(
"hColRowTCardCorrelNCellLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3988 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3989 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3995 (Form(
"hColRowTCardCorrelNCellHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3996 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3997 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
4003 (Form(
"hColRowTCardCorrelNCellExoticLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4004 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4005 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
4011 (Form(
"hColRowTCardCorrelNCellExoticHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4012 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4013 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
4224 (Form(
"hLambda0Exoticity_EBin%d%s",i,add[tm].
Data()),
4225 Form(
"#lambda^{2}_{0} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
4226 200,-1,1,ssbins,ssmin,ssmax);
4232 (Form(
"hLambda1Exoticity_EBin%d%s",i,add[tm].
Data()),
4233 Form(
"#lambda^{2}_{1} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
4234 200,-1,1,ssbins,ssmin,ssmax);
4248 (Form(
"hNCellsExoticity_EBin%d%s",i,add[tm].
Data()),
4249 Form(
"#it{n}_{cells} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
4250 200,-1,1,nceclbins,nceclmin,nceclmax);
4251 fhNCellsExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
4256 (Form(
"hTimeExoticity_EBin%d%s",i,add[tm].
Data()),
4258 200,-1,1,100,-25,25);
4259 fhTimeExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
4264 (Form(
"hLambda0Lambda1_EBin%d%s",i,add[tm].
Data()),
4265 Form(
"#lambda^{2}_{0} vs #lambda^{2}_{1}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
4266 ssbins,ssmin,ssmax,ssbins,ssmin,ssmax);
4312 (Form(
"hNCellsTCardSameAndDiff_EBin%d%s",i,add[tm].
Data()),
4313 Form(
"#it{n}_{cells} same TCard vs diff TCard, w > 0.01, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
4314 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
4320 (Form(
"hNCellsTCardSameAndDiff_Exotic_EBin%d%s",i,add[tm].
Data()),
4321 Form(
"#it{n}_{cells} same TCard vs diff TCard, w > 0.01, %2.2f<#it{E}<%2.2f GeV, exo > 0.97 %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
4322 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
4328 for(
Int_t j = 0; j < 6; j++)
4330 for(
Int_t k = 0; k < 6; k++)
4333 (Form(
"hLambda0Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
4334 Form(
"#lambda^{2}_{0} vs #it{exoticity}, #it{n}_{cell} TCard same = %d, diff =%d, #it{E}>8 GeV %s",j,k,add[tm].
Data()),
4335 200,-1,1,ssbins,ssmin,ssmax);
4341 (Form(
"hLambda1Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
4342 Form(
"#lambda^{2}_{1} vs #it{exoticity}, #it{n}_{cell} TCard same = %d, diff =%d, #it{E}>8 GeV %s",j,k,add[tm].
Data()),
4343 200,-1,1,ssbins,ssmin,ssmax);
4440 for(
Int_t i = 0; i < 12; i++)
4443 (Form(
"hTCardCorrECellMaxDiff_Case%d%s",i,add[tm].
Data()),
4444 Form(
"#it{E}_{cell}^{max}-#it{E}_{cell} vs #it{E}_{cluster}, for (un)correlated cells in TCard, case %d %s",i,add[tm].
Data()),
4445 nptbins,ptmin,ptmax,210,-1,20);
4451 (Form(
"hTCardCorrEClusterDiff_Case%d%s",i,add[tm].
Data()),
4452 Form(
"#it{E}_{cluster}-#it{E}_{cell} vs #it{E}_{cluster}, for (un)correlated cells in TCard, case %d %s",i,add[tm].
Data()),
4453 nptbins,ptmin,ptmax,210,-1,20);
4475 (Form(
"hTCardCorrTCellMaxDiff_Case%d%s",i,add[tm].
Data()),
4476 Form(
"#it{t}_{cell}^{max}-#it{t}_{cell} vs #it{E}_{cluster}, for (un)correlated cells in TCard, case %d %s",i,add[tm].
Data()),
4477 nptbins,ptmin,ptmax,1000,-100,100);
4483 (Form(
"hTCardCorrECellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
4484 Form(
"#it{E}_{cell}^{max}-#it{E}_{cell} vs #it{E}_{cluster}, for (un)correlated cells in TCard, exoticity > 0.97, case %d %s",i,add[tm].
Data()),
4485 nptbins,ptmin,ptmax,210,-1,20);
4491 (Form(
"hTCardCorrEClusterDiffExo_Case%d%s",i,add[tm].
Data()),
4492 Form(
"#it{E}_{cluster}-#it{E}_{cell} vs #it{E}_{cluster}, for (un)correlated cells in TCard, exoticity > 0.97, case %d %s",i,add[tm].
Data()),
4493 nptbins,ptmin,ptmax,210,-1,20);
4515 (Form(
"hTCardCorrTCellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
4516 Form(
"#it{t}_{cell}^{max}-#it{t}_{cell} vs #it{E}_{cluster}, for (un)correlated cells in TCard, exoticity > 0.97, case %d %s",i,add[tm].
Data()),
4517 nptbins,ptmin,ptmax,1000,-100,100);
4525 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
4531 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
4537 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
4543 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
4551 (Form(
"hTMPhiResidual_EBin%d",i),
4552 Form(
"#Delta #varphi_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
4553 200,-1,1,nresphibins,resphimin,resphimax);
4559 (Form(
"hTMEtaResidual_EBin%d",i),
4560 Form(
"#Delta #eta_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
4561 200,-1,1,nresetabins,resetamin,resetamax);
4586 nptbins,ptmin,ptmax, 100,0,1.);
4592 nptbins,ptmin,ptmax, 500,0,100.);
4597 fhClusterMaxCellDiff =
new TH2F (
"hClusterMaxCellDiff",
"energy vs difference of cluster energy - max cell energy / cluster energy, good clusters",
4598 nptbins,ptmin,ptmax, 500,0,1.);
4605 fhClusterMaxCellDiffNoCut =
new TH2F (
"hClusterMaxCellDiffNoCut",
"energy vs difference of cluster energy - max cell energy / cluster energy",
4606 nptbins,ptmin,ptmax, 500,0,1.);
4612 fhClusterMaxCellECross =
new TH2F (
"hClusterMaxCellECross",
"1 - Energy in cross around max energy cell / max energy cell vs cluster energy, good clusters",
4613 nptbins,ptmin,ptmax, 400,-1,1.);
4621 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
4627 fhNCellsPerCluster =
new TH2F (
"hNCellsPerCluster",
"# cells per cluster vs energy",nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
4647 fhNClusters =
new TH1F (
"hNClusters",
"# clusters", nclbins,nclmin,nclmax);
4653 fhBadClusterEnergy =
new TH1F (
"hBadClusterEnergy",
"Bad cluster energy", nptbins,ptmin,ptmax);
4658 netabins,etamin,etamax,nphibins,phimin,phimax);
4664 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4669 fhBadClusterLambda1 =
new TH2F (
"hBadClusterLambda1",
"Bad cluster,shower shape, #lambda^{2}_{1} vs E for bad cluster ",
4670 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4676 nptbins,ptmin,ptmax, 100,0,1.);
4682 nptbins,ptmin,ptmax, 500,0,100);
4687 fhBadClusterMaxCellDiff =
new TH2F (
"hBadClusterMaxCellDiff",
"energy vs difference of cluster energy - max cell energy / cluster energy for bad clusters",
4688 nptbins,ptmin,ptmax, 500,0,1.);
4694 nptbins,ptmin,ptmax, ntimebins,timemin,timemax);
4699 fhBadClusterPairDiffTimeE =
new TH2F(
"hBadClusterPairDiffTimeE",
"cluster pair time difference (bad - good) vs E from bad cluster",nptbins,ptmin,ptmax, tdbins,tdmin,tdmax);
4704 fhBadClusterMaxCellECross =
new TH2F (
"hBadClusterMaxCellECross",
"1 - #it{E}_{+} around max energy cell / max energy cell vs cluster energy, bad clusters",
4705 nptbins,ptmin,ptmax, 400,-1,1.);
4734 "#lambda^{2}_{0} vs 1-#it{E}_{+}/#it{E}_{max} for E > 5 GeV",
4735 400,0,1,ssbins,ssmin,ssmax);
4736 fhExoL0ECross ->SetXTitle(
"1-#it{E}_{+}/#it{E}_{cell max}");
4741 "#lambda^{2}_{1} vs 1-#it{E}_{+}/#it{E}_{max} for E > 5 GeV",
4742 400,0,1,ssbins,ssmin,ssmax);
4743 fhExoL1ECross ->SetXTitle(
"1-#it{E}_{+}/#it{E}_{cell max}");
4751 Form(
"#Delta time = t_{max}-t_{cells} vs #it{E}_{cluster} for exotic, 1-#it{E}_{+}/#it{E}_{max} < %2.2f",
fExoECrossCuts[ie]),
4752 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
4753 fhExoDTime[ie] ->SetYTitle(
"#Delta #it{t} (ns)");
4759 fhExoNCell[ie][idt] =
new TH2F(Form(
"hExoNCell_ECross%d_DT%d",ie,idt),
4760 Form(
"N cells per cluster vs E cluster, 1-#it{E}_{+}/#it{E}_{max} < %2.2f, #Delta t < %2.0f",
fExoECrossCuts[ie],
fExoDTimeCuts[idt]),
4761 nptbins,ptmin,ptmax,nceclbins,nceclmin,nceclmax);
4762 fhExoNCell[ie][idt] ->SetYTitle(
"#it{n}_cells");
4763 fhExoNCell[ie][idt] ->SetXTitle(
"#it{E} (GeV)");
4766 fhExoL0 [ie][idt] =
new TH2F(Form(
"hExoL0_ECross%d_DT%d",ie,idt),
4767 Form(
"#lambda^{2}_{0} vs E cluster for exotic, 1-#it{E}_{+}/#it{E}_{max} < %2.2f, #Delta t = %2.0f",
fExoECrossCuts[ie],
fExoDTimeCuts[idt]),
4768 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4769 fhExoL0 [ie][idt] ->SetYTitle(
"#lambda^{2}_{0}");
4770 fhExoL0 [ie][idt] ->SetXTitle(
"#it{E} (GeV)");
4771 outputContainer->Add(
fhExoL0[ie][idt]) ;
4773 fhExoL1 [ie][idt] =
new TH2F(Form(
"hExoL1_ECross%d_DT%d",ie,idt),
4774 Form(
"#lambda^{2}_{1} vs E cluster for exotic, 1-#it{E}_{+}/#it{E}_{max} < %2.2f, #Delta t = %2.0f",
fExoECrossCuts[ie],
fExoDTimeCuts[idt]),
4775 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4776 fhExoL1 [ie][idt] ->SetYTitle(
"#lambda^{2}_{1}");
4777 fhExoL1 [ie][idt] ->SetXTitle(
"#it{E} (GeV)");
4778 outputContainer->Add(
fhExoL1[ie][idt]) ;
4780 fhExoECross[ie][idt] =
new TH2F(Form(
"hExoECross_ECross%d_DT%d",ie,idt),
4781 Form(
"#it{E} cross for cells vs E cell, 1-#it{E}_{+}/#it{E}_{max} < %2.2f, #Delta t < %2.0f",
fExoECrossCuts[ie],
fExoDTimeCuts[idt]),
4782 nptbins,ptmin,ptmax,400,0,1);
4783 fhExoECross[ie][idt] ->SetYTitle(
"1-#it{E}_{+}/#it{E}_{cell max}");
4784 fhExoECross[ie][idt] ->SetXTitle(
"#it{E}_{cell} (GeV)");
4787 fhExoTime [ie][idt] =
new TH2F(Form(
"hExoTime_ECross%d_DT%d",ie,idt),
4788 Form(
"Time of cluster (max cell) vs E cluster for exotic, 1-#it{E}_{+}/#it{E}_{max} < %2.2f, #Delta t = %2.0f",
fExoECrossCuts[ie],
fExoDTimeCuts[idt]),
4789 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
4790 fhExoTime [ie][idt] ->SetYTitle(
"#it{t}_{max} (ns)");
4791 fhExoTime [ie][idt] ->SetXTitle(
"#it{E} (GeV)");
4792 outputContainer->Add(
fhExoTime[ie][idt]) ;
4795 Form(
"#lambda^{2}_{0} vs N cells per clusters for E > 5 GeV, 1-#it{E}_{+}/#it{E}_{max} < %2.2f, #Delta t = %2.0f",
fExoECrossCuts[ie],
fExoDTimeCuts[idt]),
4796 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
4802 Form(
"#lambda^{2}_{1} vs N cells per clusters for E > 5 GeV, 1-#it{E}_{+}/#it{E}_{max} < %2.2f, #Delta t = %2.0f",
fExoECrossCuts[ie],
fExoDTimeCuts[idt]),
4803 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
4815 fhDeltaIEtaDeltaIPhiE0[0] =
new TH2F (
"hDeltaIEtaDeltaIPhiE0",
" Cluster size in columns vs rows for E < 2 GeV, #it{n}_{cells} > 3",
4821 fhDeltaIEtaDeltaIPhiE2[0] =
new TH2F (
"hDeltaIEtaDeltaIPhiE2",
" Cluster size in columns vs rows for 2 <E < 6 GeV, #it{n}_{cells} > 3",
4827 fhDeltaIEtaDeltaIPhiE6[0] =
new TH2F (
"hDeltaIEtaDeltaIPhiE6",
" Cluster size in columns vs rows for E > 6 GeV, #it{n}_{cells} > 3",
4833 fhDeltaIA[0] =
new TH2F (
"hDeltaIA",
" Cluster *asymmetry* in cell units vs E",
4834 nptbins,ptmin,ptmax,21,-1.05,1.05);
4835 fhDeltaIA[0]->SetXTitle(
"#it{E}_{cluster}");
4836 fhDeltaIA[0]->SetYTitle(
"#it{A}_{cell in cluster}");
4839 fhDeltaIAL0[0] =
new TH2F (
"hDeltaIAL0",
" Cluster *asymmetry* in cell units vs #lambda^{2}_{0}",
4840 ssbins,ssmin,ssmax,21,-1.05,1.05);
4842 fhDeltaIAL0[0]->SetYTitle(
"#it{A}_{cell in cluster}");
4845 fhDeltaIAL1[0] =
new TH2F (
"hDeltaIAL1",
" Cluster *asymmetry* in cell units vs #lambda^{2}_{1}",
4846 ssbins,ssmin,ssmax,21,-1.05,1.05);
4848 fhDeltaIAL1[0]->SetYTitle(
"#it{A}_{cell in cluster}");
4851 fhDeltaIANCells[0] =
new TH2F (
"hDeltaIANCells",
" Cluster *asymmetry* in cell units vs N cells in cluster",
4852 nceclbins,nceclmin,nceclmax,21,-1.05,1.05);
4858 fhDeltaIEtaDeltaIPhiE0[1] =
new TH2F (
"hDeltaIEtaDeltaIPhiE0Charged",
" Cluster size in columns vs rows for E < 2 GeV, #it{n}_{cells} > 3, matched with track",
4864 fhDeltaIEtaDeltaIPhiE2[1] =
new TH2F (
"hDeltaIEtaDeltaIPhiE2Charged",
" Cluster size in columns vs rows for 2 <E < 6 GeV, #it{n}_{cells} > 3, matched with track",
4870 fhDeltaIEtaDeltaIPhiE6[1] =
new TH2F (
"hDeltaIEtaDeltaIPhiE6Charged",
" Cluster size in columns vs rows for E > 6 GeV, #it{n}_{cells} > 3, matched with track",
4876 fhDeltaIA[1] =
new TH2F (
"hDeltaIACharged",
" Cluster *asymmetry* in cell units vs E, matched with track",
4877 nptbins,ptmin,ptmax,21,-1.05,1.05);
4878 fhDeltaIA[1]->SetXTitle(
"#it{E}_{cluster}");
4879 fhDeltaIA[1]->SetYTitle(
"#it{A}_{cell in cluster}");
4882 fhDeltaIAL0[1] =
new TH2F (
"hDeltaIAL0Charged",
" Cluster *asymmetry* in cell units vs #lambda^{2}_{0}, matched with track",
4883 ssbins,ssmin,ssmax,21,-1.05,1.05);
4885 fhDeltaIAL0[1]->SetYTitle(
"#it{A}_{cell in cluster}");
4888 fhDeltaIAL1[1] =
new TH2F (
"hDeltaIAL1Charged",
" Cluster *asymmetry* in cell units vs #lambda^{2}_{1}, matched with track",
4889 ssbins,ssmin,ssmax,21,-1.05,1.05);
4891 fhDeltaIAL1[1]->SetYTitle(
"#it{A}_{cell in cluster}");
4894 fhDeltaIANCells[1] =
new TH2F (
"hDeltaIANCellsCharged",
" Cluster *asymmetry* in cell units vs N cells in cluster, matched with track",
4895 nceclbins,nceclmin,nceclmax,21,-1.05,1.05);
4901 TString particle[]={
"Photon",
"Electron",
"Conversion",
"Hadron"};
4902 for (
Int_t iPart = 0; iPart < 4; iPart++) {
4904 fhDeltaIAMC[iPart] =
new TH2F (Form(
"hDeltaIA_MC%s",particle[iPart].
Data()),Form(
" Cluster *asymmetry* in cell units vs E, from %s",particle[iPart].
Data()),
4905 nptbins,ptmin,ptmax,21,-1.05,1.05);
4906 fhDeltaIAMC[iPart]->SetXTitle(
"#it{E}_{cluster}");
4907 fhDeltaIAMC[iPart]->SetYTitle(
"#it{A}_{cell in cluster}");