//set config file name
TString configFile("$TRAIN_ROOT/jbook_jpsi/ConfigJpsi_jb_PbPb.C");
TString trainRoot=gSystem->Getenv("TRAIN_ROOT");
- if (trainRoot.IsNull()) configFile="../PWGDQ/dielectron/macrosJPSI/ConfigJpsi_jb_PbPb.C";
+ if (trainRoot.IsNull()) configFile="$ALICE_ROOT/PWGDQ/dielectron/macrosJPSI/ConfigJpsi_jb_PbPb.C";
//create task and add it to the manager
--- /dev/null
+AliAnalysisTask *AddTask_jbook_v2(){
+ //get the current analysis manager
+ AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
+ if (!mgr) {
+ Error("AddTask_jbook_v2", "No analysis manager found.");
+ return 0;
+ }
+
+ //Do we have an MC handler?
+ Bool_t hasMC=(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()!=0x0);
+ //Do we have an AOD handler?
+ Bool_t isAOD=(mgr->GetInputEventHandler()->IsA()==AliAODInputHandler::Class() ? kTRUE : kFALSE);
+
+ //set config file name
+ TString configFile("$TRAIN_ROOT/jbook_jpsi/ConfigJpsi_jb_PbPb.C");
+ TString trainRoot=gSystem->Getenv("TRAIN_ROOT");
+ if (trainRoot.IsNull()) configFile="$ALICE_ROOT/PWGDQ/dielectron/macrosJPSI/ConfigJpsi_jb_PbPb.C";
+
+
+ //create task and add it to the manager
+ AliAnalysisTaskMultiDielectron *task=new AliAnalysisTaskMultiDielectron("MultiDieData");
+ task->SetTriggerMask(AliVEvent::kMB+AliVEvent::kCentral+AliVEvent::kSemiCentral);
+ if (!hasMC) task->UsePhysicsSelection();
+ mgr->AddTask(task);
+
+ //load dielectron configuration file
+ TString checkconfig="ConfigJpsi_jb_PbPb";
+ if (!gROOT->GetListOfGlobalFunctions()->FindObject(checkconfig.Data()))
+ gROOT->LoadMacro(configFile.Data());
+
+ //add dielectron analysis with different cuts to the task
+ for (Int_t i=16; i<nDie; ++i){ //nDie defined in config file
+ AliDielectron *jpsi=ConfigJpsi_jb_PbPb(i);
+ if (jpsi ) task->AddDielectron(jpsi);
+ if (jpsi ) printf("add: %s\n",jpsi->GetName());
+ }
+
+ //Add event filter
+ AliDielectronEventCuts *eventCuts=new AliDielectronEventCuts("eventCuts","Vertex Track && |vtxZ|<10 && ncontrib>0");
+ if (isAOD) eventCuts->SetVertexType(AliDielectronEventCuts::kVtxAny);
+ eventCuts->SetRequireVertex();
+ eventCuts->SetMinVtxContributors(1);
+ eventCuts->SetVertexZ(-10.,10.);
+ eventCuts->SetCentralityRange(0.0,80.0);
+ task->SetEventFilter(eventCuts);
+
+
+ //create output container
+ AliAnalysisDataContainer *coutput1 =
+ mgr->CreateContainer("jbook_v2_tree",
+ TTree::Class(),
+ AliAnalysisManager::kExchangeContainer,
+ "jbook_v2_default");
+
+ AliAnalysisDataContainer *cOutputHist1 =
+ mgr->CreateContainer("jbook_v2_QA",
+ TList::Class(),
+ AliAnalysisManager::kOutputContainer,
+ "jbookv2.root");
+
+ AliAnalysisDataContainer *cOutputHist2 =
+ mgr->CreateContainer("jbook_v2_CF",
+ TList::Class(),
+ AliAnalysisManager::kOutputContainer,
+ "jbookv2.root");
+
+ AliAnalysisDataContainer *cOutputHist3 =
+ mgr->CreateContainer("jbook_v2_EventStat",
+ TH1D::Class(),
+ AliAnalysisManager::kOutputContainer,
+ "jbookv2.root");
+
+ mgr->ConnectInput(task, 0, mgr->GetCommonInputContainer());
+ mgr->ConnectOutput(task, 0, coutput1 );
+ mgr->ConnectOutput(task, 1, cOutputHist1);
+ mgr->ConnectOutput(task, 2, cOutputHist2);
+ mgr->ConnectOutput(task, 3, cOutputHist3);
+
+ return task;
+}
//Do we have an MC handler?
Bool_t hasMC=(mgr->GetMCtruthEventHandler()!=0x0);
- TString configFile("ConfigJpsi_mf_pp.C");
- if (hasMC) configFile="../PWGDQ/dielectron/macros/ConfigJpsi2eeEff.C";
+ TString configFile("$ALICE_ROOT/PWGDQ/dielectron/macrosJPSI/ConfigJpsi_mf_pp.C");
+ if (hasMC) configFile="$ALICE_ROOT/PWGDQ/dielectron/macros/ConfigJpsi2eeEff.C";
Bool_t isAOD=mgr->GetInputEventHandler()->IsA()==AliAODInputHandler::Class();
//create task and add it to the manager
void SetEtaCorrection();
TVectorD *GetRunNumbers();
-AliESDtrackCuts *SetupESDtrackCuts(Int_t cutDefinition);
-
TString names=("TPC;TOF;TRD;rec;TOFTRD;TOFTRD2;ITScls;ITSamy;dca;chi2;Gam0;Gam01;Gam05;Gam10;Gam15;Gam20;EtaGap01;EtaGap02;EtaGap03;EtaGap04;EtaGap05;SubLS;SubRndm");
enum { kTPC=0, kTOF, kTRD, krec, kTOFTRD, kTOFTRD2, kITScls, kITSamy, kDCA, kChi, kGam0, kGam01, kGam05, kGam10, kGam15, kGam20, kEtaGap01, kEtaGap02, kEtaGap03, kEtaGap04, kEtaGap05, kSubLS, kSubRndm };
// trd tables
TString pidTab="$TRAIN_ROOT/util/dielectron/dielectron/TRDpidEff_eleProb07_TRDntr4_6.root";
TString trainRoot=gSystem->Getenv("TRAIN_ROOT");
- if (trainRoot.IsNull()) pidTab="../PWGDQ/dielectron/files/TRDpidEff_eleProb07_TRDntr4_6.root";
+ if (trainRoot.IsNull()) pidTab="$ALICE_ROOT/PWGDQ/dielectron/files/TRDpidEff_eleProb07_TRDntr4_6.root";
if (gSystem->AccessPathName(gSystem->ExpandPathName(pidTab.Data())))
Error("ConfigPbPb","PID table not found: %s",pidTab.Data());
}
AliDielectronVarCuts *poi = new AliDielectronVarCuts("PoI","PoI");
- poi->AddCut(AliDielectronVarManager::kM,2.92,3.16); // particles of interest, jpsi mass window
+ poi->AddCut(AliDielectronVarManager::kM,2.92,3.20); // particles of interest, jpsi mass window
die->GetEventPlanePOIPreFilter().AddCuts(poi);
- AliDielectronVarCuts *poiTrk = new AliDielectronVarCuts("PoItracks","PoItracks");
- poiTrk->AddCut(AliDielectronVarManager::kPt,0.8,1e30);
- die->GetEventPlanePreFilter().AddCuts(poiTrk);
-
if(cutDefinition >= kEtaGap01 &&
cutDefinition < kSubLS ) {
AliDielectronVarCuts *etaGap = new AliDielectronVarCuts(AliDielectronVarManager::GetValueName(AliDielectronVarManager::kEta),"etaGap");
// setup eta correction
if(isESD) SetEtaCorrection();
+ // VZERO calibration
+ TString trainRoot=gSystem->Getenv("TRAIN_ROOT");
+ if (!trainRoot.IsNull()) {
+ die->SetVZEROCalibrationFilename("$TRAIN_ROOT/util/dielectron/dielectron/VzeroCalibrationLHC10h.root");
+ die->SetVZERORecenteringFilename("$TRAIN_ROOT/util/dielectron/dielectron/VzeroRecenteringLHC10h.root");
+ }
+
return die;
}
die->GetTrackFilter().AddCuts(cuts);
//Pt cut, should make execution a bit faster
- AliDielectronVarCuts *pt = new AliDielectronVarCuts("Pt>.8","Pt>.8");
- pt->AddCut(AliDielectronVarManager::kPt,0.8,1e30);
+ AliDielectronVarCuts *pt = new AliDielectronVarCuts("PtCut","PtCut");
+ if(cutDefinition >= kEtaGap01 )
+ pt->AddCut(AliDielectronVarManager::kPt,1.1,1e30);
+ else
+ pt->AddCut(AliDielectronVarManager::kPt,0.8,1e30);
cuts->AddCut(pt);
- //ESD track cuts
- Bool_t isESD=(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()->IsA()==AliESDInputHandler::Class());
- if(isESD) cuts->AddCut(SetupESDtrackCuts(cutDefinition));
-
- // track cuts ESD and AOD
+ // track cuts ESD and AOD
AliDielectronVarCuts *varCuts = new AliDielectronVarCuts("VarCuts","VarCuts");
varCuts->AddCut(AliDielectronVarManager::kImpactParXY, -1.0, 1.0);
varCuts->AddCut(AliDielectronVarManager::kImpactParZ, -3.0, 3.0);
varCuts->AddCut(AliDielectronVarManager::kEta, -0.9, 0.9);
varCuts->AddCut(AliDielectronVarManager::kTPCchi2Cl, 0.0, 4.0);
varCuts->AddCut(AliDielectronVarManager::kNclsTPC, 70.0, 160.0);
+ varCuts->AddCut(AliDielectronVarManager::kKinkIndex0, 0.0);
switch(cutDefinition) {
case kTOFTRD2: varCuts->AddCut(AliDielectronVarManager::kITSLayerFirstCls,-0.01,0.5); //ITS(0) = SPDfirst
break;
////////////////////////////////// DATA
if(!hasMC) {
- pid->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-3.,3.); // merge with below ones vv
pid->AddCut(AliDielectronPID::kTPC,AliPID::kPion,-100.,3.5,0.,0.,kTRUE);
- //if((cutDefinition!=kTOF && cutDefinition<kTOFTRD) || cutDefinition>=kGam0)
pid->AddCut(AliDielectronPID::kTPC,AliPID::kProton,-100.,3.5,0.,0.,kTRUE);
if(cutDefinition==kTRD || cutDefinition>=kTOFTRD || cutDefinition>=kTOFTRD2)
if(hasMC) {
// electron
- pid->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-3.,3.); //merge with data only ^^
Double_t nSigmaPi = 3.5; Double_t nSigmaP = 3.5;
Double_t resolution=0.0549;
Double_t BBpro[5] = {0};
TString list=gSystem->Getenv("LIST");
//LHC11a10b
- if (list.Contains("LHC11a10b")) {
+ if (list.Contains("LHC11a10b") || list.IsNull()) {
printf("LHC11a10b parameters\n");
ffPro->SetParameters(BBpro[0],BBpro[1],BBpro[2],BBpro[3],BBpro[4]);
ffPio->SetParameters(BBpio[0],BBpio[1],BBpio[2],BBpio[3],BBpio[4]);
// shifts for the nSigma electrons
TGraph* nSigmaCorrection = new TGraph();
// LHC11a10b
- if (list.Contains("LHC11a10b")) {
+ if (list.Contains("LHC11a10b") || list.IsNull()) {
nSigmaCorrection->SetPoint(0, 137161., -0.50-(0.28));
nSigmaCorrection->SetPoint(1, 139510., -0.50-(0.28));
pid->SetCorrGraph(nSigmaCorrection);
} //hasMC
////////////////////////////////// DATA + MC
- // pid cuts TOF & TRD
+ // pid cuts TPC + TOF & TRD
+ pid->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-3.,3.);
if(cutDefinition==kTOF || cutDefinition>=kTOFTRD || cutDefinition>=kTOFTRD2)
pid->AddCut(AliDielectronPID::kTOF,AliPID::kElectron,-3,3.,0.,0.,kFALSE,AliDielectronPID::kIfAvailable);
}
-//______________________________________________________________________________________
-AliESDtrackCuts *SetupESDtrackCuts(Int_t cutDefinition)
-{
- //
- // Setup default AliESDtrackCuts
- //
- AliESDtrackCuts *esdTrackCuts = new AliESDtrackCuts;
- // reconstruction cuts
- esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
-
- /*
- // basic track quality cuts (basicQ)
- esdTrackCuts->SetMaxDCAToVertexZ(3.0);
- esdTrackCuts->SetMaxDCAToVertexXY(1.0);
-
- // acceptance cuts
- esdTrackCuts->SetEtaRange( -0.9 , 0.9 );
- //esdTrackCuts->SetPtRange( 0.8, 1e30 );
-
- // reconstruction cuts
- esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
- esdTrackCuts->SetRequireITSRefit(kTRUE);
- esdTrackCuts->SetRequireTPCRefit(kTRUE);
-
- esdTrackCuts->SetMinNClustersTPC(70);
- // if(cutDefinition<=kChi || cutDefinition>=kEtaGap01)
- // else esdTrackCuts->SetMinNClustersTPC(120);
- esdTrackCuts->SetMaxChi2PerClusterTPC(4);
- // if(cutDefinition!=kITSamy) esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst); // SPD first
- // esdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
-
- switch(cutDefinition) {
- case kTOFTRD2:
- esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst); // SPDfirst
- break;
- default:
- esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); // SPDany
- break;
- }
- */
- return esdTrackCuts;
-}
-
//______________________________________________________________________________________
void InitHistograms(AliDielectron *die, Int_t cutDefinition)
{
////// FLOW //////
- if(0) {
if(cutDefinition == kTOFTRD || cutDefinition >= kEtaGap01) {
histos->UserHistogram("Event","TPCxH2","TPC Qx component;TPCxH2",
100,-1500.,1500.,
AliDielectronVarManager::kCentrality,AliDielectronVarManager::kv0Av0CDiffH2);
}
- }
////// MONTE CARLO //////
/*
}
//// FLOW results use tprofiles
- if(0) {
if(cutDefinition == kTOFTRD || cutDefinition == kTOFTRD2 || cutDefinition >= kEtaGap01) {
histos->UserProfile("Pair","M_Cent_Pt_V0ACrpH2FlowV2",
125,0.,125*.04,
AliDielectronVarManager::kM);
}
- }
die->SetHistogramManager(histos);
}
TString etaMap="$TRAIN_ROOT/jpsi_JPSI/EtaCorrMaps.root";
TString trainRoot=gSystem->Getenv("TRAIN_ROOT");
- if (trainRoot.IsNull()) etaMap="../PWGDQ/dielectron/files/EtaCorrMaps.root";
+ if (trainRoot.IsNull()) etaMap="$ALICE_ROOT/PWGDQ/dielectron/files/EtaCorrMaps.root";
if (gSystem->AccessPathName(gSystem->ExpandPathName(etaMap.Data()))){
Error("ConfigPbPb","Eta map not found: %s",etaMap.Data());
return;
AliESDtrackCuts *SetupESDtrackCutsDieleData(Int_t cutDefinition);
-TString namesDieleData=("basicQ+SPDAny+pt>1+PID;basicQ+EMCal+SPDAny+pt6_20;basicQ+EMCal+SPDAny+pt8_20;EMCal+SPDAny+pt6_20;EMCal+SPDAny+pt8_20;EMCal+SPDFirst+pt6_20;EMCal+SPDFirst+pt8_20;EMCal+SPDAny");
+TString namesDieleData=("basicQ+SPDAny+pt>1+PID;EMCal+SPDAny+pt6_100;EMCal+SPDFirst+pt6_100;EMCal+SPDAny;EMCal+SPDFirst");
TObjArray *arrNamesDieleData=namesDieleData.Tokenize("; ");
const Int_t nDie=arrNamesDieleData->GetEntries();
InitHistogramsDieleData(diele, cutDefinition, isAOD);
// the last definition uses no cuts and only the QA histograms should be filled!
- if (cutDefinition==nDie-1) InitCFDieleData(diele, cutDefinition, isAOD);
+ if (cutDefinition>=nDie-2) InitCFDieleData(diele, cutDefinition, isAOD);
AliDielectronTrackRotator *rot=new AliDielectronTrackRotator;
rot->SetConeAnglePhi(TMath::Pi());
Centrality->AddCut(AliDielectronVarManager::kCentrality,0.,40.);
diele->GetEventFilter().AddCuts(vtxZ);
- if(cutDefinition<nDie-1) diele->GetEventFilter().AddCuts(Centrality);
-
+ if(cutDefinition<=2)
+ diele->GetEventFilter().AddCuts(Centrality);
}
diele->GetTrackFilter().AddCuts(SetupESDtrackCutsDieleData(cutDefinition));
} else {
AliDielectronTrackCuts *trackCuts=new AliDielectronTrackCuts("trackCuts","trackCuts");
- if (cutDefinition==5||cutDefinition==6)
+ if (cutDefinition==2||cutDefinition==4)
trackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
else
trackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
pt->AddCut(AliDielectronVarManager::kPt,1.,1e30);
pt->AddCut(AliDielectronVarManager::kKinkIndex0,0.);
//AOD additions since there are no AliESDtrackCuts -----------------
- pt->AddCut(AliDielectronVarManager::kNclsTPC,60.,160.);
pt->AddCut(AliDielectronVarManager::kEta,-0.9,0.9);
- if(cutDefinition>2&&cutDefinition!=nDie-1) pt->AddCut(AliDielectronVarManager::kTPCsignal,70.,120.);
- if(cutDefinition==nDie-1) pt->AddCut(AliDielectronVarManager::kTPCsignal,50.,150.);
+
+ if(cutDefinition==0) pt->AddCut(AliDielectronVarManager::kTPCsignal,55.,150.);
+ if(cutDefinition==1||cutDefinition==2) pt->AddCut(AliDielectronVarManager::kTPCsignal,70.,120.);
+ if(cutDefinition<=2) pt->AddCut(AliDielectronVarManager::kNclsTPC,80.,160.);
+ if(cutDefinition>2) {
+ pt->AddCut(AliDielectronVarManager::kTPCsignal,55.,150.);
+ pt->AddCut(AliDielectronVarManager::kNclsTPC,60.,160.);
+ }
+
diele->GetTrackFilter().AddCuts(pt);
// PID cuts --------------------------------------------------------
pid->AddCut(AliDielectronPID::kTPC,AliPID::kPion,-3.,3.,0.,0.,kTRUE);
pid->AddCut(AliDielectronPID::kTPC,AliPID::kProton,-20.,3.,0.,0.,kTRUE);
if(cutDefinition==0) diele->GetTrackFilter().AddCuts(pid);
- // PID cuts 2 ---- For EMCal purposes (all electron band) ----------------------------------------------------
- AliDielectronPID *pid2 = new AliDielectronPID("PIDforemcal","TPC nSigma |e|<3");
- pid2->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-3.,3.);
- if(cutDefinition==1||cutDefinition==2) diele->GetTrackFilter().AddCuts(pid2);
-
+
}
//______________________________________________________________________________________
AliDielectronVarCuts *pairCut=new AliDielectronVarCuts("0<M<5+|Y|<.9+PtEMCalleg","0<M<5 + |Y|<.9+PtEMCalleg");
pairCut->AddCut(AliDielectronVarManager::kM,0.,5.);
pairCut->AddCut(AliDielectronVarManager::kY,-0.9,0.9);
- if(cutDefinition==1&&cutDefinition==3&&cutDefinition==5) pairCut->AddCut(AliDielectronVarManager::kPt,6.,20.);
- if(cutDefinition==2&&cutDefinition==4&&cutDefinition==6) pairCut->AddCut(AliDielectronVarManager::kPt,8.,20.);
- if(cutDefinition==nDie-1) pairCut->AddCut(AliDielectronVarManager::kPt,2.,100.);
+ if(cutDefinition==1||cutDefinition==2) pairCut->AddCut(AliDielectronVarManager::kPt,6.,100.);
+ if(cutDefinition>2) pairCut->AddCut(AliDielectronVarManager::kPt,2.,100.);
diele->GetPairFilter().AddCuts(pairCut);
- if(cutDefinition>1){
- AliDielectronVarCuts *mycut = new AliDielectronVarCuts("ptCutEMCAL","cut for EMCal");
- mycut->AddCut(AliDielectronVarManager::kEMCALnSigmaEle,-3.,3.);
- mycut->AddCut(AliDielectronVarManager::kP,5.,1e30);
- AliDielectronPairLegCuts *varpair=new AliDielectronPairLegCuts();
- varpair->GetLeg1Filter().AddCuts(mycut);
- diele->GetPairFilter().AddCuts(varpair);
+ if(cutDefinition==1||cutDefinition==2){
+ AliDielectronVarCuts *mycut = new AliDielectronVarCuts("ptCutEMCAL","cut for EMCal");
+ mycut->AddCut(AliDielectronVarManager::kEMCALnSigmaEle,-3.,3.);
+ mycut->AddCut(AliDielectronVarManager::kP,6.,1e30);
+ AliDielectronPairLegCuts *varpair=new AliDielectronPairLegCuts();
+ varpair->GetLeg1Filter().AddCuts(mycut);
+ diele->GetPairFilter().AddCuts(varpair);
}
}
AliDielectronCF *cf=new AliDielectronCF(diele->GetName(),diele->GetTitle());
//pair variables
- cf->AddVariable(AliDielectronVarManager::kPt,"2.0, 4.0, 6., 8.0, 10.0, 100.0");
+ cf->AddVariable(AliDielectronVarManager::kPt,"2.0, 4.0, 6.0,7.0, 8.0, 10.0, 100.0");
cf->AddVariable(AliDielectronVarManager::kY,"-5,-1,-0.9,-0.8,-0.7,-0.5,-0.3,0.3,0.5,0.7,0.8,0.9,1.0,5");
cf->AddVariable(AliDielectronVarManager::kM,125,0.,125*.04); //40Mev Steps
cf->AddVariable(AliDielectronVarManager::kPairType,11,0,11);
- cf->AddVariable(AliDielectronVarManager::kOpeningAngle,"0.,0.2,0.5,1.,2.,3.14");
+ cf->AddVariable(AliDielectronVarManager::kOpeningAngle,"0.,0.2,0.5,0.8,1.,1.5,2.,3.14");
//leg variables
cf->AddVariable(AliDielectronVarManager::kPt,"0.0, 0.8, 1.0, 2.0, 3.0, 4.0, 5.0, 100.0",kTRUE);
+ cf->AddVariable(AliDielectronVarManager::kP,"0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 8.0, 10.0, 100.0",kTRUE);
cf->AddVariable(AliDielectronVarManager::kNclsTPC,"0, 70, 75, 80, 85, 90, 100, 120, 160",kTRUE);
cf->AddVariable(AliDielectronVarManager::kEta,"-5,-1,-0.9,-0.85,-0.8,-0.75,0.75,0.8,0.85,0.9,1.0,5",kTRUE);
- cf->AddVariable(AliDielectronVarManager::kTPCnSigmaEle,"-2.5,-2,-1.5,-1,-0.5,4.",kTRUE);
+ cf->AddVariable(AliDielectronVarManager::kEMCALnSigmaEle,"-3.,-2,-1.,1.,2.,3.,4.",kTRUE);
+ cf->AddVariable(AliDielectronVarManager::kTPCnSigmaEle,"-3.,-2.5,-2,-1.5,-1,-0.5,2.,3.,4.",kTRUE);
cf->AddVariable(AliDielectronVarManager::kTPCnSigmaPio,"3.,3.5,4.,100",kTRUE);
cf->AddVariable(AliDielectronVarManager::kTPCnSigmaPro,"3.,3.5,4.,100",kTRUE);
- cf->AddVariable(AliDielectronVarManager::kTPCsignal,"60.,65.,68.,70.,80.,90.,100.,120.,150.",kTRUE);
+ cf->AddVariable(AliDielectronVarManager::kTPCsignal,"60.,65.,68.,70.,72.,80.,90.,100.,120.,150.",kTRUE);
//event variables
cf->AddVariable(AliDielectronVarManager::kNaccTrcklts,"0.0, 9.0, 17.0, 25.0, 36.0, 55.0, 500.0");
- cf->AddVariable(AliDielectronVarManager::kCentrality,"0.0,20.0,40.0,80.0");
+ cf->AddVariable(AliDielectronVarManager::kCentrality,"0.0,10.0,20.0,40.0,60.0,80.0");
+ cf->AddVariable(AliDielectronVarManager::kNclsITS,"1,2,3,4,5,6",kTRUE);
if (!isAOD){
Bool_t hasMC=(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()!=0x0);
if (hasMC){