PWG4/macros/ConfigJetCorrel.C

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