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), fStudyModuleCells(kTRUE),
50 fFillNLocMaxHistos(kTRUE),
53 fM02Min(0), fNCellMin(0),
54 fMinDistToBad(0), fNEBinCuts(0),
55 fEMinShape(0), fEMaxShape(100),
56 fNCellMinShape(-1), fNCellsBins(3),
59 fdEdXMinEle(0), fdEdXMaxEle(0),
60 fdEdXMinHad(0), fdEdXMaxHad(0),
63 fInvMassMinECut(0), fInvMassMaxECut(0),
64 fInvMassMinM02Cut(0), fInvMassMaxM02Cut(0),
65 fInvMassMaxOpenAngle(0), fInvMassMaxTimeDifference(0),
67 fConstantTimeShift(0),
69 fClusterMomentum(), fClusterMomentum2(),
70 fCaloCellList(NULL), fCaloClusList(NULL),
75 fhEnergyTMEtaResidual1Cell(0), fhEnergyTMPhiResidual1Cell(0),
76 fhColRowExoticHighE1CellPosTime(0), fhColRowExoticHighE1CellNegTime(0), fhColRowExoticHighE1CellNulTime(0),
77 fhEnergyTMEtaResidualExotic(0), fhEnergyTMPhiResidualExotic(0),
78 fhColRowExoticHighEPosTime(0), fhColRowExoticHighENegTime(0), fhColRowExoticHighENulTime(0),
79 fhColRowHighEPosTime(0), fhColRowHighENegTime(0), fhColRowHighENulTime(0),
80 fhEnergyTMEtaResidualTCardCorrNoSelection1Cell(0), fhEnergyTMPhiResidualTCardCorrNoSelection1Cell(0),
81 fhEnergyTMEtaResidualTCardCorrNoSelectionExotic(0), fhEnergyTMPhiResidualTCardCorrNoSelectionExotic(0),
86 fhClusterMaxCellCloseCellRatioM02(0), fhClusterMaxCellECrossM02(0),
87 fhInvMassNCellSM(0), fhInvMassNCellSMSame(0),
88 fhColRowM02(0), fhColRowM02NCellCut(0),
89 fhEMaxCellTimeM02SM(0), fhEMaxCellTimeNCellSM(0), fhESecCellTimeNCellSM(0),
91 fhSMNCellModuleMax(0), fhSMNCellModuleOut(0),
92 fhSMNCellModuleOutModSame(0), fhSMNCellModuleOutModDiff(0),
93 fhSMECellModuleMax(0), fhSMECellModuleMaxTot(0),
94 fhSMECellModuleOut(0), fhSMECellModuleOutModSame(0), fhSMECellModuleOutModDiff(0),
95 fhSMNCellModuleMaxOutRat(0), fhSMNCellModuleMaxOutModSameRat(0), fhSMNCellModuleMaxOutModDiffRat(0),
96 fhSMECellModuleMaxRat(0), fhSMECellModuleMaxTotRat(0), fhSMECellModuleMaxTotRatClus(0),
97 fhSMECellModuleMaxOutRat(0), fhSMECellModuleMaxOutModSameRat(0), fhSMECellModuleMaxOutModDiffRat(0),
99 fhSMNCellModuleMaxLowM02(0), fhSMNCellModuleOutLowM02(0),
100 fhSMNCellModuleOutModSameLowM02(0), fhSMNCellModuleOutModDiffLowM02(0),
101 fhSMECellModuleMaxLowM02(0), fhSMECellModuleMaxTotLowM02(0),
102 fhSMECellModuleOutLowM02(0), fhSMECellModuleOutModSameLowM02(0), fhSMECellModuleOutModDiffLowM02(0),
103 fhSMNCellModuleMaxOutRatLowM02(0), fhSMNCellModuleMaxOutModSameRatLowM02(0),fhSMNCellModuleMaxOutModDiffRatLowM02(0),
104 fhSMECellModuleMaxRatLowM02(0), fhSMECellModuleMaxTotRatLowM02(0), fhSMECellModuleMaxTotRatClusLowM02(0),
105 fhSMECellModuleMaxOutRatLowM02(0), fhSMECellModuleMaxOutModSameRatLowM02(0),fhSMECellModuleMaxOutModDiffRatLowM02(0),
107 fhSMNCellModuleMaxHighM02(0), fhSMNCellModuleOutHighM02(0),
108 fhSMNCellModuleOutModSameHighM02(0), fhSMNCellModuleOutModDiffHighM02(0),
109 fhSMECellModuleMaxHighM02(0), fhSMECellModuleMaxTotHighM02(0),
110 fhSMECellModuleOutHighM02(0), fhSMECellModuleOutModSameHighM02(0), fhSMECellModuleOutModDiffHighM02(0),
111 fhSMNCellModuleMaxOutRatHighM02(0), fhSMNCellModuleMaxOutModSameRatHighM02(0),fhSMNCellModuleMaxOutModDiffRatHighM02(0),
112 fhSMECellModuleMaxRatHighM02(0), fhSMECellModuleMaxTotRatHighM02(0), fhSMECellModuleMaxTotRatClusHighM02(0),
113 fhSMECellModuleMaxOutRatHighM02(0), fhSMECellModuleMaxOutModSameRatHighM02(0),fhSMECellModuleMaxOutModDiffRatHighM02(0),
115 fhSMEMaxEClusterRat(0),
116 fhSMEMaxEClusterRatLowM02(0), fhSMEMaxEClusterRatHighM02(0),
118 fhTCardChannelM02(0), fhTCardChannelM02NoCut(0),
119 fhTCardChannelNCellModMax(0), fhTCardChannelNCell(0),
122 fhECellClusterRatio(0), fhECellClusterLogRatio(0),
123 fhEMaxCellClusterRatio(0), fhEMaxCellClusterLogRatio(0),
124 fhECellTotalRatio(0), fhECellTotalLogRatio(0),
125 fhECellTotalRatioMod(0), fhECellTotalLogRatioMod(0)
127 for(
Int_t i=0; i < 3; i++)
192 for(
Int_t i =0; i < 12; i++)
194 for(
Int_t j = 0; j < 4; j++)
196 for(
Int_t k = 0; k < 3; k++)
203 for(
Int_t j = 0; j < 5; j++)
211 for(
Int_t tm = 0; tm < 2; tm++)
307 for(
Int_t i = 0; i < 6; i++)
309 for(
Int_t j = 0; j < 6; j++)
382 for(
Int_t i = 0; i < 14; i++)
401 for(
Int_t j = 0; j < 6; j++)
403 for(
Int_t k = 0; k < 6; k++)
424 for(
Int_t i = 0; i < 12; i++)
450 for(
Int_t i = 0; i < 20; i++)
479 for(
Int_t j = 0; j < 4; j++)
523 Int_t icol = -1, irow = -1, iRCU = -1, icolAbs = -1, irowAbs = -1;
533 Int_t ncells = clus->GetNCells();
534 Double_t time = clus->GetTOF()*1.e9;
537 Float_t deta = clus->GetTrackDz();
538 Float_t dphi = clus->GetTrackDx();
546 if(energy >= 5 && energy < 8)
561 else if(exoticity > 0.97)
563 if(energy >= 5 && energy < 8)
578 for (
Int_t ipos = 0; ipos < ncells; ipos++)
580 Int_t absId = clus->GetCellsAbsId()[ipos];
588 if( absId == absIdMax || weight < 0.01 )
continue;
590 Int_t rowDiff = -100, colDiff = -100;
594 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
599 if(energy >= 5 && energy < 8)
606 if(energy >= 5 && energy < 8)
614 else if ( energy > 8 )
627 else if ( energy > 8 && ncells > 1 && exoticity < 0.97)
637 Int_t etaRegion = -1, phiRegion = -1;
640 if ( etaRegion !=0 ) return ;
650 else if(exoticity > 0.97)
661 if ( m02 > 0.001 ) lamR = m20/m02;
663 Int_t absIdList[ncells];
671 Int_t nCellWithWeight = 1;
703 Bool_t sameTCard2ndMax = kFALSE;
707 Float_t eCellSameRowSameTCardNearCol = 0.;
708 Float_t eCellSameRowDiffTCardNearCol = 0.;
709 Double_t tCellSameRowSameTCardNearCol = 0.;
710 Double_t tCellSameRowDiffTCardNearCol = 0.;
717 for (
Int_t ipos = 0; ipos < ncells; ipos++)
719 Int_t absId = clus->GetCellsAbsId()[ipos];
736 if( absId == absIdMax || weight < 0.01 )
continue;
760 Int_t rowDiff = -100, colDiff = -100;
791 Int_t indexType = -1;
796 if(TMath::Abs(rowDiff) == 1) nearRow = kTRUE;
797 if(TMath::Abs(colDiff) == 1) nearCol = kTRUE;
799 if ( rowDiff == 0 && colDiff != 0 )
801 if ( nearCol ) indexType = 6;
807 else if ( rowDiff != 0 && colDiff == 0 )
809 if ( nearRow ) indexType = 8;
817 if ( nearRow && nearCol) indexType = 10;
828 if ( rowDiff == 0 && colDiff != 0 )
830 if ( nearCol ) indexType = 0;
834 else if ( rowDiff != 0 && colDiff == 0 )
836 if ( nearRow ) indexType = 2;
841 if ( nearCol && nearRow ) indexType = 4;
846 if ( rowDiff == 0 && TMath::Abs(colDiff) == 1 )
850 eCellSameRowSameTCardNearCol = eCell;
851 tCellSameRowSameTCardNearCol = tCell;
855 eCellSameRowDiffTCardNearCol = eCell;
856 tCellSameRowDiffTCardNearCol = tCell;
862 Float_t eCellDiff = eCellMax - eCell;
863 Float_t eClusDiff = energy - eCell;
866 Float_t tCellDiff = tCellMax - tCell;
887 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
892 if(energy >= 5 && energy < 8)
899 if(energy >= 5 && energy < 8)
906 if ( eCell > emax2nd )
912 sameTCard2ndMax = kTRUE;
918 sameTCard2ndMax = kFALSE;
953 Float_t ratioNcells = nCellWithWeight/(ncells*1.);
958 if ( nCellWithWeight <= 1 )
return;
963 printf(
"AliAnaClusterShapeCorrelStudies: M02 %f, M20 %f, E %2.3f, ncell %d, n with weight %d; max cell E %2.3f\n",
964 m02,m20,energy,ncells,nCellWithWeight,eCellMax);
973 if(energy >= 5 && energy < 8)
980 if(energy >= 5 && energy < 8)
1022 if ( sameTCard2ndMax )
1045 if(eCellSameRowSameTCardNearCol > 0 && eCellSameRowDiffTCardNearCol > 0)
1047 Float_t eDiff = eCellSameRowSameTCardNearCol - eCellSameRowDiffTCardNearCol ;
1048 Float_t tDiff = tCellSameRowSameTCardNearCol - tCellSameRowDiffTCardNearCol ;
1065 Int_t nCorrInd = nCorr;
1066 if(nCorr > 4) nCorrInd = 5;
1068 Int_t nCorrNoInd = nCorrNo;
1069 if(nCorrNoInd > 4) nCorrNoInd = 5;
1118 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1120 Int_t absId = clus->GetCellsAbsId()[ipos];
1139 if( absId == absIdMax || weight < 0.01 )
continue;
1141 Float_t tDiffMaxSecondary = tCellMax - tCell;
1148 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
1151 Int_t rowDiff = -100, colDiff = -100;
1156 if ( energy >= 5 && energy < 8 && nCorr == 0 )
1158 else if ( energy >= 8 && nCorr == 0 )
1163 if ( energy >= 5 && energy < 8 && nCorrNo == 0 )
1165 else if ( energy >= 8 && nCorrNo == 0 )
1199 AliVCluster* clus2 = (AliVCluster*)
fCaloClusList->At(jclus);
1204 Double_t tof2 = clus2->GetTOF()*1.e9;
1212 if( absIdMax == absIdMax2
1213 || !
IsGoodCluster(absIdMax2, clus2->GetM02(), clus2->GetNCells())
1245 if ( energy >= 5 && energy < 8)
1247 else if ( energy >= 8 )
1265 if ( energy >= 5 && energy < 8)
1267 else if ( energy >= 8 )
1333 if ( exoticity > 0.97 )
1348 Float_t nCellRat = nCorr*1. / ((nCorr+nCorrNo)*1.);
1450 Int_t nCaloCellsPerCluster = clus->GetNCells();
1451 if ( nCaloCellsPerCluster < 2 )
return;
1468 for (
Int_t ipos = 0; ipos < nCaloCellsPerCluster; ipos++)
1470 Int_t absId = clus ->GetCellsAbsId()[ipos];
1475 if( absId == absIdMax || eCell < 0.01 )
continue;
1479 if( weight < 0.01 )
continue;
1484 if ( nCell < 1 )
return;
1487 Int_t nCellBin = -1;
1502 Int_t ietaMax=-1, iphiMax = 0, rcuMax = 0, icolAbs = -1, irowAbs = -1;
1504 ietaMax, iphiMax, rcuMax, icolAbs, irowAbs);
1506 Int_t iTowerMax = -1, iIphiMax = -1, iIetaMax = -1;
1508 GetEMCALGeometry()->GetCellIndex(absIdMax,smMax,iTowerMax,iIphiMax,iIetaMax);
1523 Int_t etaRegion = -1, phiRegion = -1;
1526 if ( etaRegion !=0 ) return ;
1530 Int_t nPerModule = 0;
1531 Int_t nPerModuleOut = 0;
1532 Int_t nPerModuleOutD = 0;
1533 Int_t nPerModuleOutS = 0;
1538 for (
Int_t ipos = 0; ipos < nCaloCellsPerCluster; ipos++)
1540 Int_t absId = clus ->GetCellsAbsId()[ipos];
1545 if( absId == absIdMax || eCell < 0.01 )
continue;
1554 Int_t iTower = -1, iIphi = -1, iIeta = -1;
1559 if ( iTowerMax == iTower )
1569 ePerModuleOut+=eCell;
1571 if(TMath::Abs(iphiMax-iphi) == 1 || TMath::Abs(ietaMax-ieta) == 1)
1573 Int_t rowDiff = -100, colDiff = -100;
1578 ePerModuleOutS+=eCell;
1583 ePerModuleOutD+=eCell;
1593 if ( weight < 0.01 )
continue;
1595 Float_t fracCell = eCell/eCellMax;
1598 Float_t fracClus = (energy-eCell)/energy;
1644 if(dIphi < TMath::Abs(iphi-iphiMax)) dIphi = TMath::Abs(iphi-iphiMax);
1645 if(iphi-iphiMax < 0 && dIphiNeg > iphi-iphiMax) dIphiNeg = iphi-iphiMax;
1646 if(iphi-iphiMax > 0 && dIphiPos < iphi-iphiMax) dIphiPos = iphi-iphiMax;
1648 if(dIeta < TMath::Abs(ieta-ietaMax)) dIeta = TMath::Abs(ieta-ietaMax);
1649 if(ieta-ietaMax < 0 && dIetaNeg > ieta-ietaMax) dIetaNeg = ieta-ietaMax;
1650 if(ieta-ietaMax > 0 && dIetaPos < ieta-ietaMax) dIetaPos = ieta-ietaMax;
1654 Int_t ietaShift = ieta;
1655 Int_t ietaMaxShift = ietaMax;
1657 if (ieta > ietaMax) ietaMaxShift+=48;
1658 else ietaShift +=48;
1660 if(dIeta < TMath::Abs(ietaShift-ietaMaxShift)) dIeta = TMath::Abs(ietaShift-ietaMaxShift);
1661 if(ietaShift-ietaMaxShift < 0 && dIetaNeg > ietaShift-ietaMaxShift) dIetaNeg = ietaShift-ietaMaxShift;
1662 if(ietaShift-ietaMaxShift > 0 && dIetaPos < ietaShift-ietaMaxShift) dIetaPos = ietaShift-ietaMaxShift;
1675 if ( nCellBin >= 0 && nCellBin < fNCellsBins )
1687 if ( nCellBin >= 0 && nCellBin < fNCellsBins )
1704 Int_t iTCard = rowTCard+8*colTCard;
1712 Float_t ePerModuleTot = ePerModule+eCellMax;
1857 if ( matchedPID == 0 )
1891 Float_t dIA = 1.*(dIphi-dIeta)/(dIeta+dIphi);
1892 Float_t dIATot = 1.*((dIphiPos-dIphiNeg)-(dIetaPos-dIetaNeg))/((dIetaPos-dIetaNeg)+(dIphiPos-dIphiNeg));
1901 AliDebug(1,Form(
"E %2.2f, nCell %d, dPhi %d, dEta %d, dIA %2.2f, match %d",energy,nCell, dIphi,dIeta,dIA,matchedPID));
1934 Float_t dispp= 0., dEta = 0., dPhi = 0.;
1935 Float_t sEta = 0., sPhi = 0., sEtaPhi = 0.;
1939 l0, l1, dispp, dEta, dPhi, sEta, sPhi, sEtaPhi);
1942 if(sEta+sPhi>0.0001) sEtaPhiA = (sPhi-sEta)/(sEta+sPhi);
1944 AliDebug(2,Form(
"Recalculate shower shape org: m02 %2.2f, m20 %2.2f, disp %2.2f;" 1945 " new: m02 %2.2f, m20 %2.2f, disp %2.2f; " 1946 "mEta %2.2f, mPhi %2.2f, mEtaPhi %2.2f, A_EtaPhi %2.2f; dEta %2.2f dPhi %2.2f",
1947 m02,m20,clus->GetDispersion(),l0,l1,dispp,sEta,sPhi,sEtaPhi,sEtaPhiA,dEta,dPhi));
1960 if ( nCell > fNCellMinShape )
1979 AliVCluster* clus2 = (AliVCluster*)
fCaloClusList->At(jclus);
1984 Double_t tof2 = clus2->GetTOF()*1.e9;
1992 if( absIdMax == absIdMax2
1993 || !
IsGoodCluster(absIdMax2, clus2->GetM02(), clus2->GetNCells())
2030 if (
IsDataMC() && mcIndex > -1 && mcIndex < 10)
2072 AliDebug(1,Form(
"No track found with matches %d",clus->GetNTracksMatched()));
2078 Double_t dedx = track->GetTPCsignal();
2098 else matchedPID =-1;
2100 AliDebug(1,Form(
"matches %d; dEdX %2.2f; matched PID %d",clus->GetNTracksMatched(),dedx,matchedPID));
2116 Int_t nCaloCellsPerCluster = 0 ;
2127 for(
Int_t iclus = 0; iclus < nCaloClusters; iclus++)
2129 AliDebug(1,Form(
"Cluster: %d/%d, data %d",iclus+1,nCaloClusters,
GetReader()->GetDataType()));
2131 AliVCluster* clus = (AliVCluster*)
fCaloClusList->At(iclus);
2134 if ( clus->GetDistanceToBadChannel() > -1 &&
2137 AliDebug(1,Form(
"Small distance to bad channel %2.2f < %2.2f",clus->GetDistanceToBadChannel(),
fMinDistToBad));
2143 if ( nModule < fFirstModule || nModule >
fLastModule )
2145 AliDebug(1,Form(
"Cluster module out of range %d",nModule));
2155 Double_t tof = clus->GetTOF()*1.e9;
2173 AliDebug(1,Form(
"Remove cluster with phi %2.2f and eta %2.2f", phi*TMath::RadToDeg(), eta));
2177 AliDebug(1,Form(
"cluster: E %2.3f, pT %2.3f, eta %2.3f, phi %2.3f",e,pt,eta,phi*TMath::RadToDeg()));
2181 nCaloCellsPerCluster = clus->GetNCells();
2185 AliDebug(1,Form(
"Accept cluster? %d",goodCluster));
2187 if(!goodCluster)
continue;
2213 const UInt_t nlabels = clus->GetNLabels();
2214 Int_t overpdg[nlabels];
2215 Int_t overlab[nlabels];
2217 GetMC(), overpdg, overlab);
2224 Int_t matchedPID = -1;
2247 if ( ampMax > 0.01 )
2255 if (
fStudyShape && matchedPID >= 0 && matchedPID < 3 )
2256 ClusterShapeHistograms(clus, absIdMax, maxCellFraction, eCrossFrac, ampMax, tmax, matchedPID, mcIndex, noverlaps);
2276 const Int_t buffersize = 255;
2277 char onePar[buffersize] ;
2279 snprintf(onePar,buffersize,
"--- AliAnaClusterShapeCorrelStudies ---:") ;
2300 return new TObjString(parList) ;
2310 outputContainer->SetName(
"ClusterShapeStudies") ;
2321 TString mcParticleStringLabel[] = {
"Merged #gamma#gamma",
"Decay #gamma",
"Direct #gamma",
"e^{#pm}",
"#pi^{#pm}",
"k^{#pm}",
"p",
"#bar{p}",
"n",
"#bar{n}"};
2325 TString matchCase[] = {
"Neutral",
"Electron",
"Hadron"};
2353 Float_t rowcellmin = fNMaxRowsFullMin-1.5;
2362 Int_t nShShBins = 200;
2372 Int_t cellBins = 15;
2382 (
"hColRowHighEPosTime",
2383 "column vs row, exo < 0.97, E > 8 GeV, t > 5 ns",
2384 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2390 (
"hColRowHighENegTime",
2391 "column vs row, exo < 0.97, E > 8 GeV, t < -5 ns",
2392 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2398 (
"hColRowHighENulTime",
2399 "column vs row, exo < 0.97, E > 8 GeV, -5 < t < 5 ns",
2400 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2408 (
"hColRowExoticHighE1CellPosTime",
2409 "column vs row, 1 cell, E > 8 GeV, t > 5 ns",
2410 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2416 (
"hColRowExoticHighEPosTime",
2417 "column vs row, exo > 0.97, E > 8 GeV, t > 5 ns",
2418 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2424 (
"hColRowExoticHighE1CellNegTime",
2425 "column vs row, 1 cell, E > 8 GeV, t < -5 ns",
2426 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2432 (
"hColRowExoticHighENegTime",
2433 "column vs row, exo > 0.97, E > 8 GeV, t < -5 ns",
2434 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2440 (
"hColRowExoticHighE1CellNulTime",
2441 "column vs row, 1 cell, E > 8 GeV, -5 < t < 5 ns",
2442 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2448 (
"hColRowExoticHighENulTime",
2449 "column vs row, exo > 0.97, E > 8 GeV, -5 < t < 5 ns",
2450 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2456 TString add[] = {
"",
"TrackMatched"};
2457 for(
Int_t tm = 0; tm < 2; tm++)
2460 (Form(
"hColRowTCardCorrNoSelectionLowE%s",add[tm].
Data()),
2461 Form(
"column vs row, max E cell for TCard correlation selected clusters, 5 < E < 8 GeV %s",add[tm].
Data()),
2462 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2468 (Form(
"hColRowTCardCorrNoSelectionHighE%s",add[tm].
Data()),
2469 Form(
"column vs row, max E cell for TCard correlation selected clusters, E > 8 GeV %s",add[tm].
Data()),
2470 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2478 (Form(
"hNCellsTCardCorrNoSelection%s",add[tm].
Data()),
2479 Form(
"# custer # cells vs #it{E} %s",add[tm].
Data()),
2480 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
2486 (Form(
"hNCellsTCardCorrWithWeightNoSelection%s",add[tm].
Data()),
2487 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
2488 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
2494 (Form(
"hNCellsTCardCorrRatioWithWeightNoSelection%s",add[tm].
Data()),
2495 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
2496 nptbins,ptmin,ptmax, 100,0,1);
2502 (Form(
"hTimeTCardCorrNoSelection%s",add[tm].
Data()),
2503 Form(
"#it{time} vs #it{E} %s",add[tm].
Data()),
2504 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2510 (Form(
"hLambda0TCardCorrNoSelection%s",add[tm].
Data()),
2511 Form(
"#sigma^{2}_{long} vs #it{E} %s",add[tm].
Data()),
2512 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2518 (Form(
"hLambda1TCardCorrNoSelection%s",add[tm].
Data()),
2519 Form(
"#sigma^{2}_{short} vs #it{E} %s",add[tm].
Data()),
2520 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2528 (Form(
"hLambda0NLM1TCardCorrNoSelection%s",add[tm].
Data()),
2529 Form(
"#sigma^{2}_{long} vs #it{E}, nlm=1 %s",add[tm].
Data()),
2530 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2536 (Form(
"hLambda1NLM1TCardCorrNoSelection%s",add[tm].
Data()),
2537 Form(
"#sigma^{2}_{short} vs #it{E}, nlm=1 %s",add[tm].
Data()),
2538 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2544 (Form(
"hLambda0NLM2TCardCorrNoSelection%s",add[tm].
Data()),
2545 Form(
"#sigma^{2}_{long} vs #it{E}, nlm=2 %s",add[tm].
Data()),
2546 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2552 (Form(
"hLambda1NLM2TCardCorrNoSelection%s",add[tm].
Data()),
2553 Form(
"#sigma^{2}_{short} vs #it{E}, nlm=2 %s",add[tm].
Data()),
2554 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2561 (Form(
"hLambdaRTCardCorrNoSelection%s",add[tm].
Data()),
2562 Form(
"#sigma^{2}_{short}/#sigma^{2}_{long} vs #it{E} %s",add[tm].
Data()),
2563 nptbins,ptmin,ptmax,110,0,1.1);
2569 (Form(
"hNLocMaxTCardCorrNoSelection%s",add[tm].
Data()),
2570 Form(
"#it{n}_{LM} vs E %s",add[tm].
Data()),
2571 nptbins,ptmin,ptmax,10,0,10);
2579 (Form(
"hEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2580 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
2581 nptbins,ptmin,ptmax,100,0,1);
2587 (Form(
"hEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2588 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2589 nptbins,ptmin,ptmax,100,0,1);
2595 (Form(
"hEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2596 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2597 nptbins,ptmin,ptmax,100,0,1);
2603 (Form(
"hE2ndRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2604 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
2605 nptbins,ptmin,ptmax,100,0,1);
2611 (Form(
"hE2ndRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2612 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2613 nptbins,ptmin,ptmax,100,0,1);
2619 (Form(
"hE2ndRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2620 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2621 nptbins,ptmin,ptmax,100,0,1);
2628 (Form(
"hE2ndEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2629 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
2630 nptbins,ptmin,ptmax,100,0,1);
2636 (Form(
"hE2ndEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2637 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2638 nptbins,ptmin,ptmax,100,0,1);
2644 (Form(
"hE2ndEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2645 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2646 nptbins,ptmin,ptmax,100,0,1);
2654 (Form(
"hE2ndSameRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2655 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()),
2656 nptbins,ptmin,ptmax,100,0,1);
2662 (Form(
"hE2ndSameRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2663 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()),
2664 nptbins,ptmin,ptmax,100,0,1);
2670 (Form(
"hE2ndSameRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2671 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()),
2672 nptbins,ptmin,ptmax,100,0,1);
2679 (Form(
"hE2ndSameEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2680 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()),
2681 nptbins,ptmin,ptmax,100,0,1);
2687 (Form(
"hE2ndSameEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2688 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()),
2689 nptbins,ptmin,ptmax,100,0,1);
2695 (Form(
"hE2ndSameEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2696 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()),
2697 nptbins,ptmin,ptmax,100,0,1);
2705 (Form(
"hECellClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2706 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
2707 nptbins,ptmin,ptmax,100,0,1);
2713 (Form(
"hECellClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2714 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2715 nptbins,ptmin,ptmax,100,0,1);
2721 (Form(
"hECellClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2722 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2723 nptbins,ptmin,ptmax,100,0,1);
2729 (Form(
"hLogECellNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2730 Form(
"log(#it{E}_{cell}) vs E_{cluster}, #it{n}_{LM}=1, w > 0.01 %s",add[tm].
Data()),
2731 nptbins,ptmin,ptmax,150,-3,3);
2737 (Form(
"hLogECellNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2738 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}=2, w > 0.01 %s",add[tm].
Data()),
2739 nptbins,ptmin,ptmax,150,-3,3);
2745 (Form(
"hLogECellNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2746 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}>2, w > 0.01 %s",add[tm].
Data()),
2747 nptbins,ptmin,ptmax,150,-3,3);
2754 (Form(
"hECellWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2755 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
2756 nptbins,ptmin,ptmax,90,0,4.5);
2762 (Form(
"hECellWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2763 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2764 nptbins,ptmin,ptmax,90,0,4.5);
2770 (Form(
"hECellWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2771 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster}))vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2772 nptbins,ptmin,ptmax,90,0,4.5);
2779 (Form(
"hECellSameClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2780 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()),
2781 nptbins,ptmin,ptmax,100,0,1);
2787 (Form(
"hECellSameClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2788 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, cell from same T-Card as leading, #it{n}_{LM}=2 %s",add[tm].
Data()),
2789 nptbins,ptmin,ptmax,100,0,1);
2795 (Form(
"hECellSameClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2796 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, cell from same T-Card as leading, #it{n}_{LM}>2 %s",add[tm].
Data()),
2797 nptbins,ptmin,ptmax,100,0,1);
2803 (Form(
"hLogECellSameNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2804 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()),
2805 nptbins,ptmin,ptmax,150,-3,3);
2811 (Form(
"hLogECellSameNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2812 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()),
2813 nptbins,ptmin,ptmax,150,-3,3);
2819 (Form(
"hLogECellSameNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2820 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()),
2821 nptbins,ptmin,ptmax,150,-3,3);
2828 (Form(
"hECellSameWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2829 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()),
2830 nptbins,ptmin,ptmax,90,0,4.5);
2836 (Form(
"hECellSameWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2837 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()),
2838 nptbins,ptmin,ptmax,90,0,4.5);
2844 (Form(
"hECellSameWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2845 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()),
2846 nptbins,ptmin,ptmax,90,0,4.5);
2853 (Form(
"hExoticTCardCorrNoSelection%s",add[tm].
Data()),
2854 Form(
"exoticity vs #it{E} %s",add[tm].
Data()),
2855 nptbins,ptmin,ptmax,200,-1,1);
2861 (Form(
"hNCellsTCardSameAndDiffFraction%s",add[tm].
Data()),
2862 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, %s",add[tm].
Data()),
2863 nptbins,ptmin,ptmax,100,0,1);
2869 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterE%s",add[tm].
Data()),
2870 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()),
2871 nptbins,ptmin,ptmax,200,-10,10);
2877 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterE%s",add[tm].
Data()),
2878 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()),
2879 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2885 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxE%s",add[tm].
Data()),
2886 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()),
2887 nptbins,ptmin,ptmax,200,-10,10);
2893 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxE%s",add[tm].
Data()),
2894 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()),
2895 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2903 (Form(
"hEnergyTime1Cell%s",add[tm].
Data()),
2904 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
2905 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2911 (Form(
"hEnergyTimeExotic%s",add[tm].
Data()),
2912 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
2913 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2919 (Form(
"hEnergyTimeTCardCorrNoSelection1Cell%s",add[tm].
Data()),
2920 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
2921 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2927 (Form(
"hEnergyTimeTCardCorrNoSelectionExotic%s",add[tm].
Data()),
2928 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
2929 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2935 (Form(
"hColRowExoticLowE1Cell%s",add[tm].
Data()),
2936 Form(
"column vs row, 1 cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2937 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2943 (Form(
"hColRowExoticHighE1Cell%s",add[tm].
Data()),
2944 Form(
"column vs row, 1 cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2945 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2951 (Form(
"hColRowExoticLowE%s",add[tm].
Data()),
2952 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2953 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2959 (Form(
"hColRowExoticHighE%s",add[tm].
Data()),
2960 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2961 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2967 (Form(
"hColRowExotic2ndCellDiffLowE%s",add[tm].
Data()),
2968 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2969 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2975 (Form(
"hColRowExotic2ndCellDiffHighE%s",add[tm].
Data()),
2976 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2977 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2983 (Form(
"hColRowExotic2ndCellSameLowE%s",add[tm].
Data()),
2984 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2985 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2991 (Form(
"hColRowExotic2ndCellSameHighE%s",add[tm].
Data()),
2992 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2993 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2999 (Form(
"hColRowTCardCorrNoSelectionExoticLowE%s",add[tm].
Data()),
3000 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3001 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3007 (Form(
"hColRowTCardCorrNoSelectionExoticHighE%s",add[tm].
Data()),
3008 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3009 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3015 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffLowE%s",add[tm].
Data()),
3016 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3017 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3023 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffHighE%s",add[tm].
Data()),
3024 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3025 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3031 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameLowE%s",add[tm].
Data()),
3032 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3033 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3039 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameHighE%s",add[tm].
Data()),
3040 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3041 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3047 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameLowE%s",add[tm].
Data()),
3048 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3049 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3055 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameHighE%s",add[tm].
Data()),
3056 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3057 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3063 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffLowE%s",add[tm].
Data()),
3064 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
3065 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3071 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffHighE%s",add[tm].
Data()),
3072 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
3073 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3079 (Form(
"hNCellsTCardSameAndDiffFraction_Exotic%s",add[tm].
Data()),
3080 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, exo > 0.97 %s",add[tm].
Data()),
3081 nptbins,ptmin,ptmax,100,0,1);
3087 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterEExo%s",add[tm].
Data()),
3088 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()),
3089 nptbins,ptmin,ptmax,200,-10,10);
3095 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterEExo%s",add[tm].
Data()),
3096 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()),
3097 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
3103 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxEExo%s",add[tm].
Data()),
3104 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()),
3105 nptbins,ptmin,ptmax,200,-10,10);
3111 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxEExo%s",add[tm].
Data()),
3112 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()),
3113 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
3119 for(
Int_t i = 0; i < 6; i++)
3121 for(
Int_t j = 0; j < 6; j++)
3124 (Form(
"hLambda0TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3125 Form(
"#sigma^{2}_{long} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3126 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3132 (Form(
"hLambda1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3133 Form(
"#sigma^{2}_{short} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3134 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3150 (Form(
"hLambda0NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3151 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()),
3152 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3158 (Form(
"hLambda1NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3159 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()),
3160 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3166 (Form(
"hLambda0NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3167 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()),
3168 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3174 (Form(
"hLambda1NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3175 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()),
3176 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3182 (Form(
"hNLocMaxTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3183 Form(
"#it{n}_{LM} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3184 nptbins,ptmin,ptmax,10,0,10);
3190 (Form(
"hEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3191 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()),
3192 nptbins,ptmin,ptmax,100,0,1);
3198 (Form(
"hEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3199 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()),
3200 nptbins,ptmin,ptmax,100,0,1);
3206 (Form(
"hEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3207 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()),
3208 nptbins,ptmin,ptmax,100,0,1);
3214 (Form(
"hE2ndRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3215 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()),
3216 nptbins,ptmin,ptmax,100,0,1);
3222 (Form(
"hE2ndRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3223 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()),
3224 nptbins,ptmin,ptmax,100,0,1);
3230 (Form(
"hE2ndRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3231 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()),
3232 nptbins,ptmin,ptmax,100,0,1);
3238 (Form(
"hE2ndEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3239 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()),
3240 nptbins,ptmin,ptmax,100,0,1);
3246 (Form(
"hE2ndEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3247 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()),
3248 nptbins,ptmin,ptmax,100,0,1);
3254 (Form(
"hE2ndEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3255 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()),
3256 nptbins,ptmin,ptmax,100,0,1);
3263 (Form(
"hECellClusRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3264 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()),
3265 nptbins,ptmin,ptmax,100,0,1);
3271 (Form(
"hECellClusRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3272 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()),
3273 nptbins,ptmin,ptmax,100,0,1);
3279 (Form(
"hECellClusRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3280 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()),
3281 nptbins,ptmin,ptmax,100,0,1);
3287 (Form(
"hLogECellNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3288 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()),
3289 nptbins,ptmin,ptmax,150,-3,3);
3295 (Form(
"hLogECellNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3296 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()),
3297 nptbins,ptmin,ptmax,150,-3,3);
3303 (Form(
"hLogECellNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3304 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()),
3305 nptbins,ptmin,ptmax,150,-3,3);
3311 (Form(
"hECellWeightNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3312 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()),
3313 nptbins,ptmin,ptmax,90,0,4.5);
3319 (Form(
"hECellWeightNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3320 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()),
3321 nptbins,ptmin,ptmax,90,0,4.5);
3327 (Form(
"hECellWeightNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3328 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()),
3329 nptbins,ptmin,ptmax,90,0,4.5);
3338 (Form(
"hMassEClusTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3339 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()),
3340 nptbins,ptmin,ptmax,nmassbins,massmin,massmax);
3355 (Form(
"hTimeDiffTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3356 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()),
3357 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
3363 (Form(
"hColRowTCardCorrelNCellLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3364 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3365 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3371 (Form(
"hColRowTCardCorrelNCellHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3372 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3373 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3381 (Form(
"hExoticTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3382 Form(
"exoticity vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3383 nptbins,ptmin,ptmax,200,-1,1);
3389 (Form(
"hTimeDiffExoTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3390 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()),
3391 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
3397 (Form(
"hColRowTCardCorrelNCellExoticLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3398 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3399 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3405 (Form(
"hColRowTCardCorrelNCellExoticHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3406 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3407 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3622 (Form(
"hLambda0Exoticity_EBin%d%s",i,add[tm].
Data()),
3623 Form(
"#sigma^{2}_{long} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3624 200,-1,1,ssbins,ssmin,ssmax);
3630 (Form(
"hLambda1Exoticity_EBin%d%s",i,add[tm].
Data()),
3631 Form(
"#sigma^{2}_{short} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3632 200,-1,1,ssbins,ssmin,ssmax);
3646 (Form(
"hNCellsExoticity_EBin%d%s",i,add[tm].
Data()),
3647 Form(
"#it{n}_{cells} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3648 200,-1,1,nceclbins,nceclmin,nceclmax);
3649 fhNCellsExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
3654 (Form(
"hTimeExoticity_EBin%d%s",i,add[tm].
Data()),
3656 200,-1,1,100,-25,25);
3657 fhTimeExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
3662 (Form(
"hNCellsTCardSameAndDiff_Exotic_EBin%d%s",i,add[tm].
Data()),
3663 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()),
3664 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
3704 (Form(
"hNCellsTCardSameAndDiff_EBin%d%s",i,add[tm].
Data()),
3705 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()),
3706 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
3712 (Form(
"hLambda0Lambda1_EBin%d%s",i,add[tm].
Data()),
3713 Form(
"#sigma^{2}_{long} vs #sigma^{2}_{short}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3714 ssbins,ssmin,ssmax,ssbins,ssmin,ssmax);
3730 for(
Int_t j = 0; j < 6; j++)
3732 for(
Int_t k = 0; k < 6; k++)
3735 (Form(
"hLambda0Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
3736 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()),
3737 200,-1,1,ssbins,ssmin,ssmax);
3743 (Form(
"hLambda1Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
3744 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()),
3745 200,-1,1,ssbins,ssmin,ssmax);
3842 for(
Int_t i = 0; i < 12; i++)
3845 (Form(
"hTCardCorrECellMaxDiff_Case%d%s",i,add[tm].
Data()),
3846 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()),
3847 nptbins,ptmin,ptmax,210,-1,20);
3853 (Form(
"hTCardCorrEClusterDiff_Case%d%s",i,add[tm].
Data()),
3854 Form(
"#it{E}_{cluster}-#it{E}_{cell} vs #it{E}_{cluster}, for (un)correlated cells in TCard, case %d %s",i,add[tm].
Data()),
3855 nptbins,ptmin,ptmax,210,-1,20);
3877 (Form(
"hTCardCorrTCellMaxDiff_Case%d%s",i,add[tm].
Data()),
3878 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()),
3879 nptbins,ptmin,ptmax,1000,-100,100);
3887 (Form(
"hTCardCorrECellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
3888 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()),
3889 nptbins,ptmin,ptmax,210,-1,20);
3895 (Form(
"hTCardCorrEClusterDiffExo_Case%d%s",i,add[tm].
Data()),
3896 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()),
3897 nptbins,ptmin,ptmax,210,-1,20);
3919 (Form(
"hTCardCorrTCellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
3920 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()),
3921 nptbins,ptmin,ptmax,1000,-100,100);
3932 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3938 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3944 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3950 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3956 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3962 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3968 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3974 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3982 (Form(
"hTMPhiResidual_EBin%d",i),
3983 Form(
"#Delta #varphi_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3984 200,-1,1,nresphibins,resphimin,resphimax);
3990 (Form(
"hTMEtaResidual_EBin%d",i),
3991 Form(
"#Delta #eta_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3992 200,-1,1,nresetabins,resetamin,resetamax);
3998 (Form(
"hTMPhiResidual_LooseCut_EBin%d",i),
3999 Form(
"#Delta #varphi_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
4000 200,-1,1,nresphibins,resphimin,resphimax);
4006 (Form(
"hTMEtaResidual_LooseCut_EBin%d",i),
4007 Form(
"#Delta #eta_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
4008 200,-1,1,nresetabins,resetamin,resetamax);
4036 (Form(
"hColRowM02"),
4038 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax,40,0.,2.);
4045 (Form(
"hColRowM02NCellCut"),
4046 Form(
"column vs row vs M02, %2.2f<#it{E}<%2.2f GeV #it{n}_{cells}^{w>0.01} > %d for Neutral",
4048 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax,40,0.,2.);
4058 Form(
"%2.2f<#it{E}_{1}<%2.2f GeV, %2.2f<#it{E}_{2}<%2.2f GeV, %2.2f<#sigma^{2}_{long}<%2.2f" 4059 "#it{M}_{#gamma #gamma} vs #it{n}_{1, cells}^{w>0.01} vs SM number trig cluster",
4068 (
"hInvMassNCellSMSame",
4069 Form(
"%2.2f<#it{E}_{1}<%2.2f GeV, %2.2f<#it{E}_{2}<%2.2f GeV, %2.2f<#sigma^{2}_{long}<%2.2f" 4070 "#it{M}_{#gamma #gamma} vs #it{n}_{1, cells}^{w>0.01} vs SM number both cluster",
4084 (
"hEMaxCellTimeM02SM",
4085 Form(
"time vs #sigma^{2}_{long} vs SM number, " 4094 (
"hEMaxCellTimeNCellSM",
4095 Form(
"time vs #it{n}_{cells}^{w>0.01} vs SM number, " 4104 (
"hESecCellTimeNCellSM",
4105 Form(
"secondary cell time vs #it{n}_{cells}^{w>0.01} vs SM number, " 4116 Int_t nbinsdeltacells = 19 ;
4122 for(
Int_t col = 0; col < 2; col++)
4125 (Form(
"hColRowFromCellMaxLowM02_Col%d",col),
4126 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, " 4128 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells);
4134 (Form(
"hColRowFromCellMaxHighM02_Col%d",col),
4135 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, " 4137 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells);
4148 (Form(
"hNCellsPerClusterM02M20_SM%d",i),
4149 Form(
" vs #sigma^{2}_{long} vs #it{n}_{cells}^{w>0.01} vs #sigma^{2}_{short}, " 4151 (
Int_t)nShShBins/1.5,0,maxShSh/1.5, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4160 (Form(
"hNCellsPerClusterM02NLM_SM%d",i),
4161 Form(
" vs #sigma^{2}_{long} vs #it{n}_{cells}^{w>0.01} vs N_{lm}, " 4163 10,0,10, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4171 (Form(
"hEMaxCellEClusterM02NCell_SM%d",i),
4172 Form(
"(#it{E}_{cluster} - #it{E}_{max cell})/#it{E}_{cluster} vs #sigma^{2}_{long} vs #it{n}_{cells}^{w>0.01}, " 4174 50,0,1., cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4191 (Form(
"hESecCellEMaxCellM02NCell_SM%d",i),
4192 Form(
"#it{E}_{cell}/#it{E}_{cell max} vs #sigma^{2}_{long} vs #it{n}_{cells}^{w>0.01}, " 4194 50,0,1., cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4201 (Form(
"hESecCellEClusterM02NCell_SM%d",i),
4202 Form(
"(#it{E}_{cluster} - #it{E}_{cell})/#it{E}_{cluster} vs #sigma^{2}_{long} vs #it{n}_{cells}^{w>0.01}, " 4204 50,0,1., cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4221 (Form(
"hEMaxESecCellNCellLowM02_SM%d",i),
4222 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, " 4224 24,0,12, 48,0,12,cellBins,cellMin,cellMax);
4231 (Form(
"hEMaxECrossNCellLowM02_SM%d",i),
4232 Form(
"#it{E}_{cell}^{max} vs exoticity vs #it{n}_{cells}^{w>0.01} vs %2.2f < #sigma^{2}_{long} < %2.2f, " 4234 24,0,12, 40,0.6,1.,cellBins,cellMin,cellMax);
4241 (Form(
"hEMaxESecCellNCellHighM02_SM%d",i),
4242 Form(
"#it{E}_{cell}^{max} vs #it{E}_{cell}^{secondary} vs #it{n}_{cells}^{w>0.01} vs " 4243 "%2.2f < #sigma^{2}_{long} < %2.2f, %2.2f<#it{E}<%2.2f GeV, SM=%d",
4245 24,0,12, 48,0,12,cellBins,cellMin,cellMax);
4252 (Form(
"hEMaxECrossNCellHighM02_SM%d",i),
4253 Form(
"#it{E}_{cell}^{max} vs exoticity vs #it{n}_{cells}^{w>0.01} " 4254 "for %2.2f < #sigma^{2}_{long} < %2.2f, %2.2f<#it{E}<%2.2f GeV, SM=%d",
4256 24,0,12, 40,0.6,1.,cellBins,cellMin,cellMax);
4264 Int_t nbinsdeltacells = 19 ;
4267 for(
Int_t col = 0; col < 2; col++)
4270 (Form(
"hNCellsPerClusterM02_SM%d_Col%d",i,col),
4271 Form(
"#it{E} vs #it{n}_{cells} vs #sigma^{2}_{long}, SM=%d, colum %d",i,col),
4272 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4279 (Form(
"hNCellsPerClusterM20_SM%d_Col%d",i,col),
4280 Form(
"#it{E} vs #it{n}_{cells} vs #sigma^{2}_{short}, SM=%d, colum %d",i,col),
4281 nEbins,minE,maxE,cellBins,cellMin,cellMax,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4290 (Form(
"hNCellsPerClusterNLM_SM%d_Col%d",i,col),
4291 Form(
"#it{E} vs #it{n}_{cells} vs n_{lm}, SM=%d, colum %d",i,col),
4292 nEbins,minE,maxE,cellBins,cellMin,cellMax,10,0,10);
4300 (Form(
"hColRowFromCellMaxLowM02_SM%d_Col%d",i,col),
4301 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, " 4302 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d",
4304 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,cellBins,cellMin,cellMax);
4311 (Form(
"hColRowFromCellMaxHighM02_SM%d_Col%d",i,col),
4312 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #it{n}_{cells}^{w>0.01} " 4313 "for %2.2f < #sigma^{2}_{long} < %2.2f, %2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d",
4315 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,cellBins,cellMin,cellMax);
4324 (Form(
"hColRowFromCellMaxEMaxSecDiffLowM02_SM%d_Col%d_NCellBin%d",i,col,j),
4325 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #Delta #it{E}_{max-secondary} " 4326 "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",
4328 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,120,0,12);
4335 (Form(
"hColRowFromCellMaxEMaxSecDiffHighM02_SM%d_Col%d_NCellBin%d",i,col,j),
4336 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #Delta #it{E}_{max-secondary} " 4337 "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",
4339 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,120,0,12);
4346 (Form(
"hColRowFromCellMaxEMaxSecDiffFracLowM02_SM%d_Col%d_NCellBin%d",i,col,j),
4347 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #Delta #it{E}_{max-secondary}/#it{E}_{max} " 4348 "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",
4350 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,40,0,1);
4357 (Form(
"hColRowFromCellMaxEMaxSecDiffFracHighM02_SM%d_Col%d_NCellBin%d",i,col,j),
4358 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #Delta #it{E}_{max-secondary}/#it{E}_{max}" 4359 " 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",
4361 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,40,0,1);
4368 (Form(
"hColRowFromCellMaxECellClusterRatLowM02_SM%d_Col%d_NCellBin%d",i,col,j),
4369 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #it{E}_{cell}/#it{E}_{cluster} " 4370 "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",
4372 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,50,0,1);
4379 (Form(
"hColRowFromCellMaxECellClusterRatHighM02_SM%d_Col%d_NCellBin%d",i,col,j),
4380 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #it{E}_{cell}/#it{E}_{cluster} " 4381 "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",
4383 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,50,0,1);
4400 (Form(
"hDeltaIEtaDeltaIPhi_%s",matchCase[imatch].
Data()),
4401 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),
4402 nEbins,minE,maxE,cellBins,cellMin,cellMax,cellBins,cellMin,cellMax);
4409 (Form(
"hDeltaIEtaDeltaIPhiTot_%s",matchCase[imatch].
Data()),
4410 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),
4411 nEbins,minE,maxE,cellBins,cellMin,cellMax,cellBins,cellMin,cellMax);
4418 (Form(
"hDeltaIA_%s",matchCase[imatch].
Data()),
4419 Form(
"Cluster *asymmetry* in cell units vs E for %s",matchCase[imatch].
Data()),
4420 nptbins,ptmin,ptmax,asyBins,asyMin,asyMax);
4421 fhDeltaIA[imatch]->SetXTitle(
"#it{E}_{cluster}");
4422 fhDeltaIA[imatch]->SetYTitle(
"#it{A}_{cell in cluster}");
4423 outputContainer->Add(
fhDeltaIA[imatch]);
4426 (Form(
"hDeltaIATot_%s",matchCase[imatch].
Data()),
4427 Form(
"Cluster *total asymmetry* in cell units vs E for %s",matchCase[imatch].
Data()),
4428 nptbins,ptmin,ptmax,asyBins,asyMin,asyMax);
4430 fhDeltaIATot[imatch]->SetYTitle(
"#it{A}_{cell in cluster}^{total}");
4434 (Form(
"hDeltaIAM02_%s",matchCase[imatch].
Data()),
4435 Form(
"Cluster *asymmetry* in cell units vs #sigma^{2}_{long} for %s",matchCase[imatch].
Data()),
4436 nEbins,minE,maxE,nShShBins,minShSh,maxShSh,asyBins,asyMin,asyMax);
4439 fhDeltaIAM02[imatch]->SetZTitle(
"#it{A}_{cell in cluster}");
4443 (Form(
"hDeltaIATotM02_%s",matchCase[imatch].
Data()),
4444 Form(
"Cluster *total asymmetry* in cell units vs #sigma^{2}_{long} for %s",matchCase[imatch].
Data()),
4445 nEbins,minE,maxE,nShShBins,minShSh,maxShSh,asyBins,asyMin,asyMax);
4448 fhDeltaIATotM02[imatch]->SetZTitle(
"#it{A}_{cell in cluster}^{total}");
4452 (Form(
"hDeltaIAM20_%s",matchCase[imatch].
Data()),
4453 Form(
"Cluster *asymmetry* in cell units vs #sigma^{2}_{short} for %s",matchCase[imatch].
Data()),
4454 nEbins,minE,maxE,(
Int_t)nShShBins/1.5,0,maxShSh/1.5,asyBins,asyMin,asyMax);
4457 fhDeltaIAM20[imatch]->SetZTitle(
"#it{A}_{cell in cluster}");
4461 (Form(
"hDeltaIATotM20_%s",matchCase[imatch].
Data()),
4462 Form(
"Cluster *total asymmetry* in cell units vs #sigma^{2}_{short} for %s",matchCase[imatch].
Data()),
4463 nEbins,minE,maxE,(
Int_t)nShShBins/1.5,0,maxShSh/1.5,asyBins,asyMin,asyMax);
4466 fhDeltaIATotM20[imatch]->SetZTitle(
"#it{A}_{cell in cluster}^{total}");
4470 (Form(
"hDeltaIANCells_%s",matchCase[imatch].
Data()),
4471 Form(
"Cluster *asymmetry* in cell units vs N cells in cluster for %s",matchCase[imatch].
Data()),
4472 nEbins,minE,maxE,cellBins,cellMin,cellMax,asyBins,asyMin,asyMax);
4479 (Form(
"hDeltaIATotNCells_%s",matchCase[imatch].
Data()),
4480 Form(
"Cluster *total asymmetry* in cell units vs N cells in cluster for %s",matchCase[imatch].
Data()),
4481 nEbins,minE,maxE,cellBins,cellMin,cellMax,asyBins,asyMin,asyMax);
4490 (Form(
"hDeltaIAOrigin_%s",matchCase[imatch].
Data()),
4491 Form(
"Cluster *asymmetry* in cell units vs E vs primary ID, for ID %s",matchCase[imatch].
Data()),
4492 nEbins,minE,maxE,10,-0.5,9.5,asyBins,asyMin,asyMax);
4496 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4497 fhDeltaIAOrigin[imatch]->GetYaxis()->SetBinLabel(ilabel,mcParticleStringLabel[ilabel-1]);
4501 (Form(
"hDeltaIATotOrigin_%s",matchCase[imatch].
Data()),
4502 Form(
"Cluster *total asymmetry* in cell units vs E vs primary ID, for ID %s",matchCase[imatch].
Data()),
4503 nEbins,minE,maxE,10,-0.5,9.5,asyBins,asyMin,asyMax);
4507 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4508 fhDeltaIATotOrigin[imatch]->GetYaxis()->SetBinLabel(ilabel,mcParticleStringLabel[ilabel-1]);
4514 (Form(
"hClusterMaxCellDiffM02_%s",matchCase[imatch].
Data()),
4515 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()),
4516 nEbins,minE,maxE, 20,0,1.,nShShBins,minShSh,maxShSh);
4523 (Form(
"hClusterTimeEnergyM02_%s",matchCase[imatch].
Data()),
4524 Form(
"#it{E} vs TOF vs #sigma^{2}_{long} for ID %s",matchCase[imatch].
Data()),
4525 nEbins,minE,maxE,45,-25.5,20.5,nShShBins,minShSh,maxShSh);
4532 (Form(
"hNCellsPerClusterM02_%s",matchCase[imatch].
Data()),
4533 Form(
"#it{E} vs #it{n}_{cells} vs #sigma^{2}_{long} for ID %s",matchCase[imatch].
Data()),
4534 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4541 (Form(
"hNCellsPerClusterM20_%s",matchCase[imatch].
Data()),
4542 Form(
"#it{E} vs #it{n}_{cells} vs #sigma^{2}_{short} for ID %s",matchCase[imatch].
Data()),
4543 nEbins,minE,maxE,cellBins,cellMin,cellMax,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4552 (Form(
"hNCellsPerClusterNLM_%s",matchCase[imatch].
Data()),
4553 Form(
"#it{E} vs #it{n}_{cells} vs n_{lm} for ID %s",matchCase[imatch].
Data()),
4554 nEbins,minE,maxE,cellBins,cellMin,cellMax,10,0,10);
4562 (Form(
"hSMM02_%s",matchCase[imatch].
Data()),
4563 Form(
"#it{E} vs SM number vs #sigma^{2}_{long}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4564 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,minShSh,maxShSh);
4565 fhSMM02[imatch]->SetXTitle(
"#it{E} (GeV)");
4566 fhSMM02[imatch]->SetYTitle(
"SM number");
4567 fhSMM02[imatch]->SetZTitle(
"#sigma^{2}_{long}");
4568 outputContainer->Add(
fhSMM02[imatch]);
4571 (Form(
"hSMM02NoCut_%s",matchCase[imatch].
Data()),
4572 Form(
"#it{E} vs SM number vs #sigma^{2}_{long} for ID %s",matchCase[imatch].
Data()),
4573 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,minShSh,maxShSh);
4580 (Form(
"hSMM20LowM02_%s",matchCase[imatch].
Data()),
4581 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",
4583 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4590 (Form(
"hSMM20LowM02NoCut_%s",matchCase[imatch].
Data()),
4591 Form(
"#it{E} vs SM number vs #sigma^{2}_{short}, for %2.2f<#sigma^{2}_{long}<%2.2f, for ID %s",
4593 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4600 (Form(
"hSMM20HighM02_%s",matchCase[imatch].
Data()),
4601 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",
4603 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4610 (Form(
"hSMM20HighM02NoCut_%s",matchCase[imatch].
Data()),
4611 Form(
"#it{E} vs SM number vs #sigma^{2}_{short}, for %2.2f<#sigma^{2}_{long}<%2.2f, for ID %s",
4613 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);