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++)
133 for(
Int_t i =0; i < 12; i++)
135 for(
Int_t j = 0; j < 4; j++)
137 for(
Int_t k = 0; k < 3; k++)
144 for(
Int_t j = 0; j < 5; j++)
152 for(
Int_t tm = 0; tm < 2; tm++)
248 for(
Int_t i = 0; i < 6; i++)
250 for(
Int_t j = 0; j < 6; j++)
323 for(
Int_t i = 0; i < 14; i++)
342 for(
Int_t j = 0; j < 6; j++)
344 for(
Int_t k = 0; k < 6; k++)
365 for(
Int_t i = 0; i < 12; i++)
391 for(
Int_t i = 0; i < 20; i++)
412 for(
Int_t j = 0; j < 3; j++)
456 Int_t icol = -1, irow = -1, iRCU = -1, icolAbs = -1, irowAbs = -1;
457 GetModuleNumberCellIndexesAbsCaloMap(absIdMax,GetCalorimeter(), icol, irow, iRCU, icolAbs, irowAbs);
461 for(
Int_t ie = 0; ie < fNEBinCuts; ie++)
463 if( energy >= fEBinCuts[ie] && energy < fEBinCuts[ie+1] ) ebin = ie;
466 Int_t ncells = clus->GetNCells();
467 Double_t time = clus->GetTOF()*1.e9;
468 if(time>400) time-=fConstantTimeShift;
470 Float_t deta = clus->GetTrackDz();
471 Float_t dphi = clus->GetTrackDx();
479 if(energy >= 5 && energy < 8)
480 fhColRowExoticLowE1Cell [matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
483 fhColRowExoticHighE1Cell[matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
485 if ( time > 5) fhColRowExoticHighE1CellPosTime->Fill(icolAbs,irowAbs,GetEventWeight());
486 else if( time < -5) fhColRowExoticHighE1CellNegTime->Fill(icolAbs,irowAbs,GetEventWeight());
487 else fhColRowExoticHighE1CellNulTime->Fill(icolAbs,irowAbs,GetEventWeight());
490 fhEnergyTime1Cell[matched]->Fill(energy,time,GetEventWeight());
491 fhEnergyTMEtaResidual1Cell->Fill(energy,deta,GetEventWeight());
492 fhEnergyTMPhiResidual1Cell->Fill(energy,dphi,GetEventWeight());
494 else if(exoticity > 0.97)
496 if(energy >= 5 && energy < 8)
497 fhColRowExoticLowE [matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
500 fhColRowExoticHighE[matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
502 if ( time > 5) fhColRowExoticHighEPosTime->Fill(icolAbs,irowAbs,GetEventWeight());
503 else if( time < -5) fhColRowExoticHighENegTime->Fill(icolAbs,irowAbs,GetEventWeight());
504 else fhColRowExoticHighENulTime->Fill(icolAbs,irowAbs,GetEventWeight());
507 fhEnergyTimeExotic[matched]->Fill(energy,time,GetEventWeight());
508 fhEnergyTMEtaResidualExotic->Fill(energy,deta,GetEventWeight());
509 fhEnergyTMPhiResidualExotic->Fill(energy,dphi,GetEventWeight());
511 for (
Int_t ipos = 0; ipos < ncells; ipos++)
513 Int_t absId = clus->GetCellsAbsId()[ipos];
515 Float_t eCell = fCaloCellList->GetCellAmplitude(absId);
516 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
519 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
521 if( absId == absIdMax || weight < 0.01 )
continue;
523 Int_t rowDiff = -100, colDiff = -100;
524 Bool_t sameTCard = GetCaloUtils()->IsAbsIDsFromTCard(absIdMax,absId,rowDiff,colDiff);
527 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
528 GetModuleNumberCellIndexesAbsCaloMap(absId,GetCalorimeter(), icol2, irow2, iRCU2, icolAbs2, irowAbs2);
532 if(energy >= 5 && energy < 8)
533 fhColRowExotic2ndCellDiffLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
535 fhColRowExotic2ndCellDiffHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
539 if(energy >= 5 && energy < 8)
540 fhColRowExotic2ndCellSameLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
542 fhColRowExotic2ndCellSameHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
547 else if ( energy > 8 )
549 if ( time > 5) fhColRowHighEPosTime->Fill(icolAbs,irowAbs,GetEventWeight());
550 else if( time < -5) fhColRowHighENegTime->Fill(icolAbs,irowAbs,GetEventWeight());
551 else fhColRowHighENulTime->Fill(icolAbs,irowAbs,GetEventWeight());
556 fhTMPhiResidualExoticityLooseCut[ebin]->Fill(exoticity,dphi);
557 fhTMEtaResidualExoticityLooseCut[ebin]->Fill(exoticity,deta);
560 else if ( energy > 8 && ncells > 1 && exoticity < 0.97)
562 if ( time > 5) fhColRowHighEPosTime->Fill(icolAbs,irowAbs,GetEventWeight());
563 else if( time < -5) fhColRowHighENegTime->Fill(icolAbs,irowAbs,GetEventWeight());
564 else fhColRowHighENulTime->Fill(icolAbs,irowAbs,GetEventWeight());
570 Int_t etaRegion = -1, phiRegion = -1;
571 GetCaloUtils()->GetEMCALSubregion(clus,fCaloCellList,etaRegion,phiRegion);
573 if ( etaRegion !=0 ) return ;
579 fhEnergyTimeTCardCorrNoSelection1Cell[matched]->Fill(energy,time,GetEventWeight());
580 fhEnergyTMEtaResidualTCardCorrNoSelection1Cell->Fill(energy,deta,GetEventWeight());
581 fhEnergyTMPhiResidualTCardCorrNoSelection1Cell->Fill(energy,dphi,GetEventWeight());
583 else if(exoticity > 0.97)
585 fhEnergyTimeTCardCorrNoSelectionExotic[matched]->Fill(energy,time,GetEventWeight());
586 fhEnergyTMEtaResidualTCardCorrNoSelectionExotic->Fill(energy,deta,GetEventWeight());
587 fhEnergyTMPhiResidualTCardCorrNoSelectionExotic->Fill(energy,dphi,GetEventWeight());
594 if ( m02 > 0.001 ) lamR = m20/m02;
596 Int_t absIdList[ncells];
598 Int_t nlm = GetCaloUtils()->GetNumberOfLocalMaxima(clus, fCaloCellList, absIdList, maxEList) ;
604 Int_t nCellWithWeight = 1;
612 Float_t eCellMax = fCaloCellList->GetCellAmplitude(absIdMax);
613 Double_t tCellMax = fCaloCellList->GetCellTime(absIdMax);
615 GetCaloUtils()->RecalibrateCellAmplitude(eCellMax, GetCalorimeter(), absIdMax);
616 GetCaloUtils()->RecalibrateCellTime(tCellMax, GetCalorimeter(), absIdMax, GetReader()->GetInputEvent()->GetBunchCrossNumber());
620 if(tCellMax>400) tCellMax-=fConstantTimeShift;
631 Bool_t sameTCard2ndMax = kFALSE;
635 Float_t eCellSameRowSameTCardNearCol = 0.;
636 Float_t eCellSameRowDiffTCardNearCol = 0.;
637 Double_t tCellSameRowSameTCardNearCol = 0.;
638 Double_t tCellSameRowDiffTCardNearCol = 0.;
645 for (
Int_t ipos = 0; ipos < ncells; ipos++)
647 Int_t absId = clus->GetCellsAbsId()[ipos];
649 Float_t eCell = fCaloCellList->GetCellAmplitude(absId);
650 Double_t tCell = fCaloCellList->GetCellTime(absId);
652 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
653 GetCaloUtils()->RecalibrateCellTime(tCell, GetCalorimeter(), absId, GetReader()->GetInputEvent()->GetBunchCrossNumber());
655 if(tCell>400) tCell-=fConstantTimeShift;
658 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
660 if( absId == absIdMax || weight < 0.01 )
continue;
664 fhECellClusRatNLM1TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
665 fhLogECellNLM1TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
666 fhECellWeightNLM1TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
670 fhECellClusRatNLM2TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
671 fhLogECellNLM2TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
672 fhECellWeightNLM2TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
676 fhECellClusRatNLM3TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
677 fhLogECellNLM3TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
678 fhECellWeightNLM3TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
681 Int_t rowDiff = -100, colDiff = -100;
682 Bool_t sameTCard = GetCaloUtils()->IsAbsIDsFromTCard(absIdMax,absId,rowDiff,colDiff);
688 fhECellSameClusRatNLM1TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
689 fhLogECellSameNLM1TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
690 fhECellSameWeightNLM1TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
694 fhECellSameClusRatNLM2TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
695 fhLogECellSameNLM2TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
696 fhECellSameWeightNLM2TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
700 fhECellSameClusRatNLM3TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
701 fhLogECellSameNLM3TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
702 fhECellSameWeightNLM3TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
712 Int_t indexType = -1;
717 if(TMath::Abs(rowDiff) == 1) nearRow = kTRUE;
718 if(TMath::Abs(colDiff) == 1) nearCol = kTRUE;
720 if ( rowDiff == 0 && colDiff != 0 )
722 if ( nearCol ) indexType = 6;
728 else if ( rowDiff != 0 && colDiff == 0 )
730 if ( nearRow ) indexType = 8;
738 if ( nearRow && nearCol) indexType = 10;
749 if ( rowDiff == 0 && colDiff != 0 )
751 if ( nearCol ) indexType = 0;
755 else if ( rowDiff != 0 && colDiff == 0 )
757 if ( nearRow ) indexType = 2;
762 if ( nearCol && nearRow ) indexType = 4;
767 if ( rowDiff == 0 && TMath::Abs(colDiff) == 1 )
771 eCellSameRowSameTCardNearCol = eCell;
772 tCellSameRowSameTCardNearCol = tCell;
776 eCellSameRowDiffTCardNearCol = eCell;
777 tCellSameRowDiffTCardNearCol = tCell;
783 Float_t eCellDiff = eCellMax - eCell;
784 Float_t eClusDiff = energy - eCell;
787 Float_t tCellDiff = tCellMax - tCell;
789 fhTCardCorrECellMaxDiff[indexType][matched]->Fill(energy, eCellDiff, GetEventWeight());
790 fhTCardCorrEClusterDiff[indexType][matched]->Fill(energy, eClusDiff, GetEventWeight());
793 fhTCardCorrTCellMaxDiff[indexType][matched]->Fill(energy, tCellDiff, GetEventWeight());
795 if ( fStudyExotic && exoticity > 0.97 )
797 fhTCardCorrECellMaxDiffExo[indexType][matched]->Fill(energy, eCellDiff, GetEventWeight());
798 fhTCardCorrEClusterDiffExo[indexType][matched]->Fill(energy, eClusDiff, GetEventWeight());
801 fhTCardCorrTCellMaxDiffExo[indexType][matched]->Fill(energy, tCellDiff, GetEventWeight());
805 if ( fStudyExotic && exoticity > 0.97 )
808 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
809 GetModuleNumberCellIndexesAbsCaloMap(absId,GetCalorimeter(), icol2, irow2, iRCU2, icolAbs2, irowAbs2);
813 if(energy >= 5 && energy < 8)
814 fhColRowTCardCorrNoSelectionExotic2ndCellDiffLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
816 fhColRowTCardCorrNoSelectionExotic2ndCellDiffHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
820 if(energy >= 5 && energy < 8)
821 fhColRowTCardCorrNoSelectionExotic2ndCellSameLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
823 fhColRowTCardCorrNoSelectionExotic2ndCellSameHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
827 if ( eCell > emax2nd )
833 sameTCard2ndMax = kTRUE;
839 sameTCard2ndMax = kFALSE;
874 Float_t ratioNcells = nCellWithWeight/(ncells*1.);
875 fhNCellsTCardCorrRatioWithWeightNoSelection[matched]->Fill(energy, ratioNcells, GetEventWeight());
879 if ( nCellWithWeight <= 1 )
return;
884 printf(
"AliAnaClusterShapeCorrelStudies: M02 %f, M20 %f, E %2.3f, ncell %d, n with weight %d; max cell E %2.3f\n",
885 m02,m20,energy,ncells,nCellWithWeight,eCellMax);
894 if(energy >= 5 && energy < 8)
895 fhColRowTCardCorrNoSelectionLowE [matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
897 fhColRowTCardCorrNoSelectionHighE[matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
899 if ( fStudyExotic && exoticity > 0.97 )
901 if(energy >= 5 && energy < 8)
902 fhColRowTCardCorrNoSelectionExoticLowE [matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
904 fhColRowTCardCorrNoSelectionExoticHighE[matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
907 fhTimeTCardCorrNoSelection [matched]->Fill(energy,time, GetEventWeight());
908 fhLambda0TCardCorrNoSelection[matched]->Fill(energy, m02, GetEventWeight());
909 fhLambda1TCardCorrNoSelection[matched]->Fill(energy, m20, GetEventWeight());
913 fhLambda0NLM1TCardCorrNoSelection[matched]->Fill(energy, m02, GetEventWeight());
914 fhLambda1NLM1TCardCorrNoSelection[matched]->Fill(energy, m20, GetEventWeight());
918 fhLambda0NLM2TCardCorrNoSelection[matched]->Fill(energy, m02, GetEventWeight());
919 fhLambda1NLM2TCardCorrNoSelection[matched]->Fill(energy, m20, GetEventWeight());
922 fhLambdaRTCardCorrNoSelection[matched]->Fill(energy,lamR, GetEventWeight());
923 fhNLocMaxTCardCorrNoSelection[matched]->Fill(energy, nlm, GetEventWeight());
924 fhExoticTCardCorrNoSelection [matched]->Fill(energy, exoticity, GetEventWeight());
928 fhEMaxRatNLM1TCardCorrNoSelection [matched]->Fill(energy, eCellMax/energy , GetEventWeight());
929 fhE2ndRatNLM1TCardCorrNoSelection [matched]->Fill(energy, emax2nd/energy , GetEventWeight());
930 fhE2ndEMaxRatNLM1TCardCorrNoSelection[matched]->Fill(energy, emax2nd/eCellMax, GetEventWeight());
934 fhEMaxRatNLM2TCardCorrNoSelection [matched]->Fill(energy, eCellMax/energy , GetEventWeight());
935 fhE2ndRatNLM2TCardCorrNoSelection [matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
936 fhE2ndEMaxRatNLM2TCardCorrNoSelection[matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
940 fhEMaxRatNLM3TCardCorrNoSelection [matched]->Fill(energy, eCellMax/energy , GetEventWeight());
941 fhE2ndRatNLM3TCardCorrNoSelection [matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
942 fhE2ndEMaxRatNLM3TCardCorrNoSelection[matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
949 fhE2ndSameRatNLM1TCardCorrNoSelection [matched]->Fill(energy, emax2nd/energy , GetEventWeight());
950 fhE2ndSameEMaxRatNLM1TCardCorrNoSelection[matched]->Fill(energy, emax2nd/eCellMax, GetEventWeight());
954 fhE2ndSameRatNLM2TCardCorrNoSelection [matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
955 fhE2ndSameEMaxRatNLM2TCardCorrNoSelection[matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
959 fhE2ndSameRatNLM3TCardCorrNoSelection [matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
960 fhE2ndSameEMaxRatNLM3TCardCorrNoSelection[matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
964 fhNCellsTCardCorrNoSelection [matched]->Fill(energy, ncells, GetEventWeight());
965 fhNCellsTCardCorrWithWeightNoSelection[matched]->Fill(energy, nCellWithWeight, GetEventWeight());
967 if(eCellSameRowSameTCardNearCol > 0 && eCellSameRowDiffTCardNearCol > 0)
969 Float_t eDiff = eCellSameRowSameTCardNearCol - eCellSameRowDiffTCardNearCol ;
970 Float_t tDiff = tCellSameRowSameTCardNearCol - tCellSameRowDiffTCardNearCol ;
972 fhSameRowDiffColAndTCardCellsEnergyDiffClusterE[matched]->Fill(energy , eDiff, GetEventWeight());
973 fhSameRowDiffColAndTCardCellsTimeDiffClusterE [matched]->Fill(energy , tDiff, GetEventWeight());
974 fhSameRowDiffColAndTCardCellsEnergyDiffCellMaxE[matched]->Fill(eCellMax, eDiff, GetEventWeight());
975 fhSameRowDiffColAndTCardCellsTimeDiffCellMaxE [matched]->Fill(eCellMax, tDiff, GetEventWeight());
977 if ( fStudyExotic && exoticity > 0.97 )
979 fhSameRowDiffColAndTCardCellsEnergyDiffClusterEExo[matched]->Fill(energy , eDiff, GetEventWeight());
980 fhSameRowDiffColAndTCardCellsTimeDiffClusterEExo [matched]->Fill(energy , tDiff, GetEventWeight());
981 fhSameRowDiffColAndTCardCellsEnergyDiffCellMaxEExo[matched]->Fill(eCellMax, eDiff, GetEventWeight());
982 fhSameRowDiffColAndTCardCellsTimeDiffCellMaxEExo [matched]->Fill(eCellMax, tDiff, GetEventWeight());
987 Int_t nCorrInd = nCorr;
988 if(nCorr > 4) nCorrInd = 5;
990 Int_t nCorrNoInd = nCorrNo;
991 if(nCorrNoInd > 4) nCorrNoInd = 5;
1004 fhLambda0NLM1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m02, GetEventWeight());
1005 fhLambda1NLM1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m20, GetEventWeight());
1009 fhLambda0NLM2TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m02, GetEventWeight());
1010 fhLambda1NLM2TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m20, GetEventWeight());
1013 fhLambda0TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m02, GetEventWeight());
1014 fhLambda1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m20, GetEventWeight());
1016 fhNLocMaxTCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, nlm, GetEventWeight());
1019 fhExoticTCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, exoticity, GetEventWeight());
1023 fhEMaxRatNLM1TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, eCellMax/energy , GetEventWeight());
1024 fhE2ndRatNLM1TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, emax2nd/energy , GetEventWeight());
1025 fhE2ndEMaxRatNLM1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, emax2nd/eCellMax, GetEventWeight());
1029 fhEMaxRatNLM2TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, eCellMax/energy, GetEventWeight());
1030 fhE2ndRatNLM2TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
1031 fhE2ndEMaxRatNLM2TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
1035 fhEMaxRatNLM3TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, eCellMax/energy, GetEventWeight());
1036 fhE2ndRatNLM3TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
1037 fhE2ndEMaxRatNLM3TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
1041 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1043 Int_t absId = clus->GetCellsAbsId()[ipos];
1045 Float_t eCell = fCaloCellList->GetCellAmplitude(absId);
1046 Double_t tCell = fCaloCellList->GetCellTime(absId);
1048 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
1049 GetCaloUtils()->RecalibrateCellTime(tCell, GetCalorimeter(), absId, GetReader()->GetInputEvent()->GetBunchCrossNumber());
1051 if(tCell>400) tCell-=fConstantTimeShift;
1054 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
1056 if( absId == absIdMax || weight < 0.01 )
continue;
1058 Float_t tDiffMaxSecondary = tCellMax - tCell;
1060 fhTimeDiffTCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, tDiffMaxSecondary, GetEventWeight());
1061 if ( fStudyExotic && exoticity > 0.97 )
1063 fhTimeDiffExoTCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, tDiffMaxSecondary, GetEventWeight());
1065 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
1066 GetModuleNumberCellIndexesAbsCaloMap(absId,GetCalorimeter(), icol2, irow2, iRCU2, icolAbs2, irowAbs2);
1068 Int_t rowDiff = -100, colDiff = -100;
1069 Bool_t sameTCard = GetCaloUtils()->IsAbsIDsFromTCard(absIdMax,absId,rowDiff,colDiff);
1073 if ( energy >= 5 && energy < 8 && nCorr == 0 )
1074 fhColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
1075 else if ( energy >= 8 && nCorr == 0 )
1076 fhColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
1080 if ( energy >= 5 && energy < 8 && nCorrNo == 0 )
1081 fhColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
1082 else if ( energy >= 8 && nCorrNo == 0 )
1083 fhColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
1089 fhECellClusRatNLM1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, eCell/energy, GetEventWeight());
1090 fhLogECellNLM1TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
1091 fhECellWeightNLM1TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, weight, GetEventWeight());
1095 fhECellClusRatNLM2TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, eCell/energy, GetEventWeight());
1096 fhLogECellNLM2TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
1097 fhECellWeightNLM2TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, weight, GetEventWeight());
1101 fhECellClusRatNLM3TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, eCell/energy, GetEventWeight());
1102 fhLogECellNLM3TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
1103 fhECellWeightNLM3TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, weight, GetEventWeight());
1108 if(m02 > fInvMassMinM02Cut && m02 < fInvMassMaxM02Cut)
1110 for(
Int_t jclus = 0 ; jclus < fCaloClusList->GetEntriesFast() ; jclus++)
1112 AliVCluster* clus2 = (AliVCluster*) fCaloClusList->At(jclus);
1115 Int_t absIdMax2 = GetCaloUtils()->GetMaxEnergyCell(fCaloCellList, clus2, maxCellFraction);
1117 Double_t tof2 = clus2->GetTOF()*1.e9;
1118 if(tof2>400) tof2-=fConstantTimeShift;
1125 if( absIdMax == absIdMax2
1126 || !IsGoodCluster(absIdMax2, clus2->GetM02(), clus2->GetNCells())
1127 || GetCaloUtils()->GetNumberOfLocalMaxima(clus2,fCaloCellList) > 1
1128 || clus2->GetM02() > fInvMassMaxM02Cut
1129 || clus2->GetM02() < fInvMassMinM02Cut
1130 || clus2->E() < fInvMassMinECut
1131 || clus2->E() > fInvMassMaxECut
1132 || TMath::Abs(diffTof) > fInvMassMaxTimeDifference
1137 clus2->GetMomentum(fClusterMomentum2,v);
1141 if(IsFiducialCutOn()) in2 = GetFiducialCut()->IsInFiducialCut(fClusterMomentum2.Eta(),fClusterMomentum2.Phi(),GetCalorimeter()) ;
1147 Double_t angle = fClusterMomentum.Angle(fClusterMomentum2.Vect());
1149 if( angle > fInvMassMaxOpenAngle )
continue;
1152 Float_t mass = (fClusterMomentum+fClusterMomentum2).M ();
1153 fhMassEClusTCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, mass, GetEventWeight());
1158 if ( energy >= 5 && energy < 8)
1159 fhColRowTCardCorrelNCellLowE [nCorrInd][nCorrNoInd][matched]->Fill(icolAbs, irowAbs, GetEventWeight());
1160 else if ( energy >= 8 )
1161 fhColRowTCardCorrelNCellHighE[nCorrInd][nCorrNoInd][matched]->Fill(icolAbs, irowAbs, GetEventWeight());
1176 if ( fStudyExotic && exoticity > 0.97 )
1178 if ( energy >= 5 && energy < 8)
1179 fhColRowTCardCorrelNCellExoticLowE [nCorrInd][nCorrNoInd][matched]->Fill(icolAbs, irowAbs, GetEventWeight());
1180 else if ( energy >= 8 )
1181 fhColRowTCardCorrelNCellExoticHighE[nCorrInd][nCorrNoInd][matched]->Fill(icolAbs, irowAbs, GetEventWeight());
1218 fhLambda0Lambda1 [ebin][matched]->Fill(m20, m02, GetEventWeight());
1219 fhNCellsTCardSameAndDiff[ebin][matched]->Fill(nCorrNo, nCorr, GetEventWeight());
1225 fhLambda0Exoticity[ebin][matched]->Fill(exoticity, m02, GetEventWeight());
1226 fhLambda1Exoticity[ebin][matched]->Fill(exoticity, m20, GetEventWeight());
1228 fhNCellsExoticity [ebin][matched]->Fill(exoticity, nCellWithWeight, GetEventWeight());
1229 fhTimeExoticity [ebin][matched]->Fill(exoticity, tCellMax, GetEventWeight());
1233 fhLambda0ExoticityPerNCell[nCorrInd][nCorrNoInd][matched]->Fill(exoticity, m02, GetEventWeight());
1234 fhLambda1ExoticityPerNCell[nCorrInd][nCorrNoInd][matched]->Fill(exoticity, m20, GetEventWeight());
1246 if ( exoticity > 0.97 )
1247 fhNCellsTCardSameAndDiffExotic[ebin][matched]->Fill(nCorrNo, nCorr, GetEventWeight());
1250 fhTMPhiResidualExoticity[ebin]->Fill(exoticity,dphi);
1251 fhTMEtaResidualExoticity[ebin]->Fill(exoticity,deta);
1261 Float_t nCellRat = nCorr*1. / ((nCorr+nCorrNo)*1.);
1262 fhNCellsTCardSameAndDiffFraction[matched]->Fill(energy, nCellRat, GetEventWeight());
1263 if ( fStudyExotic && exoticity > 0.97 )
1264 fhNCellsTCardSameAndDiffFractionExotic[matched]->Fill(energy, nCellRat, GetEventWeight());
1363 Int_t nCaloCellsPerCluster = clus->GetNCells();
1364 if ( nCaloCellsPerCluster < 2 )
return;
1379 for (
Int_t ipos = 0; ipos < nCaloCellsPerCluster; ipos++)
1381 Int_t absId = clus ->GetCellsAbsId()[ipos];
1382 Float_t eCell = fCaloCellList->GetCellAmplitude(absId) ;
1384 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
1386 if( absId == absIdMax || eCell < 0.01 )
continue;
1388 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
1390 if( weight < 0.01 )
continue;
1395 if ( nCell < 1 )
return;
1400 Int_t ietaMax=-1, iphiMax = 0, rcuMax = 0, icolAbs = -1, irowAbs = -1;
1401 Int_t smMax = GetModuleNumberCellIndexesAbsCaloMap(absIdMax, GetCalorimeter(),
1402 ietaMax, iphiMax, rcuMax, icolAbs, irowAbs);
1404 if ( matchedPID == 0 && energy > fEMinShape && energy < fEMaxShape )
1406 fhColRowM02->Fill(icolAbs,irowAbs,m02,GetEventWeight()) ;
1408 if ( nCell > fNCellMinShape )
1409 fhColRowM02NCellCut->Fill(icolAbs,irowAbs,m02,GetEventWeight()) ;
1417 Int_t etaRegion = -1, phiRegion = -1;
1418 GetCaloUtils()->GetEMCALSubregion(clus,fCaloCellList,etaRegion,phiRegion);
1420 if ( etaRegion !=0 ) return ;
1424 for (
Int_t ipos = 0; ipos < nCaloCellsPerCluster; ipos++)
1426 Int_t absId = clus ->GetCellsAbsId()[ipos];
1427 Float_t eCell = fCaloCellList->GetCellAmplitude(absId) ;
1429 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
1431 if( absId == absIdMax || eCell < 0.01 )
continue;
1433 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
1435 if( weight < 0.01 )
continue;
1437 Float_t fracCell = eCell/eCellMax;
1438 fhClusterMaxCellCloseCellRatioM02->Fill(energy, fracCell, m02, GetEventWeight());
1440 Float_t fracClus = (energy-eCell)/energy;
1442 Float_t logECell = TMath::Log(eCell);
1453 if ( energy > fEMinShape && energy < fEMaxShape && matchedPID == 0 )
1455 fhESecCellEMaxCellM02NCellPerSM[smMax]->Fill(fracCell, nCell, m02, GetEventWeight());
1456 fhESecCellEClusterM02NCellPerSM[smMax]->Fill(fracClus, nCell, m02, GetEventWeight());
1457 fhESecCellLogM02NCellPerSM [smMax]->Fill(logECell, nCell, m02, GetEventWeight());
1459 if ( m02 > 0.1 && m02 < 0.3 )
1461 fhEMaxESecCellNCellLowM02PerSM [smMax]->Fill(eCellMax, eCell , nCell, GetEventWeight());
1462 fhEMaxECrossNCellLowM02PerSM [smMax]->Fill(eCellMax, eCrossFrac, nCell, GetEventWeight());
1464 Double_t time = fCaloCellList->GetCellTime(absId);
1465 GetCaloUtils()->RecalibrateCellTime(time, GetCalorimeter(), absId,GetReader()->GetInputEvent()->GetBunchCrossNumber());
1467 if(time>400) time-=fConstantTimeShift;
1469 fhESecCellTimeNCellSM->Fill(time, smMax, nCell, GetEventWeight());
1471 else if ( m02 > 0.5 && m02 < 2 )
1473 fhEMaxESecCellNCellHighM02PerSM[smMax]->Fill(eCellMax, eCell , nCell, GetEventWeight());
1474 fhEMaxECrossNCellHighM02PerSM [smMax]->Fill(eCellMax, eCrossFrac, nCell, GetEventWeight());
1484 Int_t sm = GetModuleNumberCellIndexes(absId,GetCalorimeter(), ieta, iphi, rcu);
1486 if(dIphi < TMath::Abs(iphi-iphiMax)) dIphi = TMath::Abs(iphi-iphiMax);
1487 if(iphi-iphiMax < 0 && dIphiNeg > iphi-iphiMax) dIphiNeg = iphi-iphiMax;
1488 if(iphi-iphiMax > 0 && dIphiPos < iphi-iphiMax) dIphiPos = iphi-iphiMax;
1492 if(dIeta < TMath::Abs(ieta-ietaMax)) dIeta = TMath::Abs(ieta-ietaMax);
1493 if(ieta-ietaMax < 0 && dIetaNeg > ieta-ietaMax) dIetaNeg = ieta-ietaMax;
1494 if(ieta-ietaMax > 0 && dIetaPos < ieta-ietaMax) dIetaPos = ieta-ietaMax;
1496 if ( energy > fEMinShape && energy < fEMaxShape && matchedPID == 0 )
1499 if ( nCell == 2 || nCell == 3 ) nCellBin = 0;
1500 else if ( nCell == 4 || nCell == 5 ) nCellBin = 1;
1501 if ( m02 > 0.1 && m02 < 0.3 )
1503 fhColRowFromCellMaxLowM02PerSM[smMax][ietaMax%2]->Fill(ietaMax-ieta, iphiMax-iphi, nCell, GetEventWeight());
1504 fhColRowFromCellMaxLowM02 [ietaMax%2]->Fill(ietaMax-ieta, iphiMax-iphi, GetEventWeight());
1506 fhColRowFromCellMaxEMaxSecDiffFracLowM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi,(eCellMax-eCell)/eCellMax, GetEventWeight());
1507 fhColRowFromCellMaxEMaxSecDiffLowM02PerSM [smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi, eCellMax-eCell, GetEventWeight());
1508 fhColRowFromCellMaxECellClusterRatLowM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi, eCell/energy , GetEventWeight());
1510 else if (m02 > 0.5 && m02 < 2)
1512 fhColRowFromCellMaxHighM02PerSM[smMax][ietaMax%2]->Fill(ietaMax-ieta, iphiMax-iphi, nCell, GetEventWeight());
1513 fhColRowFromCellMaxHighM02 [ietaMax%2]->Fill(ietaMax-ieta, iphiMax-iphi, GetEventWeight());
1515 fhColRowFromCellMaxEMaxSecDiffFracHighM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi,(eCellMax-eCell)/eCellMax, GetEventWeight());
1516 fhColRowFromCellMaxEMaxSecDiffHighM02PerSM [smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi, eCellMax-eCell, GetEventWeight());
1517 fhColRowFromCellMaxECellClusterRatHighM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi, eCell/energy , GetEventWeight());
1524 Int_t ietaShift = ieta;
1525 Int_t ietaMaxShift = ietaMax;
1527 if (ieta > ietaMax) ietaMaxShift+=48;
1528 else ietaShift +=48;
1530 if(dIeta < TMath::Abs(ietaShift-ietaMaxShift)) dIeta = TMath::Abs(ietaShift-ietaMaxShift);
1531 if(ietaShift-ietaMaxShift < 0 && dIetaNeg > ietaShift-ietaMaxShift) dIetaNeg = ietaShift-ietaMaxShift;
1532 if(ietaShift-ietaMaxShift > 0 && dIetaPos < ietaShift-ietaMaxShift) dIetaPos = ietaShift-ietaMaxShift;
1539 if ( energy > fEMinShape && energy < fEMaxShape && matchedPID == 0 )
1541 fhNCellsPerClusterM02M20PerSM[smMax]->Fill(m20, nCell, m02, GetEventWeight());
1544 if ( nCell == 2 || nCell == 3 ) nCellBin = 0;
1545 else if ( nCell == 4 || nCell == 5 ) nCellBin = 1;
1546 if ( m02 > 0.1 && m02 < 0.3 )
1547 fhColRowFromCellMaxECellClusterRatLowM02PerSM [smMax][ietaMax%2][nCellBin]->Fill(0., 0., eCellMax/energy, GetEventWeight());
1548 else if (m02 > 0.5 && m02 < 2)
1549 fhColRowFromCellMaxECellClusterRatHighM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(0., 0., eCellMax/energy, GetEventWeight());
1553 fhClusterMaxCellECrossM02->Fill(energy, eCrossFrac, m02, GetEventWeight());
1558 fhClusterMaxCellDiffM02[matchedPID]->Fill(energy, maxFrac, m02, GetEventWeight());
1559 fhClusterTimeEnergyM02 [matchedPID]->Fill(energy, tmax , m02, GetEventWeight());
1560 fhNCellsPerClusterM02 [matchedPID]->Fill(energy, nCell , m02, GetEventWeight());
1561 fhNCellsPerClusterM20 [matchedPID]->Fill(energy, nCell , m20, GetEventWeight());
1563 fhSMNCell [matchedPID]->Fill(energy, smMax , nCell, GetEventWeight());
1564 fhSMM02NoCut [matchedPID]->Fill(energy, smMax , m02 , GetEventWeight());
1566 if ( energy > fEMinShape && energy < fEMaxShape )
1568 fhSMNCellM02[matchedPID]->Fill(smMax , nCell, m02, GetEventWeight());
1570 if ( matchedPID == 0 )
1572 if ( m02 > 0.1 && m02 < 0.3 )
1573 fhEMaxCellTimeNCellSM->Fill(tmax, smMax, nCell, GetEventWeight());
1575 fhEMaxCellTimeM02SM->Fill(tmax, smMax, m02, GetEventWeight());
1577 fhEMaxCellEClusterM02NCellPerSM[smMax]->Fill(maxFrac, nCell, m02, GetEventWeight());
1578 fhEMaxCellLogM02NCellPerSM [smMax]->Fill(TMath::Log(eCellMax), nCell, m02, GetEventWeight());
1582 if ( nCell > fNCellMinShape )
1584 fhSMM02 [matchedPID]->Fill(energy, smMax , m02, GetEventWeight());
1585 fhColM02[matchedPID]->Fill(energy, ietaMax, m02, GetEventWeight());
1586 fhRowM02[matchedPID]->Fill(energy, iphiMax, m02, GetEventWeight());
1590 if ( fStudyShapeParam )
1593 Float_t dIA = 1.*(dIphi-dIeta)/(dIeta+dIphi);
1594 Float_t dIATot = 1.*((dIphiPos-dIphiNeg)-(dIetaPos-dIetaNeg))/((dIetaPos-dIetaNeg)+(dIphiPos-dIphiNeg));
1603 AliDebug(1,Form(
"E %2.2f, nCell %d, dPhi %d, dEta %d, dIA %2.2f, match %d",energy,nCell, dIphi,dIeta,dIA,matchedPID));
1613 fhDeltaIANCells [matchedPID]->Fill(energy, nCell , dIA , GetEventWeight());
1614 fhDeltaIATotNCells [matchedPID]->Fill(energy, nCell , dIATot, GetEventWeight());
1616 if ( nCell > fNCellMinShape )
1618 fhDeltaIEtaDeltaIPhi[matchedPID]->Fill(energy, dIeta, dIphi, GetEventWeight());
1619 fhDeltaIA [matchedPID]->Fill(energy, dIA , GetEventWeight());
1620 fhDeltaIAM02 [matchedPID]->Fill(energy, m02 , dIA , GetEventWeight());
1621 fhDeltaIAM20 [matchedPID]->Fill(energy, m20 , dIA , GetEventWeight());
1623 fhDeltaIEtaDeltaIPhiTot[matchedPID]->Fill(energy,dIetaPos-dIetaNeg, dIphiPos-dIphiNeg, GetEventWeight());
1624 fhDeltaIATot [matchedPID]->Fill(energy, dIATot , GetEventWeight());
1625 fhDeltaIATotM02 [matchedPID]->Fill(energy, m02 , dIATot, GetEventWeight());
1626 fhDeltaIATotM20 [matchedPID]->Fill(energy, m20 , dIATot, GetEventWeight());
1629 if ( IsDataMC() && mcIndex > -1 && mcIndex < 10 && nCell > fNCellMinShape )
1631 fhDeltaIAOrigin [matchedPID]->Fill(energy, mcIndex, dIA , GetEventWeight());
1632 fhDeltaIATotOrigin[matchedPID]->Fill(energy, mcIndex, dIATot, GetEventWeight());
1636 Float_t dispp= 0., dEta = 0., dPhi = 0.;
1637 Float_t sEta = 0., sPhi = 0., sEtaPhi = 0.;
1638 if ( GetCalorimeter() ==
kEMCAL )
1640 GetCaloUtils()->GetEMCALRecoUtils()->RecalculateClusterShowerShapeParameters(GetEMCALGeometry(), fCaloCellList, clus,
1641 l0, l1, dispp, dEta, dPhi, sEta, sPhi, sEtaPhi);
1644 if(sEta+sPhi>0.0001) sEtaPhiA = (sPhi-sEta)/(sEta+sPhi);
1646 AliDebug(2,Form(
"Recalculate shower shape org: m02 %2.2f, m20 %2.2f, disp %2.2f;"
1647 " new: m02 %2.2f, m20 %2.2f, disp %2.2f; "
1648 "mEta %2.2f, mPhi %2.2f, mEtaPhi %2.2f, A_EtaPhi %2.2f; dEta %2.2f dPhi %2.2f",
1649 m02,m20,clus->GetDispersion(),l0,l1,dispp,sEta,sPhi,sEtaPhi,sEtaPhiA,dEta,dPhi));
1651 fhNCellsPerClusterMEta [matchedPID]->Fill(energy, nCell, sEta , GetEventWeight());
1652 fhNCellsPerClusterMPhi [matchedPID]->Fill(energy, nCell, sPhi , GetEventWeight());
1653 fhNCellsPerClusterMEtaPhi [matchedPID]->Fill(energy, nCell, sEtaPhi , GetEventWeight());
1654 fhNCellsPerClusterMEtaPhiA[matchedPID]->Fill(energy, nCell, sEtaPhiA, GetEventWeight());
1659 if ( matchedPID == 0 && energy > fEMinShape && energy < fEMaxShape &&
1660 m02 > fInvMassMinM02Cut && m02 < fInvMassMaxM02Cut )
1662 for(
Int_t jclus = 0 ; jclus < fCaloClusList->GetEntriesFast() ; jclus++)
1664 AliVCluster* clus2 = (AliVCluster*) fCaloClusList->At(jclus);
1667 Int_t absIdMax2 = GetCaloUtils()->GetMaxEnergyCell(fCaloCellList, clus2, maxCellFraction);
1669 Double_t tof2 = clus2->GetTOF()*1.e9;
1670 if(tof2>400) tof2-=fConstantTimeShift;
1677 if( absIdMax == absIdMax2
1678 || !IsGoodCluster(absIdMax2, clus2->GetM02(), clus2->GetNCells())
1679 || GetCaloUtils()->GetNumberOfLocalMaxima(clus2,fCaloCellList) > 1
1680 || clus2->GetM02() > fInvMassMaxM02Cut
1681 || clus2->GetM02() < fInvMassMinM02Cut
1682 || clus2->E() < fInvMassMinECut
1683 || clus2->E() > fInvMassMaxECut
1684 || TMath::Abs(diffTof) > fInvMassMaxTimeDifference
1689 clus2->GetMomentum(fClusterMomentum2,v);
1693 if(IsFiducialCutOn()) in2 = GetFiducialCut()->IsInFiducialCut(fClusterMomentum2.Eta(),fClusterMomentum2.Phi(),GetCalorimeter()) ;
1699 Double_t angle = fClusterMomentum.Angle(fClusterMomentum2.Vect());
1701 if( angle > fInvMassMaxOpenAngle )
continue;
1704 Float_t mass = (fClusterMomentum+fClusterMomentum2).M ();
1705 fhInvMassNCellSM->Fill(mass, nCell, smMax, GetEventWeight());
1707 Int_t smMax2 = GetModuleNumber(clus2);
1708 if(smMax == smMax2) fhInvMassNCellSMSame->Fill(mass, nCell, smMax, GetEventWeight());
1713 if ( IsDataMC() && mcIndex > -1 && mcIndex < 10)
1715 fhOriginE [matchedPID]->Fill(energy, mcIndex, GetEventWeight());
1716 fhOriginM02[matchedPID]->Fill(energy, mcIndex, m02, GetEventWeight());
1732 AliVTrack *track = GetCaloUtils()->GetMatchedTrack(clus, GetReader()->GetInputEvent());
1742 Double_t dedx = track->GetTPCsignal();
1757 if(fdEdXMinEle == 0 || fdEdXMaxEle == 0 || fdEdXMinHad == 0 || fdEdXMaxHad == 0)
1758 InitdEdXParameters();
1760 if ( dedx >= fdEdXMinEle && dedx < fdEdXMaxEle ) matchedPID = 1;
1761 else if ( dedx >= fdEdXMinHad && dedx < fdEdXMaxHad ) matchedPID = 2;
1764 AliDebug(1,Form(
"dEdX out of range %2.2f",dedx));
1782 Int_t nCaloCellsPerCluster = 0 ;
1793 for(
Int_t iclus = 0; iclus < nCaloClusters; iclus++)
1795 AliDebug(1,Form(
"Cluster: %d/%d, data %d",iclus+1,nCaloClusters,
GetReader()->GetDataType()));
1797 AliVCluster* clus = (AliVCluster*)
fCaloClusList->At(iclus);
1800 if ( clus->GetDistanceToBadChannel() <
fMinDistToBad )
return ;
1804 if ( nModule < fFirstModule || nModule >
fLastModule )
1806 AliDebug(1,Form(
"Cluster module out of range %d",nModule));
1816 Double_t tof = clus->GetTOF()*1.e9;
1834 AliDebug(1,Form(
"Remove cluster with phi %2.2f and eta %2.2f", phi*TMath::RadToDeg(), eta));
1838 AliDebug(1,Form(
"cluster: E %2.3f, pT %2.3f, eta %2.3f, phi %2.3f",e,pt,eta,phi*TMath::RadToDeg()));
1842 nCaloCellsPerCluster = clus->GetNCells();
1846 AliDebug(1,Form(
"Accept cluster? %d",goodCluster));
1848 if(!goodCluster)
continue;
1877 Int_t matchedPID = 0;
1894 if ( ampMax > 0.01 )
1902 if (
fStudyShape && matchedPID >= 0 && matchedPID < 3 )
1923 const Int_t buffersize = 255;
1924 char onePar[buffersize] ;
1926 snprintf(onePar,buffersize,
"--- AliAnaClusterShapeCorrelStudies ---:") ;
1947 return new TObjString(parList) ;
1957 outputContainer->SetName(
"ClusterShapeStudies") ;
1968 TString mcParticleStringLabel[] = {
"Merged #gamma#gamma",
"Decay #gamma",
"Direct #gamma",
"e^{#pm}",
"#pi^{#pm}",
"k^{#pm}",
"p",
"#bar{p}",
"n",
"#bar{n}"};
1972 TString matchCase[] = {
"Neutral",
"Electron",
"Hadron"};
2000 Float_t rowcellmin = fNMaxRowsFullMin-1.5;
2009 Int_t nShShBins = 200;
2019 Int_t cellBins = 15;
2029 (
"hColRowHighEPosTime",
2030 "column vs row, exo < 0.97, E > 8 GeV, t > 5 ns",
2031 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2037 (
"hColRowHighENegTime",
2038 "column vs row, exo < 0.97, E > 8 GeV, t < -5 ns",
2039 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2045 (
"hColRowHighENulTime",
2046 "column vs row, exo < 0.97, E > 8 GeV, -5 < t < 5 ns",
2047 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2055 (
"hColRowExoticHighE1CellPosTime",
2056 "column vs row, 1 cell, E > 8 GeV, t > 5 ns",
2057 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2063 (
"hColRowExoticHighEPosTime",
2064 "column vs row, exo > 0.97, E > 8 GeV, t > 5 ns",
2065 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2071 (
"hColRowExoticHighE1CellNegTime",
2072 "column vs row, 1 cell, E > 8 GeV, t < -5 ns",
2073 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2079 (
"hColRowExoticHighENegTime",
2080 "column vs row, exo > 0.97, E > 8 GeV, t < -5 ns",
2081 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2087 (
"hColRowExoticHighE1CellNulTime",
2088 "column vs row, 1 cell, E > 8 GeV, -5 < t < 5 ns",
2089 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2095 (
"hColRowExoticHighENulTime",
2096 "column vs row, exo > 0.97, E > 8 GeV, -5 < t < 5 ns",
2097 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2103 TString add[] = {
"",
"TrackMatched"};
2104 for(
Int_t tm = 0; tm < 2; tm++)
2107 (Form(
"hColRowTCardCorrNoSelectionLowE%s",add[tm].
Data()),
2108 Form(
"column vs row, max E cell for TCard correlation selected clusters, 5 < E < 8 GeV %s",add[tm].
Data()),
2109 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2115 (Form(
"hColRowTCardCorrNoSelectionHighE%s",add[tm].
Data()),
2116 Form(
"column vs row, max E cell for TCard correlation selected clusters, E > 8 GeV %s",add[tm].
Data()),
2117 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2125 (Form(
"hNCellsTCardCorrNoSelection%s",add[tm].
Data()),
2126 Form(
"# custer # cells vs #it{E} %s",add[tm].
Data()),
2127 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
2133 (Form(
"hNCellsTCardCorrWithWeightNoSelection%s",add[tm].
Data()),
2134 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
2135 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
2141 (Form(
"hNCellsTCardCorrRatioWithWeightNoSelection%s",add[tm].
Data()),
2142 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
2143 nptbins,ptmin,ptmax, 100,0,1);
2149 (Form(
"hTimeTCardCorrNoSelection%s",add[tm].
Data()),
2150 Form(
"#it{time} vs #it{E} %s",add[tm].
Data()),
2151 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2157 (Form(
"hLambda0TCardCorrNoSelection%s",add[tm].
Data()),
2158 Form(
"#lambda^{2}_{0} vs #it{E} %s",add[tm].
Data()),
2159 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2165 (Form(
"hLambda1TCardCorrNoSelection%s",add[tm].
Data()),
2166 Form(
"#lambda^{2}_{1} vs #it{E} %s",add[tm].
Data()),
2167 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2173 (Form(
"hLambda0NLM1TCardCorrNoSelection%s",add[tm].
Data()),
2174 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=1 %s",add[tm].
Data()),
2175 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2181 (Form(
"hLambda1NLM1TCardCorrNoSelection%s",add[tm].
Data()),
2182 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=1 %s",add[tm].
Data()),
2183 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2189 (Form(
"hLambda0NLM2TCardCorrNoSelection%s",add[tm].
Data()),
2190 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=2 %s",add[tm].
Data()),
2191 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2197 (Form(
"hLambda1NLM2TCardCorrNoSelection%s",add[tm].
Data()),
2198 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=2 %s",add[tm].
Data()),
2199 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2206 (Form(
"hLambdaRTCardCorrNoSelection%s",add[tm].
Data()),
2207 Form(
"#lambda^{1}_{0}/#lambda^{2}_{0} vs #it{E} %s",add[tm].
Data()),
2208 nptbins,ptmin,ptmax,110,0,1.1);
2214 (Form(
"hNLocMaxTCardCorrNoSelection%s",add[tm].
Data()),
2215 Form(
"#it{n}_{LM} vs E %s",add[tm].
Data()),
2216 nptbins,ptmin,ptmax,10,0,10);
2222 (Form(
"hEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2223 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
2224 nptbins,ptmin,ptmax,100,0,1);
2230 (Form(
"hEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2231 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2232 nptbins,ptmin,ptmax,100,0,1);
2238 (Form(
"hEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2239 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2240 nptbins,ptmin,ptmax,100,0,1);
2246 (Form(
"hE2ndRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2247 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
2248 nptbins,ptmin,ptmax,100,0,1);
2254 (Form(
"hE2ndRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2255 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2256 nptbins,ptmin,ptmax,100,0,1);
2262 (Form(
"hE2ndRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2263 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2264 nptbins,ptmin,ptmax,100,0,1);
2271 (Form(
"hE2ndEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2272 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
2273 nptbins,ptmin,ptmax,100,0,1);
2279 (Form(
"hE2ndEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2280 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2281 nptbins,ptmin,ptmax,100,0,1);
2287 (Form(
"hE2ndEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2288 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2289 nptbins,ptmin,ptmax,100,0,1);
2297 (Form(
"hE2ndSameRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2298 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()),
2299 nptbins,ptmin,ptmax,100,0,1);
2305 (Form(
"hE2ndSameRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2306 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()),
2307 nptbins,ptmin,ptmax,100,0,1);
2313 (Form(
"hE2ndSameRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2314 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()),
2315 nptbins,ptmin,ptmax,100,0,1);
2322 (Form(
"hE2ndSameEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2323 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()),
2324 nptbins,ptmin,ptmax,100,0,1);
2330 (Form(
"hE2ndSameEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2331 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()),
2332 nptbins,ptmin,ptmax,100,0,1);
2338 (Form(
"hE2ndSameEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2339 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()),
2340 nptbins,ptmin,ptmax,100,0,1);
2348 (Form(
"hECellClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2349 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
2350 nptbins,ptmin,ptmax,100,0,1);
2356 (Form(
"hECellClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2357 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2358 nptbins,ptmin,ptmax,100,0,1);
2364 (Form(
"hECellClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2365 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2366 nptbins,ptmin,ptmax,100,0,1);
2372 (Form(
"hLogECellNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2373 Form(
"log(#it{E}_{cell}) vs E_{cluster}, #it{n}_{LM}=1, w > 0.01 %s",add[tm].
Data()),
2374 nptbins,ptmin,ptmax,150,-3,3);
2380 (Form(
"hLogECellNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2381 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}=2, w > 0.01 %s",add[tm].
Data()),
2382 nptbins,ptmin,ptmax,150,-3,3);
2388 (Form(
"hLogECellNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2389 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}>2, w > 0.01 %s",add[tm].
Data()),
2390 nptbins,ptmin,ptmax,150,-3,3);
2397 (Form(
"hECellWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2398 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
2399 nptbins,ptmin,ptmax,90,0,4.5);
2405 (Form(
"hECellWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2406 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2407 nptbins,ptmin,ptmax,90,0,4.5);
2413 (Form(
"hECellWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2414 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster}))vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2415 nptbins,ptmin,ptmax,90,0,4.5);
2422 (Form(
"hECellSameClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2423 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()),
2424 nptbins,ptmin,ptmax,100,0,1);
2430 (Form(
"hECellSameClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2431 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, cell from same T-Card as leading, #it{n}_{LM}=2 %s",add[tm].
Data()),
2432 nptbins,ptmin,ptmax,100,0,1);
2438 (Form(
"hECellSameClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2439 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, cell from same T-Card as leading, #it{n}_{LM}>2 %s",add[tm].
Data()),
2440 nptbins,ptmin,ptmax,100,0,1);
2446 (Form(
"hLogECellSameNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2447 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()),
2448 nptbins,ptmin,ptmax,150,-3,3);
2454 (Form(
"hLogECellSameNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2455 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()),
2456 nptbins,ptmin,ptmax,150,-3,3);
2462 (Form(
"hLogECellSameNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2463 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()),
2464 nptbins,ptmin,ptmax,150,-3,3);
2471 (Form(
"hECellSameWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2472 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()),
2473 nptbins,ptmin,ptmax,90,0,4.5);
2479 (Form(
"hECellSameWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2480 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()),
2481 nptbins,ptmin,ptmax,90,0,4.5);
2487 (Form(
"hECellSameWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2488 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()),
2489 nptbins,ptmin,ptmax,90,0,4.5);
2495 (Form(
"hExoticTCardCorrNoSelection%s",add[tm].
Data()),
2496 Form(
"exoticity vs #it{E} %s",add[tm].
Data()),
2497 nptbins,ptmin,ptmax,200,-1,1);
2503 (Form(
"hNCellsTCardSameAndDiffFraction%s",add[tm].
Data()),
2504 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, %s",add[tm].
Data()),
2505 nptbins,ptmin,ptmax,100,0,1);
2511 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterE%s",add[tm].
Data()),
2512 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()),
2513 nptbins,ptmin,ptmax,200,-10,10);
2519 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterE%s",add[tm].
Data()),
2520 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()),
2521 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2527 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxE%s",add[tm].
Data()),
2528 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()),
2529 nptbins,ptmin,ptmax,200,-10,10);
2535 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxE%s",add[tm].
Data()),
2536 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()),
2537 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2545 (Form(
"hEnergyTime1Cell%s",add[tm].
Data()),
2546 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
2547 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2553 (Form(
"hEnergyTimeExotic%s",add[tm].
Data()),
2554 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
2555 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2561 (Form(
"hEnergyTimeTCardCorrNoSelection1Cell%s",add[tm].
Data()),
2562 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
2563 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2569 (Form(
"hEnergyTimeTCardCorrNoSelectionExotic%s",add[tm].
Data()),
2570 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
2571 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2577 (Form(
"hColRowExoticLowE1Cell%s",add[tm].
Data()),
2578 Form(
"column vs row, 1 cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2579 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2585 (Form(
"hColRowExoticHighE1Cell%s",add[tm].
Data()),
2586 Form(
"column vs row, 1 cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2587 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2593 (Form(
"hColRowExoticLowE%s",add[tm].
Data()),
2594 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2595 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2601 (Form(
"hColRowExoticHighE%s",add[tm].
Data()),
2602 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2603 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2609 (Form(
"hColRowExotic2ndCellDiffLowE%s",add[tm].
Data()),
2610 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2611 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2617 (Form(
"hColRowExotic2ndCellDiffHighE%s",add[tm].
Data()),
2618 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2619 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2625 (Form(
"hColRowExotic2ndCellSameLowE%s",add[tm].
Data()),
2626 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2627 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2633 (Form(
"hColRowExotic2ndCellSameHighE%s",add[tm].
Data()),
2634 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2635 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2641 (Form(
"hColRowTCardCorrNoSelectionExoticLowE%s",add[tm].
Data()),
2642 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2643 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2649 (Form(
"hColRowTCardCorrNoSelectionExoticHighE%s",add[tm].
Data()),
2650 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2651 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2657 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffLowE%s",add[tm].
Data()),
2658 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2659 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2665 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffHighE%s",add[tm].
Data()),
2666 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2667 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2673 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameLowE%s",add[tm].
Data()),
2674 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2675 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2681 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameHighE%s",add[tm].
Data()),
2682 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2683 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2689 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameLowE%s",add[tm].
Data()),
2690 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2691 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2697 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameHighE%s",add[tm].
Data()),
2698 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2699 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2705 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffLowE%s",add[tm].
Data()),
2706 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2707 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2713 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffHighE%s",add[tm].
Data()),
2714 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2715 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2721 (Form(
"hNCellsTCardSameAndDiffFraction_Exotic%s",add[tm].
Data()),
2722 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, exo > 0.97 %s",add[tm].
Data()),
2723 nptbins,ptmin,ptmax,100,0,1);
2729 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterEExo%s",add[tm].
Data()),
2730 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()),
2731 nptbins,ptmin,ptmax,200,-10,10);
2737 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterEExo%s",add[tm].
Data()),
2738 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()),
2739 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2745 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxEExo%s",add[tm].
Data()),
2746 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()),
2747 nptbins,ptmin,ptmax,200,-10,10);
2753 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxEExo%s",add[tm].
Data()),
2754 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()),
2755 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2761 for(
Int_t i = 0; i < 6; i++)
2763 for(
Int_t j = 0; j < 6; j++)
2766 (Form(
"hLambda0TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2767 Form(
"#lambda^{2}_{0} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
2768 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2774 (Form(
"hLambda1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2775 Form(
"#lambda^{2}_{1} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
2776 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2782 (Form(
"hLambda0NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2783 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()),
2784 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2790 (Form(
"hLambda1NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2791 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()),
2792 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2798 (Form(
"hLambda0NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2799 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()),
2800 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2806 (Form(
"hLambda1NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2807 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()),
2808 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2823 (Form(
"hNLocMaxTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2824 Form(
"#it{n}_{LM} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
2825 nptbins,ptmin,ptmax,10,0,10);
2831 (Form(
"hEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2832 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()),
2833 nptbins,ptmin,ptmax,100,0,1);
2839 (Form(
"hEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2840 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()),
2841 nptbins,ptmin,ptmax,100,0,1);
2847 (Form(
"hEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2848 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()),
2849 nptbins,ptmin,ptmax,100,0,1);
2855 (Form(
"hE2ndRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2856 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()),
2857 nptbins,ptmin,ptmax,100,0,1);
2863 (Form(
"hE2ndRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2864 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()),
2865 nptbins,ptmin,ptmax,100,0,1);
2871 (Form(
"hE2ndRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2872 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()),
2873 nptbins,ptmin,ptmax,100,0,1);
2879 (Form(
"hE2ndEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2880 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()),
2881 nptbins,ptmin,ptmax,100,0,1);
2887 (Form(
"hE2ndEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2888 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()),
2889 nptbins,ptmin,ptmax,100,0,1);
2895 (Form(
"hE2ndEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2896 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()),
2897 nptbins,ptmin,ptmax,100,0,1);
2904 (Form(
"hECellClusRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2905 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()),
2906 nptbins,ptmin,ptmax,100,0,1);
2912 (Form(
"hECellClusRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2913 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()),
2914 nptbins,ptmin,ptmax,100,0,1);
2920 (Form(
"hECellClusRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2921 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()),
2922 nptbins,ptmin,ptmax,100,0,1);
2928 (Form(
"hLogECellNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2929 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()),
2930 nptbins,ptmin,ptmax,150,-3,3);
2936 (Form(
"hLogECellNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2937 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()),
2938 nptbins,ptmin,ptmax,150,-3,3);
2944 (Form(
"hLogECellNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2945 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()),
2946 nptbins,ptmin,ptmax,150,-3,3);
2952 (Form(
"hECellWeightNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2953 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()),
2954 nptbins,ptmin,ptmax,90,0,4.5);
2960 (Form(
"hECellWeightNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2961 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()),
2962 nptbins,ptmin,ptmax,90,0,4.5);
2968 (Form(
"hECellWeightNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2969 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()),
2970 nptbins,ptmin,ptmax,90,0,4.5);
2976 (Form(
"hMassEClusTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2977 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()),
2978 nptbins,ptmin,ptmax,nmassbins,massmin,massmax);
2992 (Form(
"hTimeDiffTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2993 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()),
2994 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
3000 (Form(
"hColRowTCardCorrelNCellLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3001 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3002 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3008 (Form(
"hColRowTCardCorrelNCellHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3009 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3010 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3018 (Form(
"hExoticTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3019 Form(
"exoticity vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3020 nptbins,ptmin,ptmax,200,-1,1);
3026 (Form(
"hTimeDiffExoTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3027 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()),
3028 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
3034 (Form(
"hColRowTCardCorrelNCellExoticLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3035 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3036 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3042 (Form(
"hColRowTCardCorrelNCellExoticHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3043 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3044 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3259 (Form(
"hLambda0Exoticity_EBin%d%s",i,add[tm].
Data()),
3260 Form(
"#lambda^{2}_{0} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3261 200,-1,1,ssbins,ssmin,ssmax);
3267 (Form(
"hLambda1Exoticity_EBin%d%s",i,add[tm].
Data()),
3268 Form(
"#lambda^{2}_{1} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3269 200,-1,1,ssbins,ssmin,ssmax);
3283 (Form(
"hNCellsExoticity_EBin%d%s",i,add[tm].
Data()),
3284 Form(
"#it{n}_{cells} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3285 200,-1,1,nceclbins,nceclmin,nceclmax);
3286 fhNCellsExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
3291 (Form(
"hTimeExoticity_EBin%d%s",i,add[tm].
Data()),
3293 200,-1,1,100,-25,25);
3294 fhTimeExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
3299 (Form(
"hNCellsTCardSameAndDiff_Exotic_EBin%d%s",i,add[tm].
Data()),
3300 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()),
3301 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
3341 (Form(
"hNCellsTCardSameAndDiff_EBin%d%s",i,add[tm].
Data()),
3342 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()),
3343 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
3349 (Form(
"hLambda0Lambda1_EBin%d%s",i,add[tm].
Data()),
3350 Form(
"#lambda^{2}_{0} vs #lambda^{2}_{1}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3351 ssbins,ssmin,ssmax,ssbins,ssmin,ssmax);
3367 for(
Int_t j = 0; j < 6; j++)
3369 for(
Int_t k = 0; k < 6; k++)
3372 (Form(
"hLambda0Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
3373 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()),
3374 200,-1,1,ssbins,ssmin,ssmax);
3380 (Form(
"hLambda1Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
3381 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()),
3382 200,-1,1,ssbins,ssmin,ssmax);
3479 for(
Int_t i = 0; i < 12; i++)
3482 (Form(
"hTCardCorrECellMaxDiff_Case%d%s",i,add[tm].
Data()),
3483 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()),
3484 nptbins,ptmin,ptmax,210,-1,20);
3490 (Form(
"hTCardCorrEClusterDiff_Case%d%s",i,add[tm].
Data()),
3491 Form(
"#it{E}_{cluster}-#it{E}_{cell} vs #it{E}_{cluster}, for (un)correlated cells in TCard, case %d %s",i,add[tm].
Data()),
3492 nptbins,ptmin,ptmax,210,-1,20);
3514 (Form(
"hTCardCorrTCellMaxDiff_Case%d%s",i,add[tm].
Data()),
3515 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()),
3516 nptbins,ptmin,ptmax,1000,-100,100);
3524 (Form(
"hTCardCorrECellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
3525 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()),
3526 nptbins,ptmin,ptmax,210,-1,20);
3532 (Form(
"hTCardCorrEClusterDiffExo_Case%d%s",i,add[tm].
Data()),
3533 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()),
3534 nptbins,ptmin,ptmax,210,-1,20);
3556 (Form(
"hTCardCorrTCellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
3557 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()),
3558 nptbins,ptmin,ptmax,1000,-100,100);
3569 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3575 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3581 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3587 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3593 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3599 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3605 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3611 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3619 (Form(
"hTMPhiResidual_EBin%d",i),
3620 Form(
"#Delta #varphi_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3621 200,-1,1,nresphibins,resphimin,resphimax);
3627 (Form(
"hTMEtaResidual_EBin%d",i),
3628 Form(
"#Delta #eta_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3629 200,-1,1,nresetabins,resetamin,resetamax);
3635 (Form(
"hTMPhiResidual_LooseCut_EBin%d",i),
3636 Form(
"#Delta #varphi_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3637 200,-1,1,nresphibins,resphimin,resphimax);
3643 (Form(
"hTMEtaResidual_LooseCut_EBin%d",i),
3644 Form(
"#Delta #eta_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3645 200,-1,1,nresetabins,resetamin,resetamax);
3673 (Form(
"hColRowM02"),
3675 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax,40,0.,2.);
3682 (Form(
"hColRowM02NCellCut"),
3683 Form(
"column vs row vs M02, %2.2f<#it{E}<%2.2f GeV #it{n}_{cells}^{w>0.01} > %d for Neutral",
3685 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax,40,0.,2.);
3693 Form(
"%2.2f<#it{E}_{1}<%2.2f GeV, %2.2f<#it{E}_{2}<%2.2f GeV, %2.2f<#lambda^{2}_{0}<%2.2f"
3694 "#it{M}_{#gamma #gamma} vs #it{n}_{1, cells}^{w>0.01} vs SM number trig cluster",
3703 (
"hInvMassNCellSMSame",
3704 Form(
"%2.2f<#it{E}_{1}<%2.2f GeV, %2.2f<#it{E}_{2}<%2.2f GeV, %2.2f<#lambda^{2}_{0}<%2.2f"
3705 "#it{M}_{#gamma #gamma} vs #it{n}_{1, cells}^{w>0.01} vs SM number both cluster",
3716 (
"hEMaxCellTimeM02SM",
3717 Form(
"time vs #lambda_{0}^{2} vs SM number, "
3726 (
"hEMaxCellTimeNCellSM",
3727 Form(
"time vs #it{n}_{cells}^{w>0.01} vs SM number, "
3736 (
"hESecCellTimeNCellSM",
3737 Form(
"secondary cell time vs #it{n}_{cells}^{w>0.01} vs SM number, "
3747 Int_t nbinsdeltacells = 19 ;
3750 for(
Int_t col = 0; col < 2; col++)
3753 (Form(
"hColRowFromCellMaxLowM02_Col%d",col),
3754 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, "
3756 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells);
3762 (Form(
"hColRowFromCellMaxHighM02_Col%d",col),
3763 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, "
3765 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells);
3777 (Form(
"hNCellsPerClusterM02M20_SM%d",i),
3778 Form(
" vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01} vs #lambda_{1}^{2}, "
3780 nShShBins/2,minShSh,maxShSh/2, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3787 (Form(
"hEMaxCellEClusterM02NCell_SM%d",i),
3788 Form(
"(#it{E}_{cluster} - #it{E}_{max cell})/#it{E}_{cluster} vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01}, "
3790 50,0,1., cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3797 (Form(
"hEMaxCellLogM02NCell_SM%d",i),
3798 Form(
"log(#it{E}_{max cell}) vs #lambda_{0}^{2} vs SM number, "
3800 150,-3,3, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3807 (Form(
"hESecCellEMaxCellM02NCell_SM%d",i),
3808 Form(
"#it{E}_{cell}/#it{E}_{cell max} vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01}, "
3810 50,0,1., cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3817 (Form(
"hESecCellEClusterM02NCell_SM%d",i),
3818 Form(
"(#it{E}_{cluster} - #it{E}_{cell})/#it{E}_{cluster} vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01}, "
3820 50,0,1., cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3827 (Form(
"hESecCellLogM02NCell_SM%d",i),
3828 Form(
"log(#it{E}_{cell}) vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01}, "
3830 150,-3,3, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3837 (Form(
"hEMaxESecCellNCellLowM02_SM%d",i),
3838 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, "
3840 24,0,12, 48,0,12,cellBins,cellMin,cellMax);
3847 (Form(
"hEMaxECrossNCellLowM02_SM%d",i),
3848 Form(
"#it{E}_{cell}^{max} vs exoticity vs #it{n}_{cells}^{w>0.01} vs 0.1 < #lambda_{0}^{2} < 0.3, "
3850 24,0,12, 40,0.6,1.,cellBins,cellMin,cellMax);
3857 (Form(
"hEMaxESecCellNCellHighM02_SM%d",i),
3858 Form(
"#it{E}_{cell}^{max} vs #it{E}_{cell}^{secondary} vs #it{n}_{cells}^{w>0.01} vs 0.5 < #lambda_{0}^{2} < 2, "
3860 24,0,12, 48,0,12,cellBins,cellMin,cellMax);
3867 (Form(
"hEMaxECrossNCellHighM02_SM%d",i),
3868 Form(
"#it{E}_{cell}^{max} vs exoticity vs #it{n}_{cells}^{w>0.01} for 0.5 < #lambda_{0}^{2} < 2, "
3870 24,0,12, 40,0.6,1.,cellBins,cellMin,cellMax);
3876 Int_t nbinsdeltacells = 19 ;
3879 for(
Int_t col = 0; col < 2; col++)
3882 (Form(
"hColRowFromCellMaxLowM02_SM%d_Col%d",i,col),
3883 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, "
3885 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,cellBins,cellMin,cellMax);
3892 (Form(
"hColRowFromCellMaxHighM02_SM%d_Col%d",i,col),
3893 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, "
3895 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,cellBins,cellMin,cellMax);
3901 for(
Int_t j = 0; j < 3; j++)
3904 (Form(
"hColRowFromCellMaxEMaxSecDiffLowM02_SM%d_Col%d_NCellBin%d",i,col,j),
3905 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, "
3906 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
3907 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,120,0,12);
3914 (Form(
"hColRowFromCellMaxEMaxSecDiffHighM02_SM%d_Col%d_NCellBin%d",i,col,j),
3915 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, "
3916 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
3917 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,120,0,12);
3924 (Form(
"hColRowFromCellMaxEMaxSecDiffFracLowM02_SM%d_Col%d_NCellBin%d",i,col,j),
3925 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, "
3926 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
3927 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,40,0,1);
3934 (Form(
"hColRowFromCellMaxEMaxSecDiffFracHighM02_SM%d_Col%d_NCellBin%d",i,col,j),
3935 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, "
3936 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
3937 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,40,0,1);
3944 (Form(
"hColRowFromCellMaxECellClusterRatLowM02_SM%d_Col%d_NCellBin%d",i,col,j),
3945 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, "
3946 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
3947 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,50,0,1);
3954 (Form(
"hColRowFromCellMaxECellClusterRatHighM02_SM%d_Col%d_NCellBin%d",i,col,j),
3955 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, "
3956 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
3957 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,50,0,1);
3968 for(
Int_t imatch = 0; imatch < 3; imatch++)
3973 (Form(
"hDeltaIEtaDeltaIPhi_%s",matchCase[imatch].
Data()),
3974 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),
3975 nEbins,minE,maxE,cellBins,cellMin,cellMax,cellBins,cellMin,cellMax);
3982 (Form(
"hDeltaIEtaDeltaIPhiTot_%s",matchCase[imatch].
Data()),
3983 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),
3984 nEbins,minE,maxE,cellBins,cellMin,cellMax,cellBins,cellMin,cellMax);
3991 (Form(
"hDeltaIA_%s",matchCase[imatch].
Data()),
3992 Form(
"Cluster *asymmetry* in cell units vs E for %s",matchCase[imatch].
Data()),
3993 nptbins,ptmin,ptmax,asyBins,asyMin,asyMax);
3994 fhDeltaIA[imatch]->SetXTitle(
"#it{E}_{cluster}");
3995 fhDeltaIA[imatch]->SetYTitle(
"#it{A}_{cell in cluster}");
3996 outputContainer->Add(
fhDeltaIA[imatch]);
3999 (Form(
"hDeltaIATot_%s",matchCase[imatch].
Data()),
4000 Form(
"Cluster *total asymmetry* in cell units vs E for %s",matchCase[imatch].
Data()),
4001 nptbins,ptmin,ptmax,asyBins,asyMin,asyMax);
4003 fhDeltaIATot[imatch]->SetYTitle(
"#it{A}_{cell in cluster}^{total}");
4007 (Form(
"hDeltaIAM02_%s",matchCase[imatch].
Data()),
4008 Form(
"Cluster *asymmetry* in cell units vs #lambda^{2}_{0} for %s",matchCase[imatch].
Data()),
4009 nEbins,minE,maxE,nShShBins,minShSh,maxShSh,asyBins,asyMin,asyMax);
4012 fhDeltaIAM02[imatch]->SetZTitle(
"#it{A}_{cell in cluster}");
4016 (Form(
"hDeltaIATotM02_%s",matchCase[imatch].
Data()),
4017 Form(
"Cluster *total asymmetry* in cell units vs #lambda^{2}_{0} for %s",matchCase[imatch].
Data()),
4018 nEbins,minE,maxE,nShShBins,minShSh,maxShSh,asyBins,asyMin,asyMax);
4021 fhDeltaIATotM02[imatch]->SetZTitle(
"#it{A}_{cell in cluster}^{total}");
4025 (Form(
"hDeltaIAM20_%s",matchCase[imatch].
Data()),
4026 Form(
"Cluster *asymmetry* in cell units vs #lambda^{2}_{1} for %s",matchCase[imatch].
Data()),
4027 nEbins,minE,maxE,(
Int_t)nShShBins/1.5,minShSh,(
Int_t)maxShSh/1.5,asyBins,asyMin,asyMax);
4030 fhDeltaIAM20[imatch]->SetZTitle(
"#it{A}_{cell in cluster}");
4034 (Form(
"hDeltaIATotM20_%s",matchCase[imatch].
Data()),
4035 Form(
"Cluster *total asymmetry* in cell units vs #lambda^{2}_{1} for %s",matchCase[imatch].
Data()),
4036 nEbins,minE,maxE,(
Int_t)nShShBins/1.5,minShSh,(
Int_t)maxShSh/1.5,asyBins,asyMin,asyMax);
4039 fhDeltaIATotM20[imatch]->SetZTitle(
"#it{A}_{cell in cluster}^{total}");
4043 (Form(
"hDeltaIANCells_%s",matchCase[imatch].
Data()),
4044 Form(
"Cluster *asymmetry* in cell units vs N cells in cluster for %s",matchCase[imatch].
Data()),
4045 nEbins,minE,maxE,cellBins,cellMin,cellMax,asyBins,asyMin,asyMax);
4052 (Form(
"hDeltaIATotNCells_%s",matchCase[imatch].
Data()),
4053 Form(
"Cluster *total asymmetry* in cell units vs N cells in cluster for %s",matchCase[imatch].
Data()),
4054 nEbins,minE,maxE,cellBins,cellMin,cellMax,asyBins,asyMin,asyMax);
4063 (Form(
"hDeltaIAOrigin_%s",matchCase[imatch].
Data()),
4064 Form(
"Cluster *asymmetry* in cell units vs E vs primary ID, for ID %s",matchCase[imatch].
Data()),
4065 nEbins,minE,maxE,10,-0.5,9.5,asyBins,asyMin,asyMax);
4069 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4070 fhDeltaIAOrigin[imatch]->GetYaxis()->SetBinLabel(ilabel,mcParticleStringLabel[ilabel-1]);
4074 (Form(
"hDeltaIATotOrigin_%s",matchCase[imatch].
Data()),
4075 Form(
"Cluster *total asymmetry* in cell units vs E vs primary ID, for ID %s",matchCase[imatch].
Data()),
4076 nEbins,minE,maxE,10,-0.5,9.5,asyBins,asyMin,asyMax);
4080 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4081 fhDeltaIATotOrigin[imatch]->GetYaxis()->SetBinLabel(ilabel,mcParticleStringLabel[ilabel-1]);
4087 (Form(
"hClusterMaxCellDiffM02_%s",matchCase[imatch].
Data()),
4088 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()),
4089 nEbins,minE,maxE, 20,0,1.,nShShBins,minShSh,maxShSh);
4096 (Form(
"hClusterTimeEnergyM02_%s",matchCase[imatch].
Data()),
4097 Form(
"#it{E} vs TOF vs #lambda_{0}^{2} for ID %s",matchCase[imatch].
Data()),
4098 nEbins,minE,maxE,45,-25.5,20.5,nShShBins,minShSh,maxShSh);
4105 (Form(
"hNCellsPerClusterM02_%s",matchCase[imatch].
Data()),
4106 Form(
"#it{E} vs #it{n}_{cells} vs #lambda_{0}^{2} for ID %s",matchCase[imatch].
Data()),
4107 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4114 (Form(
"hNCellsPerClusterM20_%s",matchCase[imatch].
Data()),
4115 Form(
"#it{E} vs #it{n}_{cells} vs #lambda_{1}^{2} for ID %s",matchCase[imatch].
Data()),
4116 nEbins,minE,maxE,cellBins,cellMin,cellMax,(
Int_t)nShShBins/1.5,minShSh,(
Int_t)maxShSh/1.5);
4123 (Form(
"hSMM02_%s",matchCase[imatch].
Data()),
4124 Form(
"#it{E} vs SM number vs #lambda_{0}^{2}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4125 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,minShSh,maxShSh);
4126 fhSMM02[imatch]->SetXTitle(
"#it{E} (GeV)");
4127 fhSMM02[imatch]->SetYTitle(
"SM number");
4128 fhSMM02[imatch]->SetZTitle(
"#lambda_{0}^{2}");
4129 outputContainer->Add(
fhSMM02[imatch]);
4132 (Form(
"hSMM02NoCut_%s",matchCase[imatch].
Data()),
4133 Form(
"#it{E} vs SM number vs #lambda_{0}^{2} for ID %s",matchCase[imatch].
Data()),
4134 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,minShSh,maxShSh);
4141 (Form(
"hSMNCell_%s",matchCase[imatch].
Data()),
4142 Form(
"#it{E} vs SM number vs #it{n}_{cells}^{w>0.01} for ID %s",matchCase[imatch].
Data()),
4143 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,cellBins,cellMin,cellMax);
4144 fhSMNCell[imatch]->SetXTitle(
"#it{E} (GeV)");
4145 fhSMNCell[imatch]->SetYTitle(
"SM number");
4146 fhSMNCell[imatch]->SetZTitle(
"#it{n}_{cells}^{w>0.01}");
4147 outputContainer->Add(
fhSMNCell[imatch]);
4150 (Form(
"hSMNCellM02_%s",matchCase[imatch].
Data()),
4151 Form(
"SM number vs #it{n}_{cells}^{w>0.01} vs #lambda_{0}^{2}, "
4153 fNModules,-0.5,fNModules-0.5,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4156 fhSMNCellM02[imatch]->SetYTitle(
"#it{n}_{cells}^{w>0.01}");
4160 (Form(
"hColM02_%s",matchCase[imatch].
Data()),
4161 Form(
"#it{E} vs column number vs #lambda_{0}^{2}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4162 nEbins,minE,maxE,48,-0.5,47.5,nShShBins,minShSh,maxShSh);
4163 fhColM02[imatch]->SetXTitle(
"#it{E} (GeV)");
4164 fhColM02[imatch]->SetYTitle(
"column number");
4165 fhColM02[imatch]->SetZTitle(
"#lambda_{0}^{2}");
4166 outputContainer->Add(
fhColM02[imatch]);
4169 (Form(
"hRowM02_%s",matchCase[imatch].
Data()),
4170 Form(
"#it{E} vs row number vs #lambda_{0}^{2}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4171 nEbins,minE,maxE,24,-0.5,23.5,nShShBins,minShSh,maxShSh);
4172 fhRowM02[imatch]->SetXTitle(
"#it{E} (GeV)");
4173 fhRowM02[imatch]->SetYTitle(
"row number");
4174 fhRowM02[imatch]->SetZTitle(
"#lambda_{0}^{2}");
4175 outputContainer->Add(
fhRowM02[imatch]);
4180 (Form(
"hNCellsPerClusterMEta_%s",matchCase[imatch].
Data()),
4181 Form(
"#it{E} vs #it{n}_{cells} vs #sigma_{#eta}^{2} for ID %s",matchCase[imatch].
Data()),
4182 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4189 (Form(
"hNCellsPerClusterMPhi_%s",matchCase[imatch].
Data()),
4190 Form(
"#it{E} vs #it{n}_{cells} vs #sigma_{#varphi}^{2} for ID %s",matchCase[imatch].
Data()),
4191 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4198 (Form(
"hNCellsPerClusterMEtaPhi_%s",matchCase[imatch].
Data()),
4199 Form(
"#it{E} vs #it{n}_{cells} vs #sigma_{#eta#varphi}^{2} for ID %s",matchCase[imatch].
Data()),
4200 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,-1*maxShSh,maxShSh);
4207 (Form(
"hNCellsPerClusterMEtaPhiA_%s",matchCase[imatch].
Data()),
4208 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()),
4209 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,-1*maxShSh,maxShSh);
4212 fhNCellsPerClusterMEtaPhiA[imatch]->SetZTitle(
"(#sigma_{#varphi}^{2}-#sigma_{#eta}^{2})/(#sigma_{#varphi}^{2}+#sigma_{#eta}^{2})");
4219 (Form(
"hOrigin_%s",matchCase[imatch].
Data()),
4220 Form(
"#it{E} vs origin for ID %s",matchCase[imatch].
Data()),
4221 nEbins,minE,maxE, 10,-0.5,9.5);
4222 fhOriginE[imatch]->SetXTitle(
"#it{E} (GeV)");
4223 fhOriginE[imatch]->SetYTitle(
"particle");
4224 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4225 fhOriginE[imatch]->GetYaxis()->SetBinLabel(ilabel ,mcParticleStringLabel[ilabel-1]);
4226 outputContainer->Add(
fhOriginE[imatch]);
4229 (Form(
"hOriginM02_%s",matchCase[imatch].
Data()),
4230 Form(
"#it{E} vs origin vs #lambda_{0}^{2} for ID %s",matchCase[imatch].
Data()),
4231 nEbins,minE,maxE,10,-0.5,9.5,nShShBins,minShSh,maxShSh);
4234 fhOriginM02[imatch]->SetZTitle(
"#lambda_{0}^{2}");
4235 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4236 fhOriginM02[imatch]->GetYaxis()->SetBinLabel(ilabel,mcParticleStringLabel[ilabel-1]);
4251 (
"hClusterMaxCellCloseCellRatioM02",
"#it{E} vs #it{E}_{cell-i}/#it{E}_{cell max} vs #lambda_{0}^{2}",
4252 nEbins,minE,maxE, 20,0,1.,nShShBins,minShSh,maxShSh);
4270 (
"hClusterMaxCellECrossM02",
4271 "#it{E} vs exoticity vs #lambda_{0}^{2}",
4272 nEbins,minE,maxE, 40,0.6,1.,nShShBins,minShSh,maxShSh);
4283 nptbins,ptmin,ptmax, 100,0,1.);
4289 nptbins,ptmin,ptmax, 100,-10,0);
4295 nptbins,ptmin,ptmax, 100,0,1.);
4301 nptbins,ptmin,ptmax, 100,-10,0);
4307 nptbins*2,ptmin,ptmax*2, 100,0,1.);
4313 nptbins*2,ptmin,ptmax*2, 100,-10,0);
4323 if(imod < fFirstModule || imod >
fLastModule)
continue;
4326 Form(
"#cell energy / sum all energy vs all energy in Module %d",imod),
4327 nptbins*2,ptmin,ptmax*2, 100,0,1.);
4333 Form(
"Log(cell energy / sum all energy) vs all energy in Module %d",imod),
4334 nptbins*2,ptmin,ptmax*2, 100,-10,0);
4341 Float_t cellEmin [] = {0.05,0.1,0.15,0.2};
4342 Float_t cellTmin [] = {50.,100.,1000000.};
4344 for(
Int_t iw = 0; iw < 12; iw++)
4347 for(
Int_t iEmin = 0; iEmin < 4; iEmin++)
4349 for(
Int_t iTmin = 0; iTmin < 3; iTmin++)
4352 Form(
"#lambda^{2}_{0} vs E, w0=%1.2f, cell E>%2.2f MeV, |t|<%2.0f ns",
4353 w0, cellEmin[iEmin], cellTmin[iTmin]),
4354 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4369 Form(
"#lambda^{2}_{0} vs E, w0=%1.2f, cell E>%2.2f MeV, |t|<%2.0f ns, |#eta| < 0.15",
4370 w0, cellEmin[iEmin], cellTmin[iTmin]),
4371 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4380 TString mcnames[] = {
"Photon",
"Electron",
"Conversion",
"Pi0",
"Hadron"};
4381 for(
Int_t imc = 0; imc < 5; imc++)
4384 Form(
"shower shape, #lambda^{2}_{0} vs E, w0 = %1.1f, for MC %s",w0,mcnames[imc].
Data()),
4385 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4404 return outputContainer;
4465 AliInfo(Form(
"dEdX cuts init for MC: %d<dedx_ele<%d; %d<dedx_had<%d",
4473 if ( runNumber < 146861 )
4480 else if ( runNumber < 156000 )
4487 else if ( runNumber < 165000 )
4502 AliInfo(Form(
"dEdX cuts init for run %d: %d<dedx_ele<%d; %d<dedx_had<%d",
4516 Int_t nCellsPerCluster)
4534 if(ampMax < 0.01)
return kFALSE;
4550 printf(
"**** Print %s %s ****\n", GetName(), GetTitle() ) ;
4554 printf(
"Min n cells : %d\n" ,
fNCellMin) ;
4556 printf(
"Min M02 : %1.2f\n",
fM02Min) ;
4572 AliDebug(1,
"Start");
4588 AliFatal(Form(
"AliAnaClusterShapeCorrelStudies::MakeAnalysisFillHistograms() - Wrong calorimeter name <%s>, END",
GetCalorimeterString().
Data()));
4592 AliFatal(Form(
"AliAnaClusterShapeCorrelStudies::MakeAnalysisFillHistograms() - No CaloClusters or CaloCells available"));
4622 if(clus->GetNCells() < 3)
return ;
4624 for (
Int_t ipos = 0; ipos < clus->GetNCells(); ipos++)
4626 Int_t id = clus->GetCellsAbsId()[ipos];
4641 AliWarning(Form(
"Wrong calculated energy %f",energy));
4652 for (
Int_t ipos = 0; ipos < clus->GetNCells(); ipos++)
4654 Int_t id = clus->GetCellsAbsId()[ipos];
4669 Float_t l0org = clus->GetM02();
4670 Float_t l1org = clus->GetM20();
4671 Float_t dorg = clus->GetDispersion();
4675 if(
IsDataMC() && clus->GetNLabels() > 0)
4700 Float_t cellEmin [] = {0.05,0.1,0.15,0.2};
4701 Float_t cellTmin [] = {50.,100.,1000000.};
4703 for(
Int_t iw = 0; iw < 12; iw++)
4705 for(
Int_t iEmin = 0; iEmin < 4; iEmin++)
4707 for(
Int_t iTmin = 0; iTmin < 3; iTmin++)
4713 Float_t disp = 0, dEta = 0, dPhi = 0;
4714 Float_t sEta = 0, sPhi = 0, sEtaPhi = 0;
4718 cellEmin[iEmin], cellTmin[iTmin], bc,
4719 newEnergy, l0, l1, disp, dEta, dPhi,
4720 sEta, sPhi, sEtaPhi);
4741 clus->SetM02(l0org);
4742 clus->SetM20(l1org);
4743 clus->SetDispersion(dorg);
TH2F * fhNCellsTCardSameAndDiffFraction[2]
! Cluster fraction of NCells in same TCard as leading vs energy
virtual Bool_t IsFiducialCutOn() const
Float_t GetHistoPtMax() const
TH2F * fhECellClusRatNLM1TCardCorrNoSelection[2]
! Cluster E cell / E cluster for NLM=1 vs E for clusters selected for TCard correlation studies ...
TH2F * fhLambda1NLM2TCardCorrNoSelection[2]
! Cluster m20 vs E for clusters selected for TCard correlation studies, nlm=2
TH2F * fhEMaxCellClusterRatio
! e max cell / e cluster vs e cluster
TH2F * fhLambda0NLM1TCardCorrelNCell[6][6][2]
! Cluster m02 vs E, cluster contains 0 to more than 6 cells with w > 0.01 in same TCard or diff TCard...
TH2F * fhE2ndRatNLM1TCardCorrelNCell[6][6][2]
! Cluster E cell second max / E cluster for NLM=1 vs E, cluster contains 0 to more than 6 cells with ...
void ClusterShapeHistograms(AliVCluster *cluster, Int_t absIdMax, Double_t maxCellFraction, Float_t eCrossFrac, Float_t ampMax, Double_t tmax, Int_t matchedPID, Int_t mcIndex)
TH2F * fhColRowTCardCorrNoSelectionExoticLowE[2]
! col-row cluster cell max for those selected for TCard correlation studies, 5 < E < 8 GeV...
TH2F * fhColRowTCardCorrelNCellExoticHighE[6][6][2]
! Cluster max cell col vs row, E > 8 GeV, cluster contains 0 to more than 6 cells with w > 0...
TH3F * fhClusterMaxCellCloseCellRatioM02
! Ratio between max cell energy and cell energy of the same cluster
Float_t fM02Min
Minimum M02 on clusters.
TH2F * fhECellTotalRatio
! e cell / e total vs e total
Int_t GetHistoNClusterCellMin() const
TH2F * fhExoticTCardCorrNoSelection[2]
! exoticity per cluster vs cluster energy, clusters selected for TCard correlation studies ...
TH2F * fhTMPhiResidualExoticity[14]
! Cluster-track matching residual in phi vs exoticity
Float_t GetHistoPtMin() const
TH2F * fhLogECellNLM2TCardCorrNoSelection[2]
! Cluster Log E cell for NLM=2 vs E for clusters selected for TCard correlation studies ...
TH2F * fhTimeTCardCorrNoSelection[2]
! Cluster time vs E for clusters selected for TCard correlation studies
TH2F * fhColRowTCardCorrelNCellExoticLowE[6][6][2]
! Cluster max cell col vs row, E > 2 GeV, cluster contains 0 to more than 6 cells with w > 0...
TH2F * fhEnergyTimeTCardCorrNoSelectionExotic[2]
! cluster energy vs time, exo > 0.97, T-Card strict cuts
Int_t fNMaxColsFull
Number of EMCAL/PHOS columns full detector.
TH2F * fhTCardCorrEClusterDiffExo[12][2]
! Cluster energy - secondary cell energy in cluster vs cluster energy, different secondary cell selec...
TH3F * fhClusterMaxCellDiffM02[3]
! Difference between cluster energy and energy of cell with more energy, vs m02
Int_t GetHistoShowerShapeBins() const
TH2F * fhSameRowDiffColAndTCardCellsTimeDiffCellMaxE[2]
! Secondary cell energy difference vs leading cell energy, one in same TCard as cell max...
virtual Int_t GetCalorimeter() const
TH2F * fhSameRowDiffColAndTCardCellsTimeDiffCellMaxEExo[2]
! Secondary cell energy difference vs leading cell energy, one in same TCard as cell max...
virtual void AddToHistogramsName(TString add)
virtual AliVCaloCells * GetEMCALCells() const
TH3F * fhNCellsPerClusterMEtaPhi[3]
! N cells per cluster vs cluster energy vs shape in eta*phi direction
TH3F * fhDeltaIAOrigin[3]
! Cluster "asymmetry" vs E vs origin
TH2F * fhECellClusRatNLM2TCardCorrelNCell[6][6][2]
! Cluster E cell / E cluster for NLM=2 vs E, cluster contains 0 to more than 6 cells with w > 0...
TH2F ** fhECellTotalRatioMod
! e cell / e total vs e total, per SM
TH3F * fhColRowM02
! Cluster position in cell max col-row, 8 < E < 12 GeV, neutral only
virtual Bool_t IsDataMC() const
TH3F * fhColRowFromCellMaxEMaxSecDiffLowM02PerSM[20][2][3]
! secondary cell distance to cell max in col vs row vs energy difference to cell max, (2-3), (4-5) (>5) n cells, for clusters 8 < E < 12 GeV, 0.1 < M02 < 0.3, per SM number, per odd/pair column
Int_t fNEBinCuts
Number of energy bin cuts for fStudyTCardCorrelation.
TH2F * fhNLocMaxTCardCorrelNCell[6][6][2]
! Cluster nlocmax vs E, cluster contains 0 to more than 6 cells with w > 0.01 in same TCard or diff T...
void WeightHistograms(AliVCluster *clus, Int_t mcTag)
TH3F * fhColM02[3]
! main cell column vs m02, n cell > 4
Int_t fLastModule
Last EMCAL/PHOS module, set in CaloUtils or depending fidutial cuts.
TH2F * fhNCellsTCardCorrWithWeightNoSelection[2]
! Ncells per cluster vs cluster energy, select cells with w>0.01, clusters selected for TCard correla...
TH3F * fhNCellsPerClusterMPhi[3]
! N cells per cluster vs cluster energy vs shape in phi direction
TH2F * fhE2ndEMaxRatNLM3TCardCorrNoSelection[2]
! Cluster E cell second loc max / E Max for NLM>2 vs E for clusters selected for TCard correlation st...
AliEMCALRecoUtils * GetEMCALRecoUtils() const
TH3F * fhInvMassNCellSMSame
! Invariant mass vs first selected cluster N cell with 8 < E < 12 GeV and per SM of both cluster ...
TH3F * fhDeltaIEtaDeltaIPhiTot[3]
! Difference between max cell index and farthest cells left/up to right/down, eta vs phi vs E ...
TH2F * fhColRowHighEPosTime
! cluster col-row cluster cell max, E > 8 GeV, t > 5 ns exo < 0.97, n cell > 1
TH2F * fhTCardCorrECellMaxDiff[12][2]
! Cell max energy - secondary cell energy in cluster vs cluster energy, different secondary cell sele...
Float_t GetECross(Int_t absId, AliVCaloCells *cells, Int_t bc)
TH2F * fhColRowHighENegTime
! cluster col-row cluster cell max, E > 8 GeV, t < -5 ns exo < 0.97, n cell > 1
TH2F * fhLogECellNLM3TCardCorrelNCell[6][6][2]
! Cluster log E cell vs E cluster for NLM>2 vs E, cluster contains 0 to more than 6 cells with w > 0...
Float_t GetPhi(Float_t phi) const
Shift phi angle in case of negative value 360 degrees. Example TLorenzVector::Phi defined in -pi to p...
TH2F * fhE2ndEMaxRatNLM1TCardCorrNoSelection[2]
! Cluster E cell second loc max / E Max for NLM=1 vs E for clusters selected for TCard correlation st...
TH2F * fhLogECellSameNLM2TCardCorrNoSelection[2]
! Cluster Log E cell for NLM=2 vs E for clusters selected for TCard correlation studies, same Tcard as leading
Float_t fEBinCuts[15]
Energy bins cut for fStudyTCardCorrelation.
virtual AliVEvent * GetInputEvent() const
Float_t fInvMassMaxM02Cut
Maximum M02 shower shape cut value for clusters entering the invariant mass calculation.
TH3F * fhNCellsPerClusterMEtaPhiA[3]
! N cells per cluster vs cluster energy vs shape in (phi-eta)/(eta+phi)
TH2F * fhLogECellSameNLM1TCardCorrNoSelection[2]
! Cluster Log E cell for NLM=1 vs E for clusters selected for TCard correlation studies, same Tcard as leading
Int_t fNModules
Number of EMCAL/PHOS modules to use in analysis, set in CaloUtils.
TH2F * fhLambda1TCardCorrNoSelection[2]
! Cluster m20 vs E for clusters selected for TCard correlation studies
TH2F * fhEMaxRatNLM3TCardCorrelNCell[6][6][2]
! Cluster E cell max / E cluster for NLM>2 vs E, cluster contains 0 to more than 6 cells with w > 0...
TH2F * fhEMaxRatNLM2TCardCorrelNCell[6][6][2]
! Cluster E cell max / E cluster for NLM=2 vs E, cluster contains 0 to more than 6 cells with w > 0...
TH2F * fhE2ndSameEMaxRatNLM2TCardCorrNoSelection[2]
! Cluster E cell second loc max / E Max for NLM=2 vs E for clusters selected for TCard correlation st...
TH2F * fhColRowExoticHighE1CellNegTime
! cluster col-row cluster cell max, E > 8 GeV, t < -5 ns
TH3F * fhEMaxCellTimeM02SM
! m02 vs SM number vs max cell time
Bool_t fStudyShapeParam
Study not only M02 but other kind of shape param.
TH2F * fhEnergyTMEtaResidualTCardCorrNoSelectionExotic
! cluster energy vs eta track-cluster residual, exo > 0.97, T-Card strict cuts