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),
83 fClusterMomentum(), fClusterMomentum2(),
85 fConstantTimeShift(0),
90 fhEtaPhi(0), fhEtaPhiE(0),
91 fhECharged(0), fhPtCharged(0),
92 fhPhiCharged(0), fhEtaCharged(0),
93 fhEtaPhiCharged(0), fhEtaPhiECharged(0),
96 fhIM(0), fhIMSame(0), fhIMDiff(0),
97 fhIMDCAL(0), fhIMDCALSame(0), fhIMDCALDiff(0),
98 fhIMDCALPHOS(0), fhIMDCALPHOSSame(0),
99 fhIMEMCALPHOS(0), fhIMEMCALPHOSSame(0),
101 fhOpAngle(0), fhIMvsOpAngle(0),
102 fhNCellsPerCluster(0), fhNCellsPerClusterNoCut(0),
106 fhClusterTimeEnergy(0), fhCellTimeSpreadRespectToCellMax(0),
107 fhCellIdCellLargeTimeSpread(0), fhClusterPairDiffTimeE(0), fhClusterPairDiffTimeESameMod(0),
108 fhClusterMaxCellCloseCellRatio(0), fhClusterMaxCellCloseCellDiff(0),
109 fhClusterMaxCellDiff(0), fhClusterMaxCellDiffNoCut(0),
111 fhClusterMaxCellECross(0),
113 fhLambda0(0), fhLambda1(0), fhNLocMax(0),
115 fhEnergyTMEtaResidual1Cell(0), fhEnergyTMPhiResidual1Cell(0),
116 fhColRowExoticHighE1CellPosTime(0), fhColRowExoticHighE1CellNegTime(0), fhColRowExoticHighE1CellNulTime(0),
117 fhEnergyTMEtaResidualExotic(0), fhEnergyTMPhiResidualExotic(0),
118 fhColRowExoticHighEPosTime(0), fhColRowExoticHighENegTime(0), fhColRowExoticHighENulTime(0),
119 fhColRowHighEPosTime(0), fhColRowHighENegTime(0), fhColRowHighENulTime(0),
120 fhEnergyTMEtaResidualTCardCorrNoSelection1Cell(0), fhEnergyTMPhiResidualTCardCorrNoSelection1Cell(0),
121 fhEnergyTMEtaResidualTCardCorrNoSelectionExotic(0), fhEnergyTMPhiResidualTCardCorrNoSelectionExotic(0),
125 fhBadClusterEnergy(0), fhBadClusterTimeEnergy(0), fhBadClusterEtaPhi(0),
126 fhBadClusterPairDiffTimeE(0), fhBadCellTimeSpreadRespectToCellMax(0),
127 fhBadClusterMaxCellCloseCellRatio(0), fhBadClusterMaxCellCloseCellDiff(0), fhBadClusterMaxCellDiff(0),
128 fhBadClusterMaxCellDiffAverageTime(0), fhBadClusterMaxCellDiffWeightedTime(0),
129 fhBadClusterMaxCellECross(0),
130 fhBadClusterLambda0(0), fhBadClusterLambda1(0),
131 fhBadClusterDeltaIEtaDeltaIPhiE0(0), fhBadClusterDeltaIEtaDeltaIPhiE2(0),
132 fhBadClusterDeltaIEtaDeltaIPhiE6(0), fhBadClusterDeltaIA(0),
135 fhRNCells(0), fhXNCells(0),
136 fhYNCells(0), fhZNCells(0),
140 fhRCellE(0), fhXCellE(0),
141 fhYCellE(0), fhZCellE(0),
143 fhDeltaCellClusterRNCells(0), fhDeltaCellClusterXNCells(0),
144 fhDeltaCellClusterYNCells(0), fhDeltaCellClusterZNCells(0),
145 fhDeltaCellClusterRE(0), fhDeltaCellClusterXE(0),
146 fhDeltaCellClusterYE(0), fhDeltaCellClusterZE(0),
149 fhNCells(0), fhNCellsCutAmpMin(0),
150 fhAmplitude(0), fhAmpId(0),
151 fhEtaPhiAmpCell(0), fhEtaPhiCell(0),
153 fhTimeId(0), fhTimeL1UnCorrId(0), fhTimeAmp(0),
154 fhAmpIdLowGain(0), fhTimeIdLowGain(0), fhTimeAmpLowGain(0),
158 fhEMCALPHOSCorrNClusters(0), fhEMCALPHOSCorrEClusters(0),
159 fhEMCALPHOSCorrNCells(0), fhEMCALPHOSCorrECells(0),
160 fhEMCALDCALCorrNClusters(0), fhEMCALDCALCorrEClusters(0),
161 fhEMCALDCALCorrNCells(0), fhEMCALDCALCorrECells(0),
162 fhDCALPHOSCorrNClusters(0), fhDCALPHOSCorrEClusters(0),
163 fhDCALPHOSCorrNCells(0), fhDCALPHOSCorrECells(0),
164 fhCaloV0SCorrNClusters(0), fhCaloV0SCorrEClusters(0),
165 fhCaloV0SCorrNCells(0), fhCaloV0SCorrECells(0),
166 fhCaloV0MCorrNClusters(0), fhCaloV0MCorrEClusters(0),
167 fhCaloV0MCorrNCells(0), fhCaloV0MCorrECells(0),
168 fhCaloTrackMCorrNClusters(0), fhCaloTrackMCorrEClusters(0),
169 fhCaloTrackMCorrNCells(0), fhCaloTrackMCorrECells(0),
170 fhCaloCenNClusters(0), fhCaloCenEClusters(0),
171 fhCaloCenNCells(0), fhCaloCenECells(0),
172 fhCaloEvPNClusters(0), fhCaloEvPEClusters(0),
173 fhCaloEvPNCells(0), fhCaloEvPECells(0),
176 fhEMod(0), fhAmpMod(0),
177 fhEWeirdMod(0), fhAmpWeirdMod(0),
179 fhNClustersMod(0), fhNCellsMod(0),
180 fhNCellsSumAmpPerMod(0), fhNClustersSumEnergyPerMod(0),
181 fhNCellsPerClusterMod(0), fhNCellsPerClusterModNoCut(0),
182 fhNCellsPerClusterWeirdMod(0), fhNCellsPerClusterWeirdModNoCut(0),
184 fhGridCells(0), fhGridCellsE(0), fhGridCellsTime(0),
185 fhGridCellsLowGain(0), fhGridCellsELowGain(0), fhGridCellsTimeLowGain(0),
186 fhTimeAmpPerRCU(0), fhIMMod(0),
189 fhECellClusterRatio(0), fhECellClusterLogRatio(0),
190 fhEMaxCellClusterRatio(0), fhEMaxCellClusterLogRatio(0),
191 fhECellTotalRatio(0), fhECellTotalLogRatio(0),
192 fhECellTotalRatioMod(0), fhECellTotalLogRatioMod(0),
197 fhRecoMCE(), fhRecoMCPhi(), fhRecoMCEta(),
198 fhRecoMCDeltaE(), fhRecoMCRatioE(),
199 fhRecoMCDeltaPhi(), fhRecoMCDeltaEta(),
202 fhGenMCE(), fhGenMCPt(), fhGenMCEtaPhi(),
203 fhGenMCAccE(), fhGenMCAccPt(), fhGenMCAccEtaPhi(),
206 fhEMVxyz(0), fhEMR(0),
207 fhHaVxyz(0), fhHaR(0),
208 fh1EOverP(0), fh2dR(0),
209 fh2EledEdx(0), fh2MatchdEdx(0),
210 fh1EOverPR02(0), fh1EleEOverP(0),
211 fhMCEle1EOverP(0), fhMCEle1dR(0), fhMCEle2MatchdEdx(0),
212 fhMCChHad1EOverP(0), fhMCChHad1dR(0), fhMCChHad2MatchdEdx(0),
213 fhMCNeutral1EOverP(0), fhMCNeutral1dR(0), fhMCNeutral2MatchdEdx(0),
214 fhMCEle1EOverPR02(0), fhMCChHad1EOverPR02(0), fhMCNeutral1EOverPR02(0),
215 fhMCEle1EleEOverP(0), fhMCChHad1EleEOverP(0), fhMCNeutral1EleEOverP(0),
216 fhTrackMatchedDEtaNeg(0), fhTrackMatchedDPhiNeg(0), fhTrackMatchedDEtaDPhiNeg(0),
217 fhTrackMatchedDEtaPos(0), fhTrackMatchedDPhiPos(0), fhTrackMatchedDEtaDPhiPos(0),
218 fhTrackMatchedDEtaNegMod(0), fhTrackMatchedDPhiNegMod(0),
219 fhTrackMatchedDEtaPosMod(0), fhTrackMatchedDPhiPosMod(0),
220 fhClusterTimeEnergyM02(0), fhCellTimeSpreadRespectToCellMaxM02(0),
221 fhClusterMaxCellCloseCellRatioM02(0), fhClusterMaxCellCloseCellDiffM02(0),
222 fhClusterMaxCellDiffM02(0), fhClusterMaxCellECrossM02(0), fhNCellsPerClusterM02(0)
225 for(
Int_t i =0; i < 12; i++)
227 for(
Int_t j = 0; j < 4; j++)
229 for(
Int_t k = 0; k < 3; k++)
236 for(
Int_t j = 0; j < 5; j++)
270 for(
Int_t i = 0; i < 7; i++)
281 for(
Int_t i = 0; i < 4; i++)
291 for(
Int_t i = 0; i < 14; i++)
301 for(
Int_t tm = 0; tm < 2; tm++)
397 for(
Int_t i = 0; i < 6; i++)
399 for(
Int_t j = 0; j < 6; j++)
472 for(
Int_t i = 0; i < 14; i++)
491 for(
Int_t j = 0; j < 6; j++)
493 for(
Int_t k = 0; k < 6; k++)
514 for(
Int_t i = 0; i < 12; i++)
573 if(phi < 0) phi += TMath::TwoPi();
584 for(
Int_t iclus2 = 0; iclus2 < caloClusters->GetEntriesFast(); iclus2++ )
586 AliVCluster* clus2 = (AliVCluster*) caloClusters->At(iclus2);
588 if(clus->GetID() == clus2->GetID())
continue;
593 if(
IsGoodCluster(absIdMax2, clus->GetM02(), clus->GetNCells(), cells) && clus2->GetM02() > 0.1 )
595 Double_t tof2 = clus2->GetTOF()*1.e9;
613 for (
Int_t ipos = 0; ipos < clus->GetNCells(); ipos++)
615 Int_t absId = clus->GetCellsAbsId()[ipos];
616 if(absId!=absIdMax && cells->GetCellAmplitude(absIdMax) > 0.01)
618 Float_t frac = cells->GetCellAmplitude(absId)/cells->GetCellAmplitude(absIdMax);
625 Double_t time = cells->GetCellTime(absId);
643 AliVCaloCells* cells,
651 for (
Int_t ipos = 0; ipos < clus->GetNCells(); ipos++)
653 Int_t id = clus->GetCellsAbsId()[ipos];
656 amp = cells->GetCellAmplitude(
id);
675 Int_t ncells = clus->GetNCells();
677 for (
Int_t ipos = 0; ipos < ncells; ipos++)
679 id = clus ->GetCellsAbsId()[ipos];
680 amp = cells->GetCellAmplitude(
id);
681 time = cells->GetCellTime(
id);
689 wTime += time*1e9 * w;
693 if(ncells > 0) aTime /= ncells;
696 if(wTot > 0) wTime /= wTot;
699 timeAverages[0] = aTime;
700 timeAverages[1] = wTime;
708 Int_t ncells = cells->GetNumberOfCells();
722 nCellsInModule[imod] = 0 ;
723 eCellsInModule[imod] = 0.;
726 Int_t icol = -1, icolAbs = -1;
727 Int_t irow = -1, irowAbs = -1;
737 for (
Int_t iCell = 0; iCell < cells->GetNumberOfCells(); iCell++)
739 if ( cells->GetCellNumber(iCell) < 0 )
continue;
741 AliDebug(2,Form(
"Cell : amp %f, absId %d", cells->GetAmplitude(iCell), cells->GetCellNumber(iCell)));
747 AliDebug(2,Form(
"\t module %d, column %d (%d), row %d (%d)", nModule,icolAbs,icol,irowAbs,irow));
749 if(nModule < fNModules)
756 if(
GetCaloUtils()->GetEMCALChannelStatus(nModule,icol,irow))
continue;
760 if(
GetCaloUtils()->GetPHOSChannelStatus(nModule,icol,irow) )
continue;
764 amp = cells->GetAmplitude(iCell)*recalF;
765 time = cells->GetTime(iCell);
766 id = cells->GetCellNumber(iCell);
767 highG = cells->GetCellHighGain(
id);
781 timeL1UnCorr *= 1.0e9;
786 AliDebug(1,Form(
"Remove cell with Time %f",time));
792 GetCaloUtils()->GetEMCALRecoUtils()->IsExoticCell(
id, cells, bc))
continue;
795 if(amp>=10)binWidthCorrection=1.0/4;
796 if(amp>=20)binWidthCorrection=1.0/10;
825 nCellsInModule[nModule]++ ;
877 if ( cellphi < 0 ) cellphi+=TMath::TwoPi();
891 Float_t rcell = TMath::Sqrt(cellpos[0]*cellpos[0]+cellpos[1]*cellpos[1]);
899 Int_t relId[4], module;
907 Float_t rcell = TMath::Sqrt(xyz.X()*xyz.X()+xyz.Y()*xyz.Y());
927 AliDebug(1,Form(
"Module %d, calo %s, N cells %d, sum Amp %f", imod,
GetCalorimeterString().
Data(), nCellsInModule[imod], eCellsInModule[imod]));
938 for (
Int_t iCell = 0; iCell < cells->GetNumberOfCells(); iCell++)
940 if ( cells->GetCellNumber(iCell) < 0 )
continue;
942 AliDebug(2,Form(
"Cell : amp %f, absId %d", cells->GetAmplitude(iCell), cells->GetCellNumber(iCell)));
946 AliDebug(2,Form(
"\t module %d, column %d, row %d", nModule,icol,irow));
948 if(nModule < fNModules)
955 if(
GetCaloUtils()->GetEMCALChannelStatus(nModule, icol, irow))
continue;
959 if(
GetCaloUtils()->GetPHOSChannelStatus (nModule, icol, irow) )
continue;
963 amp = cells->GetAmplitude(iCell)*recalF;
964 time = cells->GetTime(iCell);
965 id = cells->GetCellNumber(iCell);
972 GetReader()->GetInputEvent()->GetBunchCrossNumber());
980 AliDebug(1,Form(
"Remove cell with Time %f",time));
986 GetCaloUtils()->GetEMCALRecoUtils()->IsExoticCell(
id, cells, bc))
continue;
998 if(eCellsInModule[nModule] > 0)
1000 Float_t ratioMod = amp/eCellsInModule[nModule];
1009 delete [] nCellsInModule;
1010 delete [] eCellsInModule;
1018 Int_t nCaloCellsPerCluster = clus->GetNCells();
1020 UShort_t * indexList = clus->GetCellsAbsId();
1023 clus->GetPosition(pos);
1028 for (
Int_t ipos = 0; ipos < nCaloCellsPerCluster; ipos++)
1031 Int_t absId = indexList[ipos];
1048 Float_t r = TMath::Sqrt(pos[0] *pos[0] + pos[1] * pos[1] );
1049 Float_t rcell = TMath::Sqrt(cellpos[0]*cellpos[0] + cellpos[1]* cellpos[1]);
1057 Int_t relId[4], module;
1073 Float_t r = TMath::Sqrt(pos[0] * pos[0] + pos[1] * pos[1] );
1074 Float_t rcell = TMath::Sqrt(xyz.X() * xyz.X() + xyz.Y() * xyz.Y());
1098 Int_t icol = -1, irow = -1, iRCU = -1, icolAbs = -1, irowAbs = -1;
1108 Int_t ncells = clus->GetNCells();
1109 Double_t time = clus->GetTOF()*1.e9;
1111 Float_t deta = clus->GetTrackDz();
1112 Float_t dphi = clus->GetTrackDx();
1120 if(energy >= 5 && energy < 8)
1122 else if(energy >= 8)
1135 else if(exoticity > 0.97)
1137 if(energy >= 5 && energy < 8)
1139 else if(energy >= 8)
1152 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1154 Int_t absId = clus->GetCellsAbsId()[ipos];
1156 Float_t eCell = cells->GetCellAmplitude(absId);
1162 if( absId == absIdMax || weight < 0.01 )
continue;
1164 Int_t rowDiff = -100, colDiff = -100;
1168 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
1173 if(energy >= 5 && energy < 8)
1175 else if(energy >= 8)
1180 if(energy >= 5 && energy < 8)
1182 else if(energy >= 8)
1188 else if ( energy > 8 )
1201 else if ( energy > 8 && ncells > 1 && exoticity < 0.97)
1211 if ( clus->GetDistanceToBadChannel() < 5 )
return ;
1214 Int_t etaRegion = -1, phiRegion = -1;
1217 if ( etaRegion !=0 ) return ;
1227 else if(exoticity > 0.97)
1237 Int_t absIdList[ncells];
1245 Int_t nCellWithWeight = 1;
1253 Float_t eCellMax = cells->GetCellAmplitude(absIdMax);
1254 Double_t tCellMax = cells->GetCellTime(absIdMax);
1272 Bool_t sameTCard2ndMax = kFALSE;
1276 Float_t eCellSameRowSameTCardNearCol = 0.;
1277 Float_t eCellSameRowDiffTCardNearCol = 0.;
1278 Double_t tCellSameRowSameTCardNearCol = 0.;
1279 Double_t tCellSameRowDiffTCardNearCol = 0.;
1286 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1288 Int_t absId = clus->GetCellsAbsId()[ipos];
1290 Float_t eCell = cells->GetCellAmplitude(absId);
1291 Double_t tCell = cells->GetCellTime(absId);
1301 if( absId == absIdMax || weight < 0.01 )
continue;
1322 Int_t rowDiff = -100, colDiff = -100;
1353 Int_t indexType = -1;
1358 if(TMath::Abs(rowDiff) == 1) nearRow = kTRUE;
1359 if(TMath::Abs(colDiff) == 1) nearCol = kTRUE;
1361 if ( rowDiff == 0 && colDiff != 0 )
1363 if ( nearCol ) indexType = 6;
1369 else if ( rowDiff != 0 && colDiff == 0 )
1371 if ( nearRow ) indexType = 8;
1379 if ( nearRow && nearCol) indexType = 10;
1380 else indexType = 11;
1390 if ( rowDiff == 0 && colDiff != 0 )
1392 if ( nearCol ) indexType = 0;
1396 else if ( rowDiff != 0 && colDiff == 0 )
1398 if ( nearRow ) indexType = 2;
1403 if ( nearCol && nearRow ) indexType = 4;
1408 if ( rowDiff == 0 && TMath::Abs(colDiff) == 1 )
1412 eCellSameRowSameTCardNearCol = eCell;
1413 tCellSameRowSameTCardNearCol = tCell;
1417 eCellSameRowDiffTCardNearCol = eCell;
1418 tCellSameRowDiffTCardNearCol = tCell;
1424 Float_t eCellDiff = eCellMax - eCell;
1425 Float_t eClusDiff = energy - eCell;
1428 Float_t tCellDiff = tCellMax - tCell;
1449 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
1454 if(energy >= 5 && energy < 8)
1456 else if(energy >= 8)
1461 if(energy >= 5 && energy < 8)
1463 else if(energy >= 8)
1468 if ( eCell > emax2nd )
1474 sameTCard2ndMax = kTRUE;
1480 sameTCard2ndMax = kFALSE;
1515 Float_t ratioNcells = nCellWithWeight/(ncells*1.);
1520 if ( nCellWithWeight <= 1 )
return;
1528 if(energy >= 5 && energy < 8)
1530 else if(energy >= 8)
1535 if(energy >= 5 && energy < 8)
1537 else if(energy >= 8)
1600 if(eCellSameRowSameTCardNearCol > 0 && eCellSameRowDiffTCardNearCol > 0)
1602 Float_t eDiff = eCellSameRowSameTCardNearCol - eCellSameRowDiffTCardNearCol ;
1603 Float_t tDiff = tCellSameRowSameTCardNearCol - tCellSameRowDiffTCardNearCol ;
1620 Int_t nCorrInd = nCorr;
1621 if(nCorr > 4) nCorrInd = 5;
1623 Int_t nCorrNoInd = nCorrNo;
1624 if(nCorrNoInd > 4) nCorrNoInd = 5;
1674 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1676 Int_t absId = clus->GetCellsAbsId()[ipos];
1678 Float_t eCell = cells->GetCellAmplitude(absId);
1679 Double_t tCell = cells->GetCellTime(absId);
1689 if( absId == absIdMax || weight < 0.01 )
continue;
1691 Float_t tDiffMaxSecondary = tCellMax - tCell;
1698 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
1701 Int_t rowDiff = -100, colDiff = -100;
1706 if ( energy >= 5 && energy < 8 && nCorr == 0 )
1708 else if ( energy >= 8 && nCorr == 0 )
1713 if ( energy >= 5 && energy < 8 && nCorrNo == 0 )
1715 else if ( energy >= 8 && nCorrNo == 0 )
1750 Double_t tof2 = clus2->GetTOF()*1.e9;
1758 if( absIdMax == absIdMax2
1759 || !
IsGoodCluster(absIdMax2, clus2->GetM02(), clus2->GetNCells(), cells)
1791 if ( energy >= 5 && energy < 8)
1793 else if ( energy >= 8 )
1811 if ( energy >= 5 && energy < 8)
1813 else if ( energy >= 8 )
1879 if ( exoticity > 0.97 )
1894 Float_t nCellRat = nCorr*1. / ((nCorr+nCorrNo)*1.);
1985 if( clus->GetNCells() <= 3 )
return;
1993 for (
Int_t ipos = 0; ipos < clus->GetNCells(); ipos++)
1995 Int_t absId = clus->GetCellsAbsId()[ipos];
2000 if(dIphi < TMath::Abs(iphi-iphiMax)) dIphi = TMath::Abs(iphi-iphiMax);
2004 if(dIeta < TMath::Abs(ieta-ietaMax)) dIeta = TMath::Abs(ieta-ietaMax);
2008 Int_t ietaShift = ieta;
2009 Int_t ietaMaxShift = ietaMax;
2010 if (ieta > ietaMax) ietaMaxShift+=48;
2011 else ietaShift +=48;
2012 if(dIeta < TMath::Abs(ietaShift-ietaMaxShift)) dIeta = TMath::Abs(ietaShift-ietaMaxShift);
2016 Float_t dIA = 1.*(dIphi-dIeta)/(dIeta+dIphi);
2037 Int_t nLabel = clus->GetNLabels();
2038 Int_t* labels = clus->GetLabels();
2087 Double_t tof = clus->GetTOF()*1.e9;
2110 Int_t nCaloCellsPerCluster = clus->GetNCells();
2113 for(
Int_t iclus2 = 0; iclus2 < caloClusters->GetEntriesFast(); iclus2++ )
2115 AliVCluster* clus2 = (AliVCluster*) caloClusters->At(iclus2);
2117 if( clus->GetID() == clus2->GetID() )
continue;
2119 if( clus->GetM02() > 0.01 && clus2->GetM02() > 0.01 )
2123 Double_t tof2 = clus2->GetTOF()*1.e9;
2128 if ( nModule2 == nModule )
2133 if(nCaloCellsPerCluster > 1)
2147 for (
Int_t ipos = 0; ipos < nCaloCellsPerCluster; ipos++)
2149 Int_t absId = clus->GetCellsAbsId()[ipos];
2150 if( absId == absIdMax || cells->GetCellAmplitude(absIdMax) < 0.01 )
continue;
2152 Float_t frac = cells->GetCellAmplitude(absId)/cells->GetCellAmplitude(absIdMax);
2164 Double_t time = cells->GetCellTime(absId);
2184 AliDebug(1,Form(
"cluster: E %2.3f, pT %2.3f, eta %2.3f, phi %2.3f",e,pt,eta,phi*TMath::RadToDeg()));
2199 Int_t icol = -1, irow = -1, iRCU = -1, icolAbs = -1, irowAbs = -1;
2228 clus->GetPosition(pos);
2240 Float_t rxyz = TMath::Sqrt(pos[0]*pos[0]+pos[1]*pos[1]);
2263 AliVCaloCells* cells)
2266 Int_t *labels = 0x0;
2267 Int_t nCaloClusters = caloClusters->GetEntriesFast() ;
2268 Int_t nCaloClustersAccepted = 0 ;
2269 Int_t nCaloCellsPerCluster = 0 ;
2281 nClustersInModule[imod] = 0;
2282 energyInModule [imod] = 0;
2288 for(
Int_t iclus = 0; iclus < nCaloClusters; iclus++)
2290 AliDebug(1,Form(
"Cluster: %d/%d, data %d",iclus+1,nCaloClusters,
GetReader()->GetDataType()));
2292 AliVCluster* clus = (AliVCluster*) caloClusters->At(iclus);
2299 Double_t tof = clus->GetTOF()*1.e9;
2304 AliDebug(1,Form(
"Remove cluster with TOF %f",tof));
2317 nLabel = clus->GetNLabels();
2318 labels = clus->GetLabels();
2324 nCaloCellsPerCluster = clus->GetNCells();
2330 Double_t tmax = cells->GetCellTime(absIdMax);
2343 if(nModule >=0 && nModule < fNModules)
2352 Float_t ampMax = cells->GetCellAmplitude(absIdMax);
2358 Bool_t goodCluster =
IsGoodCluster(absIdMax, clus->GetM02(), nCaloCellsPerCluster, cells);
2361 if(ampMax > 0.01) eCrossFrac = 1-
GetECross(absIdMax,cells)/ampMax;
2363 AliDebug(1,Form(
"Accept cluster? %d",goodCluster));
2376 maxCellFraction, eCrossFrac, tmax);
2378 nCaloClustersAccepted++;
2387 nClustersInModule[nModule]++;
2389 energyInModule [nModule] += clus->E();
2399 Int_t mcOK = kFALSE;
2401 if(
IsDataMC() && nLabel > 0 && labels)
2413 && nCaloClusters > 1
2414 && nCaloCellsPerCluster > 1
2415 &&
GetCaloUtils()->GetNumberOfLocalMaxima(clus,cells) == 1
2431 AliDebug(1,Form(
"Module %d calo %s clusters %d, sum E %f", imod,
GetCalorimeterString().
Data(), nClustersInModule[imod], energyInModule[imod]));
2439 delete [] nClustersInModule;
2440 delete [] energyInModule;
2450 if(!labels || nLabels<=0)
2452 AliWarning(Form(
"Strange, labels array %p, n labels %d", labels,nLabels));
2456 AliDebug(1,Form(
"Primaries: nlabels %d",nLabels));
2461 if(phi < 0) phi +=TMath::TwoPi();
2463 AliAODMCParticle * aodprimary = 0x0;
2464 TParticle * primary = 0x0;
2467 Int_t label = labels[0];
2471 AliDebug(1,Form(
" *** bad label ***: label %d", label));
2490 AliDebug(1,Form(
"*** large label ***: label %d, n tracks %d", label,
GetMCStack()->GetNtrack()));
2496 pdg0 = TMath::Abs(primary->GetPdgCode());
2498 status = primary->GetStatusCode();
2499 vxMC = primary->Vx();
2500 vyMC = primary->Vy();
2501 iParent = primary->GetFirstMother();
2503 AliDebug(1,
"Cluster most contributing mother:");
2504 AliDebug(1,Form(
"\t Mother label %d, pdg %d, %s, status %d, parent %d",iMother, pdg0, primary->GetName(),status, iParent));
2512 pdg = TMath::Abs(primary->GetPdgCode());
2514 AliDebug(2,
"Converted cluster!. Find before conversion:");
2516 while((pdg == 22 || pdg == 11) && status != 1)
2518 Int_t iMotherOrg = iMother;
2521 status = primary->GetStatusCode();
2522 pdg = TMath::Abs(primary->GetPdgCode());
2523 iParent = primary->GetFirstMother();
2527 if(pdg==111 || pdg==221)
2529 primary =
GetMCStack()->Particle(iMotherOrg);
2539 AliDebug(2,Form(
"\t pdg %d, index %d, %s, status %d",pdg, iMother, primary->GetName(),status));
2542 AliDebug(1,
"Converted Cluster mother before conversion:");
2543 AliDebug(1,Form(
"\t Mother label %d, pdg %d, %s, status %d, parent %d",iMother, pdg, primary->GetName(), status, iParent));
2551 AliDebug(2,
"Overlapped Meson decay!, Find it:");
2553 while(pdg != 111 && pdg != 221)
2558 status = primary->GetStatusCode();
2559 pdg = TMath::Abs(primary->GetPdgCode());
2560 iParent = primary->GetFirstMother();
2562 if( iParent < 0 )
break;
2564 AliDebug(2,Form(
"\t pdg %d, %s, index %d",pdg, primary->GetName(),iMother));
2568 AliWarning(
"Tagged as Overlapped photon but meson not found, why?");
2573 AliDebug(2,Form(
"Overlapped %s decay, label %d",primary->GetName(),iMother));
2576 eMC = primary->Energy();
2578 phiMC = primary->Phi();
2579 etaMC = primary->Eta();
2580 pdg = TMath::Abs(primary->GetPdgCode());
2581 charge = (
Int_t) TDatabasePDG::Instance()->GetParticle(pdg)->Charge();
2584 else if(
GetReader()->ReadAODMCParticles() &&
2590 AliFatal(
"MCParticles not available!");
2594 pdg0 = TMath::Abs(aodprimary->GetPdgCode());
2596 status = aodprimary->IsPrimary();
2597 vxMC = aodprimary->Xv();
2598 vyMC = aodprimary->Yv();
2599 iParent = aodprimary->GetMother();
2601 AliDebug(1,
"Cluster most contributing mother:");
2602 AliDebug(1,Form(
"\t Mother label %d, pdg %d, Primary? %d, Physical Primary? %d, parent %d",
2603 iMother, pdg0, aodprimary->IsPrimary(), aodprimary->IsPhysicalPrimary(), iParent));
2608 AliDebug(2,
"Converted cluster!. Find before conversion:");
2612 pdg = TMath::Abs(aodprimary->GetPdgCode());
2614 while ((pdg == 22 || pdg == 11) && !aodprimary->IsPhysicalPrimary())
2616 Int_t iMotherOrg = iMother;
2619 status = aodprimary->IsPrimary();
2620 iParent = aodprimary->GetMother();
2621 pdg = TMath::Abs(aodprimary->GetPdgCode());
2625 if( pdg == 111 || pdg == 221 )
2637 AliDebug(2,Form(
"\t pdg %d, index %d, Primary? %d, Physical Primary? %d",
2638 pdg, iMother, aodprimary->IsPrimary(), aodprimary->IsPhysicalPrimary()));
2641 AliDebug(1,
"Converted Cluster mother before conversion:");
2642 AliDebug(1,Form(
"\t Mother label %d, pdg %d, parent %d, Primary? %d, Physical Primary? %d",
2643 iMother, pdg, iParent, aodprimary->IsPrimary(), aodprimary->IsPhysicalPrimary()));
2651 AliDebug(2,Form(
"Overlapped Meson decay!, Find it: PDG %d, mom %d",pdg, iMother));
2653 while(pdg != 111 && pdg != 221)
2657 status = aodprimary->IsPrimary();
2658 iParent = aodprimary->GetMother();
2659 pdg = TMath::Abs(aodprimary->GetPdgCode());
2661 if( iParent < 0 )
break;
2663 AliDebug(2,Form(
"\t pdg %d, index %d",pdg, iMother));
2667 AliWarning(
"Tagged as Overlapped photon but meson not found, why?");
2672 AliDebug(2,Form(
"Overlapped %s decay, label %d",aodprimary->GetName(),iMother));
2675 status = aodprimary->IsPrimary();
2676 eMC = aodprimary->E();
2678 phiMC = aodprimary->Phi();
2679 etaMC = aodprimary->Eta();
2680 pdg = TMath::Abs(aodprimary->GetPdgCode());
2681 charge = aodprimary->Charge();
2685 Float_t rVMC = TMath::Sqrt(vxMC*vxMC + vyMC*vyMC);
2686 if( ( pdg == 22 || TMath::Abs(pdg) == 11 ) && status != 1 )
2719 else if(charge == 0)
2734 if( mcIndex >= 0 && mcIndex < 7 && e > 0.5 && eMC > 0.5)
2745 if( primary || aodprimary )
return kTRUE ;
2761 if(phi < 0) phi +=TMath::TwoPi();
2779 Double_t dedx = track->GetTPCsignal();
2780 Int_t nITS = track->GetNcls(0);
2781 Int_t nTPC = track->GetNcls(1);
2782 Bool_t positive = kFALSE;
2783 if(track) positive = (track->Charge()>0);
2786 Float_t deta = clus->GetTrackDz();
2787 Float_t dphi = clus->GetTrackDx();
2788 Double_t dR = TMath::Sqrt(dphi*dphi + deta*deta);
2791 if( TMath::Abs(dphi) < 999 )
2828 if(dedx > 60 && dedx < 100)
2838 if(e > 0.5 && tmom > 0.5)
2844 if(dR < 0.02 && eOverP > 0.6 && eOverP < 1.2
2845 && clus->GetNCells() > 1 && nITS > 3 && nTPC > 20)
2853 Double_t charge = TDatabasePDG::Instance()->GetParticle(pdg)->Charge();
2855 if(TMath::Abs(pdg) == 11)
2879 else if(charge == 0)
2904 if(!caloClustersEMCAL || !caloClustersPHOS)
2906 AliDebug(1,Form(
"PHOS (%p) or EMCAL (%p) clusters array not available, do not correlate",caloClustersPHOS,caloClustersEMCAL));
2914 if(!cellsEMCAL || !cellsPHOS)
2916 AliDebug(1,Form(
"PHOS (%p) or EMCAL (%p) cells array ot available, do not correlate",cellsPHOS,cellsEMCAL));
2925 Float_t sumClusterEnergyEMCAL = 0;
2926 Float_t sumClusterEnergyDCAL = 0;
2927 Float_t sumClusterEnergyPHOS = 0;
2931 AliVCluster* cluster = 0;
2932 for(iclus = 0 ; iclus < caloClustersEMCAL->GetEntriesFast() ; iclus++)
2934 cluster = (AliVCluster*)caloClustersEMCAL->At(iclus);
2935 Float_t energy = cluster->E();
2937 if( energy < 0.5 )
continue;
2939 if(cluster->GetCellsAbsId()[0] < 12288)
2942 sumClusterEnergyEMCAL +=
energy;
2947 sumClusterEnergyDCAL +=
energy;
2951 for(iclus = 0 ; iclus < caloClustersPHOS ->GetEntriesFast(); iclus++)
2953 cluster = (AliVCluster*) caloClustersPHOS->At(iclus);
2955 energy = cluster->E();
2957 if( energy < 0.5 )
continue;
2960 sumClusterEnergyPHOS +=
energy;
2964 Int_t ncellsEMCAL = 0 ;
2965 Int_t ncellsDCAL = 0 ;
2966 Int_t ncellsPHOS = 0;
2968 Float_t sumCellEnergyEMCAL = 0;
2969 Float_t sumCellEnergyDCAL = 0;
2970 Float_t sumCellEnergyPHOS = 0;
2972 for(icell = 0 ; icell < cellsEMCAL->GetNumberOfCells() ; icell++)
2974 Float_t amp = cellsEMCAL->GetAmplitude(icell);
2975 Int_t cellId = cellsEMCAL->GetCellNumber(icell);
2982 sumCellEnergyEMCAL += amp;
2987 sumCellEnergyDCAL += amp;
2991 for(icell = 0 ; icell < cellsPHOS->GetNumberOfCells(); icell++)
2993 Float_t amp = cellsPHOS->GetAmplitude(icell);
2994 Int_t cellId = cellsPHOS->GetCellNumber(icell);
2996 if ( cellId < 0 ) continue ;
3001 sumCellEnergyPHOS += amp;
3027 Int_t ncl = nclPHOS;
3028 Float_t sumClusterEnergy = sumClusterEnergyPHOS;
3029 Int_t ncells = ncellsPHOS;
3030 Float_t sumCellEnergy = sumCellEnergyPHOS;
3034 ncl = nclEMCAL + nclDCAL;
3035 sumClusterEnergy = sumClusterEnergyEMCAL + sumClusterEnergyDCAL;
3036 ncells = ncellsEMCAL + ncellsDCAL;
3037 sumCellEnergy = sumCellEnergyEMCAL + sumCellEnergyDCAL;
3065 AliDebug(1,
"Correlate():");
3066 AliDebug(1,Form(
"\t EMCAL: N cells %d, N clusters %d, summed E cells %f, summed E clusters %f",
3067 ncellsEMCAL,nclEMCAL, sumCellEnergyEMCAL,sumClusterEnergyEMCAL));
3068 AliDebug(1,Form(
"\t DCAL : N cells %d, N clusters %d, summed E cells %f, summed E clusters %f",
3069 ncellsDCAL,nclDCAL, sumCellEnergyDCAL,sumClusterEnergyDCAL));
3070 AliDebug(1,Form(
"\t PHOS : N cells %d, N clusters %d, summed E cells %f, summed E clusters %f",
3071 ncellsPHOS,nclPHOS,sumCellEnergyPHOS,sumClusterEnergyPHOS));
3072 AliDebug(1,Form(
"\t V0 : Signal %d, Multiplicity %d, Track Multiplicity %d", v0S,v0M,trM));
3073 AliDebug(1,Form(
"\t centrality : %f, Event plane angle %f", cen,ep));
3082 const Int_t buffersize = 255;
3083 char onePar[buffersize] ;
3085 snprintf(onePar,buffersize,
"--- AliAnaCalorimeterQA ---:") ;
3112 return new TObjString(parList) ;
3192 outputContainer->SetName(
"QAHistos") ;
3251 fhE =
new TH1F (
"hE",
"#it{E} reconstructed clusters ", nptbins*5,ptmin,ptmax*5);
3252 fhE->SetXTitle(
"#it{E} (GeV)");
3253 outputContainer->Add(
fhE);
3255 fhPt =
new TH1F (
"hPt",
"#it{p}_{T} reconstructed clusters", nptbins,ptmin,ptmax);
3256 fhPt->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3257 outputContainer->Add(
fhPt);
3259 fhPhi =
new TH1F (
"hPhi",
"#varphi reconstructed clusters ",nphibins,phimin,phimax);
3260 fhPhi->SetXTitle(
"#varphi (rad)");
3261 outputContainer->Add(
fhPhi);
3263 fhEta =
new TH1F (
"hEta",
"#eta reconstructed clusters ",netabins,etamin,etamax);
3264 fhEta->SetXTitle(
"#eta ");
3265 outputContainer->Add(
fhEta);
3271 fhEtaPhiE =
new TH3F (
"hEtaPhiE",
"#eta vs #varphi vs energy, reconstructed clusters",
3272 netabins,etamin,etamax,nphibins,phimin,phimax,nptbins,ptmin,ptmax);
3280 fhEtaPhi =
new TH2F (
"hEtaPhi",
"#eta vs #varphi for #it{E} > 0.5 GeV, reconstructed clusters",
3281 netabins,etamin,etamax,nphibins,phimin,phimax);
3283 fhEtaPhi->SetYTitle(
"#varphi (rad)");
3289 nptbins,ptmin,ptmax, ntimebins,timemin,timemax);
3295 nptbins,ptmin,ptmax, tdbins,tdmin,tdmax);
3301 nptbins,ptmin,ptmax, tdbins,tdmin,tdmax);
3312 fhLambda0 =
new TH2F (
"hLambda0",
"shower shape, #lambda^{2}_{0} vs E",
3313 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3314 fhLambda0->SetXTitle(
"#it{E}_{cluster} (GeV)");
3315 fhLambda0->SetYTitle(
"#lambda^{2}_{0}");
3318 fhLambda1 =
new TH2F (
"hLambda1",
"shower shape, #lambda^{2}_{1} vs E",
3319 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3320 fhLambda1->SetXTitle(
"#it{E}_{cluster} (GeV)");
3321 fhLambda1->SetYTitle(
"#lambda^{2}_{1}");
3325 nptbins,ptmin,ptmax,10,0,10);
3326 fhNLocMax->SetXTitle(
"#it{E}_{cluster} (GeV)");
3333 (
"hColRowHighEPosTime",
3334 "column vs row, exo < 0.97, E > 8 GeV, t > 5 ns",
3335 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3341 (
"hColRowHighENegTime",
3342 "column vs row, exo < 0.97, E > 8 GeV, t < -5 ns",
3343 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3349 (
"hColRowHighENulTime",
3350 "column vs row, exo < 0.97, E > 8 GeV, -5 < t < 5 ns",
3351 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3359 (
"hColRowExoticHighE1CellPosTime",
3360 "column vs row, 1 cell, E > 8 GeV, t > 5 ns",
3361 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3367 (
"hColRowExoticHighEPosTime",
3368 "column vs row, exo > 0.97, E > 8 GeV, t > 5 ns",
3369 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3375 (
"hColRowExoticHighE1CellNegTime",
3376 "column vs row, 1 cell, E > 8 GeV, t < -5 ns",
3377 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3383 (
"hColRowExoticHighENegTime",
3384 "column vs row, exo > 0.97, E > 8 GeV, t < -5 ns",
3385 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3391 (
"hColRowExoticHighE1CellNulTime",
3392 "column vs row, 1 cell, E > 8 GeV, -5 < t < 5 ns",
3393 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3399 (
"hColRowExoticHighENulTime",
3400 "column vs row, exo > 0.97, E > 8 GeV, -5 < t < 5 ns",
3401 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3407 TString add[] = {
"",
"TrackMatched"};
3408 for(
Int_t tm = 0; tm < 2; tm++)
3411 (Form(
"hColRowTCardCorrNoSelectionLowE%s",add[tm].
Data()),
3412 Form(
"column vs row, max E cell for TCard correlation selected clusters, 5 < E < 8 GeV %s",add[tm].
Data()),
3413 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3419 (Form(
"hColRowTCardCorrNoSelectionHighE%s",add[tm].
Data()),
3420 Form(
"column vs row, max E cell for TCard correlation selected clusters, E > 8 GeV %s",add[tm].
Data()),
3421 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3429 (Form(
"hNCellsTCardCorrNoSelection%s",add[tm].
Data()),
3430 Form(
"# custer # cells vs #it{E} %s",add[tm].
Data()),
3431 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
3437 (Form(
"hNCellsTCardCorrWithWeightNoSelection%s",add[tm].
Data()),
3438 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
3439 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
3445 (Form(
"hNCellsTCardCorrRatioWithWeightNoSelection%s",add[tm].
Data()),
3446 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
3447 nptbins,ptmin,ptmax, 100,0,1);
3453 (Form(
"hLambda0TCardCorrNoSelection%s",add[tm].
Data()),
3454 Form(
"#lambda^{2}_{0} vs #it{E} %s",add[tm].
Data()),
3455 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3461 (Form(
"hLambda1TCardCorrNoSelection%s",add[tm].
Data()),
3462 Form(
"#lambda^{2}_{1} vs #it{E} %s",add[tm].
Data()),
3463 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3469 (Form(
"hLambda0NLM1TCardCorrNoSelection%s",add[tm].
Data()),
3470 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=1 %s",add[tm].
Data()),
3471 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3477 (Form(
"hLambda1NLM1TCardCorrNoSelection%s",add[tm].
Data()),
3478 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=1 %s",add[tm].
Data()),
3479 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3485 (Form(
"hLambda0NLM2TCardCorrNoSelection%s",add[tm].
Data()),
3486 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=2 %s",add[tm].
Data()),
3487 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3493 (Form(
"hLambda1NLM2TCardCorrNoSelection%s",add[tm].
Data()),
3494 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=2 %s",add[tm].
Data()),
3495 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3502 (Form(
"hLambdaRTCardCorrNoSelection%s",add[tm].
Data()),
3503 Form(
"#lambda^{1}_{0}/#lambda^{2}_{0} vs #it{E} %s",add[tm].
Data()),
3504 nptbins,ptmin,ptmax,110,0,1.1);
3510 (Form(
"hNLocMaxTCardCorrNoSelection%s",add[tm].
Data()),
3511 Form(
"#it{n}_{LM} vs E %s",add[tm].
Data()),
3512 nptbins,ptmin,ptmax,10,0,10);
3518 (Form(
"hEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3519 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
3520 nptbins,ptmin,ptmax,100,0,1);
3526 (Form(
"hEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3527 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3528 nptbins,ptmin,ptmax,100,0,1);
3534 (Form(
"hEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3535 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3536 nptbins,ptmin,ptmax,100,0,1);
3542 (Form(
"hE2ndRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3543 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
3544 nptbins,ptmin,ptmax,100,0,0.5);
3550 (Form(
"hE2ndRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3551 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3552 nptbins,ptmin,ptmax,100,0,0.5);
3558 (Form(
"hE2ndRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3559 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3560 nptbins,ptmin,ptmax,100,0,0.5);
3567 (Form(
"hE2ndEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3568 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
3569 nptbins,ptmin,ptmax,100,0,1);
3575 (Form(
"hE2ndEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3576 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3577 nptbins,ptmin,ptmax,100,0,1);
3583 (Form(
"hE2ndEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3584 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3585 nptbins,ptmin,ptmax,100,0,1);
3593 (Form(
"hE2ndSameRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3594 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()),
3595 nptbins,ptmin,ptmax,100,0,0.5);
3601 (Form(
"hE2ndSameRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3602 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()),
3603 nptbins,ptmin,ptmax,100,0,0.5);
3609 (Form(
"hE2ndSameRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3610 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()),
3611 nptbins,ptmin,ptmax,100,0,0.5);
3618 (Form(
"hE2ndSameEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3619 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()),
3620 nptbins,ptmin,ptmax,100,0,1);
3626 (Form(
"hE2ndSameEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3627 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()),
3628 nptbins,ptmin,ptmax,100,0,1);
3634 (Form(
"hE2ndSameEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3635 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()),
3636 nptbins,ptmin,ptmax,100,0,1);
3644 (Form(
"hECellClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3645 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
3646 nptbins,ptmin,ptmax,100,0,0.5);
3652 (Form(
"hECellClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3653 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3654 nptbins,ptmin,ptmax,100,0,0.5);
3660 (Form(
"hECellClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3661 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3662 nptbins,ptmin,ptmax,100,0,0.5);
3668 (Form(
"hLogECellNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3669 Form(
"log(#it{E}_{cell}) vs E_{cluster}, #it{n}_{LM}=1, w > 0.01 %s",add[tm].
Data()),
3670 nptbins,ptmin,ptmax,150,-3,3);
3676 (Form(
"hLogECellNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3677 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}=2, w > 0.01 %s",add[tm].
Data()),
3678 nptbins,ptmin,ptmax,150,-3,3);
3684 (Form(
"hLogECellNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3685 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}>2, w > 0.01 %s",add[tm].
Data()),
3686 nptbins,ptmin,ptmax,150,-3,3);
3693 (Form(
"hECellWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3694 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
3695 nptbins,ptmin,ptmax,90,0,4.5);
3701 (Form(
"hECellWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3702 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3703 nptbins,ptmin,ptmax,90,0,4.5);
3709 (Form(
"hECellWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3710 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster}))vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3711 nptbins,ptmin,ptmax,90,0,4.5);
3718 (Form(
"hECellSameClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3719 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()),
3720 nptbins,ptmin,ptmax,100,0,0.5);
3726 (Form(
"hECellSameClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3727 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, cell from same T-Card as leading, #it{n}_{LM}=2 %s",add[tm].
Data()),
3728 nptbins,ptmin,ptmax,100,0,0.5);
3734 (Form(
"hECellSameClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3735 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, cell from same T-Card as leading, #it{n}_{LM}>2 %s",add[tm].
Data()),
3736 nptbins,ptmin,ptmax,100,0,0.5);
3742 (Form(
"hLogECellSameNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3743 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()),
3744 nptbins,ptmin,ptmax,150,-3,3);
3750 (Form(
"hLogECellSameNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3751 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()),
3752 nptbins,ptmin,ptmax,150,-3,3);
3758 (Form(
"hLogECellSameNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3759 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()),
3760 nptbins,ptmin,ptmax,150,-3,3);
3767 (Form(
"hECellSameWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3768 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()),
3769 nptbins,ptmin,ptmax,90,0,4.5);
3775 (Form(
"hECellSameWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3776 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()),
3777 nptbins,ptmin,ptmax,90,0,4.5);
3783 (Form(
"hECellSameWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3784 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()),
3785 nptbins,ptmin,ptmax,90,0,4.5);
3791 (Form(
"hExoticTCardCorrNoSelection%s",add[tm].
Data()),
3792 Form(
"exoticity vs #it{E} %s",add[tm].
Data()),
3793 nptbins,ptmin,ptmax,200,-1,1);
3799 (Form(
"hNCellsTCardSameAndDiffFraction%s",add[tm].
Data()),
3800 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, %s",add[tm].
Data()),
3801 nptbins,ptmin,ptmax,100,0,1);
3807 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterE%s",add[tm].
Data()),
3808 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()),
3809 nptbins,ptmin,ptmax,200,-10,10);
3815 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterE%s",add[tm].
Data()),
3816 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()),
3817 nptbins,ptmin,ptmax,200,-100,100);
3823 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxE%s",add[tm].
Data()),
3824 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()),
3825 nptbins,ptmin,ptmax,200,-10,10);
3831 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxE%s",add[tm].
Data()),
3832 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()),
3833 nptbins,ptmin,ptmax,200,-100,100);
3841 (Form(
"hEnergyTime1Cell%s",add[tm].
Data()),
3842 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
3843 nptbins,ptmin,ptmax,300,-150,150);
3849 (Form(
"hEnergyTimeExotic%s",add[tm].
Data()),
3850 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
3851 nptbins,ptmin,ptmax,300,-150,150);
3857 (Form(
"hEnergyTimeTCardCorrNoSelection1Cell%s",add[tm].
Data()),
3858 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
3859 nptbins,ptmin,ptmax,300,-150,150);
3865 (Form(
"hEnergyTimeTCardCorrNoSelectionExotic%s",add[tm].
Data()),
3866 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
3867 nptbins,ptmin,ptmax,300,-150,150);
3873 (Form(
"hColRowExoticLowE1Cell%s",add[tm].
Data()),
3874 Form(
"column vs row, 1 cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3875 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3881 (Form(
"hColRowExoticHighE1Cell%s",add[tm].
Data()),
3882 Form(
"column vs row, 1 cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3883 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3889 (Form(
"hColRowExoticLowE%s",add[tm].
Data()),
3890 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3891 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3897 (Form(
"hColRowExoticHighE%s",add[tm].
Data()),
3898 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3899 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3905 (Form(
"hColRowExotic2ndCellDiffLowE%s",add[tm].
Data()),
3906 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3907 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3913 (Form(
"hColRowExotic2ndCellDiffHighE%s",add[tm].
Data()),
3914 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3915 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3921 (Form(
"hColRowExotic2ndCellSameLowE%s",add[tm].
Data()),
3922 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3923 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3929 (Form(
"hColRowExotic2ndCellSameHighE%s",add[tm].
Data()),
3930 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3931 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3937 (Form(
"hColRowTCardCorrNoSelectionExoticLowE%s",add[tm].
Data()),
3938 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3939 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3945 (Form(
"hColRowTCardCorrNoSelectionExoticHighE%s",add[tm].
Data()),
3946 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3947 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3953 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffLowE%s",add[tm].
Data()),
3954 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3955 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3961 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffHighE%s",add[tm].
Data()),
3962 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3963 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3969 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameLowE%s",add[tm].
Data()),
3970 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3971 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3977 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameHighE%s",add[tm].
Data()),
3978 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3979 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3985 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameLowE%s",add[tm].
Data()),
3986 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3987 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
3993 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameHighE%s",add[tm].
Data()),
3994 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3995 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
4001 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffLowE%s",add[tm].
Data()),
4002 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
4003 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
4009 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffHighE%s",add[tm].
Data()),
4010 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
4011 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
4017 (Form(
"hNCellsTCardSameAndDiffFraction_Exotic%s",add[tm].
Data()),
4018 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, exo > 0.97 %s",add[tm].
Data()),
4019 nptbins,ptmin,ptmax,100,0,1);
4025 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterEExo%s",add[tm].
Data()),
4026 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()),
4027 nptbins,ptmin,ptmax,200,-10,10);
4033 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterEExo%s",add[tm].
Data()),
4034 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()),
4035 nptbins,ptmin,ptmax,200,-100,100);
4041 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxEExo%s",add[tm].
Data()),
4042 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()),
4043 nptbins,ptmin,ptmax,200,-10,10);
4049 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxEExo%s",add[tm].
Data()),
4050 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()),
4051 nptbins,ptmin,ptmax,200,-100,100);
4057 for(
Int_t i = 0; i < 6; i++)
4059 for(
Int_t j = 0; j < 6; j++)
4062 (Form(
"hLambda0TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4063 Form(
"#lambda^{2}_{0} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4064 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4070 (Form(
"hLambda1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4071 Form(
"#lambda^{2}_{1} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4072 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4078 (Form(
"hLambda0NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4079 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()),
4080 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4086 (Form(
"hLambda1NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4087 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()),
4088 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4094 (Form(
"hLambda0NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4095 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()),
4096 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4102 (Form(
"hLambda1NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4103 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()),
4104 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4119 (Form(
"hNLocMaxTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4120 Form(
"#it{n}_{LM} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4121 nptbins,ptmin,ptmax,10,0,10);
4127 (Form(
"hEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4128 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()),
4129 nptbins,ptmin,ptmax,100,0,1);
4135 (Form(
"hEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4136 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()),
4137 nptbins,ptmin,ptmax,100,0,1);
4143 (Form(
"hEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4144 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()),
4145 nptbins,ptmin,ptmax,100,0,1);
4151 (Form(
"hE2ndRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4152 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()),
4153 nptbins,ptmin,ptmax,100,0,0.5);
4159 (Form(
"hE2ndRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4160 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()),
4161 nptbins,ptmin,ptmax,100,0,0.5);
4167 (Form(
"hE2ndRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4168 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()),
4169 nptbins,ptmin,ptmax,100,0,0.5);
4175 (Form(
"hE2ndEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4176 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()),
4177 nptbins,ptmin,ptmax,100,0,1);
4183 (Form(
"hE2ndEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4184 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()),
4185 nptbins,ptmin,ptmax,100,0,1);
4191 (Form(
"hE2ndEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4192 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()),
4193 nptbins,ptmin,ptmax,100,0,1);
4200 (Form(
"hECellClusRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4201 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()),
4202 nptbins,ptmin,ptmax,100,0,0.5);
4208 (Form(
"hECellClusRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4209 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()),
4210 nptbins,ptmin,ptmax,100,0,0.5);
4216 (Form(
"hECellClusRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4217 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()),
4218 nptbins,ptmin,ptmax,100,0,0.5);
4224 (Form(
"hLogECellNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4225 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()),
4226 nptbins,ptmin,ptmax,150,-3,3);
4232 (Form(
"hLogECellNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4233 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()),
4234 nptbins,ptmin,ptmax,150,-3,3);
4240 (Form(
"hLogECellNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4241 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()),
4242 nptbins,ptmin,ptmax,150,-3,3);
4248 (Form(
"hECellWeightNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4249 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()),
4250 nptbins,ptmin,ptmax,90,0,4.5);
4256 (Form(
"hECellWeightNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4257 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()),
4258 nptbins,ptmin,ptmax,90,0,4.5);
4264 (Form(
"hECellWeightNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4265 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()),
4266 nptbins,ptmin,ptmax,90,0,4.5);
4272 (Form(
"hMassEClusTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4273 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()),
4274 nptbins,ptmin,ptmax,nmassbins,massmin,massmax);
4288 (Form(
"hTimeDiffTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4289 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()),
4290 nptbins,ptmin,ptmax,300,-150,150);
4296 (Form(
"hColRowTCardCorrelNCellLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4297 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4298 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
4304 (Form(
"hColRowTCardCorrelNCellHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4305 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4306 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
4314 (Form(
"hExoticTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4315 Form(
"exoticity vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4316 nptbins,ptmin,ptmax,200,-1,1);
4322 (Form(
"hTimeDiffExoTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4323 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()),
4324 nptbins,ptmin,ptmax,300,-150,150);
4330 (Form(
"hColRowTCardCorrelNCellExoticLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4331 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4332 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
4338 (Form(
"hColRowTCardCorrelNCellExoticHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4339 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4340 96+2,-1.5,96+0.5,(8*24+2*8)+2,-1.5,(8*24+2*8)+0.5);
4555 (Form(
"hLambda0Exoticity_EBin%d%s",i,add[tm].
Data()),
4556 Form(
"#lambda^{2}_{0} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
4557 200,-1,1,ssbins,ssmin,ssmax);
4563 (Form(
"hLambda1Exoticity_EBin%d%s",i,add[tm].
Data()),
4564 Form(
"#lambda^{2}_{1} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
4565 200,-1,1,ssbins,ssmin,ssmax);
4579 (Form(
"hNCellsExoticity_EBin%d%s",i,add[tm].
Data()),
4580 Form(
"#it{n}_{cells} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
4581 200,-1,1,nceclbins,nceclmin,nceclmax);
4582 fhNCellsExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
4587 (Form(
"hTimeExoticity_EBin%d%s",i,add[tm].
Data()),
4589 200,-1,1,100,-25,25);
4590 fhTimeExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
4595 (Form(
"hNCellsTCardSameAndDiff_Exotic_EBin%d%s",i,add[tm].
Data()),
4596 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()),
4597 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
4637 (Form(
"hNCellsTCardSameAndDiff_EBin%d%s",i,add[tm].
Data()),
4638 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()),
4639 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
4645 (Form(
"hLambda0Lambda1_EBin%d%s",i,add[tm].
Data()),
4646 Form(
"#lambda^{2}_{0} vs #lambda^{2}_{1}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
4647 ssbins,ssmin,ssmax,ssbins,ssmin,ssmax);
4663 for(
Int_t j = 0; j < 6; j++)
4665 for(
Int_t k = 0; k < 6; k++)
4668 (Form(
"hLambda0Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
4669 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()),
4670 200,-1,1,ssbins,ssmin,ssmax);
4676 (Form(
"hLambda1Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
4677 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()),
4678 200,-1,1,ssbins,ssmin,ssmax);
4775 for(
Int_t i = 0; i < 12; i++)
4778 (Form(
"hTCardCorrECellMaxDiff_Case%d%s",i,add[tm].
Data()),
4779 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()),
4780 nptbins,ptmin,ptmax,210,-1,20);
4786 (Form(
"hTCardCorrEClusterDiff_Case%d%s",i,add[tm].
Data()),
4787 Form(
"#it{E}_{cluster}-#it{E}_{cell} vs #it{E}_{cluster}, for (un)correlated cells in TCard, case %d %s",i,add[tm].
Data()),
4788 nptbins,ptmin,ptmax,210,-1,20);
4810 (Form(
"hTCardCorrTCellMaxDiff_Case%d%s",i,add[tm].
Data()),
4811 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()),
4812 nptbins,ptmin,ptmax,1000,-100,100);
4820 (Form(
"hTCardCorrECellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
4821 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()),
4822 nptbins,ptmin,ptmax,210,-1,20);
4828 (Form(
"hTCardCorrEClusterDiffExo_Case%d%s",i,add[tm].
Data()),
4829 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()),
4830 nptbins,ptmin,ptmax,210,-1,20);
4852 (Form(
"hTCardCorrTCellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
4853 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()),
4854 nptbins,ptmin,ptmax,1000,-100,100);
4865 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
4871 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
4877 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
4883 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
4889 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
4895 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
4901 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
4907 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
4915 (Form(
"hTMPhiResidual_EBin%d",i),
4916 Form(
"#Delta #varphi_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
4917 200,-1,1,nresphibins,resphimin,resphimax);
4923 (Form(
"hTMEtaResidual_EBin%d",i),
4924 Form(
"#Delta #eta_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
4925 200,-1,1,nresetabins,resetamin,resetamax);
4931 (Form(
"hTMPhiResidual_LooseCut_EBin%d",i),
4932 Form(
"#Delta #varphi_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
4933 200,-1,1,nresphibins,resphimin,resphimax);
4939 (Form(
"hTMEtaResidual_LooseCut_EBin%d",i),
4940 Form(
"#Delta #eta_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
4941 200,-1,1,nresetabins,resetamin,resetamax);