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") ;
3263 fhE =
new TH1F (
"hE",
"#it{E} reconstructed clusters ", nptbins*5,ptmin,ptmax*5);
3264 fhE->SetXTitle(
"#it{E} (GeV)");
3265 outputContainer->Add(
fhE);
3267 fhPt =
new TH1F (
"hPt",
"#it{p}_{T} reconstructed clusters", nptbins,ptmin,ptmax);
3268 fhPt->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3269 outputContainer->Add(
fhPt);
3271 fhPhi =
new TH1F (
"hPhi",
"#varphi reconstructed clusters ",nphibins,phimin,phimax);
3272 fhPhi->SetXTitle(
"#varphi (rad)");
3273 outputContainer->Add(
fhPhi);
3275 fhEta =
new TH1F (
"hEta",
"#eta reconstructed clusters ",netabins,etamin,etamax);
3276 fhEta->SetXTitle(
"#eta ");
3277 outputContainer->Add(
fhEta);
3283 fhEtaPhiE =
new TH3F (
"hEtaPhiE",
"#eta vs #varphi vs energy, reconstructed clusters",
3284 netabins,etamin,etamax,nphibins,phimin,phimax,nptbins,ptmin,ptmax);
3292 fhEtaPhi =
new TH2F (
"hEtaPhi",
"#eta vs #varphi for #it{E} > 0.5 GeV, reconstructed clusters",
3293 netabins,etamin,etamax,nphibins,phimin,phimax);
3295 fhEtaPhi->SetYTitle(
"#varphi (rad)");
3301 nptbins,ptmin,ptmax, ntimebins,timemin,timemax);
3307 nptbins,ptmin,ptmax, tdbins,tdmin,tdmax);
3313 nptbins,ptmin,ptmax, tdbins,tdmin,tdmax);
3324 fhLambda0 =
new TH2F (
"hLambda0",
"shower shape, #lambda^{2}_{0} vs E",
3325 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3326 fhLambda0->SetXTitle(
"#it{E}_{cluster} (GeV)");
3327 fhLambda0->SetYTitle(
"#lambda^{2}_{0}");
3330 fhLambda1 =
new TH2F (
"hLambda1",
"shower shape, #lambda^{2}_{1} vs E",
3331 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3332 fhLambda1->SetXTitle(
"#it{E}_{cluster} (GeV)");
3333 fhLambda1->SetYTitle(
"#lambda^{2}_{1}");
3337 nptbins,ptmin,ptmax,10,0,10);
3338 fhNLocMax->SetXTitle(
"#it{E}_{cluster} (GeV)");
3345 (
"hColRowHighEPosTime",
3346 "column vs row, exo < 0.97, E > 8 GeV, t > 5 ns",
3347 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3353 (
"hColRowHighENegTime",
3354 "column vs row, exo < 0.97, E > 8 GeV, t < -5 ns",
3355 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3361 (
"hColRowHighENulTime",
3362 "column vs row, exo < 0.97, E > 8 GeV, -5 < t < 5 ns",
3363 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3371 (
"hColRowExoticHighE1CellPosTime",
3372 "column vs row, 1 cell, E > 8 GeV, t > 5 ns",
3373 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3379 (
"hColRowExoticHighEPosTime",
3380 "column vs row, exo > 0.97, E > 8 GeV, t > 5 ns",
3381 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3387 (
"hColRowExoticHighE1CellNegTime",
3388 "column vs row, 1 cell, E > 8 GeV, t < -5 ns",
3389 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3395 (
"hColRowExoticHighENegTime",
3396 "column vs row, exo > 0.97, E > 8 GeV, t < -5 ns",
3397 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3403 (
"hColRowExoticHighE1CellNulTime",
3404 "column vs row, 1 cell, E > 8 GeV, -5 < t < 5 ns",
3405 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3411 (
"hColRowExoticHighENulTime",
3412 "column vs row, exo > 0.97, E > 8 GeV, -5 < t < 5 ns",
3413 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3419 TString add[] = {
"",
"TrackMatched"};
3420 for(
Int_t tm = 0; tm < 2; tm++)
3423 (Form(
"hColRowTCardCorrNoSelectionLowE%s",add[tm].
Data()),
3424 Form(
"column vs row, max E cell for TCard correlation selected clusters, 5 < E < 8 GeV %s",add[tm].
Data()),
3425 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3431 (Form(
"hColRowTCardCorrNoSelectionHighE%s",add[tm].
Data()),
3432 Form(
"column vs row, max E cell for TCard correlation selected clusters, E > 8 GeV %s",add[tm].
Data()),
3433 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3441 (Form(
"hNCellsTCardCorrNoSelection%s",add[tm].
Data()),
3442 Form(
"# custer # cells vs #it{E} %s",add[tm].
Data()),
3443 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
3449 (Form(
"hNCellsTCardCorrWithWeightNoSelection%s",add[tm].
Data()),
3450 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
3451 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
3457 (Form(
"hNCellsTCardCorrRatioWithWeightNoSelection%s",add[tm].
Data()),
3458 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
3459 nptbins,ptmin,ptmax, 100,0,1);
3465 (Form(
"hLambda0TCardCorrNoSelection%s",add[tm].
Data()),
3466 Form(
"#lambda^{2}_{0} vs #it{E} %s",add[tm].
Data()),
3467 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3473 (Form(
"hLambda1TCardCorrNoSelection%s",add[tm].
Data()),
3474 Form(
"#lambda^{2}_{1} vs #it{E} %s",add[tm].
Data()),
3475 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3481 (Form(
"hLambda0NLM1TCardCorrNoSelection%s",add[tm].
Data()),
3482 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=1 %s",add[tm].
Data()),
3483 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3489 (Form(
"hLambda1NLM1TCardCorrNoSelection%s",add[tm].
Data()),
3490 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=1 %s",add[tm].
Data()),
3491 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3497 (Form(
"hLambda0NLM2TCardCorrNoSelection%s",add[tm].
Data()),
3498 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=2 %s",add[tm].
Data()),
3499 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3505 (Form(
"hLambda1NLM2TCardCorrNoSelection%s",add[tm].
Data()),
3506 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=2 %s",add[tm].
Data()),
3507 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3514 (Form(
"hLambdaRTCardCorrNoSelection%s",add[tm].
Data()),
3515 Form(
"#lambda^{1}_{0}/#lambda^{2}_{0} vs #it{E} %s",add[tm].
Data()),
3516 nptbins,ptmin,ptmax,110,0,1.1);
3522 (Form(
"hNLocMaxTCardCorrNoSelection%s",add[tm].
Data()),
3523 Form(
"#it{n}_{LM} vs E %s",add[tm].
Data()),
3524 nptbins,ptmin,ptmax,10,0,10);
3530 (Form(
"hEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3531 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
3532 nptbins,ptmin,ptmax,100,0,1);
3538 (Form(
"hEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3539 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3540 nptbins,ptmin,ptmax,100,0,1);
3546 (Form(
"hEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3547 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3548 nptbins,ptmin,ptmax,100,0,1);
3554 (Form(
"hE2ndRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3555 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
3556 nptbins,ptmin,ptmax,100,0,1);
3562 (Form(
"hE2ndRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3563 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3564 nptbins,ptmin,ptmax,100,0,1);
3570 (Form(
"hE2ndRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3571 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3572 nptbins,ptmin,ptmax,100,0,1);
3579 (Form(
"hE2ndEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3580 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
3581 nptbins,ptmin,ptmax,100,0,1);
3587 (Form(
"hE2ndEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3588 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3589 nptbins,ptmin,ptmax,100,0,1);
3595 (Form(
"hE2ndEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3596 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3597 nptbins,ptmin,ptmax,100,0,1);
3605 (Form(
"hE2ndSameRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3606 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()),
3607 nptbins,ptmin,ptmax,100,0,1);
3613 (Form(
"hE2ndSameRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3614 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()),
3615 nptbins,ptmin,ptmax,100,0,1);
3621 (Form(
"hE2ndSameRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3622 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()),
3623 nptbins,ptmin,ptmax,100,0,1);
3630 (Form(
"hE2ndSameEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3631 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()),
3632 nptbins,ptmin,ptmax,100,0,1);
3638 (Form(
"hE2ndSameEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3639 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()),
3640 nptbins,ptmin,ptmax,100,0,1);
3646 (Form(
"hE2ndSameEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3647 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()),
3648 nptbins,ptmin,ptmax,100,0,1);
3656 (Form(
"hECellClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3657 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
3658 nptbins,ptmin,ptmax,100,0,1);
3664 (Form(
"hECellClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3665 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3666 nptbins,ptmin,ptmax,100,0,1);
3672 (Form(
"hECellClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3673 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3674 nptbins,ptmin,ptmax,100,0,1);
3680 (Form(
"hLogECellNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3681 Form(
"log(#it{E}_{cell}) vs E_{cluster}, #it{n}_{LM}=1, w > 0.01 %s",add[tm].
Data()),
3682 nptbins,ptmin,ptmax,150,-3,3);
3688 (Form(
"hLogECellNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3689 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}=2, w > 0.01 %s",add[tm].
Data()),
3690 nptbins,ptmin,ptmax,150,-3,3);
3696 (Form(
"hLogECellNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3697 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}>2, w > 0.01 %s",add[tm].
Data()),
3698 nptbins,ptmin,ptmax,150,-3,3);
3705 (Form(
"hECellWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3706 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
3707 nptbins,ptmin,ptmax,90,0,4.5);
3713 (Form(
"hECellWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3714 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
3715 nptbins,ptmin,ptmax,90,0,4.5);
3721 (Form(
"hECellWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3722 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster}))vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
3723 nptbins,ptmin,ptmax,90,0,4.5);
3730 (Form(
"hECellSameClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3731 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()),
3732 nptbins,ptmin,ptmax,100,0,1);
3738 (Form(
"hECellSameClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3739 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, cell from same T-Card as leading, #it{n}_{LM}=2 %s",add[tm].
Data()),
3740 nptbins,ptmin,ptmax,100,0,1);
3746 (Form(
"hECellSameClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3747 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, cell from same T-Card as leading, #it{n}_{LM}>2 %s",add[tm].
Data()),
3748 nptbins,ptmin,ptmax,100,0,1);
3754 (Form(
"hLogECellSameNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3755 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()),
3756 nptbins,ptmin,ptmax,150,-3,3);
3762 (Form(
"hLogECellSameNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3763 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()),
3764 nptbins,ptmin,ptmax,150,-3,3);
3770 (Form(
"hLogECellSameNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3771 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()),
3772 nptbins,ptmin,ptmax,150,-3,3);
3779 (Form(
"hECellSameWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
3780 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()),
3781 nptbins,ptmin,ptmax,90,0,4.5);
3787 (Form(
"hECellSameWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
3788 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()),
3789 nptbins,ptmin,ptmax,90,0,4.5);
3795 (Form(
"hECellSameWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
3796 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()),
3797 nptbins,ptmin,ptmax,90,0,4.5);
3803 (Form(
"hExoticTCardCorrNoSelection%s",add[tm].
Data()),
3804 Form(
"exoticity vs #it{E} %s",add[tm].
Data()),
3805 nptbins,ptmin,ptmax,200,-1,1);
3811 (Form(
"hNCellsTCardSameAndDiffFraction%s",add[tm].
Data()),
3812 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, %s",add[tm].
Data()),
3813 nptbins,ptmin,ptmax,100,0,1);
3819 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterE%s",add[tm].
Data()),
3820 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()),
3821 nptbins,ptmin,ptmax,200,-10,10);
3827 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterE%s",add[tm].
Data()),
3828 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()),
3829 nptbins,ptmin,ptmax,200,-100,100);
3835 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxE%s",add[tm].
Data()),
3836 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()),
3837 nptbins,ptmin,ptmax,200,-10,10);
3843 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxE%s",add[tm].
Data()),
3844 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()),
3845 nptbins,ptmin,ptmax,200,-100,100);
3853 (Form(
"hEnergyTime1Cell%s",add[tm].
Data()),
3854 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
3855 nptbins,ptmin,ptmax,300,-150,150);
3861 (Form(
"hEnergyTimeExotic%s",add[tm].
Data()),
3862 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
3863 nptbins,ptmin,ptmax,300,-150,150);
3869 (Form(
"hEnergyTimeTCardCorrNoSelection1Cell%s",add[tm].
Data()),
3870 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
3871 nptbins,ptmin,ptmax,300,-150,150);
3877 (Form(
"hEnergyTimeTCardCorrNoSelectionExotic%s",add[tm].
Data()),
3878 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
3879 nptbins,ptmin,ptmax,300,-150,150);
3885 (Form(
"hColRowExoticLowE1Cell%s",add[tm].
Data()),
3886 Form(
"column vs row, 1 cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3887 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3893 (Form(
"hColRowExoticHighE1Cell%s",add[tm].
Data()),
3894 Form(
"column vs row, 1 cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3895 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3901 (Form(
"hColRowExoticLowE%s",add[tm].
Data()),
3902 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3903 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3909 (Form(
"hColRowExoticHighE%s",add[tm].
Data()),
3910 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3911 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3917 (Form(
"hColRowExotic2ndCellDiffLowE%s",add[tm].
Data()),
3918 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3919 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3925 (Form(
"hColRowExotic2ndCellDiffHighE%s",add[tm].
Data()),
3926 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3927 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3933 (Form(
"hColRowExotic2ndCellSameLowE%s",add[tm].
Data()),
3934 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3935 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3941 (Form(
"hColRowExotic2ndCellSameHighE%s",add[tm].
Data()),
3942 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3943 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3949 (Form(
"hColRowTCardCorrNoSelectionExoticLowE%s",add[tm].
Data()),
3950 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3951 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3957 (Form(
"hColRowTCardCorrNoSelectionExoticHighE%s",add[tm].
Data()),
3958 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3959 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3965 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffLowE%s",add[tm].
Data()),
3966 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3967 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3973 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffHighE%s",add[tm].
Data()),
3974 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3975 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3981 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameLowE%s",add[tm].
Data()),
3982 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3983 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3989 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameHighE%s",add[tm].
Data()),
3990 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3991 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
3997 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameLowE%s",add[tm].
Data()),
3998 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3999 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4005 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameHighE%s",add[tm].
Data()),
4006 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
4007 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4013 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffLowE%s",add[tm].
Data()),
4014 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
4015 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4021 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffHighE%s",add[tm].
Data()),
4022 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
4023 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4029 (Form(
"hNCellsTCardSameAndDiffFraction_Exotic%s",add[tm].
Data()),
4030 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, exo > 0.97 %s",add[tm].
Data()),
4031 nptbins,ptmin,ptmax,100,0,1);
4037 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterEExo%s",add[tm].
Data()),
4038 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()),
4039 nptbins,ptmin,ptmax,200,-10,10);
4045 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterEExo%s",add[tm].
Data()),
4046 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()),
4047 nptbins,ptmin,ptmax,200,-100,100);
4053 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxEExo%s",add[tm].
Data()),
4054 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()),
4055 nptbins,ptmin,ptmax,200,-10,10);
4061 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxEExo%s",add[tm].
Data()),
4062 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()),
4063 nptbins,ptmin,ptmax,200,-100,100);
4069 for(
Int_t i = 0; i < 6; i++)
4071 for(
Int_t j = 0; j < 6; j++)
4074 (Form(
"hLambda0TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4075 Form(
"#lambda^{2}_{0} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4076 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4082 (Form(
"hLambda1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4083 Form(
"#lambda^{2}_{1} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4084 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4090 (Form(
"hLambda0NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4091 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()),
4092 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4098 (Form(
"hLambda1NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4099 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()),
4100 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4106 (Form(
"hLambda0NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4107 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()),
4108 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4114 (Form(
"hLambda1NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4115 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()),
4116 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4131 (Form(
"hNLocMaxTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4132 Form(
"#it{n}_{LM} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4133 nptbins,ptmin,ptmax,10,0,10);
4139 (Form(
"hEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4140 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()),
4141 nptbins,ptmin,ptmax,100,0,1);
4147 (Form(
"hEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4148 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()),
4149 nptbins,ptmin,ptmax,100,0,1);
4155 (Form(
"hEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4156 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()),
4157 nptbins,ptmin,ptmax,100,0,1);
4163 (Form(
"hE2ndRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4164 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()),
4165 nptbins,ptmin,ptmax,100,0,1);
4171 (Form(
"hE2ndRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4172 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()),
4173 nptbins,ptmin,ptmax,100,0,1);
4179 (Form(
"hE2ndRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4180 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()),
4181 nptbins,ptmin,ptmax,100,0,1);
4187 (Form(
"hE2ndEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4188 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()),
4189 nptbins,ptmin,ptmax,100,0,1);
4195 (Form(
"hE2ndEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4196 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()),
4197 nptbins,ptmin,ptmax,100,0,1);
4203 (Form(
"hE2ndEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4204 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()),
4205 nptbins,ptmin,ptmax,100,0,1);
4212 (Form(
"hECellClusRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4213 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()),
4214 nptbins,ptmin,ptmax,100,0,1);
4220 (Form(
"hECellClusRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4221 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()),
4222 nptbins,ptmin,ptmax,100,0,1);
4228 (Form(
"hECellClusRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4229 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()),
4230 nptbins,ptmin,ptmax,100,0,1);
4236 (Form(
"hLogECellNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4237 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()),
4238 nptbins,ptmin,ptmax,150,-3,3);
4244 (Form(
"hLogECellNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4245 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()),
4246 nptbins,ptmin,ptmax,150,-3,3);
4252 (Form(
"hLogECellNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4253 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()),
4254 nptbins,ptmin,ptmax,150,-3,3);
4260 (Form(
"hECellWeightNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4261 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()),
4262 nptbins,ptmin,ptmax,90,0,4.5);
4268 (Form(
"hECellWeightNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4269 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()),
4270 nptbins,ptmin,ptmax,90,0,4.5);
4276 (Form(
"hECellWeightNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4277 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()),
4278 nptbins,ptmin,ptmax,90,0,4.5);
4284 (Form(
"hMassEClusTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4285 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()),
4286 nptbins,ptmin,ptmax,nmassbins,massmin,massmax);
4300 (Form(
"hTimeDiffTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4301 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()),
4302 nptbins,ptmin,ptmax,300,-150,150);
4308 (Form(
"hColRowTCardCorrelNCellLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4309 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4310 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4316 (Form(
"hColRowTCardCorrelNCellHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4317 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4318 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4326 (Form(
"hExoticTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4327 Form(
"exoticity vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4328 nptbins,ptmin,ptmax,200,-1,1);
4334 (Form(
"hTimeDiffExoTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4335 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()),
4336 nptbins,ptmin,ptmax,300,-150,150);
4342 (Form(
"hColRowTCardCorrelNCellExoticLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4343 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4344 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4350 (Form(
"hColRowTCardCorrelNCellExoticHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
4351 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
4352 colmaxs+2,-1.5,colmaxs+0.5, rowmaxs+2,-1.5,rowmaxs+0.5);
4567 (Form(
"hLambda0Exoticity_EBin%d%s",i,add[tm].
Data()),
4568 Form(
"#lambda^{2}_{0} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
4569 200,-1,1,ssbins,ssmin,ssmax);
4575 (Form(
"hLambda1Exoticity_EBin%d%s",i,add[tm].
Data()),
4576 Form(
"#lambda^{2}_{1} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
4577 200,-1,1,ssbins,ssmin,ssmax);
4591 (Form(
"hNCellsExoticity_EBin%d%s",i,add[tm].
Data()),
4592 Form(
"#it{n}_{cells} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
4593 200,-1,1,nceclbins,nceclmin,nceclmax);
4594 fhNCellsExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
4599 (Form(
"hTimeExoticity_EBin%d%s",i,add[tm].
Data()),
4601 200,-1,1,100,-25,25);
4602 fhTimeExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
4607 (Form(
"hNCellsTCardSameAndDiff_Exotic_EBin%d%s",i,add[tm].
Data()),
4608 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()),
4609 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
4649 (Form(
"hNCellsTCardSameAndDiff_EBin%d%s",i,add[tm].
Data()),
4650 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()),
4651 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
4657 (Form(
"hLambda0Lambda1_EBin%d%s",i,add[tm].
Data()),
4658 Form(
"#lambda^{2}_{0} vs #lambda^{2}_{1}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
4659 ssbins,ssmin,ssmax,ssbins,ssmin,ssmax);
4675 for(
Int_t j = 0; j < 6; j++)
4677 for(
Int_t k = 0; k < 6; k++)
4680 (Form(
"hLambda0Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
4681 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()),
4682 200,-1,1,ssbins,ssmin,ssmax);
4688 (Form(
"hLambda1Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
4689 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()),
4690 200,-1,1,ssbins,ssmin,ssmax);
4787 for(
Int_t i = 0; i < 12; i++)
4790 (Form(
"hTCardCorrECellMaxDiff_Case%d%s",i,add[tm].
Data()),
4791 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()),
4792 nptbins,ptmin,ptmax,210,-1,20);
4798 (Form(
"hTCardCorrEClusterDiff_Case%d%s",i,add[tm].
Data()),
4799 Form(
"#it{E}_{cluster}-#it{E}_{cell} vs #it{E}_{cluster}, for (un)correlated cells in TCard, case %d %s",i,add[tm].
Data()),
4800 nptbins,ptmin,ptmax,210,-1,20);
4822 (Form(
"hTCardCorrTCellMaxDiff_Case%d%s",i,add[tm].
Data()),
4823 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()),
4824 nptbins,ptmin,ptmax,1000,-100,100);
4832 (Form(
"hTCardCorrECellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
4833 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()),
4834 nptbins,ptmin,ptmax,210,-1,20);
4840 (Form(
"hTCardCorrEClusterDiffExo_Case%d%s",i,add[tm].
Data()),
4841 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()),
4842 nptbins,ptmin,ptmax,210,-1,20);
4864 (Form(
"hTCardCorrTCellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
4865 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()),
4866 nptbins,ptmin,ptmax,1000,-100,100);
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);
4913 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
4919 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
4927 (Form(
"hTMPhiResidual_EBin%d",i),
4928 Form(
"#Delta #varphi_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
4929 200,-1,1,nresphibins,resphimin,resphimax);
4935 (Form(
"hTMEtaResidual_EBin%d",i),
4936 Form(
"#Delta #eta_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
4937 200,-1,1,nresetabins,resetamin,resetamax);
4943 (Form(
"hTMPhiResidual_LooseCut_EBin%d",i),
4944 Form(
"#Delta #varphi_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
4945 200,-1,1,nresphibins,resphimin,resphimax);