\r
void AliAnalysisTaskEventMixingBF::UserExecMix(Option_t *)\r
{\r
+ // Main loop for event mixing\r
\r
TString gAnalysisLevel = fBalance->GetAnalysisLevel();\r
\r
chargeVector[i] = new vector<Double_t>;\r
}\r
\r
- Double_t v_charge;\r
- Double_t v_y;\r
- Double_t v_eta;\r
- Double_t v_phi;\r
- Double_t v_p[3];\r
- Double_t v_pt;\r
- Double_t v_E;\r
+ Double_t vCharge;\r
+ Double_t vY;\r
+ Double_t vEta;\r
+ Double_t vPhi;\r
+ Double_t vP[3];\r
+ Double_t vPt;\r
+ Double_t vE;\r
\r
Int_t iMainTrackUsed = -1;\r
\r
fHistTrackStats->Fill(aodTrackMain->GetFilterMap());\r
if(!aodTrackMain->TestFilterBit(nAODtrackCutBit)) continue;\r
\r
- v_charge = aodTrackMain->Charge();\r
- v_y = aodTrackMain->Y();\r
- v_eta = aodTrackMain->Eta();\r
- v_phi = aodTrackMain->Phi() * TMath::RadToDeg();\r
- v_E = aodTrackMain->E();\r
- v_pt = aodTrackMain->Pt();\r
- aodTrackMain->PxPyPz(v_p);\r
+ vCharge = aodTrackMain->Charge();\r
+ vY = aodTrackMain->Y();\r
+ vEta = aodTrackMain->Eta();\r
+ vPhi = aodTrackMain->Phi() * TMath::RadToDeg();\r
+ vE = aodTrackMain->E();\r
+ vPt = aodTrackMain->Pt();\r
+ aodTrackMain->PxPyPz(vP);\r
\r
- Float_t DCAxyMain = aodTrackMain->DCA(); // this is the DCA from global track (not exactly what is cut on)\r
- Float_t DCAzMain = aodTrackMain->ZAtDCA(); // this is the DCA from global track (not exactly what is cut on)\r
+ Float_t dcaXYMain = aodTrackMain->DCA(); // this is the DCA from global track (not exactly what is cut on)\r
+ Float_t dcaZMain = aodTrackMain->ZAtDCA(); // this is the DCA from global track (not exactly what is cut on)\r
\r
\r
// Kinematics cuts from ESD track cuts\r
- if( v_pt < fPtMin || v_pt > fPtMax) continue;\r
- if( v_eta < fEtaMin || v_eta > fEtaMax) continue;\r
+ if( vPt < fPtMin || vPt > fPtMax) continue;\r
+ if( vEta < fEtaMin || vEta > fEtaMax) continue;\r
\r
// Extra DCA cuts (for systematic studies [!= -1])\r
if( fDCAxyCut != -1 && fDCAzCut != -1){\r
- if(TMath::Sqrt((DCAxyMain*DCAxyMain)/(fDCAxyCut*fDCAxyCut)+(DCAzMain*DCAzMain)/(fDCAzCut*fDCAzCut)) > 1 ){\r
+ if(TMath::Sqrt((dcaXYMain*dcaXYMain)/(fDCAxyCut*fDCAxyCut)+(dcaZMain*dcaZMain)/(fDCAzCut*fDCAzCut)) > 1 ){\r
continue; // 2D cut\r
}\r
}\r
\r
// fill QA histograms\r
fHistClus->Fill(aodTrackMain->GetITSNcls(),aodTrackMain->GetTPCNcls());\r
- fHistDCA->Fill(DCAzMain,DCAxyMain);\r
+ fHistDCA->Fill(dcaZMain,dcaXYMain);\r
fHistChi2->Fill(aodTrackMain->Chi2perNDF());\r
- fHistPt->Fill(v_pt);\r
- fHistEta->Fill(v_eta);\r
- fHistPhi->Fill(v_phi);\r
+ fHistPt->Fill(vPt);\r
+ fHistEta->Fill(vEta);\r
+ fHistPhi->Fill(vPhi);\r
\r
// fill charge vector\r
- chargeVector[0]->push_back(v_charge);\r
- chargeVector[1]->push_back(v_y);\r
- chargeVector[2]->push_back(v_eta);\r
- chargeVector[3]->push_back(v_phi);\r
- chargeVector[4]->push_back(v_p[0]);\r
- chargeVector[5]->push_back(v_p[1]);\r
- chargeVector[6]->push_back(v_p[2]);\r
- chargeVector[7]->push_back(v_pt);\r
- chargeVector[8]->push_back(v_E);\r
+ chargeVector[0]->push_back(vCharge);\r
+ chargeVector[1]->push_back(vY);\r
+ chargeVector[2]->push_back(vEta);\r
+ chargeVector[3]->push_back(vPhi);\r
+ chargeVector[4]->push_back(vP[0]);\r
+ chargeVector[5]->push_back(vP[1]);\r
+ chargeVector[6]->push_back(vP[2]);\r
+ chargeVector[7]->push_back(vPt);\r
+ chargeVector[8]->push_back(vE);\r
\r
// ------------------------------------------------------------- \r
// for each track in main event loop over all tracks in mix event\r
fHistTrackStats->Fill(aodTrackMix->GetFilterMap());\r
if(!aodTrackMix->TestFilterBit(nAODtrackCutBit)) continue;\r
\r
- v_charge = aodTrackMix->Charge();\r
- v_y = aodTrackMix->Y();\r
- v_eta = aodTrackMix->Eta();\r
- v_phi = aodTrackMix->Phi() * TMath::RadToDeg();\r
- v_E = aodTrackMix->E();\r
- v_pt = aodTrackMix->Pt();\r
- aodTrackMix->PxPyPz(v_p);\r
+ vCharge = aodTrackMix->Charge();\r
+ vY = aodTrackMix->Y();\r
+ vEta = aodTrackMix->Eta();\r
+ vPhi = aodTrackMix->Phi() * TMath::RadToDeg();\r
+ vE = aodTrackMix->E();\r
+ vPt = aodTrackMix->Pt();\r
+ aodTrackMix->PxPyPz(vP);\r
\r
- Float_t DCAxyMix = aodTrackMix->DCA(); // this is the DCA from global track (not exactly what is cut on)\r
- Float_t DCAzMix = aodTrackMix->ZAtDCA(); // this is the DCA from global track (not exactly what is cut on)\r
+ Float_t dcaXYMix = aodTrackMix->DCA(); // this is the DCA from global track (not exactly what is cut on)\r
+ Float_t dcaZMix = aodTrackMix->ZAtDCA(); // this is the DCA from global track (not exactly what is cut on)\r
\r
\r
// Kinematics cuts from ESD track cuts\r
- if( v_pt < fPtMin || v_pt > fPtMax) continue;\r
- if( v_eta < fEtaMin || v_eta > fEtaMax) continue;\r
+ if( vPt < fPtMin || vPt > fPtMax) continue;\r
+ if( vEta < fEtaMin || vEta > fEtaMax) continue;\r
\r
// Extra DCA cuts (for systematic studies [!= -1])\r
if( fDCAxyCut != -1 && fDCAxyCut != -1){\r
- if(TMath::Sqrt((DCAxyMix*DCAxyMix)/(fDCAxyCut*fDCAxyCut)+(DCAzMix*DCAzMix)/(fDCAzCut*fDCAzCut)) > 1 ){\r
+ if(TMath::Sqrt((dcaXYMix*dcaXYMix)/(fDCAxyCut*fDCAxyCut)+(dcaZMix*dcaZMix)/(fDCAzCut*fDCAzCut)) > 1 ){\r
continue; // 2D cut\r
}\r
}\r
\r
// fill QA histograms\r
fHistClus->Fill(aodTrackMix->GetITSNcls(),aodTrackMix->GetTPCNcls());\r
- fHistDCA->Fill(DCAzMix,DCAxyMix);\r
+ fHistDCA->Fill(dcaZMix,dcaXYMix);\r
fHistChi2->Fill(aodTrackMix->Chi2perNDF());\r
- fHistPt->Fill(v_pt);\r
- fHistEta->Fill(v_eta);\r
- fHistPhi->Fill(v_phi);\r
+ fHistPt->Fill(vPt);\r
+ fHistEta->Fill(vEta);\r
+ fHistPhi->Fill(vPhi);\r
\r
// fill charge vector\r
- chargeVector[0]->push_back(v_charge);\r
- chargeVector[1]->push_back(v_y);\r
- chargeVector[2]->push_back(v_eta);\r
- chargeVector[3]->push_back(v_phi);\r
- chargeVector[4]->push_back(v_p[0]);\r
- chargeVector[5]->push_back(v_p[1]);\r
- chargeVector[6]->push_back(v_p[2]);\r
- chargeVector[7]->push_back(v_pt);\r
- chargeVector[8]->push_back(v_E);\r
+ chargeVector[0]->push_back(vCharge);\r
+ chargeVector[1]->push_back(vY);\r
+ chargeVector[2]->push_back(vEta);\r
+ chargeVector[3]->push_back(vPhi);\r
+ chargeVector[4]->push_back(vP[0]);\r
+ chargeVector[5]->push_back(vP[1]);\r
+ chargeVector[6]->push_back(vP[2]);\r
+ chargeVector[7]->push_back(vPt);\r
+ chargeVector[8]->push_back(vE);\r
\r
\r
\r
}\r
\r
AliMixInputEventHandler *AliAnalysisTaskEventMixingBF::SetupEventsForMixing() {\r
+ //sets the input handlers for event mixing\r
\r
AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();\r
AliMultiInputEventHandler *inEvHMain = dynamic_cast<AliMultiInputEventHandler *>(mgr->GetInputEventHandler());\r