PWGDQ/dielectron/AliDielectronVarManager.cxx

   1 /*************************************************************************
   2 * Copyright(c) 1998-2009, ALICE Experiment at CERN, All rights reserved. *
   3 *                                                                        *
   4 * Author: The ALICE Off-line Project.                                    *
   5 * Contributors are mentioned in the code where appropriate.              *
   6 *                                                                        *
   7 * Permission to use, copy, modify and distribute this software and its   *
   8 * documentation strictly for non-commercial purposes is hereby granted   *
   9 * without fee, provided that the above copyright notice appears in all   *
  10 * copies and that both the copyright notice and this permission notice   *
  11 * appear in the supporting documentation. The authors make no claims     *
  12 * about the suitability of this software for any purpose. It is          *
  13 * provided "as is" without express or implied warranty.                  *
  14 **************************************************************************/
  15
  16 ///////////////////////////////////////////////////////////////////////////
  17 //                Dielectron Variables Manager class                     //
  18 //                                                                       //
  19 /*
  20
  21 */
  22 //                                                                       //
  23 ///////////////////////////////////////////////////////////////////////////
  24
  25 #include "AliDielectronVarManager.h"
  26
  27 ClassImp(AliDielectronVarManager)
  28
  29 const char* AliDielectronVarManager::fgkParticleNames[AliDielectronVarManager::kNMaxValues] = {
  30   "Px",
  31   "Py",
  32   "Pz",
  33   "Pt",
  34   "P",
  35   "Xv",
  36   "Yv",
  37   "Zv",
  38   "OneOverPt",
  39   "Phi",
  40   "Theta",
  41   "Eta",
  42   "Y",
  43   "E",
  44   "M",
  45   "Charge",
  46   "NclsITS",
  47   "ITSchi2PerCluster",
  48   "NclsTPC",
  49   "NclsSTPC",
  50   "NclsSFracTPC",
  51   "NclsTPCiter1",
  52   "NFclsTPC",
  53   "NFclsTPCrobust",
  54   "NFclsTPCrobustFraction",
  55   "TPCsignalN",
  56   "TPCsignalNfrac",
  57   "TPCchi2PerCluster",
  58   "TPCclsDiff",
  59   "TrackStatus",
  60
  61   "NclsTRD",
  62   "TRDntracklets",
  63   "TRDpidQuality",
  64   "TRDpidProb_Electrons",
  65   "TRDpidProb_Pions",
  66   "TRDphi",
  67   "TRDpidEffLeg",
  68
  69   "ImpactParXY",
  70   "ImpactParZ",
  71   "TrackLength",
  72
  73   "PdgCode",
  74   "PdgCodeMother",
  75   "PdgCodeGrandMother",
  76   "NumberOfDaughters",
  77   "HaveSameMother",
  78   "IsJpsiPrimary",
  79   "NumberOfJPsisIncl",
  80   "NumberOfJPsisPrompt",
  81   "NumberOfJPsisNPrompt",
  82
  83   "ITS_signal",
  84   "SSD1_signal",
  85   "SSD2_signal",
  86   "SDD1_signal",
  87   "SDD2_signal",
  88   "ITS_clusterMap",
  89   "ITSLayerFirstCls",
  90   "ITS_nSigma_Electrons",
  91   "ITS_nSigma_Pions",
  92   "ITS_nSigma_Muons",
  93   "ITS_nSigma_Kaons",
  94   "ITS_nSigma_Protons",
  95
  96   "P_InnerParam",
  97   "P_OuterParam",
  98   "Y_signed_InnerParam",
  99   "TPC_signal",
 100   "TOF_signal",
 101   "TOF_beta",
 102   "TOF_PIDbit",
 103
 104   "TPC_nSigma_Electrons",
 105   "TPC_nSigma_Pions",
 106   "TPC_nSigma_Muons",
 107   "TPC_nSigma_Kaons",
 108   "TPC_nSigma_Protons",
 109
 110   "TOF_nSigma_Electrons",
 111   "TOF_nSigma_Pions",
 112   "TOF_nSigma_Muons",
 113   "TOF_nSigma_Kaons",
 114   "TOF_nSigma_Protons",
 115
 116   "EMCAL_nSigma_Electrons",
 117   "EMCAL_EoverP",
 118   "EMCAL_NCells",
 119   "EMCAL_M02",
 120   "EMCAL_M20",
 121   "EMCAL_Dispersion",
 122
 123   "KinkIndex0",
 124   //
 125   "Chi2NDF",
 126   "DecayLength",
 127   "R",
 128   "OpeningAngle",
 129   "ThetaHE",
 130   "PhiHE",
 131   "Cos2PhiHE",
 132   "CosTilPhiHE",
 133   "ThetaCS",
 134   "PhiCS",
 135   "Cos2PhiCS",
 136   "CosTilPhiCS",
 137   "DeltaPhiV0ArpH2",
 138   "DeltaPhiV0CrpH2",
 139   "DeltaPhiV0ACrpH2",
 140   "V0ArpH2FlowV2",
 141   "V0CrpH2FlowV2",
 142   "V0ACrpH2FlowV2",
 143   "LegDistance",
 144   "LegDistanceXY",
 145   "DeltaEta",
 146   "DeltaPhi",
 147   "Merr",
 148   "DCA",
 149   "PairType",
 150   "PseudoProperTime",
 151   "PseudoProperTimeResolution",
 152   "PseudoProperTimePull",
 153   "TRDpidEffPair",
 154   //
 155   "X",
 156   "Y",
 157   "Z",
 158   "XRes",
 159   "YRes",
 160   "ZRes",
 161
 162   "v0ArpH2",
 163   "v0CrpH2",
 164   "v0ACrpH2",
 165   "v0ATPCDiffH2",
 166   "v0CTPCDiffH2",
 167   "v0Av0CDiffH2",
 168
 169   "MultV0A",
 170   "MultV0C",
 171   "MultV0",
 172   "AdcV0A",
 173   "AdcV0C",
 174   "AdcV0",
 175   "VZERO_ch0",  "VZERO_ch1",  "VZERO_ch2",  "VZERO_ch3",  "VZERO_ch4",  "VZERO_ch5",  "VZERO_ch6",  "VZERO_ch7",  "VZERO_ch8",  "VZERO_ch9",
 176   "VZERO_ch10", "VZERO_ch11", "VZERO_ch12", "VZERO_ch13", "VZERO_ch14", "VZERO_ch15", "VZERO_ch16", "VZERO_ch17", "VZERO_ch18", "VZERO_ch19",
 177   "VZERO_ch20", "VZERO_ch21", "VZERO_ch22", "VZERO_ch23", "VZERO_ch24", "VZERO_ch25", "VZERO_ch26", "VZERO_ch27", "VZERO_ch28", "VZERO_ch29",
 178   "VZERO_ch30", "VZERO_ch31", "VZERO_ch32", "VZERO_ch33", "VZERO_ch34", "VZERO_ch35", "VZERO_ch36", "VZERO_ch37", "VZERO_ch38", "VZERO_ch39",
 179   "VZERO_ch40", "VZERO_ch41", "VZERO_ch42", "VZERO_ch43", "VZERO_ch44", "VZERO_ch45", "VZERO_ch46", "VZERO_ch47", "VZERO_ch48", "VZERO_ch49",
 180   "VZERO_ch50", "VZERO_ch51", "VZERO_ch52", "VZERO_ch53", "VZERO_ch54", "VZERO_ch55", "VZERO_ch56", "VZERO_ch57", "VZERO_ch58", "VZERO_ch59",
 181   "VZERO_ch60", "VZERO_ch61", "VZERO_ch62", "VZERO_ch63",
 182   "V0AxH2",
 183   "V0AyH2",
 184   "V0ArpH2",
 185   "V0CxH2",
 186   "V0CyH2",
 187   "V0CrpH2",
 188   "V0ACxH2",
 189   "V0ACyH2",
 190   "V0ACrpH2",
 191   "V0ArpResH2",
 192   "V0CrpResH2",
 193   "V0ACrpResH2",
 194   "V0XaXcH2",
 195   "V0XaYaH2",
 196   "V0XaYcH2",
 197   "V0YaXcH2",
 198   "V0YaYcH2",
 199   "V0XcYcH2",
 200   "V0ATPCDiffH2",
 201   "V0CTPCDiffH2",
 202   "V0AV0CDiffH2",
 203   "TPCxH2",
 204   "TPCyH2",
 205   "TPCrpH2",
 206   "TPCsub1xH2",
 207   "TPCsub1yH2",
 208   "TPCsub1rpH2",
 209   "TPCsub2xH2",
 210   "TPCsub2yH2",
 211   "TPCsub2rpH2",
 212   "TPCsub12DiffH2",
 213   "TPCsub12DiffH2Sin",
 214
 215   "TPCxH2uc",
 216   "TPCyH2uc",
 217   "TPCrpH2uc",
 218   "TPCsub1xH2uc",
 219   "TPCsub1yH2uc",
 220   "TPCsub1rpH2uc",
 221   "TPCsub2xH2uc",
 222   "TPCsub2yH2uc",
 223   "TPCsub2rpH2uc",
 224   "TPCsub12DiffH2uc",
 225
 226   "NTrk",
 227   "Tracks",
 228   "NVtxContrib",
 229   "Nacc",
 230   "NaccTrcklts",
 231   "NaccTrcklts0916",
 232
 233   "NaccTrckltsEsd05",
 234   "NaccTrckltsEsd10",
 235   "NaccTrckltsEsd16",
 236   "NaccTrckltsEsd05Corr",
 237   "NaccTrckltsEsd10Corr",
 238   "NaccTrckltsEsd16Corr",
 239   "NaccItsTpcEsd05",
 240   "NaccItsTpcEsd10",
 241   "NaccItsTpcEsd16",
 242   "NaccItsTpcEsd05Corr",
 243   "NaccItsTpcEsd10Corr",
 244   "NaccItsTpcEsd16Corr",
 245   "NaccItsPureEsd05",
 246   "NaccItsPureEsd10",
 247   "NaccItsPureEsd16",
 248   "NaccItsPureEsd05Corr",
 249   "NaccItsPureEsd10Corr",
 250   "NaccItsPureEsd16Corr",
 251
 252   "Nch",
 253   "Nch05",
 254   "Nch10",
 255   "Centrality",
 256   "CentralitySPD",
 257   "Nevents",
 258   "RunNumber",
 259   "MixingBin"
 260 };
 261
 262 AliPIDResponse* AliDielectronVarManager::fgPIDResponse      = 0x0;
 263 AliVEvent*      AliDielectronVarManager::fgEvent            = 0x0;
 264 AliEventplane*  AliDielectronVarManager::fgTPCEventPlane    = 0x0;
 265 AliKFVertex*    AliDielectronVarManager::fgKFVertex         = 0x0;
 266 TProfile*       AliDielectronVarManager::fgMultEstimatorAvg[4][9] = {{0x0}};
 267 TH3D*           AliDielectronVarManager::fgTRDpidEff[10][4] = {{0x0}};
 268 Double_t        AliDielectronVarManager::fgTRDpidEffCentRanges[10][4] = {{0.0}};
 269 TString         AliDielectronVarManager::fgVZEROCalibrationFile = "";
 270 TString         AliDielectronVarManager::fgVZERORecenteringFile = "";
 271 TProfile2D*     AliDielectronVarManager::fgVZEROCalib[64] = {0x0};
 272 TProfile2D*     AliDielectronVarManager::fgVZERORecentering[2][2] = {{0x0,0x0},{0x0,0x0}};
 273 Int_t           AliDielectronVarManager::fgCurrentRun = -1;
 274 Double_t        AliDielectronVarManager::fgData[AliDielectronVarManager::kNMaxValues] = {0.};
 275 //________________________________________________________________
 276 AliDielectronVarManager::AliDielectronVarManager() :
 277   TNamed("AliDielectronVarManager","AliDielectronVarManager")
 278 {
 279   //
 280   // Default constructor
 281   //
 282   for(Int_t i=0; i<4; ++i)
 283     for(Int_t j=0; j<9; ++j)
 284       fgMultEstimatorAvg[i][j] = 0x0;
 285   for(Int_t i=0; i<10; ++i)
 286     for(Int_t j=0; j<4; ++j)
 287       fgTRDpidEff[i][j] = 0x0;
 288   for(Int_t i=0; i<64; ++i) fgVZEROCalib[i] = 0x0;
 289   for(Int_t i=0; i<2; ++i) {
 290     for(Int_t j=0; j<2; ++j) fgVZERORecentering[i][j] = 0x0;
 291   }
 292 }
 293
 294 //________________________________________________________________
 295 AliDielectronVarManager::AliDielectronVarManager(const char* name, const char* title) :
 296   TNamed(name,title)
 297 {
 298   //
 299   // Named constructor
 300   //
 301   for(Int_t i=0; i<4; ++i)
 302     for(Int_t j=0; j<9; ++j)
 303       fgMultEstimatorAvg[i][j] = 0x0;
 304   for(Int_t i=0; i<10; ++i)
 305     for(Int_t j=0; j<4; ++j)
 306       fgTRDpidEff[i][j] = 0x0;
 307   for(Int_t i=0; i<64; ++i) fgVZEROCalib[i] = 0x0;
 308   for(Int_t i=0; i<2; ++i)
 309     for(Int_t j=0; j<2; ++j)
 310       fgVZERORecentering[i][j] = 0x0;
 311 }
 312
 313 //________________________________________________________________
 314 AliDielectronVarManager::~AliDielectronVarManager()
 315 {
 316   //
 317   // Default destructor
 318   //
 319   for(Int_t i=0; i<4; ++i)
 320     for(Int_t j=0; j<9; ++j)
 321       if(fgMultEstimatorAvg[i][j]) delete fgMultEstimatorAvg[i][j];
 322   for(Int_t i=0; i<10; ++i)
 323     for(Int_t j=0; j<4; ++j)
 324       if(fgTRDpidEff[i][j]) delete fgTRDpidEff[i][j];
 325   for(Int_t i=0; i<64; ++i)
 326     if(fgVZEROCalib[i]) delete fgVZEROCalib[i];
 327   for(Int_t i=0; i<2; ++i)
 328     for(Int_t j=0; j<2; ++j)
 329       if(fgVZERORecentering[i][j]) delete fgVZERORecentering[i][j];
 330 }
 331