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), fFillClusterMaxCellHisto(kFALSE),
56 fFillAllPi0Histo(kTRUE), fFillInvMassOpenAngle(kFALSE),
57 fFillPi0PairDiffTime(kFALSE), fFillInvMassInEMCALWithPHOSDCalAcc(kFALSE),
58 fFillAllCellHistograms(kTRUE), fFillAllCellAbsIdHistograms(kTRUE),
59 fFillEBinAcceptanceHisto(kFALSE),
60 fCorrelate(kTRUE), fStudyBadClusters(kFALSE),
61 fStudyClustersAsymmetry(kFALSE), fStudyExotic(kFALSE),
62 fStudyWeight(kFALSE), fStudyTCardCorrelation(kFALSE),
63 fStudyM02Dependence (kFALSE),
66 fNModules(12), fNRCU(2),
67 fNMaxCols(48), fNMaxRows(24),
68 fTimeCutMin(-10000), fTimeCutMax(10000),
70 fEMCALCellAmpMin(0), fPHOSCellAmpMin(0),
71 fEMCALClusterM02Min(0),
72 fEMCALClusterNCellMin(0), fPHOSClusterNCellMin(0),
76 fInvMassMinECut(0), fInvMassMaxECut(0),
77 fInvMassMinM02Cut(0), fInvMassMaxM02Cut(0),
78 fInvMassMaxOpenAngle(0), fInvMassMaxTimeDifference(0),
84 fClusterMomentum(), fClusterMomentum2(),
86 fConstantTimeShift(0),
91 fhEtaPhi(0), fhEtaPhiE(0),
92 fhECharged(0), fhPtCharged(0),
93 fhPhiCharged(0), fhEtaCharged(0),
94 fhEtaPhiCharged(0), fhEtaPhiECharged(0),
97 fhIM(0), fhIMSame(0), fhIMDiff(0),
98 fhIMDCAL(0), fhIMDCALSame(0), fhIMDCALDiff(0),
99 fhIMDCALPHOS(0), fhIMDCALPHOSSame(0),
100 fhIMEMCALPHOS(0), fhIMEMCALPHOSSame(0),
102 fhOpAngle(0), fhIMvsOpAngle(0),
103 fhNCellsPerCluster(0), fhNCellsPerClusterNoCut(0),
107 fhClusterTimeEnergy(0), fhCellTimeSpreadRespectToCellMax(0),
108 fhCellIdCellLargeTimeSpread(0), fhClusterPairDiffTimeE(0), fhClusterPairDiffTimeESameMod(0),
109 fhClusterMaxCellCloseCellRatio(0), fhClusterMaxCellCloseCellDiff(0),
110 fhClusterMaxCellDiff(0), fhClusterMaxCellDiffNoCut(0),
112 fhClusterMaxCellECross(0),
114 fhLambda0(0), fhLambda1(0), fhNLocMax(0),
116 fhEnergyTMEtaResidual1Cell(0), fhEnergyTMPhiResidual1Cell(0),
117 fhColRowExoticHighE1CellPosTime(0), fhColRowExoticHighE1CellNegTime(0), fhColRowExoticHighE1CellNulTime(0),
118 fhEnergyTMEtaResidualExotic(0), fhEnergyTMPhiResidualExotic(0),
119 fhColRowExoticHighEPosTime(0), fhColRowExoticHighENegTime(0), fhColRowExoticHighENulTime(0),
120 fhColRowHighEPosTime(0), fhColRowHighENegTime(0), fhColRowHighENulTime(0),
121 fhEnergyTMEtaResidualTCardCorrNoSelection1Cell(0), fhEnergyTMPhiResidualTCardCorrNoSelection1Cell(0),
122 fhEnergyTMEtaResidualTCardCorrNoSelectionExotic(0), fhEnergyTMPhiResidualTCardCorrNoSelectionExotic(0),
126 fhBadClusterEnergy(0), fhBadClusterTimeEnergy(0), fhBadClusterEtaPhi(0),
127 fhBadClusterPairDiffTimeE(0), fhBadCellTimeSpreadRespectToCellMax(0),
128 fhBadClusterMaxCellCloseCellRatio(0), fhBadClusterMaxCellCloseCellDiff(0), fhBadClusterMaxCellDiff(0),
130 fhBadClusterMaxCellECross(0),
131 fhBadClusterLambda0(0), fhBadClusterLambda1(0),
132 fhBadClusterDeltaIEtaDeltaIPhiE0(0), fhBadClusterDeltaIEtaDeltaIPhiE2(0),
133 fhBadClusterDeltaIEtaDeltaIPhiE6(0), fhBadClusterDeltaIA(0),
136 fhRNCells(0), fhXNCells(0),
137 fhYNCells(0), fhZNCells(0),
141 fhRCellE(0), fhXCellE(0),
142 fhYCellE(0), fhZCellE(0),
144 fhDeltaCellClusterRNCells(0), fhDeltaCellClusterXNCells(0),
145 fhDeltaCellClusterYNCells(0), fhDeltaCellClusterZNCells(0),
146 fhDeltaCellClusterRE(0), fhDeltaCellClusterXE(0),
147 fhDeltaCellClusterYE(0), fhDeltaCellClusterZE(0),
150 fhNCells(0), fhNCellsCutAmpMin(0),
151 fhAmplitude(0), fhAmpId(0),
152 fhEtaPhiAmpCell(0), fhEtaPhiCell(0),
154 fhTimeId(0), fhTimeL1UnCorrId(0), fhTimeAmp(0),
155 fhAmpIdLowGain(0), fhTimeIdLowGain(0), fhTimeAmpLowGain(0),
159 fhEMCALPHOSCorrNClusters(0), fhEMCALPHOSCorrEClusters(0),
160 fhEMCALPHOSCorrNCells(0), fhEMCALPHOSCorrECells(0),
161 fhEMCALDCALCorrNClusters(0), fhEMCALDCALCorrEClusters(0),
162 fhEMCALDCALCorrNCells(0), fhEMCALDCALCorrECells(0),
163 fhDCALPHOSCorrNClusters(0), fhDCALPHOSCorrEClusters(0),
164 fhDCALPHOSCorrNCells(0), fhDCALPHOSCorrECells(0),
165 fhCaloV0SCorrNClusters(0), fhCaloV0SCorrEClusters(0),
166 fhCaloV0SCorrNCells(0), fhCaloV0SCorrECells(0),
167 fhCaloV0MCorrNClusters(0), fhCaloV0MCorrEClusters(0),
168 fhCaloV0MCorrNCells(0), fhCaloV0MCorrECells(0),
169 fhCaloTrackMCorrNClusters(0), fhCaloTrackMCorrEClusters(0),
170 fhCaloTrackMCorrNCells(0), fhCaloTrackMCorrECells(0),
171 fhCaloCenNClusters(0), fhCaloCenEClusters(0),
172 fhCaloCenNCells(0), fhCaloCenECells(0),
173 fhCaloEvPNClusters(0), fhCaloEvPEClusters(0),
174 fhCaloEvPNCells(0), fhCaloEvPECells(0),
177 fhEMod(0), fhAmpMod(0),
178 fhEWeirdMod(0), fhAmpWeirdMod(0),
180 fhNClustersMod(0), fhNCellsMod(0),
181 fhNCellsSumAmpPerMod(0), fhNClustersSumEnergyPerMod(0),
182 fhNCellsPerClusterMod(0), fhNCellsPerClusterModNoCut(0),
183 fhNCellsPerClusterWeirdMod(0), fhNCellsPerClusterWeirdModNoCut(0),
185 fhGridCells(0), fhGridCellsE(0), fhGridCellsTime(0),
186 fhGridCellsLowGain(0), fhGridCellsELowGain(0), fhGridCellsTimeLowGain(0),
187 fhTimeAmpPerRCU(0), fhIMMod(0),
190 fhECellClusterRatio(0), fhECellClusterLogRatio(0),
191 fhEMaxCellClusterRatio(0), fhEMaxCellClusterLogRatio(0),
192 fhECellTotalRatio(0), fhECellTotalLogRatio(0),
193 fhECellTotalRatioMod(0), fhECellTotalLogRatioMod(0),
198 fhRecoMCE(), fhRecoMCPhi(), fhRecoMCEta(),
199 fhRecoMCDeltaE(), fhRecoMCRatioE(),
200 fhRecoMCDeltaPhi(), fhRecoMCDeltaEta(),
203 fhGenMCE(), fhGenMCPt(), fhGenMCEtaPhi(),
204 fhGenMCAccE(), fhGenMCAccPt(), fhGenMCAccEtaPhi(),
207 fhEMVxyz(0), fhEMR(0),
208 fhHaVxyz(0), fhHaR(0),
209 fh1EOverP(0), fh2dR(0),
210 fh2EledEdx(0), fh2MatchdEdx(0),
211 fh1EOverPR02(0), fh1EleEOverP(0),
212 fhMCEle1EOverP(0), fhMCEle1dR(0), fhMCEle2MatchdEdx(0),
213 fhMCChHad1EOverP(0), fhMCChHad1dR(0), fhMCChHad2MatchdEdx(0),
214 fhMCNeutral1EOverP(0), fhMCNeutral1dR(0), fhMCNeutral2MatchdEdx(0),
215 fhMCEle1EOverPR02(0), fhMCChHad1EOverPR02(0), fhMCNeutral1EOverPR02(0),
216 fhMCEle1EleEOverP(0), fhMCChHad1EleEOverP(0), fhMCNeutral1EleEOverP(0),
217 fhTrackMatchedDEtaNeg(0), fhTrackMatchedDPhiNeg(0), fhTrackMatchedDEtaDPhiNeg(0),
218 fhTrackMatchedDEtaPos(0), fhTrackMatchedDPhiPos(0), fhTrackMatchedDEtaDPhiPos(0),
219 fhTrackMatchedDEtaNegMod(0), fhTrackMatchedDPhiNegMod(0),
220 fhTrackMatchedDEtaPosMod(0), fhTrackMatchedDPhiPosMod(0),
221 fhClusterTimeEnergyM02(0), fhCellTimeSpreadRespectToCellMaxM02(0),
222 fhClusterMaxCellCloseCellRatioM02(0), fhClusterMaxCellCloseCellDiffM02(0),
223 fhClusterMaxCellDiffM02(0), fhClusterMaxCellECrossM02(0), fhNCellsPerClusterM02(0)
226 for(
Int_t i =0; i < 12; i++)
228 for(
Int_t j = 0; j < 4; j++)
230 for(
Int_t k = 0; k < 3; k++)
237 for(
Int_t j = 0; j < 5; j++)
271 for(
Int_t i = 0; i < 7; i++)
282 for(
Int_t i = 0; i < 4; i++)
292 for(
Int_t i = 0; i < 14; i++)
302 for(
Int_t tm = 0; tm < 2; tm++)
398 for(
Int_t i = 0; i < 6; i++)
400 for(
Int_t j = 0; j < 6; j++)
473 for(
Int_t i = 0; i < 14; i++)
492 for(
Int_t j = 0; j < 6; j++)
494 for(
Int_t k = 0; k < 6; k++)
515 for(
Int_t i = 0; i < 12; i++)
578 if(phi < 0) phi += TMath::TwoPi();
591 for(
Int_t iclus2 = 0; iclus2 < caloClusters->GetEntriesFast(); iclus2++ )
593 AliVCluster* clus2 = (AliVCluster*) caloClusters->At(iclus2);
595 if(clus->GetID() == clus2->GetID())
continue;
600 if(
IsGoodCluster(absIdMax2, clus->GetM02(), clus->GetNCells(), cells) && clus2->GetM02() > 0.1 )
602 Double_t tof2 = clus2->GetTOF()*1.e9;
623 for (
Int_t ipos = 0; ipos < clus->GetNCells(); ipos++)
625 Int_t absId = clus->GetCellsAbsId()[ipos];
626 if(absId!=absIdMax && cells->GetCellAmplitude(absIdMax) > 0.01)
628 Float_t frac = cells->GetCellAmplitude(absId)/cells->GetCellAmplitude(absIdMax);
635 Double_t time = cells->GetCellTime(absId);
654 AliVCaloCells* cells,
662 for (
Int_t ipos = 0; ipos < clus->GetNCells(); ipos++)
664 Int_t id = clus->GetCellsAbsId()[ipos];
667 amp = cells->GetCellAmplitude(
id);
686 Int_t ncells = clus->GetNCells();
688 for (
Int_t ipos = 0; ipos < ncells; ipos++)
690 id = clus ->GetCellsAbsId()[ipos];
691 amp = cells->GetCellAmplitude(
id);
692 time = cells->GetCellTime(
id);
700 wTime += time*1e9 * w;
704 if(ncells > 0) aTime /= ncells;
707 if(wTot > 0) wTime /= wTot;
710 timeAverages[0] = aTime;
711 timeAverages[1] = wTime;
721 Int_t ncells = cells->GetNumberOfCells();
735 nCellsInModule[imod] = 0 ;
736 eCellsInModule[imod] = 0.;
739 Int_t icol = -1, icolAbs = -1;
740 Int_t irow = -1, irowAbs = -1;
750 for (
Int_t iCell = 0; iCell < cells->GetNumberOfCells(); iCell++)
752 if ( cells->GetCellNumber(iCell) < 0 )
continue;
754 AliDebug(2,Form(
"Cell : amp %f, absId %d", cells->GetAmplitude(iCell), cells->GetCellNumber(iCell)));
760 AliDebug(2,Form(
"\t module %d, column %d (%d), row %d (%d)", nModule,icolAbs,icol,irowAbs,irow));
762 if(nModule < fNModules)
769 if(
GetCaloUtils()->GetEMCALChannelStatus(nModule,icol,irow))
continue;
773 if(
GetCaloUtils()->GetPHOSChannelStatus(nModule,icol,irow) )
continue;
777 amp = cells->GetAmplitude(iCell)*recalF;
778 time = cells->GetTime(iCell);
779 id = cells->GetCellNumber(iCell);
780 highG = cells->GetCellHighGain(
id);
794 timeL1UnCorr *= 1.0e9;
799 AliDebug(1,Form(
"Remove cell with Time %f",time));
805 GetCaloUtils()->GetEMCALRecoUtils()->IsExoticCell(
id, cells, bc))
continue;
808 if(amp>=10)binWidthCorrection=1.0/4;
809 if(amp>=20)binWidthCorrection=1.0/10;
843 nCellsInModule[nModule]++ ;
902 if ( cellphi < 0 ) cellphi+=TMath::TwoPi();
916 Float_t rcell = TMath::Sqrt(cellpos[0]*cellpos[0]+cellpos[1]*cellpos[1]);
924 Int_t relId[4], module;
932 Float_t rcell = TMath::Sqrt(xyz.X()*xyz.X()+xyz.Y()*xyz.Y());
952 AliDebug(1,Form(
"Module %d, calo %s, N cells %d, sum Amp %f", imod,
GetCalorimeterString().
Data(), nCellsInModule[imod], eCellsInModule[imod]));
963 for (
Int_t iCell = 0; iCell < cells->GetNumberOfCells(); iCell++)
965 if ( cells->GetCellNumber(iCell) < 0 )
continue;
967 AliDebug(2,Form(
"Cell : amp %f, absId %d", cells->GetAmplitude(iCell), cells->GetCellNumber(iCell)));
971 AliDebug(2,Form(
"\t module %d, column %d, row %d", nModule,icol,irow));
973 if(nModule < fNModules)
980 if(
GetCaloUtils()->GetEMCALChannelStatus(nModule, icol, irow))
continue;
984 if(
GetCaloUtils()->GetPHOSChannelStatus (nModule, icol, irow) )
continue;
988 amp = cells->GetAmplitude(iCell)*recalF;
989 time = cells->GetTime(iCell);
990 id = cells->GetCellNumber(iCell);
997 GetReader()->GetInputEvent()->GetBunchCrossNumber());
1005 AliDebug(1,Form(
"Remove cell with Time %f",time));
1011 GetCaloUtils()->GetEMCALRecoUtils()->IsExoticCell(
id, cells, bc))
continue;
1018 Float_t ratio = amp/ecellsCut;
1023 if(eCellsInModule[nModule] > 0)
1025 Float_t ratioMod = amp/eCellsInModule[nModule];
1034 delete [] nCellsInModule;
1035 delete [] eCellsInModule;
1043 Int_t nCaloCellsPerCluster = clus->GetNCells();
1045 UShort_t * indexList = clus->GetCellsAbsId();
1048 clus->GetPosition(pos);
1053 for (
Int_t ipos = 0; ipos < nCaloCellsPerCluster; ipos++)
1056 Int_t absId = indexList[ipos];
1073 Float_t r = TMath::Sqrt(pos[0] *pos[0] + pos[1] * pos[1] );
1074 Float_t rcell = TMath::Sqrt(cellpos[0]*cellpos[0] + cellpos[1]* cellpos[1]);
1082 Int_t relId[4], module;
1098 Float_t r = TMath::Sqrt(pos[0] * pos[0] + pos[1] * pos[1] );
1099 Float_t rcell = TMath::Sqrt(xyz.X() * xyz.X() + xyz.Y() * xyz.Y());
1123 Int_t icol = -1, irow = -1, iRCU = -1, icolAbs = -1, irowAbs = -1;
1133 Int_t ncells = clus->GetNCells();
1134 Double_t time = clus->GetTOF()*1.e9;
1136 Float_t deta = clus->GetTrackDz();
1137 Float_t dphi = clus->GetTrackDx();
1145 if(energy >= 5 && energy < 8)
1147 else if(energy >= 8)
1160 else if(exoticity > 0.97)
1162 if(energy >= 5 && energy < 8)
1164 else if(energy >= 8)
1177 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1179 Int_t absId = clus->GetCellsAbsId()[ipos];
1181 Float_t eCell = cells->GetCellAmplitude(absId);
1187 if( absId == absIdMax || weight < 0.01 )
continue;
1189 Int_t rowDiff = -100, colDiff = -100;
1193 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
1198 if(energy >= 5 && energy < 8)
1200 else if(energy >= 8)
1205 if(energy >= 5 && energy < 8)
1207 else if(energy >= 8)
1213 else if ( energy > 8 )
1226 else if ( energy > 8 && ncells > 1 && exoticity < 0.97)
1236 if ( clus->GetDistanceToBadChannel() < 5 )
return ;
1239 Int_t etaRegion = -1, phiRegion = -1;
1242 if ( etaRegion !=0 ) return ;
1252 else if(exoticity > 0.97)
1262 Int_t absIdList[ncells];
1270 Int_t nCellWithWeight = 1;
1278 Float_t eCellMax = cells->GetCellAmplitude(absIdMax);
1279 Double_t tCellMax = cells->GetCellTime(absIdMax);
1297 Bool_t sameTCard2ndMax = kFALSE;
1301 Float_t eCellSameRowSameTCardNearCol = 0.;
1302 Float_t eCellSameRowDiffTCardNearCol = 0.;
1303 Double_t tCellSameRowSameTCardNearCol = 0.;
1304 Double_t tCellSameRowDiffTCardNearCol = 0.;
1311 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1313 Int_t absId = clus->GetCellsAbsId()[ipos];
1315 Float_t eCell = cells->GetCellAmplitude(absId);
1316 Double_t tCell = cells->GetCellTime(absId);
1326 if( absId == absIdMax || weight < 0.01 )
continue;
1347 Int_t rowDiff = -100, colDiff = -100;
1378 Int_t indexType = -1;
1383 if(TMath::Abs(rowDiff) == 1) nearRow = kTRUE;
1384 if(TMath::Abs(colDiff) == 1) nearCol = kTRUE;
1386 if ( rowDiff == 0 && colDiff != 0 )
1388 if ( nearCol ) indexType = 6;
1394 else if ( rowDiff != 0 && colDiff == 0 )
1396 if ( nearRow ) indexType = 8;
1404 if ( nearRow && nearCol) indexType = 10;
1405 else indexType = 11;
1415 if ( rowDiff == 0 && colDiff != 0 )
1417 if ( nearCol ) indexType = 0;
1421 else if ( rowDiff != 0 && colDiff == 0 )
1423 if ( nearRow ) indexType = 2;
1428 if ( nearCol && nearRow ) indexType = 4;
1433 if ( rowDiff == 0 && TMath::Abs(colDiff) == 1 )
1437 eCellSameRowSameTCardNearCol = eCell;
1438 tCellSameRowSameTCardNearCol = tCell;
1442 eCellSameRowDiffTCardNearCol = eCell;
1443 tCellSameRowDiffTCardNearCol = tCell;
1449 Float_t eCellDiff = eCellMax - eCell;
1450 Float_t eClusDiff = energy - eCell;
1453 Float_t tCellDiff = tCellMax - tCell;
1474 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
1479 if(energy >= 5 && energy < 8)
1481 else if(energy >= 8)
1486 if(energy >= 5 && energy < 8)
1488 else if(energy >= 8)
1493 if ( eCell > emax2nd )
1499 sameTCard2ndMax = kTRUE;
1505 sameTCard2ndMax = kFALSE;
1540 Float_t ratioNcells = nCellWithWeight/(ncells*1.);
1545 if ( nCellWithWeight <= 1 )
return;
1553 if(energy >= 5 && energy < 8)
1555 else if(energy >= 8)
1560 if(energy >= 5 && energy < 8)
1562 else if(energy >= 8)
1625 if(eCellSameRowSameTCardNearCol > 0 && eCellSameRowDiffTCardNearCol > 0)
1627 Float_t eDiff = eCellSameRowSameTCardNearCol - eCellSameRowDiffTCardNearCol ;
1628 Float_t tDiff = tCellSameRowSameTCardNearCol - tCellSameRowDiffTCardNearCol ;
1645 Int_t nCorrInd = nCorr;
1646 if(nCorr > 4) nCorrInd = 5;
1648 Int_t nCorrNoInd = nCorrNo;
1649 if(nCorrNoInd > 4) nCorrNoInd = 5;
1699 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1701 Int_t absId = clus->GetCellsAbsId()[ipos];
1703 Float_t eCell = cells->GetCellAmplitude(absId);
1704 Double_t tCell = cells->GetCellTime(absId);
1714 if( absId == absIdMax || weight < 0.01 )
continue;
1716 Float_t tDiffMaxSecondary = tCellMax - tCell;
1723 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
1726 Int_t rowDiff = -100, colDiff = -100;
1731 if ( energy >= 5 && energy < 8 && nCorr == 0 )
1733 else if ( energy >= 8 && nCorr == 0 )
1738 if ( energy >= 5 && energy < 8 && nCorrNo == 0 )
1740 else if ( energy >= 8 && nCorrNo == 0 )
1775 Double_t tof2 = clus2->GetTOF()*1.e9;
1783 if( absIdMax == absIdMax2
1784 || !
IsGoodCluster(absIdMax2, clus2->GetM02(), clus2->GetNCells(), cells)
1816 if ( energy >= 5 && energy < 8)
1818 else if ( energy >= 8 )
1836 if ( energy >= 5 && energy < 8)
1838 else if ( energy >= 8 )
1904 if ( exoticity > 0.97 )
1919 Float_t nCellRat = nCorr*1. / ((nCorr+nCorrNo)*1.);
2010 if( clus->GetNCells() <= 3 )
return;
2018 for (
Int_t ipos = 0; ipos < clus->GetNCells(); ipos++)
2020 Int_t absId = clus->GetCellsAbsId()[ipos];
2025 if(dIphi < TMath::Abs(iphi-iphiMax)) dIphi = TMath::Abs(iphi-iphiMax);
2029 if(dIeta < TMath::Abs(ieta-ietaMax)) dIeta = TMath::Abs(ieta-ietaMax);
2033 Int_t ietaShift = ieta;
2034 Int_t ietaMaxShift = ietaMax;
2035 if (ieta > ietaMax) ietaMaxShift+=48;
2036 else ietaShift +=48;
2037 if(dIeta < TMath::Abs(ietaShift-ietaMaxShift)) dIeta = TMath::Abs(ietaShift-ietaMaxShift);
2041 Float_t dIA = 1.*(dIphi-dIeta)/(dIeta+dIphi);
2062 Int_t nLabel = clus->GetNLabels();
2063 Int_t* labels = clus->GetLabels();
2112 Double_t tof = clus->GetTOF()*1.e9;
2137 Int_t nCaloCellsPerCluster = clus->GetNCells();
2142 for(
Int_t iclus2 = 0; iclus2 < caloClusters->GetEntriesFast(); iclus2++ )
2144 AliVCluster* clus2 = (AliVCluster*) caloClusters->At(iclus2);
2146 if( clus->GetID() == clus2->GetID() )
continue;
2148 if( clus->GetM02() > 0.01 && clus2->GetM02() > 0.01 )
2152 Double_t tof2 = clus2->GetTOF()*1.e9;
2158 if ( nModule2 == nModule )
2164 if(nCaloCellsPerCluster > 1 &&
2179 for (
Int_t ipos = 0; ipos < nCaloCellsPerCluster; ipos++)
2181 Int_t absId = clus->GetCellsAbsId()[ipos];
2182 if( absId == absIdMax || cells->GetCellAmplitude(absIdMax) < 0.01 )
continue;
2184 Float_t frac = cells->GetCellAmplitude(absId)/cells->GetCellAmplitude(absIdMax);
2200 Double_t time = cells->GetCellTime(absId);
2212 TMath::Abs(TMath::Abs(diff) > 100) && clus->E() > 1 )
2224 AliDebug(1,Form(
"cluster: E %2.3f, pT %2.3f, eta %2.3f, phi %2.3f",e,pt,eta,phi*TMath::RadToDeg()));
2239 Int_t icol = -1, irow = -1, iRCU = -1, icolAbs = -1, irowAbs = -1;
2268 clus->GetPosition(pos);
2280 Float_t rxyz = TMath::Sqrt(pos[0]*pos[0]+pos[1]*pos[1]);
2303 AliVCaloCells* cells)
2306 Int_t *labels = 0x0;
2307 Int_t nCaloClusters = caloClusters->GetEntriesFast() ;
2308 Int_t nCaloClustersAccepted = 0 ;
2309 Int_t nCaloCellsPerCluster = 0 ;
2321 nClustersInModule[imod] = 0;
2322 energyInModule [imod] = 0;
2328 for(
Int_t iclus = 0; iclus < nCaloClusters; iclus++)
2330 AliDebug(1,Form(
"Cluster: %d/%d, data %d",iclus+1,nCaloClusters,
GetReader()->GetDataType()));
2332 AliVCluster* clus = (AliVCluster*) caloClusters->At(iclus);
2339 Double_t tof = clus->GetTOF()*1.e9;
2344 AliDebug(1,Form(
"Remove cluster with TOF %f",tof));
2357 nLabel = clus->GetNLabels();
2358 labels = clus->GetLabels();
2364 nCaloCellsPerCluster = clus->GetNCells();
2370 Double_t tmax = cells->GetCellTime(absIdMax);
2383 if(nModule >=0 && nModule < fNModules)
2393 Float_t ampMax = cells->GetCellAmplitude(absIdMax);
2399 Bool_t goodCluster =
IsGoodCluster(absIdMax, clus->GetM02(), nCaloCellsPerCluster, cells);
2402 if(ampMax > 0.01) eCrossFrac = 1-
GetECross(absIdMax,cells)/ampMax;
2404 AliDebug(1,Form(
"Accept cluster? %d",goodCluster));
2417 maxCellFraction, eCrossFrac, tmax);
2419 nCaloClustersAccepted++;
2428 nClustersInModule[nModule]++;
2430 energyInModule [nModule] += clus->E();
2440 Int_t mcOK = kFALSE;
2442 if(
IsDataMC() && nLabel > 0 && labels)
2454 && nCaloClusters > 1
2455 && nCaloCellsPerCluster > 1
2456 &&
GetCaloUtils()->GetNumberOfLocalMaxima(clus,cells) == 1
2472 AliDebug(1,Form(
"Module %d calo %s clusters %d, sum E %f", imod,
GetCalorimeterString().
Data(), nClustersInModule[imod], energyInModule[imod]));
2480 delete [] nClustersInModule;
2481 delete [] energyInModule;
2491 if(!labels || nLabels<=0)
2493 AliWarning(Form(
"Strange, labels array %p, n labels %d", labels,nLabels));
2497 AliDebug(1,Form(
"Primaries: nlabels %d",nLabels));
2502 if(phi < 0) phi +=TMath::TwoPi();
2504 AliAODMCParticle * aodprimary = 0x0;
2505 TParticle * primary = 0x0;
2508 Int_t label = labels[0];
2512 AliDebug(1,Form(
" *** bad label ***: label %d", label));
2531 AliDebug(1,Form(
"*** large label ***: label %d, n tracks %d", label,
GetMCStack()->GetNtrack()));
2537 pdg0 = TMath::Abs(primary->GetPdgCode());
2539 status = primary->GetStatusCode();
2540 vxMC = primary->Vx();
2541 vyMC = primary->Vy();
2542 iParent = primary->GetFirstMother();
2544 AliDebug(1,
"Cluster most contributing mother:");
2545 AliDebug(1,Form(
"\t Mother label %d, pdg %d, %s, status %d, parent %d",iMother, pdg0, primary->GetName(),status, iParent));
2553 pdg = TMath::Abs(primary->GetPdgCode());
2555 AliDebug(2,
"Converted cluster!. Find before conversion:");
2557 while((pdg == 22 || pdg == 11) && status != 1)
2559 Int_t iMotherOrg = iMother;
2562 status = primary->GetStatusCode();
2563 pdg = TMath::Abs(primary->GetPdgCode());
2564 iParent = primary->GetFirstMother();
2568 if(pdg==111 || pdg==221)
2570 primary =
GetMCStack()->Particle(iMotherOrg);
2580 AliDebug(2,Form(
"\t pdg %d, index %d, %s, status %d",pdg, iMother, primary->GetName(),status));
2583 AliDebug(1,
"Converted Cluster mother before conversion:");
2584 AliDebug(1,Form(
"\t Mother label %d, pdg %d, %s, status %d, parent %d",iMother, pdg, primary->GetName(), status, iParent));
2592 AliDebug(2,
"Overlapped Meson decay!, Find it:");
2594 while(pdg != 111 && pdg != 221)
2599 status = primary->GetStatusCode();
2600 pdg = TMath::Abs(primary->GetPdgCode());
2601 iParent = primary->GetFirstMother();
2603 if( iParent < 0 )
break;
2605 AliDebug(2,Form(
"\t pdg %d, %s, index %d",pdg, primary->GetName(),iMother));
2609 AliWarning(
"Tagged as Overlapped photon but meson not found, why?");
2614 AliDebug(2,Form(
"Overlapped %s decay, label %d",primary->GetName(),iMother));
2617 eMC = primary->Energy();
2619 phiMC = primary->Phi();
2620 etaMC = primary->Eta();
2621 pdg = TMath::Abs(primary->GetPdgCode());
2622 charge = (
Int_t) TDatabasePDG::Instance()->GetParticle(pdg)->Charge();
2625 else if(
GetReader()->ReadAODMCParticles() &&
2631 AliFatal(
"MCParticles not available!");
2635 pdg0 = TMath::Abs(aodprimary->GetPdgCode());
2637 status = aodprimary->IsPrimary();
2638 vxMC = aodprimary->Xv();
2639 vyMC = aodprimary->Yv();
2640 iParent = aodprimary->GetMother();
2642 AliDebug(1,
"Cluster most contributing mother:");
2643 AliDebug(1,Form(
"\t Mother label %d, pdg %d, Primary? %d, Physical Primary? %d, parent %d",
2644 iMother, pdg0, aodprimary->IsPrimary(), aodprimary->IsPhysicalPrimary(), iParent));
2649 AliDebug(2,
"Converted cluster!. Find before conversion:");
2653 pdg = TMath::Abs(aodprimary->GetPdgCode());
2655 while ((pdg == 22 || pdg == 11) && !aodprimary->IsPhysicalPrimary())
2657 Int_t iMotherOrg = iMother;
2660 status = aodprimary->IsPrimary();
2661 iParent = aodprimary->GetMother();
2662 pdg = TMath::Abs(aodprimary->GetPdgCode());
2666 if( pdg == 111 || pdg == 221 )
2678 AliDebug(2,Form(
"\t pdg %d, index %d, Primary? %d, Physical Primary? %d",
2679 pdg, iMother, aodprimary->IsPrimary(), aodprimary->IsPhysicalPrimary()));
2682 AliDebug(1,
"Converted Cluster mother before conversion:");
2683 AliDebug(1,Form(
"\t Mother label %d, pdg %d, parent %d, Primary? %d, Physical Primary? %d",
2684 iMother, pdg, iParent, aodprimary->IsPrimary(), aodprimary->IsPhysicalPrimary()));
2692 AliDebug(2,Form(
"Overlapped Meson decay!, Find it: PDG %d, mom %d",pdg, iMother));
2694 while(pdg != 111 && pdg != 221)
2698 status = aodprimary->IsPrimary();
2699 iParent = aodprimary->GetMother();
2700 pdg = TMath::Abs(aodprimary->GetPdgCode());
2702 if( iParent < 0 )
break;
2704 AliDebug(2,Form(
"\t pdg %d, index %d",pdg, iMother));
2708 AliWarning(
"Tagged as Overlapped photon but meson not found, why?");
2713 AliDebug(2,Form(
"Overlapped %s decay, label %d",aodprimary->GetName(),iMother));
2716 status = aodprimary->IsPrimary();
2717 eMC = aodprimary->E();
2719 phiMC = aodprimary->Phi();
2720 etaMC = aodprimary->Eta();
2721 pdg = TMath::Abs(aodprimary->GetPdgCode());
2722 charge = aodprimary->Charge();
2726 Float_t rVMC = TMath::Sqrt(vxMC*vxMC + vyMC*vyMC);
2727 if( ( pdg == 22 || TMath::Abs(pdg) == 11 ) && status != 1 )
2760 else if(charge == 0)
2775 if( mcIndex >= 0 && mcIndex < 7 && e > 0.5 && eMC > 0.5)
2786 if( primary || aodprimary )
return kTRUE ;
2802 if(phi < 0) phi +=TMath::TwoPi();
2820 Double_t dedx = track->GetTPCsignal();
2821 Int_t nITS = track->GetNcls(0);
2822 Int_t nTPC = track->GetNcls(1);
2823 Bool_t positive = kFALSE;
2824 if(track) positive = (track->Charge()>0);
2827 Float_t deta = clus->GetTrackDz();
2828 Float_t dphi = clus->GetTrackDx();
2829 Double_t dR = TMath::Sqrt(dphi*dphi + deta*deta);
2832 if( TMath::Abs(dphi) < 999 )
2869 if(dedx > 60 && dedx < 100)
2879 if(e > 0.5 && tmom > 0.5)
2885 if(dR < 0.02 && eOverP > 0.6 && eOverP < 1.2
2886 && clus->GetNCells() > 1 && nITS > 3 && nTPC > 20)
2894 Double_t charge = TDatabasePDG::Instance()->GetParticle(pdg)->Charge();
2896 if(TMath::Abs(pdg) == 11)
2920 else if(charge == 0)
2945 if(!caloClustersEMCAL || !caloClustersPHOS)
2947 AliDebug(1,Form(
"PHOS (%p) or EMCAL (%p) clusters array not available, do not correlate",caloClustersPHOS,caloClustersEMCAL));
2955 if(!cellsEMCAL || !cellsPHOS)
2957 AliDebug(1,Form(
"PHOS (%p) or EMCAL (%p) cells array ot available, do not correlate",cellsPHOS,cellsEMCAL));
2966 Float_t sumClusterEnergyEMCAL = 0;
2967 Float_t sumClusterEnergyDCAL = 0;
2968 Float_t sumClusterEnergyPHOS = 0;
2972 AliVCluster* cluster = 0;
2973 for(iclus = 0 ; iclus < caloClustersEMCAL->GetEntriesFast() ; iclus++)
2975 cluster = (AliVCluster*)caloClustersEMCAL->At(iclus);
2976 Float_t energy = cluster->E();
2978 if( energy < 0.5 )
continue;
2980 if(cluster->GetCellsAbsId()[0] < 12288)
2983 sumClusterEnergyEMCAL +=
energy;
2988 sumClusterEnergyDCAL +=
energy;
2992 for(iclus = 0 ; iclus < caloClustersPHOS ->GetEntriesFast(); iclus++)
2994 cluster = (AliVCluster*) caloClustersPHOS->At(iclus);
2996 energy = cluster->E();
2998 if( energy < 0.5 )
continue;
3001 sumClusterEnergyPHOS +=
energy;
3005 Int_t ncellsEMCAL = 0 ;
3006 Int_t ncellsDCAL = 0 ;
3007 Int_t ncellsPHOS = 0;
3009 Float_t sumCellEnergyEMCAL = 0;
3010 Float_t sumCellEnergyDCAL = 0;
3011 Float_t sumCellEnergyPHOS = 0;
3013 for(icell = 0 ; icell < cellsEMCAL->GetNumberOfCells() ; icell++)
3015 Float_t amp = cellsEMCAL->GetAmplitude(icell);
3016 Int_t cellId = cellsEMCAL->GetCellNumber(icell);
3023 sumCellEnergyEMCAL += amp;
3028 sumCellEnergyDCAL += amp;
3032 for(icell = 0 ; icell < cellsPHOS->GetNumberOfCells(); icell++)
3034 Float_t amp = cellsPHOS->GetAmplitude(icell);
3035 Int_t cellId = cellsPHOS->GetCellNumber(icell);
3037 if ( cellId < 0 ) continue ;
3042 sumCellEnergyPHOS += amp;
3068 Int_t ncl = nclPHOS;
3069 Float_t sumClusterEnergy = sumClusterEnergyPHOS;
3070 Int_t ncells = ncellsPHOS;
3071 Float_t sumCellEnergy = sumCellEnergyPHOS;
3075 ncl = nclEMCAL + nclDCAL;
3076 sumClusterEnergy = sumClusterEnergyEMCAL + sumClusterEnergyDCAL;
3077 ncells = ncellsEMCAL + ncellsDCAL;
3078 sumCellEnergy = sumCellEnergyEMCAL + sumCellEnergyDCAL;
3106 AliDebug(1,
"Correlate():");
3107 AliDebug(1,Form(
"\t EMCAL: N cells %d, N clusters %d, summed E cells %f, summed E clusters %f",
3108 ncellsEMCAL,nclEMCAL, sumCellEnergyEMCAL,sumClusterEnergyEMCAL));
3109 AliDebug(1,Form(
"\t DCAL : N cells %d, N clusters %d, summed E cells %f, summed E clusters %f",
3110 ncellsDCAL,nclDCAL, sumCellEnergyDCAL,sumClusterEnergyDCAL));
3111 AliDebug(1,Form(
"\t PHOS : N cells %d, N clusters %d, summed E cells %f, summed E clusters %f",
3112 ncellsPHOS,nclPHOS,sumCellEnergyPHOS,sumClusterEnergyPHOS));
3113 AliDebug(1,Form(
"\t V0 : Signal %d, Multiplicity %d, Track Multiplicity %d", v0S,v0M,trM));
3114 AliDebug(1,Form(
"\t centrality : %f, Event plane angle %f", cen,ep));
3123 const Int_t buffersize = 255;
3124 char onePar[buffersize] ;
3126 snprintf(onePar,buffersize,
"--- AliAnaCalorimeterQA ---:") ;
3153 return new TObjString(parList) ;
3233 outputContainer->SetName(
"QAHistos") ;
3304 fhE =
new TH1F (
"hE",
"#it{E} reconstructed clusters ", nptbins*5,ptmin,ptmax*5);
3305 fhE->SetXTitle(
"#it{E} (GeV)");
3306 outputContainer->Add(
fhE);
3308 fhPt =
new TH1F (
"hPt",
"#it{p}_{T} reconstructed clusters", nptbins,ptmin,ptmax);
3309 fhPt->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3310 outputContainer->Add(
fhPt);
3312 fhPhi =
new TH1F (
"hPhi",
"#varphi reconstructed clusters ",nphibins,phimin,phimax);
3313 fhPhi->SetXTitle(
"#varphi (rad)");
3314 outputContainer->Add(
fhPhi);
3316 fhEta =
new TH1F (
"hEta",
"#eta reconstructed clusters ",netabins,etamin,etamax);
3317 fhEta->SetXTitle(
"#eta ");
3318 outputContainer->Add(
fhEta);
3324 fhEtaPhiE =
new TH3F (
"hEtaPhiE",
"#eta vs #varphi vs energy, reconstructed clusters",
3325 netabins,etamin,etamax,nphibins,phimin,phimax,nptbins,ptmin,ptmax);
3333 fhEtaPhi =
new TH2F (
"hEtaPhi",
"#eta vs #varphi for #it{E} > 0.5 GeV, reconstructed clusters",
3334 netabins,etamin,etamax,nphibins,phimin,phimax);
3336 fhEtaPhi->SetYTitle(
"#varphi (rad)");
3342 nptbins,ptmin,ptmax, ntimebins,timemin,timemax);
3350 nptbins,ptmin,ptmax, tdbins,tdmin,tdmax);
3356 nptbins,ptmin,ptmax, tdbins,tdmin,tdmax);
3368 fhLambda0 =
new TH2F (
"hLambda0",
"shower shape, #lambda^{2}_{0} vs E",
3369 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3370 fhLambda0->SetXTitle(
"#it{E}_{cluster} (GeV)");
3371 fhLambda0->SetYTitle(
"#lambda^{2}_{0}");
3374 fhLambda1 =
new TH2F (
"hLambda1",
"shower shape, #lambda^{2}_{1} vs E",
3375 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3376 fhLambda1->SetXTitle(
"#it{E}_{cluster} (GeV)");
3377 fhLambda1->SetYTitle(
"#lambda^{2}_{1}");
3381 nptbins,ptmin,ptmax,10,0,10);
3382 fhNLocMax->SetXTitle(
"#it{E}_{cluster} (GeV)");
3389 (
"hColRowHighEPosTime",
3390 "column vs row, exo < 0.97, E > 8 GeV, t > 5 ns",
3391 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3397 (
"hColRowHighENegTime",
3398 "column vs row, exo < 0.97, E > 8 GeV, t < -5 ns",
3399 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3405 (
"hColRowHighENulTime",
3406 "column vs row, exo < 0.97, E > 8 GeV, -5 < t < 5 ns",
3407 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3415 (
"hColRowExoticHighE1CellPosTime",
3416 "column vs row, 1 cell, E > 8 GeV, t > 5 ns",
3417 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3423 (
"hColRowExoticHighEPosTime",
3424 "column vs row, exo > 0.97, E > 8 GeV, t > 5 ns",
3425 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3431 (
"hColRowExoticHighE1CellNegTime",
3432 "column vs row, 1 cell, E > 8 GeV, t < -5 ns",
3433 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3439 (
"hColRowExoticHighENegTime",
3440 "column vs row, exo > 0.97, E > 8 GeV, t < -5 ns",
3441 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3447 (
"hColRowExoticHighE1CellNulTime",
3448 "column vs row, 1 cell, E > 8 GeV, -5 < t < 5 ns",
3449 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3455 (
"hColRowExoticHighENulTime",
3456 "column vs row, exo > 0.97, E > 8 GeV, -5 < t < 5 ns",
3457 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3463 TString add[] = {
"",
"TrackMatched"};
3464 for(
Int_t tm = 0; tm < 2; tm++)
3467 (Form(
"hColRowTCardCorrNoSelectionLowE%s",add[tm].
Data()),
3468 Form(
"column vs row, max E cell for TCard correlation selected clusters, 5 < E < 8 GeV %s",add[tm].
Data()),
3469 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3475 (Form(
"hColRowTCardCorrNoSelectionHighE%s",add[tm].
Data()),
3476 Form(
"column vs row, max E cell for TCard correlation selected clusters, E > 8 GeV %s",add[tm].
Data()),
3477 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3485 (Form(
"hNCellsTCardCorrNoSelection%s",add[tm].
Data()),
3486 Form(
"# custer # cells vs #it{E} %s",add[tm].
Data()),
3487 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
3493 (Form(
"hNCellsTCardCorrWithWeightNoSelection%s",add[tm].
Data()),
3494 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
3495 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
3501 (Form(
"hNCellsTCardCorrRatioWithWeightNoSelection%s",add[tm].
Data()),
3502 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
3503 nptbins,ptmin,ptmax, 100,0,1);
3509 (Form(
"hLambda0TCardCorrNoSelection%s",add[tm].
Data()),
3510 Form(
"#lambda^{2}_{0} vs #it{E} %s",add[tm].
Data()),
3511 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3517 (Form(
"hLambda1TCardCorrNoSelection%s",add[tm].
Data()),
3518 Form(
"#lambda^{2}_{1} vs #it{E} %s",add[tm].
Data()),
3519 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3525 (Form(
"hLambda0NLM1TCardCorrNoSelection%s",add[tm].
Data()),
3526 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=1 %s",add[tm].
Data()),
3527 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3533 (Form(
"hLambda1NLM1TCardCorrNoSelection%s",add[tm].
Data()),
3534 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=1 %s",add[tm].
Data()),
3535 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3541 (Form(
"hLambda0NLM2TCardCorrNoSelection%s",add[tm].
Data()),
3542 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=2 %s",add[tm].
Data()),
3543 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3549 (Form(
"hLambda1NLM2TCardCorrNoSelection%s",add[tm].
Data()),
3550 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=2 %s",add[tm].
Data()),
3551 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3558 (Form(
"hLambdaRTCardCorrNoSelection%s",add[tm].
Data()),
3559 Form(
"#lambda^{1}_{0}/#lambda^{2}_{0} vs #it{E} %s",add[tm].
Data()),
3560 nptbins,ptmin,ptmax,110,0,1.1);
3566 (Form(
"hNLocMaxTCardCorrNoSelection%s",add[tm].
Data()),
3567 Form(
"#it{n}_{LM} vs E %s",add[tm].
Data()),
3568 nptbins,ptmin,ptmax,10,0,10);
3574 (Form(
"hEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3575 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
3576 nptbins,ptmin,ptmax,100,0,1);
3582 (Form(
"hEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3583 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3584 nptbins,ptmin,ptmax,100,0,1);
3590 (Form(
"hEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3591 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3592 nptbins,ptmin,ptmax,100,0,1);
3598 (Form(
"hE2ndRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3599 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
3600 nptbins,ptmin,ptmax,100,0,1);
3606 (Form(
"hE2ndRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3607 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3608 nptbins,ptmin,ptmax,100,0,1);
3614 (Form(
"hE2ndRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3615 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3616 nptbins,ptmin,ptmax,100,0,1);
3623 (Form(
"hE2ndEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3624 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
3625 nptbins,ptmin,ptmax,100,0,1);
3631 (Form(
"hE2ndEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3632 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3633 nptbins,ptmin,ptmax,100,0,1);
3639 (Form(
"hE2ndEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3640 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3641 nptbins,ptmin,ptmax,100,0,1);
3649 (Form(
"hE2ndSameRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3650 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cluster} vs E, 2nd in same TCard as leading, #it{n}_{LM}=1 %s",add[tm].
Data()),
3651 nptbins,ptmin,ptmax,100,0,1);
3657 (Form(
"hE2ndSameRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3658 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, 2nd in same TCard as leading, #it{n}_{LM}=2 %s",add[tm].
Data()),
3659 nptbins,ptmin,ptmax,100,0,1);
3665 (Form(
"hE2ndSameRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3666 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, 2nd in same TCard as leading, #it{n}_{LM}>2 %s",add[tm].
Data()),
3667 nptbins,ptmin,ptmax,100,0,1);
3674 (Form(
"hE2ndSameEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3675 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cell}^{max} vs E, 2nd in same TCard as leading, #it{n}_{LM}=1 %s",add[tm].
Data()),
3676 nptbins,ptmin,ptmax,100,0,1);
3682 (Form(
"hE2ndSameEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3683 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, 2nd in same TCard as leading, #it{n}_{LM}=2 %s",add[tm].
Data()),
3684 nptbins,ptmin,ptmax,100,0,1);
3690 (Form(
"hE2ndSameEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3691 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, 2nd in same TCard as leading, #it{n}_{LM}>2 %s",add[tm].
Data()),
3692 nptbins,ptmin,ptmax,100,0,1);
3700 (Form(
"hECellClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3701 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
3702 nptbins,ptmin,ptmax,100,0,1);
3708 (Form(
"hECellClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3709 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3710 nptbins,ptmin,ptmax,100,0,1);
3716 (Form(
"hECellClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3717 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3718 nptbins,ptmin,ptmax,100,0,1);
3724 (Form(
"hLogECellNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3725 Form(
"log(#it{E}_{cell}) vs E_{cluster}, #it{n}_{LM}=1, w > 0.01 %s",add[tm].
Data()),
3726 nptbins,ptmin,ptmax,150,-3,3);
3732 (Form(
"hLogECellNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3733 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}=2, w > 0.01 %s",add[tm].
Data()),
3734 nptbins,ptmin,ptmax,150,-3,3);
3740 (Form(
"hLogECellNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3741 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}>2, w > 0.01 %s",add[tm].
Data()),
3742 nptbins,ptmin,ptmax,150,-3,3);
3749 (Form(
"hECellWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3750 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
3751 nptbins,ptmin,ptmax,90,0,4.5);
3757 (Form(
"hECellWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3758 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3759 nptbins,ptmin,ptmax,90,0,4.5);
3765 (Form(
"hECellWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3766 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster}))vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3767 nptbins,ptmin,ptmax,90,0,4.5);
3774 (Form(
"hECellSameClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3775 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E_{cluster}, cell from same T-Card as leading, #it{n}_{LM}=1 %s",add[tm].
Data()),
3776 nptbins,ptmin,ptmax,100,0,1);
3782 (Form(
"hECellSameClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3783 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, cell from same T-Card as leading, #it{n}_{LM}=2 %s",add[tm].
Data()),
3784 nptbins,ptmin,ptmax,100,0,1);
3790 (Form(
"hECellSameClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3791 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, cell from same T-Card as leading, #it{n}_{LM}>2 %s",add[tm].
Data()),
3792 nptbins,ptmin,ptmax,100,0,1);
3798 (Form(
"hLogECellSameNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3799 Form(
"log(#it{E}_{cell}) vs E_{cluster}, cell from same T-Card as leading, #it{n}_{LM}=1, w > 0.01 %s",add[tm].
Data()),
3800 nptbins,ptmin,ptmax,150,-3,3);
3806 (Form(
"hLogECellSameNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3807 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}=2, cell from same T-Card as leading, w > 0.01 %s",add[tm].
Data()),
3808 nptbins,ptmin,ptmax,150,-3,3);
3814 (Form(
"hLogECellSameNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3815 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}>2, cell from same T-Card as leading, w > 0.01 %s",add[tm].
Data()),
3816 nptbins,ptmin,ptmax,150,-3,3);
3823 (Form(
"hECellSameWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3824 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E_{cluster}, cell from same T-Card as leading, #it{n}_{LM}=1 %s",add[tm].
Data()),
3825 nptbins,ptmin,ptmax,90,0,4.5);
3831 (Form(
"hECellSameWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3832 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E, cell from same T-Card as leading, #it{n}_{LM}=2 %s",add[tm].
Data()),
3833 nptbins,ptmin,ptmax,90,0,4.5);
3839 (Form(
"hECellSameWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3840 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster}))vs E, cell from same T-Card as leading, #it{n}_{LM}>2 %s",add[tm].
Data()),
3841 nptbins,ptmin,ptmax,90,0,4.5);
3847 (Form(
"hExoticTCardCorrNoSelection%s",add[tm].
Data()),
3848 Form(
"exoticity vs #it{E} %s",add[tm].
Data()),
3849 nptbins,ptmin,ptmax,200,-1,1);
3855 (Form(
"hNCellsTCardSameAndDiffFraction%s",add[tm].
Data()),
3856 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, %s",add[tm].
Data()),
3857 nptbins,ptmin,ptmax,100,0,1);
3863 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterE%s",add[tm].
Data()),
3864 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()),
3865 nptbins,ptmin,ptmax,200,-10,10);
3871 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterE%s",add[tm].
Data()),
3872 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()),
3873 nptbins,ptmin,ptmax,200,-100,100);
3879 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxE%s",add[tm].
Data()),
3880 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()),
3881 nptbins,ptmin,ptmax,200,-10,10);
3887 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxE%s",add[tm].
Data()),
3888 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()),
3889 nptbins,ptmin,ptmax,200,-100,100);
3897 (Form(
"hEnergyTime1Cell%s",add[tm].
Data()),
3898 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
3899 nptbins,ptmin,ptmax,300,-150,150);
3905 (Form(
"hEnergyTimeExotic%s",add[tm].
Data()),
3906 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
3907 nptbins,ptmin,ptmax,300,-150,150);
3913 (Form(
"hEnergyTimeTCardCorrNoSelection1Cell%s",add[tm].
Data()),
3914 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
3915 nptbins,ptmin,ptmax,300,-150,150);
3921 (Form(
"hEnergyTimeTCardCorrNoSelectionExotic%s",add[tm].
Data()),
3922 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
3923 nptbins,ptmin,ptmax,300,-150,150);
3929 (Form(
"hColRowExoticLowE1Cell%s",add[tm].
Data()),
3930 Form(
"column vs row, 1 cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3931 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3937 (Form(
"hColRowExoticHighE1Cell%s",add[tm].
Data()),
3938 Form(
"column vs row, 1 cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3939 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3945 (Form(
"hColRowExoticLowE%s",add[tm].
Data()),
3946 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3947 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3953 (Form(
"hColRowExoticHighE%s",add[tm].
Data()),
3954 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3955 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3961 (Form(
"hColRowExotic2ndCellDiffLowE%s",add[tm].
Data()),
3962 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3963 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3969 (Form(
"hColRowExotic2ndCellDiffHighE%s",add[tm].
Data()),
3970 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3971 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3977 (Form(
"hColRowExotic2ndCellSameLowE%s",add[tm].
Data()),
3978 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3979 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3985 (Form(
"hColRowExotic2ndCellSameHighE%s",add[tm].
Data()),
3986 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3987 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3993 (Form(
"hColRowTCardCorrNoSelectionExoticLowE%s",add[tm].
Data()),
3994 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3995 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4001 (Form(
"hColRowTCardCorrNoSelectionExoticHighE%s",add[tm].
Data()),
4002 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
4003 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4009 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffLowE%s",add[tm].
Data()),
4010 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
4011 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4017 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffHighE%s",add[tm].
Data()),
4018 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
4019 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4025 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameLowE%s",add[tm].
Data()),
4026 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
4027 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4033 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameHighE%s",add[tm].
Data()),
4034 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
4035 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4041 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameLowE%s",add[tm].
Data()),
4042 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
4043 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4049 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameHighE%s",add[tm].
Data()),
4050 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
4051 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4057 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffLowE%s",add[tm].
Data()),
4058 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
4059 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4065 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffHighE%s",add[tm].
Data()),
4066 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
4067 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4073 (Form(
"hNCellsTCardSameAndDiffFraction_Exotic%s",add[tm].
Data()),
4074 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, exo > 0.97 %s",add[tm].
Data()),
4075 nptbins,ptmin,ptmax,100,0,1);
4081 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterEExo%s",add[tm].
Data()),
4082 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()),
4083 nptbins,ptmin,ptmax,200,-10,10);
4089 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterEExo%s",add[tm].
Data()),
4090 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()),
4091 nptbins,ptmin,ptmax,200,-100,100);
4097 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxEExo%s",add[tm].
Data()),
4098 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()),
4099 nptbins,ptmin,ptmax,200,-10,10);
4105 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxEExo%s",add[tm].
Data()),
4106 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()),
4107 nptbins,ptmin,ptmax,200,-100,100);
4113 for(
Int_t i = 0; i < 6; i++)
4115 for(
Int_t j = 0; j < 6; j++)
4118 (Form(
"hLambda0TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4119 Form(
"#lambda^{2}_{0} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4120 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4126 (Form(
"hLambda1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4127 Form(
"#lambda^{2}_{1} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4128 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4134 (Form(
"hLambda0NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4135 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()),
4136 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4142 (Form(
"hLambda1NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4143 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()),
4144 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4150 (Form(
"hLambda0NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4151 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()),
4152 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4158 (Form(
"hLambda1NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4159 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()),
4160 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4175 (Form(
"hNLocMaxTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4176 Form(
"#it{n}_{LM} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4177 nptbins,ptmin,ptmax,10,0,10);
4183 (Form(
"hEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4184 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()),
4185 nptbins,ptmin,ptmax,100,0,1);
4191 (Form(
"hEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4192 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()),
4193 nptbins,ptmin,ptmax,100,0,1);
4199 (Form(
"hEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4200 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()),
4201 nptbins,ptmin,ptmax,100,0,1);
4207 (Form(
"hE2ndRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4208 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()),
4209 nptbins,ptmin,ptmax,100,0,1);
4215 (Form(
"hE2ndRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4216 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()),
4217 nptbins,ptmin,ptmax,100,0,1);
4223 (Form(
"hE2ndRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4224 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()),
4225 nptbins,ptmin,ptmax,100,0,1);
4231 (Form(
"hE2ndEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4232 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()),
4233 nptbins,ptmin,ptmax,100,0,1);
4239 (Form(
"hE2ndEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4240 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()),
4241 nptbins,ptmin,ptmax,100,0,1);
4247 (Form(
"hE2ndEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4248 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()),
4249 nptbins,ptmin,ptmax,100,0,1);
4256 (Form(
"hECellClusRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4257 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()),
4258 nptbins,ptmin,ptmax,100,0,1);
4264 (Form(
"hECellClusRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4265 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()),
4266 nptbins,ptmin,ptmax,100,0,1);
4272 (Form(
"hECellClusRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4273 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()),
4274 nptbins,ptmin,ptmax,100,0,1);
4280 (Form(
"hLogECellNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4281 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()),
4282 nptbins,ptmin,ptmax,150,-3,3);
4288 (Form(
"hLogECellNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4289 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()),
4290 nptbins,ptmin,ptmax,150,-3,3);
4296 (Form(
"hLogECellNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4297 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()),
4298 nptbins,ptmin,ptmax,150,-3,3);
4304 (Form(
"hECellWeightNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4305 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()),
4306 nptbins,ptmin,ptmax,90,0,4.5);
4312 (Form(
"hECellWeightNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4313 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()),
4314 nptbins,ptmin,ptmax,90,0,4.5);
4320 (Form(
"hECellWeightNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4321 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()),
4322 nptbins,ptmin,ptmax,90,0,4.5);
4328 (Form(
"hMassEClusTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4329 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()),
4330 nptbins,ptmin,ptmax,nmassbins,massmin,massmax);
4344 (Form(
"hTimeDiffTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4345 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()),
4346 nptbins,ptmin,ptmax,300,-150,150);
4352 (Form(
"hColRowTCardCorrelNCellLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4353 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4354 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4360 (Form(
"hColRowTCardCorrelNCellHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4361 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4362 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4370 (Form(
"hExoticTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4371 Form(
"exoticity vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4372 nptbins,ptmin,ptmax,200,-1,1);
4378 (Form(
"hTimeDiffExoTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4379 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()),
4380 nptbins,ptmin,ptmax,300,-150,150);
4386 (Form(
"hColRowTCardCorrelNCellExoticLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4387 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4388 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4394 (Form(
"hColRowTCardCorrelNCellExoticHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4395 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4396 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4611 (Form(
"hLambda0Exoticity_EBin%d%s",i,add[tm].
Data()),
4612 Form(
"#lambda^{2}_{0} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
4613 200,-1,1,ssbins,ssmin,ssmax);
4619 (Form(
"hLambda1Exoticity_EBin%d%s",i,add[tm].
Data()),
4620 Form(
"#lambda^{2}_{1} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
4621 200,-1,1,ssbins,ssmin,ssmax);
4635 (Form(
"hNCellsExoticity_EBin%d%s",i,add[tm].
Data()),
4636 Form(
"#it{n}_{cells} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
4637 200,-1,1,nceclbins,nceclmin,nceclmax);
4638 fhNCellsExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
4643 (Form(
"hTimeExoticity_EBin%d%s",i,add[tm].
Data()),
4645 200,-1,1,100,-25,25);
4646 fhTimeExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
4651 (Form(
"hNCellsTCardSameAndDiff_Exotic_EBin%d%s",i,add[tm].
Data()),
4652 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()),
4653 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
4693 (Form(
"hNCellsTCardSameAndDiff_EBin%d%s",i,add[tm].
Data()),
4694 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()),
4695 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
4701 (Form(
"hLambda0Lambda1_EBin%d%s",i,add[tm].
Data()),
4702 Form(
"#lambda^{2}_{0} vs #lambda^{2}_{1}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
4703 ssbins,ssmin,ssmax,ssbins,ssmin,ssmax);
4719 for(
Int_t j = 0; j < 6; j++)
4721 for(
Int_t k = 0; k < 6; k++)
4724 (Form(
"hLambda0Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
4725 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()),
4726 200,-1,1,ssbins,ssmin,ssmax);
4732 (Form(
"hLambda1Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
4733 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()),
4734 200,-1,1,ssbins,ssmin,ssmax);
4831 for(
Int_t i = 0; i < 12; i++)
4834 (Form(
"hTCardCorrECellMaxDiff_Case%d%s",i,add[tm].
Data()),
4835 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()),
4836 nptbins,ptmin,ptmax,210,-1,20);
4842 (Form(
"hTCardCorrEClusterDiff_Case%d%s",i,add[tm].
Data()),
4843 Form(
"#it{E}_{cluster}-#it{E}_{cell} vs #it{E}_{cluster}, for (un)correlated cells in TCard, case %d %s",i,add[tm].
Data()),
4844 nptbins,ptmin,ptmax,210,-1,20);
4866 (Form(
"hTCardCorrTCellMaxDiff_Case%d%s",i,add[tm].
Data()),
4867 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()),
4868 nptbins,ptmin,ptmax,1000,-100,100);
4876 (Form(
"hTCardCorrECellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
4877 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()),
4878 nptbins,ptmin,ptmax,210,-1,20);
4884 (Form(
"hTCardCorrEClusterDiffExo_Case%d%s",i,add[tm].
Data()),
4885 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()),
4886 nptbins,ptmin,ptmax,210,-1,20);
4908 (Form(
"hTCardCorrTCellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
4909 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()),
4910 nptbins,ptmin,ptmax,1000,-100,100);
4921 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
4927 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
4933 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
4939 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
4945 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
4951 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
4957 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
4963 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);