PWGDQ/dielectron/macrosLMEE/Config_mk_CCbarPbPb.C

   1 void SetupTrackCutsDieleData(AliDielectron *diele, Int_t cutDefinition, Bool_t isAOD);
   2 void SetupPairCutsDieleData(AliDielectron *diele, Int_t cutDefinition, Bool_t isAOD);
   3 void InitHistogramsDieleData(AliDielectron *diele, Int_t cutDefinition, Bool_t isAOD);
   4 void InitCFDieleData(AliDielectron *diele, Int_t cutDefinition, Bool_t isAOD);
   5 void AddMCSignals(AliDielectron *diele);
   6 AliESDtrackCuts *SetupESDtrackCutsDieleData(Int_t cutDefinition);
   7
   8 //TString namesDieleData=("EMCal+QA;TPC+QA;EMCal+ITSfirstclass");
   9 TString namesDieleData=("EMCal+QA;EMCal+QA+MixPlane;EMCal+QA+MixPlane+Rot");
  10 TObjArray *arrNamesDieleData=namesDieleData.Tokenize(";");
  11 Bool_t hasMC=kFALSE;
  12 TString list  = gSystem->Getenv("LIST");
  13 const Int_t nDie=arrNamesDieleData->GetEntries();
  14
  15 AliDielectron* Config_mk_CCbarPbPb(Int_t cutDefinition, Bool_t isAOD=kFALSE)
  16 {
  17   //
  18   // Setup the instance of AliDielectron
  19   //
  20     // find mc or not?
  21 //   if( list.IsNull()) list=prod;
  22   if( list.Contains("LHC10h")   || list.Contains("LHC11h")   ) hasMC=kFALSE;
  23   if( list.Contains("LHC11a10") || list.Contains("LHC12a17") ) hasMC=kTRUE;
  24
  25   // create the actual framework object
  26   TString name=Form("%02d",cutDefinition);
  27   if (cutDefinition<arrNamesDieleData->GetEntriesFast()){
  28     name=arrNamesDieleData->At(cutDefinition)->GetName();
  29   }
  30
  31
  32   AliDielectron *diele = new AliDielectron(Form("%s",name.Data()),
  33                              Form("Track cuts: %s",name.Data()));
  34
  35   diele->SetHasMC(hasMC);
  36
  37 //   printf(" Add %s %s config %s for %s \n",(isESD?"ESD":"AOD"),(hasMC?"MC":""),name.Data(),list.Data());
  38
  39   // Monte Carlo Signals and TRD efficiency tables
  40   if(hasMC) {
  41     AddMCSignals(diele);
  42     printf(" Add %d MC signals \n",diele->GetMCSignals()->GetEntriesFast());
  43   }
  44
  45
  46   // cut setup
  47   SetupEventCutsDieleFilter(diele, cutDefinition, isAOD);
  48   SetupTrackCutsDieleData(diele, cutDefinition, isAOD);
  49   SetupPairCutsDieleData(diele, cutDefinition, isAOD);
  50
  51   //
  52   // histogram setup
  53   // only if an AliDielectronHistos object is attached to the
  54   // dielelectron framework histograms will be filled
  55   //
  56   InitHistogramsDieleData(diele, cutDefinition, isAOD);
  57
  58   // the last definition uses no cuts and only the QA histograms should be filled!
  59   InitCFDieleData(diele, cutDefinition, isAOD);
  60
  61   if(!hasMC){
  62
  63       AliDielectronTrackRotator *rot=new AliDielectronTrackRotator;
  64       rot->SetIterations(10);
  65       rot->SetConeAnglePhi(TMath::Pi());
  66       rot->SetStartAnglePhi(TMath::Pi());
  67       if(cutDefinition==2)  diele->SetTrackRotator(rot);
  68
  69
  70       AliDielectronMixingHandler *mix=new AliDielectronMixingHandler;
  71       mix->AddVariable(AliDielectronVarManager::kZvPrim, 10,-10.,10.);
  72       mix->AddVariable(AliDielectronVarManager::kCentrality, "0,6,12,17,23,27,33,39,47,57,90");
  73       if(cutDefinition==1)mix->AddVariable(AliDielectronVarManager::kV0ArpH2,10,-1.571,1.571);
  74       mix->SetMixType(AliDielectronMixingHandler::kAll);
  75       mix->SetDepth(20);
  76       diele->SetMixingHandler(mix);
  77   }
  78
  79   return diele;
  80 }
  81
  82 //______________________________________________________________________________________
  83 void SetupEventCutsDieleFilter(AliDielectron *diele, Int_t cutDefinition, Bool_t isAOD=kFALSE)
  84 {
  85   //
  86   // Setup the event cuts
  87   //
  88   AliDielectronVarCuts *vtxZ = new AliDielectronVarCuts("vtxZ","Vertex z cut");
  89   vtxZ->AddCut(AliDielectronVarManager::kZvPrim,-10.,10.);
  90
  91   AliDielectronVarCuts *Centrality = new AliDielectronVarCuts("Centrality","Centrality Percentile");
  92   Centrality->AddCut(AliDielectronVarManager::kCentrality,0.,90.);
  93
  94   diele->GetEventFilter().AddCuts(vtxZ);
  95   diele->GetEventFilter().AddCuts(Centrality);
  96
  97 }
  98
  99
 100 //______________________________________________________________________________________
 101 void SetupTrackCutsDieleData(AliDielectron *diele, Int_t cutDefinition, Bool_t isAOD)
 102 {
 103   //
 104   // Setup the track cuts
 105   //
 106
 107
 108   //ESD quality cuts DielectronTrackCuts
 109   if (!isAOD) {
 110     diele->GetTrackFilter().AddCuts(SetupESDtrackCutsDieleData(cutDefinition));
 111   } else {
 112     AliDielectronTrackCuts *trackCuts=new AliDielectronTrackCuts("trackCuts","trackCuts");
 113     trackCuts->SetRequireTPCRefit(kTRUE);
 114     trackCuts->SetRequireITSRefit(kTRUE);
 115     diele->GetTrackFilter().AddCuts(trackCuts);
 116
 117   }
 118
 119
 120
 121   //Pt cut ----------------------------------------------------------
 122   AliDielectronVarCuts *pt = new AliDielectronVarCuts("ptCut","pt cut");
 123   pt->AddCut(AliDielectronVarManager::kPt,1.,100.);
 124   pt->AddCut(AliDielectronVarManager::kKinkIndex0,0.);
 125   //AOD additions since there are no AliESDtrackCuts -----------------
 126   pt->AddCut(AliDielectronVarManager::kEta,-0.9,0.9);
 127   pt->AddCut(AliDielectronVarManager::kImpactParXY, -1.0,   1.0);
 128   pt->AddCut(AliDielectronVarManager::kImpactParZ,  -3.0,   3.0);
 129   pt->AddCut(AliDielectronVarManager::AliDielectronVarManager::kTPCnSigmaEle,-2,3.);
 130   pt->AddCut(AliDielectronVarManager::kNclsTPC,80.,160.);
 131
 132
 133   diele->GetTrackFilter().AddCuts(pt);
 134 }
 135
 136 //______________________________________________________________________________________
 137 void SetupPairCutsDieleData(AliDielectron *diele, Int_t cutDefinition, Bool_t isAOD)
 138 {
 139   //
 140   // Setup the pair cuts
 141   //
 142
 143 AliDielectronVarCuts *mycut = new AliDielectronVarCuts("ptCutEMCAL","cut for EMCal");
 144 mycut->AddCut(AliDielectronVarManager::kEMCALE,5.,100.);
 145 mycut->AddCut(AliDielectronVarManager::kEMCALEoverP,0.75,1.25);
 146 AliDielectronPairLegCuts *varpair=new AliDielectronPairLegCuts();
 147 varpair->GetLeg1Filter().AddCuts(mycut);
 148 varpair->GetLeg2Filter().AddCuts(mycut);
 149 varpair->SetCutType(AliDielectronPairLegCuts::kAnyLeg);
 150 diele->GetPairFilter().AddCuts(varpair);
 151
 152 }
 153
 154
 155 //______________________________________________________________________________________
 156 void InitHistogramsDieleData(AliDielectron *diele, Int_t cutDefinition, Bool_t isAOD)
 157 {
 158   //
 159   // Initialise the histograms
 160   //
 161
 162   //Setup histogram Manager
 163   AliDielectronHistos *histos=new AliDielectronHistos(diele->GetName(),diele->GetTitle());
 164
 165   //Initialise histogram classes
 166   histos->SetReservedWords("Track;Pair");
 167
 168   //Track classes
 169   //to fill also track info from 2nd event loop until 2
 170   for (Int_t i=0; i<2; ++i){
 171     histos->AddClass(Form("Track_%s",AliDielectron::TrackClassName(i)));
 172   }
 173
 174   //Pair classes
 175   // to fill also mixed event histograms loop until 10
 176    if (cutDefinition==2) {
 177      for (Int_t i=0; i<10; ++i){
 178         histos->AddClass(Form("Pair_%s",AliDielectron::PairClassName(i)));
 179      }
 180    }
 181
 182   //legs from pair
 183   for (Int_t i=0; i<3; ++i){
 184     histos->AddClass(Form("Track_Legs_%s",AliDielectron::PairClassName(i)));
 185   }
 186
 187   //track rotation
 188   if (cutDefinition==2) {
 189       histos->AddClass(Form("Pair_%s",PairClassName(AliDielectron::kEv1PMRot)));
 190       histos->AddClass(Form("Track_Legs_%s",PairClassName(AliDielectron::kEv1PMRot)));
 191   }
 192
 193   //add histograms to event class
 194   if (cutDefinition==0) {
 195     histos->AddClass("Event");
 196     histos->UserHistogram("Event","VtxZ","Vertex Z;Z[cm]",300,-15.,15.,AliDielectronVarManager::kZvPrim);
 197     histos->UserHistogram("Event","Centrality","Centrality;Cent(%)",100,0.,100.,AliDielectronVarManager::kCentrality);
 198     histos->UserHistogram("Event","Multiplicity","Multiplicity V0;Multiplicity V0",500,0.,25000.,AliDielectronVarManager::kMultV0);
 199     histos->UserHistogram("Event","Multiplicity_cent","Multiplicity V0 x Cent;Cent(%);Multiplicity V0",10,0,100.,500,0.,25000.,AliDielectronVarManager::kCentrality,AliDielectronVarManager::kMultV0);
 200   }
 201
 202    //add histograms to Track classes
 203   histos->UserHistogram("Track","Pt","Pt;Pt [GeV];#tracks",100,0,20.,AliDielectronVarManager::kPt);
 204   histos->UserHistogram("Track","TPCnCls","Number of Clusters TPC;TPC number clusteres;#tracks",160,-0.5,159.5,AliDielectronVarManager::kNclsTPC);
 205
 206   histos->UserHistogram("Track","dXY","dXY;dXY [cm];#tracks",200,-1.,1.,AliDielectronVarManager::kImpactParXY);
 207   histos->UserHistogram("Track","dZ","dZ;dZ [cm];#tracks",200,-3.,3.,AliDielectronVarManager::kImpactParZ);
 208   histos->UserHistogram("Track","Eta_Phi","Eta Phi Map; Eta; Phi;#tracks",
 209                         100,-1,1,144,0,6.285,AliDielectronVarManager::kEta,AliDielectronVarManager::kPhi);
 210
 211   histos->UserHistogram("Track","dEdx_P","dEdx;P [GeV];TPC signal (arb units);#tracks",
 212                         200,0.2,20.,800,20.,200.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTPCsignal,kTRUE);
 213
 214   histos->UserHistogram("Track","dEdx_Pt","dEdx;Pt [GeV];TPC signal (arb units);#tracks",
 215                         200,0.2,20.,800,20.,200.,AliDielectronVarManager::kPt,AliDielectronVarManager::kTPCsignal,kTRUE);
 216
 217   histos->UserHistogram("Track","TPCnSigmaPion_P","TPCnSigmaPion;P [GeV];TPCnSigmaEle;#tracks",
 218                         200,0.2,20.,800,-12.,12.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTPCnSigmaPio,kTRUE);
 219
 220   histos->UserHistogram("Track","TPCnSigmaEle_P","TPCnSigmaEle;P [GeV];TPCnSigmaEle;#tracks",
 221                         200,0.2,20.,800,-12.,12.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTPCnSigmaEle,kTRUE);
 222
 223   histos->UserHistogram("Track","TPCnSigmaEle_Pt","TPCnSigmaEle;Pt [GeV];TPCnSigmaEle;#tracks",
 224                         200,0.2,20.,800,-12.,12.,AliDielectronVarManager::kPt,AliDielectronVarManager::kTPCnSigmaEle,kTRUE);
 225
 226   histos->UserHistogram("Track","TPCnSigmaEle_Phi","TPCnSigmaEle;#phi [rad];TPCnSigmaEle;#tracks",
 227                         200,0.,2*TMath::Pi(),800,-12.,12.,AliDielectronVarManager::kPhi,AliDielectronVarManager::kTPCnSigmaEle,kTRUE);
 228
 229   histos->UserHistogram("Track","TPCnSigmaEle_Eta","TPCnSigmaEle;#eta;TPCnSigmaEle;#tracks",
 230                         200,-1.,1.,800,-12.,12.,AliDielectronVarManager::kEta,AliDielectronVarManager::kTPCnSigmaEle,kTRUE);
 231
 232   histos->UserHistogram("Track","kNclsSTPC_kDeltaEta","kNclsSTPC vs kDeltaEta;kNclsSTPC;kDeltaEta;#tracks",
 233                         160,0.,159.,160,0.,1.6,AliDielectronVarManager::kNclsSTPC,AliDielectronVarManager::kDeltaEta,kTRUE);
 234   histos->UserHistogram("Track","kNclsSTPC_kDeltaPhi","kNclsSTPC vs kDeltaPhi;kNclsSTPC;kDeltaPhi;#tracks",
 235                         160,0.,159.,160,0.,6.4,AliDielectronVarManager::kNclsSTPC,AliDielectronVarManager::kDeltaPhi,kTRUE);
 236   histos->UserHistogram("Track","kDeltaPhi_kDeltaEta","kDeltaPhi vs kDeltaEta;kDeltaPhi;kDeltaEta;#tracks",
 237                         160,0.,6.4,160,0.,1.6,AliDielectronVarManager::kDeltaPhi,AliDielectronVarManager::kDeltaEta,kTRUE);
 238
 239   histos->UserHistogram("Track","dEdx_Phi","dEdx;#phi [rad];TPC signal (arb units);#tracks",
 240                         200,0.,2*TMath::Pi(),800,20.,200.,AliDielectronVarManager::kPhi,AliDielectronVarManager::kTPCsignal,kTRUE);
 241
 242   histos->UserHistogram("Track","dEdx_Eta","dEdx;#eta;TPC signal (arb units);#tracks",
 243                         200,-1.,1.,800,20.,200.,AliDielectronVarManager::kEta,AliDielectronVarManager::kTPCsignal,kTRUE);
 244
 245
 246
 247   histos->UserHistogram("Track","dEdx_nSigmaEMCal","dEdx;NsigmaEmcal;TPC signal (arb units);NSigmaEMCAL",
 248                         200,-5.,5.,800,20.,200.,AliDielectronVarManager::kEMCALnSigmaEle,AliDielectronVarManager::kTPCsignal,kTRUE);
 249
 250   histos->UserHistogram("Track","dEdx_TPCnSigmaEle","dEdx;TPC signal (arbunits);TPC number of sigmas Electrons;TPC signal (a.u.);#tracks",
 251                         100,-10.,10.,800,20.,200.,AliDielectronVarManager::kTPCnSigmaEle,AliDielectronVarManager::kTPCsignal,kTRUE);
 252
 253   histos->UserHistogram("Track","dEdx_EoverP","dEdx;EoverP;TPC signal (arbunits);E/P",100,0.,5.,800,20.,200.,AliDielectronVarManager::kEMCALEoverP,AliDielectronVarManager::kTPCsignal,kTRUE);
 254
 255   histos->UserHistogram("Track","nSigmaEMCal_EoverP","NsigmaEmcal;EoverP;NSigmaEMCAL;E/P",100,0.,5.,200,-5.,5.,AliDielectronVarManager::kEMCALEoverP,AliDielectronVarManager::kEMCALnSigmaEle,kTRUE);
 256
 257   histos->UserHistogram("Track","EMCal_E","EmcalE;Cluster Energy [GeV];#Clusters",200,0.,40.,AliDielectronVarManager::kEMCALE,kTRUE);
 258
 259   histos->UserHistogram("Track","EMCal_E_nCells","EmcalE_nCells;# cells;Cluster Energy [GeV]",200,0.,40.,AliDielectronVarManager::kEMCALNCells,AliDielectronVarManager::kEMCALE,kTRUE);
 260
 261   histos->UserHistogram("Track","ITS_FirstCls","ITS First Layer;ITS First Layer;#Entries",6,0.,6.,AliDielectronVarManager::kITSLayerFirstCls,kTRUE);
 262
 263   //add histograms to Pair classes
 264   histos->UserHistogram("Pair","InvMass","Inv.Mass;Inv. Mass [GeV];#pairs",
 265                         251,-.01,5.01,AliDielectronVarManager::kM);
 266   histos->UserHistogram("Pair","InvMass2D","Inv.Mass;Pt [GeV]; Inv. Mass [GeV]",
 267                         20,0.,20.,251,-.01,5.01,AliDielectronVarManager::kPt,AliDielectronVarManager::kM);
 268   histos->UserHistogram("Pair","Rapidity","Rapidity;Rapidity;#pairs",
 269                         100,-1.,1.,AliDielectronVarManager::kY);
 270   histos->UserHistogram("Pair","OpeningAngle","Opening angle;angle",
 271                         100,0.,3.15,AliDielectronVarManager::kOpeningAngle);
 272
 273   diele->SetHistogramManager(histos);
 274 }
 275
 276 void InitCFDieleData(AliDielectron *diele, Int_t cutDefinition, Bool_t isAOD)
 277 {
 278   //
 279   // Setupd the CF Manager if needed
 280   //
 281
 282   AliDielectronCF *cf=new AliDielectronCF(diele->GetName(),diele->GetTitle());
 283
 284   //pair variables
 285   cf->AddVariable(AliDielectronVarManager::kPt,"0., 2., 4.0, 4.5, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0, 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 20.0,25.0,30.0,40.0,50.0");
 286   cf->AddVariable(AliDielectronVarManager::kM,750,0.,15.);
 287   cf->AddVariable(AliDielectronVarManager::kPairType,12,0,12);
 288 //  cf->AddVariable(AliDielectronVarManager::kOpeningAngle,"0.,0.1,0.2,0.4,0.6,0.8,1.0,2.0,3.14");
 289 //  cf->AddVariable(AliDielectronVarManager::kY,20,-1.,1.);
 290   //leg variables
 291   cf->AddVariable(AliDielectronVarManager::kPt,"1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2., 3., 5., 10., 15., 20., 30., 40., 50.",kTRUE);
 292 //   cf->AddVariable(AliDielectronVarManager::kP,"0.0, 1.0, 1.2, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 5.0, 10.0, 100.0",kTRUE);
 293   cf->AddVariable(AliDielectronVarManager::kNclsTPC,"80, 85, 90, 95, 100, 110, 120, 160",kTRUE);
 294   cf->AddVariable(AliDielectronVarManager::kEta,"-5.0,-1.0,-0.9,-0.7,0.7,0.9,1.0,5.0",kTRUE);
 295    cf->AddVariable(AliDielectronVarManager::kPhi,"0.0,1.0,2.0,3.0,4.0,5.0,6.0",kTRUE);
 296 //  cf->AddVariable(AliDielectronVarManager::kEMCALE,"1.0,2.0, 3.0, 4.0, 4.5, 5.0, 5.5, 6.0, 7.0, 8.0, 9.0, 10.0,11.0,12.0,16.0,25.0",kTRUE);
 297 //  cf->AddVariable(AliDielectronVarManager::kEMCALnSigmaEle,"-3.5,-3.0,-2.0,-1.0,1.0,2.0,3.0,4.0,5.0",kTRUE);
 298 //  cf->AddVariable(AliDielectronVarManager::kEMCALNCells,25,0,25,kTRUE);
 299 //  cf->AddVariable(AliDielectronVarManager::kEMCALEoverP,"0.6, 0.65,0.7,0.8,0.9,1.1,1.2,1.3,1.4,1.8,2.0,2.1",kTRUE);
 300     cf->AddVariable(AliDielectronVarManager::kNclsSTPC,20,0.,1.,kTRUE);//shared cluster
 301     //ITS
 302   cf->AddVariable(AliDielectronVarManager::kNclsITS,6,0.,6.,kTRUE);
 303
 304   cf->AddVariable(AliDielectronVarManager::kTPCnSigmaEle,60,-3.,3.,kTRUE);
 305   cf->AddVariable(AliDielectronVarManager::kTPCnSigmaPio,60,-3.,3.,kTRUE);
 306   //event variables
 307   cf->AddVariable(AliDielectronVarManager::kCentrality,20,0.,100.);
 308 //  cf->AddVariable(AliDielectronVarManager::kMultV0,26,0.,26000.);
 309 //  cf->AddVariable(AliDielectronVarManager::kZvPrim,20,-10.,10.);
 310
 311 //  cf->AddVariable(AliDielectronVarManager::kITSLayerFirstCls,6,0.,6.,kTRUE);
 312
 313   if (hasMC){
 314       cf->AddVariable(AliDielectronVarManager::kPdgCode,10000,-5000.5,4999.5,kTRUE);
 315       cf->AddVariable(AliDielectronVarManager::kPdgCodeMother,10000,-5000.5,4999.5,kTRUE);
 316       cf->AddVariable(AliDielectronVarManager::kPdgCodeGrandMother,10000,-5000.5,4999.5,kTRUE);
 317   }
 318     //only in this case write MC truth info
 319   if (cutDefinition==0){
 320     cf->SetStepForMCtruth();
 321 //     cf->SetStepsForMCtruthOnly();
 322  }
 323
 324   diele->SetCFManagerPair(cf);
 325
 326 }
 327
 328 void AddMCSignals(AliDielectron *diele){
 329   //Do we have an MC handler?
 330   //Bool_t hasMC=(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()!=0x0);
 331   if (!hasMC) return;
 332
 333   AliDielectronSignalMC* inclusiveJpsi = new AliDielectronSignalMC("inclusiveJpsi","Inclusive J/psi");
 334   inclusiveJpsi->SetLegPDGs(11,-11);
 335   inclusiveJpsi->SetMotherPDGs(443,443);
 336   inclusiveJpsi->SetMothersRelation(AliDielectronSignalMC::kSame);
 337   inclusiveJpsi->SetFillPureMCStep(kTRUE);
 338   inclusiveJpsi->SetCheckBothChargesLegs(kTRUE,kTRUE);
 339   inclusiveJpsi->SetCheckBothChargesMothers(kTRUE,kTRUE);
 340   diele->AddSignalMC(inclusiveJpsi);
 341
 342   AliDielectronSignalMC* promptJpsi = new AliDielectronSignalMC("promptJpsi","Prompt J/psi");   // prompt J/psi (not from beauty decays)
 343   promptJpsi->SetLegPDGs(11,-11);
 344   promptJpsi->SetMotherPDGs(443,443);
 345   promptJpsi->SetGrandMotherPDGs(503,503,kTRUE,kTRUE);   // not from beauty hadrons
 346   promptJpsi->SetMothersRelation(AliDielectronSignalMC::kSame);
 347   promptJpsi->SetFillPureMCStep(kTRUE);
 348   promptJpsi->SetLegSources(AliDielectronSignalMC::kFinalState, AliDielectronSignalMC::kFinalState);
 349   promptJpsi->SetCheckBothChargesLegs(kTRUE,kTRUE);
 350   promptJpsi->SetCheckBothChargesMothers(kTRUE,kTRUE);
 351   promptJpsi->SetCheckBothChargesGrandMothers(kTRUE,kTRUE);
 352   diele->AddSignalMC(promptJpsi);
 353
 354   AliDielectronSignalMC* beautyJpsi = new AliDielectronSignalMC("beautyJpsi","Beauty J/psi");
 355   beautyJpsi->SetLegPDGs(11,-11);
 356   beautyJpsi->SetMotherPDGs(443,443);
 357   beautyJpsi->SetMothersRelation(AliDielectronSignalMC::kSame);
 358   beautyJpsi->SetGrandMotherPDGs(500,500);
 359   beautyJpsi->SetFillPureMCStep(kTRUE);
 360   beautyJpsi->SetCheckBothChargesLegs(kTRUE,kTRUE);
 361   beautyJpsi->SetCheckBothChargesMothers(kTRUE,kTRUE);
 362   beautyJpsi->SetCheckBothChargesGrandMothers(kTRUE,kTRUE);
 363   diele->AddSignalMC(beautyJpsi);
 364
 365   AliDielectronSignalMC* directJpsi = new AliDielectronSignalMC("directJpsi","Direct J/psi");   // embedded J/psi
 366   directJpsi->SetLegPDGs(11,-11);
 367   directJpsi->SetMotherPDGs(443,443);
 368   directJpsi->SetMothersRelation(AliDielectronSignalMC::kSame);
 369   directJpsi->SetFillPureMCStep(kTRUE);
 370   directJpsi->SetLegSources(AliDielectronSignalMC::kFinalState, AliDielectronSignalMC::kFinalState);
 371   directJpsi->SetMotherSources(AliDielectronSignalMC::kDirect, AliDielectronSignalMC::kDirect);
 372   directJpsi->SetCheckBothChargesLegs(kTRUE,kTRUE);
 373   directJpsi->SetCheckBothChargesMothers(kTRUE,kTRUE);
 374   diele->AddSignalMC(directJpsi);
 375 }
 376
 377
 378
 379
 380 /*
 381 //______________________________________________________________________________________
 382 AliESDtrackCuts *SetupESDtrackCutsDieleData(Int_t cutDefinition)
 383 {
 384   //
 385   // Setup default AliESDtrackCuts
 386   //
 387   AliESDtrackCuts *esdTrackCuts = new AliESDtrackCuts;
 388
 389   // basic track quality cuts  (basicQ)
 390   esdTrackCuts->SetMaxDCAToVertexZ(3.0);
 391   esdTrackCuts->SetMaxDCAToVertexXY(1.0);
 392   esdTrackCuts->SetEtaRange( -0.9 , 0.9 );
 393   esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
 394   esdTrackCuts->SetRequireITSRefit(kTRUE);
 395   esdTrackCuts->SetRequireTPCRefit(kTRUE);
 396   esdTrackCuts->SetPtRange(1.,1e30);
 397   esdTrackCuts->SetMinNClustersTPC(70);
 398   esdTrackCuts->SetMaxChi2PerClusterTPC(4);
 399   esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
 400
 401   return esdTrackCuts;
 402 }
 403 */