PWG4/macros/ConfigJetCorrel.C

   1 AliJetCorrelSelector* ConfigJetCorrel(){
   2
   3   ///////////////////////////////////
   4   // set correlation input parameters
   5   ///////////////////////////////////
   6   // set generic selections:
   7   Bool_t kGenQA = kFALSE;      // generate QA histos
   8   UInt_t kDPhiNumBins = 30;    // number of bins in DeltaPhi
   9   UInt_t kDEtaNumBins = 16;    // number of bins in DeltaEta
  10   Float_t kPoutBinWidth = 0.2; // number of bins in Pout = max(assocBins)/kPoutBinWidth
  11   Bool_t kUseAliKF = kFALSE;   // use AliKF or TLorentzVector for parent reconstruction
  12   UInt_t poolDepth = 5;
  13   UInt_t correlTypes[] = {0};  // 0=dihadron, 1=pi0-hadron, 2=photon-hadron
  14   Float_t centrBins[] = {1,500};
  15   Float_t zVertBins[] = {-5,5};
  16   Float_t triggBins[] = {6,8,10,20};
  17   Float_t assocBins[] = {2,10};
  18   TString sTrigg[] = {"ALL"}; // selects events where one of the strings is matched; "ALL"=no cut
  19   //TString sTrigg[] = {"CINT1B-"};
  20   // set track selections:
  21   Bool_t itsRefit = kTRUE;
  22   Bool_t tpcRefit = kTRUE;
  23   Bool_t trdRefit = kTRUE;         // used only for electron tracks
  24   Float_t maxEta = 0.8;
  25   UInt_t minNClusTPC = 70;
  26   Float_t maxTPCChi2 = 4.0;        // max track Chi2 per TPC cluster
  27   Bool_t rejectKinkChild = kTRUE;  // reject track comming from a kink
  28   Float_t trkPairCut = 0.;         // track pair proximity cut (dist at TPC entrance)
  29   // code that applies next 3 cuts (NClusITS,ITSChi2,NsigVtx) currently commented out
  30   UInt_t minNClusITS = 0;
  31   Float_t maxITSChi2 = 35;        // max track Chi2 per ITS cluster
  32   Float_t maxNsigVtx = 35;        // max dist to primary vertex (sigma)
  33   Float_t maxTrkVtx  = 2.4;       // max dist to primary vertex (absolute) - temporarily instead of sigma
  34
  35   //////////////////////////////////
  36   // load them into selector object:
  37   //////////////////////////////////
  38   AliJetCorrelSelector* selector = new AliJetCorrelSelector();
  39   selector->SetQA(kGenQA);
  40   selector->SetUseAliKF(kUseAliKF);
  41   selector->SetDPhiNumBins(kDPhiNumBins);
  42   selector->SetDEtaNumBins(kDEtaNumBins);
  43   selector->SetPoutBinWidth(kPoutBinWidth);
  44   selector->SetPoolDepth(poolDepth);
  45   selector->SetCorrelTypes(sizeof(correlTypes)/sizeof(UInt_t),correlTypes);
  46   selector->SetBinningCentr(sizeof(centrBins)/sizeof(Float_t),centrBins);
  47   selector->SetBinningZvert(sizeof(zVertBins)/sizeof(Float_t),zVertBins);
  48   selector->SetBinningTrigg(sizeof(triggBins)/sizeof(Float_t),triggBins);
  49   selector->SetBinningAssoc(sizeof(assocBins)/sizeof(Float_t),assocBins);
  50   selector->SetTriggers(sizeof(sTrigg)/sizeof(TString),sTrigg);
  51   selector->SetITSRefit(itsRefit);
  52   selector->SetTPCRefit(tpcRefit);
  53   selector->SetTRDRefit(trdRefit);
  54   selector->SetMaxEta(maxEta);
  55   selector->SetMinNClusITS(minNClusITS);
  56   selector->SetMinNClusTPC(minNClusTPC);
  57   selector->SetMaxITSChi2(maxITSChi2);
  58   selector->SetMaxTPCChi2(maxTPCChi2);
  59   selector->SetMaxNsigmaVtx(maxNsigVtx);
  60   selector->SetMaxTrkVtx(maxTrkVtx);
  61   selector->SetRejectKinkChild(rejectKinkChild);
  62   selector->SetTrkProximityCut(trkPairCut);
  63   selector->Show();
  64
  65   return selector;
  66 }