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++)
149 for(
Int_t i =0; i < 12; i++)
151 for(
Int_t j = 0; j < 4; j++)
153 for(
Int_t k = 0; k < 3; k++)
160 for(
Int_t j = 0; j < 5; j++)
168 for(
Int_t tm = 0; tm < 2; tm++)
264 for(
Int_t i = 0; i < 6; i++)
266 for(
Int_t j = 0; j < 6; j++)
339 for(
Int_t i = 0; i < 14; i++)
358 for(
Int_t j = 0; j < 6; j++)
360 for(
Int_t k = 0; k < 6; k++)
381 for(
Int_t i = 0; i < 12; i++)
407 for(
Int_t i = 0; i < 20; i++)
436 for(
Int_t j = 0; j < 3; j++)
480 Int_t icol = -1, irow = -1, iRCU = -1, icolAbs = -1, irowAbs = -1;
481 GetModuleNumberCellIndexesAbsCaloMap(absIdMax,GetCalorimeter(), icol, irow, iRCU, icolAbs, irowAbs);
485 for(
Int_t ie = 0; ie < fNEBinCuts; ie++)
487 if( energy >= fEBinCuts[ie] && energy < fEBinCuts[ie+1] ) ebin = ie;
490 Int_t ncells = clus->GetNCells();
491 Double_t time = clus->GetTOF()*1.e9;
492 if(time>400) time-=fConstantTimeShift;
494 Float_t deta = clus->GetTrackDz();
495 Float_t dphi = clus->GetTrackDx();
503 if(energy >= 5 && energy < 8)
504 fhColRowExoticLowE1Cell [matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
507 fhColRowExoticHighE1Cell[matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
509 if ( time > 5) fhColRowExoticHighE1CellPosTime->Fill(icolAbs,irowAbs,GetEventWeight());
510 else if( time < -5) fhColRowExoticHighE1CellNegTime->Fill(icolAbs,irowAbs,GetEventWeight());
511 else fhColRowExoticHighE1CellNulTime->Fill(icolAbs,irowAbs,GetEventWeight());
514 fhEnergyTime1Cell[matched]->Fill(energy,time,GetEventWeight());
515 fhEnergyTMEtaResidual1Cell->Fill(energy,deta,GetEventWeight());
516 fhEnergyTMPhiResidual1Cell->Fill(energy,dphi,GetEventWeight());
518 else if(exoticity > 0.97)
520 if(energy >= 5 && energy < 8)
521 fhColRowExoticLowE [matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
524 fhColRowExoticHighE[matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
526 if ( time > 5) fhColRowExoticHighEPosTime->Fill(icolAbs,irowAbs,GetEventWeight());
527 else if( time < -5) fhColRowExoticHighENegTime->Fill(icolAbs,irowAbs,GetEventWeight());
528 else fhColRowExoticHighENulTime->Fill(icolAbs,irowAbs,GetEventWeight());
531 fhEnergyTimeExotic[matched]->Fill(energy,time,GetEventWeight());
532 fhEnergyTMEtaResidualExotic->Fill(energy,deta,GetEventWeight());
533 fhEnergyTMPhiResidualExotic->Fill(energy,dphi,GetEventWeight());
535 for (
Int_t ipos = 0; ipos < ncells; ipos++)
537 Int_t absId = clus->GetCellsAbsId()[ipos];
539 Float_t eCell = fCaloCellList->GetCellAmplitude(absId);
540 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
543 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
545 if( absId == absIdMax || weight < 0.01 )
continue;
547 Int_t rowDiff = -100, colDiff = -100;
548 Bool_t sameTCard = GetCaloUtils()->IsAbsIDsFromTCard(absIdMax,absId,rowDiff,colDiff);
551 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
552 GetModuleNumberCellIndexesAbsCaloMap(absId,GetCalorimeter(), icol2, irow2, iRCU2, icolAbs2, irowAbs2);
556 if(energy >= 5 && energy < 8)
557 fhColRowExotic2ndCellDiffLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
559 fhColRowExotic2ndCellDiffHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
563 if(energy >= 5 && energy < 8)
564 fhColRowExotic2ndCellSameLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
566 fhColRowExotic2ndCellSameHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
571 else if ( energy > 8 )
573 if ( time > 5) fhColRowHighEPosTime->Fill(icolAbs,irowAbs,GetEventWeight());
574 else if( time < -5) fhColRowHighENegTime->Fill(icolAbs,irowAbs,GetEventWeight());
575 else fhColRowHighENulTime->Fill(icolAbs,irowAbs,GetEventWeight());
580 fhTMPhiResidualExoticityLooseCut[ebin]->Fill(exoticity,dphi);
581 fhTMEtaResidualExoticityLooseCut[ebin]->Fill(exoticity,deta);
584 else if ( energy > 8 && ncells > 1 && exoticity < 0.97)
586 if ( time > 5) fhColRowHighEPosTime->Fill(icolAbs,irowAbs,GetEventWeight());
587 else if( time < -5) fhColRowHighENegTime->Fill(icolAbs,irowAbs,GetEventWeight());
588 else fhColRowHighENulTime->Fill(icolAbs,irowAbs,GetEventWeight());
594 Int_t etaRegion = -1, phiRegion = -1;
595 GetCaloUtils()->GetEMCALSubregion(clus,fCaloCellList,etaRegion,phiRegion);
597 if ( etaRegion !=0 ) return ;
603 fhEnergyTimeTCardCorrNoSelection1Cell[matched]->Fill(energy,time,GetEventWeight());
604 fhEnergyTMEtaResidualTCardCorrNoSelection1Cell->Fill(energy,deta,GetEventWeight());
605 fhEnergyTMPhiResidualTCardCorrNoSelection1Cell->Fill(energy,dphi,GetEventWeight());
607 else if(exoticity > 0.97)
609 fhEnergyTimeTCardCorrNoSelectionExotic[matched]->Fill(energy,time,GetEventWeight());
610 fhEnergyTMEtaResidualTCardCorrNoSelectionExotic->Fill(energy,deta,GetEventWeight());
611 fhEnergyTMPhiResidualTCardCorrNoSelectionExotic->Fill(energy,dphi,GetEventWeight());
618 if ( m02 > 0.001 ) lamR = m20/m02;
620 Int_t absIdList[ncells];
622 Int_t nlm = GetCaloUtils()->GetNumberOfLocalMaxima(clus, fCaloCellList, absIdList, maxEList) ;
628 Int_t nCellWithWeight = 1;
636 Float_t eCellMax = fCaloCellList->GetCellAmplitude(absIdMax);
637 Double_t tCellMax = fCaloCellList->GetCellTime(absIdMax);
639 GetCaloUtils()->RecalibrateCellAmplitude(eCellMax, GetCalorimeter(), absIdMax);
640 GetCaloUtils()->RecalibrateCellTime(tCellMax, GetCalorimeter(), absIdMax, GetReader()->GetInputEvent()->GetBunchCrossNumber());
644 if(tCellMax>400) tCellMax-=fConstantTimeShift;
655 Bool_t sameTCard2ndMax = kFALSE;
659 Float_t eCellSameRowSameTCardNearCol = 0.;
660 Float_t eCellSameRowDiffTCardNearCol = 0.;
661 Double_t tCellSameRowSameTCardNearCol = 0.;
662 Double_t tCellSameRowDiffTCardNearCol = 0.;
669 for (
Int_t ipos = 0; ipos < ncells; ipos++)
671 Int_t absId = clus->GetCellsAbsId()[ipos];
673 Float_t eCell = fCaloCellList->GetCellAmplitude(absId);
674 Double_t tCell = fCaloCellList->GetCellTime(absId);
676 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
677 GetCaloUtils()->RecalibrateCellTime(tCell, GetCalorimeter(), absId, GetReader()->GetInputEvent()->GetBunchCrossNumber());
679 if(tCell>400) tCell-=fConstantTimeShift;
682 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
684 if( absId == absIdMax || weight < 0.01 )
continue;
688 fhECellClusRatNLM1TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
689 fhLogECellNLM1TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
690 fhECellWeightNLM1TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
694 fhECellClusRatNLM2TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
695 fhLogECellNLM2TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
696 fhECellWeightNLM2TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
700 fhECellClusRatNLM3TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
701 fhLogECellNLM3TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
702 fhECellWeightNLM3TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
705 Int_t rowDiff = -100, colDiff = -100;
706 Bool_t sameTCard = GetCaloUtils()->IsAbsIDsFromTCard(absIdMax,absId,rowDiff,colDiff);
712 fhECellSameClusRatNLM1TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
713 fhLogECellSameNLM1TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
714 fhECellSameWeightNLM1TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
718 fhECellSameClusRatNLM2TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
719 fhLogECellSameNLM2TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
720 fhECellSameWeightNLM2TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
724 fhECellSameClusRatNLM3TCardCorrNoSelection[matched]->Fill(energy, eCell/energy, GetEventWeight());
725 fhLogECellSameNLM3TCardCorrNoSelection [matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
726 fhECellSameWeightNLM3TCardCorrNoSelection [matched]->Fill(energy, weight, GetEventWeight());
736 Int_t indexType = -1;
741 if(TMath::Abs(rowDiff) == 1) nearRow = kTRUE;
742 if(TMath::Abs(colDiff) == 1) nearCol = kTRUE;
744 if ( rowDiff == 0 && colDiff != 0 )
746 if ( nearCol ) indexType = 6;
752 else if ( rowDiff != 0 && colDiff == 0 )
754 if ( nearRow ) indexType = 8;
762 if ( nearRow && nearCol) indexType = 10;
773 if ( rowDiff == 0 && colDiff != 0 )
775 if ( nearCol ) indexType = 0;
779 else if ( rowDiff != 0 && colDiff == 0 )
781 if ( nearRow ) indexType = 2;
786 if ( nearCol && nearRow ) indexType = 4;
791 if ( rowDiff == 0 && TMath::Abs(colDiff) == 1 )
795 eCellSameRowSameTCardNearCol = eCell;
796 tCellSameRowSameTCardNearCol = tCell;
800 eCellSameRowDiffTCardNearCol = eCell;
801 tCellSameRowDiffTCardNearCol = tCell;
807 Float_t eCellDiff = eCellMax - eCell;
808 Float_t eClusDiff = energy - eCell;
811 Float_t tCellDiff = tCellMax - tCell;
813 fhTCardCorrECellMaxDiff[indexType][matched]->Fill(energy, eCellDiff, GetEventWeight());
814 fhTCardCorrEClusterDiff[indexType][matched]->Fill(energy, eClusDiff, GetEventWeight());
817 fhTCardCorrTCellMaxDiff[indexType][matched]->Fill(energy, tCellDiff, GetEventWeight());
819 if ( fStudyExotic && exoticity > 0.97 )
821 fhTCardCorrECellMaxDiffExo[indexType][matched]->Fill(energy, eCellDiff, GetEventWeight());
822 fhTCardCorrEClusterDiffExo[indexType][matched]->Fill(energy, eClusDiff, GetEventWeight());
825 fhTCardCorrTCellMaxDiffExo[indexType][matched]->Fill(energy, tCellDiff, GetEventWeight());
829 if ( fStudyExotic && exoticity > 0.97 )
832 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
833 GetModuleNumberCellIndexesAbsCaloMap(absId,GetCalorimeter(), icol2, irow2, iRCU2, icolAbs2, irowAbs2);
837 if(energy >= 5 && energy < 8)
838 fhColRowTCardCorrNoSelectionExotic2ndCellDiffLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
840 fhColRowTCardCorrNoSelectionExotic2ndCellDiffHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
844 if(energy >= 5 && energy < 8)
845 fhColRowTCardCorrNoSelectionExotic2ndCellSameLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
847 fhColRowTCardCorrNoSelectionExotic2ndCellSameHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
851 if ( eCell > emax2nd )
857 sameTCard2ndMax = kTRUE;
863 sameTCard2ndMax = kFALSE;
898 Float_t ratioNcells = nCellWithWeight/(ncells*1.);
899 fhNCellsTCardCorrRatioWithWeightNoSelection[matched]->Fill(energy, ratioNcells, GetEventWeight());
903 if ( nCellWithWeight <= 1 )
return;
908 printf(
"AliAnaClusterShapeCorrelStudies: M02 %f, M20 %f, E %2.3f, ncell %d, n with weight %d; max cell E %2.3f\n",
909 m02,m20,energy,ncells,nCellWithWeight,eCellMax);
918 if(energy >= 5 && energy < 8)
919 fhColRowTCardCorrNoSelectionLowE [matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
921 fhColRowTCardCorrNoSelectionHighE[matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
923 if ( fStudyExotic && exoticity > 0.97 )
925 if(energy >= 5 && energy < 8)
926 fhColRowTCardCorrNoSelectionExoticLowE [matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
928 fhColRowTCardCorrNoSelectionExoticHighE[matched]->Fill(icolAbs,irowAbs,GetEventWeight()) ;
931 fhTimeTCardCorrNoSelection [matched]->Fill(energy,time, GetEventWeight());
932 fhLambda0TCardCorrNoSelection[matched]->Fill(energy, m02, GetEventWeight());
933 fhLambda1TCardCorrNoSelection[matched]->Fill(energy, m20, GetEventWeight());
937 fhLambda0NLM1TCardCorrNoSelection[matched]->Fill(energy, m02, GetEventWeight());
938 fhLambda1NLM1TCardCorrNoSelection[matched]->Fill(energy, m20, GetEventWeight());
942 fhLambda0NLM2TCardCorrNoSelection[matched]->Fill(energy, m02, GetEventWeight());
943 fhLambda1NLM2TCardCorrNoSelection[matched]->Fill(energy, m20, GetEventWeight());
946 fhLambdaRTCardCorrNoSelection[matched]->Fill(energy,lamR, GetEventWeight());
947 fhNLocMaxTCardCorrNoSelection[matched]->Fill(energy, nlm, GetEventWeight());
948 fhExoticTCardCorrNoSelection [matched]->Fill(energy, exoticity, GetEventWeight());
952 fhEMaxRatNLM1TCardCorrNoSelection [matched]->Fill(energy, eCellMax/energy , GetEventWeight());
953 fhE2ndRatNLM1TCardCorrNoSelection [matched]->Fill(energy, emax2nd/energy , GetEventWeight());
954 fhE2ndEMaxRatNLM1TCardCorrNoSelection[matched]->Fill(energy, emax2nd/eCellMax, GetEventWeight());
958 fhEMaxRatNLM2TCardCorrNoSelection [matched]->Fill(energy, eCellMax/energy , GetEventWeight());
959 fhE2ndRatNLM2TCardCorrNoSelection [matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
960 fhE2ndEMaxRatNLM2TCardCorrNoSelection[matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
964 fhEMaxRatNLM3TCardCorrNoSelection [matched]->Fill(energy, eCellMax/energy , GetEventWeight());
965 fhE2ndRatNLM3TCardCorrNoSelection [matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
966 fhE2ndEMaxRatNLM3TCardCorrNoSelection[matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
973 fhE2ndSameRatNLM1TCardCorrNoSelection [matched]->Fill(energy, emax2nd/energy , GetEventWeight());
974 fhE2ndSameEMaxRatNLM1TCardCorrNoSelection[matched]->Fill(energy, emax2nd/eCellMax, GetEventWeight());
978 fhE2ndSameRatNLM2TCardCorrNoSelection [matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
979 fhE2ndSameEMaxRatNLM2TCardCorrNoSelection[matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
983 fhE2ndSameRatNLM3TCardCorrNoSelection [matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
984 fhE2ndSameEMaxRatNLM3TCardCorrNoSelection[matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
988 fhNCellsTCardCorrNoSelection [matched]->Fill(energy, ncells, GetEventWeight());
989 fhNCellsTCardCorrWithWeightNoSelection[matched]->Fill(energy, nCellWithWeight, GetEventWeight());
991 if(eCellSameRowSameTCardNearCol > 0 && eCellSameRowDiffTCardNearCol > 0)
993 Float_t eDiff = eCellSameRowSameTCardNearCol - eCellSameRowDiffTCardNearCol ;
994 Float_t tDiff = tCellSameRowSameTCardNearCol - tCellSameRowDiffTCardNearCol ;
996 fhSameRowDiffColAndTCardCellsEnergyDiffClusterE[matched]->Fill(energy , eDiff, GetEventWeight());
997 fhSameRowDiffColAndTCardCellsTimeDiffClusterE [matched]->Fill(energy , tDiff, GetEventWeight());
998 fhSameRowDiffColAndTCardCellsEnergyDiffCellMaxE[matched]->Fill(eCellMax, eDiff, GetEventWeight());
999 fhSameRowDiffColAndTCardCellsTimeDiffCellMaxE [matched]->Fill(eCellMax, tDiff, GetEventWeight());
1001 if ( fStudyExotic && exoticity > 0.97 )
1003 fhSameRowDiffColAndTCardCellsEnergyDiffClusterEExo[matched]->Fill(energy , eDiff, GetEventWeight());
1004 fhSameRowDiffColAndTCardCellsTimeDiffClusterEExo [matched]->Fill(energy , tDiff, GetEventWeight());
1005 fhSameRowDiffColAndTCardCellsEnergyDiffCellMaxEExo[matched]->Fill(eCellMax, eDiff, GetEventWeight());
1006 fhSameRowDiffColAndTCardCellsTimeDiffCellMaxEExo [matched]->Fill(eCellMax, tDiff, GetEventWeight());
1011 Int_t nCorrInd = nCorr;
1012 if(nCorr > 4) nCorrInd = 5;
1014 Int_t nCorrNoInd = nCorrNo;
1015 if(nCorrNoInd > 4) nCorrNoInd = 5;
1028 fhLambda0NLM1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m02, GetEventWeight());
1029 fhLambda1NLM1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m20, GetEventWeight());
1033 fhLambda0NLM2TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m02, GetEventWeight());
1034 fhLambda1NLM2TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m20, GetEventWeight());
1037 fhLambda0TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m02, GetEventWeight());
1038 fhLambda1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, m20, GetEventWeight());
1040 fhNLocMaxTCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, nlm, GetEventWeight());
1043 fhExoticTCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, exoticity, GetEventWeight());
1047 fhEMaxRatNLM1TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, eCellMax/energy , GetEventWeight());
1048 fhE2ndRatNLM1TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, emax2nd/energy , GetEventWeight());
1049 fhE2ndEMaxRatNLM1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, emax2nd/eCellMax, GetEventWeight());
1053 fhEMaxRatNLM2TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, eCellMax/energy, GetEventWeight());
1054 fhE2ndRatNLM2TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
1055 fhE2ndEMaxRatNLM2TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
1059 fhEMaxRatNLM3TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, eCellMax/energy, GetEventWeight());
1060 fhE2ndRatNLM3TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, maxEList[1]/energy , GetEventWeight());
1061 fhE2ndEMaxRatNLM3TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, maxEList[1]/eCellMax, GetEventWeight());
1065 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1067 Int_t absId = clus->GetCellsAbsId()[ipos];
1069 Float_t eCell = fCaloCellList->GetCellAmplitude(absId);
1070 Double_t tCell = fCaloCellList->GetCellTime(absId);
1072 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
1073 GetCaloUtils()->RecalibrateCellTime(tCell, GetCalorimeter(), absId, GetReader()->GetInputEvent()->GetBunchCrossNumber());
1075 if(tCell>400) tCell-=fConstantTimeShift;
1078 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
1080 if( absId == absIdMax || weight < 0.01 )
continue;
1082 Float_t tDiffMaxSecondary = tCellMax - tCell;
1084 fhTimeDiffTCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, tDiffMaxSecondary, GetEventWeight());
1085 if ( fStudyExotic && exoticity > 0.97 )
1087 fhTimeDiffExoTCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, tDiffMaxSecondary, GetEventWeight());
1089 Int_t icol2 = -1, irow2 = -1, iRCU2 = -1, icolAbs2 = -1, irowAbs2 = -1;
1090 GetModuleNumberCellIndexesAbsCaloMap(absId,GetCalorimeter(), icol2, irow2, iRCU2, icolAbs2, irowAbs2);
1092 Int_t rowDiff = -100, colDiff = -100;
1093 Bool_t sameTCard = GetCaloUtils()->IsAbsIDsFromTCard(absIdMax,absId,rowDiff,colDiff);
1097 if ( energy >= 5 && energy < 8 && nCorr == 0 )
1098 fhColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
1099 else if ( energy >= 8 && nCorr == 0 )
1100 fhColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
1104 if ( energy >= 5 && energy < 8 && nCorrNo == 0 )
1105 fhColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffLowE [matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
1106 else if ( energy >= 8 && nCorrNo == 0 )
1107 fhColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffHighE[matched]->Fill(icolAbs2,irowAbs2,GetEventWeight()) ;
1113 fhECellClusRatNLM1TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, eCell/energy, GetEventWeight());
1114 fhLogECellNLM1TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
1115 fhECellWeightNLM1TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, weight, GetEventWeight());
1119 fhECellClusRatNLM2TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, eCell/energy, GetEventWeight());
1120 fhLogECellNLM2TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
1121 fhECellWeightNLM2TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, weight, GetEventWeight());
1125 fhECellClusRatNLM3TCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, eCell/energy, GetEventWeight());
1126 fhLogECellNLM3TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, TMath::Log(eCell), GetEventWeight());
1127 fhECellWeightNLM3TCardCorrelNCell [nCorrInd][nCorrNoInd][matched]->Fill(energy, weight, GetEventWeight());
1132 if(m02 > fInvMassMinM02Cut && m02 < fInvMassMaxM02Cut)
1134 for(
Int_t jclus = 0 ; jclus < fCaloClusList->GetEntriesFast() ; jclus++)
1136 AliVCluster* clus2 = (AliVCluster*) fCaloClusList->At(jclus);
1139 Int_t absIdMax2 = GetCaloUtils()->GetMaxEnergyCell(fCaloCellList, clus2, maxCellFraction);
1141 Double_t tof2 = clus2->GetTOF()*1.e9;
1142 if(tof2>400) tof2-=fConstantTimeShift;
1149 if( absIdMax == absIdMax2
1150 || !IsGoodCluster(absIdMax2, clus2->GetM02(), clus2->GetNCells())
1151 || GetCaloUtils()->GetNumberOfLocalMaxima(clus2,fCaloCellList) > 1
1152 || clus2->GetM02() > fInvMassMaxM02Cut
1153 || clus2->GetM02() < fInvMassMinM02Cut
1154 || clus2->E() < fInvMassMinECut
1155 || clus2->E() > fInvMassMaxECut
1156 || TMath::Abs(diffTof) > fInvMassMaxTimeDifference
1161 clus2->GetMomentum(fClusterMomentum2,v);
1165 if(IsFiducialCutOn()) in2 = GetFiducialCut()->IsInFiducialCut(fClusterMomentum2.Eta(),fClusterMomentum2.Phi(),GetCalorimeter()) ;
1171 Double_t angle = fClusterMomentum.Angle(fClusterMomentum2.Vect());
1173 if( angle > fInvMassMaxOpenAngle )
continue;
1176 Float_t mass = (fClusterMomentum+fClusterMomentum2).M ();
1177 fhMassEClusTCardCorrelNCell[nCorrInd][nCorrNoInd][matched]->Fill(energy, mass, GetEventWeight());
1182 if ( energy >= 5 && energy < 8)
1183 fhColRowTCardCorrelNCellLowE [nCorrInd][nCorrNoInd][matched]->Fill(icolAbs, irowAbs, GetEventWeight());
1184 else if ( energy >= 8 )
1185 fhColRowTCardCorrelNCellHighE[nCorrInd][nCorrNoInd][matched]->Fill(icolAbs, irowAbs, GetEventWeight());
1200 if ( fStudyExotic && exoticity > 0.97 )
1202 if ( energy >= 5 && energy < 8)
1203 fhColRowTCardCorrelNCellExoticLowE [nCorrInd][nCorrNoInd][matched]->Fill(icolAbs, irowAbs, GetEventWeight());
1204 else if ( energy >= 8 )
1205 fhColRowTCardCorrelNCellExoticHighE[nCorrInd][nCorrNoInd][matched]->Fill(icolAbs, irowAbs, GetEventWeight());
1242 fhLambda0Lambda1 [ebin][matched]->Fill(m20, m02, GetEventWeight());
1243 fhNCellsTCardSameAndDiff[ebin][matched]->Fill(nCorrNo, nCorr, GetEventWeight());
1249 fhLambda0Exoticity[ebin][matched]->Fill(exoticity, m02, GetEventWeight());
1250 fhLambda1Exoticity[ebin][matched]->Fill(exoticity, m20, GetEventWeight());
1252 fhNCellsExoticity [ebin][matched]->Fill(exoticity, nCellWithWeight, GetEventWeight());
1253 fhTimeExoticity [ebin][matched]->Fill(exoticity, tCellMax, GetEventWeight());
1257 fhLambda0ExoticityPerNCell[nCorrInd][nCorrNoInd][matched]->Fill(exoticity, m02, GetEventWeight());
1258 fhLambda1ExoticityPerNCell[nCorrInd][nCorrNoInd][matched]->Fill(exoticity, m20, GetEventWeight());
1270 if ( exoticity > 0.97 )
1271 fhNCellsTCardSameAndDiffExotic[ebin][matched]->Fill(nCorrNo, nCorr, GetEventWeight());
1274 fhTMPhiResidualExoticity[ebin]->Fill(exoticity,dphi);
1275 fhTMEtaResidualExoticity[ebin]->Fill(exoticity,deta);
1285 Float_t nCellRat = nCorr*1. / ((nCorr+nCorrNo)*1.);
1286 fhNCellsTCardSameAndDiffFraction[matched]->Fill(energy, nCellRat, GetEventWeight());
1287 if ( fStudyExotic && exoticity > 0.97 )
1288 fhNCellsTCardSameAndDiffFractionExotic[matched]->Fill(energy, nCellRat, GetEventWeight());
1387 Int_t nCaloCellsPerCluster = clus->GetNCells();
1388 if ( nCaloCellsPerCluster < 2 )
return;
1393 Int_t nlm = GetCaloUtils()->GetNumberOfLocalMaxima(clus, fCaloCellList) ;
1405 for (
Int_t ipos = 0; ipos < nCaloCellsPerCluster; ipos++)
1407 Int_t absId = clus ->GetCellsAbsId()[ipos];
1408 Float_t eCell = fCaloCellList->GetCellAmplitude(absId) ;
1410 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
1412 if( absId == absIdMax || eCell < 0.01 )
continue;
1414 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
1416 if( weight < 0.01 )
continue;
1421 if ( nCell < 1 )
return;
1426 Int_t ietaMax=-1, iphiMax = 0, rcuMax = 0, icolAbs = -1, irowAbs = -1;
1427 Int_t smMax = GetModuleNumberCellIndexesAbsCaloMap(absIdMax, GetCalorimeter(),
1428 ietaMax, iphiMax, rcuMax, icolAbs, irowAbs);
1430 if ( matchedPID == 0 && energy > fEMinShape && energy < fEMaxShape )
1432 fhColRowM02->Fill(icolAbs,irowAbs,m02,GetEventWeight()) ;
1434 if ( nCell > fNCellMinShape )
1435 fhColRowM02NCellCut->Fill(icolAbs,irowAbs,m02,GetEventWeight()) ;
1443 Int_t etaRegion = -1, phiRegion = -1;
1444 GetCaloUtils()->GetEMCALSubregion(clus,fCaloCellList,etaRegion,phiRegion);
1446 if ( etaRegion !=0 ) return ;
1450 for (
Int_t ipos = 0; ipos < nCaloCellsPerCluster; ipos++)
1452 Int_t absId = clus ->GetCellsAbsId()[ipos];
1453 Float_t eCell = fCaloCellList->GetCellAmplitude(absId) ;
1455 GetCaloUtils()->RecalibrateCellAmplitude(eCell, GetCalorimeter(), absId);
1457 if( absId == absIdMax || eCell < 0.01 )
continue;
1459 Float_t weight = GetCaloUtils()->GetEMCALRecoUtils()->GetCellWeight(eCell, energy);
1461 if( weight < 0.01 )
continue;
1463 Float_t fracCell = eCell/eCellMax;
1464 fhClusterMaxCellCloseCellRatioM02->Fill(energy, fracCell, m02, GetEventWeight());
1466 Float_t fracClus = (energy-eCell)/energy;
1468 Float_t logECell = TMath::Log(eCell);
1479 if ( energy > fEMinShape && energy < fEMaxShape && matchedPID == 0 )
1481 fhESecCellEMaxCellM02NCellPerSM[smMax]->Fill(fracCell, nCell, m02, GetEventWeight());
1482 fhESecCellEClusterM02NCellPerSM[smMax]->Fill(fracClus, nCell, m02, GetEventWeight());
1483 fhESecCellLogM02NCellPerSM [smMax]->Fill(logECell, nCell, m02, GetEventWeight());
1485 if ( m02 > 0.1 && m02 < 0.3 )
1487 fhEMaxESecCellNCellLowM02PerSM [smMax]->Fill(eCellMax, eCell , nCell, GetEventWeight());
1488 fhEMaxECrossNCellLowM02PerSM [smMax]->Fill(eCellMax, eCrossFrac, nCell, GetEventWeight());
1490 Double_t time = fCaloCellList->GetCellTime(absId);
1491 GetCaloUtils()->RecalibrateCellTime(time, GetCalorimeter(), absId,GetReader()->GetInputEvent()->GetBunchCrossNumber());
1493 if(time>400) time-=fConstantTimeShift;
1495 fhESecCellTimeNCellSM->Fill(time, smMax, nCell, GetEventWeight());
1497 else if ( m02 > 0.5 && m02 < 2 )
1499 fhEMaxESecCellNCellHighM02PerSM[smMax]->Fill(eCellMax, eCell , nCell, GetEventWeight());
1500 fhEMaxECrossNCellHighM02PerSM [smMax]->Fill(eCellMax, eCrossFrac, nCell, GetEventWeight());
1510 Int_t sm = GetModuleNumberCellIndexes(absId,GetCalorimeter(), ieta, iphi, rcu);
1512 if(dIphi < TMath::Abs(iphi-iphiMax)) dIphi = TMath::Abs(iphi-iphiMax);
1513 if(iphi-iphiMax < 0 && dIphiNeg > iphi-iphiMax) dIphiNeg = iphi-iphiMax;
1514 if(iphi-iphiMax > 0 && dIphiPos < iphi-iphiMax) dIphiPos = iphi-iphiMax;
1518 if(dIeta < TMath::Abs(ieta-ietaMax)) dIeta = TMath::Abs(ieta-ietaMax);
1519 if(ieta-ietaMax < 0 && dIetaNeg > ieta-ietaMax) dIetaNeg = ieta-ietaMax;
1520 if(ieta-ietaMax > 0 && dIetaPos < ieta-ietaMax) dIetaPos = ieta-ietaMax;
1522 if ( energy > fEMinShape && energy < fEMaxShape && matchedPID == 0 )
1525 if ( nCell == 2 || nCell == 3 ) nCellBin = 0;
1526 else if ( nCell == 4 || nCell == 5 ) nCellBin = 1;
1527 if ( m02 > 0.1 && m02 < 0.3 )
1529 fhColRowFromCellMaxLowM02PerSM[smMax][ietaMax%2]->Fill(ietaMax-ieta, iphiMax-iphi, nCell, GetEventWeight());
1530 fhColRowFromCellMaxLowM02 [ietaMax%2]->Fill(ietaMax-ieta, iphiMax-iphi, GetEventWeight());
1532 fhColRowFromCellMaxEMaxSecDiffFracLowM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi,(eCellMax-eCell)/eCellMax, GetEventWeight());
1533 fhColRowFromCellMaxEMaxSecDiffLowM02PerSM [smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi, eCellMax-eCell, GetEventWeight());
1534 fhColRowFromCellMaxECellClusterRatLowM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi, eCell/energy , GetEventWeight());
1536 else if (m02 > 0.5 && m02 < 2)
1538 fhColRowFromCellMaxHighM02PerSM[smMax][ietaMax%2]->Fill(ietaMax-ieta, iphiMax-iphi, nCell, GetEventWeight());
1539 fhColRowFromCellMaxHighM02 [ietaMax%2]->Fill(ietaMax-ieta, iphiMax-iphi, GetEventWeight());
1541 fhColRowFromCellMaxEMaxSecDiffFracHighM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi,(eCellMax-eCell)/eCellMax, GetEventWeight());
1542 fhColRowFromCellMaxEMaxSecDiffHighM02PerSM [smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi, eCellMax-eCell, GetEventWeight());
1543 fhColRowFromCellMaxECellClusterRatHighM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(ietaMax-ieta, iphiMax-iphi, eCell/energy , GetEventWeight());
1550 Int_t ietaShift = ieta;
1551 Int_t ietaMaxShift = ietaMax;
1553 if (ieta > ietaMax) ietaMaxShift+=48;
1554 else ietaShift +=48;
1556 if(dIeta < TMath::Abs(ietaShift-ietaMaxShift)) dIeta = TMath::Abs(ietaShift-ietaMaxShift);
1557 if(ietaShift-ietaMaxShift < 0 && dIetaNeg > ietaShift-ietaMaxShift) dIetaNeg = ietaShift-ietaMaxShift;
1558 if(ietaShift-ietaMaxShift > 0 && dIetaPos < ietaShift-ietaMaxShift) dIetaPos = ietaShift-ietaMaxShift;
1564 if ( energy > fEMinShape && energy < fEMaxShape && matchedPID == 0 )
1566 fhNCellsPerClusterM02M20PerSM[smMax]->Fill(m20, nCell, m02, GetEventWeight());
1567 fhNCellsPerClusterM02NLMPerSM[smMax]->Fill(nlm, nCell, m02, GetEventWeight());
1571 if ( nCell == 2 || nCell == 3 ) nCellBin = 0;
1572 else if ( nCell == 4 || nCell == 5 ) nCellBin = 1;
1573 if ( m02 > 0.1 && m02 < 0.3 )
1574 fhColRowFromCellMaxECellClusterRatLowM02PerSM [smMax][ietaMax%2][nCellBin]->Fill(0., 0., eCellMax/energy, GetEventWeight());
1575 else if (m02 > 0.5 && m02 < 2)
1576 fhColRowFromCellMaxECellClusterRatHighM02PerSM[smMax][ietaMax%2][nCellBin]->Fill(0., 0., eCellMax/energy, GetEventWeight());
1580 fhClusterMaxCellECrossM02->Fill(energy, eCrossFrac, m02, GetEventWeight());
1585 fhClusterMaxCellDiffM02[matchedPID]->Fill(energy, maxFrac, m02, GetEventWeight());
1586 fhClusterTimeEnergyM02 [matchedPID]->Fill(energy, tmax , m02, GetEventWeight());
1587 fhNCellsPerClusterM02 [matchedPID]->Fill(energy, nCell , m02, GetEventWeight());
1588 fhNCellsPerClusterM20 [matchedPID]->Fill(energy, nCell , m20, GetEventWeight());
1589 fhNCellsPerClusterNLM [matchedPID]->Fill(energy, nCell , nlm, GetEventWeight());
1591 if ( matchedPID == 0 )
1593 fhNCellsPerClusterM02PerSMPerCol[smMax][ietaMax%2]->Fill(energy, nCell, m02, GetEventWeight());
1594 fhNCellsPerClusterM20PerSMPerCol[smMax][ietaMax%2]->Fill(energy, nCell, m20, GetEventWeight());
1595 fhNCellsPerClusterNLMPerSMPerCol[smMax][ietaMax%2]->Fill(energy, nCell, nlm, GetEventWeight());
1598 fhSMNCell [matchedPID]->Fill(energy, smMax , nCell, GetEventWeight());
1599 fhSMM02NoCut [matchedPID]->Fill(energy, smMax , m02 , GetEventWeight());
1600 if ( m02 > 0.1 && m02 < 0.3 )
1601 fhSMM20LowM02NoCut [matchedPID]->Fill(energy, smMax , m20 , GetEventWeight());
1602 else if ( m02 > 0.5 && m02 < 2 )
1603 fhSMM20HighM02NoCut[matchedPID]->Fill(energy, smMax , m20 , GetEventWeight());
1606 if ( energy > fEMinShape && energy < fEMaxShape )
1608 fhSMNCellM02[matchedPID]->Fill(smMax , nCell, m02, GetEventWeight());
1610 if ( matchedPID == 0 )
1612 if ( m02 > 0.1 && m02 < 0.3 )
1613 fhEMaxCellTimeNCellSM->Fill(tmax, smMax, nCell, GetEventWeight());
1615 fhEMaxCellTimeM02SM->Fill(tmax, smMax, m02, GetEventWeight());
1617 fhEMaxCellEClusterM02NCellPerSM[smMax]->Fill(maxFrac, nCell, m02, GetEventWeight());
1618 fhEMaxCellLogM02NCellPerSM [smMax]->Fill(TMath::Log(eCellMax), nCell, m02, GetEventWeight());
1623 if ( nCell > fNCellMinShape )
1625 if ( m02 > 0.1 && m02 < 0.3 )
1626 fhSMM20LowM02 [matchedPID]->Fill(energy, smMax , m20 , GetEventWeight());
1627 else if ( m02 > 0.5 && m02 < 2 )
1628 fhSMM20HighM02[matchedPID]->Fill(energy, smMax , m20 , GetEventWeight());
1630 fhSMM02 [matchedPID]->Fill(energy, smMax , m02, GetEventWeight());
1631 fhColM02[matchedPID]->Fill(energy, ietaMax, m02, GetEventWeight());
1632 fhRowM02[matchedPID]->Fill(energy, iphiMax, m02, GetEventWeight());
1638 if ( fStudyShapeParam )
1641 Float_t dIA = 1.*(dIphi-dIeta)/(dIeta+dIphi);
1642 Float_t dIATot = 1.*((dIphiPos-dIphiNeg)-(dIetaPos-dIetaNeg))/((dIetaPos-dIetaNeg)+(dIphiPos-dIphiNeg));
1651 AliDebug(1,Form(
"E %2.2f, nCell %d, dPhi %d, dEta %d, dIA %2.2f, match %d",energy,nCell, dIphi,dIeta,dIA,matchedPID));
1661 fhDeltaIANCells [matchedPID]->Fill(energy, nCell , dIA , GetEventWeight());
1662 fhDeltaIATotNCells [matchedPID]->Fill(energy, nCell , dIATot, GetEventWeight());
1664 if ( nCell > fNCellMinShape )
1666 fhDeltaIEtaDeltaIPhi[matchedPID]->Fill(energy, dIeta, dIphi, GetEventWeight());
1667 fhDeltaIA [matchedPID]->Fill(energy, dIA , GetEventWeight());
1668 fhDeltaIAM02 [matchedPID]->Fill(energy, m02 , dIA , GetEventWeight());
1669 fhDeltaIAM20 [matchedPID]->Fill(energy, m20 , dIA , GetEventWeight());
1671 fhDeltaIEtaDeltaIPhiTot[matchedPID]->Fill(energy,dIetaPos-dIetaNeg, dIphiPos-dIphiNeg, GetEventWeight());
1672 fhDeltaIATot [matchedPID]->Fill(energy, dIATot , GetEventWeight());
1673 fhDeltaIATotM02 [matchedPID]->Fill(energy, m02 , dIATot, GetEventWeight());
1674 fhDeltaIATotM20 [matchedPID]->Fill(energy, m20 , dIATot, GetEventWeight());
1677 if ( IsDataMC() && mcIndex > -1 && mcIndex < 10 && nCell > fNCellMinShape )
1679 fhDeltaIAOrigin [matchedPID]->Fill(energy, mcIndex, dIA , GetEventWeight());
1680 fhDeltaIATotOrigin[matchedPID]->Fill(energy, mcIndex, dIATot, GetEventWeight());
1684 Float_t dispp= 0., dEta = 0., dPhi = 0.;
1685 Float_t sEta = 0., sPhi = 0., sEtaPhi = 0.;
1686 if ( GetCalorimeter() ==
kEMCAL )
1688 GetCaloUtils()->GetEMCALRecoUtils()->RecalculateClusterShowerShapeParameters(GetEMCALGeometry(), fCaloCellList, clus,
1689 l0, l1, dispp, dEta, dPhi, sEta, sPhi, sEtaPhi);
1692 if(sEta+sPhi>0.0001) sEtaPhiA = (sPhi-sEta)/(sEta+sPhi);
1694 AliDebug(2,Form(
"Recalculate shower shape org: m02 %2.2f, m20 %2.2f, disp %2.2f;"
1695 " new: m02 %2.2f, m20 %2.2f, disp %2.2f; "
1696 "mEta %2.2f, mPhi %2.2f, mEtaPhi %2.2f, A_EtaPhi %2.2f; dEta %2.2f dPhi %2.2f",
1697 m02,m20,clus->GetDispersion(),l0,l1,dispp,sEta,sPhi,sEtaPhi,sEtaPhiA,dEta,dPhi));
1699 fhNCellsPerClusterMEta [matchedPID]->Fill(energy, nCell, sEta , GetEventWeight());
1700 fhNCellsPerClusterMPhi [matchedPID]->Fill(energy, nCell, sPhi , GetEventWeight());
1701 fhNCellsPerClusterMEtaPhi [matchedPID]->Fill(energy, nCell, sEtaPhi , GetEventWeight());
1702 fhNCellsPerClusterMEtaPhiA[matchedPID]->Fill(energy, nCell, sEtaPhiA, GetEventWeight());
1704 fhSMMEtaNoCut [matchedPID]->Fill(energy, smMax, sEta , GetEventWeight());
1705 fhSMMPhiNoCut [matchedPID]->Fill(energy, smMax, sPhi , GetEventWeight());
1706 fhSMMEtaPhiNoCut [matchedPID]->Fill(energy, smMax, sEtaPhi , GetEventWeight());
1707 fhSMMEtaPhiANoCut[matchedPID]->Fill(energy, smMax, sEtaPhiA, GetEventWeight());
1710 if ( nCell > fNCellMinShape )
1712 fhSMMEta [matchedPID]->Fill(energy, smMax, sEta , GetEventWeight());
1713 fhSMMPhi [matchedPID]->Fill(energy, smMax, sPhi , GetEventWeight());
1714 fhSMMEtaPhi [matchedPID]->Fill(energy, smMax, sEtaPhi , GetEventWeight());
1715 fhSMMEtaPhiA[matchedPID]->Fill(energy, smMax, sEtaPhiA, GetEventWeight());
1723 if ( matchedPID == 0 && energy > fEMinShape && energy < fEMaxShape &&
1724 m02 > fInvMassMinM02Cut && m02 < fInvMassMaxM02Cut )
1726 for(
Int_t jclus = 0 ; jclus < fCaloClusList->GetEntriesFast() ; jclus++)
1728 AliVCluster* clus2 = (AliVCluster*) fCaloClusList->At(jclus);
1731 Int_t absIdMax2 = GetCaloUtils()->GetMaxEnergyCell(fCaloCellList, clus2, maxCellFraction);
1733 Double_t tof2 = clus2->GetTOF()*1.e9;
1734 if(tof2>400) tof2-=fConstantTimeShift;
1741 if( absIdMax == absIdMax2
1742 || !IsGoodCluster(absIdMax2, clus2->GetM02(), clus2->GetNCells())
1743 || GetCaloUtils()->GetNumberOfLocalMaxima(clus2,fCaloCellList) > 1
1744 || clus2->GetM02() > fInvMassMaxM02Cut
1745 || clus2->GetM02() < fInvMassMinM02Cut
1746 || clus2->E() < fInvMassMinECut
1747 || clus2->E() > fInvMassMaxECut
1748 || TMath::Abs(diffTof) > fInvMassMaxTimeDifference
1753 clus2->GetMomentum(fClusterMomentum2,v);
1757 if(IsFiducialCutOn()) in2 = GetFiducialCut()->IsInFiducialCut(fClusterMomentum2.Eta(),fClusterMomentum2.Phi(),GetCalorimeter()) ;
1763 Double_t angle = fClusterMomentum.Angle(fClusterMomentum2.Vect());
1765 if( angle > fInvMassMaxOpenAngle )
continue;
1768 Float_t mass = (fClusterMomentum+fClusterMomentum2).M ();
1769 fhInvMassNCellSM->Fill(mass, nCell, smMax, GetEventWeight());
1771 Int_t smMax2 = GetModuleNumber(clus2);
1772 if(smMax == smMax2) fhInvMassNCellSMSame->Fill(mass, nCell, smMax, GetEventWeight());
1779 if ( IsDataMC() && mcIndex > -1 && mcIndex < 10)
1781 fhOriginE [matchedPID]->Fill(energy, mcIndex, GetEventWeight());
1782 fhOriginM02[matchedPID]->Fill(energy, mcIndex, m02, GetEventWeight());
1798 AliVTrack *track = GetCaloUtils()->GetMatchedTrack(clus, GetReader()->GetInputEvent());
1803 AliDebug(1,Form(
"No track found with matches %d",clus->GetNTracksMatched()));
1809 Double_t dedx = track->GetTPCsignal();
1824 if(fdEdXMinEle == 0 || fdEdXMaxEle == 0 || fdEdXMinHad == 0 || fdEdXMaxHad == 0)
1825 InitdEdXParameters();
1827 if ( dedx >= fdEdXMinEle && dedx < fdEdXMaxEle ) matchedPID = 1;
1828 else if ( dedx >= fdEdXMinHad && dedx < fdEdXMaxHad ) matchedPID = 2;
1829 else matchedPID =-1;
1831 AliDebug(1,Form(
"matches %d; dEdX %2.2f; matched PID %d",clus->GetNTracksMatched(),dedx,matchedPID));
1847 Int_t nCaloCellsPerCluster = 0 ;
1858 for(
Int_t iclus = 0; iclus < nCaloClusters; iclus++)
1860 AliDebug(1,Form(
"Cluster: %d/%d, data %d",iclus+1,nCaloClusters,
GetReader()->GetDataType()));
1862 AliVCluster* clus = (AliVCluster*)
fCaloClusList->At(iclus);
1867 AliDebug(1,Form(
"Small distance to bad channel %2.2f < %2.2f",clus->GetDistanceToBadChannel(),
fMinDistToBad));
1873 if ( nModule < fFirstModule || nModule >
fLastModule )
1875 AliDebug(1,Form(
"Cluster module out of range %d",nModule));
1885 Double_t tof = clus->GetTOF()*1.e9;
1903 AliDebug(1,Form(
"Remove cluster with phi %2.2f and eta %2.2f", phi*TMath::RadToDeg(), eta));
1907 AliDebug(1,Form(
"cluster: E %2.3f, pT %2.3f, eta %2.3f, phi %2.3f",e,pt,eta,phi*TMath::RadToDeg()));
1911 nCaloCellsPerCluster = clus->GetNCells();
1915 AliDebug(1,Form(
"Accept cluster? %d",goodCluster));
1917 if(!goodCluster)
continue;
1946 Int_t matchedPID = 0;
1963 if ( ampMax > 0.01 )
1971 if (
fStudyShape && matchedPID >= 0 && matchedPID < 3 )
1992 const Int_t buffersize = 255;
1993 char onePar[buffersize] ;
1995 snprintf(onePar,buffersize,
"--- AliAnaClusterShapeCorrelStudies ---:") ;
2016 return new TObjString(parList) ;
2026 outputContainer->SetName(
"ClusterShapeStudies") ;
2037 TString mcParticleStringLabel[] = {
"Merged #gamma#gamma",
"Decay #gamma",
"Direct #gamma",
"e^{#pm}",
"#pi^{#pm}",
"k^{#pm}",
"p",
"#bar{p}",
"n",
"#bar{n}"};
2041 TString matchCase[] = {
"Neutral",
"Electron",
"Hadron"};
2069 Float_t rowcellmin = fNMaxRowsFullMin-1.5;
2078 Int_t nShShBins = 200;
2088 Int_t cellBins = 15;
2098 (
"hColRowHighEPosTime",
2099 "column vs row, exo < 0.97, E > 8 GeV, t > 5 ns",
2100 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2106 (
"hColRowHighENegTime",
2107 "column vs row, exo < 0.97, E > 8 GeV, t < -5 ns",
2108 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2114 (
"hColRowHighENulTime",
2115 "column vs row, exo < 0.97, E > 8 GeV, -5 < t < 5 ns",
2116 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2124 (
"hColRowExoticHighE1CellPosTime",
2125 "column vs row, 1 cell, E > 8 GeV, t > 5 ns",
2126 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2132 (
"hColRowExoticHighEPosTime",
2133 "column vs row, exo > 0.97, E > 8 GeV, t > 5 ns",
2134 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2140 (
"hColRowExoticHighE1CellNegTime",
2141 "column vs row, 1 cell, E > 8 GeV, t < -5 ns",
2142 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2148 (
"hColRowExoticHighENegTime",
2149 "column vs row, exo > 0.97, E > 8 GeV, t < -5 ns",
2150 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2156 (
"hColRowExoticHighE1CellNulTime",
2157 "column vs row, 1 cell, E > 8 GeV, -5 < t < 5 ns",
2158 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2164 (
"hColRowExoticHighENulTime",
2165 "column vs row, exo > 0.97, E > 8 GeV, -5 < t < 5 ns",
2166 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2172 TString add[] = {
"",
"TrackMatched"};
2173 for(
Int_t tm = 0; tm < 2; tm++)
2176 (Form(
"hColRowTCardCorrNoSelectionLowE%s",add[tm].
Data()),
2177 Form(
"column vs row, max E cell for TCard correlation selected clusters, 5 < E < 8 GeV %s",add[tm].
Data()),
2178 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2184 (Form(
"hColRowTCardCorrNoSelectionHighE%s",add[tm].
Data()),
2185 Form(
"column vs row, max E cell for TCard correlation selected clusters, E > 8 GeV %s",add[tm].
Data()),
2186 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2194 (Form(
"hNCellsTCardCorrNoSelection%s",add[tm].
Data()),
2195 Form(
"# custer # cells vs #it{E} %s",add[tm].
Data()),
2196 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
2202 (Form(
"hNCellsTCardCorrWithWeightNoSelection%s",add[tm].
Data()),
2203 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
2204 nptbins,ptmin,ptmax, nceclbins,nceclmin,nceclmax);
2210 (Form(
"hNCellsTCardCorrRatioWithWeightNoSelection%s",add[tm].
Data()),
2211 Form(
"custer # cells vs #it{E}, w > 0.01 %s",add[tm].
Data()),
2212 nptbins,ptmin,ptmax, 100,0,1);
2218 (Form(
"hTimeTCardCorrNoSelection%s",add[tm].
Data()),
2219 Form(
"#it{time} vs #it{E} %s",add[tm].
Data()),
2220 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2226 (Form(
"hLambda0TCardCorrNoSelection%s",add[tm].
Data()),
2227 Form(
"#lambda^{2}_{0} vs #it{E} %s",add[tm].
Data()),
2228 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2234 (Form(
"hLambda1TCardCorrNoSelection%s",add[tm].
Data()),
2235 Form(
"#lambda^{2}_{1} vs #it{E} %s",add[tm].
Data()),
2236 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2242 (Form(
"hLambda0NLM1TCardCorrNoSelection%s",add[tm].
Data()),
2243 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=1 %s",add[tm].
Data()),
2244 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2250 (Form(
"hLambda1NLM1TCardCorrNoSelection%s",add[tm].
Data()),
2251 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=1 %s",add[tm].
Data()),
2252 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2258 (Form(
"hLambda0NLM2TCardCorrNoSelection%s",add[tm].
Data()),
2259 Form(
"#lambda^{2}_{0} vs #it{E}, nlm=2 %s",add[tm].
Data()),
2260 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2266 (Form(
"hLambda1NLM2TCardCorrNoSelection%s",add[tm].
Data()),
2267 Form(
"#lambda^{2}_{1} vs #it{E}, nlm=2 %s",add[tm].
Data()),
2268 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2275 (Form(
"hLambdaRTCardCorrNoSelection%s",add[tm].
Data()),
2276 Form(
"#lambda^{1}_{0}/#lambda^{2}_{0} vs #it{E} %s",add[tm].
Data()),
2277 nptbins,ptmin,ptmax,110,0,1.1);
2283 (Form(
"hNLocMaxTCardCorrNoSelection%s",add[tm].
Data()),
2284 Form(
"#it{n}_{LM} vs E %s",add[tm].
Data()),
2285 nptbins,ptmin,ptmax,10,0,10);
2291 (Form(
"hEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2292 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
2293 nptbins,ptmin,ptmax,100,0,1);
2299 (Form(
"hEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2300 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2301 nptbins,ptmin,ptmax,100,0,1);
2307 (Form(
"hEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2308 Form(
"#it{E}_{cell}^{max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2309 nptbins,ptmin,ptmax,100,0,1);
2315 (Form(
"hE2ndRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2316 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cluster} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
2317 nptbins,ptmin,ptmax,100,0,1);
2323 (Form(
"hE2ndRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2324 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2325 nptbins,ptmin,ptmax,100,0,1);
2331 (Form(
"hE2ndRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2332 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2333 nptbins,ptmin,ptmax,100,0,1);
2340 (Form(
"hE2ndEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2341 Form(
"#it{E}_{cell}^{2nd max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=1 %s",add[tm].
Data()),
2342 nptbins,ptmin,ptmax,100,0,1);
2348 (Form(
"hE2ndEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2349 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2350 nptbins,ptmin,ptmax,100,0,1);
2356 (Form(
"hE2ndEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2357 Form(
"#it{E}_{cell}^{2nd loc max}/#it{E}_{cell}^{max} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2358 nptbins,ptmin,ptmax,100,0,1);
2366 (Form(
"hE2ndSameRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2367 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()),
2368 nptbins,ptmin,ptmax,100,0,1);
2374 (Form(
"hE2ndSameRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2375 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()),
2376 nptbins,ptmin,ptmax,100,0,1);
2382 (Form(
"hE2ndSameRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2383 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()),
2384 nptbins,ptmin,ptmax,100,0,1);
2391 (Form(
"hE2ndSameEMaxRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2392 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()),
2393 nptbins,ptmin,ptmax,100,0,1);
2399 (Form(
"hE2ndSameEMaxRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2400 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()),
2401 nptbins,ptmin,ptmax,100,0,1);
2407 (Form(
"hE2ndSameEMaxRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2408 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()),
2409 nptbins,ptmin,ptmax,100,0,1);
2417 (Form(
"hECellClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2418 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
2419 nptbins,ptmin,ptmax,100,0,1);
2425 (Form(
"hECellClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2426 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2427 nptbins,ptmin,ptmax,100,0,1);
2433 (Form(
"hECellClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2434 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2435 nptbins,ptmin,ptmax,100,0,1);
2441 (Form(
"hLogECellNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2442 Form(
"log(#it{E}_{cell}) vs E_{cluster}, #it{n}_{LM}=1, w > 0.01 %s",add[tm].
Data()),
2443 nptbins,ptmin,ptmax,150,-3,3);
2449 (Form(
"hLogECellNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2450 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}=2, w > 0.01 %s",add[tm].
Data()),
2451 nptbins,ptmin,ptmax,150,-3,3);
2457 (Form(
"hLogECellNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2458 Form(
"log(#it{E}_{cell}) vs E, #it{n}_{LM}>2, w > 0.01 %s",add[tm].
Data()),
2459 nptbins,ptmin,ptmax,150,-3,3);
2466 (Form(
"hECellWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2467 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E_{cluster}, #it{n}_{LM}=1 %s",add[tm].
Data()),
2468 nptbins,ptmin,ptmax,90,0,4.5);
2474 (Form(
"hECellWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2475 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster})) vs E, #it{n}_{LM}=2 %s",add[tm].
Data()),
2476 nptbins,ptmin,ptmax,90,0,4.5);
2482 (Form(
"hECellWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2483 Form(
"#it{w}=Max(4,5+log(#it{E}_{cell}/#it{E}_{cluster}))vs E, #it{n}_{LM}>2 %s",add[tm].
Data()),
2484 nptbins,ptmin,ptmax,90,0,4.5);
2491 (Form(
"hECellSameClusRatNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2492 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()),
2493 nptbins,ptmin,ptmax,100,0,1);
2499 (Form(
"hECellSameClusRatNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2500 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, cell from same T-Card as leading, #it{n}_{LM}=2 %s",add[tm].
Data()),
2501 nptbins,ptmin,ptmax,100,0,1);
2507 (Form(
"hECellSameClusRatNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2508 Form(
"#it{E}_{cell}/#it{E}_{cluster} vs E, cell from same T-Card as leading, #it{n}_{LM}>2 %s",add[tm].
Data()),
2509 nptbins,ptmin,ptmax,100,0,1);
2515 (Form(
"hLogECellSameNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2516 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()),
2517 nptbins,ptmin,ptmax,150,-3,3);
2523 (Form(
"hLogECellSameNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2524 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()),
2525 nptbins,ptmin,ptmax,150,-3,3);
2531 (Form(
"hLogECellSameNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2532 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()),
2533 nptbins,ptmin,ptmax,150,-3,3);
2540 (Form(
"hECellSameWeightNLM1TCardCorrNoSelection%s",add[tm].
Data()),
2541 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()),
2542 nptbins,ptmin,ptmax,90,0,4.5);
2548 (Form(
"hECellSameWeightNLM2TCardCorrNoSelection%s",add[tm].
Data()),
2549 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()),
2550 nptbins,ptmin,ptmax,90,0,4.5);
2556 (Form(
"hECellSameWeightNLM3TCardCorrNoSelection%s",add[tm].
Data()),
2557 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()),
2558 nptbins,ptmin,ptmax,90,0,4.5);
2564 (Form(
"hExoticTCardCorrNoSelection%s",add[tm].
Data()),
2565 Form(
"exoticity vs #it{E} %s",add[tm].
Data()),
2566 nptbins,ptmin,ptmax,200,-1,1);
2572 (Form(
"hNCellsTCardSameAndDiffFraction%s",add[tm].
Data()),
2573 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, %s",add[tm].
Data()),
2574 nptbins,ptmin,ptmax,100,0,1);
2580 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterE%s",add[tm].
Data()),
2581 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()),
2582 nptbins,ptmin,ptmax,200,-10,10);
2588 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterE%s",add[tm].
Data()),
2589 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()),
2590 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2596 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxE%s",add[tm].
Data()),
2597 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()),
2598 nptbins,ptmin,ptmax,200,-10,10);
2604 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxE%s",add[tm].
Data()),
2605 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()),
2606 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2614 (Form(
"hEnergyTime1Cell%s",add[tm].
Data()),
2615 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
2616 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2622 (Form(
"hEnergyTimeExotic%s",add[tm].
Data()),
2623 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
2624 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2630 (Form(
"hEnergyTimeTCardCorrNoSelection1Cell%s",add[tm].
Data()),
2631 Form(
"#it{t} vs #it{E}, 1 cells cluster %s",add[tm].
Data()),
2632 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2638 (Form(
"hEnergyTimeTCardCorrNoSelectionExotic%s",add[tm].
Data()),
2639 Form(
"#it{t} vs #it{E}, exo > 0.97, %s",add[tm].
Data()),
2640 nptbins,ptmin,ptmax,ntimebins,timemin,timemax);
2646 (Form(
"hColRowExoticLowE1Cell%s",add[tm].
Data()),
2647 Form(
"column vs row, 1 cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2648 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2654 (Form(
"hColRowExoticHighE1Cell%s",add[tm].
Data()),
2655 Form(
"column vs row, 1 cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2656 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2662 (Form(
"hColRowExoticLowE%s",add[tm].
Data()),
2663 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2664 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2670 (Form(
"hColRowExoticHighE%s",add[tm].
Data()),
2671 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2672 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2678 (Form(
"hColRowExotic2ndCellDiffLowE%s",add[tm].
Data()),
2679 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2680 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2686 (Form(
"hColRowExotic2ndCellDiffHighE%s",add[tm].
Data()),
2687 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2688 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2694 (Form(
"hColRowExotic2ndCellSameLowE%s",add[tm].
Data()),
2695 Form(
"column vs row, max E cell, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2696 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2702 (Form(
"hColRowExotic2ndCellSameHighE%s",add[tm].
Data()),
2703 Form(
"column vs row, max E cell, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2704 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2710 (Form(
"hColRowTCardCorrNoSelectionExoticLowE%s",add[tm].
Data()),
2711 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2712 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2718 (Form(
"hColRowTCardCorrNoSelectionExoticHighE%s",add[tm].
Data()),
2719 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2720 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2726 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffLowE%s",add[tm].
Data()),
2727 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2728 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2734 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffHighE%s",add[tm].
Data()),
2735 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2736 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2742 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameLowE%s",add[tm].
Data()),
2743 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2744 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2750 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameHighE%s",add[tm].
Data()),
2751 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2752 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2758 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameLowE%s",add[tm].
Data()),
2759 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2760 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2766 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellDiffNoSameHighE%s",add[tm].
Data()),
2767 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2768 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2774 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffLowE%s",add[tm].
Data()),
2775 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, 5 < E < 8 GeV %s",add[tm].
Data()),
2776 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2782 (Form(
"hColRowTCardCorrNoSelectionExotic2ndCellSameNoDiffHighE%s",add[tm].
Data()),
2783 Form(
"column vs row, max E cell for TCard correlation selected clusters, exo > 0.97, E > 8 GeV %s",add[tm].
Data()),
2784 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
2790 (Form(
"hNCellsTCardSameAndDiffFraction_Exotic%s",add[tm].
Data()),
2791 Form(
"#it{n}_{cells} same TCard vs diff TCard fraction, w > 0.01, exo > 0.97 %s",add[tm].
Data()),
2792 nptbins,ptmin,ptmax,100,0,1);
2798 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffClusterEExo%s",add[tm].
Data()),
2799 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()),
2800 nptbins,ptmin,ptmax,200,-10,10);
2806 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffClusterEExo%s",add[tm].
Data()),
2807 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()),
2808 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2814 (Form(
"hSameRowDiffColAndTCardCellsEnergyDiffCellMaxEExo%s",add[tm].
Data()),
2815 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()),
2816 nptbins,ptmin,ptmax,200,-10,10);
2822 (Form(
"hSameRowDiffColAndTCardCellsTimeDiffCellMaxEExo%s",add[tm].
Data()),
2823 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()),
2824 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
2830 for(
Int_t i = 0; i < 6; i++)
2832 for(
Int_t j = 0; j < 6; j++)
2835 (Form(
"hLambda0TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2836 Form(
"#lambda^{2}_{0} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
2837 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2843 (Form(
"hLambda1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2844 Form(
"#lambda^{2}_{1} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
2845 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2851 (Form(
"hLambda0NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2852 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()),
2853 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2859 (Form(
"hLambda1NLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2860 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()),
2861 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2867 (Form(
"hLambda0NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2868 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()),
2869 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2875 (Form(
"hLambda1NLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2876 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()),
2877 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2892 (Form(
"hNLocMaxTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2893 Form(
"#it{n}_{LM} vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
2894 nptbins,ptmin,ptmax,10,0,10);
2900 (Form(
"hEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2901 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()),
2902 nptbins,ptmin,ptmax,100,0,1);
2908 (Form(
"hEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2909 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()),
2910 nptbins,ptmin,ptmax,100,0,1);
2916 (Form(
"hEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2917 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()),
2918 nptbins,ptmin,ptmax,100,0,1);
2924 (Form(
"hE2ndRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2925 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()),
2926 nptbins,ptmin,ptmax,100,0,1);
2932 (Form(
"hE2ndRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2933 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()),
2934 nptbins,ptmin,ptmax,100,0,1);
2940 (Form(
"hE2ndRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2941 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()),
2942 nptbins,ptmin,ptmax,100,0,1);
2948 (Form(
"hE2ndEMaxRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2949 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()),
2950 nptbins,ptmin,ptmax,100,0,1);
2956 (Form(
"hE2ndEMaxRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2957 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()),
2958 nptbins,ptmin,ptmax,100,0,1);
2964 (Form(
"hE2ndEMaxRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2965 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()),
2966 nptbins,ptmin,ptmax,100,0,1);
2973 (Form(
"hECellClusRatNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2974 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()),
2975 nptbins,ptmin,ptmax,100,0,1);
2981 (Form(
"hECellClusRatNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2982 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()),
2983 nptbins,ptmin,ptmax,100,0,1);
2989 (Form(
"hECellClusRatNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2990 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()),
2991 nptbins,ptmin,ptmax,100,0,1);
2997 (Form(
"hLogECellNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
2998 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()),
2999 nptbins,ptmin,ptmax,150,-3,3);
3005 (Form(
"hLogECellNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3006 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()),
3007 nptbins,ptmin,ptmax,150,-3,3);
3013 (Form(
"hLogECellNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3014 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()),
3015 nptbins,ptmin,ptmax,150,-3,3);
3021 (Form(
"hECellWeightNLM1TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3022 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()),
3023 nptbins,ptmin,ptmax,90,0,4.5);
3029 (Form(
"hECellWeightNLM2TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3030 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()),
3031 nptbins,ptmin,ptmax,90,0,4.5);
3037 (Form(
"hECellWeightNLM3TCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3038 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()),
3039 nptbins,ptmin,ptmax,90,0,4.5);
3045 (Form(
"hMassEClusTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3046 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()),
3047 nptbins,ptmin,ptmax,nmassbins,massmin,massmax);
3061 (Form(
"hTimeDiffTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3062 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()),
3063 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
3069 (Form(
"hColRowTCardCorrelNCellLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3070 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3071 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3077 (Form(
"hColRowTCardCorrelNCellHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3078 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3079 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3087 (Form(
"hExoticTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3088 Form(
"exoticity vs #it{E}, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3089 nptbins,ptmin,ptmax,200,-1,1);
3095 (Form(
"hTimeDiffExoTCardCorrelNCell_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3096 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()),
3097 nptbins,ptmin,ptmax,tdbins,tdmin,tdmax);
3103 (Form(
"hColRowTCardCorrelNCellExoticLowE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3104 Form(
"column vs row, N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3105 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3111 (Form(
"hColRowTCardCorrelNCellExoticHighE_Same%d_Diff%d%s",i,j,add[tm].
Data()),
3112 Form(
"column vs row,N cells with w > 0.01, TCard same = %d, diff =%d %s",i,j,add[tm].
Data()),
3113 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax);
3328 (Form(
"hLambda0Exoticity_EBin%d%s",i,add[tm].
Data()),
3329 Form(
"#lambda^{2}_{0} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3330 200,-1,1,ssbins,ssmin,ssmax);
3336 (Form(
"hLambda1Exoticity_EBin%d%s",i,add[tm].
Data()),
3337 Form(
"#lambda^{2}_{1} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3338 200,-1,1,ssbins,ssmin,ssmax);
3352 (Form(
"hNCellsExoticity_EBin%d%s",i,add[tm].
Data()),
3353 Form(
"#it{n}_{cells} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3354 200,-1,1,nceclbins,nceclmin,nceclmax);
3355 fhNCellsExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
3360 (Form(
"hTimeExoticity_EBin%d%s",i,add[tm].
Data()),
3362 200,-1,1,100,-25,25);
3363 fhTimeExoticity[i][tm]->SetXTitle(
"#it{F}_{+}=1-#it{E}_{+}/#it{E}_{lead cell}");
3368 (Form(
"hNCellsTCardSameAndDiff_Exotic_EBin%d%s",i,add[tm].
Data()),
3369 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()),
3370 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
3410 (Form(
"hNCellsTCardSameAndDiff_EBin%d%s",i,add[tm].
Data()),
3411 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()),
3412 nceclbins,nceclmin,nceclmax,nceclbins,nceclmin,nceclmax);
3418 (Form(
"hLambda0Lambda1_EBin%d%s",i,add[tm].
Data()),
3419 Form(
"#lambda^{2}_{0} vs #lambda^{2}_{1}, %2.2f<#it{E}<%2.2f GeV %s",
fEBinCuts[i],
fEBinCuts[i+1],add[tm].
Data()),
3420 ssbins,ssmin,ssmax,ssbins,ssmin,ssmax);
3436 for(
Int_t j = 0; j < 6; j++)
3438 for(
Int_t k = 0; k < 6; k++)
3441 (Form(
"hLambda0Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
3442 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()),
3443 200,-1,1,ssbins,ssmin,ssmax);
3449 (Form(
"hLambda1Exoticity_NCell_Same%d_Diff%d%s",j,k,add[tm].
Data()),
3450 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()),
3451 200,-1,1,ssbins,ssmin,ssmax);
3548 for(
Int_t i = 0; i < 12; i++)
3551 (Form(
"hTCardCorrECellMaxDiff_Case%d%s",i,add[tm].
Data()),
3552 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()),
3553 nptbins,ptmin,ptmax,210,-1,20);
3559 (Form(
"hTCardCorrEClusterDiff_Case%d%s",i,add[tm].
Data()),
3560 Form(
"#it{E}_{cluster}-#it{E}_{cell} vs #it{E}_{cluster}, for (un)correlated cells in TCard, case %d %s",i,add[tm].
Data()),
3561 nptbins,ptmin,ptmax,210,-1,20);
3583 (Form(
"hTCardCorrTCellMaxDiff_Case%d%s",i,add[tm].
Data()),
3584 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()),
3585 nptbins,ptmin,ptmax,1000,-100,100);
3593 (Form(
"hTCardCorrECellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
3594 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()),
3595 nptbins,ptmin,ptmax,210,-1,20);
3601 (Form(
"hTCardCorrEClusterDiffExo_Case%d%s",i,add[tm].
Data()),
3602 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()),
3603 nptbins,ptmin,ptmax,210,-1,20);
3625 (Form(
"hTCardCorrTCellMaxDiffExo_Case%d%s",i,add[tm].
Data()),
3626 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()),
3627 nptbins,ptmin,ptmax,1000,-100,100);
3638 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3644 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3650 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3656 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3662 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3668 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3674 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3680 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3688 (Form(
"hTMPhiResidual_EBin%d",i),
3689 Form(
"#Delta #varphi_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3690 200,-1,1,nresphibins,resphimin,resphimax);
3696 (Form(
"hTMEtaResidual_EBin%d",i),
3697 Form(
"#Delta #eta_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3698 200,-1,1,nresetabins,resetamin,resetamax);
3704 (Form(
"hTMPhiResidual_LooseCut_EBin%d",i),
3705 Form(
"#Delta #varphi_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3706 200,-1,1,nresphibins,resphimin,resphimax);
3712 (Form(
"hTMEtaResidual_LooseCut_EBin%d",i),
3713 Form(
"#Delta #eta_{cluster-track} vs #it{exoticity}, %2.2f<#it{E}<%2.2f GeV",
fEBinCuts[i],
fEBinCuts[i+1]),
3714 200,-1,1,nresetabins,resetamin,resetamax);
3742 (Form(
"hColRowM02"),
3744 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax,40,0.,2.);
3751 (Form(
"hColRowM02NCellCut"),
3752 Form(
"column vs row vs M02, %2.2f<#it{E}<%2.2f GeV #it{n}_{cells}^{w>0.01} > %d for Neutral",
3754 ncolcell,colcellmin,colcellmax,nrowcell,rowcellmin,rowcellmax,40,0.,2.);
3762 Form(
"%2.2f<#it{E}_{1}<%2.2f GeV, %2.2f<#it{E}_{2}<%2.2f GeV, %2.2f<#lambda^{2}_{0}<%2.2f"
3763 "#it{M}_{#gamma #gamma} vs #it{n}_{1, cells}^{w>0.01} vs SM number trig cluster",
3772 (
"hInvMassNCellSMSame",
3773 Form(
"%2.2f<#it{E}_{1}<%2.2f GeV, %2.2f<#it{E}_{2}<%2.2f GeV, %2.2f<#lambda^{2}_{0}<%2.2f"
3774 "#it{M}_{#gamma #gamma} vs #it{n}_{1, cells}^{w>0.01} vs SM number both cluster",
3785 (
"hEMaxCellTimeM02SM",
3786 Form(
"time vs #lambda_{0}^{2} vs SM number, "
3795 (
"hEMaxCellTimeNCellSM",
3796 Form(
"time vs #it{n}_{cells}^{w>0.01} vs SM number, "
3805 (
"hESecCellTimeNCellSM",
3806 Form(
"secondary cell time vs #it{n}_{cells}^{w>0.01} vs SM number, "
3816 Int_t nbinsdeltacells = 19 ;
3819 for(
Int_t col = 0; col < 2; col++)
3822 (Form(
"hColRowFromCellMaxLowM02_Col%d",col),
3823 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, "
3825 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells);
3831 (Form(
"hColRowFromCellMaxHighM02_Col%d",col),
3832 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, "
3834 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells);
3845 (Form(
"hNCellsPerClusterM02M20_SM%d",i),
3846 Form(
" vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01} vs #lambda_{1}^{2}, "
3848 (
Int_t)nShShBins/1.5,0,maxShSh/1.5, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3855 (Form(
"hNCellsPerClusterM02NLM_SM%d",i),
3856 Form(
" vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01} vs N_{lm}, "
3858 10,0,10, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3865 (Form(
"hEMaxCellEClusterM02NCell_SM%d",i),
3866 Form(
"(#it{E}_{cluster} - #it{E}_{max cell})/#it{E}_{cluster} vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01}, "
3868 50,0,1., cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3875 (Form(
"hEMaxCellLogM02NCell_SM%d",i),
3876 Form(
"log(#it{E}_{max cell}) vs #lambda_{0}^{2} vs SM number, "
3878 150,-3,3, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3885 (Form(
"hESecCellEMaxCellM02NCell_SM%d",i),
3886 Form(
"#it{E}_{cell}/#it{E}_{cell max} vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01}, "
3888 50,0,1., cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3895 (Form(
"hESecCellEClusterM02NCell_SM%d",i),
3896 Form(
"(#it{E}_{cluster} - #it{E}_{cell})/#it{E}_{cluster} vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01}, "
3898 50,0,1., cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3905 (Form(
"hESecCellLogM02NCell_SM%d",i),
3906 Form(
"log(#it{E}_{cell}) vs #lambda_{0}^{2} vs #it{n}_{cells}^{w>0.01}, "
3908 150,-3,3, cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3915 (Form(
"hEMaxESecCellNCellLowM02_SM%d",i),
3916 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, "
3918 24,0,12, 48,0,12,cellBins,cellMin,cellMax);
3925 (Form(
"hEMaxECrossNCellLowM02_SM%d",i),
3926 Form(
"#it{E}_{cell}^{max} vs exoticity vs #it{n}_{cells}^{w>0.01} vs 0.1 < #lambda_{0}^{2} < 0.3, "
3928 24,0,12, 40,0.6,1.,cellBins,cellMin,cellMax);
3935 (Form(
"hEMaxESecCellNCellHighM02_SM%d",i),
3936 Form(
"#it{E}_{cell}^{max} vs #it{E}_{cell}^{secondary} vs #it{n}_{cells}^{w>0.01} vs 0.5 < #lambda_{0}^{2} < 2, "
3938 24,0,12, 48,0,12,cellBins,cellMin,cellMax);
3945 (Form(
"hEMaxECrossNCellHighM02_SM%d",i),
3946 Form(
"#it{E}_{cell}^{max} vs exoticity vs #it{n}_{cells}^{w>0.01} for 0.5 < #lambda_{0}^{2} < 2, "
3948 24,0,12, 40,0.6,1.,cellBins,cellMin,cellMax);
3954 Int_t nbinsdeltacells = 19 ;
3957 for(
Int_t col = 0; col < 2; col++)
3960 (Form(
"hNCellsPerClusterM02_SM%d_Col%d",i,col),
3961 Form(
"#it{E} vs #it{n}_{cells} vs #lambda_{0}^{2}, SM=%d, colum %d",i,col),
3962 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
3969 (Form(
"hNCellsPerClusterM20_SM%d_Col%d",i,col),
3970 Form(
"#it{E} vs #it{n}_{cells} vs #lambda_{1}^{2}, SM=%d, colum %d",i,col),
3971 nEbins,minE,maxE,cellBins,cellMin,cellMax,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
3978 (Form(
"hNCellsPerClusterNLM_SM%d_Col%d",i,col),
3979 Form(
"#it{E} vs #it{n}_{cells} vs n_{lm}, SM=%d, colum %d",i,col),
3980 nEbins,minE,maxE,cellBins,cellMin,cellMax,10,0,10);
3987 (Form(
"hColRowFromCellMaxLowM02_SM%d_Col%d",i,col),
3988 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, "
3990 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,cellBins,cellMin,cellMax);
3997 (Form(
"hColRowFromCellMaxHighM02_SM%d_Col%d",i,col),
3998 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, "
4000 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,cellBins,cellMin,cellMax);
4006 for(
Int_t j = 0; j < 3; j++)
4009 (Form(
"hColRowFromCellMaxEMaxSecDiffLowM02_SM%d_Col%d_NCellBin%d",i,col,j),
4010 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, "
4011 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
4012 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,120,0,12);
4019 (Form(
"hColRowFromCellMaxEMaxSecDiffHighM02_SM%d_Col%d_NCellBin%d",i,col,j),
4020 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, "
4021 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
4022 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,120,0,12);
4029 (Form(
"hColRowFromCellMaxEMaxSecDiffFracLowM02_SM%d_Col%d_NCellBin%d",i,col,j),
4030 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, "
4031 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
4032 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,40,0,1);
4039 (Form(
"hColRowFromCellMaxEMaxSecDiffFracHighM02_SM%d_Col%d_NCellBin%d",i,col,j),
4040 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, "
4041 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
4042 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,40,0,1);
4049 (Form(
"hColRowFromCellMaxECellClusterRatLowM02_SM%d_Col%d_NCellBin%d",i,col,j),
4050 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, "
4051 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
4052 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,50,0,1);
4059 (Form(
"hColRowFromCellMaxECellClusterRatHighM02_SM%d_Col%d_NCellBin%d",i,col,j),
4060 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, "
4061 "%2.2f<#it{E}<%2.2f GeV, SM=%d, colum %d, #it{n}_{cells}^{w>0.01} bin %d",
fEMinShape,
fEMaxShape,i,col,j),
4062 nbinsdeltacells,mindeltacells,maxdeltacells,nbinsdeltacells,mindeltacells,maxdeltacells,50,0,1);
4073 for(
Int_t imatch = 0; imatch < 3; imatch++)
4078 (Form(
"hDeltaIEtaDeltaIPhi_%s",matchCase[imatch].
Data()),
4079 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),
4080 nEbins,minE,maxE,cellBins,cellMin,cellMax,cellBins,cellMin,cellMax);
4087 (Form(
"hDeltaIEtaDeltaIPhiTot_%s",matchCase[imatch].
Data()),
4088 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),
4089 nEbins,minE,maxE,cellBins,cellMin,cellMax,cellBins,cellMin,cellMax);
4096 (Form(
"hDeltaIA_%s",matchCase[imatch].
Data()),
4097 Form(
"Cluster *asymmetry* in cell units vs E for %s",matchCase[imatch].
Data()),
4098 nptbins,ptmin,ptmax,asyBins,asyMin,asyMax);
4099 fhDeltaIA[imatch]->SetXTitle(
"#it{E}_{cluster}");
4100 fhDeltaIA[imatch]->SetYTitle(
"#it{A}_{cell in cluster}");
4101 outputContainer->Add(
fhDeltaIA[imatch]);
4104 (Form(
"hDeltaIATot_%s",matchCase[imatch].
Data()),
4105 Form(
"Cluster *total asymmetry* in cell units vs E for %s",matchCase[imatch].
Data()),
4106 nptbins,ptmin,ptmax,asyBins,asyMin,asyMax);
4108 fhDeltaIATot[imatch]->SetYTitle(
"#it{A}_{cell in cluster}^{total}");
4112 (Form(
"hDeltaIAM02_%s",matchCase[imatch].
Data()),
4113 Form(
"Cluster *asymmetry* in cell units vs #lambda^{2}_{0} for %s",matchCase[imatch].
Data()),
4114 nEbins,minE,maxE,nShShBins,minShSh,maxShSh,asyBins,asyMin,asyMax);
4117 fhDeltaIAM02[imatch]->SetZTitle(
"#it{A}_{cell in cluster}");
4121 (Form(
"hDeltaIATotM02_%s",matchCase[imatch].
Data()),
4122 Form(
"Cluster *total asymmetry* in cell units vs #lambda^{2}_{0} for %s",matchCase[imatch].
Data()),
4123 nEbins,minE,maxE,nShShBins,minShSh,maxShSh,asyBins,asyMin,asyMax);
4126 fhDeltaIATotM02[imatch]->SetZTitle(
"#it{A}_{cell in cluster}^{total}");
4130 (Form(
"hDeltaIAM20_%s",matchCase[imatch].
Data()),
4131 Form(
"Cluster *asymmetry* in cell units vs #lambda^{2}_{1} for %s",matchCase[imatch].
Data()),
4132 nEbins,minE,maxE,(
Int_t)nShShBins/1.5,0,maxShSh/1.5,asyBins,asyMin,asyMax);
4135 fhDeltaIAM20[imatch]->SetZTitle(
"#it{A}_{cell in cluster}");
4139 (Form(
"hDeltaIATotM20_%s",matchCase[imatch].
Data()),
4140 Form(
"Cluster *total asymmetry* in cell units vs #lambda^{2}_{1} for %s",matchCase[imatch].
Data()),
4141 nEbins,minE,maxE,(
Int_t)nShShBins/1.5,0,maxShSh/1.5,asyBins,asyMin,asyMax);
4144 fhDeltaIATotM20[imatch]->SetZTitle(
"#it{A}_{cell in cluster}^{total}");
4148 (Form(
"hDeltaIANCells_%s",matchCase[imatch].
Data()),
4149 Form(
"Cluster *asymmetry* in cell units vs N cells in cluster for %s",matchCase[imatch].
Data()),
4150 nEbins,minE,maxE,cellBins,cellMin,cellMax,asyBins,asyMin,asyMax);
4157 (Form(
"hDeltaIATotNCells_%s",matchCase[imatch].
Data()),
4158 Form(
"Cluster *total asymmetry* in cell units vs N cells in cluster for %s",matchCase[imatch].
Data()),
4159 nEbins,minE,maxE,cellBins,cellMin,cellMax,asyBins,asyMin,asyMax);
4168 (Form(
"hDeltaIAOrigin_%s",matchCase[imatch].
Data()),
4169 Form(
"Cluster *asymmetry* in cell units vs E vs primary ID, for ID %s",matchCase[imatch].
Data()),
4170 nEbins,minE,maxE,10,-0.5,9.5,asyBins,asyMin,asyMax);
4174 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4175 fhDeltaIAOrigin[imatch]->GetYaxis()->SetBinLabel(ilabel,mcParticleStringLabel[ilabel-1]);
4179 (Form(
"hDeltaIATotOrigin_%s",matchCase[imatch].
Data()),
4180 Form(
"Cluster *total asymmetry* in cell units vs E vs primary ID, for ID %s",matchCase[imatch].
Data()),
4181 nEbins,minE,maxE,10,-0.5,9.5,asyBins,asyMin,asyMax);
4185 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4186 fhDeltaIATotOrigin[imatch]->GetYaxis()->SetBinLabel(ilabel,mcParticleStringLabel[ilabel-1]);
4192 (Form(
"hClusterMaxCellDiffM02_%s",matchCase[imatch].
Data()),
4193 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()),
4194 nEbins,minE,maxE, 20,0,1.,nShShBins,minShSh,maxShSh);
4201 (Form(
"hClusterTimeEnergyM02_%s",matchCase[imatch].
Data()),
4202 Form(
"#it{E} vs TOF vs #lambda_{0}^{2} for ID %s",matchCase[imatch].
Data()),
4203 nEbins,minE,maxE,45,-25.5,20.5,nShShBins,minShSh,maxShSh);
4210 (Form(
"hNCellsPerClusterM02_%s",matchCase[imatch].
Data()),
4211 Form(
"#it{E} vs #it{n}_{cells} vs #lambda_{0}^{2} for ID %s",matchCase[imatch].
Data()),
4212 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4219 (Form(
"hNCellsPerClusterM20_%s",matchCase[imatch].
Data()),
4220 Form(
"#it{E} vs #it{n}_{cells} vs #lambda_{1}^{2} for ID %s",matchCase[imatch].
Data()),
4221 nEbins,minE,maxE,cellBins,cellMin,cellMax,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4228 (Form(
"hNCellsPerClusterNLM_%s",matchCase[imatch].
Data()),
4229 Form(
"#it{E} vs #it{n}_{cells} vs n_{lm} for ID %s",matchCase[imatch].
Data()),
4230 nEbins,minE,maxE,cellBins,cellMin,cellMax,10,0,10);
4237 (Form(
"hSMM02_%s",matchCase[imatch].
Data()),
4238 Form(
"#it{E} vs SM number vs #lambda_{0}^{2}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4239 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,minShSh,maxShSh);
4240 fhSMM02[imatch]->SetXTitle(
"#it{E} (GeV)");
4241 fhSMM02[imatch]->SetYTitle(
"SM number");
4242 fhSMM02[imatch]->SetZTitle(
"#lambda_{0}^{2}");
4243 outputContainer->Add(
fhSMM02[imatch]);
4246 (Form(
"hSMM02NoCut_%s",matchCase[imatch].
Data()),
4247 Form(
"#it{E} vs SM number vs #lambda_{0}^{2} for ID %s",matchCase[imatch].
Data()),
4248 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,minShSh,maxShSh);
4255 (Form(
"hSMM20LowM02_%s",matchCase[imatch].
Data()),
4256 Form(
"#it{E} vs SM number vs #lambda_{1}^{2}, #it{n}_{cells}^{w>0.01}>%d, for 0.1<#lambda_{0}^{2}<0.3, for ID %s",
4258 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4265 (Form(
"hSMM20LowM02NoCut_%s",matchCase[imatch].
Data()),
4266 Form(
"#it{E} vs SM number vs #lambda_{1}^{2}, for 0.1<#lambda_{0}^{2}<0.3, for ID %s",matchCase[imatch].
Data()),
4267 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4274 (Form(
"hSMM20HighM02_%s",matchCase[imatch].
Data()),
4275 Form(
"#it{E} vs SM number vs #lambda_{1}^{2}, #it{n}_{cells}^{w>0.01}>%d, for 0.5<#lambda_{0}^{2}<2, for ID %s",
4277 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4284 (Form(
"hSMM20HighM02NoCut_%s",matchCase[imatch].
Data()),
4285 Form(
"#it{E} vs SM number vs #lambda_{1}^{2}, for 0.5<#lambda_{0}^{2}<2, for ID %s",matchCase[imatch].
Data()),
4286 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,(
Int_t)nShShBins/1.5,0,maxShSh/1.5);
4293 (Form(
"hSMNCell_%s",matchCase[imatch].
Data()),
4294 Form(
"#it{E} vs SM number vs #it{n}_{cells}^{w>0.01} for ID %s",matchCase[imatch].
Data()),
4295 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,cellBins,cellMin,cellMax);
4296 fhSMNCell[imatch]->SetXTitle(
"#it{E} (GeV)");
4297 fhSMNCell[imatch]->SetYTitle(
"SM number");
4298 fhSMNCell[imatch]->SetZTitle(
"#it{n}_{cells}^{w>0.01}");
4299 outputContainer->Add(
fhSMNCell[imatch]);
4302 (Form(
"hSMNCellM02_%s",matchCase[imatch].
Data()),
4303 Form(
"SM number vs #it{n}_{cells}^{w>0.01} vs #lambda_{0}^{2}, "
4305 fNModules,-0.5,fNModules-0.5,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4308 fhSMNCellM02[imatch]->SetYTitle(
"#it{n}_{cells}^{w>0.01}");
4312 (Form(
"hColM02_%s",matchCase[imatch].
Data()),
4313 Form(
"#it{E} vs column number vs #lambda_{0}^{2}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4314 nEbins,minE,maxE,48,-0.5,47.5,nShShBins,minShSh,maxShSh);
4315 fhColM02[imatch]->SetXTitle(
"#it{E} (GeV)");
4316 fhColM02[imatch]->SetYTitle(
"column number");
4317 fhColM02[imatch]->SetZTitle(
"#lambda_{0}^{2}");
4318 outputContainer->Add(
fhColM02[imatch]);
4321 (Form(
"hRowM02_%s",matchCase[imatch].
Data()),
4322 Form(
"#it{E} vs row number vs #lambda_{0}^{2}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4323 nEbins,minE,maxE,24,-0.5,23.5,nShShBins,minShSh,maxShSh);
4324 fhRowM02[imatch]->SetXTitle(
"#it{E} (GeV)");
4325 fhRowM02[imatch]->SetYTitle(
"row number");
4326 fhRowM02[imatch]->SetZTitle(
"#lambda_{0}^{2}");
4327 outputContainer->Add(
fhRowM02[imatch]);
4332 (Form(
"hNCellsPerClusterMEta_%s",matchCase[imatch].
Data()),
4333 Form(
"#it{E} vs #it{n}_{cells} vs #sigma_{#eta}^{2} for ID %s",matchCase[imatch].
Data()),
4334 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4341 (Form(
"hNCellsPerClusterMPhi_%s",matchCase[imatch].
Data()),
4342 Form(
"#it{E} vs #it{n}_{cells} vs #sigma_{#varphi}^{2} for ID %s",matchCase[imatch].
Data()),
4343 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,minShSh,maxShSh);
4350 (Form(
"hNCellsPerClusterMEtaPhi_%s",matchCase[imatch].
Data()),
4351 Form(
"#it{E} vs #it{n}_{cells} vs #sigma_{#eta#varphi}^{2} for ID %s",matchCase[imatch].
Data()),
4352 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,-1*maxShSh,maxShSh);
4359 (Form(
"hNCellsPerClusterMEtaPhiA_%s",matchCase[imatch].
Data()),
4360 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()),
4361 nEbins,minE,maxE,cellBins,cellMin,cellMax,nShShBins,-1*maxShSh,maxShSh);
4364 fhNCellsPerClusterMEtaPhiA[imatch]->SetZTitle(
"(#sigma_{#varphi}^{2}-#sigma_{#eta}^{2})/(#sigma_{#varphi}^{2}+#sigma_{#eta}^{2})");
4368 (Form(
"hSMMEta_%s",matchCase[imatch].
Data()),
4369 Form(
"#it{E} vs SM number vs #sigma_{#eta}^{2}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4370 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,minShSh,maxShSh);
4371 fhSMMEta[imatch]->SetXTitle(
"#it{E} (GeV)");
4372 fhSMMEta[imatch]->SetYTitle(
"SM number");
4373 fhSMMEta[imatch]->SetZTitle(
"#sigma_{#eta}^{2}");
4374 outputContainer->Add(
fhSMMEta[imatch]);
4377 (Form(
"hSMMEtaNoCut_%s",matchCase[imatch].
Data()),
4378 Form(
"#it{E} vs SM number vs #sigma_{#eta}^{2} for ID %s",matchCase[imatch].
Data()),
4379 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,minShSh,maxShSh);
4386 (Form(
"hSMMPhi_%s",matchCase[imatch].
Data()),
4387 Form(
"#it{E} vs SM number vs #sigma_{#varphi}^{2}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4388 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,minShSh,maxShSh);
4389 fhSMMPhi[imatch]->SetXTitle(
"#it{E} (GeV)");
4390 fhSMMPhi[imatch]->SetYTitle(
"SM number");
4391 fhSMMPhi[imatch]->SetZTitle(
"#sigma_{#varphi}^{2}");
4392 outputContainer->Add(
fhSMMPhi[imatch]);
4395 (Form(
"hSMMPhiNoCut_%s",matchCase[imatch].
Data()),
4396 Form(
"#it{E} vs SM number vs #sigma_{#varphi}^{2} for ID %s",matchCase[imatch].
Data()),
4397 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,minShSh,maxShSh);
4404 (Form(
"hSMMEtaPhi_%s",matchCase[imatch].
Data()),
4405 Form(
"#it{E} vs SM number vs #sigma_{#eta#varphi}^{2}, #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4406 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,minShSh,maxShSh);
4409 fhSMMEtaPhi[imatch]->SetZTitle(
"#sigma_{#eta#varphi}^{2}");
4413 (Form(
"hSMMEtaPhiNoCut_%s",matchCase[imatch].
Data()),
4414 Form(
"#it{E} vs SM number vs #sigma_{#eta#varphi}^{2} for ID %s",matchCase[imatch].
Data()),
4415 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,minShSh,maxShSh);
4422 (Form(
"hSMMEtaPhiA_%s",matchCase[imatch].
Data()),
4423 Form(
"#it{E} vs SM number vs (#sigma_{#varphi}^{2}-#sigma_{#eta}^{2})/(#sigma_{#varphi}^{2}+#sigma_{#eta}^{2}), #it{n}_{cells}^{w>0.01}>%d, for ID %s",
fNCellMinShape,matchCase[imatch].
Data()),
4424 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,-1*maxShSh,maxShSh);
4427 fhSMMEtaPhiA[imatch]->SetZTitle(
"(#sigma_{#varphi}^{2}-#sigma_{#eta}^{2})/(#sigma_{#varphi}^{2}+#sigma_{#eta}^{2})");
4431 (Form(
"hSMMEtaPhiANoCut_%s",matchCase[imatch].
Data()),
4432 Form(
"#it{E} vs SM number vs (#sigma_{#varphi}^{2}-#sigma_{#eta}^{2})/(#sigma_{#varphi}^{2}+#sigma_{#eta}^{2}) for ID %s",matchCase[imatch].
Data()),
4433 nEbins,minE,maxE,fNModules,-0.5,fNModules-0.5,nShShBins,-1*maxShSh,maxShSh);
4436 fhSMMEtaPhiANoCut[imatch]->SetZTitle(
"(#sigma_{#varphi}^{2}-#sigma_{#eta}^{2})/(#sigma_{#varphi}^{2}+#sigma_{#eta}^{2})");
4443 (Form(
"hOrigin_%s",matchCase[imatch].
Data()),
4444 Form(
"#it{E} vs origin for ID %s",matchCase[imatch].
Data()),
4445 nEbins,minE,maxE, 10,-0.5,9.5);
4446 fhOriginE[imatch]->SetXTitle(
"#it{E} (GeV)");
4447 fhOriginE[imatch]->SetYTitle(
"particle");
4448 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4449 fhOriginE[imatch]->GetYaxis()->SetBinLabel(ilabel ,mcParticleStringLabel[ilabel-1]);
4450 outputContainer->Add(
fhOriginE[imatch]);
4453 (Form(
"hOriginM02_%s",matchCase[imatch].
Data()),
4454 Form(
"#it{E} vs origin vs #lambda_{0}^{2} for ID %s",matchCase[imatch].
Data()),
4455 nEbins,minE,maxE,10,-0.5,9.5,nShShBins,minShSh,maxShSh);
4458 fhOriginM02[imatch]->SetZTitle(
"#lambda_{0}^{2}");
4459 for(
Int_t ilabel = 1; ilabel <=10; ilabel++)
4460 fhOriginM02[imatch]->GetYaxis()->SetBinLabel(ilabel,mcParticleStringLabel[ilabel-1]);
4475 (
"hClusterMaxCellCloseCellRatioM02",
"#it{E} vs #it{E}_{cell-i}/#it{E}_{cell max} vs #lambda_{0}^{2}",
4476 nEbins,minE,maxE, 20,0,1.,nShShBins,minShSh,maxShSh);
4494 (
"hClusterMaxCellECrossM02",
4495 "#it{E} vs exoticity vs #lambda_{0}^{2}",
4496 nEbins,minE,maxE, 40,0.6,1.,nShShBins,minShSh,maxShSh);
4507 nptbins,ptmin,ptmax, 100,0,1.);
4513 nptbins,ptmin,ptmax, 100,-10,0);
4519 nptbins,ptmin,ptmax, 100,0,1.);
4525 nptbins,ptmin,ptmax, 100,-10,0);
4531 nptbins*2,ptmin,ptmax*2, 100,0,1.);
4537 nptbins*2,ptmin,ptmax*2, 100,-10,0);
4547 if(imod < fFirstModule || imod >
fLastModule)
continue;
4550 Form(
"#cell energy / sum all energy vs all energy in Module %d",imod),
4551 nptbins*2,ptmin,ptmax*2, 100,0,1.);
4557 Form(
"Log(cell energy / sum all energy) vs all energy in Module %d",imod),
4558 nptbins*2,ptmin,ptmax*2, 100,-10,0);
4565 Float_t cellEmin [] = {0.05,0.1,0.15,0.2};
4566 Float_t cellTmin [] = {50.,100.,1000000.};
4568 for(
Int_t iw = 0; iw < 12; iw++)
4571 for(
Int_t iEmin = 0; iEmin < 4; iEmin++)
4573 for(
Int_t iTmin = 0; iTmin < 3; iTmin++)
4576 Form(
"#lambda^{2}_{0} vs E, w0=%1.2f, cell E>%2.2f MeV, |t|<%2.0f ns",
4577 w0, cellEmin[iEmin], cellTmin[iTmin]),
4578 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4593 Form(
"#lambda^{2}_{0} vs E, w0=%1.2f, cell E>%2.2f MeV, |t|<%2.0f ns, |#eta| < 0.15",
4594 w0, cellEmin[iEmin], cellTmin[iTmin]),
4595 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4604 TString mcnames[] = {
"Photon",
"Electron",
"Conversion",
"Pi0",
"Hadron"};
4605 for(
Int_t imc = 0; imc < 5; imc++)
4608 Form(
"shower shape, #lambda^{2}_{0} vs E, w0 = %1.1f, for MC %s",w0,mcnames[imc].
Data()),
4609 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4628 return outputContainer;
4689 AliInfo(Form(
"dEdX cuts init for MC: %d<dedx_ele<%d; %d<dedx_had<%d",
4697 if ( runNumber < 146861 )
4704 else if ( runNumber < 156000 )
4711 else if ( runNumber < 165000 )
4726 AliInfo(Form(
"dEdX cuts init for run %d: %d<dedx_ele<%d; %d<dedx_had<%d",
4740 Int_t nCellsPerCluster)
4758 if(ampMax < 0.01)
return kFALSE;
4774 printf(
"**** Print %s %s ****\n", GetName(), GetTitle() ) ;
4778 printf(
"Min n cells : %d\n" ,
fNCellMin) ;
4780 printf(
"Min M02 : %1.2f\n",
fM02Min) ;
4796 AliDebug(1,
"Start");
4812 AliFatal(Form(
"AliAnaClusterShapeCorrelStudies::MakeAnalysisFillHistograms() - Wrong calorimeter name <%s>, END",
GetCalorimeterString().
Data()));
4816 AliFatal(Form(
"AliAnaClusterShapeCorrelStudies::MakeAnalysisFillHistograms() - No CaloClusters or CaloCells available"));
4846 if(clus->GetNCells() < 3)
return ;
4848 for (
Int_t ipos = 0; ipos < clus->GetNCells(); ipos++)
4850 Int_t id = clus->GetCellsAbsId()[ipos];
4865 AliWarning(Form(
"Wrong calculated energy %f",energy));
4876 for (
Int_t ipos = 0; ipos < clus->GetNCells(); ipos++)
4878 Int_t id = clus->GetCellsAbsId()[ipos];
4893 Float_t l0org = clus->GetM02();
4894 Float_t l1org = clus->GetM20();
4895 Float_t dorg = clus->GetDispersion();
4899 if(
IsDataMC() && clus->GetNLabels() > 0)
4924 Float_t cellEmin [] = {0.05,0.1,0.15,0.2};
4925 Float_t cellTmin [] = {50.,100.,1000000.};
4927 for(
Int_t iw = 0; iw < 12; iw++)
4929 for(
Int_t iEmin = 0; iEmin < 4; iEmin++)
4931 for(
Int_t iTmin = 0; iTmin < 3; iTmin++)