void InitHistograms(AliDielectron *die, Int_t cutDefinition);
+void AddQAHistsPID(AliDielectron *die);
+void AddQAHistsEP(AliDielectron *die);
+void AddQAHistsEff(AliDielectron *die);
+void AddHistsEleEff(AliDielectron *die);
+
void InitCF(AliDielectron* die, Int_t cutDefinition);
+void InitHF(AliDielectron* die, Int_t cutDefinition);
+void SetupEventCuts(AliDielectron *die, ULong64_t triggers, Int_t cutDefinition);
void SetupTrackCuts(AliDielectron *die, Int_t cutDefinition);
-void SetupPairCuts(AliDielectron *die, Int_t cutDefinition);
+void SetupV0Cuts( AliDielectron *die, Int_t cutDefinition);
+void SetupPairCuts( AliDielectron *die, Int_t cutDefinition);
+
+void ConfigEvtPlane(AliDielectron *die, Int_t cutDefinition);
+void ConfigBgrd( AliDielectron *die, Int_t cutDefinition);
-void AddMCSignals(AliDielectron *die);
-void SetEtaCorrection();
-TVectorD *GetRunNumbers();
+void AddMCSignals(AliDielectron *die, Int_t cutDefinition);
+void SetEtaCorrection(AliDielectron *die);
+TVectorD *GetRunNumbers2011();
+TVectorD *GetPDGcodes();
+TVectorD *GetDeltaPhiBins();
-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 };
+TString names=("Event;NoCut;PIDqa;LegEff;Std;SysMC;Flow;avgPt;Rec;TPC;TOF;TRD;TPC-TOF-TRD;NOPID;SysPt;SysEta;SysEle;SysPro;SysSPD;SysMCele");
+enum { kEvent, // event quantities (mult, ep, trigger, ...)
+ kNoCut, // pure event quantities (mult, ep, trigger, ...)
+ kPIDqa, // post calibration and validation of TPC PID
+ kLegEff, // single electron efficiency calculation
+ kStd, // standard Raa analysis
+ kSysMC,
+ kFlow, // flow calculation
+ kAvgPt, // mean pt and pt^2 analysis
+ kRec, // to calculate partial efficiencies, in particular the geom. acceptance
+ kTPC,
+ kTOF,
+ kTRD,
+ kTPCTOFTRD,
+ kNoPID,
+ kSysPt,
+ kSysEta,
+ kSysEle,
+ kSysPro,
+ kSysSPD,
+ kSysMCele,
+ kPIDQA
+ };
TObjArray *arrNames=names.Tokenize(";");
const Int_t nDie=arrNames->GetEntries();
-Bool_t hasMC = kFALSE;
+Bool_t isESD = kTRUE;
TString list = gSystem->Getenv("LIST");
-AliDielectron* ConfigJpsi_jb_PbPb(Int_t cutDefinition, TString prod="")
+AliDielectron* ConfigJpsi_jb_PbPb(Int_t cutDefinition, Bool_t hasMC=kFALSE, ULong64_t triggers=AliVEvent::kCentral | AliVEvent::kSemiCentral | AliVEvent::kMB)
{
//
// Setup the instance of AliDielectron
//
-
- // find mc or not?
- if( list.IsNull()) list=prod;
- if( list.Contains("LHC10h") || list.Contains("LHC11h") ) hasMC=kFALSE;
- if( list.Contains("LHC11a10") || list.Contains("LHC12a17") ) hasMC=kTRUE;
-
- // MC event handler?
- // hasMC=isMC;
- //(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()!=0x0);
+ // gsi train?
+ TString trainRoot = gSystem->Getenv("TRAIN_ROOT");
+ Bool_t isGSItrain = (trainRoot.IsNull()?kFALSE:kTRUE);
+ if(isGSItrain) {
+ // find mc or not?
+ if( list.IsNull()) list=prod;
+ if( list.Contains("LHC10h") || list.Contains("LHC11h") ) hasMC=kFALSE;
+ if( list.Contains("LHC11a10") || list.Contains("LHC12a17") || list.Contains("LHC13c7")) hasMC=kTRUE;
+ }
//ESD handler?
- Bool_t isESD=(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()->IsA()==AliESDInputHandler::Class());
-
- // switch off some configurations
- switch(cutDefinition) {
- case kTPC:
- case kTOF:
- case kTRD:
- return 0x0;
- break;
- //case kTOFTRD:
- case krec:
- if(!hasMC) return 0x0;
- break;
- //case kTOFTRD2:
- case kITScls:
- case kITSamy:
- case kDCA:
- case kChi:
- return 0x0;
- break;
- // case kGam0:
- // case kGam01:
- // case kGam05:
- // case kGam10:
- // case kGam15:
- // case kGam20:
- case kEtaGap01:
- case kEtaGap02:
- case kEtaGap03:
- case kEtaGap04:
- case kEtaGap05:
- case kSubLS:
- case kSubRndm:
- if( hasMC) return 0x0;
- break;
- }
-
- // create the actual framework object
+ isESD=(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()->IsA()==AliESDInputHandler::Class());
+
+ // task name
TString name=Form("%02d",cutDefinition);
- if (cutDefinition<arrNames->GetEntriesFast()){
- name=arrNames->At(cutDefinition)->GetName();
- }
- AliDielectron *die = new AliDielectron(Form("%s",name.Data()),
- Form("Track cuts: %s",name.Data()));
-
- // Monte Carlo Signals and TRD efficiency tables
- if(hasMC) {
- AddMCSignals(die);
-
- // trd tables
- TString pidTab="$TRAIN_ROOT/util/dielectron/dielectron/TRDpidEff_eleProb07_TRDntr4_6.root";
- TString trainRoot=gSystem->Getenv("TRAIN_ROOT");
- if (trainRoot.IsNull()) pidTab="$ALICE_ROOT/PWGDQ/dielectron/files/TRDpidEff_eleProb07_TRDntr4_6.root";
+ if (cutDefinition<arrNames->GetEntriesFast()) name=arrNames->At(cutDefinition)->GetName();
+ printf(" Adding %s%s config %s for %s \n",(isESD?"ESD":"AOD"),(hasMC?" MC":""),name.Data(),list.Data());
- if (gSystem->AccessPathName(gSystem->ExpandPathName(pidTab.Data())))
- Error("ConfigPbPb","PID table not found: %s",pidTab.Data());
- else
- die->SetTRDcorrectionFilename(pidTab.Data());
- }
-
- // cut setup
+ // init AliDielectron
+ AliDielectron *die = new AliDielectron(Form("%s",name.Data()), Form("%s",name.Data()));
+ die->SetHasMC(hasMC);
+
+ /* vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv CUTS vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv */
+ SetupEventCuts(die,triggers,cutDefinition);
SetupTrackCuts(die,cutDefinition);
+ SetupV0Cuts(die,cutDefinition);
SetupPairCuts(die,cutDefinition);
-
+
+ /* vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv MISC vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv */
+ // Monte Carlo Signals
+ AddMCSignals(die, cutDefinition);
+ // PID post calibartion
+ // if(cutDefinition!=kEvent && cutDefinition!=kLegEff) SetEtaCorrection(die);
+ // bgrd estimators
+ //ConfigBgrd(die,cutDefinition);
+ // tpc event plane configuration
+ //ConfigEvtPlane(die,cutDefinition);
+ // prefilter settings NEW
+ // if(cutDefinition==kNoCut ||cutDefinition==kEvent || cutDefinition==kNoPID || cutDefinition==kPIDqa || cutDefinition==kLegEff)
+ die->SetNoPairing();
+ // else
+ // die->SetPreFilterAllSigns();
+ //die->SetPreFilterUnlikeOnly();
+ // load single electron effieciency map ATTENTION
+ if(cutDefinition==kAvgPt) die->InitLegEffMap("alien:///alice/cern.ch/user/j/jbook/PWGDQ/dielectron/files/effMap.root");
+
+
+ /* vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv OUTPUT vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv */
// histogram setup
- if(cutDefinition == kTOFTRD ||
- cutDefinition == kGam0 ||
- cutDefinition == kTOFTRD2 ||
- cutDefinition >= kEtaGap01 )
- InitHistograms(die,cutDefinition);
-
+ InitHistograms(die,cutDefinition);
+ // histogram grid setup
+ InitHF(die,cutDefinition);
// CF container setup
- if(cutDefinition < kEtaGap01 ||
- cutDefinition == kSubRndm )
- InitCF(die,cutDefinition);
-
- // bgrd estimators
- if(!hasMC) {
-
- if(cutDefinition == kTOFTRD) {
- // rotations
- AliDielectronTrackRotator *rot=new AliDielectronTrackRotator;
- rot->SetIterations(10);
- rot->SetConeAnglePhi(TMath::Pi());
- rot->SetStartAnglePhi(TMath::Pi());
- die->SetTrackRotator(rot);
- // mixing
- AliDielectronMixingHandler *mix=new AliDielectronMixingHandler;
- mix->AddVariable(AliDielectronVarManager::kZvPrim,20,-10.,10.);
- mix->AddVariable(AliDielectronVarManager::kCentrality,"0,5,10,20,50,80");
- mix->SetMixType(AliDielectronMixingHandler::kAll);
- mix->SetDepth(50);
- die->SetMixingHandler(mix);
- }
-
-
- // TPC event plane configurations
- Double_t gGap;
- switch(cutDefinition) {
- case kEtaGap01: gGap=0.1; break;
- case kEtaGap02: gGap=0.2; break;
- case kEtaGap03: gGap=0.3; break;
- case kEtaGap04: gGap=0.4; break;
- case kEtaGap05: gGap=0.5; break;
- default: gGap=0.0;
- }
-
- AliDielectronVarCuts *poi = new AliDielectronVarCuts("PoI","PoI");
- poi->AddCut(AliDielectronVarManager::kM,2.92,3.20); // particles of interest, jpsi mass window
- die->GetEventPlanePOIPreFilter().AddCuts(poi);
-
- if(cutDefinition >= kEtaGap01 &&
- cutDefinition < kSubLS ) {
- AliDielectronVarCuts *etaGap = new AliDielectronVarCuts(AliDielectronVarManager::GetValueName(AliDielectronVarManager::kEta),"etaGap");
- etaGap->AddCut(AliDielectronVarManager::kEta,-1*gGap,gGap,kTRUE);
- die->GetEventPlanePreFilter().AddCuts(etaGap);
+ InitCF(die,cutDefinition);
+
+ // cut QA
+ die->SetCutQA();
+
+ return die;
+}
+
+//______________________________________________________________________________________
+void SetupEventCuts(AliDielectron *die, ULong64_t triggers, Int_t cutDefinition)
+{
+ //
+ // Setup the event cuts
+ //
+
+ // trigger specific centrality cuts (reject trigger inefficiencies)
+ Double_t minCent=0.0, maxCent=100.;
+ if(!die->GetHasMC()) {
+ switch(triggers) {
+ // case AliVEvent::kCentral: minCent= 0.; maxCent=10.; break; //0-9
+ // case AliVEvent::kSemiCentral: minCent=10.; maxCent=50.; break; //12-53
+ // case AliVEvent::kMB: minCent=10.; maxCent=90.; break;
+ default: minCent= 0.; maxCent=90.; break;
}
-
- if(cutDefinition==kSubLS) die->SetLikeSignSubEvents();
- die->SetPreFilterEventPlane();
- }
-
- // prefilter settings
- if(cutDefinition == kTOFTRD2)
- die->SetPreFilterAllSigns();
- else
- die->SetPreFilterUnlikeOnly();
-
- // 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;
+ // if(cutDefinition >= kEtaGap01) {minCent=20.; maxCent=50.;} // v2 analysis
+ if(cutDefinition == kSysMC) {minCent=0.; maxCent=50.;}
+ if(cutDefinition == kSysMCele) {minCent=0.; maxCent=50.;}
+
+ // VZERO multiplicity vs. number ob global tracks cut
+ TF1 *fMean = new TF1("fMean", "pol1", 0,25e+3);
+ fMean->SetParameters(691.633, 1.4892);
+ TF1 *fSigma = new TF1("fSigma","[0]+sqrt([1]*x+[2])",0,25e+3);
+ fSigma->SetParameters(-83.6599, 36.7677, 69530.7);
+
+ AliDielectronEventCuts *eventCuts=new AliDielectronEventCuts("vertex","vertex");
+ if(!isESD) eventCuts->SetVertexType(AliDielectronEventCuts::kVtxAny);
+ eventCuts->SetRequireVertex();
+ eventCuts->SetMinVtxContributors(1);
+ if(cutDefinition!=kNoCut) eventCuts->SetVertexZ(-10.,+10.);
+ eventCuts->SetCentralityRange(minCent,maxCent);
+ // eventCuts->SetCutOnV0MultipicityNTrks(fMean, fSigma, 4.0);
+ // eventCuts->SetRunRejection(AliDielectronHelper::MakeArbitraryBinning("170592,170593,170594"));
+ eventCuts->Print();
+ die->GetEventFilter().AddCuts(eventCuts);
+
}
//______________________________________________________________________________________
//
// Setup the track cuts
//
- //ESD handler?
- Bool_t isESD=(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()->IsA()==AliESDInputHandler::Class());
+ Bool_t hasMC=die->GetHasMC();
// Quality cuts
- AliDielectronCutGroup* cuts = new AliDielectronCutGroup("cuts","cuts",AliDielectronCutGroup::kCompAND);
- die->GetTrackFilter().AddCuts(cuts);
-
- //Pt cut, should make execution a bit faster
- 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);
-
- // 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);
- cuts->AddCut(varCuts);
-
- AliDielectronTrackCuts *trkCuts = new AliDielectronTrackCuts("TrkCuts","TrkCuts");
+ // die->GetTrackFilter().AddCuts(cuts);
+
+ /* vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv FILTER CUTS vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv */
+ // AOD track filter (needs to be first cut to speed up)
+ AliDielectronTrackCuts *trkFilter = new AliDielectronTrackCuts("filter","filter");
+ // config specific cuts
switch(cutDefinition) {
- case kTOFTRD2: trkCuts->SetClusterRequirementITS(AliDielectronTrackCuts::kSPD,AliDielectronTrackCuts::kFirst);
+ case kPIDqa: trkFilter->SetAODFilterBit(AliDielectronTrackCuts::kTPCqual); break;
+ default: trkFilter->SetAODFilterBit(AliDielectronTrackCuts::kTPCqualSPDany);
+ }
+ // trkFilter->SetMinNCrossedRowsOverFindable(0.6);
+
+ /* vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv TRACK CUTS vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv */
+ AliDielectronVarCuts *varAccCuts = new AliDielectronVarCuts("acc","acc");
+ AliDielectronCutGroup *grpRecCuts = new AliDielectronCutGroup("rec","rec",AliDielectronCutGroup::kCompAND);
+ AliDielectronVarCuts *varRecCuts = new AliDielectronVarCuts("VarRecCuts","VarRecCuts");
+ AliDielectronTrackCuts *trkRecCuts = new AliDielectronTrackCuts("TrkRecCuts","TrkRecCuts");
+
+ // config specific cuts
+ switch(cutDefinition) {
+ case kPIDqa:
+ varAccCuts->AddCut(AliDielectronVarManager::kPt, 0.5, 1e30); // 0.85 ATTENTION
+ varRecCuts->AddCut(AliDielectronVarManager::kTPCclsSegments,7., 8.0);
+ break;
+ case kEvent:
+ case kNoCut:
+ case kFlow:
+ case kRec:
+ case kStd:
+ case kTRD:
+ case kTPCTOFTRD:
+ case kAvgPt:
+ case kNoPID:
+ case kSysMCele:
+ varAccCuts->AddCut(AliDielectronVarManager::kPt, 0.85, 1e30);
+ varAccCuts->AddCut(AliDielectronVarManager::kEta, -0.8, 0.8);
+ case kLegEff:
+ varRecCuts->AddCut(AliDielectronVarManager::kTPCclsSegments,7., 8.0);
+ trkRecCuts->SetITSclusterCut(AliDielectronTrackCuts::kOneOf, 3); // SPD any
+ //if(hasMC) varAccCuts->AddCut(AliDielectronVarManager::kEta, -0.9, 0.9);
+ //else
+ ///
+ break;
+ case kTOF:
+ case kTPC:
+ varAccCuts->AddCut(AliDielectronVarManager::kPt, 0.95, 1e30); //0.8
+ varAccCuts->AddCut(AliDielectronVarManager::kEta, -0.9, 0.9);
+ varRecCuts->AddCut(AliDielectronVarManager::kNclsTPC, 70.0, 160.0);
+ trkRecCuts->SetITSclusterCut(AliDielectronTrackCuts::kOneOf, 7); // ITS-3 = 1+2+4
+ break;
+ ///////////////////////////////////////////////////////////////////////////////////////////// systematics
+ case kSysPt:
+ varAccCuts->AddCut(AliDielectronVarManager::kPt, 1.1, 20./*1e30*/); ///ATTENTION
+ varAccCuts->AddCut(AliDielectronVarManager::kEta, -0.8, 0.8);
+ varRecCuts->AddCut(AliDielectronVarManager::kTPCclsSegments,7., 8.0);
+ trkRecCuts->SetITSclusterCut(AliDielectronTrackCuts::kOneOf, 3); // SPD any
break;
- default: trkCuts->SetClusterRequirementITS(AliDielectronTrackCuts::kSPD,AliDielectronTrackCuts::kAny);
+ case kSysEta:
+ varAccCuts->AddCut(AliDielectronVarManager::kPt, 0.85, 20./*1e30*/);
+ varAccCuts->AddCut(AliDielectronVarManager::kEta, -0.9, 0.9); ////ATTENTION
+ varRecCuts->AddCut(AliDielectronVarManager::kTPCclsSegments,7., 8.0);
+ trkRecCuts->SetITSclusterCut(AliDielectronTrackCuts::kOneOf, 3); // SPD any
+ break;
+ case kSysEle:
+ case kSysPro:
+ varAccCuts->AddCut(AliDielectronVarManager::kPt, 0.85, 20./*1e30*/);
+ varAccCuts->AddCut(AliDielectronVarManager::kEta, -0.8, 0.8);
+ varRecCuts->AddCut(AliDielectronVarManager::kTPCclsSegments,7., 8.0);
+ trkRecCuts->SetITSclusterCut(AliDielectronTrackCuts::kOneOf, 3); // SPD any
+ break;
+ case kSysSPD:
+ varAccCuts->AddCut(AliDielectronVarManager::kPt, 0.85, 20./*1e30*/);
+ varAccCuts->AddCut(AliDielectronVarManager::kEta, -0.8, 0.8);
+ varRecCuts->AddCut(AliDielectronVarManager::kTPCclsSegments,7., 8.0);
+ trkRecCuts->SetITSclusterCut(AliDielectronTrackCuts::kOneOf, 1); // SPD any
+ break;
+ case kSysMC:
+ varAccCuts->AddCut(AliDielectronVarManager::kPt, 0.85, 1e30);
+ varAccCuts->AddCut(AliDielectronVarManager::kEta, -0.9, 0.9); ////ATTENTION
+ varRecCuts->AddCut(AliDielectronVarManager::kTPCclsSegments,7., 8.0);
+ trkRecCuts->SetITSclusterCut(AliDielectronTrackCuts::kOneOf, 3); // SPD any
break;
}
- trkCuts->SetRequireITSRefit(kTRUE);
- trkCuts->SetRequireTPCRefit(kTRUE);
- cuts->AddCut(trkCuts);
-
- //Do we have an MC handler?
- // Bool_t hasMC=(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()!=0x0);
-
+
+ // standrad reconstruction cuts
+ varRecCuts->AddCut(AliDielectronVarManager::kImpactParXY, -1.0, 1.0);
+ varRecCuts->AddCut(AliDielectronVarManager::kImpactParZ, -3.0, 3.0);
+ varRecCuts->AddCut(AliDielectronVarManager::kTPCchi2Cl, 0.0, 4.0);
+ varRecCuts->AddCut(AliDielectronVarManager::kKinkIndex0, 0.0);
+ // varRecCuts->AddCut(AliDielectronVarManager::kV0Index0, 0.0);
+ if(cutDefinition!=kPIDqa) trkRecCuts->SetRequireITSRefit(kTRUE);
+ trkRecCuts->SetRequireTPCRefit(kTRUE);
+
/* vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv PID CUTS vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv */
- AliDielectronPID *pid = new AliDielectronPID("PID","PID");
-
- ////////////////////////////////// DATA
- if(!hasMC) {
- pid->AddCut(AliDielectronPID::kTPC,AliPID::kPion,-100.,3.5,0.,0.,kTRUE);
- pid->AddCut(AliDielectronPID::kTPC,AliPID::kProton,-100.,3.5,0.,0.,kTRUE);
-
- if(cutDefinition==kTRD || cutDefinition>=kTOFTRD || cutDefinition>=kTOFTRD2)
- pid->AddCut(AliDielectronPID::kTRDeleEff,AliPID::kElectron,.8,1.,3.5.,6.,kFALSE,
- AliDielectronPID::kIfAvailable,AliDielectronVarManager::kTRDpidQuality);
+ AliDielectronCutGroup *grpPIDCuts = new AliDielectronCutGroup("PID","PID",AliDielectronCutGroup::kCompAND);
+ AliDielectronVarCuts *pidSelCuts = new AliDielectronVarCuts("selPIDCuts","selPIDCuts");
+ pidSelCuts->AddCut(AliDielectronVarManager::kTRDpidQuality, 4.0, 6.0);
+ pidSelCuts->AddCut(AliDielectronVarManager::kTRDchi2, 0.0, 2.0);
+ AliDielectronVarCuts *pidVarCuts = new AliDielectronVarCuts("varPIDCuts","varPIDCuts");
+ AliDielectronVarCuts *pidMCCuts = new AliDielectronVarCuts("mcPIDCuts","mcPIDCuts");
+ AliDielectronPID *pidCuts = new AliDielectronPID("PIDCuts","PIDCuts");
+
+ switch(cutDefinition) {
+ case kPIDqa:
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-10.,10.);
+ pidCuts->AddCut(AliDielectronPID::kTRDeleEff,AliPID::kElectron,.95,1.,pidSelCuts,
+ kFALSE, AliDielectronPID::kIfAvailable);
+ break;
+ case kTOF:
+ pidVarCuts->AddCut(AliDielectronVarManager::kTOFbeta, 0.2, 0.9, kTRUE);
+ pidCuts->AddCut(AliDielectronPID::kTOF,AliPID::kElectron,-3,3.,0.,0.,kFALSE,
+ AliDielectronPID::kIfAvailable);
+ case kTPC:
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-3.,4.);
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kPion,-100.,4.0,0.,0.,kTRUE);
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kProton,-100.,3.5,0.,0.,kTRUE);
+ break;
+ case kTRD:
+ pidCuts->AddCut(AliDielectronPID::kTRDeleEff,AliPID::kElectron,.95,1.,pidSelCuts,
+ kFALSE, AliDielectronPID::kIfAvailable);
+ case kEvent:
+ case kNoCut:
+ case kFlow:
+ case kStd:
+ case kLegEff:
+ case kAvgPt:
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-1.5.,3.);
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-0.9,3.0,-0.5,+0.5,kFALSE,
+ AliDielectronPID::kRequire,AliDielectronVarManager::kEta);
+ // pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-0.65,3.0,-0.3,+0.3,kFALSE,
+ // AliDielectronPID::kRequire,AliDielectronVarManager::kEta);
+ // pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-0.4,3.0,-0.1,+0.1,kFALSE,
+ // AliDielectronPID::kRequire,AliDielectronVarManager::kEta);
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kProton,-100.,4.0,0.,0.,kTRUE);
+ // tof heavy particle exclusion
+ // pidVarCuts->AddCut(AliDielectronVarManager::kTOFbeta, 0.3, 0.7, kTRUE);
+ //pidCuts->AddCut(AliDielectronPID::kTOF,AliPID::kElectron,-3,3.,0.,0.,kFALSE,
+ // AliDielectronPID::kIfAvailable);
+ break;
+ case kTPCTOFTRD:
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-1.3.,3.);
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-0.9,3.0,-0.5,+0.5,kFALSE,
+ AliDielectronPID::kRequire,AliDielectronVarManager::kEta);
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-0.65,3.0,-0.3,+0.3,kFALSE,
+ AliDielectronPID::kRequire,AliDielectronVarManager::kEta);
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-0.4,3.0,-0.1,+0.1,kFALSE,
+ AliDielectronPID::kRequire,AliDielectronVarManager::kEta);
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kProton,-100.,4.0,0.,0.,kTRUE);
+ pidVarCuts->AddCut(AliDielectronVarManager::kTOFbeta, 0.2, 0.9, kTRUE);
+ pidCuts->AddCut(AliDielectronPID::kTOF,AliPID::kElectron,-3,3.,0.,0.,kFALSE,
+ AliDielectronPID::kIfAvailable);
+ pidCuts->AddCut(AliDielectronPID::kTRDeleEff,AliPID::kElectron,.95,1.,pidSelCuts,
+ kFALSE, AliDielectronPID::kIfAvailable);
+ break;
+ ///////////////////////////////////////////////////////////////////////////////////////////// systematics
+ case kSysPt:
+ case kSysEta:
+ case kSysSPD:
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-1.5.,3.);
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-0.9,3.0,-0.5,+0.5,kFALSE,
+ AliDielectronPID::kRequire,AliDielectronVarManager::kEta);
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kProton,-100.,4.0,0.,0.,kTRUE);
+ break;
+ case kSysEle:
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,+0.7.,3.);
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kProton,-100.,4.0,0.,0.,kTRUE);
+ break;
+ case kSysPro:
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-1.5.,3.);
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-0.9,3.0,-0.5,+0.5,kFALSE,
+ AliDielectronPID::kRequire,AliDielectronVarManager::kEta);
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kProton,-100.,3.5,0.,0.,kTRUE);
+ break;
+ case kSysMC:
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-1.5.,3.);
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-0.9,3.0,-0.5,+0.5,kFALSE,
+ AliDielectronPID::kRequire,AliDielectronVarManager::kEta);
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kProton,-100.,3.5,0.,0.,kTRUE);
+ break;
+ case kSysMCele:
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-1.5.,3.);
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,-0.9,3.0,-0.5,+0.5,kFALSE,
+ AliDielectronPID::kRequire,AliDielectronVarManager::kEta);
+ pidCuts->AddCut(AliDielectronPID::kTPC,AliPID::kProton,-100.,4.0,0.,0.,kTRUE);
+ break;
}
-
- ////////////////////////////////// MC
- if(hasMC) {
-
- // electron
- Double_t nSigmaPi = 3.5; Double_t nSigmaP = 3.5;
- Double_t resolution=0.0549;
- Double_t BBpro[5] = {0};
- Double_t BBpio[5] = {0};
-
- for(Int_t icent=0; icent<8; icent++) {
-
- switch (icent) {
- case 0: // 0-10%
- BBpro[0] = 0.031555; BBpro[1] = 26.0595; BBpro[2] = 3.02422e-11; BBpro[3] = 2.05594; BBpro[4] = 5.99848;
- BBpio[0] = 0.0252122; BBpio[1] = 38.8991; BBpio[2] = 4.0901e-11; BBpio[3] = 5.27988; BBpio[4] = 4.3108;
- break;
- case 1: // 10-20%
- BBpro[0] = 0.0315171; BBpro[1] = 25.8656; BBpro[2] = 3.03896e-11; BBpro[3] = 2.05802; BBpro[4] = 5.99999;
- BBpio[0] = 0.0252127; BBpio[1] = 33.8617; BBpio[2] = 3.56866e-11; BBpio[3] = 5.24831; BBpio[4] = 4.31093;
- break;
- case 2: // 20-30%
- BBpro[0] = 0.0315171; BBpro[1] = 25.8656; BBpro[2] = 3.03896e-11; BBpro[3] = 2.05802; BBpro[4] = 5.99999;
- BBpio[0] = 0.0263205; BBpio[1] = 37.9307; BBpio[2] = 4.29724e-11; BBpio[3] = 5.74458; BBpio[4] = 4.32459;
- break;
- case 3: // 30-40%
- BBpro[0] = 0.0315171; BBpro[1] = 25.8656; BBpro[2] = 3.03896e-11; BBpro[3] = 2.05802; BBpro[4] = 5.99999;
- BBpio[0] = 0.026294; BBpio[1] = 39.0346; BBpio[2] = 4.12261e-11; BBpio[3] = 5.28808; BBpio[4] = 4.31301;
- break;
- case 4: // 40-50%
- BBpro[0] = 0.0315171; BBpro[1] = 25.8656; BBpro[2] = 3.03896e-11; BBpro[3] = 2.05802; BBpro[4] = 5.99999;
- BBpio[0] = 0.0263134; BBpio[1] = 38.2084; BBpio[2] = 3.75159e-11; BBpio[3] = 5.78125; BBpio[4] = 4.31363;
- break;
- case 5: // 50-60%
- BBpro[0] = 0.0315171; BBpro[1] = 25.8656; BBpro[2] = 3.03896e-11; BBpro[3] = 2.05802; BBpro[4] = 5.99999;
- BBpio[0] = 0.0263134; BBpio[1] = 38.2084; BBpio[2] = 3.75159e-11; BBpio[3] = 5.78125; BBpio[4] = 4.31363;
- break;
- case 6: // 60-70%
- BBpro[0] = 0.031555; BBpro[1] = 26.0595; BBpro[2] = 3.02422e-11; BBpro[3] = 2.05594; BBpro[4] = 5.99848;
- BBpio[0] = 0.026302; BBpio[1] = 38.6888; BBpio[2] = 3.56792e-11; BBpio[3] = 5.2465; BBpio[4] = 4.31094;
- break;
- case 7: // 70-80%
- BBpro[0] = 0.0315171; BBpro[1] = 25.8656; BBpro[2] = 3.03896e-11; BBpro[3] = 2.05802; BBpro[4] = 5.99999;
- BBpio[0] = 0.0263134; BBpio[1] = 38.2084; BBpio[2] = 3.75159e-11; BBpio[3] = 5.78125; BBpio[4] = 4.31363;
- break;
- case 8: // 80-90%
- BBpro[0] = 0.0313438; BBpro[1] = 25.8666; BBpro[2] = 4.5457e-11; BBpro[3] = 2.07912; BBpro[4] = 5.99986;
- BBpio[0] = 0.0252127; BBpio[1] = 33.8617; BBpio[2] = 3.56866e-11; BBpio[3] = 5.24831; BBpio[4] = 4.31093;
- break;
- case 9: // 90-100%
- BBpro[0] = 0.0319126; BBpro[1] = 36.8784; BBpro[2] = 3.4274e-11; BBpro[3] = 3.2431; BBpro[4] = 5.93388;
- BBpio[0] = 0.027079; BBpio[1] = 67.5936; BBpio[2] = 9.72548e-11; BBpio[3] = 9.61382; BBpio[4] = 5.99372;
- break;
- }
-
-
- TF1 *ffPro=new TF1(Form("fBethe%d_c%d",AliPID::kProton,icent), Form("(%f*%f+(AliExternalTrackParam::BetheBlochAleph(x/%f,[0],[1],[2],[3],[4])-AliExternalTrackParam::BetheBlochAleph(x/%f,[0],[1],[2],[3],[4])))/%f", nSigmaP,resolution, AliPID::ParticleMass(AliPID::kProton), AliPID::ParticleMass(AliPID::kElectron), resolution), 0.05,200.);
-
- TF1 *ffPio=new TF1(Form("fBethe%d_c%d",AliPID::kPion,icent), Form("(%f*%f+(AliExternalTrackParam::BetheBlochAleph(x/%f,[0],[1],[2],[3],[4])-AliExternalTrackParam::BetheBlochAleph(x/%f,[0],[1],[2],[3],[4])))/%f", nSigmaPi,resolution, AliPID::ParticleMass(AliPID::kPion), AliPID::ParticleMass(AliPID::kElectron), resolution), 0.05,200.);
-
- //LHC11a10b
- if (list.Contains("LHC11a10b")) {
- 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]);
-
- // proton cut
- pid->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,ffPro,10,((double)icent)*10.,((double)icent+1)*10,
- kFALSE,AliDielectronPID::kRequire,AliDielectronVarManager::kCentrality);
- // pion cut
- pid->AddCut(AliDielectronPID::kTPC,AliPID::kElectron,ffPio,10,((double)icent)*10.,((double)icent+1)*10,
- kFALSE,AliDielectronPID::kRequire,AliDielectronVarManager::kCentrality);
- }
- }
-
- // shifts for the nSigma electrons
- TGraph* nSigmaCorrection = new TGraph();
- // LHC11a10b
- if (list.Contains("LHC11a10b")) {
- nSigmaCorrection->SetPoint(0, 137161., -0.50-(0.28));
- nSigmaCorrection->SetPoint(1, 139510., -0.50-(0.28));
- pid->SetCorrGraph(nSigmaCorrection);
- }
-
- } //hasMC
-
- ////////////////////////////////// DATA + MC
- // 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);
-
- if(cutDefinition!=krec) cuts->AddCut(pid);
- /* ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ PID CUTS ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^*/
-
-
+
+ // mc identification
+ if(cutDefinition==kPIDqa && hasMC) {
+ pidMCCuts->SetCutType(AliDielectronVarCuts::kAny);//only apply any of the two cuts
+ pidMCCuts->AddCut(AliDielectronVarManager::kPdgCode,-11.5,-10.5 );
+ pidMCCuts->AddCut(AliDielectronVarManager::kPdgCode,10.5,11.5 );
+ }
+
+ /* vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv PID POST CORRECTION vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv */
+ // if(cutDefinition!=kEvent) SetEtaCorrection(pidCuts,hasMC);
+
+
+ /* vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv TENDER CUTS vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv */
// exclude conversion electrons selected by the tender
AliDielectronTrackCuts *noconv=new AliDielectronTrackCuts("noConv","noConv");
noconv->SetV0DaughterCut(AliPID::kElectron,kTRUE);
- cuts->AddCut(noconv);
-
+
+ // activate the cut sets (order might be CPU timewise important)
+ switch(cutDefinition) {
+ case kNoPID:
+ case kPIDqa:
+ if(!isESD) die->GetTrackFilter().AddCuts(trkFilter);
+ die-> GetTrackFilter().AddCuts(varAccCuts);
+ grpRecCuts->AddCut(trkRecCuts);
+ grpRecCuts->AddCut(varRecCuts);
+ die-> GetTrackFilter().AddCuts(grpRecCuts);
+ if(hasMC) grpPIDCuts->AddCut(pidMCCuts);
+ else grpPIDCuts->AddCut(pidCuts);
+ die-> GetTrackFilter().AddCuts(grpPIDCuts);
+ break;
+ case kRec:
+ die->GetTrackFilter().AddCuts(varAccCuts);
+ grpRecCuts->AddCut(trkRecCuts);
+ grpRecCuts->AddCut(varRecCuts);
+ die->GetTrackFilter().AddCuts(grpRecCuts);
+ break;
+ case kEvent:
+ case kNoCut:
+ case kFlow:
+ case kStd:
+ case kLegEff:
+ case kTRD:
+ case kTPCTOFTRD:
+ case kAvgPt:
+ case kTOF:
+ case kTPC:
+ case kSysPt:
+ case kSysEta:
+ case kSysEle:
+ case kSysPro:
+ case kSysSPD:
+ case kSysMC:
+ case kSysMCele:
+ if(!isESD) die->GetTrackFilter().AddCuts(trkFilter);
+ die-> GetTrackFilter().AddCuts(varAccCuts);
+ grpRecCuts->AddCut(trkRecCuts);
+ grpRecCuts->AddCut(varRecCuts);
+ die-> GetTrackFilter().AddCuts(grpRecCuts);
+ grpPIDCuts->AddCut(pidCuts);
+ // grpPIDCuts->AddCut(pidVarCuts);
+ die-> GetTrackFilter().AddCuts(grpPIDCuts);
+ //cuts->AddCut(noconv);
+ // debug
+ trkFilter->Print();
+ varAccCuts->Print();
+ grpRecCuts->Print();
+ grpPIDCuts->Print();
+ }
+
+}
+
+//______________________________________________________________________________________
+void SetupV0Cuts(AliDielectron *die, Int_t cutDefinition)
+{
+ //
+ // Setup the V0 cuts
+ //
+
+ switch(cutDefinition) {
+ case kEvent: return;
+ case kNoCut: return;
+ // case kLegEff: return;
+ }
+
+ Bool_t bRej = kTRUE;
+ Int_t defPID = 16;
+ if(cutDefinition==kPIDqa) {
+ bRej = kFALSE;
+ defPID = 13;//13
+ }
+
+ AliDielectronV0Cuts *gammaV0Cuts = new AliDielectronV0Cuts("V0","V0");
+ gammaV0Cuts->SetV0finder(AliDielectronV0Cuts::kOnTheFly);
+ gammaV0Cuts->SetPdgCodes(22,11,11);
+ gammaV0Cuts->SetDefaultPID(defPID);
+ gammaV0Cuts->AddCut(AliDielectronVarManager::kCosPointingAngle, TMath::Cos(0.05), 1.0, kFALSE); //0.02 -- 0.05
+ gammaV0Cuts->AddCut(AliDielectronVarManager::kArmPt, 0.0, 0.05, kFALSE);
+ gammaV0Cuts->AddCut(AliDielectronVarManager::kArmAlpha, -0.35, 0.35, kFALSE);
+ gammaV0Cuts->AddCut(AliDielectronVarManager::kM, 0.0, 0.1, kFALSE); //0.05 -- 0.1
+ gammaV0Cuts->AddCut(AliDielectronVarManager::kChi2NDF, 0.0, 10.0, kFALSE);
+ gammaV0Cuts->AddCut(AliDielectronVarManager::kLegDist, 0.0, 0.25, kFALSE);
+ gammaV0Cuts->AddCut(AliDielectronVarManager::kR, 3.0, 90.0, kFALSE);
+ gammaV0Cuts->AddCut(AliDielectronVarManager::kPsiPair, 0.0, 0.2, kFALSE); //0.05 -- 0.2
+ if(cutDefinition==kPIDqa && 0) {
+ gammaV0Cuts->AddCut(AliDielectronVarManager::kImpactParXY, 0.0, 1.0, kFALSE);
+ gammaV0Cuts->AddCut(AliDielectronVarManager::kImpactParZ, 0.0, 1.0, kFALSE);
+ }
+ // gammaV0Cuts->AddCut(AliDielectronVarManager::kOpeningAngle, 0.0, 0.1, kFALSE);
+ gammaV0Cuts->SetExcludeTracks(bRej);
+ gammaV0Cuts->Print();
+
+ // const Double_t |cutAlphaG| < 0.35; && const Double_t cutQTG < 0.05;
+ // const Double_t |cutAlphaG2|[2] = {0.6, 0.8}; && const Double_t cutQTG2 < 0.04;
+
+ // if(cuts)
+ // ((AliDielectronCutGroup*)cuts)->AddCut(gammaV0Cuts);
+ // else
+ die->GetTrackFilter().AddCuts(gammaV0Cuts);
}
//______________________________________________________________________________________
//
// Setup the pair cuts
//
-
- // conversion rejection
- Double_t gCut;
+ Bool_t hasMC=die->GetHasMC();
+
+ // rap and mass rejection
+ Double_t gCut=0.05, yCut=0.9, eCut=0.01;
switch(cutDefinition) {
- case kTPC: gCut=0.05; break;
- case krec: gCut=0.05; break;
- case kGam10: gCut=0.1; break;
- case kGam15: gCut=0.15; break;
- case kGam20: gCut=0.2; break;
- case kGam05: gCut=0.05; break;
- case kGam01: gCut=0.01; break;
- case kGam0: gCut=0.0; break;
- default: gCut=0.05; // default
+ case kNoCut:
+ case kEvent: yCut=0.0; break;
+ case kPIDqa: return;
+ case kFlow: yCut=0.8; break;
+ case kRec: yCut=0.9; break;
+ case kNoPID: return;
+ case kTPC: yCut=0.9; break;
+ case kTOF: yCut=0.9; break;
+ case kTRD: yCut=0.8; break;
+ case kStd: yCut=0.8; break;
+ case kLegEff: return;
+ case kTPCTOFTRD: yCut=0.8; break;
+ case kAvgPt: yCut=0.8; break;
+ case kSysMC: yCut=0.9; break;
+ case kSysMCele: yCut=0.8; break;
+ case kSysEta: yCut=0.9; break;
+ case kSysPt:
+ case kSysEle:
+ case kSysPro:
+ case kSysSPD:
+ yCut=0.8; break;
+ // default: gCut=0.05; // default
}
-
+
+ // MC
+ //if(hasMC) yCut=0.9;
+
+ // rapidity selection
+ AliDielectronVarCuts *rapCut=new AliDielectronVarCuts(Form("|Y|<%.1f",yCut),Form("|Y|<%.1f",yCut));
+ rapCut->AddCut(AliDielectronVarManager::kY,-1.*yCut,yCut);
+ die->GetPairFilter().AddCuts(rapCut);
+
+ // gamma rejection
AliDielectronVarCuts *gammaCuts = new AliDielectronVarCuts("GammaCuts","GammaCuts");
-// gammaCuts->AddCut(AliDielectronVarManager::kOpeningAngle, 0.0, 0.1, kTRUE);
-// gammaCuts->AddCut(AliDielectronVarManager::kLegDist, 0.0, 0.25, kTRUE);
-// gammaCuts->AddCut(AliDielectronVarManager::kR, 3.0, 90.0, kTRUE);
-// gammaCuts->AddCut(AliDielectronVarManager::kPsiPair, 0.0, 0.05, kTRUE);
-// gammaCuts->AddCut(AliDielectronVarManager::kChi2NDF, 0.0, 10.0, kTRUE);
gammaCuts->AddCut(AliDielectronVarManager::kM, 0.0, gCut);
die->GetPairPreFilter().AddCuts(gammaCuts);
-
-
- // rapidity selection
- AliDielectronVarCuts *rapCut=new AliDielectronVarCuts("|Y|<.9","|Y|<.9");
- rapCut->AddCut(AliDielectronVarManager::kY,-0.9,0.9);
- die->GetPairFilter().AddCuts(rapCut);
-
+
+ // pair efficiency cut
+ // if(cutDefinition==kAvgPt && !hasMC) {
+ // AliDielectronVarCuts *effCut=new AliDielectronVarCuts(Form("(Axe)>%.2f",eCut),Form("(Axe)>%.2f",eCut));
+ // effCut->AddCut(AliDielectronVarManager::kPairEff,0.0,eCut,kTRUE);
+ // die->GetPairFilter().AddCuts(effCut);
+ // }
+
+ // random signal rejection
+ AliDielectronCutGroup *grpRNDMCuts = new AliDielectronCutGroup("RNDM","RNDM",AliDielectronCutGroup::kCompOR);
+ AliDielectronVarCuts *exclCutCEN=new AliDielectronVarCuts("exclCEN","exclCEN");
+ exclCutCEN->SetCutType(AliDielectronVarCuts::kAll); // all criteria need to be fullfilled
+ exclCutCEN->AddCut(AliDielectronVarManager::kRndmPair,0.0,1./40);
+ exclCutCEN->AddCut(AliDielectronVarManager::kPdgCode,443);
+ exclCutCEN->AddCut(AliDielectronVarManager::kCentrality,0.,10.);
+ AliDielectronVarCuts *exclCutSEMI=new AliDielectronVarCuts("exclSEMI","exclSEMI");
+ exclCutSEMI->SetCutType(AliDielectronVarCuts::kAll); // all criteria need to be fullfilled
+ exclCutSEMI->AddCut(AliDielectronVarManager::kRndmPair,0.0,1./20);
+ exclCutSEMI->AddCut(AliDielectronVarManager::kPdgCode,443);
+ exclCutSEMI->AddCut(AliDielectronVarManager::kCentrality,10.,90.);
+ AliDielectronVarCuts *inclCut=new AliDielectronVarCuts("incl","incl");
+ inclCut->AddCut(AliDielectronVarManager::kPdgCode,443,443,kTRUE);
+ grpRNDMCuts->AddCut(exclCutCEN);
+ grpRNDMCuts->AddCut(exclCutSEMI);
+ grpRNDMCuts->AddCut(inclCut);
+ if(hasMC && cutDefinition==kStd && 0) die->GetPairFilter().AddCuts(grpRNDMCuts); //ATTENTION
+
+}
+
+//______________________________________________________________________________________
+void ConfigBgrd(AliDielectron *die, Int_t cutDefinition)
+{
+ //
+ // Configurate the background estimators
+ //
+
+ // default no processing of LS
+ // die->SetProcessLS(kFALSE);
+
+ // skip config
+ switch(cutDefinition) {
+ case kEvent:
+ case kNoCut:
+ case kNoPID:
+ case kPIDqa:
+ case kRec:
+ case kLegEff:
+ return;
+ }
+
+ Bool_t hasMC=die->GetHasMC();
+ if(hasMC && cutDefinition!=kAvgPt) return;
+
+ // add track rotations
+ AliDielectronTrackRotator *rot=new AliDielectronTrackRotator;
+ rot->SetIterations(10);
+ rot->SetConeAnglePhi(TMath::Pi());
+ rot->SetStartAnglePhi(TMath::Pi());
+ // die->SetTrackRotator(rot);
+
+ // add mixed events
+ AliDielectronMixingHandler *mix=new AliDielectronMixingHandler;
+ switch(cutDefinition) {
+ // case kStd:
+ // mix->AddVariable(AliDielectronVarManager::kZvPrim, 20, -10., 10.);
+ // mix->AddVariable(AliDielectronVarManager::kCentrality, 36, 0., 90.);
+ // mix->AddVariable(AliDielectronVarManager::kv0ACrpH2, 10, 0., TMath::Pi());
+ // break;
+ default:
+ mix->AddVariable(AliDielectronVarManager::kZvPrim, "-10.,-5.,-4.,-3.,-2.,-1.,0.,1.,2.,3.,4.,5.,10.");
+ mix->AddVariable(AliDielectronVarManager::kCentrality, 9, 0.,90.);
+ mix->AddVariable(AliDielectronVarManager::kTPCrpH2, 10, TMath::Pi()/-2, TMath::Pi()/2); // max res: 10%->10bins // 8bins
+ // mix->AddVariable(AliDielectronVarManager::kTPCmagH2, "0.,20.,50.,80.,110.,150.,500.");
+ break;
+ }
+ mix->SetSkipFirstEvent(kTRUE); // needed for flow analysis
+ mix->SetMixType(AliDielectronMixingHandler::kAll);
+ mix->SetDepth(150);
+ die->SetMixingHandler(mix);
+
+}
+
+//______________________________________________________________________________________
+void ConfigEvtPlane(AliDielectron *die, Int_t cutDefinition)
+{
+ //
+ // Configurate the TPC event plane
+ //
+
+ switch(cutDefinition) {
+ // case kEvent:
+ case kNoPID:
+ case kPIDqa:
+ // case kRec: // TOTHINK: might be needed when checking efficiencies versus dPhi
+ case kLegEff:
+ return;
+ }
+
+ Bool_t hasMC=die->GetHasMC();
+ if(hasMC && cutDefinition!=kAvgPt) return;
+
+ // Double_t gGap;
+ // switch(cutDefinition) {
+ // case kEtaGap01: gGap=0.1; break;
+ // case kEtaGap02: gGap=0.2; break;
+ // case kEtaGap03: gGap=0.3; break;
+ // case kEtaGap04: gGap=0.4; break;
+ // case kEtaGap05: gGap=0.5; break;
+ // default: gGap=0.0;
+ // }
+ // eta gap in tpc event plane
+ // AliDielectronVarCuts *etaGap = new AliDielectronVarCuts(AliDielectronVarManager::GetValueName(AliDielectronVarManager::kEta),"etaGap");
+ // etaGap->AddCut(AliDielectronVarManager::kEta,-1*gGap,gGap,kTRUE);
+ // die->GetEventPlanePreFilter().AddCuts(etaGap);
+
+ AliDielectronVarCuts *poi = new AliDielectronVarCuts("PoI","PoI");
+ poi->AddCut(AliDielectronVarManager::kM,2.92,3.20); // particles of interest, jpsi mass window
+ die->GetEventPlanePOIPreFilter().AddCuts(poi);
+
+ // die->SetLikeSignSubEvents();
+ die->SetPreFilterEventPlane();
}
//______________________________________________________________________________________
//
// Initialise the histograms
//
- // Bool_t hasMC=(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()!=0x0);
-
- //Setup histogram Manager
- AliDielectronHistos *histos=new AliDielectronHistos(die->GetName(),die->GetTitle());
-
-
- //add histograms to event class
- histos->AddClass("Event");
- histos->UserHistogram("Event","VtxZ","Vertex Z;z (cm)",
- 300,-15.,15.,
- AliDielectronVarManager::kZvPrim);
- histos->UserHistogram("Event","Centrality","Centrality;centrality (%)",
- "0.,5.,10.,20.,40.,50.,60.,80.,100.",
- AliDielectronVarManager::kCentrality);
- histos->UserHistogram("Event","Multiplicity","Multiplicity V0;Multiplicity V0",
- 500,0.,25000.,
- AliDielectronVarManager::kMultV0);
- histos->UserHistogram("Event","Cent_Mult","Centrality vs. Multiplicity;centrality (%);Multiplicity V0",
- 10,0.,100., 500,0.,25000.,
- AliDielectronVarManager::kCentrality,AliDielectronVarManager::kMultV0);
- histos->UserProfile("Event","Cent_Nacc",
- "accepted tracks;centrality (%)",
- AliDielectronVarManager::kNacc,
- "0.,5.,10.,20.,40.,50.,60.,80.,100.",
- AliDielectronVarManager::kCentrality);
- histos->UserProfile("Event","Cent_NVtxContrib",
- "number of vertex contributors;centrality (%)",
- AliDielectronVarManager::kNVtxContrib,
- "0.,5.,10.,20.,40.,50.,60.,80.,100.",
- AliDielectronVarManager::kCentrality);
-
-
- ////// FLOW //////
- if(cutDefinition == kTOFTRD || cutDefinition >= kEtaGap01) {
- histos->UserHistogram("Event","TPCxH2","TPC Qx component;TPCxH2",
- 100,-1500.,1500.,
- AliDielectronVarManager::kTPCxH2);
- histos->UserHistogram("Event","TPCyH2","TPC Qy component;TPCyH2",
- 100,-1500.,1500.,
- AliDielectronVarManager::kTPCyH2);
- histos->UserHistogram("Event","TPCrpH2","TPC reaction plane; #Psi^{TPC}",
- 100,-2.,2.,
- AliDielectronVarManager::kTPCrpH2);
- histos->UserHistogram("Event","TPCsub1xH2","TPC Qx component sub1;TPCsub1xH2",
- 100,-1500.,1500.,
- AliDielectronVarManager::kTPCsub1xH2);
- histos->UserHistogram("Event","TPCsub1yH2","TPC Qy component sub1;TPCsub1yH2",
- 100,-1500.,1500.,
- AliDielectronVarManager::kTPCsub1yH2);
- histos->UserHistogram("Event","TPCsub1rpH2","TPC reaction plane sub1; #Psi^{sub1}",
- 100,-2.,2.,
- AliDielectronVarManager::kTPCsub1rpH2);
- histos->UserHistogram("Event","TPCsub2xH2","TPC Qx component sub2;TPCsub2xH2",
- 100,-1500.,1500.,
- AliDielectronVarManager::kTPCsub2xH2);
- histos->UserHistogram("Event","TPCsub2yH2","TPC Qy component sub2;TPCsub2yH2",
- 100,-1500.,1500.,
- AliDielectronVarManager::kTPCsub2yH2);
- histos->UserHistogram("Event","TPCsub2rpH2","TPC reaction plane sub2; #Psi^{sub2}",
- 100,-2.,2.,
- AliDielectronVarManager::kTPCsub2rpH2);
- histos->UserHistogram("Event","TPCsub12DiffH2","TPC reaction plane diff; cos(2(#Psi^{sub1}-#Psi^{sub2}))",
- 100,-1.,1.,
- AliDielectronVarManager::kTPCsub12DiffH2);
- /* // uncorrected eventplane
- histos->UserHistogram("Event","TPCxH2uc","TPC Qx component;TPCxH2uc",
- 100,-1500.,1500.,
- AliDielectronVarManager::kTPCxH2uc);
- histos->UserHistogram("Event","TPCyH2uc","TPC Qy component;TPCyH2uc",
- 100,-1500.,1500.,
- AliDielectronVarManager::kTPCyH2uc);
- histos->UserHistogram("Event","TPCrpH2uc","TPC reaction plane;TPCrpH2uc",
- 100,-2.,2.,
- AliDielectronVarManager::kTPCrpH2uc);
- histos->UserHistogram("Event","TPCsub1xH2uc","TPC Qx component sub1;TPCsub1xH2uc",
- 100,-1500.,1500.,
- AliDielectronVarManager::kTPCsub1xH2uc);
- histos->UserHistogram("Event","TPCsub1yH2uc","TPC Qy component sub1;TPCsub1yH2uc",
- 100,-1500.,1500.,
- AliDielectronVarManager::kTPCsub1yH2uc);
- histos->UserHistogram("Event","TPCsub1rpH2uc","TPC reaction plane sub1;TPCsub1rpH2uc",
- 100,-2.,2.,
- AliDielectronVarManager::kTPCsub1rpH2uc);
- histos->UserHistogram("Event","TPCsub2xH2uc","TPC Qx component sub2;TPCsub2xH2uc",
- 100,-1500.,1500.,
- AliDielectronVarManager::kTPCsub2xH2uc);
- histos->UserHistogram("Event","TPCsub2yH2uc","TPC Qy component sub2;TPCsub2yH2uc",
- 100,-1500.,1500.,
- AliDielectronVarManager::kTPCsub2yH2uc);
- histos->UserHistogram("Event","TPCsub2rpH2uc","TPC reaction plane sub2;TPCsub2rpH2uc",
- 100,-2.,2.,
- AliDielectronVarManager::kTPCsub2rpH2uc);
- histos->UserHistogram("Event","TPCsub12DiffH2uc","TPC reaction plane difference;TPCsub12DiffH2uc",
- 100,-1.,1.,
- AliDielectronVarManager::kTPCsub12DiffH2uc);
- */
- histos->UserHistogram("Event","V0ACrpH2","VZERO-AC RP; #Psi_{2}^{V0AC} (rad.)",
- 100,-2.0,2.0,
- AliDielectronVarManager::kV0ACrpH2);
- histos->UserHistogram("Event","V0ArpH2","VZERO-A RP; #Psi_{2}^{V0A} (rad.)",
- 100,-2.0,2.0,
- AliDielectronVarManager::kV0ArpH2);
- histos->UserHistogram("Event","V0CrpH2","VZERO-C RP; #Psi_{2}^{V0C} (rad.)",
- 100,-2.0,2.0,
- AliDielectronVarManager::kV0CrpH2);
-
- histos->UserHistogram("Event","V0ATPCDiffH2","VZERO-A TPC diff; cos(2(#Psi^{V0A}-#Psi^{TPC}))",
- 300,-1.0,1.0,
- AliDielectronVarManager::kV0ATPCDiffH2);
- histos->UserHistogram("Event","V0CTPCDiffH2","VZERO-C TPC diff; cos(2(#Psi^{V0C}-#Psi^{TPC}))",
- 300,-1.0,1.0,
- AliDielectronVarManager::kV0CTPCDiffH2);
- histos->UserHistogram("Event","V0AV0CDiffH2","VZERO-A VZERO-C diff; cos(2(#Psi^{V0A}-#Psi^{V0C}))",
- 300,-1.0,1.0,
- AliDielectronVarManager::kV0AV0CDiffH2);
-
- // centrality dependent event plane histograms
- histos->UserHistogram("Event","Cent_TPCrpH2","TPC RP;centrality (%);#Psi^{TPC} (rad.)",
- 10,0.,100.,100,-2.,2.,
- AliDielectronVarManager::kCentrality,AliDielectronVarManager::kTPCrpH2);
- histos->UserHistogram("Event","Cent_TPCsub1rpH2","TPC-1 RP;centrality (%);#Psi^{sub1} (rad.)",
- 10,0.,100.,100,-2.,2.,
- AliDielectronVarManager::kCentrality,AliDielectronVarManager::kTPCsub1rpH2);
- histos->UserHistogram("Event","Cent_TPCsub2rpH2","TPC-2 RP;centrality (%);#Psi^{sub2} (rad.)",
- 10,0.,100.,100,-2.,2.,
- AliDielectronVarManager::kCentrality,AliDielectronVarManager::kTPCsub2rpH2);
-
- histos->UserHistogram("Event","Cent_V0ACrpH2","VZERO-AC RP;centrality (%);#Psi_{2}^{V0AC} (rad.)",
- 10,0.,100.,100,-2.0,2.0,
- AliDielectronVarManager::kCentrality,AliDielectronVarManager::kV0ACrpH2);
- histos->UserHistogram("Event","Cent_V0ArpH2","VZERO-A RP;centrality (%);#Psi_{2}^{V0A} (rad.)",
- 10,0.,100.,100,-2.0,2.0,
- AliDielectronVarManager::kCentrality,AliDielectronVarManager::kV0ArpH2);
- histos->UserHistogram("Event","Cent_V0CrpH2","VZERO-C RP;centrality (%);#Psi_{2}^{V0C} (rad.)",
- 10,0.,100.,100,-2.0,2.0,
- AliDielectronVarManager::kCentrality,AliDielectronVarManager::kV0CrpH2);
- // for resolution calculation
- histos->UserHistogram("Event","Cent_V0ATPCDiffH2","VZERO-A TPC diff;centrality (%);cos(2(#Psi^{V0A}-#Psi^{TPC}))",
- 10,0.,100.,300,-1.0,1.0,
- AliDielectronVarManager::kCentrality,AliDielectronVarManager::kV0ATPCDiffH2);
- histos->UserHistogram("Event","Cent_V0CTPCDiffH2","VZERO-C TPC diff;centrality (%);cos(2(#Psi^{V0C}-#Psi^{TPC}))",
- 10,0.,100.,300,-1.0,1.0,
- AliDielectronVarManager::kCentrality,AliDielectronVarManager::kV0CTPCDiffH2);
- histos->UserHistogram("Event","Cent_V0AV0CDiffH2","VZERO-A VZERO-C diff;centrality (%);cos(2(#Psi^{V0A}-#Psi^{V0C}))",
- 10,0.,100.,300,-1.0,1.0,
- AliDielectronVarManager::kCentrality,AliDielectronVarManager::kV0AV0CDiffH2);
- histos->UserHistogram("Event","Cent_TPCsub12DiffH2","TPC-sub1 TPC-sub2 diff;centrality (%);cos(2(#Psi^{sub1}-#Psi^{sub2}))",
- 10,0.,100.,300,-1.0,1.0,
- AliDielectronVarManager::kCentrality,AliDielectronVarManager::kTPCsub12DiffH2);
- // detector effects
- histos->UserHistogram("Event","Cent_TPCsub12DiffH2Sin","TPC-sub1 TPC-sub2 diff;centrality (%);sin(2(#Psi^{sub1}-#Psi^{sub2}))",
- 10,0.,100.,300,-1.0,1.0,AliDielectronVarManager::kCentrality,AliDielectronVarManager::kTPCsub12DiffH2Sin);
- //// EPSelectionTask V0 information directly from the ESDs
- histos->UserHistogram("Event","Cent_v0ACrpH2","VZERO-AC RP;centrality (%);#Psi_{2}^{v0AC} (rad.)",
- 10,0.,100.,100,-2.0,2.0,
- AliDielectronVarManager::kCentrality,AliDielectronVarManager::kv0ACrpH2);
- histos->UserHistogram("Event","Cent_v0ArpH2","VZERO-A RP;centrality (%);#Psi_{2}^{v0A} (rad.)",
- 10,0.,100.,100,-2.0,2.0,
- AliDielectronVarManager::kCentrality,AliDielectronVarManager::kv0ArpH2);
- histos->UserHistogram("Event","Cent_v0CrpH2","VZERO-C RP;centrality (%);#Psi_{2}^{v0C} (rad.)",
- 10,0.,100.,100,-2.0,2.0,
- AliDielectronVarManager::kCentrality,AliDielectronVarManager::kv0CrpH2);
- // for resolution calculation
- histos->UserHistogram("Event","Cent_v0ATPCDiffH2","VZERO-A TPC diff;centrality (%);cos(2(#Psi^{v0A}-#Psi^{TPC}))",
- 10,0.,100.,300,-1.0,1.0,
- AliDielectronVarManager::kCentrality,AliDielectronVarManager::kv0ATPCDiffH2);
- histos->UserHistogram("Event","Cent_v0CTPCDiffH2","VZERO-C TPC diff;centrality (%);cos(2(#Psi^{v0C}-#Psi^{TPC}))",
- 10,0.,100.,300,-1.0,1.0,
- AliDielectronVarManager::kCentrality,AliDielectronVarManager::kv0CTPCDiffH2);
- histos->UserHistogram("Event","Cent_v0Av0CDiffH2","VZERO-A VZERO-C diff;centrality (%);cos(2(#Psi^{v0A}-#Psi^{v0C}))",
- 10,0.,100.,300,-1.0,1.0,
- AliDielectronVarManager::kCentrality,AliDielectronVarManager::kv0Av0CDiffH2);
-
+ Bool_t hasMC=die->GetHasMC();
+
+ // booleans for histo selection
+ Bool_t bHistTrackQA=kFALSE, bHistEvts=kFALSE, bHistPair=kFALSE, bHistTrk=kFALSE, bHistPairME=kFALSE, bHistFlow=kFALSE, bHistFlowQA=kFALSE, bHistPID=kFALSE;
+ switch (cutDefinition) {
+ case kNoCut:
+ case kEvent: bHistEvts=kTRUE; break;
+ case kPIDqa: bHistTrk=kTRUE; break;
+ case kNoPID: bHistTrk=kTRUE; break;
+ case kRec: /* */ break;
+ case kTPC:
+ case kTOF: //bHistEvts=kTRUE; //bHistFlow=kTRUE;
+ case kTRD: ///bHistEvts=kTRUE; //bHistFlow=kTRUE;
+ case kAvgPt:
+ case kSysMC:
+ case kStd: ///*bHistEvts=kTRUE;*/ bHistPair=kTRUE; break; //bHistPairME=kTRUE;
+ case kTPCTOFTRD: bHistPair=kTRUE; bHistTrk=kTRUE; bHistPID=kTRUE; break; //bHistPairME=kTRUE;
+ case kSysMCele: bHistPair=kTRUE; bHistPID=kFALSE; break; //bHistPairME=kTRUE;
+ case kLegEff: bHistTrk=kTRUE; break;
}
-
- ////// MONTE CARLO //////
- /*
- if(cutDefinition == kTOFTRD && hasMC) {
- histos->AddClass("MCEvent");
- histos->UserHistogram("MCEvent","Cent_NJPsis","Centrality vs. generated incl. J/#psi per event;centrality (%);N_{J/#psi}",
- 10,0.,100., 21,-0.5,20.5,
- AliDielectronVarManager::kCentrality,AliDielectronVarManager::kNumberOfJPsis);
+ if(hasMC) {
+ bHistPID=kFALSE;
}
- */
-
-
+
+ //Setup histogram Manager
+ AliDielectronHistos *histos=new AliDielectronHistos(die->GetName(),die->GetTitle());
+ die->SetHistogramManager(histos);
+
//Initialise histogram classes
histos->SetReservedWords("Track;Pair");
-
-
- //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)));
+
+ //add histograms to event class
+ histos->AddClass("Event");
+ Int_t maxMixBins = (die->GetMixingHandler() ? die->GetMixingHandler()->GetNumberOfBins() : 0);
+ histos->UserHistogram("Event","","", 100, 0.0, 100.0, AliDielectronVarManager::kCentrality);
+ if(die->GetMixingHandler() && maxMixBins)
+ histos->UserHistogram("Event","","", maxMixBins, 0, maxMixBins, AliDielectronVarManager::kMixingBin);
+ // candidates monitoring
+ histos->UserProfile("Event","","", AliDielectronVarManager::kTracks, GetRunNumbers2011(), AliDielectronVarManager::kRunNumber);
+ histos->UserProfile("Event","","", AliDielectronVarManager::kPairs, GetRunNumbers2011(), AliDielectronVarManager::kRunNumber);
+
+
+ //event plane histograms
+ if(cutDefinition==kFlow) {
+ if(!hasMC) AddQAHistsEP(die);
+ histos->UserHistogram("Pair","","", 100,-1*TMath::Pi(),+1*TMath::Pi(), AliDielectronVarManager::kDeltaPhiv0ArpH2);
+ histos->UserHistogram("Pair","","", 100,-1*TMath::Pi(),+1*TMath::Pi(), AliDielectronVarManager::kDeltaPhiv0CrpH2);
+ histos->UserHistogram("Pair","","", 100,-1*TMath::Pi(),+1*TMath::Pi(), AliDielectronVarManager::kDeltaPhiTPCrpH2);
}
-
- if(cutDefinition < kEtaGap01) {
-
- //legs from pair
- for (Int_t i=0; i<3; ++i){
- histos->AddClass(Form("Track_Legs_%s",AliDielectron::PairClassName(i)));
+
+ ////// EVENT HISTOS /////
+ if(bHistEvts) {
+ histos->UserHistogram("Event","","", GetRunNumbers2011(), AliDielectronVarManager::kRunNumber);
+ histos->UserHistogram("Event","","", 300,-15.,15., AliDielectronVarManager::kZvPrim);
+ // histos->UserHistogram("Event","","", 300,0.,15000., AliDielectronVarManager::kRefMult);
+ histos->UserHistogram("Event","","", 300,0.,3000., AliDielectronVarManager::kRefMultTPConly);
+
+ // histos->UserHistogram("Event","","", 100,0.,100., AliDielectronVarManager::kRefMultTPConly);
+ histos->UserHistogram("Event","","", GetRunNumbers2011(), AliDielectronHelper::MakeLinBinning(3000, 0., 3000.),
+ AliDielectronVarManager::kRunNumber, AliDielectronVarManager::kRefMultTPConly);
+ histos->UserHistogram("Event","","", AliDielectronHelper::MakeLinBinning(90, 0., 90.), AliDielectronHelper::MakeLinBinning(3000, 0., 3000.),
+ AliDielectronVarManager::kCentrality, AliDielectronVarManager::kRefMultTPConly);
+ histos->UserProfile( "Event","","", AliDielectronVarManager::kRefMultTPConly, 90, 0., 90., AliDielectronVarManager::kCentrality);
+ histos->UserProfile( "Event","","", AliDielectronVarManager::kMultV0A, 90, 0., 90., AliDielectronVarManager::kCentrality);
+ histos->UserProfile( "Event","","", AliDielectronVarManager::kMultV0C, 90, 0., 90., AliDielectronVarManager::kCentrality);
+ // histos->UserProfile( "Event","","", AliDielectronVarManager::kEqMultV0A, 90, 0., 90., AliDielectronVarManager::kCentrality);
+ // histos->UserProfile( "Event","","", AliDielectronVarManager::kEqMultV0C, 90, 0., 90., AliDielectronVarManager::kCentrality);
+ histos->UserProfile( "Event","","", AliDielectronVarManager::kNVtxContrib, 90, 0., 90., AliDielectronVarManager::kCentrality);
+ histos->UserHistogram("Event","","", 100, 0., 4000., AliDielectronVarManager::kNVtxContribTPC);
+
+ // event plane histograms
+ if(!hasMC) AddQAHistsEP(die);
+ // trigger histograms
+ if(!hasMC) {
+ histos->UserHistogram("Event","","", 29,0.,29., AliDielectronVarManager::kTriggerInclONL);
+ histos->UserHistogram("Event","","", 29,0.,29., AliDielectronVarManager::kTriggerInclOFF);
+ histos->UserHistogram("Event","","", 29,0.,29., AliDielectronVarManager::kTriggerExclOFF);
+ histos->UserHistogram("Event","","", 100,0.,100.,29,0.,29.,
+ AliDielectronVarManager::kCentrality,AliDielectronVarManager::kTriggerInclONL);
+ histos->UserHistogram("Event","","", 100,0.,100.,29,0.,29.,
+ AliDielectronVarManager::kCentrality,AliDielectronVarManager::kTriggerExclOFF);
+ histos->UserHistogram("Event","","", 300,0.,3000.,29,0.,29.,
+ AliDielectronVarManager::kRefMultTPConly,AliDielectronVarManager::kTriggerExclOFF);
+ histos->UserHistogram("Event","","", 300,0.,3000.,29,0.,29.,
+ AliDielectronVarManager::kRefMultTPConly,AliDielectronVarManager::kTriggerInclOFF);
}
-
+
+ } //hist: event
+
+ ///// PAIR HISTOS /////
+ if(bHistPair) {
+
+ //Pair classes
+ // to fill also mixed event histograms loop until 7 or 10
+ for (Int_t i=1; i<(bHistPairME ? 8 : 2); ++i){
+ histos->AddClass(Form("Pair_%s",AliDielectron::PairClassName(i)));
+ }
+
+ //add MC signal histograms to pair class
+ if(die->GetMCSignals()) {
+ for (Int_t i=0; i<die->GetMCSignals()->GetEntriesFast(); ++i) {
+ histos->AddClass(Form("Pair_%s",die->GetMCSignals()->At(i)->GetName()));
+ histos->AddClass(Form("Pair_%s_MCtruth",die->GetMCSignals()->At(i)->GetName()));
+ }
+ }
+
+ ///// Pair classes /////
+ histos->UserHistogram("Pair","","", 375,.0,375*0.04, AliDielectronVarManager::kM); // 40MeV bins, 12GeV/c2
+ histos->UserHistogram("Pair","","", 100,-1.,1., AliDielectronVarManager::kY);
+ histos->UserHistogram("Pair","","", 400,0,20., AliDielectronVarManager::kPt);
+ histos->UserHistogram("Pair","","", 125,.0,125*0.04, 400,0,20., AliDielectronVarManager::kM,AliDielectronVarManager::kPt);
+ histos->UserHistogram("Pair","","", 100,0.,3.15, AliDielectronVarManager::kOpeningAngle);
+ histos->UserHistogram("Pair","","", 100,0.,20, AliDielectronVarManager::kChi2NDF);
+ histos->UserHistogram("Pair","","", 100,0.,3.15, AliDielectronVarManager::kPsiPair);
+ histos->UserHistogram("Pair","","", 200,0.,100., AliDielectronVarManager::kR);
+ histos->UserHistogram("Pair","","", 50,0.,5., AliDielectronVarManager::kLegDist);
+ histos->UserHistogram("Pair","","", 50,0.,5., AliDielectronVarManager::kLegDistXY);
+ // histos->UserHistogram("Pair","","", 210,-1.05,1.05, 100,0.,2.5,
+ // AliDielectronVarManager::kArmAlpha,AliDielectronVarManager::kArmPt);
+ histos->UserHistogram("Pair","","", 500,-2.5,2.5, AliDielectronVarManager::kImpactParXY);
+ histos->UserHistogram("Pair","","", 600,-3.,3., AliDielectronVarManager::kImpactParZ);
+
+ if(!hasMC && die->GetMixingHandler() && maxMisBins)
+ histos->UserHistogram("Pair","","", maxMixBins, 0, maxMixBins, AliDielectronVarManager::kMixingBin);
+
+ histos->UserHistogram("Pair","","", 100,.0,1., AliDielectronVarManager::kRndmPair);
+ // histos->UserHistogram("Pair","","", GetPDGcodes(), AliDielectronVarManager::kPdgCode);
+ // histos->UserHistogram("Pair","","", GetPDGcodes(), AliDielectronVarManager::kPdgCodeMother);
+ // histos->UserHistogram("Pair","","", GetPDGcodes(), AliDielectronVarManager::kPdgCodeGrandMother);
+
+ } //hist: pairs
+
+
+ ///// TRACK HISTOS /////
+ if(bHistTrk) {
//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)));
+ //legs from pair (fill SE)
+ for (Int_t i=1; i<2; ++i){
+ if(bHistPair) histos->AddClass(Form("Track_Legs_%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 Track classes
- //histos->UserHistogram("Track","TOFbit","TOFbit;bit;#tracks",19,-9.5,9.5,AliDielectronVarManager::kTOFPIDBit);
- histos->UserHistogram("Track","Pt","Pt;Pt [GeV];#tracks",
- 400,0,20.,
- AliDielectronVarManager::kPt);
- histos->UserHistogram("Track","TPCnCls","Number of Clusters TPC;TPC number clusteres;#tracks",
- 160,-0.5,159.5,
- AliDielectronVarManager::kNclsTPC);
- histos->UserHistogram("Track","TPCsignalN","Number of Clusters TPC;TPC number clusteres;#tracks",
- 160,-0.5,159.5,
- AliDielectronVarManager::kTPCsignalN);
-
- histos->UserHistogram("Track","dXY","dXY;dXY [cm];#tracks",
- 500,-1.,1.,
- AliDielectronVarManager::kImpactParXY);
- histos->UserHistogram("Track","dZ","dZ;dZ [cm];#tracks",
- 600,-3.,3.,
- AliDielectronVarManager::kImpactParZ);
- histos->UserHistogram("Track","Eta_Phi","Eta Phi Map; Eta; Phi;#tracks",
- 200,-1,1,200,0,6.285,
- AliDielectronVarManager::kEta,AliDielectronVarManager::kPhi);
-
- histos->UserHistogram("Track","dEdx_P","dEdx;P [GeV];TPC signal (arb units);#tracks",
- 400,0.2,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.2,20.,200,-10.,10.,
- AliDielectronVarManager::kPIn,AliDielectronVarManager::kTPCnSigmaEle,kTRUE);
-
- histos->UserHistogram("Track","Ncl",";Number clusters TPC;Number clusters TPC",
- 160,-0.5,159.5,
- AliDielectronVarManager::kNclsTPC);
- histos->UserHistogram("Track","NclFr",";Number of findable clusters (robust);Number findable clusters TPC",
- 160,-0.5,159.5,
- AliDielectronVarManager::kNFclsTPCr);
- histos->UserHistogram("Track","Ncl_NclFr","Number of (findable) clusters TPC;found clusters;findable clusters",
- 160,-0.5,159.5,160,-0.5,159.5,
- AliDielectronVarManager::kNclsTPC,AliDielectronVarManager::kNFclsTPCr);
- histos->UserHistogram("Track","NtrklTRD",";Number tracklets TRD for pid;Number tracklets TRD",
- 8,-0.5,7.5,
- AliDielectronVarManager::kTRDpidQuality);
-
- //add histograms to Pair classes
- histos->UserHistogram("Pair","InvMass","Inv.Mass;Inv. Mass [GeV];#pairs",
- 300,.0,300*0.04, AliDielectronVarManager::kM); // 40MeV bins, 12GeV/c2
- histos->UserHistogram("Pair","Rapidity","Rapidity;Rapidity;#pairs",
- 100,-1.,1., AliDielectronVarManager::kY);
- histos->UserHistogram("Pair","OpeningAngle","Opening angle;angle",
- 100,0.,3.15, AliDielectronVarManager::kOpeningAngle);
- histos->UserHistogram("Pair","Chi2NDF","#chi^{2}/NDF;#chi^{2}/NDF",
- 100,0.,20, AliDielectronVarManager::kChi2NDF);
- histos->UserHistogram("Pair","PsiPair","PsiPair;#psi",
- 100,0.,3.15, AliDielectronVarManager::kPsiPair);
- histos->UserHistogram("Pair","R","Radius;r (cm)",
- 200,0.,100., AliDielectronVarManager::kR);
- histos->UserHistogram("Pair","LegDist",";dca (cm)",
- 50,0.,5., AliDielectronVarManager::kLegDist);
- histos->UserHistogram("Pair","LegDistXY",";dca_{xy} (cm)",
- 50,0.,5., AliDielectronVarManager::kLegDistXY);
-
-
+ //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)));
+ if(cutDefinition!=kLegEff)
+ histos->AddClass(Form("Track_%s", AliDielectron::PairClassName(AliDielectron::kEv1PM)));
+
+ // PID post calibration
+ if(cutDefinition==kPIDqa) AddQAHistsPID(die);
+
+ // Vertex
+ // histos->UserHistogram("Track","","", 500,-1.,1., AliDielectronVarManager::kImpactParXY);
+ // histos->UserHistogram("Track","","", 600,-3.,3., AliDielectronVarManager::kImpactParZ);
+ // Kinematics
+ histos->UserHistogram("Track","","", 400,0,20., AliDielectronVarManager::kPt);
+ histos->UserHistogram("Track","","", 200,-1,1, 200,0,6.285, AliDielectronVarManager::kEta,AliDielectronVarManager::kPhi);
+ // TPC
+ // histos->UserHistogram("Track","","", 160,-0.5,159.5, AliDielectronVarManager::kNclsTPC);
+ // histos->UserHistogram("Track","","", 160,-0.5,159.5, AliDielectronVarManager::kTPCsignalN);
+ // histos->UserHistogram("Track","","", 160,-0.5,159.5, AliDielectronVarManager::kNFclsTPCr);
+ // histos->UserHistogram("Track","","", 160,-0.5,159.5, 160,-0.5,159.5, AliDielectronVarManager::kNclsTPC,AliDielectronVarManager::kNFclsTPCr);
+ // TRD
+ histos->UserHistogram("Track","","", 8,-0.5, 7.5, AliDielectronVarManager::kTRDpidQuality);
+ histos->UserHistogram("Track","","", 101,-0.5, 100.5, AliDielectronVarManager::kTRDsignal);
+ histos->UserHistogram("Track","","", 101,-0.5, 100.5, AliDielectronVarManager::kTRDchi2);
+ // PID
+ histos->UserHistogram("Track","","", 400,0.2,20.,200,0.,200., AliDielectronVarManager::kPIn,AliDielectronVarManager::kTPCsignal,kTRUE);
+ histos->UserHistogram("Track","","", 400,0.2,20.,200,-10.,10.,AliDielectronVarManager::kPIn,AliDielectronVarManager::kTPCnSigmaEle,kTRUE);
+ histos->UserHistogram("Track","","", 250,0.2,5.0,300,0.,1.2, AliDielectronVarManager::kPIn,AliDielectronVarManager::kTOFbeta,kTRUE);
+ // histos->UserHistogram("Track","","", 100,-1.,+1.,200,-10.,10.,AliDielectronVarManager::kEta,AliDielectronVarManager::kTPCnSigmaEle);
+
+ // histos->UserHistogram("Track","","", 100,0.2,10.,100,0.,200.,
+ // AliDielectronVarManager::kPIn,AliDielectronVarManager::kTPCsignal,kTRUE);
+ // histos->UserHistogram("Track","","", 100,0.2,10.,100,-5.,+5.,
+ // AliDielectronVarManager::kPIn,AliDielectronVarManager::kTPCnSigmaEle,kTRUE);
+ histos->UserHistogram("Track","","", 100,0.,4000.,100,0.,200.,
+ AliDielectronVarManager::kRefMultTPConly,AliDielectronVarManager::kTPCsignal);
+ histos->UserHistogram("Track","","", 100,0.,4000.,100,-5.,+5.,
+ AliDielectronVarManager::kRefMultTPConly,AliDielectronVarManager::kTPCnSigmaEle);
+ histos->UserHistogram("Track","","", 100,0.,100.,100,0.,200.,
+ AliDielectronVarManager::kCentrality,AliDielectronVarManager::kTPCsignal);
+ histos->UserHistogram("Track","","", 100,0.,100.,100,-5.,+5.,
+ AliDielectronVarManager::kCentrality,AliDielectronVarManager::kTPCnSigmaEle);
+ histos->UserHistogram("Track","","", 100,-1.,+1.,100,0.,200.,
+ AliDielectronVarManager::kEta,AliDielectronVarManager::kTPCsignal);
+ histos->UserHistogram("Track","","", 100,-1.,+1.,100,-5.,+5.,
+ AliDielectronVarManager::kEta,AliDielectronVarManager::kTPCnSigmaEle);
+ histos->UserHistogram("Track","","", 100,-1.,+1.,100,-5.,+5.,
+ AliDielectronVarManager::kEta,AliDielectronVarManager::kTPCnSigmaEleRaw);
+ histos->UserHistogram("Track","","", 100,-1.,+1.,100,-5.,+5.,
+ AliDielectronVarManager::kEta,AliDielectronVarManager::kTPCnSigmaPio);
+ } //hist: tracks
+
+
+ if(cutDefinition==kAvgPt) {
+ // add single electron efficiency histograms
+ AddQAHistsEff(die);
}
-
- //// FLOW results use tprofiles
- if(cutDefinition == kTOFTRD || cutDefinition == kTOFTRD2 || cutDefinition >= kEtaGap01) {
-
- histos->UserProfile("Pair","M_Cent_Pt_V0ACrpH2FlowV2",
- "cos(2(#varphi-#Psi^{V0AC}));mass (GeV/c^{2});centrality (%);p_{T} (GeV/c)",
- AliDielectronVarManager::kV0ACrpH2FlowV2,
- 125,0.,125*.04, 10, 0.,100., 200,0.,100.,
- AliDielectronVarManager::kM, AliDielectronVarManager::kCentrality, AliDielectronVarManager::kPt);
-
- histos->UserProfile("Pair","M_Cent_Pt_V0ArpH2FlowV2",
- "cos(2(#varphi-#Psi^{V0A}));mass (GeV/c^{2});centrality (%);p_{T} (GeV/c)",
- AliDielectronVarManager::kV0ArpH2FlowV2,
- 125,0.,125*.04, 10, 0.,100., 200,0.,100.,
- AliDielectronVarManager::kM, AliDielectronVarManager::kCentrality, AliDielectronVarManager::kPt);
-
- histos->UserProfile("Pair","M_Cent_Pt_V0CrpH2FlowV2",
- "cos(2(#varphi-#Psi^{V0C}));mass (GeV/c^{2});centrality (%);p_{T} (GeV/c)",
- AliDielectronVarManager::kV0CrpH2FlowV2,
- 125,0.,125*.04, 10, 0.,100., 200,0.,100.,
- AliDielectronVarManager::kM, AliDielectronVarManager::kCentrality, AliDielectronVarManager::kPt);
-
- histos->UserProfile("Pair","M_Cent_Pt_v0ACrpH2FlowV2",
- "cos(2(#varphi-#Psi^{V0AC}));mass (GeV/c^{2});centrality (%);p_{T} (GeV/c)",
- AliDielectronVarManager::kv0ACrpH2FlowV2,
- 125,0.,125*.04, 10, 0.,100., 200,0.,100.,
- AliDielectronVarManager::kM, AliDielectronVarManager::kCentrality, AliDielectronVarManager::kPt);
-
- histos->UserProfile("Pair","M_Cent_Pt_v0ArpH2FlowV2",
- "cos(2(#varphi-#Psi^{V0A}));mass (GeV/c^{2});centrality (%);p_{T} (GeV/c)",
- AliDielectronVarManager::kv0ArpH2FlowV2,
- 125,0.,125*.04, 10, 0.,100., 200,0.,100.,
- AliDielectronVarManager::kM, AliDielectronVarManager::kCentrality, AliDielectronVarManager::kPt);
-
- histos->UserProfile("Pair","M_Cent_Pt_v0CrpH2FlowV2",
- "cos(2(#varphi-#Psi^{V0C}));mass (GeV/c^{2});centrality (%);p_{T} (GeV/c)",
- AliDielectronVarManager::kv0CrpH2FlowV2,
- 125,0.,125*.04, 10, 0.,100., 200,0.,100.,
- AliDielectronVarManager::kM, AliDielectronVarManager::kCentrality, AliDielectronVarManager::kPt);
-
- // 1D
- histos->UserProfile("Pair","M_V0ACrpH2FlowV2",
- "cos(2(#varphi-#Psi^{V0AC}));mass (GeV/c^{2})",
- AliDielectronVarManager::kV0ACrpH2FlowV2,
- 125,0.,125*.04,
- AliDielectronVarManager::kM);
-
- histos->UserProfile("Pair","M_V0ArpH2FlowV2",
- "cos(2(#varphi-#Psi^{V0AC}));mass (GeV/c^{2})",
- AliDielectronVarManager::kV0ArpH2FlowV2,
- 125,0.,125*.04,
- AliDielectronVarManager::kM);
-
- histos->UserProfile("Pair","M_V0CrpH2FlowV2",
- "cos(2(#varphi-#Psi^{V0AC}));mass (GeV/c^{2})",
- AliDielectronVarManager::kV0CrpH2FlowV2,
- 125,0.,125*.04,
- AliDielectronVarManager::kM);
- }
-
+
+ ////// MONTE CARLO //////
+ if(cutDefinition==kLegEff) {
+
+ //add MC signal histograms to track class
+ if(die->GetMCSignals()) {
+
+ TString className="";
+ for (Int_t i=0; i<die->GetMCSignals()->GetEntriesFast(); ++i) {
+ TString sigMCname = die->GetMCSignals()->At(i)->GetName();
+
+ // mc truth
+ // histos->AddClass(Form("Pair_%s_MCtruth", sigMCname.Data()));
+ // histos->AddClass(Form("Track_Legs_%s_MCtruth",sigMCname.Data()));
+ // mc reconstructed
+ // histos->AddClass(Form("Pair_%s", sigMCname.Data()));
+ // histos->AddClass(Form("Track_Legs_%s", sigMCname.Data()));
+ // tracks
+ histos->AddClass(Form("Track_%s_%s",AliDielectron::PairClassName(AliDielectron::kEv1PM),sigMCname.Data()));
+ histos->AddClass(Form("Track_%s_%s_MCtruth",AliDielectron::PairClassName(AliDielectron::kEv1PM),sigMCname.Data()));
+ } //end: loop signals
+ } //end: has signals
+
+ // add single electron histograms
+ AddHistsEleEff(die);
+ // pdg codes
+ // histos->UserHistogram("Track","","", GetPDGcodes(), AliDielectronVarManager::kPdgCode);
+ // histos->UserHistogram("Track","","", GetPDGcodes(), AliDielectronVarManager::kPdgCodeMother);
+ // histos->UserHistogram("Track","","", GetPDGcodes(), AliDielectronVarManager::kPdgCodeGrandMother);
+ // histos->UserHistogram("Track","","", GetPDGcodes(), GetPDGcodes(),
+ // AliDielectronVarManager::kPdgCodeMother, AliDielectronVarManager::kPdgCodeGrandMother);
+ }
+
+
+
die->SetHistogramManager(histos);
+
+
+ ////// LOG //////
+ TIter nextClass(histos->GetHistogramList());
+ THashList *l=0;
+ while ( (l=static_cast<THashList*>(nextClass())) ) {
+ //printf(" [D] HistogramManger: Class %s: Histograms: %04d \n", l->GetName(), l->GetEntries());
+ }
+
+} //end: init histograms
+
+void AddQAHistsPID(AliDielectron *die) {
+ //
+ // add histograms for PID validation, comparison, post calibration aso.
+ //
+
+ Bool_t hasMC=die->GetHasMC();
+ AliDielectronHistos *histos = die->GetHistoManager();
+
+ // add MC signal tracks
+ if(hasMC && die->GetMCSignals()) {
+ TString sigMCname = die->GetMCSignals()->Last()->GetName();
+ histos->AddClass(Form("Track_%s_%s",AliDielectron::PairClassName(AliDielectron::kEv1PM),sigMCname.Data()));
+ }
+
+
+ // for TPC post calibration
+
+ // arbitrary binning for variables
+ TVectorD *vcen = AliDielectronHelper::MakeLinBinning( 11, 0., 55.);
+ (*vcen)[21] = 90.;
+
+ // array of bin limits
+ TObjArray *limits = new TObjArray();
+ limits->Add(AliDielectronHelper::MakeLinBinning( 75,-10., 5.));
+ limits->Add(AliDielectronHelper::MakeLinBinning( 75,-10., 5.));
+ // limits->Add(AliDielectronHelper::MakeLinBinning(100, 0., 200.));
+ // limits->Add(AliDielectronHelper::MakeLinBinning(60, 0., 1.2));
+ limits->Add(AliDielectronHelper::MakeLinBinning( 50, 0., 4000.));
+ limits->Add(AliDielectronHelper::MakeLinBinning( 50, 0., 10.));
+ // limits->Add(vcen);
+ // limits->Add(AliDielectronHelper::MakeLinBinning( 36, 0., 90.));
+ limits->Add(AliDielectronHelper::MakeLinBinning( 20, -1., 1.));
+ limits->Add(GetRunNumbers2011());
+ UInt_t var[]={AliDielectronVarManager::kTPCnSigmaEleRaw, // NOTE: raw nsigma w/o corr
+ AliDielectronVarManager::kTPCnSigmaEle,
+ // AliDielectronVarManager::kTOFbeta,
+ // AliDielectronVarManager::kTPCsignal,
+ AliDielectronVarManager::kRefMultTPConly,
+ AliDielectronVarManager::kPIn,
+ // AliDielectronVarManager::kCentrality,
+ AliDielectronVarManager::kEta,
+ AliDielectronVarManager::kRunNumber
+ };
+ // add merged track histogram
+ histos->UserSparse("Track", limits->GetEntriesFast(), limits, var);
+
+ // run dependence of nisgma electron
+ histos->UserHistogram("Track","","", GetRunNumbers2011(),AliDielectronHelper::MakeLinBinning(100,-5.,+5.),
+ AliDielectronVarManager::kRunNumber,AliDielectronVarManager::kTPCnSigmaEleRaw);
+ histos->UserHistogram("Track","","", GetRunNumbers2011(),AliDielectronHelper::MakeLinBinning(100,-5.,+5.),
+ AliDielectronVarManager::kRunNumber,AliDielectronVarManager::kTPCnSigmaEle);
+
+ // post calibration check
+ histos->UserHistogram("Track","","", 100,0.2,10.,100,-5.,+5.,
+ AliDielectronVarManager::kPIn,AliDielectronVarManager::kTPCnSigmaEleRaw,kTRUE);
+ histos->UserHistogram("Track","","", 100,0.2,10.,100,-5.,+5.,
+ AliDielectronVarManager::kPIn,AliDielectronVarManager::kTPCnSigmaEle,kTRUE);
+
+ // TPC number of points used for nsigma calculation
+ histos->UserHistogram("Track","","", 160,-0.5,159.5, AliDielectronVarManager::kTPCsignalN);
}
+void AddQAHistsEP(AliDielectron *die) {
+ //
+ // add histograms for event plane flattening aso.
+ //
+
+ Bool_t hasMC=die->GetHasMC();
+ AliDielectronHistos *histos = die->GetHistoManager();
+
+ histos->UserHistogram("Event","","", AliDielectronHelper::MakeLinBinning(100,-1.6,1.6), AliDielectronVarManager::kTPCrpH2);
+ histos->UserHistogram("Event","","", AliDielectronHelper::MakeLinBinning(100,-1.6,1.6), AliDielectronVarManager::kTPCrpH2uc);
+
+ // event plane resolutions
+ TObjArray *limits = new TObjArray();
+ limits->Add(AliDielectronHelper::MakeLinBinning( 18, 0., 90.));
+ // limits->Add(GetRunNumbers2011());
+ limits->Add(AliDielectronHelper::MakeLinBinning(100, -1., 1.));
+ limits->Add(AliDielectronHelper::MakeLinBinning(100, -1., 1.));
+ limits->Add(AliDielectronHelper::MakeLinBinning(100, -1., 1.));
+
+ UInt_t var[]={AliDielectronVarManager::kCentrality,
+ // AliDielectronVarManager::kRunNumber,
+ AliDielectronVarManager::kv0ATPCDiffH2,
+ AliDielectronVarManager::kv0CTPCDiffH2,
+ AliDielectronVarManager::kv0Av0CDiffH2 };
+ // histos->UserSparse("Event", 4, limits, var);
+
+ // event plane angles
+ TObjArray *limits2 = new TObjArray();
+ limits2->Add(AliDielectronHelper::MakeLinBinning( 18, 0., 90.));
+ // limits2->Add(GetRunNumbers2011());
+ limits2->Add(AliDielectronHelper::MakeLinBinning(100, -1.6, 1.6));
+ limits2->Add(AliDielectronHelper::MakeLinBinning(100, -1.6, 1.6));
+ limits2->Add(AliDielectronHelper::MakeLinBinning(100, -1.6., 1.6));
+
+ UInt_t var2[]={AliDielectronVarManager::kCentrality,
+ // AliDielectronVarManager::kRunNumber,
+ AliDielectronVarManager::kv0ArpH2,
+ AliDielectronVarManager::kv0CrpH2,
+ AliDielectronVarManager::kTPCrpH2 };
+ histos->UserSparse("Event", limits2->GetEntriesFast(), limits2, var2);
+
+ // VZERO event plane angles (sub rings)
+ histos->UserHistogram("Event","","", 100,-2.,2., AliDielectronVarManager::kv0A0rpH2);
+ histos->UserHistogram("Event","","", 100,-2.,2., AliDielectronVarManager::kv0C0rpH2);
+ histos->UserHistogram("Event","","", 100,-2.,2., AliDielectronVarManager::kv0A3rpH2);
+ histos->UserHistogram("Event","","", 100,-2.,2., AliDielectronVarManager::kv0C3rpH2);
+ histos->UserHistogram("Event","","", 100,-2.,2., AliDielectronVarManager::kv0ACrpH2); // combined A+C
+
+ // run dependence of the angles
+ histos->UserHistogram("Event","","", GetRunNumbers2011(), AliDielectronHelper::MakeLinBinning(100,-1.6,1.6),
+ AliDielectronVarManager::kRunNumber, AliDielectronVarManager::kv0ArpH2);
+ histos->UserHistogram("Event","","", GetRunNumbers2011(), AliDielectronHelper::MakeLinBinning(100,-1.6,1.6),
+ AliDielectronVarManager::kRunNumber, AliDielectronVarManager::kv0CrpH2);
+ histos->UserHistogram("Event","","", GetRunNumbers2011(), AliDielectronHelper::MakeLinBinning(100,-1.6,1.6),
+ AliDielectronVarManager::kRunNumber, AliDielectronVarManager::kTPCrpH2);
+
+
+ // TPC q vector components
+ histos->UserHistogram("Event","","", 100,-1500.,1500., AliDielectronVarManager::kTPCxH2);
+ histos->UserHistogram("Event","","", 100,-1500.,1500., AliDielectronVarManager::kTPCyH2);
+ // histos->UserHistogram("Event","","", 100, -2., 2., AliDielectronVarManager::kTPCsub1rpH2);
+ // histos->UserHistogram("Event","","", 100, -2., 2., AliDielectronVarManager::kTPCsub2rpH2);
+ // histos->UserHistogram("Event","","", 100, -1., 1., AliDielectronVarManager::kTPCsub12DiffH2);
+
+ // further event plane resolutions (used in 3 sub event method)
+ histos->UserProfile("Event","","", AliDielectronVarManager::kv0ATPCDiffH2, 18, 0.,90., AliDielectronVarManager::kCentrality);
+ histos->UserProfile("Event","","", AliDielectronVarManager::kv0CTPCDiffH2, 18, 0.,90., AliDielectronVarManager::kCentrality);
+ histos->UserProfile("Event","","", AliDielectronVarManager::kv0Av0CDiffH2, 18, 0.,90., AliDielectronVarManager::kCentrality);
+ histos->UserProfile("Event","","", AliDielectronVarManager::kTPCsub12DiffH2, 18, 0.,90., AliDielectronVarManager::kCentrality);
+ histos->UserProfile("Event","","", AliDielectronVarManager::kv0Av0C0DiffH2, 18, 0.,90., AliDielectronVarManager::kCentrality);
+ histos->UserProfile("Event","","", AliDielectronVarManager::kv0Av0C3DiffH2, 18, 0.,90., AliDielectronVarManager::kCentrality);
+ histos->UserProfile("Event","","", AliDielectronVarManager::kv0Cv0A0DiffH2, 18, 0.,90., AliDielectronVarManager::kCentrality);
+ histos->UserProfile("Event","","", AliDielectronVarManager::kv0Cv0A3DiffH2, 18, 0.,90., AliDielectronVarManager::kCentrality);
+ histos->UserProfile("Event","","", AliDielectronVarManager::kv0A0v0A3DiffH2, 18, 0.,90., AliDielectronVarManager::kCentrality);
+ histos->UserProfile("Event","","", AliDielectronVarManager::kv0C0v0C3DiffH2, 18, 0.,90., AliDielectronVarManager::kCentrality);
+
+ // EP angle correlation (range of phi angle)
+ // histos->UserHistogram("Event","","", 320,-3.2.,3.2, 320,-3.2.,3.2,
+ // AliDielectronVarManager::kTPCrpH2,AliDielectronVarManager::kv0ACrpH2);
+
+ // EP Qvector magnitudes
+ histos->UserHistogram("Event","","", 200,0.,200., AliDielectronVarManager::kTPCmagH2);
+ histos->UserHistogram("Event","","", 200,0.,800., AliDielectronVarManager::kv0ACmagH2);
+ histos->UserHistogram("Event","","", 200,0.,800., AliDielectronVarManager::kv0AmagH2);
+ histos->UserHistogram("Event","","", 200,0.,800., AliDielectronVarManager::kv0CmagH2);
+
+ // detector effects checks
+ histos->UserHistogram("Event","","", 10,0.,100., 300,-1.0,1.0,
+ AliDielectronVarManager::kCentrality, AliDielectronVarManager::kTPCsub12DiffH2Sin);
+
+ // TPC recentering
+ // histos->UserProfile("Event","","", AliDielectronVarManager::kTPCxH2,
+ // AliDielectronHelper::MakeLinBinning(9, 0.,90.), GetRunNumbers2011(),
+ // AliDielectronVarManager::kCentrality, AliDielectronVarManager::kRunNumber);
+ // histos->UserProfile("Event","","", AliDielectronVarManager::kTPCyH2,
+ // AliDielectronHelper::MakeLinBinning(9, 0.,90.), GetRunNumbers2011(),
+ // AliDielectronVarManager::kCentrality, AliDielectronVarManager::kRunNumber);
+
+}
+
+void AddHistsEleEff(AliDielectron *die) {
+ //
+ // adding histograms for single electron efficiencies
+ //
+
+ Bool_t hasMC=die->GetHasMC();
+ AliDielectronHistos *histos = die->GetHistoManager();
+
+ // single electron efficiecy
+
+ // arbitrary binning for variables
+ // TVectorD *vpt = AliDielectronHelper::MakeLinBinning( 41, 0.0, 10.25);
+ // (*vpt)[41] = 20.;
+ TVectorD *vpt1 = AliDielectronHelper::MakeLinBinning( (int)(( 3. - 0.)/0.10), 0., 3.); //steps of 100MeV
+ TVectorD *vpt2 = AliDielectronHelper::MakeLinBinning( (int)(( 10.- 3.25)/0.25), 3.25, 10.); //steps of 250MeV
+ TVectorD *vpt3 = AliDielectronHelper::MakeLinBinning( (int)((100.-20.)/10.0), 20., 100.); //steps of 10GeV
+ TVectorD *vpt = new TVectorD(vpt1->GetNrows()+vpt2->GetNrows()+vpt3->GetNrows());
+ for(Int_t i=0; i<vpt1->GetNrows(); i++) (*vpt)[i] = (*vpt1)[i];
+ for(Int_t i=0; i<vpt2->GetNrows(); i++) (*vpt)[vpt1->GetNrows()+i] = (*vpt2)[i];
+ for(Int_t i=0; i<vpt3->GetNrows(); i++) (*vpt)[vpt1->GetNrows()+vpt2->GetNrows()+i] = (*vpt3)[i];
+ delete vpt1; delete vpt2; delete vpt3; //vpt->Print();
+
+ // array of bin limits
+ TObjArray *limEpm = new TObjArray();
+ // limEpm->Add(AliDielectronHelper::MakeLinBinning( 75,-10., 5.));
+ limEpm->Add(AliDielectronHelper::MakeLinBinning( 18, 0., 90.));
+ limEpm->Add(vpt);
+ limEpm->Add(AliDielectronHelper::MakeLinBinning( 20, 0., TMath::TwoPi()));
+ limEpm->Add(AliDielectronHelper::MakeLinBinning( 20, -1., +1.));
+ // limEpm->Add(AliDielectronHelper::MakeLinBinning( 1, -1., +1.));
+ // limEpm->Add(AliDielectronHelper::MakeLinBinning( 1, -1., +1.));
+ limEpm->Add(GetRunNumbers2011());
+ limEpm->Add(GetPDGcodes());
+ limEpm->Add(GetPDGcodes());
+ UInt_t varEpm[]={//AliDielectronVarManager::kTPCnSigmaEle,
+ AliDielectronVarManager::kCentrality
+ ,AliDielectronVarManager::kPt
+ ,AliDielectronVarManager::kPhi
+ ,AliDielectronVarManager::kEta
+ // ,AliDielectronVarManager::kImpactParXY
+ // ,AliDielectronVarManager::kImpactParZ
+ ,AliDielectronVarManager::kRunNumber
+ ,AliDielectronVarManager::kPdgCodeMother
+ ,AliDielectronVarManager::kPdgCodeGrandMother
+ };
+ // adding histogram
+ if(hasMC) histos->UserSparse("Track", limEpm->GetEntriesFast(), limEpm, varEpm);
+ delete limEpm;
+
+ histos->UserHistogram("Track","","", 500,-1.,1., AliDielectronVarManager::kImpactParXY);
+ histos->UserHistogram("Track","","", 600,-3.,3., AliDielectronVarManager::kImpactParZ);
+
+ // polarisation //
+ // array of bin limits
+ TVectorD *vpt = AliDielectronHelper::MakeLinBinning( 21, 0.0, 10.5);
+ (*vpt)[21] = 100.;
+ TObjArray *limPair = new TObjArray();
+ limPair->Add(AliDielectronHelper::MakeLinBinning( 18, 0., 90.));
+ limPair->Add(AliDielectronHelper::MakeLinBinning( 20, -1., +1.));
+ limPair->Add(vpt);
+ // limPair->Add(GetPDGcodes());
+ UInt_t varPair[]={
+ AliDielectronVarManager::kCentrality
+ ,AliDielectronVarManager::kThetaCS
+ ,AliDielectronVarManager::kPt
+ // ,AliDielectronVarManager::kPdgCodeMother
+ };
+ // if(hasMC) histos->UserSparse("Pair", limPair->GetEntriesFast(), limPair, varPair); // TAKES 4ever
+ delete limPair;
+
+}
+
+void AddQAHistsEff(AliDielectron *die) {
+ //
+ // adding histograms for single electron efficiencies
+ //
+
+ Bool_t hasMC=die->GetHasMC();
+ AliDielectronHistos *histos = die->GetHistoManager();
+
+ // arbitrary binning for variables
+ TVectorD *vpt1 = AliDielectronHelper::MakeLinBinning( (int)(( 3. - 0.)/0.10), 0., 3.);
+ TVectorD *vpt2 = AliDielectronHelper::MakeLinBinning( (int)(( 10.- 3.25)/0.25), 3.25, 10.);
+ TVectorD *vpt3 = AliDielectronHelper::MakeLinBinning( (int)((100.-20.)/10.0), 20., 100.);
+ TVectorD *vpt = new TVectorD(vpt1->GetNrows()+vpt2->GetNrows()+vpt3->GetNrows());
+ for(Int_t i=0; i<vpt1->GetNrows(); i++) (*vpt)[i] = (*vpt1)[i];
+ for(Int_t i=0; i<vpt2->GetNrows(); i++) (*vpt)[vpt1->GetNrows()+i] = (*vpt2)[i];
+ for(Int_t i=0; i<vpt3->GetNrows(); i++) (*vpt)[vpt1->GetNrows()+vpt2->GetNrows()+i] = (*vpt3)[i];
+ delete vpt1; delete vpt2; delete vpt3;
+
+ // single electron efficiecy
+ // applied efficiencies in collision data
+ histos->UserHistogram("Track","","", 101,-0.01,1.0, AliDielectronVarManager::kLegEff);
+ histos->UserProfile("Track","","", AliDielectronVarManager::kLegEff, GetRunNumbers2011(), AliDielectronVarManager::kRunNumber);
+ histos->UserProfile("Track","","", AliDielectronVarManager::kLegEff, 18,0.0,90.0, AliDielectronVarManager::kCentrality);
+ histos->UserProfile("Track","","", AliDielectronVarManager::kLegEff, vpt, AliDielectronVarManager::kPt);
+ histos->UserProfile("Track","","", AliDielectronVarManager::kLegEff, 20,0.0,TMath::TwoPi(), AliDielectronVarManager::kPhi);
+ histos->UserProfile("Track","","", AliDielectronVarManager::kLegEff, 20,-1.,+1., AliDielectronVarManager::kEta);
+ // histos->UserProfile("Track","","", AliDielectronVarManager::kLegEff,
+ // vpt, AliDielectronHelper::MakeLinBinning(20,0.0,TMath::TwoPi()),
+ // AliDielectronVarManager::kPt,AliDielectronVarManager::kPhi);
+
+ histos->UserHistogram("Pair","","", 101,-0.01,1.0, AliDielectronVarManager::kPairEff);
+ histos->UserProfile("Pair","","", AliDielectronVarManager::kPairEff, 125,.0,125*0.04, AliDielectronVarManager::kM);
+ histos->UserProfile("Pair","","", AliDielectronVarManager::kPairEff, 100,.0,10.0, AliDielectronVarManager::kPt);
+ histos->UserProfile("Pair","","", AliDielectronVarManager::kPairEff, 18,0.0,90.0, AliDielectronVarManager::kCentrality);
+
+ histos->UserHistogram("Pair","","", 125,.0,125*0.04, 101,-0.01,1.0,
+ AliDielectronVarManager::kM, AliDielectronVarManager::kPairEff);
+ histos->UserHistogram("Pair","","", 200,0,10., 101,-0.01,1.0,
+ AliDielectronVarManager::kPt, AliDielectronVarManager::kPairEff);
+
+ // array of bin limits
+ TObjArray *limEpm = new TObjArray();
+ limEpm->Add(AliDielectronHelper::MakeLinBinning( 125, 0., 125*0.04));
+ limEpm->Add(AliDielectronHelper::MakeLinBinning( 40, 0., 10.));
+ limEpm->Add(AliDielectronHelper::MakeLinBinning( 101, -0.01, +1.0));
+ UInt_t varEpm[]={
+ AliDielectronVarManager::kM
+ ,AliDielectronVarManager::kPt
+ ,AliDielectronVarManager::kPairEff
+ };
+ // histos->UserSparse("Pair", limEpm->GetEntriesFast(), limEpm, varEpm);
+ delete limEpm;
+
+ histos->UserProfile("Pair","","", AliDielectronVarManager::kPt,125,.0,125*0.04, AliDielectronVarManager::kM);
+ //weighted
+ histos->UserHistogram("Pair","","", AliDielectronHelper::MakeLinBinning(125,.0,125*0.04),
+ AliDielectronVarManager::kM, AliDielectronVarManager::kOneOverPairEff);
+ histos->UserProfile("Pair", "","", AliDielectronVarManager::kPt,AliDielectronHelper::MakeLinBinning(125,.0,125*0.04),
+ AliDielectronVarManager::kM, "", AliDielectronVarManager::kOneOverPairEff);
+}
+
+
+void InitHF(AliDielectron* die, Int_t cutDefinition)
+{
+ //
+ // Setup the HF arrays
+ //
+ // do not fill
+ switch(cutDefinition) {
+ case kNoCut:
+ case kEvent:
+ case kPIDqa:
+ case kLegEff:
+ case kNoPID:
+ case kRec:
+ case kSysMC:
+ return;
+ }
+
+ // has mc
+ Bool_t hasMC=die->GetHasMC();
+ // booleans for histo selection
+ Bool_t bHistEff = kFALSE; //ATTENTION
+ // mixing
+ AliDielectronMixingHandler *mixH=die->GetMixingHandler();
+ Int_t maxMixBins = (mixH ? mixH->GetNumberOfBins() : 0);
+
+ // container
+ AliDielectronHF *hf=new AliDielectronHF(die->GetName(),die->GetTitle());
+ // define pair types and sources
+ if(hasMC) hf->SetStepForMCGenerated();
+ hf->SetPairTypes(AliDielectronHF::kOSandMIX);
+ if(hasMC && cutDefinition!=kAvgPt) hf->SetPairTypes(AliDielectronHF::kMConly); // only mc truth
+ // hf->SetPairTypes(AliDielectronHF::kAll); // all pair types
+
+ //// define the grid size and granularity /////
+ hf->AddCutVariable(AliDielectronVarManager::kCentrality, AliDielectronHelper::MakeArbitraryBinning("0.,10.,40.,50.,90."));
+ hf->AddCutVariable(AliDielectronVarManager::kPt, AliDielectronHelper::MakeArbitraryBinning("0.,1.5,2.,2.5,3.,5.,6.,7.,8.,10.,100."));
+ hf->AddCutVariable(AliDielectronVarManager::kY, AliDielectronHelper::MakeArbitraryBinning("-0.9,-0.8,0.8,0.9"));
+
+ // hf->AddCutVariable(AliDielectronVarManager::kCentrality, 9, 0., 90.);
+ // TVectorD *vpt = AliDielectronHelper::MakeLinBinning( 21, 0.0, 10.5);
+ // (*vpt)[21] = 100.;
+ // hf->AddCutVariable(AliDielectronVarManager::kPt, vpt);
+ // hf->AddCutVariable(AliDielectronVarManager::kY, 1, -0.8, 0.8);
+ if(cutDefinition==kAvgPt ) //NEW
+ hf->AddCutVariable(AliDielectronVarManager::kPairEff, AliDielectronHelper::MakeArbitraryBinning(".0,.01,.02,.03,.04,.05,.1,1."));
+
+ // defaults
+ hf->UserHistogram("Pair", AliDielectronHelper::MakeLinBinning(125,.0,125*0.04), AliDielectronVarManager::kM);
+ // for mixed event weighting
+ // if(mixH) hf->UserHistogram("Pair", AliDielectronHelper::MakeLinBinning(maxMixBins,0,maxMixBins), AliDielectronVarManager::kMixingBin);
+
+ // mean pt analysis
+ if(cutDefinition==kAvgPt) {
+ hf->UserProfile("Pair", AliDielectronVarManager::kPt,
+ AliDielectronHelper::MakeLinBinning(125,.0,125*0.04), AliDielectronVarManager::kM);
+ // hf->UserProfile("Pair", AliDielectronVarManager::kPtSq,
+ // AliDielectronHelper::MakeLinBinning(125,.0,125*0.04), AliDielectronVarManager::kM);
+ // ME binning
+ // hf->UserHistogram("Pair", AliDielectronHelper::MakeLinBinning(125,.0,125*0.04),
+ // AliDielectronHelper::MakeLinBinning(maxMixBins,0,maxMixBins),
+ // AliDielectronVarManager::kM, AliDielectronVarManager::kMixingBin);
+ // hf->UserProfile("Pair", AliDielectronVarManager::kPt,
+ // AliDielectronHelper::MakeLinBinning(125,.0,125*0.04),
+ // AliDielectronHelper::MakeLinBinning(maxMixBins,0,maxMixBins),
+ // AliDielectronVarManager::kM, AliDielectronVarManager::kMixingBin);
+ }
+
+ // flow analysis
+ if(cutDefinition==kFlow && !hasMC) {
+ // flow versus minv
+ hf->UserProfile("Pair", AliDielectronVarManager::kv0ArpH2FlowV2,
+ AliDielectronHelper::MakeLinBinning(125,.0,125*0.04), AliDielectronVarManager::kM);
+ hf->UserProfile("Pair", AliDielectronVarManager::kv0CrpH2FlowV2,
+ AliDielectronHelper::MakeLinBinning(125,.0,125*0.04), AliDielectronVarManager::kM);
+ hf->UserProfile("Pair", AliDielectronVarManager::kTPCrpH2FlowV2,
+ AliDielectronHelper::MakeLinBinning(125,.0,125*0.04), AliDielectronVarManager::kM);
+ // detector effects
+ hf->UserProfile("Pair", AliDielectronVarManager::kCosTPCrpH2,
+ AliDielectronHelper::MakeLinBinning(1,.0,125*0.04), AliDielectronVarManager::kM);
+ hf->UserProfile("Pair", AliDielectronVarManager::kSinTPCrpH2,
+ AliDielectronHelper::MakeLinBinning(1,.0,125*0.04), AliDielectronVarManager::kM);
+ hf->UserProfile("Pair", AliDielectronVarManager::kCosPhiH2,
+ AliDielectronHelper::MakeLinBinning(1,.0,125*0.04), AliDielectronVarManager::kM);
+ hf->UserProfile("Pair", AliDielectronVarManager::kSinPhiH2,
+ AliDielectronHelper::MakeLinBinning(1,.0,125*0.04), AliDielectronVarManager::kM);
+ hf->UserProfile("Pair", AliDielectronVarManager::kTPCrpH2FlowV2Sin,
+ AliDielectronHelper::MakeLinBinning(125,.0,125*0.04), AliDielectronVarManager::kM);
+
+ // variables
+ hf->AddCutVariable(AliDielectronVarManager::kDeltaPhiv0ArpH2, GetDeltaPhiBins());
+ hf->AddCutVariable(AliDielectronVarManager::kDeltaPhiv0CrpH2, GetDeltaPhiBins());
+ hf->AddCutVariable(AliDielectronVarManager::kDeltaPhiTPCrpH2, GetDeltaPhiBins());
+ }
+
+ // on the fly efficienies
+ if(cutDefinition==kAvgPt/* && (!hasMC || 0)*/) {
+ hf->UserProfile("Pair", AliDielectronVarManager::kPairEff,
+ AliDielectronHelper::MakeLinBinning(125,.0,125*0.04), AliDielectronVarManager::kM);
+ hf->UserHistogram("Pair", AliDielectronHelper::MakeLinBinning(125,.0,125*0.04), AliDielectronVarManager::kM,
+ AliDielectronVarManager::kOneOverPairEff);
+ hf->UserProfile("Pair", AliDielectronVarManager::kPt,
+ AliDielectronHelper::MakeLinBinning(125,.0,125*0.04), AliDielectronVarManager::kM,
+ "", AliDielectronVarManager::kOneOverPairEff);
+ // hf->UserProfile("Pair", AliDielectronVarManager::kPtSq,
+ // AliDielectronHelper::MakeLinBinning(125,.0,125*0.04), AliDielectronVarManager::kM,
+ // "", AliDielectronVarManager::kOneOverPairEff);
+
+ // ME binning
+ if(maxMixBins) {
+ hf->UserHistogram("Pair", AliDielectronHelper::MakeLinBinning(125,.0,125*0.04),
+ AliDielectronHelper::MakeLinBinning(maxMixBins,0,maxMixBins),
+ AliDielectronVarManager::kM, AliDielectronVarManager::kMixingBin,
+ AliDielectronVarManager::kOneOverPairEff);
+ hf->UserProfile("Pair", AliDielectronVarManager::kPt,
+ AliDielectronHelper::MakeLinBinning(125,.0,125*0.04),
+ AliDielectronHelper::MakeLinBinning(maxMixBins,0,maxMixBins),
+ AliDielectronVarManager::kM, AliDielectronVarManager::kMixingBin,"",
+ AliDielectronVarManager::kOneOverPairEff);
+ }
+ }
+
+
+ /*
+//// define the grid size and granularity /////
+// event variables //
+// hf->AddCutVariable(AliDielectronVarManager::kCentrality, AliDielectronHelper::MakeArbitraryBinning("0,10,50,90")); // flow only
+if(hasMC) hf->AddCutVariable(AliDielectronVarManager::kCentrality, 18, 0., 90.);
+else hf->AddCutVariable(AliDielectronVarManager::kCentrality, 9, 0., 90.);
+// if(hasMC) hf->AddCutVariable(AliDielectronVarManager::kRunNumber, GetRunNumbers2011() );
+// if(hasMC) hf->AddCutVariable(AliDielectronVarManager::kNacc, 3000,0.,3000.);
+// if(hasMC) hf->AddCutVariable(AliDielectronVarManager::kNVtxContrib, 20,0.,4000.);
+
+ // pair variables //ATTENTION
+ if(hasMC && 0) hf->AddCutVariable(AliDielectronVarManager::kY, 18, -0.9, 0.9);
+ // hf->AddCutVariable(AliDielectronVarManager::kPt, 10, 0.0, 10.0);
+ if(!hasMC)
+ hf->AddCutVariable(AliDielectronVarManager::kPt, AliDielectronHelper::MakeArbitraryBinning("0,1,2,3,4,5,6,7,8,9,10,100"));
+ // hf->AddCutVariable(AliDielectronVarManager::kPt, 20, 0.0, 10.0);
+ else
+ hf->AddCutVariable(AliDielectronVarManager::kPt, 50, 0.0, 10.0);
+
+ // if(hasMC) hf->AddCutVariable(AliDielectronVarManager::kTRDpidEffPair,101,0.0,1.01);
+ // if(hasMC) hf->AddCutVariable(AliDielectronVarManager::kThetaCS,15,-1.,1.);
+
+ // flow variables //
+ // if(!hasMC) hf->AddCutVariable(AliDielectronVarManager::kDeltaPhiv0ArpH2, GetDeltaPhiBins());
+ // if(!hasMC) hf->AddCutVariable(AliDielectronVarManager::kDeltaPhiv0CrpH2, GetDeltaPhiBins());
+ if(!hasMC && bHistFlow) hf->AddCutVariable(AliDielectronVarManager::kDeltaPhiTPCrpH2, GetDeltaPhiBins());
+
+ // leg variables // NOTE: switched off in HF??
+ // if(hasMC) hf->AddCutVariable(AliDielectronVarManager::kPt, "0.85, 0.95, 1.1, 100.0", kTRUE, AliDielectronHF::kBinToMax);
+ // if(hasMC) hf->AddCutVariable(AliDielectronVarManager::kEta,"-0.9,-0.8,0.8,0.9", kTRUE, AliDielectronHF::kSymBin);
+ // hf->AddCutVariable(AliDielectronVarManager::kY, 1, -0.9, 0.9 );
+ // hf->AddCutVariable(AliDielectronVarManager::kPt, "0.8, 1.0, 1.1, 1.2, 1.5, 100.0", kTRUE, AliDielectronHF::kBinToMax);
+ // hf->AddCutVariable(AliDielectronVarManager::kNclsTPC, "70,90,100,120,160", kTRUE, AliDielectronHF::kBinToMax);
+ // hf->AddCutVariable(AliDielectronVarManager::kTPCnSigmaEle,"-4,-3,-2.5,-2,2,2.5,3,4", kTRUE, AliDielectronHF::kSymBin);
+ //hf->AddCutVariable(AliDielectronVarManager::kTPCnSigmaPio,"3.,3.5,4.,100.", kTRUE, AliDielectronHF::kBinToMax);
+ //hf->AddCutVariable(AliDielectronVarManager::kITSLayerFirstCls,4,0.,4., kFALSE, kTRUE, AliDielectronHF::kBinFromMin);
+ //hf->AddCutVariable(AliDielectronVarManager::kNclsITS, 5,2.,7., kFALSE, kTRUE, AliDielectronHF::kBinToMax);
+ //hf->AddCutVariable(AliDielectronVarManager::kRunNumber, GetRunNumbers2011());
+ */
+ die->SetHistogramArray(hf);
+}
void InitCF(AliDielectron* die, Int_t cutDefinition)
{
//
// Setup the CF Manager if needed
//
- // Bool_t hasMC=(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()!=0x0);
-
+ // do not fill
+ switch(cutDefinition) {
+ case kNoCut:
+ case kEvent:
+ case kPIDqa:
+ case kLegEff:
+ case kNoPID:
+ case kAvgPt: //////////////////NEW
+ return;
+ }
+
+ // has mc
+ Bool_t hasMC=die->GetHasMC();
+
+ // mixing
+ AliDielectronMixingHandler *mixH=die->GetMixingHandler();
+ Int_t maxMixBins = (mixH ? mixH->GetNumberOfBins() : 0);
+
+ // container
AliDielectronCF *cf=new AliDielectronCF(die->GetName(),die->GetTitle());
-
- // pair variables
+
+ //event variables
+ //cf->AddVariable(AliDielectronVarManager::kCentrality,"0.,5.0,10.0,15.0,20.0,25.0,30.0,35.0,40.0,50.0");
+ cf->AddVariable(AliDielectronVarManager::kCentrality,"0.,5.0,10.0,15.0,20.0,25.0,30.0,35.0,40.0,50.0,60.0,70.0,80.,90.");
+ // cf->AddVariable(AliDielectronVarManager::kRunNumber, GetRunNumbers2011() );
+ if(mixH && maxMixBins) cf->AddVariable(AliDielectronVarManager::kMixingBin, AliDielectronHelper::MakeLinBinning(maxMixBins,0,maxMixBins));
+
+ //pair variables
+ TVectorD *vpt = AliDielectronHelper::MakeLinBinning( 21, 0.0, 10.5);
+ (*vpt)[21] = 100.;
+ cf->AddVariable(AliDielectronVarManager::kPt, vpt);
+ // TVectorD *vpt = AliDielectronHelper::MakeLinBinning( 41, 0.0, 10.25);
+ // (*vpt)[41] = 100.;
+ // cf->AddVariable(AliDielectronVarManager::kPt,vpt);
+ // cf->AddVariable(AliDielectronVarManager::kPt,80,0.0,100*0.25);
+
+ // cf->AddVariable(AliDielectronVarManager::kY,"-5,-1,-0.9,-0.8,-0.5,0.5,0.8,0.9,1.0,5");
+ cf->AddVariable(AliDielectronVarManager::kY,"-0.9,-0.8,-0.7,+0.7,+0.8,+0.9");
+ // cf->AddVariable(AliDielectronVarManager::kY,18,-0.9,+0.9);
cf->AddVariable(AliDielectronVarManager::kM,125,0.,125*.04); //40Mev Steps
- if(cutDefinition!=kSubRndm) cf->AddVariable(AliDielectronVarManager::kPairType,11,0,11);
-
- if(cutDefinition < kGam0 || cutDefinition == kSubRndm) {
-
- // pair and event vars
- if(cutDefinition <= kChi || cutDefinition == kSubRndm) {
- cf->AddVariable(AliDielectronVarManager::kCentrality,"0.,5.,10.,20.,40.,50.,60.,80.");
- cf->AddVariable(AliDielectronVarManager::kPt,"0., 1., 2.5, 5., 100.0");
- if(!hasMC) cf->AddVariable(AliDielectronVarManager::kZvPrim,20, -10., 10.);
- if(hasMC) cf->AddVariable(AliDielectronVarManager::kNacc,20,0.,3000.0);
- if(hasMC) cf->AddVariable(AliDielectronVarManager::kNVtxContrib,20,0.,4000.);
- if(hasMC) cf->AddVariable(AliDielectronVarManager::kTRDpidEffPair,101,0.0,1.01);
- //cf->AddVariable(AliDielectronVarManager::kY,"-0.8,0.8");
- }
-
- //leg variables
- if(cutDefinition!=kSubRndm) {
- cf->AddVariable(AliDielectronVarManager::kPt,"0.8, 1.0, 1.1, 1.2, 1.5, 100.0",kTRUE);
- cf->AddVariable(AliDielectronVarManager::kTPCnSigmaEle,"-3,-2.5,-2,2,2.5,3",kTRUE);
- //cf->AddVariable(AliDielectronVarManager::kTPCnSigmaPio,"3.5,4.0,4.5,5.0,100",kTRUE);
- // cf->AddVariable(AliDielectronVarManager::kTPCnSigmaPro,"3.5,4.0,4.5,5.0,100",kTRUE);
- cf->AddVariable(AliDielectronVarManager::kTRDpidQuality,"3.5, 4.5, 5.5, 6.5",kTRUE);
-
- // standard vars
- if(cutDefinition<=kChi) {
- if(hasMC) cf->AddVariable(AliDielectronVarManager::kEta,"-0.9,0.9",kTRUE);
- //cf->AddVariable(AliDielectronVarManager::kNclsTPC,"70, 90, 100, 120, 160",kTRUE);
- }
- }
-
- switch(cutDefinition) {
- case kTPC:
- case krec:
- case kTOF: //cf->AddVariable(AliDielectronVarManager::kTOFnSigmaEle,"-3,-2,2,3",kTRUE); break;
- case kTRD:
- case kTOFTRD:
- // if(hasMC) cf->AddVariable(AliDielectronVarManager::kThetaCS,15,-1.,1.);
- //cf->AddVariable(AliDielectronVarManager::kITSLayerFirstCls,7,-1.5,5.5,kTRUE); break;
- break;
- case kITScls: cf->AddVariable(AliDielectronVarManager::kNclsITS,"1,2,3,4,5,6",kTRUE); break;
- case kITSamy: cf->AddVariable(AliDielectronVarManager::kITSLayerFirstCls,7,-1.5,5.5,kTRUE); break;
- case kDCA:
- cf->AddVariable(AliDielectronVarManager::kITSLayerFirstCls,7,-1.5,5.5,kTRUE);
- cf->AddVariable(AliDielectronVarManager::kImpactParXY,8,-2.,2.,kTRUE);
- cf->AddVariable(AliDielectronVarManager::kImpactParZ,8,-4.,4.,kTRUE);
- break;
- case kChi: cf->AddVariable(AliDielectronVarManager::kChi2NDF,"0,1,2,3,4,5",kTRUE); break;
- // cf->AddVariable(AliDielectronVarManager::kTOFnSigmaEle,"-3,-2,2,3",kTRUE);
- // cf->AddVariable(AliDielectronVarManager::kTOFPIDBit,"-.5,.5,1.5",kTRUE);
- }
-
+ // cf->AddVariable(AliDielectronVarManager::kPairType,1,1,1);
+ cf->AddVariable(AliDielectronVarManager::kPairType,11,0,11);
+ // cf->AddVariable(AliDielectronVarManager::kThetaCS,20,-1.,+1.);
+ // cf->AddVariable(AliDielectronVarManager::kThetaHE,20,-1.,+1.);
+ // cf->AddVariable(AliDielectronVarManager::kPhiCS,20,-3.2,+3.2);
+ // cf->AddVariable(AliDielectronVarManager::kPhiHE,20,-3.2,+3.2);
+ /// cf->AddVariable(AliDielectronVarManager::kDeltaPhiTPCrpH2,GetDeltaPhiBins());
+
+ if(cutDefinition==kAvgPt && 0) { //ATTENTION
+ cf->AddVariable(AliDielectronVarManager::kPairEff, AliDielectronHelper::MakeArbitraryBinning(".0,.01,.02,.03,.04,.05,.1,1."));
}
-
+ //leg variables
+ if(cutDefinition==kSysMC) { //ATTENTION
+ // if(cutDefinition!=kSysMCele && 0) { //ATTENTION
+ cf->AddVariable(AliDielectronVarManager::kPt,"0.85, 1.1, 100.0",kTRUE);
+ // cf->AddVariable(AliDielectronVarManager::kNclsTPC,"0, 70, 80, 90, 100, 120, 160",kTRUE);
+ cf->AddVariable(AliDielectronVarManager::kEta,"-0.9,-0.85,-0.8,-0.75,-0.70,0.70,0.75,0.8,0.85,0.9",kTRUE);
+ cf->AddVariable(AliDielectronVarManager::kITSLayerFirstCls,3,-1.5,1.5,kTRUE);
+ cf->AddVariable(AliDielectronVarManager::kTPCnSigmaPro,"3.5,4.0,4.5,100",kTRUE);
+ cf->AddVariable(AliDielectronVarManager::kTPCnSigmaEle,"-1.5,-1.0,-0.6,3.0",kTRUE);
+ }
+ else if(0){
+ cf->AddVariable(AliDielectronVarManager::kTPCnSigmaEle,"-3.0,-1.5,-1.0,-0.6,3.0",kTRUE);
+ }
+ /*
+ // event variables
+ cf->AddVariable(AliDielectronVarManager::kCentrality, 18,0., 90.);
+ if(hasMC) cf->AddVariable(AliDielectronVarManager::kRunNumber, GetRunNumbers2011() );
+ // if(hasMC) cf->AddVariable(AliDielectronVarManager::kNacc, 3000,0.,3000.);
+ // if(hasMC) cf->AddVariable(AliDielectronVarManager::kNVtxContrib, 20,0.,4000.);
+
+ // pair variables
+ if(hasMC) cf->AddVariable(AliDielectronVarManager::kPairType,1,1,1);
+ else cf->AddVariable(AliDielectronVarManager::kPairType,11,0,11);
+ cf->AddVariable(AliDielectronVarManager::kM, 125, 0.0, 5.0); //40Mev Steps
+ if(hasMC) cf->AddVariable(AliDielectronVarManager::kY, 18, -0.9, 0.9);
+ if(hasMC) cf->AddVariable(AliDielectronVarManager::kPt, 100, 0.0, 20.0);
+ else cf->AddVariable(AliDielectronVarManager::kPt, AliDielectronHelper::MakeArbitraryBinning("0,1,2,3,4,5,6,7,8,9,10,20"));
+ // if(hasMC) cf->AddVariable(AliDielectronVarManager::kTRDpidEffPair,101,0.0,1.01);
+ // if(hasMC) cf->AddVariable(AliDielectronVarManager::kThetaCS,15,-1.,1.);
+
+ // flow variables
+ //if(!hasMC) cf->AddVariable(AliDielectronVarManager::kDeltaPhiv0ArpH2, GetDeltaPhiBins());
+ //if(!hasMC) cf->AddVariable(AliDielectronVarManager::kDeltaPhiv0CrpH2, GetDeltaPhiBins());
+
+ // leg variables
+ if(hasMC) cf->AddVariable(AliDielectronVarManager::kPt,"0.0, 0.85, 0.95, 1.0, 1.1, 100.0",kTRUE);
+ if(hasMC) cf->AddVariable(AliDielectronVarManager::kEta,"-0.9,-0.8,-0.7, 0.7, 0.8, 0.9", kTRUE);
+ // cf->AddVariable(AliDielectronVarManager::kITSLayerFirstCls,7,-1.5,5.5,kTRUE);
+ // cf->AddVariable(AliDielectronVarManager::kNclsITS,"1,2,3,4,5,6",kTRUE);
+ // cf->AddVariable(AliDielectronVarManager::kTPCnSigmaEle,"-3,-2.5,-2,2,2.5,3",kTRUE);
+ // cf->AddVariable(AliDielectronVarManager::kTPCnSigmaPio,"2.5,3.0,3.5,4.0,4.5,100",kTRUE);
+ // cf->AddVariable(AliDielectronVarManager::kNclsTPC,"70, 90, 100, 120, 160",kTRUE);
+ // cf->AddVariable(AliDielectronVarManager::kTPCnSigmaPro,"3.5,4.0,4.5,5.0,100",kTRUE);
+ // cf->AddVariable(AliDielectronVarManager::kTOFnSigmaEle,"-3,-2,2,3",kTRUE); break;
+ // cf->AddVariable(AliDielectronVarManager::kTRDpidQuality,"3.5, 4.5, 5.5, 6.5",kTRUE);
+ // if(!hasMC && isESD) cf->AddVariable(AliDielectronVarManager::kTRDchi2,"-1.,0.,2.,4.",kTRUE);
+ */
+ // mc steps
if(hasMC) {
- cf->AddVariable(AliDielectronVarManager::kRunNumber, GetRunNumbers() );
- if(cutDefinition==kTOFTRD || cutDefinition==kGam0) cf->SetStepForMCtruth();
- // if(cutDefinition!=kTOFTRD)
- cf->SetStepsForMCtruthOnly();
- // cf->SetStepsForBackground();
+ cf->SetStepForMCtruth();
+ // cf->SetStepForNoCutsMCmotherPid();
+ // cf->SetStepForAfterAllCuts();
+ // cf->SetStepsForEachCut();
+ // cf->SetStepsForSignal();
+ // cf->SetStepsForBackground();
+ if(cutDefinition!=kAvgPt) cf->SetStepsForMCtruthOnly();
}
-
+ else
+ cf->SetStepsForSignal();
+
die->SetCFManagerPair(cf);
}
-void AddMCSignals(AliDielectron *die){
+void AddMCSignals(AliDielectron *die, Int_t cutDefinition){
//Do we have an MC handler?
- //Bool_t hasMC=(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()!=0x0);
- if (!hasMC) return;
-
- AliDielectronSignalMC* inclusiveJpsi = new AliDielectronSignalMC("inclusiveJpsi","Inclusive J/psi");
+ if (!die->GetHasMC()) return;
+
+ AliDielectronSignalMC* inclusiveJpsi = new AliDielectronSignalMC("inclusiveJpsi","Inclusive");
inclusiveJpsi->SetLegPDGs(11,-11);
inclusiveJpsi->SetMotherPDGs(443,443);
inclusiveJpsi->SetMothersRelation(AliDielectronSignalMC::kSame);
inclusiveJpsi->SetFillPureMCStep(kTRUE);
inclusiveJpsi->SetCheckBothChargesLegs(kTRUE,kTRUE);
inclusiveJpsi->SetCheckBothChargesMothers(kTRUE,kTRUE);
- die->AddSignalMC(inclusiveJpsi);
-
- AliDielectronSignalMC* promptJpsi = new AliDielectronSignalMC("promptJpsi","Prompt J/psi"); // prompt J/psi (not from beauty decays)
+
+ AliDielectronSignalMC* beautyJpsi = new AliDielectronSignalMC("beautyJpsi","Beauty");
+ beautyJpsi->SetLegPDGs(11,-11);
+ beautyJpsi->SetMotherPDGs(443,443);
+ beautyJpsi->SetMothersRelation(AliDielectronSignalMC::kSame);
+ beautyJpsi->SetGrandMotherPDGs(500,500);
+ beautyJpsi->SetFillPureMCStep(kTRUE);
+ beautyJpsi->SetCheckBothChargesLegs(kTRUE,kTRUE);
+ beautyJpsi->SetCheckBothChargesMothers(kTRUE,kTRUE);
+ beautyJpsi->SetCheckBothChargesGrandMothers(kTRUE,kTRUE);
+
+ AliDielectronSignalMC* promptJpsi = new AliDielectronSignalMC("promptJpsi","Prompt"); // prompt J/psi (not from beauty decays)
promptJpsi->SetLegPDGs(11,-11);
promptJpsi->SetMotherPDGs(443,443);
promptJpsi->SetGrandMotherPDGs(503,503,kTRUE,kTRUE); // not from beauty hadrons
promptJpsi->SetCheckBothChargesLegs(kTRUE,kTRUE);
promptJpsi->SetCheckBothChargesMothers(kTRUE,kTRUE);
promptJpsi->SetCheckBothChargesGrandMothers(kTRUE,kTRUE);
- die->AddSignalMC(promptJpsi);
-
- AliDielectronSignalMC* beautyJpsi = new AliDielectronSignalMC("beautyJpsi","Beauty J/psi");
- beautyJpsi->SetLegPDGs(11,-11);
- beautyJpsi->SetMotherPDGs(443,443);
- beautyJpsi->SetMothersRelation(AliDielectronSignalMC::kSame);
- beautyJpsi->SetGrandMotherPDGs(500,500);
- beautyJpsi->SetFillPureMCStep(kTRUE);
- beautyJpsi->SetCheckBothChargesLegs(kTRUE,kTRUE);
- beautyJpsi->SetCheckBothChargesMothers(kTRUE,kTRUE);
- beautyJpsi->SetCheckBothChargesGrandMothers(kTRUE,kTRUE);
- die->AddSignalMC(beautyJpsi);
-
- AliDielectronSignalMC* directJpsi = new AliDielectronSignalMC("directJpsi","Direct J/psi"); // embedded J/psi
+
+ // prompt J/psi radiative channel
+ AliDielectronSignalMC* promptJpsiRad = new AliDielectronSignalMC("promptJpsiRad","PromptRadiative"); // prompt J/psi (not from beauty decays)
+ promptJpsiRad->SetLegPDGs(11,-11);
+ promptJpsiRad->SetMotherPDGs(443,443);
+ promptJpsiRad->SetGrandMotherPDGs(503,503,kTRUE,kTRUE); // not from beauty hadrons
+ promptJpsiRad->SetMothersRelation(AliDielectronSignalMC::kSame);
+ promptJpsiRad->SetFillPureMCStep(kTRUE);
+ promptJpsiRad->SetLegSources(AliDielectronSignalMC::kFinalState, AliDielectronSignalMC::kFinalState);
+ promptJpsiRad->SetCheckBothChargesLegs(kTRUE,kTRUE);
+ promptJpsiRad->SetCheckBothChargesMothers(kTRUE,kTRUE);
+ promptJpsiRad->SetCheckBothChargesGrandMothers(kTRUE,kTRUE);
+ promptJpsiRad->SetJpsiRadiative(AliDielectronSignalMC::kIsRadiative);
+
+ // prompt J/psi Non radiative channel
+ AliDielectronSignalMC* promptJpsiNonRad = new AliDielectronSignalMC("promptJpsiNonRad","PromptNonRadiative"); // prompt J/psi (not from beauty decays)
+ promptJpsiNonRad->SetLegPDGs(11,-11);
+ promptJpsiNonRad->SetMotherPDGs(443,443);
+ promptJpsiNonRad->SetGrandMotherPDGs(503,503,kTRUE,kTRUE); // not from beauty hadrons
+ promptJpsiNonRad->SetMothersRelation(AliDielectronSignalMC::kSame);
+ promptJpsiNonRad->SetFillPureMCStep(kTRUE);
+ promptJpsiNonRad->SetLegSources(AliDielectronSignalMC::kFinalState, AliDielectronSignalMC::kFinalState);
+ promptJpsiNonRad->SetCheckBothChargesLegs(kTRUE,kTRUE);
+ promptJpsiNonRad->SetCheckBothChargesMothers(kTRUE,kTRUE);
+ promptJpsiNonRad->SetCheckBothChargesGrandMothers(kTRUE,kTRUE);
+ promptJpsiNonRad->SetJpsiRadiative(AliDielectronSignalMC::kIsNotRadiative);
+
+ AliDielectronSignalMC* directJpsi = new AliDielectronSignalMC("directJpsi","Direct"); // embedded J/psi
directJpsi->SetLegPDGs(11,-11);
directJpsi->SetMotherPDGs(443,443);
directJpsi->SetMothersRelation(AliDielectronSignalMC::kSame);
directJpsi->SetMotherSources(AliDielectronSignalMC::kDirect, AliDielectronSignalMC::kDirect);
directJpsi->SetCheckBothChargesLegs(kTRUE,kTRUE);
directJpsi->SetCheckBothChargesMothers(kTRUE,kTRUE);
- die->AddSignalMC(directJpsi);
-
- AliDielectronSignalMC* conversionElePairs = new AliDielectronSignalMC("conversionElePairs","conversion electron pairs"); // pairs made from conversion (may be also from 2 different conversions)
- conversionElePairs->SetLegPDGs(11,-11);
- conversionElePairs->SetCheckBothChargesLegs(kTRUE,kTRUE);
- conversionElePairs->SetLegSources(AliDielectronSignalMC::kSecondary, AliDielectronSignalMC::kSecondary);
- conversionElePairs->SetMotherPDGs(22,22);
- // die->AddSignalMC(conversionElePairs);
-}
-void SetEtaCorrection()
-{
- if (AliDielectronPID::GetEtaCorrFunction()) return;
-
- TString etaMap="$TRAIN_ROOT/jpsi_JPSI/EtaCorrMaps.root";
- TString trainRoot=gSystem->Getenv("TRAIN_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());
+ AliDielectronSignalMC* gammaConversion = new AliDielectronSignalMC("gammaConversion","gamma conversions");
+ gammaConversion->SetLegPDGs(11,-11);
+ gammaConversion->SetCheckBothChargesLegs(kTRUE,kTRUE);
+ gammaConversion->SetLegSources(AliDielectronSignalMC::kSecondary, AliDielectronSignalMC::kSecondary);
+ gammaConversion->SetMotherPDGs(22,22);
+ gammaConversion->SetMothersRelation(AliDielectronSignalMC::kSame);
+
+
+ AliDielectronSignalMC* electrons = new AliDielectronSignalMC("electrons","electrons");
+ electrons->SetLegPDGs(11,1); //dummy second leg (never MCtrue)
+ electrons->SetCheckBothChargesLegs(kTRUE,kTRUE);
+ electrons->SetLegSources(AliDielectronSignalMC::kFinalState, AliDielectronSignalMC::kFinalState);
+ // electrons->SetGrandMotherPDGs(500,500,kTRUE,kTRUE); // exclude non-prompt jpsi eletrons
+ electrons->SetFillPureMCStep(kTRUE);
+ // electrons->SetMothersRelation(AliDielectronSignalMC::kSame);
+
+ AliDielectronSignalMC* directElec = new AliDielectronSignalMC("directElec","directElec");
+ directElec->SetLegPDGs(11,1); //NEW
+ // directElec->SetMothersRelation(AliDielectronSignalMC::kSame);
+ // directElec->SetGrandMotherPDGs(-1103,-1103);
+ directElec->SetFillPureMCStep(kTRUE);
+ directElec->SetLegSources(AliDielectronSignalMC::kFinalState, AliDielectronSignalMC::kFinalState);
+ directElec->SetMotherSources(AliDielectronSignalMC::kDirect, AliDielectronSignalMC::kDirect);
+ directElec->SetCheckBothChargesLegs(kTRUE,kTRUE);
+ // new
+ directElec->SetCheckBothChargesGrandMothers(kTRUE,kTRUE);
+ directElec->SetGrandMotherPDGs(902,902,kTRUE,kTRUE); // exclude open charm,beauty hadrons
+
+
+ AliDielectronSignalMC* elecPrim = new AliDielectronSignalMC("elecPrim","elecPrim");
+ elecPrim->SetLegPDGs(11,1); //dummy second leg (never MCtrue)
+ elecPrim->SetCheckBothChargesLegs(kTRUE,kTRUE);
+ elecPrim->SetLegSources(AliDielectronSignalMC::kFinalState, AliDielectronSignalMC::kFinalState);
+ elecPrim->SetCheckBothChargesMothers(kTRUE,kTRUE);
+ elecPrim->SetMotherPDGs(902,902,kTRUE,kTRUE); // exclude open charm,beauty hadrons
+ elecPrim->SetCheckBothChargesGrandMothers(kTRUE,kTRUE);
+ elecPrim->SetGrandMotherPDGs(902,902,kTRUE,kTRUE); // exclude open charm,beauty hadrons
+ elecPrim->SetFillPureMCStep(kTRUE);
+
+ // add direct di lepton resonances
+ /*
+ AliDielectronSignalMC* directP[7];
+ TParticlePDG *ap;
+ Int_t pdg[] = {111, 113, 221, 223, 331, 333, 443};
+ for(Int_t i=0; i<7; i++) {
+ ap = TDatabasePDG::Instance()->GetParticle(pdg[i]);
+ directP[i] = new AliDielectronSignalMC(Form("direct%s",ap->GetName()),Form("direct%s",ap->GetName()));
+ directP[i]->SetLegPDGs(11,-11);
+ directP[i]->SetMotherPDGs(pdg[i],pdg[i]);
+ directP[i]->SetMothersRelation(AliDielectronSignalMC::kSame);
+ directP[i]->SetFillPureMCStep(kTRUE);
+ directP[i]->SetLegSources(AliDielectronSignalMC::kFinalState, AliDielectronSignalMC::kFinalState);
+ directP[i]->SetMotherSources(AliDielectronSignalMC::kDirect, AliDielectronSignalMC::kDirect);
+ // directP[i]->SetCheckBothChargesLegs(kTRUE,kTRUE);
+ // directP[i]->SetCheckBothChargesMothers(kTRUE,kTRUE);
+ }
+ */
+
+ /*
+ AliDielectronSignalMC* eleHijing = new AliDielectronSignalMC("eleHijing","eleHijing");
+ eleHijing->SetLegPDGs(11,1); //dummy second leg (never MCtrue)
+ eleHijing->SetCheckBothChargesLegs(kTRUE,kTRUE);
+ eleHijing->SetLegSources(AliDielectronSignalMC::kNoCocktail, AliDielectronSignalMC::kNoCocktail);
+ eleHijing->SetFillPureMCStep(kTRUE);
+ */
+ /*
+ AliDielectronSignalMC* electrons = new AliDielectronSignalMC("electrons","electrons");
+ electrons->SetLegPDGs(11,-11); //dummy second leg (never MCtrue)
+ electrons->SetCheckBothChargesLegs(kTRUE,kTRUE);
+ electrons->SetLegSources(AliDielectronSignalMC::kFinalState, AliDielectronSignalMC::kFinalState);
+ // electrons->SetMotherPDGs(111,111,kTRUE,kTRUE); // not from pi0
+ // electrons->SetMothersRelation(AliDielectronSignalMC::kSame);
+ electrons->SetFillPureMCStep(kTRUE);
+ */
+
+ // selection
+ switch(cutDefinition) {
+ case kPIDqa:
return;
+ //die->AddSignalMC(directElec);
+ break;
+ case kRec:
+ die->AddSignalMC(inclusiveJpsi);
+ die->AddSignalMC(directJpsi);
+ break;
+ case kStd:
+ die->AddSignalMC(inclusiveJpsi);
+ // die->AddSignalMC(beautyJpsi);
+ //die->AddSignalMC(promptJpsi);
+ //die->AddSignalMC(promptJpsiRad);
+ //die->AddSignalMC(promptJpsiNonRad);
+ die->AddSignalMC(directJpsi);
+ // die->AddSignalMC(gammaConversion);
+ break;
+ case kAvgPt:
+ die->AddSignalMC(inclusiveJpsi);
+ // die->AddSignalMC(directJpsi);
+ break;
+ case kLegEff:
+ // die->AddSignalMC(directJpsi);
+ // die->AddSignalMC(inclusiveJpsi);
+ // die->AddSignalMC(electrons);
+ // die->AddSignalMC(elecPrim);
+ die->AddSignalMC(directElec);
+ //for(Int_t i=0; i<7; i++) die->AddSignalMC(directP[i]);
+ break;
+ default:
+ die->AddSignalMC(directJpsi);
+ break;
}
- TFile f(etaMap.Data());
- if (!f.IsOpen()) return;
- TList *keys=f.GetListOfKeys();
- for (Int_t i=0; i<keys->GetEntries(); ++i){
- TString kName=keys->At(i)->GetName();
- TPRegexp reg(kName);
- if (reg.MatchB(list)){
- printf("Using Eta Correction Function: %s\n",kName.Data());
- AliDielectronPID::SetEtaCorrFunction((TF1*)f.Get(kName.Data()));
- }
+}
+
+void SetEtaCorrection(AliDielectron *die) {
+
+
+ // if(cutDefinition==kLegEff) return;
+ // if (pid->GetCentroidCorrFunction()) return;
+ Bool_t hasMC=die->GetHasMC();
+ Bool_t hasTuneOnData=kFALSE;
+ // ((AliAnalysisTaskPIDResponse*)AliAnalysisManager::GetAnalysisManager()->GetTasks()->At(0))->GetTuneOnData();
+ printf("tune on data switched: %d \n",hasTuneOnData);
+ // printf("name task at 0: %s \n",AliAnalysisManager::GetAnalysisManager()->GetTasks()->At(0)->GetName());
+ // printf("input event %p \n", AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());
+ // printf("pid response %p \n",((AliInputEventHandler*)AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler())->GetPIDResponse());
+ // printf("pid response task %p \n",AliAnalysisManager::GetAnalysisManager()->GetTasks()->At(0));
+ // AliAnalysisManager::GetAnalysisManager()->GetTasks()->At(0)->Dump();;
+
+ // AliAnalysisManager* man = AliAnalysisManager::GetAnalysisManager();
+ // AliInputEventHandler* inputHandler = dynamic_cast<AliInputEventHandler*>(man->GetInputEventHandler());
+ // AliPIDResponse* pidResponse = inputHandler->GetPIDResponse();
+ // if(pidResponse) hasTuneOnData = pidResponse->IsTunedOnData();
+ // printf("man %p inp %p pid %p ====> %d \n",man,inputHandler,pidResponse,hasTuneOnData);
+
+ TF2 *fCntrdCorr=0x0;
+ TF1 *fWdthCorr=0x0;
+ /* vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv DATA vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv */
+ // either data or MC with tune on data option
+ if( !hasMC /*|| hasTuneOnData*/ ) {
+ // 2-dimensional eta correction for the centroid of electron sigmas
+ fCntrdCorr = new TF2("fCntrdCorr",
+ "[0] + [1]*y + [2]*y*y + [3]*TMath::Power(y,3) + [4]*TMath::Power(y,4) + [5]*TMath::Power(y,5) + [6]*TMath::Power(y,6) + [7]*x",
+ // 0.0, 3000.0, -0.9, +0.9);
+ 0.0, 90.0, -0.9, +0.9);
+ // fCntrdCorr->SetParameters(0.723106, 0.23958, -6.31221, -0.687976, 15.912, 0.579609, -11.6901, -0.000354381); //Nacc dep.
+ fCntrdCorr->SetParameters(+0.149002, +0.214644 , -6.034930, -0.529588, +14.97902, +0.402640, -10.890027, +0.011248); //Cent
+
+ // 1-dimensional eta correction for the width of electron sigmas
+ fWdthCorr = new TF1("fWdthCorr", "pol2", 0.0, 90.0);
+ // fWdthCorr->SetParameters(1.06108, 0.000217804,-5.80291e-08); //Nacc dep.
+ fWdthCorr->SetParameters(+1.290755, -0.005261, +0.000021); //Cent dep.
+
+ // apply corrections
+ // die->SetCentroidCorrFunction(fCntrdCorr,AliDielectronVarManager::kNacc,AliDielectronVarManager::kEta);
+ // die->SetWidthCorrFunction(fWdthCorr,AliDielectronVarManager::kNacc);
+ die->SetCentroidCorrFunction(fCntrdCorr,AliDielectronVarManager::kCentrality,AliDielectronVarManager::kEta);
+ die->SetWidthCorrFunction(fWdthCorr,AliDielectronVarManager::kCentrality);
+ printf(" DATA PID correction loaded!!!\n");
}
+ else {
+ /* vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv MONTE CARLO vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv */
+ // 2-dimensional eta correction for the centroid and width of electron sigmas
+ fCntrdCorr = new TF2("fCntrdCorr",
+ "[0] + [1]*y + [2]*y*y + [3]*TMath::Power(y,3) + [4]*TMath::Power(y,4) + [5]*TMath::Power(y,5) + [6]*TMath::Power(y,6) + [7]*x",
+ 0.0, 3000.0, -0.9, +0.9);
+ fCntrdCorr->SetParameters(+0.293718,
+ +0.010037,
+ -2.632949,
+ -0.241412,
+ +8.304244,
+ +0.525481,
+ -4.874357,
+ -0.000103); //TPCrefMult dep.
+ fWdthCorr = new TF2("fWdthCorr",
+ "[0] + [1]*y + [2]*y*y + [3]*TMath::Power(y,3) + [4]*TMath::Power(y,4) + [5]*TMath::Power(y,5) + [6]*TMath::Power(y,6) + [7]*x",
+ 0.0, 3000.0, -0.9, +0.9);
+ fWdthCorr->SetParameters(+0.917840,
+ -0.021500,
+ -0.628371,
+ +0.230847,
+ +1.434907,
+ -0.330751,
+ -0.458941,
+ +0.000036); //TPCrefMult dep.
+
+ // apply corrections
+ die->SetCentroidCorrFunction(fCntrdCorr,
+ AliDielectronVarManager::kRefMultTPConly,
+ AliDielectronVarManager::kEta);
+ die->SetWidthCorrFunction(fWdthCorr,
+ AliDielectronVarManager::kRefMultTPConly,
+ AliDielectronVarManager::kEta);
+ /*
+ // 2-dimensional eta correction for the centroid of electron sigmas
+ fCntrdCorr = new TF2("fCntrdCorr", "[0] + [1]*y + [2]*y*y + [3]*TMath::Power(y,3) + [4]*TMath::Power(y,4) + [5]*TMath::Power(y,5) + [6]*TMath::Power(y,6) + [7]*x",
+ 0.0, 3000.0, -0.9, +0.9);
+ fCntrdCorr->SetParameters(+0.378611, -0.070831, -3.076778, +0.121977, +8.576097, +0.113009, -5.001368, -0.000181);
+ // 1-dimensional eta correction for the width of electron sigmas
+ fWdthCorr = new TF1("fWdthCorr", "pol1", 0.0, 3000.0);
+ fWdthCorr->SetParameters(+0.881894, +0.000053);
+
+ // apply corrections
+ die->SetCentroidCorrFunction(fCntrdCorr,AliDielectronVarManager::kNacc,AliDielectronVarManager::kEta);
+ die->SetWidthCorrFunction(fWdthCorr,AliDielectronVarManager::kNacc);
+ // die->SetCentroidCorrFunction(fCntrdCorr,AliDielectronVarManager::kCentrality,AliDielectronVarManager::kEta);
+ // die->SetWidthCorrFunction(fWdthCorr,AliDielectronVarManager::kCentrality);
+ */
+ printf(" MC PID correction loaded!!!\n");
+ }
+
+
}
-TVectorD *GetRunNumbers() {
+TVectorD *GetRunNumbers2011() {
Double_t runLHC10h[] = { // all good runs based on RCT 29.Mai
139510, 139507, 139505, 139503, 139465, 139438, 139437, 139360, 139329, 139328, 139314, 139310, 139309, 139173, 139107, 139105, 139038, 139037, 139036, 139029, 139028, 138872, 138871, 138870, 138837, 138732, 138730, 138666, 138662, 138653, 138652, 138638, 138624, 138621, 138583, 138582, 138579, 138578, 138534, 138469, 138442, 138439, 138438, 138396, 138364, 138275, 138225, 138201, 138197, 138192, 138190, 137848, 137844, 137752, 137751, 137724, 137722, 137718, 137704, 137693, 137692, 137691, 137686, 137685, 137639, 137638, 137608, 137595, 137549, 137546, 137544, 137541, 137539, 137531, 137530, 137443, 137441, 137440, 137439, 137434, 137432, 137431, 137430, 137366, 137243, 137236, 137235, 137232, 137231, 137230, 137162, 137161, 137135
};
-
Double_t runLHC11h[] = { // all good runs based on RCT 29.Mai
- 170593, 170572, 170388, 170387, 170315, 170313, 170312, 170311, 170309, 170308, 170306, 170270, 170269, 170268, 170230, 170228, 170207, 170204, 170203, 170193, 170163, 170159, 170155, 170091, 170089, 170088, 170085, 170084, 170083, 170081, 170040, 170027, 169965, 169923, 169859, 169858, 169855, 169846, 169838, 169837, 169835, 169591, 169588, 169586, 169557, 169554, 169550, 169512, 169504, 169498, 169475, 169419, 169418, 169417, 169415, 169411, 169238, 169167, 169160, 169156, 169148, 169145, 169144, 169138, 169099, 169094, 169091, 169044, 169035, 168992, 168988, 168826, 168777, 168514, 168512, 168511, 168467, 168464, 168460, 168458, 168362, 168361, 168342, 168341, 168325, 168322, 168311, 168310, 168115, 168108, 168107, 168105, 168076, 168069, 167988, 167987, 167985, 167920, 167915
+ 167915, 167920, 167985, 167987, 167988, 168069, 168076, 168105, 168107, 168108, 168115, 168310, 168311, 168322, 168325, 168341, 168342, 168361, 168362, 168458, 168460, 168464, 168467, 168511, 168512, 168514, 168777, 168826, 168988, 168992, 169035, 169040, 169044, 169045, 169091, 169094, 169099, 169138, 169144, 169145, 169148, 169156, 169160, 169167, 169238, 169411, 169415, 169417, 169418, 169419, 169420, 169475, 169498, 169504, 169506, 169512, 169515, 169550, 169553, 169554, 169555, 169557, 169586, 169587, 169588, 169590, 169591, 169835, 169837, 169838, 169846, 169855, 169858, 169859, 169923, 169965, 170027, 170040, 170081, 170083, 170084, 170085, 170088, 170089, 170091, 170155, 170159, 170163, 170193, 170203, 170204, 170207, 170228, 170230, 170268, 170269, 170270, 170306, 170308, 170309, 170311, 170312, 170313, 170315, 170387, 170388, 170572, 170593,
+ 170593+1
};
-
- // selection via environement variable (works only for gsi trains)
-
-
+ /*
if(list.Contains("LHC10h") || list.Contains("LHC11a10")) {
- Int_t size = (int) (sizeof(runLHC10h)/sizeof(Double_t));
- TVectorD *vec = new TVectorD(size+1);
-
- (*vec)[size] = runLHC10h[0] + 1;
- for (int i = 0; i < size; i++) {
- (*vec)[i] = runLHC10h[size-1-i];
- }
- // vec->Print("");
- return vec;
- }
-
if( list.Contains("LHC11h") || list.Contains("LHC12a17") ) {
-
- Int_t size = (int) (sizeof(runLHC11h)/sizeof(Double_t));
- TVectorD *vec = new TVectorD(size+1);
-
- (*vec)[size] = runLHC11h[0] + 1;
- for (int i = 0; i < size; i++) {
- (*vec)[i] = runLHC11h[size-1-i];
+ */
+ Int_t size = (int) (sizeof(runLHC11h)/sizeof(Double_t));
+ TVectorD *vec = new TVectorD(size,runLHC11h);
+ //vec->Print("");
+ return vec;
+}
+
+TVectorD *GetPDGcodes() {
+ //
+ // array of pdgcodes stored in TDatabasePDG
+ //
+ TDatabasePDG *pdg = new TDatabasePDG();
+ pdg->ReadPDGTable();
+ TGraph *gr = new TGraph();
+ TIter next(pdg->ParticleList());
+ TParticlePDG *p;
+ Int_t i=0;
+ while ((p = (TParticlePDG *)next())) {
+ if(TMath::Abs(p->PdgCode()) < 1e+6) {
+ // printf("%s -> %d \n",p->GetName(),p->PdgCode());
+ gr->SetPoint(i++, p->PdgCode(),1.);
}
- // vec->Print("");
- return vec;
}
-
+ gr->Sort();
+ TVectorD *vec = new TVectorD(gr->GetN(), gr->GetX());
+ // vec->Print();
+ delete pdg;
+ delete gr;
+ return vec;
+
+}
+
+TVectorD *GetDeltaPhiBins() {
+ //
+ // for in and out of event plane bins
+ //
+ Double_t pi = TMath::Pi();
+ TVectorD *deltaPhi = new TVectorD(6);
+ (*deltaPhi)[0] = -1. *pi;
+ (*deltaPhi)[1] = -3./4. *pi;
+ (*deltaPhi)[2] = -1./4. *pi;
+ (*deltaPhi)[3] = +1./4. *pi;
+ (*deltaPhi)[4] = +3./4. *pi;
+ (*deltaPhi)[5] = +1. *pi;
+ return deltaPhi;
}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff