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),
50 fM02Min(0), fNCellMin(0),
51 fMinDistToBad(0), fNEBinCuts(0),
52 fEMinShape(0), fEMaxShape(100),
55 fdEdXMinEle(0), fdEdXMaxEle(0),
56 fdEdXMinHad(0), fdEdXMaxHad(0),
59 fInvMassMinECut(0), fInvMassMaxECut(0),
60 fInvMassMinM02Cut(0), fInvMassMaxM02Cut(0),
61 fInvMassMaxOpenAngle(0), fInvMassMaxTimeDifference(0),
63 fConstantTimeShift(0),
65 fClusterMomentum(), fClusterMomentum2(),
66 fCaloCellList(NULL), fCaloClusList(NULL),
71 fhEnergyTMEtaResidual1Cell(0), fhEnergyTMPhiResidual1Cell(0),
72 fhColRowExoticHighE1CellPosTime(0), fhColRowExoticHighE1CellNegTime(0), fhColRowExoticHighE1CellNulTime(0),
73 fhEnergyTMEtaResidualExotic(0), fhEnergyTMPhiResidualExotic(0),
74 fhColRowExoticHighEPosTime(0), fhColRowExoticHighENegTime(0), fhColRowExoticHighENulTime(0),
75 fhColRowHighEPosTime(0), fhColRowHighENegTime(0), fhColRowHighENulTime(0),
76 fhEnergyTMEtaResidualTCardCorrNoSelection1Cell(0), fhEnergyTMPhiResidualTCardCorrNoSelection1Cell(0),
77 fhEnergyTMEtaResidualTCardCorrNoSelectionExotic(0), fhEnergyTMPhiResidualTCardCorrNoSelectionExotic(0),
82 fhClusterMaxCellCloseCellRatioM02(0), fhClusterMaxCellECrossM02(0),
83 fhInvMassNCellSM(0), fhInvMassNCellSMSame(0),
84 fhColRowM02(0), fhColRowM02NCellCut(0),
85 fhEMaxCellTimeM02SM(0), fhEMaxCellTimeNCellSM(0), fhESecCellTimeNCellSM(0),
88 fhECellClusterRatio(0), fhECellClusterLogRatio(0),
89 fhEMaxCellClusterRatio(0), fhEMaxCellClusterLogRatio(0),
90 fhECellTotalRatio(0), fhECellTotalLogRatio(0),
91 fhECellTotalRatioMod(0), fhECellTotalLogRatioMod(0)
93 for(
Int_t i=0; i < 3; i++)
140 for(
Int_t i =0; i < 12; i++)
142 for(
Int_t j = 0; j < 4; j++)
144 for(
Int_t k = 0; k < 3; k++)
151 for(
Int_t j = 0; j < 5; j++)
159 for(
Int_t tm = 0; tm < 2; tm++)
255 for(
Int_t i = 0; i < 6; i++)
257 for(
Int_t j = 0; j < 6; j++)
330 for(
Int_t i = 0; i < 14; i++)
349 for(
Int_t j = 0; j < 6; j++)
351 for(
Int_t k = 0; k < 6; k++)
372 for(
Int_t i = 0; i < 12; i++)
398 for(
Int_t i = 0; i < 20; i++)
420 for(
Int_t j = 0; j < 3; j++)
464 Int_t icol = -1, irow = -1, iRCU = -1, icolAbs = -1, irowAbs = -1;
465 GetModuleNumberCellIndexesAbsCaloMap(absIdMax,GetCalorimeter(), icol, irow, iRCU, icolAbs, irowAbs);
469 for(
Int_t ie = 0; ie < fNEBinCuts; ie++)
471 if( energy >= fEBinCuts[ie] && energy < fEBinCuts[ie+1] ) ebin = ie;
474 Int_t ncells = clus->GetNCells();
475 Double_t time = clus->GetTOF()*1.e9;
476 if(time>400) time-=fConstantTimeShift;
478 Float_t deta = clus->GetTrackDz();
479 Float_t dphi = clus->GetTrackDx();
487 if(energy >= 5 && energy < 8)
488 fhColRowExoticLowE1Cell [matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
491 fhColRowExoticHighE1Cell[matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
493 if ( time > 5) fhColRowExoticHighE1CellPosTime->Fill(icolAbs,irowAbs,GetEventWeight());
494 else if( time < -5) fhColRowExoticHighE1CellNegTime->Fill(icolAbs,irowAbs,GetEventWeight());
495 else fhColRowExoticHighE1CellNulTime->Fill(icolAbs,irowAbs,GetEventWeight());
498 fhEnergyTime1Cell[matched]->Fill(energy,time,GetEventWeight());
499 fhEnergyTMEtaResidual1Cell->Fill(energy,deta,GetEventWeight());
500 fhEnergyTMPhiResidual1Cell->Fill(energy,dphi,GetEventWeight());
502 else if(exoticity > 0.97)
504 if(energy >= 5 && energy < 8)
505 fhColRowExoticLowE [matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
508 fhColRowExoticHighE[matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
510 if ( time > 5) fhColRowExoticHighEPosTime->Fill(icolAbs,irowAbs,GetEventWeight());
511 else if( time < -5) fhColRowExoticHighENegTime->Fill(icolAbs,irowAbs,GetEventWeight());
512 else fhColRowExoticHighENulTime->Fill(icolAbs,irowAbs,GetEventWeight());
515 fhEnergyTimeExotic[matched]->Fill(energy,time,GetEventWeight());
516 fhEnergyTMEtaResidualExotic->Fill(energy,deta,GetEventWeight());
517 fhEnergyTMPhiResidualExotic->Fill(energy,dphi,GetEventWeight());
519 for (
Int_t ipos = 0; ipos < ncells; ipos++)
521 Int_t absId = clus->GetCellsAbsId()[ipos];
523 Float_t eCell = fCaloCellList->GetCellAmplitude(absId);
524 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
527 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
529 if( absId == absIdMax || weight < 0.01 )
continue;
531 Int_t rowDiff = -100, colDiff = -100;
532 Bool_t sameTCard = GetCaloUtils()->IsAbsIDsFromTCard(absIdMax,absId,rowDiff,colDiff);
535 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
536 GetModuleNumberCellIndexesAbsCaloMap(absId,GetCalorimeter(), icol2, irow2, iRCU2, icolAbs2, irowAbs2);
540 if(energy >= 5 && energy < 8)
541 fhColRowExotic2ndCellDiffLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
543 fhColRowExotic2ndCellDiffHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
547 if(energy >= 5 && energy < 8)
548 fhColRowExotic2ndCellSameLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
550 fhColRowExotic2ndCellSameHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
555 else if ( energy > 8 )
557 if ( time > 5) fhColRowHighEPosTime->Fill(icolAbs,irowAbs,GetEventWeight());
558 else if( time < -5) fhColRowHighENegTime->Fill(icolAbs,irowAbs,GetEventWeight());
559 else fhColRowHighENulTime->Fill(icolAbs,irowAbs,GetEventWeight());
564 fhTMPhiResidualExoticityLooseCut[ebin]->Fill(exoticity,dphi);
565 fhTMEtaResidualExoticityLooseCut[ebin]->Fill(exoticity,deta);
568 else if ( energy > 8 && ncells > 1 && exoticity < 0.97)
570 if ( time > 5) fhColRowHighEPosTime->Fill(icolAbs,irowAbs,GetEventWeight());
571 else if( time < -5) fhColRowHighENegTime->Fill(icolAbs,irowAbs,GetEventWeight());
572 else fhColRowHighENulTime->Fill(icolAbs,irowAbs,GetEventWeight());
578 Int_t etaRegion = -1, phiRegion = -1;
579 GetCaloUtils()->GetEMCALSubregion(clus,fCaloCellList,etaRegion,phiRegion);
581 if ( etaRegion !=0 ) return ;
587 fhEnergyTimeTCardCorrNoSelection1Cell[matched]->Fill(energy,time,GetEventWeight());
588 fhEnergyTMEtaResidualTCardCorrNoSelection1Cell->Fill(energy,deta,GetEventWeight());
589 fhEnergyTMPhiResidualTCardCorrNoSelection1Cell->Fill(energy,dphi,GetEventWeight());
591 else if(exoticity > 0.97)
593 fhEnergyTimeTCardCorrNoSelectionExotic[matched]->Fill(energy,time,GetEventWeight());
594 fhEnergyTMEtaResidualTCardCorrNoSelectionExotic->Fill(energy,deta,GetEventWeight());
595 fhEnergyTMPhiResidualTCardCorrNoSelectionExotic->Fill(energy,dphi,GetEventWeight());
602 if ( m02 > 0.001 ) lamR = m20/m02;
604 Int_t absIdList[ncells];
606 Int_t nlm = GetCaloUtils()->GetNumberOfLocalMaxima(clus, fCaloCellList, absIdList, maxEList) ;
612 Int_t nCellWithWeight = 1;
620 Float_t eCellMax = fCaloCellList->GetCellAmplitude(absIdMax);
621 Double_t tCellMax = fCaloCellList->GetCellTime(absIdMax);
623 GetCaloUtils()->RecalibrateCellAmplitude(eCellMax, GetCalorimeter(), absIdMax);
624 GetCaloUtils()->RecalibrateCellTime(tCellMax, GetCalorimeter(), absIdMax, GetReader()->GetInputEvent()->GetBunchCrossNumber());
628 if(tCellMax>400) tCellMax-=fConstantTimeShift;
639 Bool_t sameTCard2ndMax = kFALSE;
643 Float_t eCellSameRowSameTCardNearCol = 0.;
644 Float_t eCellSameRowDiffTCardNearCol = 0.;
645 Double_t tCellSameRowSameTCardNearCol = 0.;
646 Double_t tCellSameRowDiffTCardNearCol = 0.;
653 for (
Int_t ipos = 0; ipos < ncells; ipos++)
655 Int_t absId = clus->GetCellsAbsId()[ipos];
657 Float_t eCell = fCaloCellList->GetCellAmplitude(absId);
658 Double_t tCell = fCaloCellList->GetCellTime(absId);
660 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
661 GetCaloUtils()->RecalibrateCellTime(tCell, GetCalorimeter(), absId, GetReader()->GetInputEvent()->GetBunchCrossNumber());
663 if(tCell>400) tCell-=fConstantTimeShift;
666 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
668 if( absId == absIdMax || weight < 0.01 )
continue;
672 fhECellClusRatNLM1TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
673 fhLogECellNLM1TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
674 fhECellWeightNLM1TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
678 fhECellClusRatNLM2TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
679 fhLogECellNLM2TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
680 fhECellWeightNLM2TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
684 fhECellClusRatNLM3TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
685 fhLogECellNLM3TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
686 fhECellWeightNLM3TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
689 Int_t rowDiff = -100, colDiff = -100;
690 Bool_t sameTCard = GetCaloUtils()->IsAbsIDsFromTCard(absIdMax,absId,rowDiff,colDiff);
696 fhECellSameClusRatNLM1TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
697 fhLogECellSameNLM1TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
698 fhECellSameWeightNLM1TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
702 fhECellSameClusRatNLM2TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
703 fhLogECellSameNLM2TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
704 fhECellSameWeightNLM2TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
708 fhECellSameClusRatNLM3TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
709 fhLogECellSameNLM3TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
710 fhECellSameWeightNLM3TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
720 Int_t indexType = -1;
725 if(TMath::Abs(rowDiff) == 1) nearRow = kTRUE;
726 if(TMath::Abs(colDiff) == 1) nearCol = kTRUE;
728 if ( rowDiff == 0 && colDiff != 0 )
730 if ( nearCol ) indexType = 6;
736 else if ( rowDiff != 0 && colDiff == 0 )
738 if ( nearRow ) indexType = 8;
746 if ( nearRow && nearCol) indexType = 10;
757 if ( rowDiff == 0 && colDiff != 0 )
759 if ( nearCol ) indexType = 0;
763 else if ( rowDiff != 0 && colDiff == 0 )
765 if ( nearRow ) indexType = 2;
770 if ( nearCol && nearRow ) indexType = 4;
775 if ( rowDiff == 0 && TMath::Abs(colDiff) == 1 )
779 eCellSameRowSameTCardNearCol = eCell;
780 tCellSameRowSameTCardNearCol = tCell;
784 eCellSameRowDiffTCardNearCol = eCell;
785 tCellSameRowDiffTCardNearCol = tCell;
791 Float_t eCellDiff = eCellMax - eCell;
792 Float_t eClusDiff = energy - eCell;
795 Float_t tCellDiff = tCellMax - tCell;
797 fhTCardCorrECellMaxDiff[indexType][matched]->Fill(energy, eCellDiff, GetEventWeight());
798 fhTCardCorrEClusterDiff[indexType][matched]->Fill(energy, eClusDiff, GetEventWeight());
801 fhTCardCorrTCellMaxDiff[indexType][matched]->Fill(energy, tCellDiff, GetEventWeight());
803 if ( fStudyExotic && exoticity > 0.97 )
805 fhTCardCorrECellMaxDiffExo[indexType][matched]->Fill(energy, eCellDiff, GetEventWeight());
806 fhTCardCorrEClusterDiffExo[indexType][matched]->Fill(energy, eClusDiff, GetEventWeight());
809 fhTCardCorrTCellMaxDiffExo[indexType][matched]->Fill(energy, tCellDiff, GetEventWeight());
813 if ( fStudyExotic && exoticity > 0.97 )
816 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
817 GetModuleNumberCellIndexesAbsCaloMap(absId,GetCalorimeter(), icol2, irow2, iRCU2, icolAbs2, irowAbs2);
821 if(energy >= 5 && energy < 8)
822 fhColRowTCardCorrNoSelectionExotic2ndCellDiffLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
824 fhColRowTCardCorrNoSelectionExotic2ndCellDiffHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
828 if(energy >= 5 && energy < 8)
829 fhColRowTCardCorrNoSelectionExotic2ndCellSameLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
831 fhColRowTCardCorrNoSelectionExotic2ndCellSameHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
835 if ( eCell > emax2nd )
841 sameTCard2ndMax = kTRUE;
847 sameTCard2ndMax = kFALSE;
882 Float_t ratioNcells = nCellWithWeight/(ncells*1.);
883 fhNCellsTCardCorrRatioWithWeightNoSelection[matched]->Fill(energy, ratioNcells, GetEventWeight());
887 if ( nCellWithWeight <= 1 )
return;
892 printf(
"AliAnaClusterShapeCorrelStudies: M02 %f, M20 %f, E %2.3f, ncell %d, n with weight %d; max cell E %2.3f\n",
893 m02,m20,energy,ncells,nCellWithWeight,eCellMax);
902 if(energy >= 5 && energy < 8)
903 fhColRowTCardCorrNoSelectionLowE [matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
905 fhColRowTCardCorrNoSelectionHighE[matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
907 if ( fStudyExotic && exoticity > 0.97 )
909 if(energy >= 5 && energy < 8)
910 fhColRowTCardCorrNoSelectionExoticLowE [matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
912 fhColRowTCardCorrNoSelectionExoticHighE[matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
915 fhTimeTCardCorrNoSelection [matched]->Fill(energy,time, GetEventWeight());
916 fhLambda0TCardCorrNoSelection[matched]->Fill(energy, m02, GetEventWeight());
917 fhLambda1TCardCorrNoSelection[matched]->Fill(energy, m20, GetEventWeight());
921 fhLambda0NLM1TCardCorrNoSelection[matched]->Fill(energy, m02, GetEventWeight());
922 fhLambda1NLM1TCardCorrNoSelection[matched]->Fill(energy, m20, GetEventWeight());
926 fhLambda0NLM2TCardCorrNoSelection[matched]->Fill(energy, m02, GetEventWeight());
927 fhLambda1NLM2TCardCorrNoSelection[matched]->Fill(energy, m20, GetEventWeight());
930 fhLambdaRTCardCorrNoSelection[matched]->Fill(energy,lamR, GetEventWeight());
931 fhNLocMaxTCardCorrNoSelection[matched]->Fill(energy, nlm, GetEventWeight());
932 fhExoticTCardCorrNoSelection [matched]->Fill(energy, exoticity, GetEventWeight());
936 fhEMaxRatNLM1TCardCorrNoSelection [matched]->Fill(energy, eCellMax/energy , GetEventWeight());
937 fhE2ndRatNLM1TCardCorrNoSelection [matched]->Fill(energy, emax2nd/energy , GetEventWeight());
938 fhE2ndEMaxRatNLM1TCardCorrNoSelection[matched]->Fill(energy, emax2nd/eCellMax, GetEventWeight());
942 fhEMaxRatNLM2TCardCorrNoSelection [matched]->Fill(energy, eCellMax/energy , GetEventWeight());
943 fhE2ndRatNLM2TCardCorrNoSelection [matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
944 fhE2ndEMaxRatNLM2TCardCorrNoSelection[matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
948 fhEMaxRatNLM3TCardCorrNoSelection [matched]->Fill(energy, eCellMax/energy , GetEventWeight());
949 fhE2ndRatNLM3TCardCorrNoSelection [matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
950 fhE2ndEMaxRatNLM3TCardCorrNoSelection[matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
957 fhE2ndSameRatNLM1TCardCorrNoSelection [matched]->Fill(energy, emax2nd/energy , GetEventWeight());
958 fhE2ndSameEMaxRatNLM1TCardCorrNoSelection[matched]->Fill(energy, emax2nd/eCellMax, GetEventWeight());
962 fhE2ndSameRatNLM2TCardCorrNoSelection [matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
963 fhE2ndSameEMaxRatNLM2TCardCorrNoSelection[matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
967 fhE2ndSameRatNLM3TCardCorrNoSelection [matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
968 fhE2ndSameEMaxRatNLM3TCardCorrNoSelection[matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
972 fhNCellsTCardCorrNoSelection [matched]->Fill(energy, ncells, GetEventWeight());
973 fhNCellsTCardCorrWithWeightNoSelection[matched]->Fill(energy, nCellWithWeight, GetEventWeight());
975 if(eCellSameRowSameTCardNearCol > 0 && eCellSameRowDiffTCardNearCol > 0)
977 Float_t eDiff = eCellSameRowSameTCardNearCol - eCellSameRowDiffTCardNearCol ;
978 Float_t tDiff = tCellSameRowSameTCardNearCol - tCellSameRowDiffTCardNearCol ;
980 fhSameRowDiffColAndTCardCellsEnergyDiffClusterE[matched]->Fill(energy , eDiff, GetEventWeight());
981 fhSameRowDiffColAndTCardCellsTimeDiffClusterE [matched]->Fill(energy , tDiff, GetEventWeight());
982 fhSameRowDiffColAndTCardCellsEnergyDiffCellMaxE[matched]->Fill(eCellMax, eDiff, GetEventWeight());
983 fhSameRowDiffColAndTCardCellsTimeDiffCellMaxE [matched]->Fill(eCellMax, tDiff, GetEventWeight());
985 if ( fStudyExotic && exoticity > 0.97 )
987 fhSameRowDiffColAndTCardCellsEnergyDiffClusterEExo[matched]->Fill(energy , eDiff, GetEventWeight());
988 fhSameRowDiffColAndTCardCellsTimeDiffClusterEExo [matched]->Fill(energy , tDiff, GetEventWeight());
989 fhSameRowDiffColAndTCardCellsEnergyDiffCellMaxEExo[matched]->Fill(eCellMax, eDiff, GetEventWeight());
990 fhSameRowDiffColAndTCardCellsTimeDiffCellMaxEExo [matched]->Fill(eCellMax, tDiff, GetEventWeight());
995 Int_t nCorrInd = nCorr;
996 if(nCorr > 4) nCorrInd = 5;
998 Int_t nCorrNoInd = nCorrNo;
999 if(nCorrNoInd > 4) nCorrNoInd = 5;
1012 fhLambda0NLM1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m02, GetEventWeight());
1013 fhLambda1NLM1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m20, GetEventWeight());
1017 fhLambda0NLM2TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m02, GetEventWeight());
1018 fhLambda1NLM2TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m20, GetEventWeight());
1021 fhLambda0TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m02, GetEventWeight());
1022 fhLambda1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m20, GetEventWeight());
1024 fhNLocMaxTCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, nlm, GetEventWeight());
1027 fhExoticTCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, exoticity, GetEventWeight());
1031 fhEMaxRatNLM1TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, eCellMax/energy , GetEventWeight());
1032 fhE2ndRatNLM1TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, emax2nd/energy , GetEventWeight());
1033 fhE2ndEMaxRatNLM1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, emax2nd/eCellMax, GetEventWeight());
1037 fhEMaxRatNLM2TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, eCellMax/energy, GetEventWeight());
1038 fhE2ndRatNLM2TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
1039 fhE2ndEMaxRatNLM2TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
1043 fhEMaxRatNLM3TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, eCellMax/energy, GetEventWeight());
1044 fhE2ndRatNLM3TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
1045 fhE2ndEMaxRatNLM3TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
1049 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1051 Int_t absId = clus->GetCellsAbsId()[ipos];
1053 Float_t eCell = fCaloCellList->GetCellAmplitude(absId);
1054 Double_t tCell = fCaloCellList->GetCellTime(absId);
1056 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
1057 GetCaloUtils()->RecalibrateCellTime(tCell, GetCalorimeter(), absId, GetReader()->GetInputEvent()->GetBunchCrossNumber());
1059 if(tCell>400) tCell-=fConstantTimeShift;
1062 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
1064 if( absId == absIdMax || weight < 0.01 )
continue;
1066 Float_t tDiffMaxSecondary = tCellMax - tCell;
1068 fhTimeDiffTCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, tDiffMaxSecondary, GetEventWeight());
1069 if ( fStudyExotic && exoticity > 0.97 )
1071 fhTimeDiffExoTCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, tDiffMaxSecondary, GetEventWeight());
1073 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
1074 GetModuleNumberCellIndexesAbsCaloMap(absId,GetCalorimeter(), icol2, irow2, iRCU2, icolAbs2, irowAbs2);
1076 Int_t rowDiff = -100, colDiff = -100;
1077 Bool_t sameTCard = GetCaloUtils()->IsAbsIDsFromTCard(absIdMax,absId,rowDiff,colDiff);
1081 if ( energy >= 5 && energy < 8 && nCorr == 0 )
1082 fhColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
1083 else if ( energy >= 8 && nCorr == 0 )
1084 fhColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
1088 if ( energy >= 5 && energy < 8 && nCorrNo == 0 )
1089 fhColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
1090 else if ( energy >= 8 && nCorrNo == 0 )
1091 fhColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
1097 fhECellClusRatNLM1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, eCell/energy, GetEventWeight());
1098 fhLogECellNLM1TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
1099 fhECellWeightNLM1TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, weight, GetEventWeight());
1103 fhECellClusRatNLM2TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, eCell/energy, GetEventWeight());
1104 fhLogECellNLM2TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
1105 fhECellWeightNLM2TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, weight, GetEventWeight());
1109 fhECellClusRatNLM3TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, eCell/energy, GetEventWeight());
1110 fhLogECellNLM3TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
1111 fhECellWeightNLM3TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, weight, GetEventWeight());
1116 if(m02 > fInvMassMinM02Cut && m02 < fInvMassMaxM02Cut)
1118 for(
Int_t jclus = 0 ; jclus < fCaloClusList->GetEntriesFast() ; jclus++)
1120 AliVCluster* clus2 = (AliVCluster*) fCaloClusList->At(jclus);
1123 Int_t absIdMax2 = GetCaloUtils()->GetMaxEnergyCell(fCaloCellList, clus2, maxCellFraction);
1125 Double_t tof2 = clus2->GetTOF()*1.e9;
1126 if(tof2>400) tof2-=fConstantTimeShift;
1133 if( absIdMax == absIdMax2
1134 || !IsGoodCluster(absIdMax2, clus2->GetM02(), clus2->GetNCells())
1135 || GetCaloUtils()->GetNumberOfLocalMaxima(clus2,fCaloCellList) > 1
1136 || clus2->GetM02() > fInvMassMaxM02Cut
1137 || clus2->GetM02() < fInvMassMinM02Cut
1138 || clus2->E() < fInvMassMinECut
1139 || clus2->E() > fInvMassMaxECut
1140 || TMath::Abs(diffTof) > fInvMassMaxTimeDifference
1145 clus2->GetMomentum(fClusterMomentum2,v);
1149 if(IsFiducialCutOn()) in2 = GetFiducialCut()->IsInFiducialCut(fClusterMomentum2.Eta(),fClusterMomentum2.Phi(),GetCalorimeter()) ;
1155 Double_t angle = fClusterMomentum.Angle(fClusterMomentum2.Vect());
1157 if( angle > fInvMassMaxOpenAngle )
continue;
1160 Float_t mass = (fClusterMomentum+fClusterMomentum2).M ();
1161 fhMassEClusTCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, mass, GetEventWeight());
1166 if ( energy >= 5 && energy < 8)
1167 fhColRowTCardCorrelNCellLowE [nCorrInd][nCorrNoInd][matched]->Fill(icolAbs, irowAbs, GetEventWeight());
1168 else if ( energy >= 8 )
1169 fhColRowTCardCorrelNCellHighE[nCorrInd][nCorrNoInd][matched]->Fill(icolAbs, irowAbs, GetEventWeight());
1184 if ( fStudyExotic && exoticity > 0.97 )
1186 if ( energy >= 5 && energy < 8)
1187 fhColRowTCardCorrelNCellExoticLowE [nCorrInd][nCorrNoInd][matched]->Fill(icolAbs, irowAbs, GetEventWeight());
1188 else if ( energy >= 8 )
1189 fhColRowTCardCorrelNCellExoticHighE[nCorrInd][nCorrNoInd][matched]->Fill(icolAbs, irowAbs, GetEventWeight());
1226 fhLambda0Lambda1 [ebin][matched]->Fill(m20, m02, GetEventWeight());
1227 fhNCellsTCardSameAndDiff[ebin][matched]->Fill(nCorrNo, nCorr, GetEventWeight());
1233 fhLambda0Exoticity[ebin][matched]->Fill(exoticity, m02, GetEventWeight());
1234 fhLambda1Exoticity[ebin][matched]->Fill(exoticity, m20, GetEventWeight());
1236 fhNCellsExoticity [ebin][matched]->Fill(exoticity, nCellWithWeight, GetEventWeight());
1237 fhTimeExoticity [ebin][matched]->Fill(exoticity, tCellMax, GetEventWeight());
1241 fhLambda0ExoticityPerNCell[nCorrInd][nCorrNoInd][matched]->Fill(exoticity, m02, GetEventWeight());
1242 fhLambda1ExoticityPerNCell[nCorrInd][nCorrNoInd][matched]->Fill(exoticity, m20, GetEventWeight());
1254 if ( exoticity > 0.97 )
1255 fhNCellsTCardSameAndDiffExotic[ebin][matched]->Fill(nCorrNo, nCorr, GetEventWeight());
1258 fhTMPhiResidualExoticity[ebin]->Fill(exoticity,dphi);
1259 fhTMEtaResidualExoticity[ebin]->Fill(exoticity,deta);
1269 Float_t nCellRat = nCorr*1. / ((nCorr+nCorrNo)*1.);
1270 fhNCellsTCardSameAndDiffFraction[matched]->Fill(energy, nCellRat, GetEventWeight());
1271 if ( fStudyExotic && exoticity > 0.97 )
1272 fhNCellsTCardSameAndDiffFractionExotic[matched]->Fill(energy, nCellRat, GetEventWeight());
1371 Int_t nCaloCellsPerCluster = clus->GetNCells();
1372 if ( nCaloCellsPerCluster < 2 )
return;
1377 Int_t nlm = GetCaloUtils()->GetNumberOfLocalMaxima(clus, fCaloCellList) ;
1389 for (
Int_t ipos = 0; ipos < nCaloCellsPerCluster; ipos++)
1391 Int_t absId = clus ->GetCellsAbsId()[ipos];
1392 Float_t eCell = fCaloCellList->GetCellAmplitude(absId) ;
1394 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
1396 if( absId == absIdMax || eCell < 0.01 )
continue;
1398 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
1400 if( weight < 0.01 )
continue;
1405 if ( nCell < 1 )
return;
1410 Int_t ietaMax=-1, iphiMax = 0, rcuMax = 0, icolAbs = -1, irowAbs = -1;
1411 Int_t smMax = GetModuleNumberCellIndexesAbsCaloMap(absIdMax, GetCalorimeter(),
1412 ietaMax, iphiMax, rcuMax, icolAbs, irowAbs);
1414 if ( matchedPID == 0 && energy > fEMinShape && energy < fEMaxShape )
1416 fhColRowM02->Fill(icolAbs,irowAbs,m02,GetEventWeight()) ;
1418 if ( nCell > fNCellMinShape )
1419 fhColRowM02NCellCut->Fill(icolAbs,irowAbs,m02,GetEventWeight()) ;
1427 Int_t etaRegion = -1, phiRegion = -1;
1428 GetCaloUtils()->GetEMCALSubregion(clus,fCaloCellList,etaRegion,phiRegion);
1430 if ( etaRegion !=0 ) return ;
1434 for (
Int_t ipos = 0; ipos < nCaloCellsPerCluster; ipos++)
1436 Int_t absId = clus ->GetCellsAbsId()[ipos];
1437 Float_t eCell = fCaloCellList->GetCellAmplitude(absId) ;
1439 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
1441 if( absId == absIdMax || eCell < 0.01 )
continue;
1443 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
1445 if( weight < 0.01 )
continue;
1447 Float_t fracCell = eCell/eCellMax;
1448 fhClusterMaxCellCloseCellRatioM02->Fill(energy, fracCell, m02, GetEventWeight());
1450 Float_t fracClus = (energy-eCell)/energy;
1452 Float_t logECell = TMath::Log(eCell);
1463 if ( energy > fEMinShape && energy < fEMaxShape && matchedPID == 0 )
1465 fhESecCellEMaxCellM02NCellPerSM[smMax]->Fill(fracCell, nCell, m02, GetEventWeight());
1466 fhESecCellEClusterM02NCellPerSM[smMax]->Fill(fracClus, nCell, m02, GetEventWeight());
1467 fhESecCellLogM02NCellPerSM [smMax]->Fill(logECell, nCell, m02, GetEventWeight());
1469 if ( m02 > 0.1 && m02 < 0.3 )
1471 fhEMaxESecCellNCellLowM02PerSM [smMax]->Fill(eCellMax, eCell , nCell, GetEventWeight());
1472 fhEMaxECrossNCellLowM02PerSM [smMax]->Fill(eCellMax, eCrossFrac, nCell, GetEventWeight());
1474 Double_t time = fCaloCellList->GetCellTime(absId);
1475 GetCaloUtils()->RecalibrateCellTime(time, GetCalorimeter(), absId,GetReader()->GetInputEvent()->GetBunchCrossNumber());
1477 if(time>400) time-=fConstantTimeShift;
1479 fhESecCellTimeNCellSM->Fill(time, smMax, nCell, GetEventWeight());
1481 else if ( m02 > 0.5 && m02 < 2 )
1483 fhEMaxESecCellNCellHighM02PerSM[smMax]->Fill(eCellMax, eCell , nCell, GetEventWeight());
1484 fhEMaxECrossNCellHighM02PerSM [smMax]->Fill(eCellMax, eCrossFrac, nCell, GetEventWeight());
1494 Int_t sm = GetModuleNumberCellIndexes(absId,GetCalorimeter(), ieta, iphi, rcu);
1496 if(dIphi < TMath::Abs(iphi-iphiMax)) dIphi = TMath::Abs(iphi-iphiMax);
1497 if(iphi-iphiMax < 0 && dIphiNeg > iphi-iphiMax) dIphiNeg = iphi-iphiMax;
1498 if(iphi-iphiMax > 0 && dIphiPos < iphi-iphiMax) dIphiPos = iphi-iphiMax;
1502 if(dIeta < TMath::Abs(ieta-ietaMax)) dIeta = TMath::Abs(ieta-ietaMax);
1503 if(ieta-ietaMax < 0 && dIetaNeg > ieta-ietaMax) dIetaNeg = ieta-ietaMax;
1504 if(ieta-ietaMax > 0 && dIetaPos < ieta-ietaMax) dIetaPos = ieta-ietaMax;
1506 if ( energy > fEMinShape && energy < fEMaxShape && matchedPID == 0 )
1509 if ( nCell == 2 || nCell == 3 ) nCellBin = 0;
1510 else if ( nCell == 4 || nCell == 5 ) nCellBin = 1;
1511 if ( m02 > 0.1 && m02 < 0.3 )
1513 fhColRowFromCellMaxLowM02PerSM[smMax][ietaMax%2]->Fill(ietaMax-ieta, iphiMax-iphi, nCell, GetEventWeight());
1514 fhColRowFromCellMaxLowM02 [ietaMax%2]->Fill(ietaMax-ieta, iphiMax-iphi, GetEventWeight());
1516 fhColRowFromCellMaxEMaxSecDiffFracLowM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi,(eCellMax-eCell)/eCellMax, GetEventWeight());
1517 fhColRowFromCellMaxEMaxSecDiffLowM02PerSM [smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi, eCellMax-eCell, GetEventWeight());
1518 fhColRowFromCellMaxECellClusterRatLowM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi, eCell/energy , GetEventWeight());
1520 else if (m02 > 0.5 && m02 < 2)
1522 fhColRowFromCellMaxHighM02PerSM[smMax][ietaMax%2]->Fill(ietaMax-ieta, iphiMax-iphi, nCell, GetEventWeight());
1523 fhColRowFromCellMaxHighM02 [ietaMax%2]->Fill(ietaMax-ieta, iphiMax-iphi, GetEventWeight());
1525 fhColRowFromCellMaxEMaxSecDiffFracHighM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi,(eCellMax-eCell)/eCellMax, GetEventWeight());
1526 fhColRowFromCellMaxEMaxSecDiffHighM02PerSM [smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi, eCellMax-eCell, GetEventWeight());
1527 fhColRowFromCellMaxECellClusterRatHighM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi, eCell/energy , GetEventWeight());
1534 Int_t ietaShift = ieta;
1535 Int_t ietaMaxShift = ietaMax;
1537 if (ieta > ietaMax) ietaMaxShift+=48;
1538 else ietaShift +=48;
1540 if(dIeta < TMath::Abs(ietaShift-ietaMaxShift)) dIeta = TMath::Abs(ietaShift-ietaMaxShift);
1541 if(ietaShift-ietaMaxShift < 0 && dIetaNeg > ietaShift-ietaMaxShift) dIetaNeg = ietaShift-ietaMaxShift;
1542 if(ietaShift-ietaMaxShift > 0 && dIetaPos < ietaShift-ietaMaxShift) dIetaPos = ietaShift-ietaMaxShift;
1548 if ( energy > fEMinShape && energy < fEMaxShape && matchedPID == 0 )
1550 fhNCellsPerClusterM02M20PerSM[smMax]->Fill(m20, nCell, m02, GetEventWeight());
1551 fhNCellsPerClusterM02NLMPerSM[smMax]->Fill(nlm, nCell, m02, GetEventWeight());
1555 if ( nCell == 2 || nCell == 3 ) nCellBin = 0;
1556 else if ( nCell == 4 || nCell == 5 ) nCellBin = 1;
1557 if ( m02 > 0.1 && m02 < 0.3 )
1558 fhColRowFromCellMaxECellClusterRatLowM02PerSM [smMax][ietaMax%2][nCellBin]->Fill(0., 0., eCellMax/energy, GetEventWeight());
1559 else if (m02 > 0.5 && m02 < 2)
1560 fhColRowFromCellMaxECellClusterRatHighM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(0., 0., eCellMax/energy, GetEventWeight());
1564 fhClusterMaxCellECrossM02->Fill(energy, eCrossFrac, m02, GetEventWeight());
1569 fhClusterMaxCellDiffM02[matchedPID]->Fill(energy, maxFrac, m02, GetEventWeight());
1570 fhClusterTimeEnergyM02 [matchedPID]->Fill(energy, tmax , m02, GetEventWeight());
1571 fhNCellsPerClusterM02 [matchedPID]->Fill(energy, nCell , m02, GetEventWeight());
1572 fhNCellsPerClusterM20 [matchedPID]->Fill(energy, nCell , m20, GetEventWeight());
1573 fhNCellsPerClusterNLM [matchedPID]->Fill(energy, nCell , nlm, GetEventWeight());
1575 fhSMNCell [matchedPID]->Fill(energy, smMax , nCell, GetEventWeight());
1576 fhSMM02NoCut [matchedPID]->Fill(energy, smMax , m02 , GetEventWeight());
1577 if ( m02 > 0.1 && m02 < 0.3 )
1578 fhSMM20LowM02NoCut [matchedPID]->Fill(energy, smMax , m20 , GetEventWeight());
1579 else if ( m02 > 0.5 && m02 < 2 )
1580 fhSMM20HighM02NoCut[matchedPID]->Fill(energy, smMax , m20 , GetEventWeight());
1583 if ( energy > fEMinShape && energy < fEMaxShape )
1585 fhSMNCellM02[matchedPID]->Fill(smMax , nCell, m02, GetEventWeight());
1587 if ( matchedPID == 0 )
1589 if ( m02 > 0.1 && m02 < 0.3 )
1590 fhEMaxCellTimeNCellSM->Fill(tmax, smMax, nCell, GetEventWeight());
1592 fhEMaxCellTimeM02SM->Fill(tmax, smMax, m02, GetEventWeight());
1594 fhEMaxCellEClusterM02NCellPerSM[smMax]->Fill(maxFrac, nCell, m02, GetEventWeight());
1595 fhEMaxCellLogM02NCellPerSM [smMax]->Fill(TMath::Log(eCellMax), nCell, m02, GetEventWeight());
1599 if ( nCell > fNCellMinShape )
1601 if ( m02 > 0.1 && m02 < 0.3 )
1602 fhSMM20LowM02 [matchedPID]->Fill(energy, smMax , m20 , GetEventWeight());
1603 else if ( m02 > 0.5 && m02 < 2 )
1604 fhSMM20HighM02[matchedPID]->Fill(energy, smMax , m20 , GetEventWeight());
1606 fhSMM02 [matchedPID]->Fill(energy, smMax , m02, GetEventWeight());
1607 fhColM02[matchedPID]->Fill(energy, ietaMax, m02, GetEventWeight());
1608 fhRowM02[matchedPID]->Fill(energy, iphiMax, m02, GetEventWeight());
1612 if ( fStudyShapeParam )
1615 Float_t dIA = 1.*(dIphi-dIeta)/(dIeta+dIphi);
1616 Float_t dIATot = 1.*((dIphiPos-dIphiNeg)-(dIetaPos-dIetaNeg))/((dIetaPos-dIetaNeg)+(dIphiPos-dIphiNeg));
1625 AliDebug(1,Form(
"E %2.2f, nCell %d, dPhi %d, dEta %d, dIA %2.2f, match %d",energy,nCell, dIphi,dIeta,dIA,matchedPID));
1635 fhDeltaIANCells [matchedPID]->Fill(energy, nCell , dIA , GetEventWeight());
1636 fhDeltaIATotNCells [matchedPID]->Fill(energy, nCell , dIATot, GetEventWeight());
1638 if ( nCell > fNCellMinShape )
1640 fhDeltaIEtaDeltaIPhi[matchedPID]->Fill(energy, dIeta, dIphi, GetEventWeight());
1641 fhDeltaIA [matchedPID]->Fill(energy, dIA , GetEventWeight());
1642 fhDeltaIAM02 [matchedPID]->Fill(energy, m02 , dIA , GetEventWeight());
1643 fhDeltaIAM20 [matchedPID]->Fill(energy, m20 , dIA , GetEventWeight());
1645 fhDeltaIEtaDeltaIPhiTot[matchedPID]->Fill(energy,dIetaPos-dIetaNeg, dIphiPos-dIphiNeg, GetEventWeight());
1646 fhDeltaIATot [matchedPID]->Fill(energy, dIATot , GetEventWeight());
1647 fhDeltaIATotM02 [matchedPID]->Fill(energy, m02 , dIATot, GetEventWeight());
1648 fhDeltaIATotM20 [matchedPID]->Fill(energy, m20 , dIATot, GetEventWeight());
1651 if ( IsDataMC() && mcIndex > -1 && mcIndex < 10 && nCell > fNCellMinShape )
1653 fhDeltaIAOrigin [matchedPID]->Fill(energy, mcIndex, dIA , GetEventWeight());
1654 fhDeltaIATotOrigin[matchedPID]->Fill(energy, mcIndex, dIATot, GetEventWeight());
1658 Float_t dispp= 0., dEta = 0., dPhi = 0.;
1659 Float_t sEta = 0., sPhi = 0., sEtaPhi = 0.;
1660 if ( GetCalorimeter() ==
kEMCAL )
1662 GetCaloUtils()->GetEMCALRecoUtils()->RecalculateClusterShowerShapeParameters(GetEMCALGeometry(), fCaloCellList, clus,
1663 l0, l1, dispp, dEta, dPhi, sEta, sPhi, sEtaPhi);
1666 if(sEta+sPhi>0.0001) sEtaPhiA = (sPhi-sEta)/(sEta+sPhi);
1668 AliDebug(2,Form(
"Recalculate shower shape org: m02 %2.2f, m20 %2.2f, disp %2.2f;"
1669 " new: m02 %2.2f, m20 %2.2f, disp %2.2f; "
1670 "mEta %2.2f, mPhi %2.2f, mEtaPhi %2.2f, A_EtaPhi %2.2f; dEta %2.2f dPhi %2.2f",
1671 m02,m20,clus->GetDispersion(),l0,l1,dispp,sEta,sPhi,sEtaPhi,sEtaPhiA,dEta,dPhi));
1673 fhNCellsPerClusterMEta [matchedPID]->Fill(energy, nCell, sEta , GetEventWeight());
1674 fhNCellsPerClusterMPhi [matchedPID]->Fill(energy, nCell, sPhi , GetEventWeight());
1675 fhNCellsPerClusterMEtaPhi [matchedPID]->Fill(energy, nCell, sEtaPhi , GetEventWeight());
1676 fhNCellsPerClusterMEtaPhiA[matchedPID]->Fill(energy, nCell, sEtaPhiA, GetEventWeight());
1681 if ( matchedPID == 0 && energy > fEMinShape && energy < fEMaxShape &&
1682 m02 > fInvMassMinM02Cut && m02 < fInvMassMaxM02Cut )
1684 for(
Int_t jclus = 0 ; jclus < fCaloClusList->GetEntriesFast() ; jclus++)
1686 AliVCluster* clus2 = (AliVCluster*) fCaloClusList->At(jclus);
1689 Int_t absIdMax2 = GetCaloUtils()->GetMaxEnergyCell(fCaloCellList, clus2, maxCellFraction);
1691 Double_t tof2 = clus2->GetTOF()*1.e9;
1692 if(tof2>400) tof2-=fConstantTimeShift;
1699 if( absIdMax == absIdMax2
1700 || !IsGoodCluster(absIdMax2, clus2->GetM02(), clus2->GetNCells())
1701 || GetCaloUtils()->GetNumberOfLocalMaxima(clus2,fCaloCellList) > 1
1702 || clus2->GetM02() > fInvMassMaxM02Cut
1703 || clus2->GetM02() < fInvMassMinM02Cut
1704 || clus2->E() < fInvMassMinECut
1705 || clus2->E() > fInvMassMaxECut
1706 || TMath::Abs(diffTof) > fInvMassMaxTimeDifference
1711 clus2->GetMomentum(fClusterMomentum2,v);
1715 if(IsFiducialCutOn()) in2 = GetFiducialCut()->IsInFiducialCut(fClusterMomentum2.Eta(),fClusterMomentum2.Phi(),GetCalorimeter()) ;
1721 Double_t angle = fClusterMomentum.Angle(fClusterMomentum2.Vect());
1723 if( angle > fInvMassMaxOpenAngle )
continue;
1726 Float_t mass = (fClusterMomentum+fClusterMomentum2).M ();
1727 fhInvMassNCellSM->Fill(mass, nCell, smMax, GetEventWeight());
1729 Int_t smMax2 = GetModuleNumber(clus2);
1730 if(smMax == smMax2) fhInvMassNCellSMSame->Fill(mass, nCell, smMax, GetEventWeight());
1735 if ( IsDataMC() && mcIndex > -1 && mcIndex < 10)
1737 fhOriginE [matchedPID]->Fill(energy, mcIndex, GetEventWeight());
1738 fhOriginM02[matchedPID]->Fill(energy, mcIndex, m02, GetEventWeight());
1754 AliVTrack *track = GetCaloUtils()->GetMatchedTrack(clus, GetReader()->GetInputEvent());
1764 Double_t dedx = track->GetTPCsignal();
1779 if(fdEdXMinEle == 0 || fdEdXMaxEle == 0 || fdEdXMinHad == 0 || fdEdXMaxHad == 0)
1780 InitdEdXParameters();
1782 if ( dedx >= fdEdXMinEle && dedx < fdEdXMaxEle ) matchedPID = 1;
1783 else if ( dedx >= fdEdXMinHad && dedx < fdEdXMaxHad ) matchedPID = 2;
1786 AliDebug(1,Form(
"dEdX out of range %2.2f",dedx));
1804 Int_t nCaloCellsPerCluster = 0 ;
1815 for(
Int_t iclus = 0; iclus < nCaloClusters; iclus++)
1817 AliDebug(1,Form(
"Cluster: %d/%d, data %d",iclus+1,nCaloClusters,
GetReader()->GetDataType()));
1819 AliVCluster* clus = (AliVCluster*)
fCaloClusList->At(iclus);
1822 if ( clus->GetDistanceToBadChannel() <
fMinDistToBad )
return ;
1826 if ( nModule < fFirstModule || nModule >
fLastModule )
1828 AliDebug(1,Form(
"Cluster module out of range %d",nModule));
1838 Double_t tof = clus->GetTOF()*1.e9;
1856 AliDebug(1,Form(
"Remove cluster with phi %2.2f and eta %2.2f", phi*TMath::RadToDeg(), eta));
1860 AliDebug(1,Form(
"cluster: E %2.3f, pT %2.3f, eta %2.3f, phi %2.3f",e,pt,eta,phi*TMath::RadToDeg()));
1864 nCaloCellsPerCluster = clus->GetNCells();
1868 AliDebug(1,Form(
"Accept cluster? %d",goodCluster));
1870 if(!goodCluster)
continue;
1899 Int_t matchedPID = 0;
1916 if ( ampMax > 0.01 )
1924 if (
fStudyShape && matchedPID >= 0 && matchedPID < 3 )
1945 const Int_t buffersize = 255;
1946 char onePar[buffersize] ;
1948 snprintf(onePar,buffersize,
"--- AliAnaClusterShapeCorrelStudies ---:") ;
1969 return new TObjString(parList) ;
1979 outputContainer->SetName(
"ClusterShapeStudies") ;
1990 TString mcParticleStringLabel[] = {
"Merged #gamma#gamma",
"Decay #gamma",
"Direct #gamma",
"e^{#pm}",
"#pi^{#pm}",
"k^{#pm}",
"p",
"#bar{p}",
"n",
"#bar{n}"};
1994 TString matchCase[] = {
"Neutral",
"Electron",
"Hadron"};
2022 Float_t rowcellmin = fNMaxRowsFullMin-1.5;
2031 Int_t nShShBins = 200;
2041 Int_t cellBins = 15;
2051 (
"hColRowHighEPosTime",
2052 "column vs row, exo < 0.97, E > 8 GeV, t > 5 ns",
2053 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2059 (
"hColRowHighENegTime",
2060 "column vs row, exo < 0.97, E > 8 GeV, t < -5 ns",
2061 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2067 (
"hColRowHighENulTime",
2068 "column vs row, exo < 0.97, E > 8 GeV, -5 < t < 5 ns",
2069 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2077 (
"hColRowExoticHighE1CellPosTime",
2078 "column vs row, 1 cell, E > 8 GeV, t > 5 ns",
2079 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2085 (
"hColRowExoticHighEPosTime",
2086 "column vs row, exo > 0.97, E > 8 GeV, t > 5 ns",
2087 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2093 (
"hColRowExoticHighE1CellNegTime",
2094 "column vs row, 1 cell, E > 8 GeV, t < -5 ns",
2095 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2101 (
"hColRowExoticHighENegTime",
2102 "column vs row, exo > 0.97, E > 8 GeV, t < -5 ns",
2103 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2109 (
"hColRowExoticHighE1CellNulTime",
2110 "column vs row, 1 cell, E > 8 GeV, -5 < t < 5 ns",
2111 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2117 (
"hColRowExoticHighENulTime",
2118 "column vs row, exo > 0.97, E > 8 GeV, -5 < t < 5 ns",
2119 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2125 TString add[] = {
"",
"TrackMatched"};
2126 for(
Int_t tm = 0; tm < 2; tm++)
2129 (Form(
"hColRowTCardCorrNoSelectionLowE%s",add[tm].
Data()),
2130 Form(
"column vs row, max E cell for TCard correlation selected clusters, 5 < E < 8 GeV %s",add[tm].
Data()),
2131 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2137 (Form(
"hColRowTCardCorrNoSelectionHighE%s",add[tm].
Data()),
2138 Form(
"column vs row, max E cell for TCard correlation selected clusters, E > 8 GeV %s",add[tm].
Data()),
2139 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2147 (Form(
"hNCellsTCardCorrNoSelection%s",add[tm].
Data()),
2148 Form(
"# custer # cells vs #it{E} %s",add[tm].
Data()),
2149 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
2155 (Form(
"hNCellsTCardCorrWithWeightNoSelection%s",add[tm].
Data()),
2156 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
2157 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
2163 (Form(
"hNCellsTCardCorrRatioWithWeightNoSelection%s",add[tm].
Data()),
2164 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
2165 nptbins,ptmin,ptmax, 100,0,1);
2171 (Form(
"hTimeTCardCorrNoSelection%s",add[tm].
Data()),
2172 Form(
"#it{time} vs #it{E} %s",add[tm].
Data()),
2173 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2179 (Form(
"hLambda0TCardCorrNoSelection%s",add[tm].
Data()),
2180 Form(
"#lambda^{2}_{0} vs #it{E} %s",add[tm].
Data()),
2181 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2187 (Form(
"hLambda1TCardCorrNoSelection%s",add[tm].
Data()),
2188 Form(
"#lambda^{2}_{1} vs #it{E} %s",add[tm].
Data()),
2189 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2195 (Form(
"hLambda0NLM1TCardCorrNoSelection%s",add[tm].
Data()),
2196 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=1 %s",add[tm].
Data()),
2197 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2203 (Form(
"hLambda1NLM1TCardCorrNoSelection%s",add[tm].
Data()),
2204 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=1 %s",add[tm].
Data()),
2205 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2211 (Form(
"hLambda0NLM2TCardCorrNoSelection%s",add[tm].
Data()),
2212 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=2 %s",add[tm].
Data()),
2213 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2219 (Form(
"hLambda1NLM2TCardCorrNoSelection%s",add[tm].
Data()),
2220 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=2 %s",add[tm].
Data()),
2221 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2228 (Form(
"hLambdaRTCardCorrNoSelection%s",add[tm].
Data()),
2229 Form(
"#lambda^{1}_{0}/#lambda^{2}_{0} vs #it{E} %s",add[tm].
Data()),
2230 nptbins,ptmin,ptmax,110,0,1.1);
2236 (Form(
"hNLocMaxTCardCorrNoSelection%s",add[tm].
Data()),
2237 Form(
"#it{n}_{LM} vs E %s",add[tm].
Data()),
2238 nptbins,ptmin,ptmax,10,0,10);
2244 (Form(
"hEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2245 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
2246 nptbins,ptmin,ptmax,100,0,1);
2252 (Form(
"hEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2253 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2254 nptbins,ptmin,ptmax,100,0,1);
2260 (Form(
"hEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2261 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2262 nptbins,ptmin,ptmax,100,0,1);
2268 (Form(
"hE2ndRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2269 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
2270 nptbins,ptmin,ptmax,100,0,1);
2276 (Form(
"hE2ndRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2277 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2278 nptbins,ptmin,ptmax,100,0,1);
2284 (Form(
"hE2ndRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2285 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2286 nptbins,ptmin,ptmax,100,0,1);
2293 (Form(
"hE2ndEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2294 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
2295 nptbins,ptmin,ptmax,100,0,1);
2301 (Form(
"hE2ndEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2302 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2303 nptbins,ptmin,ptmax,100,0,1);
2309 (Form(
"hE2ndEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2310 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2311 nptbins,ptmin,ptmax,100,0,1);
2319 (Form(
"hE2ndSameRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2320 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()),
2321 nptbins,ptmin,ptmax,100,0,1);
2327 (Form(
"hE2ndSameRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2328 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()),
2329 nptbins,ptmin,ptmax,100,0,1);
2335 (Form(
"hE2ndSameRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2336 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()),
2337 nptbins,ptmin,ptmax,100,0,1);
2344 (Form(
"hE2ndSameEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2345 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()),
2346 nptbins,ptmin,ptmax,100,0,1);
2352 (Form(
"hE2ndSameEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2353 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()),
2354 nptbins,ptmin,ptmax,100,0,1);
2360 (Form(
"hE2ndSameEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2361 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()),
2362 nptbins,ptmin,ptmax,100,0,1);
2370 (Form(
"hECellClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2371 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
2372 nptbins,ptmin,ptmax,100,0,1);
2378 (Form(
"hECellClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2379 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2380 nptbins,ptmin,ptmax,100,0,1);
2386 (Form(
"hECellClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2387 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2388 nptbins,ptmin,ptmax,100,0,1);
2394 (Form(
"hLogECellNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2395 Form(
"log(#it{E}_{cell}) vs E_{cluster}, #it{n}_{LM}=1, w > 0.01 %s",add[tm].
Data()),
2396 nptbins,ptmin,ptmax,150,-3,3);
2402 (Form(
"hLogECellNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2403 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}=2, w > 0.01 %s",add[tm].
Data()),
2404 nptbins,ptmin,ptmax,150,-3,3);
2410 (Form(
"hLogECellNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2411 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}>2, w > 0.01 %s",add[tm].
Data()),
2412 nptbins,ptmin,ptmax,150,-3,3);
2419 (Form(
"hECellWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2420 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
2421 nptbins,ptmin,ptmax,90,0,4.5);
2427 (Form(
"hECellWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2428 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2429 nptbins,ptmin,ptmax,90,0,4.5);
2435 (Form(
"hECellWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2436 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster}))vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2437 nptbins,ptmin,ptmax,90,0,4.5);
2444 (Form(
"hECellSameClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2445 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()),
2446 nptbins,ptmin,ptmax,100,0,1);
2452 (Form(
"hECellSameClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2453 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, cell from same T-Card as leading, #it{n}_{LM}=2 %s",add[tm].
Data()),
2454 nptbins,ptmin,ptmax,100,0,1);
2460 (Form(
"hECellSameClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2461 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, cell from same T-Card as leading, #it{n}_{LM}>2 %s",add[tm].
Data()),
2462 nptbins,ptmin,ptmax,100,0,1);
2468 (Form(
"hLogECellSameNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2469 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()),
2470 nptbins,ptmin,ptmax,150,-3,3);
2476 (Form(
"hLogECellSameNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2477 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()),
2478 nptbins,ptmin,ptmax,150,-3,3);
2484 (Form(
"hLogECellSameNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2485 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()),
2486 nptbins,ptmin,ptmax,150,-3,3);
2493 (Form(
"hECellSameWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2494 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()),
2495 nptbins,ptmin,ptmax,90,0,4.5);
2501 (Form(
"hECellSameWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2502 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()),
2503 nptbins,ptmin,ptmax,90,0,4.5);
2509 (Form(
"hECellSameWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2510 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()),
2511 nptbins,ptmin,ptmax,90,0,4.5);
2517 (Form(
"hExoticTCardCorrNoSelection%s",add[tm].
Data()),
2518 Form(
"exoticity vs #it{E} %s",add[tm].
Data()),
2519 nptbins,ptmin,ptmax,200,-1,1);
2525 (Form(
"hNCellsTCardSameAndDiffFraction%s",add[tm].
Data()),
2526 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, %s",add[tm].
Data()),
2527 nptbins,ptmin,ptmax,100,0,1);
2533 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterE%s",add[tm].
Data()),
2534 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()),
2535 nptbins,ptmin,ptmax,200,-10,10);
2541 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterE%s",add[tm].
Data()),
2542 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()),
2543 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2549 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxE%s",add[tm].
Data()),
2550 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()),
2551 nptbins,ptmin,ptmax,200,-10,10);
2557 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxE%s",add[tm].
Data()),
2558 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()),
2559 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2567 (Form(
"hEnergyTime1Cell%s",add[tm].
Data()),
2568 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
2569 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2575 (Form(
"hEnergyTimeExotic%s",add[tm].
Data()),
2576 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
2577 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2583 (Form(
"hEnergyTimeTCardCorrNoSelection1Cell%s",add[tm].
Data()),
2584 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
2585 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2591 (Form(
"hEnergyTimeTCardCorrNoSelectionExotic%s",add[tm].
Data()),
2592 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
2593 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2599 (Form(
"hColRowExoticLowE1Cell%s",add[tm].
Data()),
2600 Form(
"column vs row, 1 cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2601 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2607 (Form(
"hColRowExoticHighE1Cell%s",add[tm].
Data()),
2608 Form(
"column vs row, 1 cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2609 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2615 (Form(
"hColRowExoticLowE%s",add[tm].
Data()),
2616 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2617 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2623 (Form(
"hColRowExoticHighE%s",add[tm].
Data()),
2624 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2625 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2631 (Form(
"hColRowExotic2ndCellDiffLowE%s",add[tm].
Data()),
2632 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2633 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2639 (Form(
"hColRowExotic2ndCellDiffHighE%s",add[tm].
Data()),
2640 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2641 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2647 (Form(
"hColRowExotic2ndCellSameLowE%s",add[tm].
Data()),
2648 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2649 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2655 (Form(
"hColRowExotic2ndCellSameHighE%s",add[tm].
Data()),
2656 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2657 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2663 (Form(
"hColRowTCardCorrNoSelectionExoticLowE%s",add[tm].
Data()),
2664 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2665 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2671 (Form(
"hColRowTCardCorrNoSelectionExoticHighE%s",add[tm].
Data()),
2672 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2673 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2679 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffLowE%s",add[tm].
Data()),
2680 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2681 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2687 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffHighE%s",add[tm].
Data()),
2688 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2689 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2695 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameLowE%s",add[tm].
Data()),
2696 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2697 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2703 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameHighE%s",add[tm].
Data()),
2704 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2705 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2711 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameLowE%s",add[tm].
Data()),
2712 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2713 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2719 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameHighE%s",add[tm].
Data()),
2720 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2721 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2727 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffLowE%s",add[tm].
Data()),
2728 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2729 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2735 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffHighE%s",add[tm].
Data()),
2736 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2737 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2743 (Form(
"hNCellsTCardSameAndDiffFraction_Exotic%s",add[tm].
Data()),
2744 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, exo > 0.97 %s",add[tm].
Data()),
2745 nptbins,ptmin,ptmax,100,0,1);
2751 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterEExo%s",add[tm].
Data()),
2752 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()),
2753 nptbins,ptmin,ptmax,200,-10,10);
2759 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterEExo%s",add[tm].
Data()),
2760 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()),
2761 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2767 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxEExo%s",add[tm].
Data()),
2768 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()),
2769 nptbins,ptmin,ptmax,200,-10,10);
2775 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxEExo%s",add[tm].
Data()),
2776 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()),
2777 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2783 for(
Int_t i = 0; i < 6; i++)
2785 for(
Int_t j = 0; j < 6; j++)
2788 (Form(
"hLambda0TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2789 Form(
"#lambda^{2}_{0} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
2790 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2796 (Form(
"hLambda1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2797 Form(
"#lambda^{2}_{1} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
2798 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2804 (Form(
"hLambda0NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2805 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=1, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
2806 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2812 (Form(
"hLambda1NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2813 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=1, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
2814 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2820 (Form(
"hLambda0NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2821 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=2, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
2822 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2828 (Form(
"hLambda1NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2829 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=2, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
2830 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2845 (Form(
"hNLocMaxTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2846 Form(
"#it{n}_{LM} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
2847 nptbins,ptmin,ptmax,10,0,10);
2853 (Form(
"hEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2854 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()),
2855 nptbins,ptmin,ptmax,100,0,1);
2861 (Form(
"hEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2862 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()),
2863 nptbins,ptmin,ptmax,100,0,1);
2869 (Form(
"hEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2870 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()),
2871 nptbins,ptmin,ptmax,100,0,1);
2877 (Form(
"hE2ndRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2878 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()),
2879 nptbins,ptmin,ptmax,100,0,1);
2885 (Form(
"hE2ndRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2886 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()),
2887 nptbins,ptmin,ptmax,100,0,1);
2893 (Form(
"hE2ndRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2894 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()),
2895 nptbins,ptmin,ptmax,100,0,1);
2901 (Form(
"hE2ndEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2902 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()),
2903 nptbins,ptmin,ptmax,100,0,1);
2909 (Form(
"hE2ndEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2910 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()),
2911 nptbins,ptmin,ptmax,100,0,1);
2917 (Form(
"hE2ndEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2918 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()),
2919 nptbins,ptmin,ptmax,100,0,1);
2926 (Form(
"hECellClusRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2927 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()),
2928 nptbins,ptmin,ptmax,100,0,1);
2934 (Form(
"hECellClusRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2935 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()),
2936 nptbins,ptmin,ptmax,100,0,1);
2942 (Form(
"hECellClusRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2943 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()),
2944 nptbins,ptmin,ptmax,100,0,1);
2950 (Form(
"hLogECellNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2951 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()),
2952 nptbins,ptmin,ptmax,150,-3,3);
2958 (Form(
"hLogECellNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2959 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()),
2960 nptbins,ptmin,ptmax,150,-3,3);
2966 (Form(
"hLogECellNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2967 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()),
2968 nptbins,ptmin,ptmax,150,-3,3);
2974 (Form(
"hECellWeightNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2975 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()),
2976 nptbins,ptmin,ptmax,90,0,4.5);
2982 (Form(
"hECellWeightNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2983 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()),
2984 nptbins,ptmin,ptmax,90,0,4.5);
2990 (Form(
"hECellWeightNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2991 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()),
2992 nptbins,ptmin,ptmax,90,0,4.5);
2998 (Form(
"hMassEClusTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2999 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()),
3000 nptbins,ptmin,ptmax,nmassbins,massmin,massmax);
3014 (Form(
"hTimeDiffTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3015 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()),
3016 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
3022 (Form(
"hColRowTCardCorrelNCellLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3023 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3024 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3030 (Form(
"hColRowTCardCorrelNCellHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3031 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3032 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3040 (Form(
"hExoticTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3041 Form(
"exoticity vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3042 nptbins,ptmin,ptmax,200,-1,1);
3048 (Form(
"hTimeDiffExoTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3049 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()),
3050 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
3056 (Form(
"hColRowTCardCorrelNCellExoticLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3057 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3058 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3064 (Form(
"hColRowTCardCorrelNCellExoticHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3065 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3066 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3281 (Form(
"hLambda0Exoticity_EBin%d%s",i,add[tm].
Data()),
3282 Form(
"#lambda^{2}_{0} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3283 200,-1,1,ssbins,ssmin,ssmax);
3289 (Form(
"hLambda1Exoticity_EBin%d%s",i,add[tm].
Data()),
3290 Form(
"#lambda^{2}_{1} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3291 200,-1,1,ssbins,ssmin,ssmax);
3305 (Form(
"hNCellsExoticity_EBin%d%s",i,add[tm].
Data()),
3306 Form(
"#it{n}_{cells} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3307 200,-1,1,nceclbins,nceclmin,nceclmax);
3308 fhNCellsExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
3313 (Form(
"hTimeExoticity_EBin%d%s",i,add[tm].
Data()),
3315 200,-1,1,100,-25,25);
3316 fhTimeExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
3321 (Form(
"hNCellsTCardSameAndDiff_Exotic_EBin%d%s",i,add[tm].
Data()),
3322 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()),
3323 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
3363 (Form(
"hNCellsTCardSameAndDiff_EBin%d%s",i,add[tm].
Data()),
3364 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()),
3365 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
3371 (Form(
"hLambda0Lambda1_EBin%d%s",i,add[tm].
Data()),
3372 Form(
"#lambda^{2}_{0} vs #lambda^{2}_{1}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3373 ssbins,ssmin,ssmax,ssbins,ssmin,ssmax);
3389 for(
Int_t j = 0; j < 6; j++)
3391 for(
Int_t k = 0; k < 6; k++)
3394 (Form(
"hLambda0Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
3395 Form(
"#lambda^{2}_{0} vs #it{exoticity}, #it{n}_{cell} TCard same = %d, diff =%d, #it{E}>8 GeV %s",j,k,add[tm].
Data()),
3396 200,-1,1,ssbins,ssmin,ssmax);
3402 (Form(
"hLambda1Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
3403 Form(
"#lambda^{2}_{1} vs #it{exoticity}, #it{n}_{cell} TCard same = %d, diff =%d, #it{E}>8 GeV %s",j,k,add[tm].
Data()),
3404 200,-1,1,ssbins,ssmin,ssmax);
3501 for(
Int_t i = 0; i < 12; i++)
3504 (Form(
"hTCardCorrECellMaxDiff_Case%d%s",i,add[tm].
Data()),
3505 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()),
3506 nptbins,ptmin,ptmax,210,-1,20);
3512 (Form(
"hTCardCorrEClusterDiff_Case%d%s",i,add[tm].
Data()),
3513 Form(
"#it{E}_{cluster}-#it{E}_{cell} vs #it{E}_{cluster}, for (un)correlated cells in TCard, case %d %s",i,add[tm].
Data()),
3514 nptbins,ptmin,ptmax,210,-1,20);
3536 (Form(
"hTCardCorrTCellMaxDiff_Case%d%s",i,add[tm].
Data()),
3537 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()),
3538 nptbins,ptmin,ptmax,1000,-100,100);
3546 (Form(
"hTCardCorrECellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
3547 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()),
3548 nptbins,ptmin,ptmax,210,-1,20);
3554 (Form(
"hTCardCorrEClusterDiffExo_Case%d%s",i,add[tm].
Data()),
3555 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()),
3556 nptbins,ptmin,ptmax,210,-1,20);
3578 (Form(
"hTCardCorrTCellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
3579 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()),
3580 nptbins,ptmin,ptmax,1000,-100,100);
3591 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3597 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3603 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3609 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3615 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3621 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3627 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3633 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3641 (Form(
"hTMPhiResidual_EBin%d",i),
3642 Form(
"#Delta #varphi_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3643 200,-1,1,nresphibins,resphimin,resphimax);
3649 (Form(
"hTMEtaResidual_EBin%d",i),
3650 Form(
"#Delta #eta_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3651 200,-1,1,nresetabins,resetamin,resetamax);
3657 (Form(
"hTMPhiResidual_LooseCut_EBin%d",i),
3658 Form(
"#Delta #varphi_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3659 200,-1,1,nresphibins,resphimin,resphimax);
3665 (Form(
"hTMEtaResidual_LooseCut_EBin%d",i),
3666 Form(
"#Delta #eta_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3667 200,-1,1,nresetabins,resetamin,resetamax);
3695 (Form(
"hColRowM02"),
3697 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax,40,0.,2.);
3704 (Form(
"hColRowM02NCellCut"),
3705 Form(
"column vs row vs M02, %2.2f<#it{E}<%2.2f GeV #it{n}_{cells}^{w>0.01} > %d for Neutral",
3707 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax,40,0.,2.);
3715 Form(
"%2.2f<#it{E}_{1}<%2.2f GeV, %2.2f<#it{E}_{2}<%2.2f GeV, %2.2f<#lambda^{2}_{0}<%2.2f"
3716 "#it{M}_{#gamma #gamma} vs #it{n}_{1, cells}^{w>0.01} vs SM number trig cluster",
3725 (
"hInvMassNCellSMSame",
3726 Form(
"%2.2f<#it{E}_{1}<%2.2f GeV, %2.2f<#it{E}_{2}<%2.2f GeV, %2.2f<#lambda^{2}_{0}<%2.2f"
3727 "#it{M}_{#gamma #gamma} vs #it{n}_{1, cells}^{w>0.01} vs SM number both cluster",
3738 (
"hEMaxCellTimeM02SM",
3739 Form(
"time vs #lambda_{0}^{2} vs SM number, "
3748 (
"hEMaxCellTimeNCellSM",
3749 Form(
"time vs #it{n}_{cells}^{w>0.01} vs SM number, "
3758 (
"hESecCellTimeNCellSM",
3759 Form(
"secondary cell time vs #it{n}_{cells}^{w>0.01} vs SM number, "
3769 Int_t nbinsdeltacells = 19 ;
3772 for(
Int_t col = 0; col < 2; col++)
3775 (Form(
"hColRowFromCellMaxLowM02_Col%d",col),
3776 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #it{n}_{cells}^{w>0.01} for 0.1 < #lambda_{0}^{2} < 0.3, "
3778 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells);
3784 (Form(
"hColRowFromCellMaxHighM02_Col%d",col),
3785 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #it{n}_{cells}^{w>0.01} for 0.5 < #lambda_{0}^{2} < 2, "
3787 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells);
3799 (Form(
"hNCellsPerClusterM02M20_SM%d",i),
3800 Form(
" vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01} vs #lambda_{1}^{2}, "
3802 (
Int_t)nShShBins/1.5,0,maxShSh/1.5, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3809 (Form(
"hNCellsPerClusterM02NLM_SM%d",i),
3810 Form(
" vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01} vs N_{lm}, "
3812 10,0,10, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3819 (Form(
"hEMaxCellEClusterM02NCell_SM%d",i),
3820 Form(
"(#it{E}_{cluster} - #it{E}_{max cell})/#it{E}_{cluster} vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01}, "
3822 50,0,1., cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3829 (Form(
"hEMaxCellLogM02NCell_SM%d",i),
3830 Form(
"log(#it{E}_{max cell}) vs #lambda_{0}^{2} vs SM number, "
3832 150,-3,3, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3839 (Form(
"hESecCellEMaxCellM02NCell_SM%d",i),
3840 Form(
"#it{E}_{cell}/#it{E}_{cell max} vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01}, "
3842 50,0,1., cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3849 (Form(
"hESecCellEClusterM02NCell_SM%d",i),
3850 Form(
"(#it{E}_{cluster} - #it{E}_{cell})/#it{E}_{cluster} vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01}, "
3852 50,0,1., cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3859 (Form(
"hESecCellLogM02NCell_SM%d",i),
3860 Form(
"log(#it{E}_{cell}) vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01}, "
3862 150,-3,3, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3869 (Form(
"hEMaxESecCellNCellLowM02_SM%d",i),
3870 Form(
"#it{E}_{cell}^{max} vs #it{E}_{cell}^{secondary} vs #it{n}_{cells}^{w>0.01} vs 0.1 < #lambda_{0}^{2} < 0.3, "
3872 24,0,12, 48,0,12,cellBins,cellMin,cellMax);
3879 (Form(
"hEMaxECrossNCellLowM02_SM%d",i),
3880 Form(
"#it{E}_{cell}^{max} vs exoticity vs #it{n}_{cells}^{w>0.01} vs 0.1 < #lambda_{0}^{2} < 0.3, "
3882 24,0,12, 40,0.6,1.,cellBins,cellMin,cellMax);
3889 (Form(
"hEMaxESecCellNCellHighM02_SM%d",i),
3890 Form(
"#it{E}_{cell}^{max} vs #it{E}_{cell}^{secondary} vs #it{n}_{cells}^{w>0.01} vs 0.5 < #lambda_{0}^{2} < 2, "
3892 24,0,12, 48,0,12,cellBins,cellMin,cellMax);
3899 (Form(
"hEMaxECrossNCellHighM02_SM%d",i),
3900 Form(
"#it{E}_{cell}^{max} vs exoticity vs #it{n}_{cells}^{w>0.01} for 0.5 < #lambda_{0}^{2} < 2, "
3902 24,0,12, 40,0.6,1.,cellBins,cellMin,cellMax);
3908 Int_t nbinsdeltacells = 19 ;
3911 for(
Int_t col = 0; col < 2; col++)
3914 (Form(
"hColRowFromCellMaxLowM02_SM%d_Col%d",i,col),
3915 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #it{n}_{cells}^{w>0.01} for 0.1 < #lambda_{0}^{2} < 0.3, "
3917 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,cellBins,cellMin,cellMax);
3924 (Form(
"hColRowFromCellMaxHighM02_SM%d_Col%d",i,col),
3925 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #it{n}_{cells}^{w>0.01} for 0.5 < #lambda_{0}^{2} < 2, "
3927 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,cellBins,cellMin,cellMax);
3933 for(
Int_t j = 0; j < 3; j++)
3936 (Form(
"hColRowFromCellMaxEMaxSecDiffLowM02_SM%d_Col%d_NCellBin%d",i,col,j),
3937 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #Delta #it{E}_{max-secondary} for 0.1 < #lambda_{0}^{2} < 0.3, "
3938 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
3939 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,120,0,12);
3946 (Form(
"hColRowFromCellMaxEMaxSecDiffHighM02_SM%d_Col%d_NCellBin%d",i,col,j),
3947 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #Delta #it{E}_{max-secondary} for 0.5 < #lambda_{0}^{2} < 2, "
3948 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
3949 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,120,0,12);
3956 (Form(
"hColRowFromCellMaxEMaxSecDiffFracLowM02_SM%d_Col%d_NCellBin%d",i,col,j),
3957 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #Delta #it{E}_{max-secondary}/#it{E}_{max} for 0.1 < #lambda_{0}^{2} < 0.3, "
3958 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
3959 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,40,0,1);
3966 (Form(
"hColRowFromCellMaxEMaxSecDiffFracHighM02_SM%d_Col%d_NCellBin%d",i,col,j),
3967 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #Delta #it{E}_{max-secondary}/#it{E}_{max} for 0.5 < #lambda_{0}^{2} < 2, "
3968 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
3969 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,40,0,1);
3976 (Form(
"hColRowFromCellMaxECellClusterRatLowM02_SM%d_Col%d_NCellBin%d",i,col,j),
3977 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #it{E}_{cell}/#it{E}_{cluster} for 0.1 < #lambda_{0}^{2} < 0.3, "
3978 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
3979 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,50,0,1);
3986 (Form(
"hColRowFromCellMaxECellClusterRatHighM02_SM%d_Col%d_NCellBin%d",i,col,j),
3987 Form(
"cell col_{max}-#col_{secondary} vs cell #row_{max}-#row_{secondary} vs #it{E}_{cell}/#it{E}_{cluster} for 0.5 < #lambda_{0}^{2} < 2, "
3988 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
3989 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,50,0,1);
4000 for(
Int_t imatch = 0; imatch < 3; imatch++)
4005 (Form(
"hDeltaIEtaDeltaIPhi_%s",matchCase[imatch].
Data()),
4006 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),
4007 nEbins,minE,maxE,cellBins,cellMin,cellMax,cellBins,cellMin,cellMax);
4014 (Form(
"hDeltaIEtaDeltaIPhiTot_%s",matchCase[imatch].
Data()),
4015 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),
4016 nEbins,minE,maxE,cellBins,cellMin,cellMax,cellBins,cellMin,cellMax);
4023 (Form(
"hDeltaIA_%s",matchCase[imatch].
Data()),
4024 Form(
"Cluster *asymmetry* in cell units vs E for %s",matchCase[imatch].
Data()),
4025 nptbins,ptmin,ptmax,asyBins,asyMin,asyMax);
4026 fhDeltaIA[imatch]->SetXTitle(
"#it{E}_{cluster}");
4027 fhDeltaIA[imatch]->SetYTitle(
"#it{A}_{cell in cluster}");
4028 outputContainer->Add(
fhDeltaIA[imatch]);
4031 (Form(
"hDeltaIATot_%s",matchCase[imatch].
Data()),
4032 Form(
"Cluster *total asymmetry* in cell units vs E for %s",matchCase[imatch].
Data()),
4033 nptbins,ptmin,ptmax,asyBins,asyMin,asyMax);
4035 fhDeltaIATot[imatch]->SetYTitle(
"#it{A}_{cell in cluster}^{total}");
4039 (Form(
"hDeltaIAM02_%s",matchCase[imatch].
Data()),
4040 Form(
"Cluster *asymmetry* in cell units vs #lambda^{2}_{0} for %s",matchCase[imatch].
Data()),
4041 nEbins,minE,maxE,nShShBins,minShSh,maxShSh,asyBins,asyMin,asyMax);
4044 fhDeltaIAM02[imatch]->SetZTitle(
"#it{A}_{cell in cluster}");
4048 (Form(
"hDeltaIATotM02_%s",matchCase[imatch].
Data()),
4049 Form(
"Cluster *total asymmetry* in cell units vs #lambda^{2}_{0} for %s",matchCase[imatch].
Data()),
4050 nEbins,minE,maxE,nShShBins,minShSh,maxShSh,asyBins,asyMin,asyMax);
4053 fhDeltaIATotM02[imatch]->SetZTitle(
"#it{A}_{cell in cluster}^{total}");
4057 (Form(
"hDeltaIAM20_%s",matchCase[imatch].
Data()),
4058 Form(
"Cluster *asymmetry* in cell units vs #lambda^{2}_{1} for %s",matchCase[imatch].
Data()),
4059 nEbins,minE,maxE,(
Int_t)nShShBins/1.5,0,maxShSh/1.5,asyBins,asyMin,asyMax);
4062 fhDeltaIAM20[imatch]->SetZTitle(
"#it{A}_{cell in cluster}");
4066 (Form(
"hDeltaIATotM20_%s",matchCase[imatch].
Data()),
4067 Form(
"Cluster *total asymmetry* in cell units vs #lambda^{2}_{1} for %s",matchCase[imatch].
Data()),
4068 nEbins,minE,maxE,(
Int_t)nShShBins/1.5,0,maxShSh/1.5,asyBins,asyMin,asyMax);
4071 fhDeltaIATotM20[imatch]->SetZTitle(
"#it{A}_{cell in cluster}^{total}");
4075 (Form(
"hDeltaIANCells_%s",matchCase[imatch].
Data()),
4076 Form(
"Cluster *asymmetry* in cell units vs N cells in cluster for %s",matchCase[imatch].
Data()),
4077 nEbins,minE,maxE,cellBins,cellMin,cellMax,asyBins,asyMin,asyMax);
4084 (Form(
"hDeltaIATotNCells_%s",matchCase[imatch].
Data()),
4085 Form(
"Cluster *total asymmetry* in cell units vs N cells in cluster for %s",matchCase[imatch].
Data()),
4086 nEbins,minE,maxE,cellBins,cellMin,cellMax,asyBins,asyMin,asyMax);
4095 (Form(
"hDeltaIAOrigin_%s",matchCase[imatch].
Data()),
4096 Form(
"Cluster *asymmetry* in cell units vs E vs primary ID, for ID %s",matchCase[imatch].
Data()),
4097 nEbins,minE,maxE,10,-0.5,9.5,asyBins,asyMin,asyMax);
4101 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4102 fhDeltaIAOrigin[imatch]->GetYaxis()->SetBinLabel(ilabel,mcParticleStringLabel[ilabel-1]);
4106 (Form(
"hDeltaIATotOrigin_%s",matchCase[imatch].
Data()),
4107 Form(
"Cluster *total asymmetry* in cell units vs E vs primary ID, for ID %s",matchCase[imatch].
Data()),
4108 nEbins,minE,maxE,10,-0.5,9.5,asyBins,asyMin,asyMax);
4112 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4113 fhDeltaIATotOrigin[imatch]->GetYaxis()->SetBinLabel(ilabel,mcParticleStringLabel[ilabel-1]);
4119 (Form(
"hClusterMaxCellDiffM02_%s",matchCase[imatch].
Data()),
4120 Form(
"#it{E}_{cluster} vs (#it{E}_{cluster} - #it{E}_{cell max})/#it{E}_{cluster} vs #lambda_{0}^{2} for ID %s",matchCase[imatch].
Data()),
4121 nEbins,minE,maxE, 20,0,1.,nShShBins,minShSh,maxShSh);
4128 (Form(
"hClusterTimeEnergyM02_%s",matchCase[imatch].
Data()),
4129 Form(
"#it{E} vs TOF vs #lambda_{0}^{2} for ID %s",matchCase[imatch].
Data()),
4130 nEbins,minE,maxE,45,-25.5,20.5,nShShBins,minShSh,maxShSh);
4137 (Form(
"hNCellsPerClusterM02_%s",matchCase[imatch].
Data()),
4138 Form(
"#it{E} vs #it{n}_{cells} vs #lambda_{0}^{2} for ID %s",matchCase[imatch].
Data()),
4139 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4146 (Form(
"hNCellsPerClusterM20_%s",matchCase[imatch].
Data()),
4147 Form(
"#it{E} vs #it{n}_{cells} vs #lambda_{1}^{2} for ID %s",matchCase[imatch].
Data()),
4148 nEbins,minE,maxE,cellBins,cellMin,cellMax,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4155 (Form(
"hNCellsPerClusterNLM_%s",matchCase[imatch].
Data()),
4156 Form(
"#it{E} vs #it{n}_{cells} vs n_{lm} for ID %s",matchCase[imatch].
Data()),
4157 nEbins,minE,maxE,cellBins,cellMin,cellMax,10,0,10);
4164 (Form(
"hSMM02_%s",matchCase[imatch].
Data()),
4165 Form(
"#it{E} vs SM number vs #lambda_{0}^{2}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4166 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,minShSh,maxShSh);
4167 fhSMM02[imatch]->SetXTitle(
"#it{E} (GeV)");
4168 fhSMM02[imatch]->SetYTitle(
"SM number");
4169 fhSMM02[imatch]->SetZTitle(
"#lambda_{0}^{2}");
4170 outputContainer->Add(
fhSMM02[imatch]);
4173 (Form(
"hSMM02NoCut_%s",matchCase[imatch].
Data()),
4174 Form(
"#it{E} vs SM number vs #lambda_{0}^{2} for ID %s",matchCase[imatch].
Data()),
4175 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,minShSh,maxShSh);
4182 (Form(
"hSMM20LowM02_%s",matchCase[imatch].
Data()),
4183 Form(
"#it{E} vs SM number vs #lambda_{1}^{2}, #it{n}_{cells}^{w>0.01}>%d, for 0.1<#lambda_{0}^{2}<0.3, for ID %s",
4185 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4192 (Form(
"hSMM20LowM02NoCut_%s",matchCase[imatch].
Data()),
4193 Form(
"#it{E} vs SM number vs #lambda_{1}^{2}, for 0.1<#lambda_{0}^{2}<0.3, for ID %s",matchCase[imatch].
Data()),
4194 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4201 (Form(
"hSMM20HighM02_%s",matchCase[imatch].
Data()),
4202 Form(
"#it{E} vs SM number vs #lambda_{1}^{2}, #it{n}_{cells}^{w>0.01}>%d, for 0.5<#lambda_{0}^{2}<2, for ID %s",
4204 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4211 (Form(
"hSMM20HighM02NoCut_%s",matchCase[imatch].
Data()),
4212 Form(
"#it{E} vs SM number vs #lambda_{1}^{2}, for 0.5<#lambda_{0}^{2}<2, for ID %s",matchCase[imatch].
Data()),
4213 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4220 (Form(
"hSMNCell_%s",matchCase[imatch].
Data()),
4221 Form(
"#it{E} vs SM number vs #it{n}_{cells}^{w>0.01} for ID %s",matchCase[imatch].
Data()),
4222 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,cellBins,cellMin,cellMax);
4223 fhSMNCell[imatch]->SetXTitle(
"#it{E} (GeV)");
4224 fhSMNCell[imatch]->SetYTitle(
"SM number");
4225 fhSMNCell[imatch]->SetZTitle(
"#it{n}_{cells}^{w>0.01}");
4226 outputContainer->Add(
fhSMNCell[imatch]);
4229 (Form(
"hSMNCellM02_%s",matchCase[imatch].
Data()),
4230 Form(
"SM number vs #it{n}_{cells}^{w>0.01} vs #lambda_{0}^{2}, "
4232 fNModules,-0.5,fNModules-0.5,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4235 fhSMNCellM02[imatch]->SetYTitle(
"#it{n}_{cells}^{w>0.01}");
4239 (Form(
"hColM02_%s",matchCase[imatch].
Data()),
4240 Form(
"#it{E} vs column number vs #lambda_{0}^{2}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4241 nEbins,minE,maxE,48,-0.5,47.5,nShShBins,minShSh,maxShSh);
4242 fhColM02[imatch]->SetXTitle(
"#it{E} (GeV)");
4243 fhColM02[imatch]->SetYTitle(
"column number");
4244 fhColM02[imatch]->SetZTitle(
"#lambda_{0}^{2}");
4245 outputContainer->Add(
fhColM02[imatch]);
4248 (Form(
"hRowM02_%s",matchCase[imatch].
Data()),
4249 Form(
"#it{E} vs row number vs #lambda_{0}^{2}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4250 nEbins,minE,maxE,24,-0.5,23.5,nShShBins,minShSh,maxShSh);
4251 fhRowM02[imatch]->SetXTitle(
"#it{E} (GeV)");
4252 fhRowM02[imatch]->SetYTitle(
"row number");
4253 fhRowM02[imatch]->SetZTitle(
"#lambda_{0}^{2}");
4254 outputContainer->Add(
fhRowM02[imatch]);
4259 (Form(
"hNCellsPerClusterMEta_%s",matchCase[imatch].
Data()),
4260 Form(
"#it{E} vs #it{n}_{cells} vs #sigma_{#eta}^{2} for ID %s",matchCase[imatch].
Data()),
4261 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4268 (Form(
"hNCellsPerClusterMPhi_%s",matchCase[imatch].
Data()),
4269 Form(
"#it{E} vs #it{n}_{cells} vs #sigma_{#varphi}^{2} for ID %s",matchCase[imatch].
Data()),
4270 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4277 (Form(
"hNCellsPerClusterMEtaPhi_%s",matchCase[imatch].
Data()),
4278 Form(
"#it{E} vs #it{n}_{cells} vs #sigma_{#eta#varphi}^{2} for ID %s",matchCase[imatch].
Data()),
4279 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,-1*maxShSh,maxShSh);
4286 (Form(
"hNCellsPerClusterMEtaPhiA_%s",matchCase[imatch].
Data()),
4287 Form(
"#it{E} vs #it{n}_{cells} vs (#sigma_{#varphi}^{2}-#sigma_{#eta}^{2})/(#sigma_{#varphi}^{2}+#sigma_{#eta}^{2}) for ID %s",matchCase[imatch].
Data()),
4288 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,-1*maxShSh,maxShSh);
4291 fhNCellsPerClusterMEtaPhiA[imatch]->SetZTitle(
"(#sigma_{#varphi}^{2}-#sigma_{#eta}^{2})/(#sigma_{#varphi}^{2}+#sigma_{#eta}^{2})");
4298 (Form(
"hOrigin_%s",matchCase[imatch].
Data()),
4299 Form(
"#it{E} vs origin for ID %s",matchCase[imatch].
Data()),
4300 nEbins,minE,maxE, 10,-0.5,9.5);
4301 fhOriginE[imatch]->SetXTitle(
"#it{E} (GeV)");
4302 fhOriginE[imatch]->SetYTitle(
"particle");
4303 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4304 fhOriginE[imatch]->GetYaxis()->SetBinLabel(ilabel ,mcParticleStringLabel[ilabel-1]);
4305 outputContainer->Add(
fhOriginE[imatch]);
4308 (Form(
"hOriginM02_%s",matchCase[imatch].
Data()),
4309 Form(
"#it{E} vs origin vs #lambda_{0}^{2} for ID %s",matchCase[imatch].
Data()),
4310 nEbins,minE,maxE,10,-0.5,9.5,nShShBins,minShSh,maxShSh);
4313 fhOriginM02[imatch]->SetZTitle(
"#lambda_{0}^{2}");
4314 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4315 fhOriginM02[imatch]->GetYaxis()->SetBinLabel(ilabel,mcParticleStringLabel[ilabel-1]);
4330 (
"hClusterMaxCellCloseCellRatioM02",
"#it{E} vs #it{E}_{cell-i}/#it{E}_{cell max} vs #lambda_{0}^{2}",
4331 nEbins,minE,maxE, 20,0,1.,nShShBins,minShSh,maxShSh);
4349 (
"hClusterMaxCellECrossM02",
4350 "#it{E} vs exoticity vs #lambda_{0}^{2}",
4351 nEbins,minE,maxE, 40,0.6,1.,nShShBins,minShSh,maxShSh);
4362 nptbins,ptmin,ptmax, 100,0,1.);
4368 nptbins,ptmin,ptmax, 100,-10,0);
4374 nptbins,ptmin,ptmax, 100,0,1.);
4380 nptbins,ptmin,ptmax, 100,-10,0);
4386 nptbins*2,ptmin,ptmax*2, 100,0,1.);
4392 nptbins*2,ptmin,ptmax*2, 100,-10,0);
4402 if(imod < fFirstModule || imod >
fLastModule)
continue;
4405 Form(
"#cell energy / sum all energy vs all energy in Module %d",imod),
4406 nptbins*2,ptmin,ptmax*2, 100,0,1.);
4412 Form(
"Log(cell energy / sum all energy) vs all energy in Module %d",imod),
4413 nptbins*2,ptmin,ptmax*2, 100,-10,0);
4420 Float_t cellEmin [] = {0.05,0.1,0.15,0.2};
4421 Float_t cellTmin [] = {50.,100.,1000000.};
4423 for(
Int_t iw = 0; iw < 12; iw++)
4426 for(
Int_t iEmin = 0; iEmin < 4; iEmin++)
4428 for(
Int_t iTmin = 0; iTmin < 3; iTmin++)
4431 Form(
"#lambda^{2}_{0} vs E, w0=%1.2f, cell E>%2.2f MeV, |t|<%2.0f ns",
4432 w0, cellEmin[iEmin], cellTmin[iTmin]),
4433 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4448 Form(
"#lambda^{2}_{0} vs E, w0=%1.2f, cell E>%2.2f MeV, |t|<%2.0f ns, |#eta| < 0.15",
4449 w0, cellEmin[iEmin], cellTmin[iTmin]),
4450 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4459 TString mcnames[] = {
"Photon",
"Electron",
"Conversion",
"Pi0",
"Hadron"};
4460 for(
Int_t imc = 0; imc < 5; imc++)
4463 Form(
"shower shape, #lambda^{2}_{0} vs E, w0 = %1.1f, for MC %s",w0,mcnames[imc].
Data()),
4464 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4483 return outputContainer;
4544 AliInfo(Form(
"dEdX cuts init for MC: %d<dedx_ele<%d; %d<dedx_had<%d",
4552 if ( runNumber < 146861 )
4559 else if ( runNumber < 156000 )
4566 else if ( runNumber < 165000 )
4581 AliInfo(Form(
"dEdX cuts init for run %d: %d<dedx_ele<%d; %d<dedx_had<%d",
4595 Int_t nCellsPerCluster)
4613 if(ampMax < 0.01)
return kFALSE;
4629 printf(
"**** Print %s %s ****\n", GetName(), GetTitle() ) ;
4633 printf(
"Min n cells : %d\n" ,
fNCellMin) ;
4635 printf(
"Min M02 : %1.2f\n",
fM02Min) ;
4651 AliDebug(1,
"Start");
4667 AliFatal(Form(
"AliAnaClusterShapeCorrelStudies::MakeAnalysisFillHistograms() - Wrong calorimeter name <%s>, END",
GetCalorimeterString().
Data()));
4671 AliFatal(Form(
"AliAnaClusterShapeCorrelStudies::MakeAnalysisFillHistograms() - No CaloClusters or CaloCells available"));
4701 if(clus->GetNCells() < 3)
return ;
4703 for (
Int_t ipos = 0; ipos < clus->GetNCells(); ipos++)
4705 Int_t id = clus->GetCellsAbsId()[ipos];
4720 AliWarning(Form(
"Wrong calculated energy %f",energy));
4731 for (
Int_t ipos = 0; ipos < clus->GetNCells(); ipos++)
4733 Int_t id = clus->GetCellsAbsId()[ipos];
4748 Float_t l0org = clus->GetM02();
4749 Float_t l1org = clus->GetM20();
4750 Float_t dorg = clus->GetDispersion();
4754 if(
IsDataMC() && clus->GetNLabels() > 0)
4779 Float_t cellEmin [] = {0.05,0.1,0.15,0.2};
4780 Float_t cellTmin [] = {50.,100.,1000000.};
4782 for(
Int_t iw = 0; iw < 12; iw++)
4784 for(
Int_t iEmin = 0; iEmin < 4; iEmin++)
4786 for(
Int_t iTmin = 0; iTmin < 3; iTmin++)
4792 Float_t disp = 0, dEta = 0, dPhi = 0;
4793 Float_t sEta = 0, sPhi = 0, sEtaPhi = 0;
4797 cellEmin[iEmin], cellTmin[iTmin], bc,
4798 newEnergy, l0, l1, disp, dEta, dPhi,
4799 sEta, sPhi, sEtaPhi);
4820 clus->SetM02(l0org);
4821 clus->SetM20(l1org);
4822 clus->SetDispersion(dorg);
TH2F * fhNCellsTCardSameAndDiffFraction[2]
! Cluster fraction of NCells in same TCard as leading vs energy
virtual Bool_t IsFiducialCutOn() const
Float_t GetHistoPtMax() const
TH2F * fhECellClusRatNLM1TCardCorrNoSelection[2]
! Cluster E cell / E cluster for NLM=1 vs E for clusters selected for TCard correlation studies ...
TH2F * fhLambda1NLM2TCardCorrNoSelection[2]
! Cluster m20 vs E for clusters selected for TCard correlation studies, nlm=2
TH2F * fhEMaxCellClusterRatio
! e max cell / e cluster vs e cluster
TH2F * fhLambda0NLM1TCardCorrelNCell[6][6][2]
! Cluster m02 vs E, cluster contains 0 to more than 6 cells with w > 0.01 in same TCard or diff TCard...
TH2F * fhE2ndRatNLM1TCardCorrelNCell[6][6][2]
! Cluster E cell second max / E cluster for NLM=1 vs E, cluster contains 0 to more than 6 cells with ...
void ClusterShapeHistograms(AliVCluster *cluster, Int_t absIdMax, Double_t maxCellFraction, Float_t eCrossFrac, Float_t ampMax, Double_t tmax, Int_t matchedPID, Int_t mcIndex)
TH2F * fhColRowTCardCorrNoSelectionExoticLowE[2]
! col-row cluster cell max for those selected for TCard correlation studies, 5 < E < 8 GeV...
TH2F * fhColRowTCardCorrelNCellExoticHighE[6][6][2]
! Cluster max cell col vs row, E > 8 GeV, cluster contains 0 to more than 6 cells with w > 0...
TH3F * fhClusterMaxCellCloseCellRatioM02
! Ratio between max cell energy and cell energy of the same cluster
Float_t fM02Min
Minimum M02 on clusters.
TH2F * fhECellTotalRatio
! e cell / e total vs e total
Int_t GetHistoNClusterCellMin() const
TH2F * fhExoticTCardCorrNoSelection[2]
! exoticity per cluster vs cluster energy, clusters selected for TCard correlation studies ...
TH2F * fhTMPhiResidualExoticity[14]
! Cluster-track matching residual in phi vs exoticity
Float_t GetHistoPtMin() const
TH2F * fhLogECellNLM2TCardCorrNoSelection[2]
! Cluster Log E cell for NLM=2 vs E for clusters selected for TCard correlation studies ...
TH2F * fhTimeTCardCorrNoSelection[2]
! Cluster time vs E for clusters selected for TCard correlation studies
TH2F * fhColRowTCardCorrelNCellExoticLowE[6][6][2]
! Cluster max cell col vs row, E > 2 GeV, cluster contains 0 to more than 6 cells with w > 0...
TH2F * fhEnergyTimeTCardCorrNoSelectionExotic[2]
! cluster energy vs time, exo > 0.97, T-Card strict cuts
Int_t fNMaxColsFull
Number of EMCAL/PHOS columns full detector.
TH2F * fhTCardCorrEClusterDiffExo[12][2]
! Cluster energy - secondary cell energy in cluster vs cluster energy, different secondary cell selec...
TH3F * fhClusterMaxCellDiffM02[3]
! Difference between cluster energy and energy of cell with more energy, vs m02
Int_t GetHistoShowerShapeBins() const
TH2F * fhSameRowDiffColAndTCardCellsTimeDiffCellMaxE[2]
! Secondary cell energy difference vs leading cell energy, one in same TCard as cell max...
virtual Int_t GetCalorimeter() const
TH2F * fhSameRowDiffColAndTCardCellsTimeDiffCellMaxEExo[2]
! Secondary cell energy difference vs leading cell energy, one in same TCard as cell max...
virtual void AddToHistogramsName(TString add)
virtual AliVCaloCells * GetEMCALCells() const
TH3F * fhNCellsPerClusterMEtaPhi[3]
! N cells per cluster vs cluster energy vs shape in eta*phi direction
TH3F * fhDeltaIAOrigin[3]
! Cluster "asymmetry" vs E vs origin
TH2F * fhECellClusRatNLM2TCardCorrelNCell[6][6][2]
! Cluster E cell / E cluster for NLM=2 vs E, cluster contains 0 to more than 6 cells with w > 0...
TH2F ** fhECellTotalRatioMod
! e cell / e total vs e total, per SM
TH3F * fhColRowM02
! Cluster position in cell max col-row, 8 < E < 12 GeV, neutral only
virtual Bool_t IsDataMC() const
TH3F * fhColRowFromCellMaxEMaxSecDiffLowM02PerSM[20][2][3]
! secondary cell distance to cell max in col vs row vs energy difference to cell max, (2-3), (4-5) (>5) n cells, for clusters 8 < E < 12 GeV, 0.1 < M02 < 0.3, per SM number, per odd/pair column
Int_t fNEBinCuts
Number of energy bin cuts for fStudyTCardCorrelation.
TH2F * fhNLocMaxTCardCorrelNCell[6][6][2]
! Cluster nlocmax vs E, cluster contains 0 to more than 6 cells with w > 0.01 in same TCard or diff T...
void WeightHistograms(AliVCluster *clus, Int_t mcTag)
TH3F * fhColM02[3]
! main cell column vs m02, n cell > 4
Int_t fLastModule
Last EMCAL/PHOS module, set in CaloUtils or depending fidutial cuts.