17 #include <TObjArray.h> 19 #include <TObjString.h> 24 #include "AliVCaloCells.h" 26 #include "AliVCluster.h" 27 #include "AliVTrack.h" 28 #include "AliVEvent.h" 31 #include "AliEMCALGeometry.h" 45 fStudyShape(kFALSE), fStudyShapeParam(kFALSE),
47 fStudyTCardCorrelation(kFALSE), fStudyExotic(kFALSE),
48 fStudyColRowFromCellMax(kFALSE), fStudyInvMass(kFALSE),
49 fStudyTimeCellHisto(kFALSE),
52 fM02Min(0), fNCellMin(0),
53 fMinDistToBad(0), fNEBinCuts(0),
54 fEMinShape(0), fEMaxShape(100),
55 fNCellMinShape(-1), fNCellsBins(3),
58 fdEdXMinEle(0), fdEdXMaxEle(0),
59 fdEdXMinHad(0), fdEdXMaxHad(0),
62 fInvMassMinECut(0), fInvMassMaxECut(0),
63 fInvMassMinM02Cut(0), fInvMassMaxM02Cut(0),
64 fInvMassMaxOpenAngle(0), fInvMassMaxTimeDifference(0),
66 fConstantTimeShift(0),
68 fClusterMomentum(), fClusterMomentum2(),
69 fCaloCellList(NULL), fCaloClusList(NULL),
74 fhEnergyTMEtaResidual1Cell(0), fhEnergyTMPhiResidual1Cell(0),
75 fhColRowExoticHighE1CellPosTime(0), fhColRowExoticHighE1CellNegTime(0), fhColRowExoticHighE1CellNulTime(0),
76 fhEnergyTMEtaResidualExotic(0), fhEnergyTMPhiResidualExotic(0),
77 fhColRowExoticHighEPosTime(0), fhColRowExoticHighENegTime(0), fhColRowExoticHighENulTime(0),
78 fhColRowHighEPosTime(0), fhColRowHighENegTime(0), fhColRowHighENulTime(0),
79 fhEnergyTMEtaResidualTCardCorrNoSelection1Cell(0), fhEnergyTMPhiResidualTCardCorrNoSelection1Cell(0),
80 fhEnergyTMEtaResidualTCardCorrNoSelectionExotic(0), fhEnergyTMPhiResidualTCardCorrNoSelectionExotic(0),
85 fhClusterMaxCellCloseCellRatioM02(0), fhClusterMaxCellECrossM02(0),
86 fhInvMassNCellSM(0), fhInvMassNCellSMSame(0),
87 fhColRowM02(0), fhColRowM02NCellCut(0),
88 fhEMaxCellTimeM02SM(0), fhEMaxCellTimeNCellSM(0), fhESecCellTimeNCellSM(0),
90 fhSMNCellModuleMax(0), fhSMNCellModuleOut(0),
91 fhSMNCellModuleOutModSame(0), fhSMNCellModuleOutModDiff(0),
92 fhSMECellModuleMax(0), fhSMECellModuleMaxTot(0),
93 fhSMECellModuleOut(0), fhSMECellModuleOutModSame(0), fhSMECellModuleOutModDiff(0),
94 fhSMNCellModuleMaxOutRat(0), fhSMNCellModuleMaxOutModSameRat(0), fhSMNCellModuleMaxOutModDiffRat(0),
95 fhSMECellModuleMaxRat(0), fhSMECellModuleMaxTotRat(0), fhSMECellModuleMaxTotRatClus(0),
96 fhSMECellModuleMaxOutRat(0), fhSMECellModuleMaxOutModSameRat(0), fhSMECellModuleMaxOutModDiffRat(0),
98 fhSMNCellModuleMaxLowM02(0), fhSMNCellModuleOutLowM02(0),
99 fhSMNCellModuleOutModSameLowM02(0), fhSMNCellModuleOutModDiffLowM02(0),
100 fhSMECellModuleMaxLowM02(0), fhSMECellModuleMaxTotLowM02(0),
101 fhSMECellModuleOutLowM02(0), fhSMECellModuleOutModSameLowM02(0), fhSMECellModuleOutModDiffLowM02(0),
102 fhSMNCellModuleMaxOutRatLowM02(0), fhSMNCellModuleMaxOutModSameRatLowM02(0),fhSMNCellModuleMaxOutModDiffRatLowM02(0),
103 fhSMECellModuleMaxRatLowM02(0), fhSMECellModuleMaxTotRatLowM02(0), fhSMECellModuleMaxTotRatClusLowM02(0),
104 fhSMECellModuleMaxOutRatLowM02(0), fhSMECellModuleMaxOutModSameRatLowM02(0),fhSMECellModuleMaxOutModDiffRatLowM02(0),
106 fhSMNCellModuleMaxHighM02(0), fhSMNCellModuleOutHighM02(0),
107 fhSMNCellModuleOutModSameHighM02(0), fhSMNCellModuleOutModDiffHighM02(0),
108 fhSMECellModuleMaxHighM02(0), fhSMECellModuleMaxTotHighM02(0),
109 fhSMECellModuleOutHighM02(0), fhSMECellModuleOutModSameHighM02(0), fhSMECellModuleOutModDiffHighM02(0),
110 fhSMNCellModuleMaxOutRatHighM02(0), fhSMNCellModuleMaxOutModSameRatHighM02(0),fhSMNCellModuleMaxOutModDiffRatHighM02(0),
111 fhSMECellModuleMaxRatHighM02(0), fhSMECellModuleMaxTotRatHighM02(0), fhSMECellModuleMaxTotRatClusHighM02(0),
112 fhSMECellModuleMaxOutRatHighM02(0), fhSMECellModuleMaxOutModSameRatHighM02(0),fhSMECellModuleMaxOutModDiffRatHighM02(0),
114 fhSMEMaxEClusterRat(0),
115 fhSMEMaxEClusterRatLowM02(0), fhSMEMaxEClusterRatHighM02(0),
117 fhTCardChannelM02(0), fhTCardChannelM02NoCut(0),
118 fhTCardChannelNCellModMax(0), fhTCardChannelNCell(0),
121 fhECellClusterRatio(0), fhECellClusterLogRatio(0),
122 fhEMaxCellClusterRatio(0), fhEMaxCellClusterLogRatio(0),
123 fhECellTotalRatio(0), fhECellTotalLogRatio(0),
124 fhECellTotalRatioMod(0), fhECellTotalLogRatioMod(0)
126 for(
Int_t i=0; i < 3; i++)
182 for(
Int_t i =0; i < 12; i++)
184 for(
Int_t j = 0; j < 4; j++)
186 for(
Int_t k = 0; k < 3; k++)
193 for(
Int_t j = 0; j < 5; j++)
201 for(
Int_t tm = 0; tm < 2; tm++)
297 for(
Int_t i = 0; i < 6; i++)
299 for(
Int_t j = 0; j < 6; j++)
372 for(
Int_t i = 0; i < 14; i++)
391 for(
Int_t j = 0; j < 6; j++)
393 for(
Int_t k = 0; k < 6; k++)
414 for(
Int_t i = 0; i < 12; i++)
440 for(
Int_t i = 0; i < 20; i++)
469 for(
Int_t j = 0; j < 4; j++)
513 Int_t icol = -1, irow = -1, iRCU = -1, icolAbs = -1, irowAbs = -1;
523 Int_t ncells = clus->GetNCells();
524 Double_t time = clus->GetTOF()*1.e9;
527 Float_t deta = clus->GetTrackDz();
528 Float_t dphi = clus->GetTrackDx();
536 if(energy >= 5 && energy < 8)
551 else if(exoticity > 0.97)
553 if(energy >= 5 && energy < 8)
568 for (
Int_t ipos = 0; ipos < ncells; ipos++)
570 Int_t absId = clus->GetCellsAbsId()[ipos];
578 if( absId == absIdMax || weight < 0.01 )
continue;
580 Int_t rowDiff = -100, colDiff = -100;
584 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
589 if(energy >= 5 && energy < 8)
596 if(energy >= 5 && energy < 8)
604 else if ( energy > 8 )
617 else if ( energy > 8 && ncells > 1 && exoticity < 0.97)
627 Int_t etaRegion = -1, phiRegion = -1;
630 if ( etaRegion !=0 ) return ;
640 else if(exoticity > 0.97)
651 if ( m02 > 0.001 ) lamR = m20/m02;
653 Int_t absIdList[ncells];
661 Int_t nCellWithWeight = 1;
693 Bool_t sameTCard2ndMax = kFALSE;
697 Float_t eCellSameRowSameTCardNearCol = 0.;
698 Float_t eCellSameRowDiffTCardNearCol = 0.;
699 Double_t tCellSameRowSameTCardNearCol = 0.;
700 Double_t tCellSameRowDiffTCardNearCol = 0.;
707 for (
Int_t ipos = 0; ipos < ncells; ipos++)
709 Int_t absId = clus->GetCellsAbsId()[ipos];
726 if( absId == absIdMax || weight < 0.01 )
continue;
747 Int_t rowDiff = -100, colDiff = -100;
778 Int_t indexType = -1;
783 if(TMath::Abs(rowDiff) == 1) nearRow = kTRUE;
784 if(TMath::Abs(colDiff) == 1) nearCol = kTRUE;
786 if ( rowDiff == 0 && colDiff != 0 )
788 if ( nearCol ) indexType = 6;
794 else if ( rowDiff != 0 && colDiff == 0 )
796 if ( nearRow ) indexType = 8;
804 if ( nearRow && nearCol) indexType = 10;
815 if ( rowDiff == 0 && colDiff != 0 )
817 if ( nearCol ) indexType = 0;
821 else if ( rowDiff != 0 && colDiff == 0 )
823 if ( nearRow ) indexType = 2;
828 if ( nearCol && nearRow ) indexType = 4;
833 if ( rowDiff == 0 && TMath::Abs(colDiff) == 1 )
837 eCellSameRowSameTCardNearCol = eCell;
838 tCellSameRowSameTCardNearCol = tCell;
842 eCellSameRowDiffTCardNearCol = eCell;
843 tCellSameRowDiffTCardNearCol = tCell;
849 Float_t eCellDiff = eCellMax - eCell;
850 Float_t eClusDiff = energy - eCell;
853 Float_t tCellDiff = tCellMax - tCell;
874 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
879 if(energy >= 5 && energy < 8)
886 if(energy >= 5 && energy < 8)
893 if ( eCell > emax2nd )
899 sameTCard2ndMax = kTRUE;
905 sameTCard2ndMax = kFALSE;
940 Float_t ratioNcells = nCellWithWeight/(ncells*1.);
945 if ( nCellWithWeight <= 1 )
return;
950 printf(
"AliAnaClusterShapeCorrelStudies: M02 %f, M20 %f, E %2.3f, ncell %d, n with weight %d; max cell E %2.3f\n",
951 m02,m20,energy,ncells,nCellWithWeight,eCellMax);
960 if(energy >= 5 && energy < 8)
967 if(energy >= 5 && energy < 8)
1033 if(eCellSameRowSameTCardNearCol > 0 && eCellSameRowDiffTCardNearCol > 0)
1035 Float_t eDiff = eCellSameRowSameTCardNearCol - eCellSameRowDiffTCardNearCol ;
1036 Float_t tDiff = tCellSameRowSameTCardNearCol - tCellSameRowDiffTCardNearCol ;
1053 Int_t nCorrInd = nCorr;
1054 if(nCorr > 4) nCorrInd = 5;
1056 Int_t nCorrNoInd = nCorrNo;
1057 if(nCorrNoInd > 4) nCorrNoInd = 5;
1107 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1109 Int_t absId = clus->GetCellsAbsId()[ipos];
1128 if( absId == absIdMax || weight < 0.01 )
continue;
1130 Float_t tDiffMaxSecondary = tCellMax - tCell;
1137 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
1140 Int_t rowDiff = -100, colDiff = -100;
1145 if ( energy >= 5 && energy < 8 && nCorr == 0 )
1147 else if ( energy >= 8 && nCorr == 0 )
1152 if ( energy >= 5 && energy < 8 && nCorrNo == 0 )
1154 else if ( energy >= 8 && nCorrNo == 0 )
1184 AliVCluster* clus2 = (AliVCluster*)
fCaloClusList->At(jclus);
1189 Double_t tof2 = clus2->GetTOF()*1.e9;
1197 if( absIdMax == absIdMax2
1198 || !
IsGoodCluster(absIdMax2, clus2->GetM02(), clus2->GetNCells())
1230 if ( energy >= 5 && energy < 8)
1232 else if ( energy >= 8 )
1250 if ( energy >= 5 && energy < 8)
1252 else if ( energy >= 8 )
1318 if ( exoticity > 0.97 )
1333 Float_t nCellRat = nCorr*1. / ((nCorr+nCorrNo)*1.);
1435 Int_t nCaloCellsPerCluster = clus->GetNCells();
1436 if ( nCaloCellsPerCluster < 2 )
return;
1453 for (
Int_t ipos = 0; ipos < nCaloCellsPerCluster; ipos++)
1455 Int_t absId = clus ->GetCellsAbsId()[ipos];
1460 if( absId == absIdMax || eCell < 0.01 )
continue;
1464 if( weight < 0.01 )
continue;
1469 if ( nCell < 1 )
return;
1472 Int_t nCellBin = -1;
1487 Int_t ietaMax=-1, iphiMax = 0, rcuMax = 0, icolAbs = -1, irowAbs = -1;
1489 ietaMax, iphiMax, rcuMax, icolAbs, irowAbs);
1491 Int_t iTowerMax = -1, iIphiMax = -1, iIetaMax = -1;
1493 GetEMCALGeometry()->GetCellIndex(absIdMax,smMax,iTowerMax,iIphiMax,iIetaMax);
1508 Int_t etaRegion = -1, phiRegion = -1;
1511 if ( etaRegion !=0 ) return ;
1515 Int_t nPerModule = 0;
1516 Int_t nPerModuleOut = 0;
1517 Int_t nPerModuleOutD = 0;
1518 Int_t nPerModuleOutS = 0;
1523 for (
Int_t ipos = 0; ipos < nCaloCellsPerCluster; ipos++)
1525 Int_t absId = clus ->GetCellsAbsId()[ipos];
1530 if( absId == absIdMax || eCell < 0.01 )
continue;
1539 Int_t iTower = -1, iIphi = -1, iIeta = -1;
1544 if ( iTowerMax == iTower )
1554 ePerModuleOut+=eCell;
1556 if(TMath::Abs(iphiMax-iphi) == 1 || TMath::Abs(ietaMax-ieta) == 1)
1558 Int_t rowDiff = -100, colDiff = -100;
1563 ePerModuleOutS+=eCell;
1568 ePerModuleOutD+=eCell;
1578 if ( weight < 0.01 )
continue;
1580 Float_t fracCell = eCell/eCellMax;
1583 Float_t fracClus = (energy-eCell)/energy;
1629 if(dIphi < TMath::Abs(iphi-iphiMax)) dIphi = TMath::Abs(iphi-iphiMax);
1630 if(iphi-iphiMax < 0 && dIphiNeg > iphi-iphiMax) dIphiNeg = iphi-iphiMax;
1631 if(iphi-iphiMax > 0 && dIphiPos < iphi-iphiMax) dIphiPos = iphi-iphiMax;
1633 if(dIeta < TMath::Abs(ieta-ietaMax)) dIeta = TMath::Abs(ieta-ietaMax);
1634 if(ieta-ietaMax < 0 && dIetaNeg > ieta-ietaMax) dIetaNeg = ieta-ietaMax;
1635 if(ieta-ietaMax > 0 && dIetaPos < ieta-ietaMax) dIetaPos = ieta-ietaMax;
1639 Int_t ietaShift = ieta;
1640 Int_t ietaMaxShift = ietaMax;
1642 if (ieta > ietaMax) ietaMaxShift+=48;
1643 else ietaShift +=48;
1645 if(dIeta < TMath::Abs(ietaShift-ietaMaxShift)) dIeta = TMath::Abs(ietaShift-ietaMaxShift);
1646 if(ietaShift-ietaMaxShift < 0 && dIetaNeg > ietaShift-ietaMaxShift) dIetaNeg = ietaShift-ietaMaxShift;
1647 if(ietaShift-ietaMaxShift > 0 && dIetaPos < ietaShift-ietaMaxShift) dIetaPos = ietaShift-ietaMaxShift;
1660 if ( nCellBin >= 0 && nCellBin < fNCellsBins )
1672 if ( nCellBin >= 0 && nCellBin < fNCellsBins )
1689 Int_t iTCard = rowTCard+8*colTCard;
1697 Float_t ePerModuleTot = ePerModule+eCellMax;
1835 if ( matchedPID == 0 )
1869 Float_t dIA = 1.*(dIphi-dIeta)/(dIeta+dIphi);
1870 Float_t dIATot = 1.*((dIphiPos-dIphiNeg)-(dIetaPos-dIetaNeg))/((dIetaPos-dIetaNeg)+(dIphiPos-dIphiNeg));
1879 AliDebug(1,Form(
"E %2.2f, nCell %d, dPhi %d, dEta %d, dIA %2.2f, match %d",energy,nCell, dIphi,dIeta,dIA,matchedPID));
1912 Float_t dispp= 0., dEta = 0., dPhi = 0.;
1913 Float_t sEta = 0., sPhi = 0., sEtaPhi = 0.;
1917 l0, l1, dispp, dEta, dPhi, sEta, sPhi, sEtaPhi);
1920 if(sEta+sPhi>0.0001) sEtaPhiA = (sPhi-sEta)/(sEta+sPhi);
1922 AliDebug(2,Form(
"Recalculate shower shape org: m02 %2.2f, m20 %2.2f, disp %2.2f;" 1923 " new: m02 %2.2f, m20 %2.2f, disp %2.2f; " 1924 "mEta %2.2f, mPhi %2.2f, mEtaPhi %2.2f, A_EtaPhi %2.2f; dEta %2.2f dPhi %2.2f",
1925 m02,m20,clus->GetDispersion(),l0,l1,dispp,sEta,sPhi,sEtaPhi,sEtaPhiA,dEta,dPhi));
1938 if ( nCell > fNCellMinShape )
1957 AliVCluster* clus2 = (AliVCluster*)
fCaloClusList->At(jclus);
1962 Double_t tof2 = clus2->GetTOF()*1.e9;
1970 if( absIdMax == absIdMax2
1971 || !
IsGoodCluster(absIdMax2, clus2->GetM02(), clus2->GetNCells())
2008 if (
IsDataMC() && mcIndex > -1 && mcIndex < 10)
2032 AliDebug(1,Form(
"No track found with matches %d",clus->GetNTracksMatched()));
2038 Double_t dedx = track->GetTPCsignal();
2058 else matchedPID =-1;
2060 AliDebug(1,Form(
"matches %d; dEdX %2.2f; matched PID %d",clus->GetNTracksMatched(),dedx,matchedPID));
2076 Int_t nCaloCellsPerCluster = 0 ;
2087 for(
Int_t iclus = 0; iclus < nCaloClusters; iclus++)
2089 AliDebug(1,Form(
"Cluster: %d/%d, data %d",iclus+1,nCaloClusters,
GetReader()->GetDataType()));
2091 AliVCluster* clus = (AliVCluster*)
fCaloClusList->At(iclus);
2094 if ( clus->GetDistanceToBadChannel() > -1 &&
2097 AliDebug(1,Form(
"Small distance to bad channel %2.2f < %2.2f",clus->GetDistanceToBadChannel(),
fMinDistToBad));
2103 if ( nModule < fFirstModule || nModule >
fLastModule )
2105 AliDebug(1,Form(
"Cluster module out of range %d",nModule));
2115 Double_t tof = clus->GetTOF()*1.e9;
2133 AliDebug(1,Form(
"Remove cluster with phi %2.2f and eta %2.2f", phi*TMath::RadToDeg(), eta));
2137 AliDebug(1,Form(
"cluster: E %2.3f, pT %2.3f, eta %2.3f, phi %2.3f",e,pt,eta,phi*TMath::RadToDeg()));
2141 nCaloCellsPerCluster = clus->GetNCells();
2145 AliDebug(1,Form(
"Accept cluster? %d",goodCluster));
2147 if(!goodCluster)
continue;
2176 Int_t matchedPID = -1;
2199 if ( ampMax > 0.01 )
2207 if (
fStudyShape && matchedPID >= 0 && matchedPID < 3 )
2228 const Int_t buffersize = 255;
2229 char onePar[buffersize] ;
2231 snprintf(onePar,buffersize,
"--- AliAnaClusterShapeCorrelStudies ---:") ;
2252 return new TObjString(parList) ;
2262 outputContainer->SetName(
"ClusterShapeStudies") ;
2273 TString mcParticleStringLabel[] = {
"Merged #gamma#gamma",
"Decay #gamma",
"Direct #gamma",
"e^{#pm}",
"#pi^{#pm}",
"k^{#pm}",
"p",
"#bar{p}",
"n",
"#bar{n}"};
2277 TString matchCase[] = {
"Neutral",
"Electron",
"Hadron"};
2305 Float_t rowcellmin = fNMaxRowsFullMin-1.5;
2314 Int_t nShShBins = 200;
2324 Int_t cellBins = 15;
2334 (
"hColRowHighEPosTime",
2335 "column vs row, exo < 0.97, E > 8 GeV, t > 5 ns",
2336 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2342 (
"hColRowHighENegTime",
2343 "column vs row, exo < 0.97, E > 8 GeV, t < -5 ns",
2344 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2350 (
"hColRowHighENulTime",
2351 "column vs row, exo < 0.97, E > 8 GeV, -5 < t < 5 ns",
2352 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2360 (
"hColRowExoticHighE1CellPosTime",
2361 "column vs row, 1 cell, E > 8 GeV, t > 5 ns",
2362 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2368 (
"hColRowExoticHighEPosTime",
2369 "column vs row, exo > 0.97, E > 8 GeV, t > 5 ns",
2370 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2376 (
"hColRowExoticHighE1CellNegTime",
2377 "column vs row, 1 cell, E > 8 GeV, t < -5 ns",
2378 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2384 (
"hColRowExoticHighENegTime",
2385 "column vs row, exo > 0.97, E > 8 GeV, t < -5 ns",
2386 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2392 (
"hColRowExoticHighE1CellNulTime",
2393 "column vs row, 1 cell, E > 8 GeV, -5 < t < 5 ns",
2394 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2400 (
"hColRowExoticHighENulTime",
2401 "column vs row, exo > 0.97, E > 8 GeV, -5 < t < 5 ns",
2402 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2408 TString add[] = {
"",
"TrackMatched"};
2409 for(
Int_t tm = 0; tm < 2; tm++)
2412 (Form(
"hColRowTCardCorrNoSelectionLowE%s",add[tm].
Data()),
2413 Form(
"column vs row, max E cell for TCard correlation selected clusters, 5 < E < 8 GeV %s",add[tm].
Data()),
2414 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2420 (Form(
"hColRowTCardCorrNoSelectionHighE%s",add[tm].
Data()),
2421 Form(
"column vs row, max E cell for TCard correlation selected clusters, E > 8 GeV %s",add[tm].
Data()),
2422 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2430 (Form(
"hNCellsTCardCorrNoSelection%s",add[tm].
Data()),
2431 Form(
"# custer # cells vs #it{E} %s",add[tm].
Data()),
2432 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
2438 (Form(
"hNCellsTCardCorrWithWeightNoSelection%s",add[tm].
Data()),
2439 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
2440 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
2446 (Form(
"hNCellsTCardCorrRatioWithWeightNoSelection%s",add[tm].
Data()),
2447 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
2448 nptbins,ptmin,ptmax, 100,0,1);
2454 (Form(
"hTimeTCardCorrNoSelection%s",add[tm].
Data()),
2455 Form(
"#it{time} vs #it{E} %s",add[tm].
Data()),
2456 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2462 (Form(
"hLambda0TCardCorrNoSelection%s",add[tm].
Data()),
2463 Form(
"#sigma^{2}_{long} vs #it{E} %s",add[tm].
Data()),
2464 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2470 (Form(
"hLambda1TCardCorrNoSelection%s",add[tm].
Data()),
2471 Form(
"#sigma^{2}_{short} vs #it{E} %s",add[tm].
Data()),
2472 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2478 (Form(
"hLambda0NLM1TCardCorrNoSelection%s",add[tm].
Data()),
2479 Form(
"#sigma^{2}_{long} vs #it{E}, nlm=1 %s",add[tm].
Data()),
2480 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2486 (Form(
"hLambda1NLM1TCardCorrNoSelection%s",add[tm].
Data()),
2487 Form(
"#sigma^{2}_{short} vs #it{E}, nlm=1 %s",add[tm].
Data()),
2488 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2494 (Form(
"hLambda0NLM2TCardCorrNoSelection%s",add[tm].
Data()),
2495 Form(
"#sigma^{2}_{long} vs #it{E}, nlm=2 %s",add[tm].
Data()),
2496 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2502 (Form(
"hLambda1NLM2TCardCorrNoSelection%s",add[tm].
Data()),
2503 Form(
"#sigma^{2}_{short} vs #it{E}, nlm=2 %s",add[tm].
Data()),
2504 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2511 (Form(
"hLambdaRTCardCorrNoSelection%s",add[tm].
Data()),
2512 Form(
"#sigma^{2}_{short}/#sigma^{2}_{long} vs #it{E} %s",add[tm].
Data()),
2513 nptbins,ptmin,ptmax,110,0,1.1);
2519 (Form(
"hNLocMaxTCardCorrNoSelection%s",add[tm].
Data()),
2520 Form(
"#it{n}_{LM} vs E %s",add[tm].
Data()),
2521 nptbins,ptmin,ptmax,10,0,10);
2527 (Form(
"hEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2528 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
2529 nptbins,ptmin,ptmax,100,0,1);
2535 (Form(
"hEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2536 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2537 nptbins,ptmin,ptmax,100,0,1);
2543 (Form(
"hEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2544 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2545 nptbins,ptmin,ptmax,100,0,1);
2551 (Form(
"hE2ndRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2552 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
2553 nptbins,ptmin,ptmax,100,0,1);
2559 (Form(
"hE2ndRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2560 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2561 nptbins,ptmin,ptmax,100,0,1);
2567 (Form(
"hE2ndRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2568 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2569 nptbins,ptmin,ptmax,100,0,1);
2576 (Form(
"hE2ndEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2577 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
2578 nptbins,ptmin,ptmax,100,0,1);
2584 (Form(
"hE2ndEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2585 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2586 nptbins,ptmin,ptmax,100,0,1);
2592 (Form(
"hE2ndEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2593 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2594 nptbins,ptmin,ptmax,100,0,1);
2602 (Form(
"hE2ndSameRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2603 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()),
2604 nptbins,ptmin,ptmax,100,0,1);
2610 (Form(
"hE2ndSameRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2611 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()),
2612 nptbins,ptmin,ptmax,100,0,1);
2618 (Form(
"hE2ndSameRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2619 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()),
2620 nptbins,ptmin,ptmax,100,0,1);
2627 (Form(
"hE2ndSameEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2628 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()),
2629 nptbins,ptmin,ptmax,100,0,1);
2635 (Form(
"hE2ndSameEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2636 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()),
2637 nptbins,ptmin,ptmax,100,0,1);
2643 (Form(
"hE2ndSameEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2644 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()),
2645 nptbins,ptmin,ptmax,100,0,1);
2653 (Form(
"hECellClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2654 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
2655 nptbins,ptmin,ptmax,100,0,1);
2661 (Form(
"hECellClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2662 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2663 nptbins,ptmin,ptmax,100,0,1);
2669 (Form(
"hECellClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2670 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2671 nptbins,ptmin,ptmax,100,0,1);
2677 (Form(
"hLogECellNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2678 Form(
"log(#it{E}_{cell}) vs E_{cluster}, #it{n}_{LM}=1, w > 0.01 %s",add[tm].
Data()),
2679 nptbins,ptmin,ptmax,150,-3,3);
2685 (Form(
"hLogECellNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2686 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}=2, w > 0.01 %s",add[tm].
Data()),
2687 nptbins,ptmin,ptmax,150,-3,3);
2693 (Form(
"hLogECellNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2694 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}>2, w > 0.01 %s",add[tm].
Data()),
2695 nptbins,ptmin,ptmax,150,-3,3);
2702 (Form(
"hECellWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2703 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
2704 nptbins,ptmin,ptmax,90,0,4.5);
2710 (Form(
"hECellWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2711 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2712 nptbins,ptmin,ptmax,90,0,4.5);
2718 (Form(
"hECellWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2719 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster}))vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2720 nptbins,ptmin,ptmax,90,0,4.5);
2727 (Form(
"hECellSameClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2728 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()),
2729 nptbins,ptmin,ptmax,100,0,1);
2735 (Form(
"hECellSameClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2736 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, cell from same T-Card as leading, #it{n}_{LM}=2 %s",add[tm].
Data()),
2737 nptbins,ptmin,ptmax,100,0,1);
2743 (Form(
"hECellSameClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2744 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, cell from same T-Card as leading, #it{n}_{LM}>2 %s",add[tm].
Data()),
2745 nptbins,ptmin,ptmax,100,0,1);
2751 (Form(
"hLogECellSameNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2752 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()),
2753 nptbins,ptmin,ptmax,150,-3,3);
2759 (Form(
"hLogECellSameNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2760 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()),
2761 nptbins,ptmin,ptmax,150,-3,3);
2767 (Form(
"hLogECellSameNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2768 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()),
2769 nptbins,ptmin,ptmax,150,-3,3);
2776 (Form(
"hECellSameWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2777 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()),
2778 nptbins,ptmin,ptmax,90,0,4.5);
2784 (Form(
"hECellSameWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2785 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()),
2786 nptbins,ptmin,ptmax,90,0,4.5);
2792 (Form(
"hECellSameWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2793 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()),
2794 nptbins,ptmin,ptmax,90,0,4.5);
2800 (Form(
"hExoticTCardCorrNoSelection%s",add[tm].
Data()),
2801 Form(
"exoticity vs #it{E} %s",add[tm].
Data()),
2802 nptbins,ptmin,ptmax,200,-1,1);
2808 (Form(
"hNCellsTCardSameAndDiffFraction%s",add[tm].
Data()),
2809 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, %s",add[tm].
Data()),
2810 nptbins,ptmin,ptmax,100,0,1);
2816 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterE%s",add[tm].
Data()),
2817 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()),
2818 nptbins,ptmin,ptmax,200,-10,10);
2824 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterE%s",add[tm].
Data()),
2825 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()),
2826 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2832 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxE%s",add[tm].
Data()),
2833 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()),
2834 nptbins,ptmin,ptmax,200,-10,10);
2840 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxE%s",add[tm].
Data()),
2841 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()),
2842 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2850 (Form(
"hEnergyTime1Cell%s",add[tm].
Data()),
2851 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
2852 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2858 (Form(
"hEnergyTimeExotic%s",add[tm].
Data()),
2859 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
2860 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2866 (Form(
"hEnergyTimeTCardCorrNoSelection1Cell%s",add[tm].
Data()),
2867 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
2868 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2874 (Form(
"hEnergyTimeTCardCorrNoSelectionExotic%s",add[tm].
Data()),
2875 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
2876 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2882 (Form(
"hColRowExoticLowE1Cell%s",add[tm].
Data()),
2883 Form(
"column vs row, 1 cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2884 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2890 (Form(
"hColRowExoticHighE1Cell%s",add[tm].
Data()),
2891 Form(
"column vs row, 1 cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2892 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2898 (Form(
"hColRowExoticLowE%s",add[tm].
Data()),
2899 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2900 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2906 (Form(
"hColRowExoticHighE%s",add[tm].
Data()),
2907 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2908 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2914 (Form(
"hColRowExotic2ndCellDiffLowE%s",add[tm].
Data()),
2915 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2916 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2922 (Form(
"hColRowExotic2ndCellDiffHighE%s",add[tm].
Data()),
2923 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2924 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2930 (Form(
"hColRowExotic2ndCellSameLowE%s",add[tm].
Data()),
2931 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2932 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2938 (Form(
"hColRowExotic2ndCellSameHighE%s",add[tm].
Data()),
2939 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2940 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2946 (Form(
"hColRowTCardCorrNoSelectionExoticLowE%s",add[tm].
Data()),
2947 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2948 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2954 (Form(
"hColRowTCardCorrNoSelectionExoticHighE%s",add[tm].
Data()),
2955 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2956 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2962 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffLowE%s",add[tm].
Data()),
2963 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2964 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2970 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffHighE%s",add[tm].
Data()),
2971 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2972 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2978 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameLowE%s",add[tm].
Data()),
2979 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2980 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2986 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameHighE%s",add[tm].
Data()),
2987 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2988 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2994 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameLowE%s",add[tm].
Data()),
2995 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2996 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3002 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameHighE%s",add[tm].
Data()),
3003 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3004 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3010 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffLowE%s",add[tm].
Data()),
3011 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3012 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3018 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffHighE%s",add[tm].
Data()),
3019 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3020 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3026 (Form(
"hNCellsTCardSameAndDiffFraction_Exotic%s",add[tm].
Data()),
3027 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, exo > 0.97 %s",add[tm].
Data()),
3028 nptbins,ptmin,ptmax,100,0,1);
3034 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterEExo%s",add[tm].
Data()),
3035 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()),
3036 nptbins,ptmin,ptmax,200,-10,10);
3042 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterEExo%s",add[tm].
Data()),
3043 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()),
3044 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
3050 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxEExo%s",add[tm].
Data()),
3051 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()),
3052 nptbins,ptmin,ptmax,200,-10,10);
3058 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxEExo%s",add[tm].
Data()),
3059 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()),
3060 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
3066 for(
Int_t i = 0; i < 6; i++)
3068 for(
Int_t j = 0; j < 6; j++)
3071 (Form(
"hLambda0TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3072 Form(
"#sigma^{2}_{long} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3073 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3079 (Form(
"hLambda1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3080 Form(
"#sigma^{2}_{short} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3081 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3087 (Form(
"hLambda0NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3088 Form(
"#sigma^{2}_{long} vs #it{E}, nlm=1, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3089 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3095 (Form(
"hLambda1NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3096 Form(
"#sigma^{2}_{short} vs #it{E}, nlm=1, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3097 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3103 (Form(
"hLambda0NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3104 Form(
"#sigma^{2}_{long} vs #it{E}, nlm=2, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3105 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3111 (Form(
"hLambda1NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3112 Form(
"#sigma^{2}_{short} vs #it{E}, nlm=2, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3113 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3128 (Form(
"hNLocMaxTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3129 Form(
"#it{n}_{LM} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3130 nptbins,ptmin,ptmax,10,0,10);
3136 (Form(
"hEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3137 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()),
3138 nptbins,ptmin,ptmax,100,0,1);
3144 (Form(
"hEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3145 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()),
3146 nptbins,ptmin,ptmax,100,0,1);
3152 (Form(
"hEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3153 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()),
3154 nptbins,ptmin,ptmax,100,0,1);
3160 (Form(
"hE2ndRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3161 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()),
3162 nptbins,ptmin,ptmax,100,0,1);
3168 (Form(
"hE2ndRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3169 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()),
3170 nptbins,ptmin,ptmax,100,0,1);
3176 (Form(
"hE2ndRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3177 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()),
3178 nptbins,ptmin,ptmax,100,0,1);
3184 (Form(
"hE2ndEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3185 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()),
3186 nptbins,ptmin,ptmax,100,0,1);
3192 (Form(
"hE2ndEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3193 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()),
3194 nptbins,ptmin,ptmax,100,0,1);
3200 (Form(
"hE2ndEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3201 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()),
3202 nptbins,ptmin,ptmax,100,0,1);
3209 (Form(
"hECellClusRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3210 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()),
3211 nptbins,ptmin,ptmax,100,0,1);
3217 (Form(
"hECellClusRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3218 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()),
3219 nptbins,ptmin,ptmax,100,0,1);
3225 (Form(
"hECellClusRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3226 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()),
3227 nptbins,ptmin,ptmax,100,0,1);
3233 (Form(
"hLogECellNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3234 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()),
3235 nptbins,ptmin,ptmax,150,-3,3);
3241 (Form(
"hLogECellNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3242 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()),
3243 nptbins,ptmin,ptmax,150,-3,3);
3249 (Form(
"hLogECellNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3250 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()),
3251 nptbins,ptmin,ptmax,150,-3,3);
3257 (Form(
"hECellWeightNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3258 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()),
3259 nptbins,ptmin,ptmax,90,0,4.5);
3265 (Form(
"hECellWeightNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3266 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()),
3267 nptbins,ptmin,ptmax,90,0,4.5);
3273 (Form(
"hECellWeightNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3274 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()),
3275 nptbins,ptmin,ptmax,90,0,4.5);
3283 (Form(
"hMassEClusTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3284 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()),
3285 nptbins,ptmin,ptmax,nmassbins,massmin,massmax);
3300 (Form(
"hTimeDiffTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3301 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()),
3302 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
3308 (Form(
"hColRowTCardCorrelNCellLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3309 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3310 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3316 (Form(
"hColRowTCardCorrelNCellHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3317 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3318 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3326 (Form(
"hExoticTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3327 Form(
"exoticity vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3328 nptbins,ptmin,ptmax,200,-1,1);
3334 (Form(
"hTimeDiffExoTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3335 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()),
3336 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
3342 (Form(
"hColRowTCardCorrelNCellExoticLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3343 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3344 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3350 (Form(
"hColRowTCardCorrelNCellExoticHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3351 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3352 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3567 (Form(
"hLambda0Exoticity_EBin%d%s",i,add[tm].
Data()),
3568 Form(
"#sigma^{2}_{long} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3569 200,-1,1,ssbins,ssmin,ssmax);
3575 (Form(
"hLambda1Exoticity_EBin%d%s",i,add[tm].
Data()),
3576 Form(
"#sigma^{2}_{short} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3577 200,-1,1,ssbins,ssmin,ssmax);
3591 (Form(
"hNCellsExoticity_EBin%d%s",i,add[tm].
Data()),
3592 Form(
"#it{n}_{cells} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3593 200,-1,1,nceclbins,nceclmin,nceclmax);
3594 fhNCellsExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
3599 (Form(
"hTimeExoticity_EBin%d%s",i,add[tm].
Data()),
3601 200,-1,1,100,-25,25);
3602 fhTimeExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
3607 (Form(
"hNCellsTCardSameAndDiff_Exotic_EBin%d%s",i,add[tm].
Data()),
3608 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()),
3609 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
3649 (Form(
"hNCellsTCardSameAndDiff_EBin%d%s",i,add[tm].
Data()),
3650 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()),
3651 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
3657 (Form(
"hLambda0Lambda1_EBin%d%s",i,add[tm].
Data()),
3658 Form(
"#sigma^{2}_{long} vs #sigma^{2}_{short}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3659 ssbins,ssmin,ssmax,ssbins,ssmin,ssmax);
3675 for(
Int_t j = 0; j < 6; j++)
3677 for(
Int_t k = 0; k < 6; k++)
3680 (Form(
"hLambda0Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
3681 Form(
"#sigma^{2}_{long} vs #it{exoticity}, #it{n}_{cell} TCard same = %d, diff =%d, #it{E}>8 GeV %s",j,k,add[tm].
Data()),
3682 200,-1,1,ssbins,ssmin,ssmax);
3688 (Form(
"hLambda1Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
3689 Form(
"#sigma^{2}_{short} vs #it{exoticity}, #it{n}_{cell} TCard same = %d, diff =%d, #it{E}>8 GeV %s",j,k,add[tm].
Data()),
3690 200,-1,1,ssbins,ssmin,ssmax);
3787 for(
Int_t i = 0; i < 12; i++)
3790 (Form(
"hTCardCorrECellMaxDiff_Case%d%s",i,add[tm].
Data()),
3791 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()),
3792 nptbins,ptmin,ptmax,210,-1,20);
3798 (Form(
"hTCardCorrEClusterDiff_Case%d%s",i,add[tm].
Data()),
3799 Form(
"#it{E}_{cluster}-#it{E}_{cell} vs #it{E}_{cluster}, for (un)correlated cells in TCard, case %d %s",i,add[tm].
Data()),
3800 nptbins,ptmin,ptmax,210,-1,20);
3822 (Form(
"hTCardCorrTCellMaxDiff_Case%d%s",i,add[tm].
Data()),
3823 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()),
3824 nptbins,ptmin,ptmax,1000,-100,100);
3832 (Form(
"hTCardCorrECellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
3833 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()),
3834 nptbins,ptmin,ptmax,210,-1,20);
3840 (Form(
"hTCardCorrEClusterDiffExo_Case%d%s",i,add[tm].
Data()),
3841 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()),
3842 nptbins,ptmin,ptmax,210,-1,20);
3864 (Form(
"hTCardCorrTCellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
3865 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()),
3866 nptbins,ptmin,ptmax,1000,-100,100);
3877 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3883 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3889 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3895 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3901 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3907 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3913 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3919 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3927 (Form(
"hTMPhiResidual_EBin%d",i),
3928 Form(
"#Delta #varphi_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3929 200,-1,1,nresphibins,resphimin,resphimax);
3935 (Form(
"hTMEtaResidual_EBin%d",i),
3936 Form(
"#Delta #eta_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3937 200,-1,1,nresetabins,resetamin,resetamax);
3943 (Form(
"hTMPhiResidual_LooseCut_EBin%d",i),
3944 Form(
"#Delta #varphi_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3945 200,-1,1,nresphibins,resphimin,resphimax);
3951 (Form(
"hTMEtaResidual_LooseCut_EBin%d",i),
3952 Form(
"#Delta #eta_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3953 200,-1,1,nresetabins,resetamin,resetamax);
3981 (Form(
"hColRowM02"),
3983 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax,40,0.,2.);
3990 (Form(
"hColRowM02NCellCut"),
3991 Form(
"column vs row vs M02, %2.2f<#it{E}<%2.2f GeV #it{n}_{cells}^{w>0.01} > %d for Neutral",
3993 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax,40,0.,2.);
4003 Form(
"%2.2f<#it{E}_{1}<%2.2f GeV, %2.2f<#it{E}_{2}<%2.2f GeV, %2.2f<#sigma^{2}_{long}<%2.2f" 4004 "#it{M}_{#gamma #gamma} vs #it{n}_{1, cells}^{w>0.01} vs SM number trig cluster",
4013 (
"hInvMassNCellSMSame",
4014 Form(
"%2.2f<#it{E}_{1}<%2.2f GeV, %2.2f<#it{E}_{2}<%2.2f GeV, %2.2f<#sigma^{2}_{long}<%2.2f" 4015 "#it{M}_{#gamma #gamma} vs #it{n}_{1, cells}^{w>0.01} vs SM number both cluster",
4029 (
"hEMaxCellTimeM02SM",
4030 Form(
"time vs #sigma^{2}_{long} vs SM number, " 4039 (
"hEMaxCellTimeNCellSM",
4040 Form(
"time vs #it{n}_{cells}^{w>0.01} vs SM number, " 4049 (
"hESecCellTimeNCellSM",
4050 Form(
"secondary cell time vs #it{n}_{cells}^{w>0.01} vs SM number, " 4061 Int_t nbinsdeltacells = 19 ;
4067 for(
Int_t col = 0; col < 2; col++)
4070 (Form(
"hColRowFromCellMaxLowM02_Col%d",col),
4071 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #it{n}_{cells}^{w>0.01} for %2.2f < #sigma^{2}_{long} < %2.2f, " 4073 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells);
4079 (Form(
"hColRowFromCellMaxHighM02_Col%d",col),
4080 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #it{n}_{cells}^{w>0.01} for %2.2f < #sigma^{2}_{long} < %2.2f, " 4082 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells);
4093 (Form(
"hNCellsPerClusterM02M20_SM%d",i),
4094 Form(
" vs #sigma^{2}_{long} vs #it{n}_{cells}^{w>0.01} vs #sigma^{2}_{short}, " 4096 (
Int_t)nShShBins/1.5,0,maxShSh/1.5, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4103 (Form(
"hNCellsPerClusterM02NLM_SM%d",i),
4104 Form(
" vs #sigma^{2}_{long} vs #it{n}_{cells}^{w>0.01} vs N_{lm}, " 4106 10,0,10, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4113 (Form(
"hEMaxCellEClusterM02NCell_SM%d",i),
4114 Form(
"(#it{E}_{cluster} - #it{E}_{max cell})/#it{E}_{cluster} vs #sigma^{2}_{long} vs #it{n}_{cells}^{w>0.01}, " 4116 50,0,1., cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4133 (Form(
"hESecCellEMaxCellM02NCell_SM%d",i),
4134 Form(
"#it{E}_{cell}/#it{E}_{cell max} vs #sigma^{2}_{long} vs #it{n}_{cells}^{w>0.01}, " 4136 50,0,1., cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4143 (Form(
"hESecCellEClusterM02NCell_SM%d",i),
4144 Form(
"(#it{E}_{cluster} - #it{E}_{cell})/#it{E}_{cluster} vs #sigma^{2}_{long} vs #it{n}_{cells}^{w>0.01}, " 4146 50,0,1., cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4163 (Form(
"hEMaxESecCellNCellLowM02_SM%d",i),
4164 Form(
"#it{E}_{cell}^{max} vs #it{E}_{cell}^{secondary} vs #it{n}_{cells}^{w>0.01} vs %2.2f < #sigma^{2}_{long} < %2.2f, " 4166 24,0,12, 48,0,12,cellBins,cellMin,cellMax);
4173 (Form(
"hEMaxECrossNCellLowM02_SM%d",i),
4174 Form(
"#it{E}_{cell}^{max} vs exoticity vs #it{n}_{cells}^{w>0.01} vs %2.2f < #sigma^{2}_{long} < %2.2f, " 4176 24,0,12, 40,0.6,1.,cellBins,cellMin,cellMax);
4183 (Form(
"hEMaxESecCellNCellHighM02_SM%d",i),
4184 Form(
"#it{E}_{cell}^{max} vs #it{E}_{cell}^{secondary} vs #it{n}_{cells}^{w>0.01} vs " 4185 "%2.2f < #sigma^{2}_{long} < %2.2f, %2.2f<#it{E}<%2.2f GeV, SM=%d",
4187 24,0,12, 48,0,12,cellBins,cellMin,cellMax);
4194 (Form(
"hEMaxECrossNCellHighM02_SM%d",i),
4195 Form(
"#it{E}_{cell}^{max} vs exoticity vs #it{n}_{cells}^{w>0.01} " 4196 "for %2.2f < #sigma^{2}_{long} < %2.2f, %2.2f<#it{E}<%2.2f GeV, SM=%d",
4198 24,0,12, 40,0.6,1.,cellBins,cellMin,cellMax);
4206 Int_t nbinsdeltacells = 19 ;
4209 for(
Int_t col = 0; col < 2; col++)
4212 (Form(
"hNCellsPerClusterM02_SM%d_Col%d",i,col),
4213 Form(
"#it{E} vs #it{n}_{cells} vs #sigma^{2}_{long}, SM=%d, colum %d",i,col),
4214 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4221 (Form(
"hNCellsPerClusterM20_SM%d_Col%d",i,col),
4222 Form(
"#it{E} vs #it{n}_{cells} vs #sigma^{2}_{short}, SM=%d, colum %d",i,col),
4223 nEbins,minE,maxE,cellBins,cellMin,cellMax,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4230 (Form(
"hNCellsPerClusterNLM_SM%d_Col%d",i,col),
4231 Form(
"#it{E} vs #it{n}_{cells} vs n_{lm}, SM=%d, colum %d",i,col),
4232 nEbins,minE,maxE,cellBins,cellMin,cellMax,10,0,10);
4240 (Form(
"hColRowFromCellMaxLowM02_SM%d_Col%d",i,col),
4241 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #it{n}_{cells}^{w>0.01} for %2.2f < #sigma^{2}_{long} < %2.2f, " 4242 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d",
4244 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,cellBins,cellMin,cellMax);
4251 (Form(
"hColRowFromCellMaxHighM02_SM%d_Col%d",i,col),
4252 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #it{n}_{cells}^{w>0.01} " 4253 "for %2.2f < #sigma^{2}_{long} < %2.2f, %2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d",
4255 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,cellBins,cellMin,cellMax);
4264 (Form(
"hColRowFromCellMaxEMaxSecDiffLowM02_SM%d_Col%d_NCellBin%d",i,col,j),
4265 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #Delta #it{E}_{max-secondary} " 4266 "for %2.2f < #sigma^{2}_{long} < %2.2f, %2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d,%d <= #it{n}_{cells}^{w>0.01} < %d",
4268 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,120,0,12);
4275 (Form(
"hColRowFromCellMaxEMaxSecDiffHighM02_SM%d_Col%d_NCellBin%d",i,col,j),
4276 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #Delta #it{E}_{max-secondary} " 4277 "for %2.2f < #sigma^{2}_{long} < %2.2f, %2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, %d <= #it{n}_{cells}^{w>0.01} < %d",
4279 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,120,0,12);
4286 (Form(
"hColRowFromCellMaxEMaxSecDiffFracLowM02_SM%d_Col%d_NCellBin%d",i,col,j),
4287 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #Delta #it{E}_{max-secondary}/#it{E}_{max} " 4288 "for %2.2f < #sigma^{2}_{long} < %2.2f, %2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d,%d <= #it{n}_{cells}^{w>0.01} < %d",
4290 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,40,0,1);
4297 (Form(
"hColRowFromCellMaxEMaxSecDiffFracHighM02_SM%d_Col%d_NCellBin%d",i,col,j),
4298 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #Delta #it{E}_{max-secondary}/#it{E}_{max}" 4299 " for %2.2f < #sigma^{2}_{long} < %2.2f, %2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, %d <= #it{n}_{cells}^{w>0.01} < %d",
4301 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,40,0,1);
4308 (Form(
"hColRowFromCellMaxECellClusterRatLowM02_SM%d_Col%d_NCellBin%d",i,col,j),
4309 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #it{E}_{cell}/#it{E}_{cluster} " 4310 "for %2.2f < #sigma^{2}_{long} < %2.2f, %2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, %d <= #it{n}_{cells}^{w>0.01} < %d",
4312 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,50,0,1);
4319 (Form(
"hColRowFromCellMaxECellClusterRatHighM02_SM%d_Col%d_NCellBin%d",i,col,j),
4320 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #it{E}_{cell}/#it{E}_{cluster} " 4321 "for %2.2f < #sigma^{2}_{long} < %2.2f, %2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, %d <= #it{n}_{cells}^{w>0.01} < %d",
4323 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,50,0,1);
4340 (Form(
"hDeltaIEtaDeltaIPhi_%s",matchCase[imatch].
Data()),
4341 Form(
"Cluster max size with respect main cell in columns vs rows vs E for %s, #it{n}_{cells}^{w>0.01}>%d",matchCase[imatch].
Data(),
fNCellMinShape),
4342 nEbins,minE,maxE,cellBins,cellMin,cellMax,cellBins,cellMin,cellMax);
4349 (Form(
"hDeltaIEtaDeltaIPhiTot_%s",matchCase[imatch].
Data()),
4350 Form(
"Cluster size in columns vs rows, minus main cell, vs E for %s, #it{n}_{cells}^{w>0.01}>%d",matchCase[imatch].
Data(),
fNCellMinShape),
4351 nEbins,minE,maxE,cellBins,cellMin,cellMax,cellBins,cellMin,cellMax);
4358 (Form(
"hDeltaIA_%s",matchCase[imatch].
Data()),
4359 Form(
"Cluster *asymmetry* in cell units vs E for %s",matchCase[imatch].
Data()),
4360 nptbins,ptmin,ptmax,asyBins,asyMin,asyMax);
4361 fhDeltaIA[imatch]->SetXTitle(
"#it{E}_{cluster}");
4362 fhDeltaIA[imatch]->SetYTitle(
"#it{A}_{cell in cluster}");
4363 outputContainer->Add(
fhDeltaIA[imatch]);
4366 (Form(
"hDeltaIATot_%s",matchCase[imatch].
Data()),
4367 Form(
"Cluster *total asymmetry* in cell units vs E for %s",matchCase[imatch].
Data()),
4368 nptbins,ptmin,ptmax,asyBins,asyMin,asyMax);
4370 fhDeltaIATot[imatch]->SetYTitle(
"#it{A}_{cell in cluster}^{total}");
4374 (Form(
"hDeltaIAM02_%s",matchCase[imatch].
Data()),
4375 Form(
"Cluster *asymmetry* in cell units vs #sigma^{2}_{long} for %s",matchCase[imatch].
Data()),
4376 nEbins,minE,maxE,nShShBins,minShSh,maxShSh,asyBins,asyMin,asyMax);
4379 fhDeltaIAM02[imatch]->SetZTitle(
"#it{A}_{cell in cluster}");
4383 (Form(
"hDeltaIATotM02_%s",matchCase[imatch].
Data()),
4384 Form(
"Cluster *total asymmetry* in cell units vs #sigma^{2}_{long} for %s",matchCase[imatch].
Data()),
4385 nEbins,minE,maxE,nShShBins,minShSh,maxShSh,asyBins,asyMin,asyMax);
4388 fhDeltaIATotM02[imatch]->SetZTitle(
"#it{A}_{cell in cluster}^{total}");
4392 (Form(
"hDeltaIAM20_%s",matchCase[imatch].
Data()),
4393 Form(
"Cluster *asymmetry* in cell units vs #sigma^{2}_{short} for %s",matchCase[imatch].
Data()),
4394 nEbins,minE,maxE,(
Int_t)nShShBins/1.5,0,maxShSh/1.5,asyBins,asyMin,asyMax);
4397 fhDeltaIAM20[imatch]->SetZTitle(
"#it{A}_{cell in cluster}");
4401 (Form(
"hDeltaIATotM20_%s",matchCase[imatch].
Data()),
4402 Form(
"Cluster *total asymmetry* in cell units vs #sigma^{2}_{short} for %s",matchCase[imatch].
Data()),
4403 nEbins,minE,maxE,(
Int_t)nShShBins/1.5,0,maxShSh/1.5,asyBins,asyMin,asyMax);
4406 fhDeltaIATotM20[imatch]->SetZTitle(
"#it{A}_{cell in cluster}^{total}");
4410 (Form(
"hDeltaIANCells_%s",matchCase[imatch].
Data()),
4411 Form(
"Cluster *asymmetry* in cell units vs N cells in cluster for %s",matchCase[imatch].
Data()),
4412 nEbins,minE,maxE,cellBins,cellMin,cellMax,asyBins,asyMin,asyMax);
4419 (Form(
"hDeltaIATotNCells_%s",matchCase[imatch].
Data()),
4420 Form(
"Cluster *total asymmetry* in cell units vs N cells in cluster for %s",matchCase[imatch].
Data()),
4421 nEbins,minE,maxE,cellBins,cellMin,cellMax,asyBins,asyMin,asyMax);
4430 (Form(
"hDeltaIAOrigin_%s",matchCase[imatch].
Data()),
4431 Form(
"Cluster *asymmetry* in cell units vs E vs primary ID, for ID %s",matchCase[imatch].
Data()),
4432 nEbins,minE,maxE,10,-0.5,9.5,asyBins,asyMin,asyMax);
4436 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4437 fhDeltaIAOrigin[imatch]->GetYaxis()->SetBinLabel(ilabel,mcParticleStringLabel[ilabel-1]);
4441 (Form(
"hDeltaIATotOrigin_%s",matchCase[imatch].
Data()),
4442 Form(
"Cluster *total asymmetry* in cell units vs E vs primary ID, for ID %s",matchCase[imatch].
Data()),
4443 nEbins,minE,maxE,10,-0.5,9.5,asyBins,asyMin,asyMax);
4447 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4448 fhDeltaIATotOrigin[imatch]->GetYaxis()->SetBinLabel(ilabel,mcParticleStringLabel[ilabel-1]);
4454 (Form(
"hClusterMaxCellDiffM02_%s",matchCase[imatch].
Data()),
4455 Form(
"#it{E}_{cluster} vs (#it{E}_{cluster} - #it{E}_{cell max})/#it{E}_{cluster} vs #sigma^{2}_{long} for ID %s",matchCase[imatch].
Data()),
4456 nEbins,minE,maxE, 20,0,1.,nShShBins,minShSh,maxShSh);
4463 (Form(
"hClusterTimeEnergyM02_%s",matchCase[imatch].
Data()),
4464 Form(
"#it{E} vs TOF vs #sigma^{2}_{long} for ID %s",matchCase[imatch].
Data()),
4465 nEbins,minE,maxE,45,-25.5,20.5,nShShBins,minShSh,maxShSh);
4472 (Form(
"hNCellsPerClusterM02_%s",matchCase[imatch].
Data()),
4473 Form(
"#it{E} vs #it{n}_{cells} vs #sigma^{2}_{long} for ID %s",matchCase[imatch].
Data()),
4474 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4481 (Form(
"hNCellsPerClusterM20_%s",matchCase[imatch].
Data()),
4482 Form(
"#it{E} vs #it{n}_{cells} vs #sigma^{2}_{short} for ID %s",matchCase[imatch].
Data()),
4483 nEbins,minE,maxE,cellBins,cellMin,cellMax,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4490 (Form(
"hNCellsPerClusterNLM_%s",matchCase[imatch].
Data()),
4491 Form(
"#it{E} vs #it{n}_{cells} vs n_{lm} for ID %s",matchCase[imatch].
Data()),
4492 nEbins,minE,maxE,cellBins,cellMin,cellMax,10,0,10);
4499 (Form(
"hSMM02_%s",matchCase[imatch].
Data()),
4500 Form(
"#it{E} vs SM number vs #sigma^{2}_{long}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4501 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,minShSh,maxShSh);
4502 fhSMM02[imatch]->SetXTitle(
"#it{E} (GeV)");
4503 fhSMM02[imatch]->SetYTitle(
"SM number");
4504 fhSMM02[imatch]->SetZTitle(
"#sigma^{2}_{long}");
4505 outputContainer->Add(
fhSMM02[imatch]);
4508 (Form(
"hSMM02NoCut_%s",matchCase[imatch].
Data()),
4509 Form(
"#it{E} vs SM number vs #sigma^{2}_{long} for ID %s",matchCase[imatch].
Data()),
4510 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,minShSh,maxShSh);
4517 (Form(
"hSMM20LowM02_%s",matchCase[imatch].
Data()),
4518 Form(
"#it{E} vs SM number vs #sigma^{2}_{short}, #it{n}_{cells}^{w>0.01}>%d, for %2.2f<#sigma^{2}_{long}<%2.2f, for ID %s",
4520 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4527 (Form(
"hSMM20LowM02NoCut_%s",matchCase[imatch].
Data()),
4528 Form(
"#it{E} vs SM number vs #sigma^{2}_{short}, for %2.2f<#sigma^{2}_{long}<%2.2f, for ID %s",
4530 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4537 (Form(
"hSMM20HighM02_%s",matchCase[imatch].
Data()),
4538 Form(
"#it{E} vs SM number vs #sigma^{2}_{short}, #it{n}_{cells}^{w>0.01}>%d, for 0.5<#sigma^{2}_{long}<2, for ID %s",
4540 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4547 (Form(
"hSMM20HighM02NoCut_%s",matchCase[imatch].
Data()),
4548 Form(
"#it{E} vs SM number vs #sigma^{2}_{short}, for 0.5<#sigma^{2}_{long}<2, for ID %s",matchCase[imatch].
Data()),
4549 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4556 (Form(
"hSMNCell_%s",matchCase[imatch].
Data()),
4557 Form(
"#it{E} vs SM number vs #it{n}_{cells}^{w>0.01} for ID %s",matchCase[imatch].
Data()),
4558 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,cellBins,cellMin,cellMax);
4559 fhSMNCell[imatch]->SetXTitle(
"#it{E} (GeV)");
4560 fhSMNCell[imatch]->SetYTitle(
"SM number");
4561 fhSMNCell[imatch]->SetZTitle(
"#it{n}_{cells}^{w>0.01}");
4562 outputContainer->Add(
fhSMNCell[imatch]);
4565 (Form(
"hSMNCellM02_%s",matchCase[imatch].
Data()),
4566 Form(
"SM number vs #it{n}_{cells}^{w>0.01} vs #sigma^{2}_{long}, " 4568 fNModules,-0.5,fNModules-0.5,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4571 fhSMNCellM02[imatch]->SetYTitle(
"#it{n}_{cells}^{w>0.01}");
4575 (Form(
"hColM02_%s",matchCase[imatch].
Data()),
4576 Form(
"#it{E} vs column number vs #sigma^{2}_{long}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4577 nEbins,minE,maxE,48,-0.5,47.5,nShShBins,minShSh,maxShSh);
4578 fhColM02[imatch]->SetXTitle(
"#it{E} (GeV)");
4579 fhColM02[imatch]->SetYTitle(
"column number");
4580 fhColM02[imatch]->SetZTitle(
"#sigma^{2}_{long}");
4581 outputContainer->Add(
fhColM02[imatch]);
4584 (Form(
"hRowM02_%s",matchCase[imatch].
Data()),
4585 Form(
"#it{E} vs row number vs #sigma^{2}_{long}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4586 nEbins,minE,maxE,24,-0.5,23.5,nShShBins,minShSh,maxShSh);
4587 fhRowM02[imatch]->SetXTitle(
"#it{E} (GeV)");
4588 fhRowM02[imatch]->SetYTitle(
"row number");
4589 fhRowM02[imatch]->SetZTitle(
"#sigma^{2}_{long}");
4590 outputContainer->Add(
fhRowM02[imatch]);
4595 (Form(
"hNCellsPerClusterMEta_%s",matchCase[imatch].
Data()),
4596 Form(
"#it{E} vs #it{n}_{cells} vs #sigma_{#eta}^{2} for ID %s",matchCase[imatch].
Data()),
4597 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4604 (Form(
"hNCellsPerClusterMPhi_%s",matchCase[imatch].
Data()),
4605 Form(
"#it{E} vs #it{n}_{cells} vs #sigma_{#varphi}^{2} for ID %s",matchCase[imatch].
Data()),
4606 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4613 (Form(
"hNCellsPerClusterMEtaPhi_%s",matchCase[imatch].
Data()),
4614 Form(
"#it{E} vs #it{n}_{cells} vs #sigma_{#eta#varphi}^{2} for ID %s",matchCase[imatch].
Data()),
4615 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,-1*maxShSh,maxShSh);