Add missed macro

diff --git a/PWGDQ/dielectron/macrosLMEE/AddTask_cbaumann_LMEEnoPID.C b/PWGDQ/dielectron/macrosLMEE/AddTask_cbaumann_LMEEnoPID.C
index 8348af7a8ce9791a8e3788143ec4732676538c32..cdf5d1c2512cd85c4ca64c9653290162eb260771 100644 (file)
--- a/PWGDQ/dielectron/macrosLMEE/AddTask_cbaumann_LMEEnoPID.C
+++ b/PWGDQ/dielectron/macrosLMEE/AddTask_cbaumann_LMEEnoPID.C
@@ -55,6 +55,10 @@ AliAnalysisTask *AddTask_cbaumann_LMEEnoPID(Bool_t withMC = kFALSE,Bool_t enable
   AliDielectron *lowmass1=ConfigLMEEnoPID(1,hasMC,enableCF);
   task->AddDielectron(lowmass1);
   printf("add: %s\n",lowmass1->GetName());
+
+   AliDielectron *lowmass2=ConfigLMEEnoPID(2,hasMC,enableCF);
+  task->AddDielectron(lowmass2);
+  printf("add: %s\n",lowmass2->GetName());
    
   mgr->AddTask(task);
 

diff --git a/PWGDQ/dielectron/macrosLMEE/ConfigLMEEnoPID.C b/PWGDQ/dielectron/macrosLMEE/ConfigLMEEnoPID.C
new file mode 100644 (file)
index 0000000..ef56e97
--- /dev/null
+++ b/PWGDQ/dielectron/macrosLMEE/ConfigLMEEnoPID.C
@@ -0,0 +1,296 @@
+#include "LMEECutLib.C"
+void InitHistograms(AliDielectron *die, Int_t cutDefinition);
+void InitCF(AliDielectron* die, Int_t cutDefinition);
+void EnableMC();
+
+TString names=("noPairingNoPID;noPairingTOFonly,noPairingTPCTOFPbPb2011");
+TObjArray *arrNames=names.Tokenize(";");
+const Int_t nDie=arrNames->GetEntries();
+
+Bool_t MCenabled=kFALSE;
+
+
+AliDielectron* ConfigLMEEnoPID(Int_t cutDefinition, Bool_t hasMC=kFALSE,Bool_t CFenabled=kFALSE)
+{
+
+  Int_t selectedPID=-1;
+  Bool_t rejectionStep=kFALSE;
+  LMEECutLib*  LMCL = new LMEECutLib();
+
+  //
+  // Setup the instance of AliDielectron
+  //
+
+  MCenabled=hasMC;
+
+  // create the actual framework object
+
+  TString name=Form("%02d",cutDefinition);
+  if ((cutDefinition)<arrNames->GetEntriesFast()){
+       name=arrNames->At((cutDefinition))->GetName();
+  }
+
+  //thisCut only relevant for MC:
+  AliDielectron *die =
+       new AliDielectron(Form
+               ("%s",name.Data()),
+               Form("Track cuts: %s",name.Data()));
+
+
+  //Setup AnalysisSelection:
+  //Choose Arbitratry cuts, Basic Track Cuts shoudl be the same for all
+  if (cutDefinition==0) {
+       //not yet implemented
+       selectedPID = LMEECutLib::kpp2010TPCandTOF;
+       rejectionStep = kFALSE;
+  }
+  else if (cutDefinition==1) {
+   //TPCTOFCentnoRej =>
+       selectedPID = LMEECutLib::kpp2010TPCandTOF;
+       rejectionStep = kFALSE;
+  }
+  else if (cutDefinition==2) {
+   //TPCTOFCentnoRej =>
+       selectedPID = LMEECutLib::kPbPb2011TPCandTOF;
+       rejectionStep = kFALSE;
+  }
+  else {
+       cout << " =============================== " << endl;
+       cout << " ==== INVALID CONFIGURATION ==== " << endl;
+       cout << " =============================== " << endl;
+  }
+
+
+  //Now configure task
+
+  //Apply correct Pre-Filter Scheme, if necessary
+  die->SetPreFilterAllSigns();
+  
+   //SWITCH OFF PAIRING FOR Track-only analysis
+  die->SetNoPairing();
+
+  die->GetTrackFilter().AddCuts( LMCL->GetTrackCutsAna(selectedPID) );
+
+  if (cutDefinition==1) {
+         AliDielectronPID *pidTOFonly = new AliDielectronPID("TOFonly","TOFonly");
+         pidTOFonly->AddCut(AliDielectronPID::kTOF ,AliPID::kElectron , -3. , 3. , 0.0 , 100., kFALSE );
+         die->GetTrackFilter().AddCuts(pidTOFonly);
+       }
+
+  if (cutDefinition==2) {
+         die->GetTrackFilter().AddCuts(LMCL->GetPIDCutsAna(selectedPID));
+       }
+
+
+
+  // histogram setup
+  // only if an AliDielectronHistos object is attached to the
+  // dielectron framework histograms will be filled
+  //
+  InitHistograms(die,cutDefinition);
+
+  // the last definition uses no cuts and only the QA histograms should be filled!
+  if (CFenabled) InitCF(die,cutDefinition);
+
+  return die;
+}
+
+//______________________________________________________________________________________
+
+void InitHistograms(AliDielectron *die, Int_t cutDefinition)
+{
+  //
+  // Initialise the histograms
+  //
+
+  //Setup histogram Manager
+  AliDielectronHistos *histos=
+       new AliDielectronHistos(die->GetName(),
+               die->GetTitle());
+  //Initialise histogram classes
+  histos->SetReservedWords("Track;Pair;Pre;RejTrack;RejPair");
+
+  //Event class
+//  if (cutDefinition==nDie-1) 
+               histos->AddClass("Event");
+
+  //Track classes
+  //to fill also track info from 2nd event loop until 2
+  for (Int_t i=0; i<2; ++i){
+       histos->AddClass(Form("Track_%s",AliDielectron::TrackClassName(i)));
+  }
+
+  //Pair classes
+  // to fill also mixed event histograms loop until 10
+  for (Int_t i=0; i<3; ++i){
+       histos->AddClass(Form("Pair_%s",AliDielectron::PairClassName(i)));
+  }
+
+  //ME and track rot
+  if (die->GetMixingHandler()) {
+       histos->AddClass(Form("Pair_%s",AliDielectron::PairClassName(3)));
+       histos->AddClass(Form("Pair_%s",AliDielectron::PairClassName(4)));
+       histos->AddClass(Form("Pair_%s",AliDielectron::PairClassName(6)));
+       histos->AddClass(Form("Pair_%s",AliDielectron::PairClassName(7)));
+  }
+  if (die->GetTrackRotator()) {
+       histos->AddClass(Form("Pair_%s",AliDielectron::PairClassName(10)));
+  }
+
+  //PreFilter Classes
+  //to fill also track info from 2nd event loop until 2
+  for (Int_t i=0; i<2; ++i){
+       histos->AddClass(Form("Pre_%s",AliDielectron::TrackClassName(i)));
+  }
+
+
+  //Create Classes for Rejected Tracks/Pairs:
+  for (Int_t i=0; i<2; ++i){
+       histos->AddClass(Form("RejTrack_%s",AliDielectron::TrackClassName(i)));
+  }
+  for (Int_t i=0; i<3; ++i){
+       histos->AddClass(Form("RejPair_%s",AliDielectron::PairClassName(i)));
+  }
+
+  /*
+  //track rotation
+
+  histos->AddClass(Form("Pair_%s",AliDielectron::PairClassName(AliDielectron::kEv1PMRot)));
+  histos->AddClass(Form("Track_Legs_%s",AliDielectron::PairClassName(AliDielectron::kEv1PMRot)));
+  */
+       //add histograms to event class
+       histos->UserHistogram("Event","nEvents","Number of processed events after cuts;Number events",
+               1,0.,1.,AliDielectronVarManager::kNevents);
+       histos->UserHistogram("Event","Centrality","Centrality;Centrality [%]","0,10,20,40,80,100,101",
+               AliDielectronVarManager::kCentrality);
+
+
+  //add histograms to Track classes, also fills RejTrack
+  histos->UserHistogram("Track","Pt","Pt;Pt [GeV];#tracks",200,0,20.,AliDielectronVarManager::kPt);
+  histos->UserHistogram("Track","NclsSFracTPC","NclsSFracTPC; NclsSFracTPC;#tracks",200,0,10.,AliDielectronVarManager::kNclsSFracTPC);
+  histos->UserHistogram("Track","TPCclsDiff","TPCclsDiff; TPCclsDiff;#tracks",200,0,10.,AliDielectronVarManager::kTPCclsDiff);
+
+  histos->UserHistogram("Track","ITS_dEdx_P","ITS_dEdx;P [GeV];ITS signal (arb units);#tracks",
+         400,0.0,20.,200,0.,1000.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kITSsignal,kTRUE);
+
+  histos->UserHistogram("Track","dEdx_P","dEdx;P [GeV];TPC signal (arb units);#tracks",
+         400,0.0,20.,200,0.,200.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTPCsignal,kTRUE);
+/*
+  histos->UserHistogram("Track","TPCnSigmaEle_P","TPC number of sigmas Electrons;P [GeV];TPC number of sigmas Electrons;#tracks",
+         400,0.0,20.,100,-5.,5.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTPCnSigmaEle,kTRUE);
+  histos->UserHistogram("Track","TPCnSigmaKao_P","TPC number of sigmas Kaons;P [GeV];TPC number of sigmas Kaons;#tracks",
+         400,0.0,20.,100,-5.,5.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTPCnSigmaKao,kTRUE);
+  histos->UserHistogram("Track","TPCnSigmaPio_P","TPC number of sigmas Pions;P [GeV];TPC number of sigmas Pions;#tracks",
+         400,0.0,20.,100,-5.,5.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTPCnSigmaPio,kTRUE);
+
+  histos->UserHistogram("Track","TRDpidPobEle_P","TRD PID probability Electrons;P [GeV];TRD prob Electrons;#tracks",
+         400,0.0,20.,100,0.,1.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTRDprobEle,kTRUE);
+  histos->UserHistogram("Track","TRDpidPobPio_P","TRD PID probability Pions;P [GeV];TRD prob Pions;#tracks",
+         400,0.0,20.,100,0.,1.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTRDprobPio,kTRUE);
+
+  histos->UserHistogram("Track","TOFnSigmaKao_P","TOF number of sigmas Kaons;P [GeV];TOF number of sigmas Kaons;#tracks",
+         400,0.0,20.,100,-5.,5.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTOFnSigmaKao,kTRUE);
+  histos->UserHistogram("Track","TOFnSigmaPro_P","TOF number of sigmas Protons;P [GeV];TOF number of sigmas Protons;#tracks",
+         400,0.0,20.,100,-5.,5.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTOFnSigmaPro,kTRUE);
+*/
+  histos->UserHistogram("Track","Eta_Phi","Eta Phi Map; Eta; Phi;#tracks",
+         100,-2,2,100,0,3.15,AliDielectronVarManager::kEta,AliDielectronVarManager::kPhi);
+
+  histos->UserHistogram("Track","dXY","dXY;dXY [cm];#tracks",200,-2.,2.,AliDielectronVarManager::kImpactParXY);
+
+  histos->UserHistogram("Track","TPCnCls","Number of Clusters TPC;TPC number clusteres;#tracks",159,0.,159.,AliDielectronVarManager::kNclsTPC);
+
+  histos->UserHistogram("Track","TPCnCls_kNFclsTPCr","nTPC vs nTPCr;nTPC vs nTPCr;#tracks",159,0.,159.,159,0.,159.,AliDielectronVarManager::kNclsTPC,AliDielectronVarManager::kNFclsTPCr);
+
+  histos->UserHistogram("Track","kNFclsTPCr_pT","nTPCr vs pt;nTPCr vs pt;#tracks",159,0.,159.,200,0.,20.,AliDielectronVarManager::kNFclsTPCr,AliDielectronVarManager::kPt);
+
+  //add histograms to Pair classes, also fills RejPair
+  /*
+  histos->UserHistogram("Pair","InvMass","Inv.Mass;Inv. Mass [GeV];#pairs",
+         500,0.0,5.00,AliDielectronVarManager::kM);
+  histos->UserHistogram("Pair","Rapidity","Rapidity;Rapidity;#pairs",
+         100,-2.,2.,AliDielectronVarManager::kY);
+  histos->UserHistogram("Pair","OpeningAngle","Opening angle;angle",
+         100,0.,3.15,AliDielectronVarManager::kOpeningAngle);
+  //2D Histo Plot
+  histos->UserHistogram("Pair","InvMassPairPt","Inv.Mass vs PairPt;Inv. Mass [GeV], pT [GeV];#pairs",
+         500,0.0,5.0,100,0.,10.,AliDielectronVarManager::kM,AliDielectronVarManager::kPt);
+
+  histos->UserHistogram("Pair","InvMassOpeningAngle","Opening Angle vs Inv.Mass;Inv. Mass [GeV];#pairs",
+         500,0.0,5.0,200,0.,6.3,AliDielectronVarManager::kM,AliDielectronVarManager::kOpeningAngle);
+
+  //add histograms to PRE-Track classes
+  histos->UserHistogram("Pre","Pt","Pt;Pt [GeV];#tracks",200,0,20.,AliDielectronVarManager::kPt);
+
+  histos->UserHistogram("Pre","ITS_dEdx_P","ITS_dEdx;P [GeV];ITS signal (arb units);#tracks",
+         400,0.0,20.,200,0.,1000.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kITSsignal,kTRUE);
+
+  histos->UserHistogram("Pre","dEdx_P","dEdx;P [GeV];TPC signal (arb units);#tracks",
+         400,0.0,20.,200,0.,200.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTPCsignal,kTRUE);
+*/
+/*
+  histos->UserHistogram("Pre","TPCnSigmaEle_P","TPC number of sigmas Electrons;P [GeV];TPC number of sigmas Electrons;#tracks",
+         400,0.0,20.,100,-5.,5.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTPCnSigmaEle,kTRUE);
+  histos->UserHistogram("Pre","TPCnSigmaKao_P","TPC number of sigmas Kaons;P [GeV];TPC number of sigmas Kaons;#tracks",
+         400,0.0,20.,100,-5.,5.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTPCnSigmaKao,kTRUE);
+  histos->UserHistogram("Pre","TPCnSigmaPio_P","TPC number of sigmas Pions;P [GeV];TPC number of sigmas Pions;#tracks",
+         400,0.0,20.,100,-5.,5.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTPCnSigmaPio,kTRUE);
+
+  histos->UserHistogram("Pre","TRDpidPobEle_P","TRD PID probability Electrons;P [GeV];TRD prob Electrons;#tracks",
+         400,0.0,20.,100,0.,1.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTRDprobEle,kTRUE);
+  histos->UserHistogram("Pre","TRDpidPobPio_P","TRD PID probability Pions;P [GeV];TRD prob Pions;#tracks",
+         400,0.0,20.,100,0.,1.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTRDprobPio,kTRUE);
+
+  histos->UserHistogram("Pre","TOFnSigmaKao_P","TOF number of sigmas Kaons;P [GeV];TOF number of sigmas Kaons;#tracks",
+         400,0.0,20.,100,-5.,5.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTOFnSigmaKao,kTRUE);
+  histos->UserHistogram("Pre","TOFnSigmaPro_P","TOF number of sigmas Protons;P [GeV];TOF number of sigmas Protons;#tracks",
+         400,0.0,20.,100,-5.,5.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTOFnSigmaPro,kTRUE);
+*/
+  /*
+  histos->UserHistogram("Pre","Eta_Phi","Eta Phi Map; Eta; Phi;#tracks",
+         100,-2,2,100,0,3.15,AliDielectronVarManager::kEta,AliDielectronVarManager::kPhi);
+
+  histos->UserHistogram("Pre","dXY","dXY;dXY [cm];#tracks",200,-2.,2.,AliDielectronVarManager::kImpactParXY);
+
+  histos->UserHistogram("Pre","TPCnCls","Number of Clusters TPC;TPC number clusteres;#tracks",159,0.,159.,AliDielectronVarManager::kNclsTPC);
+*/
+  die->SetHistogramManager(histos);
+}
+
+
+void InitCF(AliDielectron* die, Int_t cutDefinition)
+{
+  //
+  // Setupd the CF Manager if needed
+  //
+  AliDielectronCF *cf=new AliDielectronCF(die->GetName(),die->GetTitle());
+
+  //pair variables
+/* cf->AddVariable(AliDielectronVarManager::kP,200,0,20);
+  cf->AddVariable(AliDielectronVarManager::kM,201,-0.01,4.01); //20Mev Steps
+*/
+   cf->AddVariable(AliDielectronVarManager::kPairType,10,0,10);
+/*
+  //leg variables
+  cf->AddVariable(AliDielectronVarManager::kP,200,0.,20.,kTRUE);
+    cf->AddVariable(AliDielectronVarManager::kITSsignal,1000,0.0.,1000.,kTRUE);
+  cf->AddVariable(AliDielectronVarManager::kTPCsignal,500,0.0.,500.,kTRUE);
+*/
+  //only in this case write MC truth info
+  if (MCenabled) {
+       cf->SetStepForMCtruth();
+       cf->SetStepsForMCtruthOnly();
+       cf->AddVariable(AliDielectronVarManager::kHaveSameMother,21,-10,10,kTRUE);
+       cf->AddVariable(AliDielectronVarManager::kPdgCode,10000,-5000.5,4999.5,kTRUE);
+       cf->AddVariable(AliDielectronVarManager::kPdgCodeMother,10000,-5000.5,4999.5,kTRUE);
+  }
+
+  cf->SetStepsForSignal();
+  die->SetCFManagerPair(cf);
+}
+
+//--------------------------------------
+void EnableMC() {
+  MCenabled=kTRUE;
+}
+