100,-3.15,3.15,AliDielectronVarManager::kDeltaPhiv0ACrpH2);
- histos->UserHistogram("Pair","PairPlanev0ACrpH2Angle","Phi;Phi[rad];#counts",
- 100,0,1.6,AliDielectronVarManager::kPairPlanev0rpH2Angle);
- histos->UserHistogram("Pair","PairPlaneMagAngle","Phi;Phi[rad];#counts",
- 100,0,1.6,AliDielectronVarManager::kPairPlaneMagAngle);
+ // histos->UserHistogram("Pair","PairPlanev0ACrpH2Angle","Phi;Phi[rad];#counts",
+ // 100,0,1.6,AliDielectronVarManager::kPairPlanev0rpH2Angle);
+ // histos->UserHistogram("Pair","PairPlaneMagAngle","Phi;Phi[rad];#counts",
+ // 100,0,1.6,AliDielectronVarManager::kPairPlaneMagAngle);
histos->UserHistogram("Pair","PairPlaneAngle","Phi;Phi[rad];#counts",
100,0,1.6,AliDielectronVarManager::kPairPlaneAngle);
histos->UserHistogram("Pair","InvMassHighOpeningAngle","Opening Angle vs Inv.Mass;Inv. Mass [GeV];#pairs",
200,0.3,0.5,200,0.,6.3,AliDielectronVarManager::kM,AliDielectronVarManager::kOpeningAngle);
- */
+
histos->UserHistogram("Pair","InvMassAllPairplaneMagAngle","ee plane and Mag Vector Angle vs Inv.Mass;Inv. Mass [GeV];Phi [rad]",
180,0.12,0.3,100,0.,3.15,AliDielectronVarManager::kM,AliDielectronVarManager::kPairPlanev0rpH2Angle);
histos->UserHistogram("Pair","InvMassHighPairplaneRPAngle","ee plane and Mag Vector Angle vs Inv.Mass;Inv. Mass [GeV];Phi [rad]",
200,0.3,0.5,100,0.,3.15,AliDielectronVarManager::kM,AliDielectronVarManager::kPairPlanev0rpH2Angle);
-
+ */
histos->UserHistogram("Pair","InvMassAllPairplaneAngle","ee plane and Mag Vector Angle vs Inv.Mass;Inv. Mass [GeV];Phi [rad]",
1000,0.0,5.0,100,0.,3.15,AliDielectronVarManager::kM,AliDielectronVarManager::kPairPlaneAngle);
histos->UserHistogram("Pair","InvMassAll2PairplaneAngle","ee plane and Mag Vector Angle vs Inv.Mass;Inv. Mass [GeV];Phi [rad]",
500,0.0,0.50,100,0.,3.15,AliDielectronVarManager::kM,AliDielectronVarManager::kPairPlaneAngle);
+ histos->UserHistogram("Pair","InvMassMiddle2PairplaneAngle","ee plane and Mag Vector Angle vs Inv.Mass;Inv. Mass [GeV];Phi [rad]",
+ 400,0.1,0.50,100,0.,3.15,AliDielectronVarManager::kM,AliDielectronVarManager::kPairPlaneAngle);
+
histos->UserHistogram("Pair","InvMassLoweeplaneAngle","ee plane and Mag Vector Angle vs Inv.Mass;Inv. Mass [GeV];Phi [rad]",
300,0.0,0.03,100,0.,3.15,AliDielectronVarManager::kM,AliDielectronVarManager::kPairPlaneAngle);
histos->UserHistogram("Pair","InvMassMiddleeelaneAngle","ee plane and Mag Vector Angle vs Inv.Mass;Inv. Mass [GeV];Phi [rad]",
180,0.12,0.3,100,0.,3.15,AliDielectronVarManager::kM,AliDielectronVarManager::kPairPlaneAngle);
histos->UserHistogram("Pair","InvMassHighPairplaneAngle","ee plane and Mag Vector Angle vs Inv.Mass;Inv. Mass [GeV];Phi [rad]",
200,0.3,0.5,100,0.,3.15,AliDielectronVarManager::kM,AliDielectronVarManager::kPairPlaneAngle);
-
+
//add histograms to Track classes
}
-
AliESDtrackCuts *SetupESDtrackCutsDieleData(Int_t cutDefinition);
-TString namesDieleData=("baseMixzVert;baseMixzVert+OnlyComplete;baseMixzVert40Bins");
+TString namesDieleData=("baseMixzVert40Bins;baseMixzVert40BinsComplete;kMixLegs");
TObjArray *arrNamesDieleData=namesDieleData.Tokenize(";");
// mixing
AliDielectronMixingHandler *mix=new AliDielectronMixingHandler;
- mix->AddVariable(AliDielectronVarManager::kZvPrim,20,-10.,10.);
- if(cutDefinition == 2)mix->AddVariable(AliDielectronVarManager::kZvPrim,40,-10.,10.);
- mix->SetDepth(100);
+ mix->AddVariable(AliDielectronVarManager::kZvPrim,40,-10.,10.);
if(cutDefinition == 1)mix->SetMixUncomplete(kFALSE);
+ mix->SetDepth(100);
mix->SetMixType(AliDielectronMixingHandler::kAll);
diele->SetMixingHandler(mix);
-// check MultV0
-
// old code
+// if(cutDefinition == 2)mix->AddVariable(AliDielectronVarManager::kZvPrim,20,-10.,10.);
// if(cutDefinition==1)mix->AddVariable(AliDielectronVarManager::kRefMult,"1,35,44,50,57,62,68,73,79,85,91,97,104,111,119,129,140,155,176,212,600");// 20 binssame integral per bin
// if(cutDefinition==2)mix->AddVariable(AliDielectronVarManager::kVZEROchMult,"1,2,3,4,5,6,7,8,9,10,12,120");
// if(cutDefinition==3){mix->SetDepth(20); mix->SetMixUncomplete(kFALSE);}
{
// Setup the pair cuts
AliDielectronVarCuts *mycut = new AliDielectronVarCuts("CutEMCAL","cut for EMCal");
-// mycut->AddCut(AliDielectronVarManager::kEMCALnSigmaEle,-2.,3.);//not needed anymore
+ mycut->AddCut(AliDielectronVarManager::kEMCALnSigmaEle,-2.,3.);//not needed anymore
mycut->AddCut(AliDielectronVarManager::kEMCALE,3.5,100.);
- mycut->AddCut(AliDielectronVarManager::kEMCALEoverP,0.75,1.25);//
+// mycut->AddCut(AliDielectronVarManager::kEMCALEoverP,0.75,1.25);//
AliDielectronPairLegCuts *varpair=new AliDielectronPairLegCuts();
varpair->GetLeg1Filter().AddCuts(mycut);
varpair->GetLeg2Filter().AddCuts(mycut);
- varpair->SetCutType(AliDielectronPairLegCuts::kAnyLeg);
+ if(cutDefinition == 0 || cutDefinition == 1){
+ varpair->SetCutType(AliDielectronPairLegCuts::kAnyLeg);
+ } else if(cutDefinition == 2){
+ varpair->SetCutType(AliDielectronPairLegCuts::kMixLegs);
+ }
diele->GetPairFilter().AddCuts(varpair);
histos->UserHistogram("Track","dEdx_Eta","dEdx vs eta;#eta;TPC signal (arb units);#tracks",
200,-1.,1.,800,20.,200.,AliDielectronVarManager::kEta,AliDielectronVarManager::kTPCsignal,kTRUE);
- /*
- histos->UserHistogram("Track","kNclsSTPC_kDeltaEta","kNclsSTPC vs kDeltaEta;kNclsSTPC;kDeltaEta;#tracks",
- 100,0.,1.,160,0.,1.6,AliDielectronVarManager::kNclsSTPC,AliDielectronVarManager::kDeltaEta,kTRUE);
- histos->UserHistogram("Track","kNclsSTPC_kDeltaPhi","kNclsSTPC vs kDeltaPhi;kNclsSTPC;kDeltaPhi;#tracks",
- 100,0.,1.,160,0.,6.4,AliDielectronVarManager::kNclsSTPC,AliDielectronVarManager::kDeltaPhi,kTRUE);
- histos->UserHistogram("Track","kDeltaPhi_kDeltaEta","kDeltaPhi vs kDeltaEta;kDeltaPhi;kDeltaEta;#tracks",
- 160,0.,6.4,160,0.,1.6,AliDielectronVarManager::kDeltaPhi,AliDielectronVarManager::kDeltaEta,kTRUE);
- */
histos->UserHistogram("Track","dEdx_nSigmaEMCal","dEdx vs nSigmaEMCal;NsigmaEmcal;TPC signal (arb units);NSigmaEMCAL",
200,-5.,5.,800,20.,200.,AliDielectronVarManager::kEMCALnSigmaEle,AliDielectronVarManager::kTPCsignal,kTRUE);
cf->AddVariable(AliDielectronVarManager::kPairType,12,0,12);
cf->AddVariable(AliDielectronVarManager::kEta,20,-1.,1.);
cf->AddVariable(AliDielectronVarManager::kY,20,-1.,1.);
- cf->AddVariable(AliDielectronVarManager::kPhi,20,0.,20*0.32);
-
+ cf->AddVariable(AliDielectronVarManager::kPhi,64,-64.*0.1,64.*0.1);
+ cf->AddVariable(AliDielectronVarManager::kOpeningAngle,64,0.,3.2);
+
//global leg variables
cf->AddVariable(AliDielectronVarManager::kPt,"0.,0.5,0.75,0.9,1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 2.0, 3.0, 4.0, 8.0,9.0,10.0,11.0,12.0,13.0,14.0,15.0,16.0,17.0,18.0,19.0,20.0,30.0,50.0,100.0",kTRUE);
cf->AddVariable(AliDielectronVarManager::kEta,44,-1.2,1.2,kTRUE);