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++)
427 for(
Int_t j = 0; j < 3; j++)
471 Int_t icol = -1, irow = -1, iRCU = -1, icolAbs = -1, irowAbs = -1;
472 GetModuleNumberCellIndexesAbsCaloMap(absIdMax,GetCalorimeter(), icol, irow, iRCU, icolAbs, irowAbs);
476 for(
Int_t ie = 0; ie < fNEBinCuts; ie++)
478 if( energy >= fEBinCuts[ie] && energy < fEBinCuts[ie+1] ) ebin = ie;
481 Int_t ncells = clus->GetNCells();
482 Double_t time = clus->GetTOF()*1.e9;
483 if(time>400) time-=fConstantTimeShift;
485 Float_t deta = clus->GetTrackDz();
486 Float_t dphi = clus->GetTrackDx();
494 if(energy >= 5 && energy < 8)
495 fhColRowExoticLowE1Cell [matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
498 fhColRowExoticHighE1Cell[matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
500 if ( time > 5) fhColRowExoticHighE1CellPosTime->Fill(icolAbs,irowAbs,GetEventWeight());
501 else if( time < -5) fhColRowExoticHighE1CellNegTime->Fill(icolAbs,irowAbs,GetEventWeight());
502 else fhColRowExoticHighE1CellNulTime->Fill(icolAbs,irowAbs,GetEventWeight());
505 fhEnergyTime1Cell[matched]->Fill(energy,time,GetEventWeight());
506 fhEnergyTMEtaResidual1Cell->Fill(energy,deta,GetEventWeight());
507 fhEnergyTMPhiResidual1Cell->Fill(energy,dphi,GetEventWeight());
509 else if(exoticity > 0.97)
511 if(energy >= 5 && energy < 8)
512 fhColRowExoticLowE [matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
515 fhColRowExoticHighE[matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
517 if ( time > 5) fhColRowExoticHighEPosTime->Fill(icolAbs,irowAbs,GetEventWeight());
518 else if( time < -5) fhColRowExoticHighENegTime->Fill(icolAbs,irowAbs,GetEventWeight());
519 else fhColRowExoticHighENulTime->Fill(icolAbs,irowAbs,GetEventWeight());
522 fhEnergyTimeExotic[matched]->Fill(energy,time,GetEventWeight());
523 fhEnergyTMEtaResidualExotic->Fill(energy,deta,GetEventWeight());
524 fhEnergyTMPhiResidualExotic->Fill(energy,dphi,GetEventWeight());
526 for (
Int_t ipos = 0; ipos < ncells; ipos++)
528 Int_t absId = clus->GetCellsAbsId()[ipos];
530 Float_t eCell = fCaloCellList->GetCellAmplitude(absId);
531 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
534 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
536 if( absId == absIdMax || weight < 0.01 )
continue;
538 Int_t rowDiff = -100, colDiff = -100;
539 Bool_t sameTCard = GetCaloUtils()->IsAbsIDsFromTCard(absIdMax,absId,rowDiff,colDiff);
542 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
543 GetModuleNumberCellIndexesAbsCaloMap(absId,GetCalorimeter(), icol2, irow2, iRCU2, icolAbs2, irowAbs2);
547 if(energy >= 5 && energy < 8)
548 fhColRowExotic2ndCellDiffLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
550 fhColRowExotic2ndCellDiffHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
554 if(energy >= 5 && energy < 8)
555 fhColRowExotic2ndCellSameLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
557 fhColRowExotic2ndCellSameHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
562 else if ( energy > 8 )
564 if ( time > 5) fhColRowHighEPosTime->Fill(icolAbs,irowAbs,GetEventWeight());
565 else if( time < -5) fhColRowHighENegTime->Fill(icolAbs,irowAbs,GetEventWeight());
566 else fhColRowHighENulTime->Fill(icolAbs,irowAbs,GetEventWeight());
571 fhTMPhiResidualExoticityLooseCut[ebin]->Fill(exoticity,dphi);
572 fhTMEtaResidualExoticityLooseCut[ebin]->Fill(exoticity,deta);
575 else if ( energy > 8 && ncells > 1 && exoticity < 0.97)
577 if ( time > 5) fhColRowHighEPosTime->Fill(icolAbs,irowAbs,GetEventWeight());
578 else if( time < -5) fhColRowHighENegTime->Fill(icolAbs,irowAbs,GetEventWeight());
579 else fhColRowHighENulTime->Fill(icolAbs,irowAbs,GetEventWeight());
585 Int_t etaRegion = -1, phiRegion = -1;
586 GetCaloUtils()->GetEMCALSubregion(clus,fCaloCellList,etaRegion,phiRegion);
588 if ( etaRegion !=0 ) return ;
594 fhEnergyTimeTCardCorrNoSelection1Cell[matched]->Fill(energy,time,GetEventWeight());
595 fhEnergyTMEtaResidualTCardCorrNoSelection1Cell->Fill(energy,deta,GetEventWeight());
596 fhEnergyTMPhiResidualTCardCorrNoSelection1Cell->Fill(energy,dphi,GetEventWeight());
598 else if(exoticity > 0.97)
600 fhEnergyTimeTCardCorrNoSelectionExotic[matched]->Fill(energy,time,GetEventWeight());
601 fhEnergyTMEtaResidualTCardCorrNoSelectionExotic->Fill(energy,deta,GetEventWeight());
602 fhEnergyTMPhiResidualTCardCorrNoSelectionExotic->Fill(energy,dphi,GetEventWeight());
609 if ( m02 > 0.001 ) lamR = m20/m02;
611 Int_t absIdList[ncells];
613 Int_t nlm = GetCaloUtils()->GetNumberOfLocalMaxima(clus, fCaloCellList, absIdList, maxEList) ;
619 Int_t nCellWithWeight = 1;
627 Float_t eCellMax = fCaloCellList->GetCellAmplitude(absIdMax);
628 Double_t tCellMax = fCaloCellList->GetCellTime(absIdMax);
630 GetCaloUtils()->RecalibrateCellAmplitude(eCellMax, GetCalorimeter(), absIdMax);
631 GetCaloUtils()->RecalibrateCellTime(tCellMax, GetCalorimeter(), absIdMax, GetReader()->GetInputEvent()->GetBunchCrossNumber());
635 if(tCellMax>400) tCellMax-=fConstantTimeShift;
646 Bool_t sameTCard2ndMax = kFALSE;
650 Float_t eCellSameRowSameTCardNearCol = 0.;
651 Float_t eCellSameRowDiffTCardNearCol = 0.;
652 Double_t tCellSameRowSameTCardNearCol = 0.;
653 Double_t tCellSameRowDiffTCardNearCol = 0.;
660 for (
Int_t ipos = 0; ipos < ncells; ipos++)
662 Int_t absId = clus->GetCellsAbsId()[ipos];
664 Float_t eCell = fCaloCellList->GetCellAmplitude(absId);
665 Double_t tCell = fCaloCellList->GetCellTime(absId);
667 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
668 GetCaloUtils()->RecalibrateCellTime(tCell, GetCalorimeter(), absId, GetReader()->GetInputEvent()->GetBunchCrossNumber());
670 if(tCell>400) tCell-=fConstantTimeShift;
673 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
675 if( absId == absIdMax || weight < 0.01 )
continue;
679 fhECellClusRatNLM1TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
680 fhLogECellNLM1TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
681 fhECellWeightNLM1TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
685 fhECellClusRatNLM2TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
686 fhLogECellNLM2TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
687 fhECellWeightNLM2TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
691 fhECellClusRatNLM3TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
692 fhLogECellNLM3TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
693 fhECellWeightNLM3TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
696 Int_t rowDiff = -100, colDiff = -100;
697 Bool_t sameTCard = GetCaloUtils()->IsAbsIDsFromTCard(absIdMax,absId,rowDiff,colDiff);
703 fhECellSameClusRatNLM1TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
704 fhLogECellSameNLM1TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
705 fhECellSameWeightNLM1TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
709 fhECellSameClusRatNLM2TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
710 fhLogECellSameNLM2TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
711 fhECellSameWeightNLM2TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
715 fhECellSameClusRatNLM3TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
716 fhLogECellSameNLM3TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
717 fhECellSameWeightNLM3TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
727 Int_t indexType = -1;
732 if(TMath::Abs(rowDiff) == 1) nearRow = kTRUE;
733 if(TMath::Abs(colDiff) == 1) nearCol = kTRUE;
735 if ( rowDiff == 0 && colDiff != 0 )
737 if ( nearCol ) indexType = 6;
743 else if ( rowDiff != 0 && colDiff == 0 )
745 if ( nearRow ) indexType = 8;
753 if ( nearRow && nearCol) indexType = 10;
764 if ( rowDiff == 0 && colDiff != 0 )
766 if ( nearCol ) indexType = 0;
770 else if ( rowDiff != 0 && colDiff == 0 )
772 if ( nearRow ) indexType = 2;
777 if ( nearCol && nearRow ) indexType = 4;
782 if ( rowDiff == 0 && TMath::Abs(colDiff) == 1 )
786 eCellSameRowSameTCardNearCol = eCell;
787 tCellSameRowSameTCardNearCol = tCell;
791 eCellSameRowDiffTCardNearCol = eCell;
792 tCellSameRowDiffTCardNearCol = tCell;
798 Float_t eCellDiff = eCellMax - eCell;
799 Float_t eClusDiff = energy - eCell;
802 Float_t tCellDiff = tCellMax - tCell;
804 fhTCardCorrECellMaxDiff[indexType][matched]->Fill(energy, eCellDiff, GetEventWeight());
805 fhTCardCorrEClusterDiff[indexType][matched]->Fill(energy, eClusDiff, GetEventWeight());
808 fhTCardCorrTCellMaxDiff[indexType][matched]->Fill(energy, tCellDiff, GetEventWeight());
810 if ( fStudyExotic && exoticity > 0.97 )
812 fhTCardCorrECellMaxDiffExo[indexType][matched]->Fill(energy, eCellDiff, GetEventWeight());
813 fhTCardCorrEClusterDiffExo[indexType][matched]->Fill(energy, eClusDiff, GetEventWeight());
816 fhTCardCorrTCellMaxDiffExo[indexType][matched]->Fill(energy, tCellDiff, GetEventWeight());
820 if ( fStudyExotic && exoticity > 0.97 )
823 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
824 GetModuleNumberCellIndexesAbsCaloMap(absId,GetCalorimeter(), icol2, irow2, iRCU2, icolAbs2, irowAbs2);
828 if(energy >= 5 && energy < 8)
829 fhColRowTCardCorrNoSelectionExotic2ndCellDiffLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
831 fhColRowTCardCorrNoSelectionExotic2ndCellDiffHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
835 if(energy >= 5 && energy < 8)
836 fhColRowTCardCorrNoSelectionExotic2ndCellSameLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
838 fhColRowTCardCorrNoSelectionExotic2ndCellSameHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
842 if ( eCell > emax2nd )
848 sameTCard2ndMax = kTRUE;
854 sameTCard2ndMax = kFALSE;
889 Float_t ratioNcells = nCellWithWeight/(ncells*1.);
890 fhNCellsTCardCorrRatioWithWeightNoSelection[matched]->Fill(energy, ratioNcells, GetEventWeight());
894 if ( nCellWithWeight <= 1 )
return;
899 printf(
"AliAnaClusterShapeCorrelStudies: M02 %f, M20 %f, E %2.3f, ncell %d, n with weight %d; max cell E %2.3f\n",
900 m02,m20,energy,ncells,nCellWithWeight,eCellMax);
909 if(energy >= 5 && energy < 8)
910 fhColRowTCardCorrNoSelectionLowE [matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
912 fhColRowTCardCorrNoSelectionHighE[matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
914 if ( fStudyExotic && exoticity > 0.97 )
916 if(energy >= 5 && energy < 8)
917 fhColRowTCardCorrNoSelectionExoticLowE [matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
919 fhColRowTCardCorrNoSelectionExoticHighE[matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
922 fhTimeTCardCorrNoSelection [matched]->Fill(energy,time, GetEventWeight());
923 fhLambda0TCardCorrNoSelection[matched]->Fill(energy, m02, GetEventWeight());
924 fhLambda1TCardCorrNoSelection[matched]->Fill(energy, m20, GetEventWeight());
928 fhLambda0NLM1TCardCorrNoSelection[matched]->Fill(energy, m02, GetEventWeight());
929 fhLambda1NLM1TCardCorrNoSelection[matched]->Fill(energy, m20, GetEventWeight());
933 fhLambda0NLM2TCardCorrNoSelection[matched]->Fill(energy, m02, GetEventWeight());
934 fhLambda1NLM2TCardCorrNoSelection[matched]->Fill(energy, m20, GetEventWeight());
937 fhLambdaRTCardCorrNoSelection[matched]->Fill(energy,lamR, GetEventWeight());
938 fhNLocMaxTCardCorrNoSelection[matched]->Fill(energy, nlm, GetEventWeight());
939 fhExoticTCardCorrNoSelection [matched]->Fill(energy, exoticity, GetEventWeight());
943 fhEMaxRatNLM1TCardCorrNoSelection [matched]->Fill(energy, eCellMax/energy , GetEventWeight());
944 fhE2ndRatNLM1TCardCorrNoSelection [matched]->Fill(energy, emax2nd/energy , GetEventWeight());
945 fhE2ndEMaxRatNLM1TCardCorrNoSelection[matched]->Fill(energy, emax2nd/eCellMax, GetEventWeight());
949 fhEMaxRatNLM2TCardCorrNoSelection [matched]->Fill(energy, eCellMax/energy , GetEventWeight());
950 fhE2ndRatNLM2TCardCorrNoSelection [matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
951 fhE2ndEMaxRatNLM2TCardCorrNoSelection[matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
955 fhEMaxRatNLM3TCardCorrNoSelection [matched]->Fill(energy, eCellMax/energy , GetEventWeight());
956 fhE2ndRatNLM3TCardCorrNoSelection [matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
957 fhE2ndEMaxRatNLM3TCardCorrNoSelection[matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
964 fhE2ndSameRatNLM1TCardCorrNoSelection [matched]->Fill(energy, emax2nd/energy , GetEventWeight());
965 fhE2ndSameEMaxRatNLM1TCardCorrNoSelection[matched]->Fill(energy, emax2nd/eCellMax, GetEventWeight());
969 fhE2ndSameRatNLM2TCardCorrNoSelection [matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
970 fhE2ndSameEMaxRatNLM2TCardCorrNoSelection[matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
974 fhE2ndSameRatNLM3TCardCorrNoSelection [matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
975 fhE2ndSameEMaxRatNLM3TCardCorrNoSelection[matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
979 fhNCellsTCardCorrNoSelection [matched]->Fill(energy, ncells, GetEventWeight());
980 fhNCellsTCardCorrWithWeightNoSelection[matched]->Fill(energy, nCellWithWeight, GetEventWeight());
982 if(eCellSameRowSameTCardNearCol > 0 && eCellSameRowDiffTCardNearCol > 0)
984 Float_t eDiff = eCellSameRowSameTCardNearCol - eCellSameRowDiffTCardNearCol ;
985 Float_t tDiff = tCellSameRowSameTCardNearCol - tCellSameRowDiffTCardNearCol ;
987 fhSameRowDiffColAndTCardCellsEnergyDiffClusterE[matched]->Fill(energy , eDiff, GetEventWeight());
988 fhSameRowDiffColAndTCardCellsTimeDiffClusterE [matched]->Fill(energy , tDiff, GetEventWeight());
989 fhSameRowDiffColAndTCardCellsEnergyDiffCellMaxE[matched]->Fill(eCellMax, eDiff, GetEventWeight());
990 fhSameRowDiffColAndTCardCellsTimeDiffCellMaxE [matched]->Fill(eCellMax, tDiff, GetEventWeight());
992 if ( fStudyExotic && exoticity > 0.97 )
994 fhSameRowDiffColAndTCardCellsEnergyDiffClusterEExo[matched]->Fill(energy , eDiff, GetEventWeight());
995 fhSameRowDiffColAndTCardCellsTimeDiffClusterEExo [matched]->Fill(energy , tDiff, GetEventWeight());
996 fhSameRowDiffColAndTCardCellsEnergyDiffCellMaxEExo[matched]->Fill(eCellMax, eDiff, GetEventWeight());
997 fhSameRowDiffColAndTCardCellsTimeDiffCellMaxEExo [matched]->Fill(eCellMax, tDiff, GetEventWeight());
1002 Int_t nCorrInd = nCorr;
1003 if(nCorr > 4) nCorrInd = 5;
1005 Int_t nCorrNoInd = nCorrNo;
1006 if(nCorrNoInd > 4) nCorrNoInd = 5;
1019 fhLambda0NLM1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m02, GetEventWeight());
1020 fhLambda1NLM1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m20, GetEventWeight());
1024 fhLambda0NLM2TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m02, GetEventWeight());
1025 fhLambda1NLM2TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m20, GetEventWeight());
1028 fhLambda0TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m02, GetEventWeight());
1029 fhLambda1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m20, GetEventWeight());
1031 fhNLocMaxTCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, nlm, GetEventWeight());
1034 fhExoticTCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, exoticity, GetEventWeight());
1038 fhEMaxRatNLM1TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, eCellMax/energy , GetEventWeight());
1039 fhE2ndRatNLM1TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, emax2nd/energy , GetEventWeight());
1040 fhE2ndEMaxRatNLM1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, emax2nd/eCellMax, GetEventWeight());
1044 fhEMaxRatNLM2TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, eCellMax/energy, GetEventWeight());
1045 fhE2ndRatNLM2TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
1046 fhE2ndEMaxRatNLM2TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
1050 fhEMaxRatNLM3TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, eCellMax/energy, GetEventWeight());
1051 fhE2ndRatNLM3TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
1052 fhE2ndEMaxRatNLM3TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
1056 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1058 Int_t absId = clus->GetCellsAbsId()[ipos];
1060 Float_t eCell = fCaloCellList->GetCellAmplitude(absId);
1061 Double_t tCell = fCaloCellList->GetCellTime(absId);
1063 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
1064 GetCaloUtils()->RecalibrateCellTime(tCell, GetCalorimeter(), absId, GetReader()->GetInputEvent()->GetBunchCrossNumber());
1066 if(tCell>400) tCell-=fConstantTimeShift;
1069 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
1071 if( absId == absIdMax || weight < 0.01 )
continue;
1073 Float_t tDiffMaxSecondary = tCellMax - tCell;
1075 fhTimeDiffTCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, tDiffMaxSecondary, GetEventWeight());
1076 if ( fStudyExotic && exoticity > 0.97 )
1078 fhTimeDiffExoTCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, tDiffMaxSecondary, GetEventWeight());
1080 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
1081 GetModuleNumberCellIndexesAbsCaloMap(absId,GetCalorimeter(), icol2, irow2, iRCU2, icolAbs2, irowAbs2);
1083 Int_t rowDiff = -100, colDiff = -100;
1084 Bool_t sameTCard = GetCaloUtils()->IsAbsIDsFromTCard(absIdMax,absId,rowDiff,colDiff);
1088 if ( energy >= 5 && energy < 8 && nCorr == 0 )
1089 fhColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
1090 else if ( energy >= 8 && nCorr == 0 )
1091 fhColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
1095 if ( energy >= 5 && energy < 8 && nCorrNo == 0 )
1096 fhColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
1097 else if ( energy >= 8 && nCorrNo == 0 )
1098 fhColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
1104 fhECellClusRatNLM1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, eCell/energy, GetEventWeight());
1105 fhLogECellNLM1TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
1106 fhECellWeightNLM1TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, weight, GetEventWeight());
1110 fhECellClusRatNLM2TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, eCell/energy, GetEventWeight());
1111 fhLogECellNLM2TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
1112 fhECellWeightNLM2TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, weight, GetEventWeight());
1116 fhECellClusRatNLM3TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, eCell/energy, GetEventWeight());
1117 fhLogECellNLM3TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
1118 fhECellWeightNLM3TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, weight, GetEventWeight());
1123 if(m02 > fInvMassMinM02Cut && m02 < fInvMassMaxM02Cut)
1125 for(
Int_t jclus = 0 ; jclus < fCaloClusList->GetEntriesFast() ; jclus++)
1127 AliVCluster* clus2 = (AliVCluster*) fCaloClusList->At(jclus);
1130 Int_t absIdMax2 = GetCaloUtils()->GetMaxEnergyCell(fCaloCellList, clus2, maxCellFraction);
1132 Double_t tof2 = clus2->GetTOF()*1.e9;
1133 if(tof2>400) tof2-=fConstantTimeShift;
1140 if( absIdMax == absIdMax2
1141 || !IsGoodCluster(absIdMax2, clus2->GetM02(), clus2->GetNCells())
1142 || GetCaloUtils()->GetNumberOfLocalMaxima(clus2,fCaloCellList) > 1
1143 || clus2->GetM02() > fInvMassMaxM02Cut
1144 || clus2->GetM02() < fInvMassMinM02Cut
1145 || clus2->E() < fInvMassMinECut
1146 || clus2->E() > fInvMassMaxECut
1147 || TMath::Abs(diffTof) > fInvMassMaxTimeDifference
1152 clus2->GetMomentum(fClusterMomentum2,v);
1156 if(IsFiducialCutOn()) in2 = GetFiducialCut()->IsInFiducialCut(fClusterMomentum2.Eta(),fClusterMomentum2.Phi(),GetCalorimeter()) ;
1162 Double_t angle = fClusterMomentum.Angle(fClusterMomentum2.Vect());
1164 if( angle > fInvMassMaxOpenAngle )
continue;
1167 Float_t mass = (fClusterMomentum+fClusterMomentum2).M ();
1168 fhMassEClusTCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, mass, GetEventWeight());
1173 if ( energy >= 5 && energy < 8)
1174 fhColRowTCardCorrelNCellLowE [nCorrInd][nCorrNoInd][matched]->Fill(icolAbs, irowAbs, GetEventWeight());
1175 else if ( energy >= 8 )
1176 fhColRowTCardCorrelNCellHighE[nCorrInd][nCorrNoInd][matched]->Fill(icolAbs, irowAbs, GetEventWeight());
1191 if ( fStudyExotic && exoticity > 0.97 )
1193 if ( energy >= 5 && energy < 8)
1194 fhColRowTCardCorrelNCellExoticLowE [nCorrInd][nCorrNoInd][matched]->Fill(icolAbs, irowAbs, GetEventWeight());
1195 else if ( energy >= 8 )
1196 fhColRowTCardCorrelNCellExoticHighE[nCorrInd][nCorrNoInd][matched]->Fill(icolAbs, irowAbs, GetEventWeight());
1233 fhLambda0Lambda1 [ebin][matched]->Fill(m20, m02, GetEventWeight());
1234 fhNCellsTCardSameAndDiff[ebin][matched]->Fill(nCorrNo, nCorr, GetEventWeight());
1240 fhLambda0Exoticity[ebin][matched]->Fill(exoticity, m02, GetEventWeight());
1241 fhLambda1Exoticity[ebin][matched]->Fill(exoticity, m20, GetEventWeight());
1243 fhNCellsExoticity [ebin][matched]->Fill(exoticity, nCellWithWeight, GetEventWeight());
1244 fhTimeExoticity [ebin][matched]->Fill(exoticity, tCellMax, GetEventWeight());
1248 fhLambda0ExoticityPerNCell[nCorrInd][nCorrNoInd][matched]->Fill(exoticity, m02, GetEventWeight());
1249 fhLambda1ExoticityPerNCell[nCorrInd][nCorrNoInd][matched]->Fill(exoticity, m20, GetEventWeight());
1261 if ( exoticity > 0.97 )
1262 fhNCellsTCardSameAndDiffExotic[ebin][matched]->Fill(nCorrNo, nCorr, GetEventWeight());
1265 fhTMPhiResidualExoticity[ebin]->Fill(exoticity,dphi);
1266 fhTMEtaResidualExoticity[ebin]->Fill(exoticity,deta);
1276 Float_t nCellRat = nCorr*1. / ((nCorr+nCorrNo)*1.);
1277 fhNCellsTCardSameAndDiffFraction[matched]->Fill(energy, nCellRat, GetEventWeight());
1278 if ( fStudyExotic && exoticity > 0.97 )
1279 fhNCellsTCardSameAndDiffFractionExotic[matched]->Fill(energy, nCellRat, GetEventWeight());
1378 Int_t nCaloCellsPerCluster = clus->GetNCells();
1379 if ( nCaloCellsPerCluster < 2 )
return;
1384 Int_t nlm = GetCaloUtils()->GetNumberOfLocalMaxima(clus, fCaloCellList) ;
1396 for (
Int_t ipos = 0; ipos < nCaloCellsPerCluster; ipos++)
1398 Int_t absId = clus ->GetCellsAbsId()[ipos];
1399 Float_t eCell = fCaloCellList->GetCellAmplitude(absId) ;
1401 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
1403 if( absId == absIdMax || eCell < 0.01 )
continue;
1405 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
1407 if( weight < 0.01 )
continue;
1412 if ( nCell < 1 )
return;
1417 Int_t ietaMax=-1, iphiMax = 0, rcuMax = 0, icolAbs = -1, irowAbs = -1;
1418 Int_t smMax = GetModuleNumberCellIndexesAbsCaloMap(absIdMax, GetCalorimeter(),
1419 ietaMax, iphiMax, rcuMax, icolAbs, irowAbs);
1421 if ( matchedPID == 0 && energy > fEMinShape && energy < fEMaxShape )
1423 fhColRowM02->Fill(icolAbs,irowAbs,m02,GetEventWeight()) ;
1425 if ( nCell > fNCellMinShape )
1426 fhColRowM02NCellCut->Fill(icolAbs,irowAbs,m02,GetEventWeight()) ;
1434 Int_t etaRegion = -1, phiRegion = -1;
1435 GetCaloUtils()->GetEMCALSubregion(clus,fCaloCellList,etaRegion,phiRegion);
1437 if ( etaRegion !=0 ) return ;
1441 for (
Int_t ipos = 0; ipos < nCaloCellsPerCluster; ipos++)
1443 Int_t absId = clus ->GetCellsAbsId()[ipos];
1444 Float_t eCell = fCaloCellList->GetCellAmplitude(absId) ;
1446 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
1448 if( absId == absIdMax || eCell < 0.01 )
continue;
1450 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
1452 if( weight < 0.01 )
continue;
1454 Float_t fracCell = eCell/eCellMax;
1455 fhClusterMaxCellCloseCellRatioM02->Fill(energy, fracCell, m02, GetEventWeight());
1457 Float_t fracClus = (energy-eCell)/energy;
1459 Float_t logECell = TMath::Log(eCell);
1470 if ( energy > fEMinShape && energy < fEMaxShape && matchedPID == 0 )
1472 fhESecCellEMaxCellM02NCellPerSM[smMax]->Fill(fracCell, nCell, m02, GetEventWeight());
1473 fhESecCellEClusterM02NCellPerSM[smMax]->Fill(fracClus, nCell, m02, GetEventWeight());
1474 fhESecCellLogM02NCellPerSM [smMax]->Fill(logECell, nCell, m02, GetEventWeight());
1476 if ( m02 > 0.1 && m02 < 0.3 )
1478 fhEMaxESecCellNCellLowM02PerSM [smMax]->Fill(eCellMax, eCell , nCell, GetEventWeight());
1479 fhEMaxECrossNCellLowM02PerSM [smMax]->Fill(eCellMax, eCrossFrac, nCell, GetEventWeight());
1481 Double_t time = fCaloCellList->GetCellTime(absId);
1482 GetCaloUtils()->RecalibrateCellTime(time, GetCalorimeter(), absId,GetReader()->GetInputEvent()->GetBunchCrossNumber());
1484 if(time>400) time-=fConstantTimeShift;
1486 fhESecCellTimeNCellSM->Fill(time, smMax, nCell, GetEventWeight());
1488 else if ( m02 > 0.5 && m02 < 2 )
1490 fhEMaxESecCellNCellHighM02PerSM[smMax]->Fill(eCellMax, eCell , nCell, GetEventWeight());
1491 fhEMaxECrossNCellHighM02PerSM [smMax]->Fill(eCellMax, eCrossFrac, nCell, GetEventWeight());
1501 Int_t sm = GetModuleNumberCellIndexes(absId,GetCalorimeter(), ieta, iphi, rcu);
1503 if(dIphi < TMath::Abs(iphi-iphiMax)) dIphi = TMath::Abs(iphi-iphiMax);
1504 if(iphi-iphiMax < 0 && dIphiNeg > iphi-iphiMax) dIphiNeg = iphi-iphiMax;
1505 if(iphi-iphiMax > 0 && dIphiPos < iphi-iphiMax) dIphiPos = iphi-iphiMax;
1509 if(dIeta < TMath::Abs(ieta-ietaMax)) dIeta = TMath::Abs(ieta-ietaMax);
1510 if(ieta-ietaMax < 0 && dIetaNeg > ieta-ietaMax) dIetaNeg = ieta-ietaMax;
1511 if(ieta-ietaMax > 0 && dIetaPos < ieta-ietaMax) dIetaPos = ieta-ietaMax;
1513 if ( energy > fEMinShape && energy < fEMaxShape && matchedPID == 0 )
1516 if ( nCell == 2 || nCell == 3 ) nCellBin = 0;
1517 else if ( nCell == 4 || nCell == 5 ) nCellBin = 1;
1518 if ( m02 > 0.1 && m02 < 0.3 )
1520 fhColRowFromCellMaxLowM02PerSM[smMax][ietaMax%2]->Fill(ietaMax-ieta, iphiMax-iphi, nCell, GetEventWeight());
1521 fhColRowFromCellMaxLowM02 [ietaMax%2]->Fill(ietaMax-ieta, iphiMax-iphi, GetEventWeight());
1523 fhColRowFromCellMaxEMaxSecDiffFracLowM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi,(eCellMax-eCell)/eCellMax, GetEventWeight());
1524 fhColRowFromCellMaxEMaxSecDiffLowM02PerSM [smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi, eCellMax-eCell, GetEventWeight());
1525 fhColRowFromCellMaxECellClusterRatLowM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi, eCell/energy , GetEventWeight());
1527 else if (m02 > 0.5 && m02 < 2)
1529 fhColRowFromCellMaxHighM02PerSM[smMax][ietaMax%2]->Fill(ietaMax-ieta, iphiMax-iphi, nCell, GetEventWeight());
1530 fhColRowFromCellMaxHighM02 [ietaMax%2]->Fill(ietaMax-ieta, iphiMax-iphi, GetEventWeight());
1532 fhColRowFromCellMaxEMaxSecDiffFracHighM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi,(eCellMax-eCell)/eCellMax, GetEventWeight());
1533 fhColRowFromCellMaxEMaxSecDiffHighM02PerSM [smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi, eCellMax-eCell, GetEventWeight());
1534 fhColRowFromCellMaxECellClusterRatHighM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi, eCell/energy , GetEventWeight());
1541 Int_t ietaShift = ieta;
1542 Int_t ietaMaxShift = ietaMax;
1544 if (ieta > ietaMax) ietaMaxShift+=48;
1545 else ietaShift +=48;
1547 if(dIeta < TMath::Abs(ietaShift-ietaMaxShift)) dIeta = TMath::Abs(ietaShift-ietaMaxShift);
1548 if(ietaShift-ietaMaxShift < 0 && dIetaNeg > ietaShift-ietaMaxShift) dIetaNeg = ietaShift-ietaMaxShift;
1549 if(ietaShift-ietaMaxShift > 0 && dIetaPos < ietaShift-ietaMaxShift) dIetaPos = ietaShift-ietaMaxShift;
1555 if ( energy > fEMinShape && energy < fEMaxShape && matchedPID == 0 )
1557 fhNCellsPerClusterM02M20PerSM[smMax]->Fill(m20, nCell, m02, GetEventWeight());
1558 fhNCellsPerClusterM02NLMPerSM[smMax]->Fill(nlm, nCell, m02, GetEventWeight());
1562 if ( nCell == 2 || nCell == 3 ) nCellBin = 0;
1563 else if ( nCell == 4 || nCell == 5 ) nCellBin = 1;
1564 if ( m02 > 0.1 && m02 < 0.3 )
1565 fhColRowFromCellMaxECellClusterRatLowM02PerSM [smMax][ietaMax%2][nCellBin]->Fill(0., 0., eCellMax/energy, GetEventWeight());
1566 else if (m02 > 0.5 && m02 < 2)
1567 fhColRowFromCellMaxECellClusterRatHighM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(0., 0., eCellMax/energy, GetEventWeight());
1571 fhClusterMaxCellECrossM02->Fill(energy, eCrossFrac, m02, GetEventWeight());
1576 fhClusterMaxCellDiffM02[matchedPID]->Fill(energy, maxFrac, m02, GetEventWeight());
1577 fhClusterTimeEnergyM02 [matchedPID]->Fill(energy, tmax , m02, GetEventWeight());
1578 fhNCellsPerClusterM02 [matchedPID]->Fill(energy, nCell , m02, GetEventWeight());
1579 fhNCellsPerClusterM20 [matchedPID]->Fill(energy, nCell , m20, GetEventWeight());
1580 fhNCellsPerClusterNLM [matchedPID]->Fill(energy, nCell , nlm, GetEventWeight());
1582 if ( matchedPID == 0 )
1584 fhNCellsPerClusterM02PerSMPerCol[smMax][ietaMax%2]->Fill(energy, nCell, m02, GetEventWeight());
1585 fhNCellsPerClusterM20PerSMPerCol[smMax][ietaMax%2]->Fill(energy, nCell, m20, GetEventWeight());
1586 fhNCellsPerClusterNLMPerSMPerCol[smMax][ietaMax%2]->Fill(energy, nCell, nlm, GetEventWeight());
1589 fhSMNCell [matchedPID]->Fill(energy, smMax , nCell, GetEventWeight());
1590 fhSMM02NoCut [matchedPID]->Fill(energy, smMax , m02 , GetEventWeight());
1591 if ( m02 > 0.1 && m02 < 0.3 )
1592 fhSMM20LowM02NoCut [matchedPID]->Fill(energy, smMax , m20 , GetEventWeight());
1593 else if ( m02 > 0.5 && m02 < 2 )
1594 fhSMM20HighM02NoCut[matchedPID]->Fill(energy, smMax , m20 , GetEventWeight());
1597 if ( energy > fEMinShape && energy < fEMaxShape )
1599 fhSMNCellM02[matchedPID]->Fill(smMax , nCell, m02, GetEventWeight());
1601 if ( matchedPID == 0 )
1603 if ( m02 > 0.1 && m02 < 0.3 )
1604 fhEMaxCellTimeNCellSM->Fill(tmax, smMax, nCell, GetEventWeight());
1606 fhEMaxCellTimeM02SM->Fill(tmax, smMax, m02, GetEventWeight());
1608 fhEMaxCellEClusterM02NCellPerSM[smMax]->Fill(maxFrac, nCell, m02, GetEventWeight());
1609 fhEMaxCellLogM02NCellPerSM [smMax]->Fill(TMath::Log(eCellMax), nCell, m02, GetEventWeight());
1613 if ( nCell > fNCellMinShape )
1615 if ( m02 > 0.1 && m02 < 0.3 )
1616 fhSMM20LowM02 [matchedPID]->Fill(energy, smMax , m20 , GetEventWeight());
1617 else if ( m02 > 0.5 && m02 < 2 )
1618 fhSMM20HighM02[matchedPID]->Fill(energy, smMax , m20 , GetEventWeight());
1620 fhSMM02 [matchedPID]->Fill(energy, smMax , m02, GetEventWeight());
1621 fhColM02[matchedPID]->Fill(energy, ietaMax, m02, GetEventWeight());
1622 fhRowM02[matchedPID]->Fill(energy, iphiMax, m02, GetEventWeight());
1626 if ( fStudyShapeParam )
1629 Float_t dIA = 1.*(dIphi-dIeta)/(dIeta+dIphi);
1630 Float_t dIATot = 1.*((dIphiPos-dIphiNeg)-(dIetaPos-dIetaNeg))/((dIetaPos-dIetaNeg)+(dIphiPos-dIphiNeg));
1639 AliDebug(1,Form(
"E %2.2f, nCell %d, dPhi %d, dEta %d, dIA %2.2f, match %d",energy,nCell, dIphi,dIeta,dIA,matchedPID));
1649 fhDeltaIANCells [matchedPID]->Fill(energy, nCell , dIA , GetEventWeight());
1650 fhDeltaIATotNCells [matchedPID]->Fill(energy, nCell , dIATot, GetEventWeight());
1652 if ( nCell > fNCellMinShape )
1654 fhDeltaIEtaDeltaIPhi[matchedPID]->Fill(energy, dIeta, dIphi, GetEventWeight());
1655 fhDeltaIA [matchedPID]->Fill(energy, dIA , GetEventWeight());
1656 fhDeltaIAM02 [matchedPID]->Fill(energy, m02 , dIA , GetEventWeight());
1657 fhDeltaIAM20 [matchedPID]->Fill(energy, m20 , dIA , GetEventWeight());
1659 fhDeltaIEtaDeltaIPhiTot[matchedPID]->Fill(energy,dIetaPos-dIetaNeg, dIphiPos-dIphiNeg, GetEventWeight());
1660 fhDeltaIATot [matchedPID]->Fill(energy, dIATot , GetEventWeight());
1661 fhDeltaIATotM02 [matchedPID]->Fill(energy, m02 , dIATot, GetEventWeight());
1662 fhDeltaIATotM20 [matchedPID]->Fill(energy, m20 , dIATot, GetEventWeight());
1665 if ( IsDataMC() && mcIndex > -1 && mcIndex < 10 && nCell > fNCellMinShape )
1667 fhDeltaIAOrigin [matchedPID]->Fill(energy, mcIndex, dIA , GetEventWeight());
1668 fhDeltaIATotOrigin[matchedPID]->Fill(energy, mcIndex, dIATot, GetEventWeight());
1672 Float_t dispp= 0., dEta = 0., dPhi = 0.;
1673 Float_t sEta = 0., sPhi = 0., sEtaPhi = 0.;
1674 if ( GetCalorimeter() ==
kEMCAL )
1676 GetCaloUtils()->GetEMCALRecoUtils()->RecalculateClusterShowerShapeParameters(GetEMCALGeometry(), fCaloCellList, clus,
1677 l0, l1, dispp, dEta, dPhi, sEta, sPhi, sEtaPhi);
1680 if(sEta+sPhi>0.0001) sEtaPhiA = (sPhi-sEta)/(sEta+sPhi);
1682 AliDebug(2,Form(
"Recalculate shower shape org: m02 %2.2f, m20 %2.2f, disp %2.2f;"
1683 " new: m02 %2.2f, m20 %2.2f, disp %2.2f; "
1684 "mEta %2.2f, mPhi %2.2f, mEtaPhi %2.2f, A_EtaPhi %2.2f; dEta %2.2f dPhi %2.2f",
1685 m02,m20,clus->GetDispersion(),l0,l1,dispp,sEta,sPhi,sEtaPhi,sEtaPhiA,dEta,dPhi));
1687 fhNCellsPerClusterMEta [matchedPID]->Fill(energy, nCell, sEta , GetEventWeight());
1688 fhNCellsPerClusterMPhi [matchedPID]->Fill(energy, nCell, sPhi , GetEventWeight());
1689 fhNCellsPerClusterMEtaPhi [matchedPID]->Fill(energy, nCell, sEtaPhi , GetEventWeight());
1690 fhNCellsPerClusterMEtaPhiA[matchedPID]->Fill(energy, nCell, sEtaPhiA, GetEventWeight());
1695 if ( matchedPID == 0 && energy > fEMinShape && energy < fEMaxShape &&
1696 m02 > fInvMassMinM02Cut && m02 < fInvMassMaxM02Cut )
1698 for(
Int_t jclus = 0 ; jclus < fCaloClusList->GetEntriesFast() ; jclus++)
1700 AliVCluster* clus2 = (AliVCluster*) fCaloClusList->At(jclus);
1703 Int_t absIdMax2 = GetCaloUtils()->GetMaxEnergyCell(fCaloCellList, clus2, maxCellFraction);
1705 Double_t tof2 = clus2->GetTOF()*1.e9;
1706 if(tof2>400) tof2-=fConstantTimeShift;
1713 if( absIdMax == absIdMax2
1714 || !IsGoodCluster(absIdMax2, clus2->GetM02(), clus2->GetNCells())
1715 || GetCaloUtils()->GetNumberOfLocalMaxima(clus2,fCaloCellList) > 1
1716 || clus2->GetM02() > fInvMassMaxM02Cut
1717 || clus2->GetM02() < fInvMassMinM02Cut
1718 || clus2->E() < fInvMassMinECut
1719 || clus2->E() > fInvMassMaxECut
1720 || TMath::Abs(diffTof) > fInvMassMaxTimeDifference
1725 clus2->GetMomentum(fClusterMomentum2,v);
1729 if(IsFiducialCutOn()) in2 = GetFiducialCut()->IsInFiducialCut(fClusterMomentum2.Eta(),fClusterMomentum2.Phi(),GetCalorimeter()) ;
1735 Double_t angle = fClusterMomentum.Angle(fClusterMomentum2.Vect());
1737 if( angle > fInvMassMaxOpenAngle )
continue;
1740 Float_t mass = (fClusterMomentum+fClusterMomentum2).M ();
1741 fhInvMassNCellSM->Fill(mass, nCell, smMax, GetEventWeight());
1743 Int_t smMax2 = GetModuleNumber(clus2);
1744 if(smMax == smMax2) fhInvMassNCellSMSame->Fill(mass, nCell, smMax, GetEventWeight());
1749 if ( IsDataMC() && mcIndex > -1 && mcIndex < 10)
1751 fhOriginE [matchedPID]->Fill(energy, mcIndex, GetEventWeight());
1752 fhOriginM02[matchedPID]->Fill(energy, mcIndex, m02, GetEventWeight());
1768 AliVTrack *track = GetCaloUtils()->GetMatchedTrack(clus, GetReader()->GetInputEvent());
1778 Double_t dedx = track->GetTPCsignal();
1793 if(fdEdXMinEle == 0 || fdEdXMaxEle == 0 || fdEdXMinHad == 0 || fdEdXMaxHad == 0)
1794 InitdEdXParameters();
1796 if ( dedx >= fdEdXMinEle && dedx < fdEdXMaxEle ) matchedPID = 1;
1797 else if ( dedx >= fdEdXMinHad && dedx < fdEdXMaxHad ) matchedPID = 2;
1800 AliDebug(1,Form(
"dEdX out of range %2.2f",dedx));
1818 Int_t nCaloCellsPerCluster = 0 ;
1829 for(
Int_t iclus = 0; iclus < nCaloClusters; iclus++)
1831 AliDebug(1,Form(
"Cluster: %d/%d, data %d",iclus+1,nCaloClusters,
GetReader()->GetDataType()));
1833 AliVCluster* clus = (AliVCluster*)
fCaloClusList->At(iclus);
1836 if ( clus->GetDistanceToBadChannel() <
fMinDistToBad )
return ;
1840 if ( nModule < fFirstModule || nModule >
fLastModule )
1842 AliDebug(1,Form(
"Cluster module out of range %d",nModule));
1852 Double_t tof = clus->GetTOF()*1.e9;
1870 AliDebug(1,Form(
"Remove cluster with phi %2.2f and eta %2.2f", phi*TMath::RadToDeg(), eta));
1874 AliDebug(1,Form(
"cluster: E %2.3f, pT %2.3f, eta %2.3f, phi %2.3f",e,pt,eta,phi*TMath::RadToDeg()));
1878 nCaloCellsPerCluster = clus->GetNCells();
1882 AliDebug(1,Form(
"Accept cluster? %d",goodCluster));
1884 if(!goodCluster)
continue;
1913 Int_t matchedPID = 0;
1930 if ( ampMax > 0.01 )
1938 if (
fStudyShape && matchedPID >= 0 && matchedPID < 3 )
1959 const Int_t buffersize = 255;
1960 char onePar[buffersize] ;
1962 snprintf(onePar,buffersize,
"--- AliAnaClusterShapeCorrelStudies ---:") ;
1983 return new TObjString(parList) ;
1993 outputContainer->SetName(
"ClusterShapeStudies") ;
2004 TString mcParticleStringLabel[] = {
"Merged #gamma#gamma",
"Decay #gamma",
"Direct #gamma",
"e^{#pm}",
"#pi^{#pm}",
"k^{#pm}",
"p",
"#bar{p}",
"n",
"#bar{n}"};
2008 TString matchCase[] = {
"Neutral",
"Electron",
"Hadron"};
2036 Float_t rowcellmin = fNMaxRowsFullMin-1.5;
2045 Int_t nShShBins = 200;
2055 Int_t cellBins = 15;
2065 (
"hColRowHighEPosTime",
2066 "column vs row, exo < 0.97, E > 8 GeV, t > 5 ns",
2067 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2073 (
"hColRowHighENegTime",
2074 "column vs row, exo < 0.97, E > 8 GeV, t < -5 ns",
2075 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2081 (
"hColRowHighENulTime",
2082 "column vs row, exo < 0.97, E > 8 GeV, -5 < t < 5 ns",
2083 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2091 (
"hColRowExoticHighE1CellPosTime",
2092 "column vs row, 1 cell, E > 8 GeV, t > 5 ns",
2093 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2099 (
"hColRowExoticHighEPosTime",
2100 "column vs row, exo > 0.97, E > 8 GeV, t > 5 ns",
2101 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2107 (
"hColRowExoticHighE1CellNegTime",
2108 "column vs row, 1 cell, E > 8 GeV, t < -5 ns",
2109 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2115 (
"hColRowExoticHighENegTime",
2116 "column vs row, exo > 0.97, E > 8 GeV, t < -5 ns",
2117 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2123 (
"hColRowExoticHighE1CellNulTime",
2124 "column vs row, 1 cell, E > 8 GeV, -5 < t < 5 ns",
2125 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2131 (
"hColRowExoticHighENulTime",
2132 "column vs row, exo > 0.97, E > 8 GeV, -5 < t < 5 ns",
2133 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2139 TString add[] = {
"",
"TrackMatched"};
2140 for(
Int_t tm = 0; tm < 2; tm++)
2143 (Form(
"hColRowTCardCorrNoSelectionLowE%s",add[tm].
Data()),
2144 Form(
"column vs row, max E cell for TCard correlation selected clusters, 5 < E < 8 GeV %s",add[tm].
Data()),
2145 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2151 (Form(
"hColRowTCardCorrNoSelectionHighE%s",add[tm].
Data()),
2152 Form(
"column vs row, max E cell for TCard correlation selected clusters, E > 8 GeV %s",add[tm].
Data()),
2153 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2161 (Form(
"hNCellsTCardCorrNoSelection%s",add[tm].
Data()),
2162 Form(
"# custer # cells vs #it{E} %s",add[tm].
Data()),
2163 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
2169 (Form(
"hNCellsTCardCorrWithWeightNoSelection%s",add[tm].
Data()),
2170 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
2171 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
2177 (Form(
"hNCellsTCardCorrRatioWithWeightNoSelection%s",add[tm].
Data()),
2178 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
2179 nptbins,ptmin,ptmax, 100,0,1);
2185 (Form(
"hTimeTCardCorrNoSelection%s",add[tm].
Data()),
2186 Form(
"#it{time} vs #it{E} %s",add[tm].
Data()),
2187 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2193 (Form(
"hLambda0TCardCorrNoSelection%s",add[tm].
Data()),
2194 Form(
"#lambda^{2}_{0} vs #it{E} %s",add[tm].
Data()),
2195 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2201 (Form(
"hLambda1TCardCorrNoSelection%s",add[tm].
Data()),
2202 Form(
"#lambda^{2}_{1} vs #it{E} %s",add[tm].
Data()),
2203 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2209 (Form(
"hLambda0NLM1TCardCorrNoSelection%s",add[tm].
Data()),
2210 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=1 %s",add[tm].
Data()),
2211 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2217 (Form(
"hLambda1NLM1TCardCorrNoSelection%s",add[tm].
Data()),
2218 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=1 %s",add[tm].
Data()),
2219 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2225 (Form(
"hLambda0NLM2TCardCorrNoSelection%s",add[tm].
Data()),
2226 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=2 %s",add[tm].
Data()),
2227 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2233 (Form(
"hLambda1NLM2TCardCorrNoSelection%s",add[tm].
Data()),
2234 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=2 %s",add[tm].
Data()),
2235 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2242 (Form(
"hLambdaRTCardCorrNoSelection%s",add[tm].
Data()),
2243 Form(
"#lambda^{1}_{0}/#lambda^{2}_{0} vs #it{E} %s",add[tm].
Data()),
2244 nptbins,ptmin,ptmax,110,0,1.1);
2250 (Form(
"hNLocMaxTCardCorrNoSelection%s",add[tm].
Data()),
2251 Form(
"#it{n}_{LM} vs E %s",add[tm].
Data()),
2252 nptbins,ptmin,ptmax,10,0,10);
2258 (Form(
"hEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2259 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
2260 nptbins,ptmin,ptmax,100,0,1);
2266 (Form(
"hEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2267 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2268 nptbins,ptmin,ptmax,100,0,1);
2274 (Form(
"hEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2275 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2276 nptbins,ptmin,ptmax,100,0,1);
2282 (Form(
"hE2ndRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2283 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
2284 nptbins,ptmin,ptmax,100,0,1);
2290 (Form(
"hE2ndRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2291 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2292 nptbins,ptmin,ptmax,100,0,1);
2298 (Form(
"hE2ndRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2299 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2300 nptbins,ptmin,ptmax,100,0,1);
2307 (Form(
"hE2ndEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2308 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
2309 nptbins,ptmin,ptmax,100,0,1);
2315 (Form(
"hE2ndEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2316 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2317 nptbins,ptmin,ptmax,100,0,1);
2323 (Form(
"hE2ndEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2324 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2325 nptbins,ptmin,ptmax,100,0,1);
2333 (Form(
"hE2ndSameRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2334 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()),
2335 nptbins,ptmin,ptmax,100,0,1);
2341 (Form(
"hE2ndSameRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2342 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()),
2343 nptbins,ptmin,ptmax,100,0,1);
2349 (Form(
"hE2ndSameRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2350 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()),
2351 nptbins,ptmin,ptmax,100,0,1);
2358 (Form(
"hE2ndSameEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2359 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()),
2360 nptbins,ptmin,ptmax,100,0,1);
2366 (Form(
"hE2ndSameEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2367 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()),
2368 nptbins,ptmin,ptmax,100,0,1);
2374 (Form(
"hE2ndSameEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2375 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()),
2376 nptbins,ptmin,ptmax,100,0,1);
2384 (Form(
"hECellClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2385 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
2386 nptbins,ptmin,ptmax,100,0,1);
2392 (Form(
"hECellClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2393 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2394 nptbins,ptmin,ptmax,100,0,1);
2400 (Form(
"hECellClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2401 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2402 nptbins,ptmin,ptmax,100,0,1);
2408 (Form(
"hLogECellNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2409 Form(
"log(#it{E}_{cell}) vs E_{cluster}, #it{n}_{LM}=1, w > 0.01 %s",add[tm].
Data()),
2410 nptbins,ptmin,ptmax,150,-3,3);
2416 (Form(
"hLogECellNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2417 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}=2, w > 0.01 %s",add[tm].
Data()),
2418 nptbins,ptmin,ptmax,150,-3,3);
2424 (Form(
"hLogECellNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2425 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}>2, w > 0.01 %s",add[tm].
Data()),
2426 nptbins,ptmin,ptmax,150,-3,3);
2433 (Form(
"hECellWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2434 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
2435 nptbins,ptmin,ptmax,90,0,4.5);
2441 (Form(
"hECellWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2442 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2443 nptbins,ptmin,ptmax,90,0,4.5);
2449 (Form(
"hECellWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2450 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster}))vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2451 nptbins,ptmin,ptmax,90,0,4.5);
2458 (Form(
"hECellSameClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2459 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()),
2460 nptbins,ptmin,ptmax,100,0,1);
2466 (Form(
"hECellSameClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2467 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, cell from same T-Card as leading, #it{n}_{LM}=2 %s",add[tm].
Data()),
2468 nptbins,ptmin,ptmax,100,0,1);
2474 (Form(
"hECellSameClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2475 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, cell from same T-Card as leading, #it{n}_{LM}>2 %s",add[tm].
Data()),
2476 nptbins,ptmin,ptmax,100,0,1);
2482 (Form(
"hLogECellSameNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2483 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()),
2484 nptbins,ptmin,ptmax,150,-3,3);
2490 (Form(
"hLogECellSameNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2491 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()),
2492 nptbins,ptmin,ptmax,150,-3,3);
2498 (Form(
"hLogECellSameNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2499 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()),
2500 nptbins,ptmin,ptmax,150,-3,3);
2507 (Form(
"hECellSameWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2508 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()),
2509 nptbins,ptmin,ptmax,90,0,4.5);
2515 (Form(
"hECellSameWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2516 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()),
2517 nptbins,ptmin,ptmax,90,0,4.5);
2523 (Form(
"hECellSameWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2524 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()),
2525 nptbins,ptmin,ptmax,90,0,4.5);
2531 (Form(
"hExoticTCardCorrNoSelection%s",add[tm].
Data()),
2532 Form(
"exoticity vs #it{E} %s",add[tm].
Data()),
2533 nptbins,ptmin,ptmax,200,-1,1);
2539 (Form(
"hNCellsTCardSameAndDiffFraction%s",add[tm].
Data()),
2540 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, %s",add[tm].
Data()),
2541 nptbins,ptmin,ptmax,100,0,1);
2547 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterE%s",add[tm].
Data()),
2548 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()),
2549 nptbins,ptmin,ptmax,200,-10,10);
2555 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterE%s",add[tm].
Data()),
2556 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()),
2557 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2563 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxE%s",add[tm].
Data()),
2564 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()),
2565 nptbins,ptmin,ptmax,200,-10,10);
2571 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxE%s",add[tm].
Data()),
2572 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()),
2573 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2581 (Form(
"hEnergyTime1Cell%s",add[tm].
Data()),
2582 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
2583 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2589 (Form(
"hEnergyTimeExotic%s",add[tm].
Data()),
2590 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
2591 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2597 (Form(
"hEnergyTimeTCardCorrNoSelection1Cell%s",add[tm].
Data()),
2598 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
2599 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2605 (Form(
"hEnergyTimeTCardCorrNoSelectionExotic%s",add[tm].
Data()),
2606 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
2607 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2613 (Form(
"hColRowExoticLowE1Cell%s",add[tm].
Data()),
2614 Form(
"column vs row, 1 cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2615 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2621 (Form(
"hColRowExoticHighE1Cell%s",add[tm].
Data()),
2622 Form(
"column vs row, 1 cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2623 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2629 (Form(
"hColRowExoticLowE%s",add[tm].
Data()),
2630 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2631 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2637 (Form(
"hColRowExoticHighE%s",add[tm].
Data()),
2638 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2639 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2645 (Form(
"hColRowExotic2ndCellDiffLowE%s",add[tm].
Data()),
2646 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2647 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2653 (Form(
"hColRowExotic2ndCellDiffHighE%s",add[tm].
Data()),
2654 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2655 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2661 (Form(
"hColRowExotic2ndCellSameLowE%s",add[tm].
Data()),
2662 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2663 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2669 (Form(
"hColRowExotic2ndCellSameHighE%s",add[tm].
Data()),
2670 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2671 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2677 (Form(
"hColRowTCardCorrNoSelectionExoticLowE%s",add[tm].
Data()),
2678 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2679 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2685 (Form(
"hColRowTCardCorrNoSelectionExoticHighE%s",add[tm].
Data()),
2686 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2687 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2693 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffLowE%s",add[tm].
Data()),
2694 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2695 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2701 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffHighE%s",add[tm].
Data()),
2702 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2703 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2709 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameLowE%s",add[tm].
Data()),
2710 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2711 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2717 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameHighE%s",add[tm].
Data()),
2718 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2719 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2725 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameLowE%s",add[tm].
Data()),
2726 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2727 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2733 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameHighE%s",add[tm].
Data()),
2734 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2735 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2741 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffLowE%s",add[tm].
Data()),
2742 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2743 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2749 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffHighE%s",add[tm].
Data()),
2750 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2751 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2757 (Form(
"hNCellsTCardSameAndDiffFraction_Exotic%s",add[tm].
Data()),
2758 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, exo > 0.97 %s",add[tm].
Data()),
2759 nptbins,ptmin,ptmax,100,0,1);
2765 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterEExo%s",add[tm].
Data()),
2766 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()),
2767 nptbins,ptmin,ptmax,200,-10,10);
2773 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterEExo%s",add[tm].
Data()),
2774 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()),
2775 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2781 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxEExo%s",add[tm].
Data()),
2782 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()),
2783 nptbins,ptmin,ptmax,200,-10,10);
2789 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxEExo%s",add[tm].
Data()),
2790 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()),
2791 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2797 for(
Int_t i = 0; i < 6; i++)
2799 for(
Int_t j = 0; j < 6; j++)
2802 (Form(
"hLambda0TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2803 Form(
"#lambda^{2}_{0} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
2804 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2810 (Form(
"hLambda1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2811 Form(
"#lambda^{2}_{1} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
2812 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2818 (Form(
"hLambda0NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2819 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()),
2820 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2826 (Form(
"hLambda1NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2827 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()),
2828 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2834 (Form(
"hLambda0NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2835 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()),
2836 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2842 (Form(
"hLambda1NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2843 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()),
2844 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2859 (Form(
"hNLocMaxTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2860 Form(
"#it{n}_{LM} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
2861 nptbins,ptmin,ptmax,10,0,10);
2867 (Form(
"hEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2868 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()),
2869 nptbins,ptmin,ptmax,100,0,1);
2875 (Form(
"hEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2876 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()),
2877 nptbins,ptmin,ptmax,100,0,1);
2883 (Form(
"hEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2884 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()),
2885 nptbins,ptmin,ptmax,100,0,1);
2891 (Form(
"hE2ndRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2892 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()),
2893 nptbins,ptmin,ptmax,100,0,1);
2899 (Form(
"hE2ndRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2900 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()),
2901 nptbins,ptmin,ptmax,100,0,1);
2907 (Form(
"hE2ndRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2908 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()),
2909 nptbins,ptmin,ptmax,100,0,1);
2915 (Form(
"hE2ndEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2916 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()),
2917 nptbins,ptmin,ptmax,100,0,1);
2923 (Form(
"hE2ndEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2924 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()),
2925 nptbins,ptmin,ptmax,100,0,1);
2931 (Form(
"hE2ndEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2932 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()),
2933 nptbins,ptmin,ptmax,100,0,1);
2940 (Form(
"hECellClusRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2941 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()),
2942 nptbins,ptmin,ptmax,100,0,1);
2948 (Form(
"hECellClusRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2949 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()),
2950 nptbins,ptmin,ptmax,100,0,1);
2956 (Form(
"hECellClusRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2957 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()),
2958 nptbins,ptmin,ptmax,100,0,1);
2964 (Form(
"hLogECellNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2965 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()),
2966 nptbins,ptmin,ptmax,150,-3,3);
2972 (Form(
"hLogECellNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2973 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()),
2974 nptbins,ptmin,ptmax,150,-3,3);
2980 (Form(
"hLogECellNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2981 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()),
2982 nptbins,ptmin,ptmax,150,-3,3);
2988 (Form(
"hECellWeightNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2989 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()),
2990 nptbins,ptmin,ptmax,90,0,4.5);
2996 (Form(
"hECellWeightNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2997 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()),
2998 nptbins,ptmin,ptmax,90,0,4.5);
3004 (Form(
"hECellWeightNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3005 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()),
3006 nptbins,ptmin,ptmax,90,0,4.5);
3012 (Form(
"hMassEClusTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3013 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()),
3014 nptbins,ptmin,ptmax,nmassbins,massmin,massmax);
3028 (Form(
"hTimeDiffTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3029 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()),
3030 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
3036 (Form(
"hColRowTCardCorrelNCellLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3037 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3038 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3044 (Form(
"hColRowTCardCorrelNCellHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3045 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3046 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3054 (Form(
"hExoticTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3055 Form(
"exoticity vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3056 nptbins,ptmin,ptmax,200,-1,1);
3062 (Form(
"hTimeDiffExoTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3063 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()),
3064 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
3070 (Form(
"hColRowTCardCorrelNCellExoticLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3071 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3072 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3078 (Form(
"hColRowTCardCorrelNCellExoticHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3079 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3080 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3295 (Form(
"hLambda0Exoticity_EBin%d%s",i,add[tm].
Data()),
3296 Form(
"#lambda^{2}_{0} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3297 200,-1,1,ssbins,ssmin,ssmax);
3303 (Form(
"hLambda1Exoticity_EBin%d%s",i,add[tm].
Data()),
3304 Form(
"#lambda^{2}_{1} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3305 200,-1,1,ssbins,ssmin,ssmax);
3319 (Form(
"hNCellsExoticity_EBin%d%s",i,add[tm].
Data()),
3320 Form(
"#it{n}_{cells} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3321 200,-1,1,nceclbins,nceclmin,nceclmax);
3322 fhNCellsExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
3327 (Form(
"hTimeExoticity_EBin%d%s",i,add[tm].
Data()),
3329 200,-1,1,100,-25,25);
3330 fhTimeExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
3335 (Form(
"hNCellsTCardSameAndDiff_Exotic_EBin%d%s",i,add[tm].
Data()),
3336 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()),
3337 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
3377 (Form(
"hNCellsTCardSameAndDiff_EBin%d%s",i,add[tm].
Data()),
3378 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()),
3379 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
3385 (Form(
"hLambda0Lambda1_EBin%d%s",i,add[tm].
Data()),
3386 Form(
"#lambda^{2}_{0} vs #lambda^{2}_{1}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3387 ssbins,ssmin,ssmax,ssbins,ssmin,ssmax);
3403 for(
Int_t j = 0; j < 6; j++)
3405 for(
Int_t k = 0; k < 6; k++)
3408 (Form(
"hLambda0Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
3409 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()),
3410 200,-1,1,ssbins,ssmin,ssmax);
3416 (Form(
"hLambda1Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
3417 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()),
3418 200,-1,1,ssbins,ssmin,ssmax);
3515 for(
Int_t i = 0; i < 12; i++)
3518 (Form(
"hTCardCorrECellMaxDiff_Case%d%s",i,add[tm].
Data()),
3519 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()),
3520 nptbins,ptmin,ptmax,210,-1,20);
3526 (Form(
"hTCardCorrEClusterDiff_Case%d%s",i,add[tm].
Data()),
3527 Form(
"#it{E}_{cluster}-#it{E}_{cell} vs #it{E}_{cluster}, for (un)correlated cells in TCard, case %d %s",i,add[tm].
Data()),
3528 nptbins,ptmin,ptmax,210,-1,20);
3550 (Form(
"hTCardCorrTCellMaxDiff_Case%d%s",i,add[tm].
Data()),
3551 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()),
3552 nptbins,ptmin,ptmax,1000,-100,100);
3560 (Form(
"hTCardCorrECellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
3561 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()),
3562 nptbins,ptmin,ptmax,210,-1,20);
3568 (Form(
"hTCardCorrEClusterDiffExo_Case%d%s",i,add[tm].
Data()),
3569 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()),
3570 nptbins,ptmin,ptmax,210,-1,20);
3592 (Form(
"hTCardCorrTCellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
3593 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()),
3594 nptbins,ptmin,ptmax,1000,-100,100);
3605 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3611 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3617 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3623 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3629 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3635 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3641 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3647 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3655 (Form(
"hTMPhiResidual_EBin%d",i),
3656 Form(
"#Delta #varphi_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3657 200,-1,1,nresphibins,resphimin,resphimax);
3663 (Form(
"hTMEtaResidual_EBin%d",i),
3664 Form(
"#Delta #eta_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3665 200,-1,1,nresetabins,resetamin,resetamax);
3671 (Form(
"hTMPhiResidual_LooseCut_EBin%d",i),
3672 Form(
"#Delta #varphi_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3673 200,-1,1,nresphibins,resphimin,resphimax);
3679 (Form(
"hTMEtaResidual_LooseCut_EBin%d",i),
3680 Form(
"#Delta #eta_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3681 200,-1,1,nresetabins,resetamin,resetamax);
3709 (Form(
"hColRowM02"),
3711 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax,40,0.,2.);
3718 (Form(
"hColRowM02NCellCut"),
3719 Form(
"column vs row vs M02, %2.2f<#it{E}<%2.2f GeV #it{n}_{cells}^{w>0.01} > %d for Neutral",
3721 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax,40,0.,2.);
3729 Form(
"%2.2f<#it{E}_{1}<%2.2f GeV, %2.2f<#it{E}_{2}<%2.2f GeV, %2.2f<#lambda^{2}_{0}<%2.2f"
3730 "#it{M}_{#gamma #gamma} vs #it{n}_{1, cells}^{w>0.01} vs SM number trig cluster",
3739 (
"hInvMassNCellSMSame",
3740 Form(
"%2.2f<#it{E}_{1}<%2.2f GeV, %2.2f<#it{E}_{2}<%2.2f GeV, %2.2f<#lambda^{2}_{0}<%2.2f"
3741 "#it{M}_{#gamma #gamma} vs #it{n}_{1, cells}^{w>0.01} vs SM number both cluster",
3752 (
"hEMaxCellTimeM02SM",
3753 Form(
"time vs #lambda_{0}^{2} vs SM number, "
3762 (
"hEMaxCellTimeNCellSM",
3763 Form(
"time vs #it{n}_{cells}^{w>0.01} vs SM number, "
3772 (
"hESecCellTimeNCellSM",
3773 Form(
"secondary cell time vs #it{n}_{cells}^{w>0.01} vs SM number, "
3783 Int_t nbinsdeltacells = 19 ;
3786 for(
Int_t col = 0; col < 2; col++)
3789 (Form(
"hColRowFromCellMaxLowM02_Col%d",col),
3790 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, "
3792 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells);
3798 (Form(
"hColRowFromCellMaxHighM02_Col%d",col),
3799 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, "
3801 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells);
3812 (Form(
"hNCellsPerClusterM02M20_SM%d",i),
3813 Form(
" vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01} vs #lambda_{1}^{2}, "
3815 (
Int_t)nShShBins/1.5,0,maxShSh/1.5, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3822 (Form(
"hNCellsPerClusterM02NLM_SM%d",i),
3823 Form(
" vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01} vs N_{lm}, "
3825 10,0,10, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3832 (Form(
"hEMaxCellEClusterM02NCell_SM%d",i),
3833 Form(
"(#it{E}_{cluster} - #it{E}_{max cell})/#it{E}_{cluster} vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01}, "
3835 50,0,1., cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3842 (Form(
"hEMaxCellLogM02NCell_SM%d",i),
3843 Form(
"log(#it{E}_{max cell}) vs #lambda_{0}^{2} vs SM number, "
3845 150,-3,3, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3852 (Form(
"hESecCellEMaxCellM02NCell_SM%d",i),
3853 Form(
"#it{E}_{cell}/#it{E}_{cell max} vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01}, "
3855 50,0,1., cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3862 (Form(
"hESecCellEClusterM02NCell_SM%d",i),
3863 Form(
"(#it{E}_{cluster} - #it{E}_{cell})/#it{E}_{cluster} vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01}, "
3865 50,0,1., cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3872 (Form(
"hESecCellLogM02NCell_SM%d",i),
3873 Form(
"log(#it{E}_{cell}) vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01}, "
3875 150,-3,3, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3882 (Form(
"hEMaxESecCellNCellLowM02_SM%d",i),
3883 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, "
3885 24,0,12, 48,0,12,cellBins,cellMin,cellMax);
3892 (Form(
"hEMaxECrossNCellLowM02_SM%d",i),
3893 Form(
"#it{E}_{cell}^{max} vs exoticity vs #it{n}_{cells}^{w>0.01} vs 0.1 < #lambda_{0}^{2} < 0.3, "
3895 24,0,12, 40,0.6,1.,cellBins,cellMin,cellMax);
3902 (Form(
"hEMaxESecCellNCellHighM02_SM%d",i),
3903 Form(
"#it{E}_{cell}^{max} vs #it{E}_{cell}^{secondary} vs #it{n}_{cells}^{w>0.01} vs 0.5 < #lambda_{0}^{2} < 2, "
3905 24,0,12, 48,0,12,cellBins,cellMin,cellMax);
3912 (Form(
"hEMaxECrossNCellHighM02_SM%d",i),
3913 Form(
"#it{E}_{cell}^{max} vs exoticity vs #it{n}_{cells}^{w>0.01} for 0.5 < #lambda_{0}^{2} < 2, "
3915 24,0,12, 40,0.6,1.,cellBins,cellMin,cellMax);
3921 Int_t nbinsdeltacells = 19 ;
3924 for(
Int_t col = 0; col < 2; col++)
3927 (Form(
"hNCellsPerClusterM02_SM%d_Col%d",i,col),
3928 Form(
"#it{E} vs #it{n}_{cells} vs #lambda_{0}^{2}, SM=%d, colum %d",i,col),
3929 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3936 (Form(
"hNCellsPerClusterM20_SM%d_Col%d",i,col),
3937 Form(
"#it{E} vs #it{n}_{cells} vs #lambda_{1}^{2}, SM=%d, colum %d",i,col),
3938 nEbins,minE,maxE,cellBins,cellMin,cellMax,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
3945 (Form(
"hNCellsPerClusterNLM_SM%d_Col%d",i,col),
3946 Form(
"#it{E} vs #it{n}_{cells} vs n_{lm}, SM=%d, colum %d",i,col),
3947 nEbins,minE,maxE,cellBins,cellMin,cellMax,10,0,10);
3954 (Form(
"hColRowFromCellMaxLowM02_SM%d_Col%d",i,col),
3955 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, "
3957 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,cellBins,cellMin,cellMax);
3964 (Form(
"hColRowFromCellMaxHighM02_SM%d_Col%d",i,col),
3965 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, "
3967 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,cellBins,cellMin,cellMax);
3973 for(
Int_t j = 0; j < 3; j++)
3976 (Form(
"hColRowFromCellMaxEMaxSecDiffLowM02_SM%d_Col%d_NCellBin%d",i,col,j),
3977 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, "
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,120,0,12);
3986 (Form(
"hColRowFromCellMaxEMaxSecDiffHighM02_SM%d_Col%d_NCellBin%d",i,col,j),
3987 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, "
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,120,0,12);
3996 (Form(
"hColRowFromCellMaxEMaxSecDiffFracLowM02_SM%d_Col%d_NCellBin%d",i,col,j),
3997 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, "
3998 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
3999 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,40,0,1);
4006 (Form(
"hColRowFromCellMaxEMaxSecDiffFracHighM02_SM%d_Col%d_NCellBin%d",i,col,j),
4007 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, "
4008 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
4009 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,40,0,1);
4016 (Form(
"hColRowFromCellMaxECellClusterRatLowM02_SM%d_Col%d_NCellBin%d",i,col,j),
4017 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, "
4018 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
4019 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,50,0,1);
4026 (Form(
"hColRowFromCellMaxECellClusterRatHighM02_SM%d_Col%d_NCellBin%d",i,col,j),
4027 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, "
4028 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
4029 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,50,0,1);
4040 for(
Int_t imatch = 0; imatch < 3; imatch++)
4045 (Form(
"hDeltaIEtaDeltaIPhi_%s",matchCase[imatch].
Data()),
4046 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),
4047 nEbins,minE,maxE,cellBins,cellMin,cellMax,cellBins,cellMin,cellMax);
4054 (Form(
"hDeltaIEtaDeltaIPhiTot_%s",matchCase[imatch].
Data()),
4055 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),
4056 nEbins,minE,maxE,cellBins,cellMin,cellMax,cellBins,cellMin,cellMax);
4063 (Form(
"hDeltaIA_%s",matchCase[imatch].
Data()),
4064 Form(
"Cluster *asymmetry* in cell units vs E for %s",matchCase[imatch].
Data()),
4065 nptbins,ptmin,ptmax,asyBins,asyMin,asyMax);
4066 fhDeltaIA[imatch]->SetXTitle(
"#it{E}_{cluster}");
4067 fhDeltaIA[imatch]->SetYTitle(
"#it{A}_{cell in cluster}");
4068 outputContainer->Add(
fhDeltaIA[imatch]);
4071 (Form(
"hDeltaIATot_%s",matchCase[imatch].
Data()),
4072 Form(
"Cluster *total asymmetry* in cell units vs E for %s",matchCase[imatch].
Data()),
4073 nptbins,ptmin,ptmax,asyBins,asyMin,asyMax);
4075 fhDeltaIATot[imatch]->SetYTitle(
"#it{A}_{cell in cluster}^{total}");
4079 (Form(
"hDeltaIAM02_%s",matchCase[imatch].
Data()),
4080 Form(
"Cluster *asymmetry* in cell units vs #lambda^{2}_{0} for %s",matchCase[imatch].
Data()),
4081 nEbins,minE,maxE,nShShBins,minShSh,maxShSh,asyBins,asyMin,asyMax);
4084 fhDeltaIAM02[imatch]->SetZTitle(
"#it{A}_{cell in cluster}");
4088 (Form(
"hDeltaIATotM02_%s",matchCase[imatch].
Data()),
4089 Form(
"Cluster *total asymmetry* in cell units vs #lambda^{2}_{0} for %s",matchCase[imatch].
Data()),
4090 nEbins,minE,maxE,nShShBins,minShSh,maxShSh,asyBins,asyMin,asyMax);
4093 fhDeltaIATotM02[imatch]->SetZTitle(
"#it{A}_{cell in cluster}^{total}");
4097 (Form(
"hDeltaIAM20_%s",matchCase[imatch].
Data()),
4098 Form(
"Cluster *asymmetry* in cell units vs #lambda^{2}_{1} for %s",matchCase[imatch].
Data()),
4099 nEbins,minE,maxE,(
Int_t)nShShBins/1.5,0,maxShSh/1.5,asyBins,asyMin,asyMax);
4102 fhDeltaIAM20[imatch]->SetZTitle(
"#it{A}_{cell in cluster}");
4106 (Form(
"hDeltaIATotM20_%s",matchCase[imatch].
Data()),
4107 Form(
"Cluster *total asymmetry* in cell units vs #lambda^{2}_{1} for %s",matchCase[imatch].
Data()),
4108 nEbins,minE,maxE,(
Int_t)nShShBins/1.5,0,maxShSh/1.5,asyBins,asyMin,asyMax);
4111 fhDeltaIATotM20[imatch]->SetZTitle(
"#it{A}_{cell in cluster}^{total}");
4115 (Form(
"hDeltaIANCells_%s",matchCase[imatch].
Data()),
4116 Form(
"Cluster *asymmetry* in cell units vs N cells in cluster for %s",matchCase[imatch].
Data()),
4117 nEbins,minE,maxE,cellBins,cellMin,cellMax,asyBins,asyMin,asyMax);
4124 (Form(
"hDeltaIATotNCells_%s",matchCase[imatch].
Data()),
4125 Form(
"Cluster *total asymmetry* in cell units vs N cells in cluster for %s",matchCase[imatch].
Data()),
4126 nEbins,minE,maxE,cellBins,cellMin,cellMax,asyBins,asyMin,asyMax);
4135 (Form(
"hDeltaIAOrigin_%s",matchCase[imatch].
Data()),
4136 Form(
"Cluster *asymmetry* in cell units vs E vs primary ID, for ID %s",matchCase[imatch].
Data()),
4137 nEbins,minE,maxE,10,-0.5,9.5,asyBins,asyMin,asyMax);
4141 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4142 fhDeltaIAOrigin[imatch]->GetYaxis()->SetBinLabel(ilabel,mcParticleStringLabel[ilabel-1]);
4146 (Form(
"hDeltaIATotOrigin_%s",matchCase[imatch].
Data()),
4147 Form(
"Cluster *total asymmetry* in cell units vs E vs primary ID, for ID %s",matchCase[imatch].
Data()),
4148 nEbins,minE,maxE,10,-0.5,9.5,asyBins,asyMin,asyMax);
4152 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4153 fhDeltaIATotOrigin[imatch]->GetYaxis()->SetBinLabel(ilabel,mcParticleStringLabel[ilabel-1]);
4159 (Form(
"hClusterMaxCellDiffM02_%s",matchCase[imatch].
Data()),
4160 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()),
4161 nEbins,minE,maxE, 20,0,1.,nShShBins,minShSh,maxShSh);
4168 (Form(
"hClusterTimeEnergyM02_%s",matchCase[imatch].
Data()),
4169 Form(
"#it{E} vs TOF vs #lambda_{0}^{2} for ID %s",matchCase[imatch].
Data()),
4170 nEbins,minE,maxE,45,-25.5,20.5,nShShBins,minShSh,maxShSh);
4177 (Form(
"hNCellsPerClusterM02_%s",matchCase[imatch].
Data()),
4178 Form(
"#it{E} vs #it{n}_{cells} vs #lambda_{0}^{2} for ID %s",matchCase[imatch].
Data()),
4179 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4186 (Form(
"hNCellsPerClusterM20_%s",matchCase[imatch].
Data()),
4187 Form(
"#it{E} vs #it{n}_{cells} vs #lambda_{1}^{2} for ID %s",matchCase[imatch].
Data()),
4188 nEbins,minE,maxE,cellBins,cellMin,cellMax,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4195 (Form(
"hNCellsPerClusterNLM_%s",matchCase[imatch].
Data()),
4196 Form(
"#it{E} vs #it{n}_{cells} vs n_{lm} for ID %s",matchCase[imatch].
Data()),
4197 nEbins,minE,maxE,cellBins,cellMin,cellMax,10,0,10);
4204 (Form(
"hSMM02_%s",matchCase[imatch].
Data()),
4205 Form(
"#it{E} vs SM number vs #lambda_{0}^{2}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4206 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,minShSh,maxShSh);
4207 fhSMM02[imatch]->SetXTitle(
"#it{E} (GeV)");
4208 fhSMM02[imatch]->SetYTitle(
"SM number");
4209 fhSMM02[imatch]->SetZTitle(
"#lambda_{0}^{2}");
4210 outputContainer->Add(
fhSMM02[imatch]);
4213 (Form(
"hSMM02NoCut_%s",matchCase[imatch].
Data()),
4214 Form(
"#it{E} vs SM number vs #lambda_{0}^{2} for ID %s",matchCase[imatch].
Data()),
4215 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,minShSh,maxShSh);
4222 (Form(
"hSMM20LowM02_%s",matchCase[imatch].
Data()),
4223 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",
4225 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4232 (Form(
"hSMM20LowM02NoCut_%s",matchCase[imatch].
Data()),
4233 Form(
"#it{E} vs SM number vs #lambda_{1}^{2}, for 0.1<#lambda_{0}^{2}<0.3, for ID %s",matchCase[imatch].
Data()),
4234 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4241 (Form(
"hSMM20HighM02_%s",matchCase[imatch].
Data()),
4242 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",
4244 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4251 (Form(
"hSMM20HighM02NoCut_%s",matchCase[imatch].
Data()),
4252 Form(
"#it{E} vs SM number vs #lambda_{1}^{2}, for 0.5<#lambda_{0}^{2}<2, for ID %s",matchCase[imatch].
Data()),
4253 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4260 (Form(
"hSMNCell_%s",matchCase[imatch].
Data()),
4261 Form(
"#it{E} vs SM number vs #it{n}_{cells}^{w>0.01} for ID %s",matchCase[imatch].
Data()),
4262 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,cellBins,cellMin,cellMax);
4263 fhSMNCell[imatch]->SetXTitle(
"#it{E} (GeV)");
4264 fhSMNCell[imatch]->SetYTitle(
"SM number");
4265 fhSMNCell[imatch]->SetZTitle(
"#it{n}_{cells}^{w>0.01}");
4266 outputContainer->Add(
fhSMNCell[imatch]);
4269 (Form(
"hSMNCellM02_%s",matchCase[imatch].
Data()),
4270 Form(
"SM number vs #it{n}_{cells}^{w>0.01} vs #lambda_{0}^{2}, "
4272 fNModules,-0.5,fNModules-0.5,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4275 fhSMNCellM02[imatch]->SetYTitle(
"#it{n}_{cells}^{w>0.01}");
4279 (Form(
"hColM02_%s",matchCase[imatch].
Data()),
4280 Form(
"#it{E} vs column number vs #lambda_{0}^{2}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4281 nEbins,minE,maxE,48,-0.5,47.5,nShShBins,minShSh,maxShSh);
4282 fhColM02[imatch]->SetXTitle(
"#it{E} (GeV)");
4283 fhColM02[imatch]->SetYTitle(
"column number");
4284 fhColM02[imatch]->SetZTitle(
"#lambda_{0}^{2}");
4285 outputContainer->Add(
fhColM02[imatch]);
4288 (Form(
"hRowM02_%s",matchCase[imatch].
Data()),
4289 Form(
"#it{E} vs row number vs #lambda_{0}^{2}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4290 nEbins,minE,maxE,24,-0.5,23.5,nShShBins,minShSh,maxShSh);
4291 fhRowM02[imatch]->SetXTitle(
"#it{E} (GeV)");
4292 fhRowM02[imatch]->SetYTitle(
"row number");
4293 fhRowM02[imatch]->SetZTitle(
"#lambda_{0}^{2}");
4294 outputContainer->Add(
fhRowM02[imatch]);
4299 (Form(
"hNCellsPerClusterMEta_%s",matchCase[imatch].
Data()),
4300 Form(
"#it{E} vs #it{n}_{cells} vs #sigma_{#eta}^{2} for ID %s",matchCase[imatch].
Data()),
4301 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4308 (Form(
"hNCellsPerClusterMPhi_%s",matchCase[imatch].
Data()),
4309 Form(
"#it{E} vs #it{n}_{cells} vs #sigma_{#varphi}^{2} for ID %s",matchCase[imatch].
Data()),
4310 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4317 (Form(
"hNCellsPerClusterMEtaPhi_%s",matchCase[imatch].
Data()),
4318 Form(
"#it{E} vs #it{n}_{cells} vs #sigma_{#eta#varphi}^{2} for ID %s",matchCase[imatch].
Data()),
4319 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,-1*maxShSh,maxShSh);
4326 (Form(
"hNCellsPerClusterMEtaPhiA_%s",matchCase[imatch].
Data()),
4327 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()),
4328 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,-1*maxShSh,maxShSh);
4331 fhNCellsPerClusterMEtaPhiA[imatch]->SetZTitle(
"(#sigma_{#varphi}^{2}-#sigma_{#eta}^{2})/(#sigma_{#varphi}^{2}+#sigma_{#eta}^{2})");
4338 (Form(
"hOrigin_%s",matchCase[imatch].
Data()),
4339 Form(
"#it{E} vs origin for ID %s",matchCase[imatch].
Data()),
4340 nEbins,minE,maxE, 10,-0.5,9.5);
4341 fhOriginE[imatch]->SetXTitle(
"#it{E} (GeV)");
4342 fhOriginE[imatch]->SetYTitle(
"particle");
4343 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4344 fhOriginE[imatch]->GetYaxis()->SetBinLabel(ilabel ,mcParticleStringLabel[ilabel-1]);
4345 outputContainer->Add(
fhOriginE[imatch]);
4348 (Form(
"hOriginM02_%s",matchCase[imatch].
Data()),
4349 Form(
"#it{E} vs origin vs #lambda_{0}^{2} for ID %s",matchCase[imatch].
Data()),
4350 nEbins,minE,maxE,10,-0.5,9.5,nShShBins,minShSh,maxShSh);
4353 fhOriginM02[imatch]->SetZTitle(
"#lambda_{0}^{2}");
4354 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4355 fhOriginM02[imatch]->GetYaxis()->SetBinLabel(ilabel,mcParticleStringLabel[ilabel-1]);
4370 (
"hClusterMaxCellCloseCellRatioM02",
"#it{E} vs #it{E}_{cell-i}/#it{E}_{cell max} vs #lambda_{0}^{2}",
4371 nEbins,minE,maxE, 20,0,1.,nShShBins,minShSh,maxShSh);
4389 (
"hClusterMaxCellECrossM02",
4390 "#it{E} vs exoticity vs #lambda_{0}^{2}",
4391 nEbins,minE,maxE, 40,0.6,1.,nShShBins,minShSh,maxShSh);
4402 nptbins,ptmin,ptmax, 100,0,1.);
4408 nptbins,ptmin,ptmax, 100,-10,0);
4414 nptbins,ptmin,ptmax, 100,0,1.);
4420 nptbins,ptmin,ptmax, 100,-10,0);
4426 nptbins*2,ptmin,ptmax*2, 100,0,1.);
4432 nptbins*2,ptmin,ptmax*2, 100,-10,0);
4442 if(imod < fFirstModule || imod >
fLastModule)
continue;
4445 Form(
"#cell energy / sum all energy vs all energy in Module %d",imod),
4446 nptbins*2,ptmin,ptmax*2, 100,0,1.);
4452 Form(
"Log(cell energy / sum all energy) vs all energy in Module %d",imod),
4453 nptbins*2,ptmin,ptmax*2, 100,-10,0);
4460 Float_t cellEmin [] = {0.05,0.1,0.15,0.2};
4461 Float_t cellTmin [] = {50.,100.,1000000.};
4463 for(
Int_t iw = 0; iw < 12; iw++)
4466 for(
Int_t iEmin = 0; iEmin < 4; iEmin++)
4468 for(
Int_t iTmin = 0; iTmin < 3; iTmin++)
4471 Form(
"#lambda^{2}_{0} vs E, w0=%1.2f, cell E>%2.2f MeV, |t|<%2.0f ns",
4472 w0, cellEmin[iEmin], cellTmin[iTmin]),
4473 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4488 Form(
"#lambda^{2}_{0} vs E, w0=%1.2f, cell E>%2.2f MeV, |t|<%2.0f ns, |#eta| < 0.15",
4489 w0, cellEmin[iEmin], cellTmin[iTmin]),
4490 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4499 TString mcnames[] = {
"Photon",
"Electron",
"Conversion",
"Pi0",
"Hadron"};
4500 for(
Int_t imc = 0; imc < 5; imc++)
4503 Form(
"shower shape, #lambda^{2}_{0} vs E, w0 = %1.1f, for MC %s",w0,mcnames[imc].
Data()),
4504 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4523 return outputContainer;
4584 AliInfo(Form(
"dEdX cuts init for MC: %d<dedx_ele<%d; %d<dedx_had<%d",
4592 if ( runNumber < 146861 )
4599 else if ( runNumber < 156000 )
4606 else if ( runNumber < 165000 )
4621 AliInfo(Form(
"dEdX cuts init for run %d: %d<dedx_ele<%d; %d<dedx_had<%d",
4635 Int_t nCellsPerCluster)
4653 if(ampMax < 0.01)
return kFALSE;
4669 printf(
"**** Print %s %s ****\n", GetName(), GetTitle() ) ;
4673 printf(
"Min n cells : %d\n" ,
fNCellMin) ;
4675 printf(
"Min M02 : %1.2f\n",
fM02Min) ;
4691 AliDebug(1,
"Start");
4707 AliFatal(Form(
"AliAnaClusterShapeCorrelStudies::MakeAnalysisFillHistograms() - Wrong calorimeter name <%s>, END",
GetCalorimeterString().
Data()));
4711 AliFatal(Form(
"AliAnaClusterShapeCorrelStudies::MakeAnalysisFillHistograms() - No CaloClusters or CaloCells available"));
4741 if(clus->GetNCells() < 3)
return ;
4743 for (
Int_t ipos = 0; ipos < clus->GetNCells(); ipos++)
4745 Int_t id = clus->GetCellsAbsId()[ipos];
4760 AliWarning(Form(
"Wrong calculated energy %f",energy));
4771 for (
Int_t ipos = 0; ipos < clus->GetNCells(); ipos++)
4773 Int_t id = clus->GetCellsAbsId()[ipos];
4788 Float_t l0org = clus->GetM02();
4789 Float_t l1org = clus->GetM20();
4790 Float_t dorg = clus->GetDispersion();
4794 if(
IsDataMC() && clus->GetNLabels() > 0)
4819 Float_t cellEmin [] = {0.05,0.1,0.15,0.2};
4820 Float_t cellTmin [] = {50.,100.,1000000.};
4822 for(
Int_t iw = 0; iw < 12; iw++)
4824 for(
Int_t iEmin = 0; iEmin < 4; iEmin++)
4826 for(
Int_t iTmin = 0; iTmin < 3; iTmin++)
4832 Float_t disp = 0, dEta = 0, dPhi = 0;
4833 Float_t sEta = 0, sPhi = 0, sEtaPhi = 0;
4837 cellEmin[iEmin], cellTmin[iTmin], bc,
4838 newEnergy, l0, l1, disp, dEta, dPhi,
4839 sEta, sPhi, sEtaPhi);
4860 clus->SetM02(l0org);
4861 clus->SetM20(l1org);
4862 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