fPtAssocMin(0),
fPtAssocMax(0),
fAlpha(alpha),
+ fUseLeadingPt(1),
fMergeCount(0)
{
// Constructor
Int_t triggerAway_entries = triggerAway->GetEntriesFast();
AliVParticle* found_particle[triggerAway_entries];
+ Bool_t do_not_use_T1 = false;
+
+ //in case only the leading pt of a jet should be used, check every particle on the trigger near side if it's closer than alpha and if it has a higher pt than trigger 1
+ if(fUseLeadingPt){
+ for (Int_t i2=0; i2<triggerNear->GetEntriesFast(); i2++){
+ if(i==i2)
+ continue;
+
+ AliVParticle* part_i2 = (AliVParticle*) triggerNear->UncheckedAt(i2);
+
+ if(part_i2->Pt()<=part_pt)
+ continue;
+
+ Double_t dphi_check = part_phi-part_i2->Phi();
+ if(dphi_check>1.5*TMath::Pi()) dphi_check -= TMath::TwoPi();
+ else if(dphi_check<-0.5*TMath::Pi()) dphi_check += TMath::TwoPi();
+
+ if(TMath::Abs(dphi_check)<fAlpha){
+ do_not_use_T1 = true;
+ break;
+ }
+ }
+ }
+
+ //if there is a particle with higher energy than T1 closer than alpha to T1, do not use this T1
+ if(do_not_use_T1)
+ continue;
+
//have to fake the away side triggers for the 1+1 analysis
if(is1plus1){
found_particle[ind_found] = part;//in 1plus1 use first trigger particle also as pseudo second trigger particle
AliVParticle* part2 = (AliVParticle*) triggerAway->UncheckedAt(j);
Double_t part2_pt = part2->Pt();
- //check if pT of trigger 2 is within the trigger range
- //also pt of trigger 2 needs to be smaller than the pt of trigger 1 (to have an ordering if both pt are close to each other)
- if(part2_pt<fTriggerPt2Min || part2_pt>fTriggerPt2Max || part2_pt>part_pt)
+ //check if pT of trigger 2 has enough energy to be a trigger
+ //maximum energy is checked later after checking this particle may have to much energy for trigger 1 or 2
+ if(part2_pt<fTriggerPt2Min)
continue;
-
+
// don't use the same particle (is in any case impossible because the Delta phi angle will be 0)
if(part==part2){
continue;
}
Double_t dphi_triggers = part_phi-part2->Phi();
-
if(dphi_triggers>1.5*TMath::Pi()) dphi_triggers -= TMath::TwoPi();
else if(dphi_triggers<-0.5*TMath::Pi()) dphi_triggers += TMath::TwoPi();
if(TMath::Abs(dphi_triggers)>fAlpha)
continue;
+ //check if pT of trigger 2 is too high
+ if(part2_pt>fTriggerPt2Max || part2_pt>part_pt){
+ //pt of trigger 2 needs to be smaller than the pt of trigger 1 (to have an ordering if both pt are close to each other)
+ if(fUseLeadingPt){
+ do_not_use_T1 = true;
+ break;
+ }else
+ continue;
+ }
+
found_particle[ind_found] = part2;
if(ind_max_found_pt==-1 || part2_pt>found_particle[ind_max_found_pt]->Pt()) ind_max_found_pt = ind_found;
ind_found++;
}//end loop to search for the second trigger particle
}
+ //if there is a particle with higher energy than T1 or max(T2) within Delta phi = pi +/- alpha to T1, do not use this T1
+ if(do_not_use_T1)
+ continue;
+
//if no second trigger particle was found continue to search for the next first trigger particle
if(ind_found==0)
continue;
+ //use only the highest energetic particle on the away side, if there is only 1 away side trigger this is already the case
+ if(fUseLeadingPt && ind_found>1){
+ found_particle[0] = found_particle[ind_max_found_pt];
+ ind_found=1;
+ ind_max_found_pt = 0;
+ }
+
//the energy of the second trigger particle is set for the near side to the maximum energy of all trigger 2 particles on the away side
// this leads to the fact that the number of accepted trigger combinations can be artificial smaller than the real number if there is a cut on the pT 2 energy from the top; cutting away the smallest energy of pT 2 is still save; this is the reason why it is not allowed to use a cut on the top pt of trigger particle 2
//fill trigger particles
if(ind_found>0){
- Double_t vars[4];
- vars[0] = part_pt;
- vars[1] = centrality;
- vars[2] = zVtx;
- vars[3] = found_particle[ind_max_found_pt]->Pt();
- if(is1plus1)
- vars[3] = (fTriggerPt2Max+fTriggerPt2Min)/2;
-
- event_hist->Fill(vars, stepUEHist, weight);//near side
-
- if(!is1plus1)
- for(Int_t k=0; k< ind_found; k++){
- vars[3] = found_particle[k]->Pt();
- event_hist->Fill(vars, stepUEHist+1, weight);//away side
- }
+ Double_t vars[4];
+ vars[0] = part_pt;
+ vars[1] = centrality;
+ vars[2] = zVtx;
+ vars[3] = found_particle[ind_max_found_pt]->Pt();
+ if(is1plus1)
+ vars[3] = (fTriggerPt2Max+fTriggerPt2Min)/2;
+
+ event_hist->Fill(vars, stepUEHist, weight);//near side
+
+ if(!is1plus1)
+ for(Int_t k=0; k< ind_found; k++){
+ vars[3] = found_particle[k]->Pt();
+ event_hist->Fill(vars, stepUEHist+1, weight);//away side
+ }
- //fill fTriggerPt only once, choosed kSameNS
- if(step==AliTwoPlusOneContainer::kSameNS)
- for(Int_t k=0; k< ind_found; k++)
- fTriggerPt->Fill(part_pt, found_particle[k]->Pt());
-
- //fill asymmetry only for kSameNS and kMixedNS
- if(step==AliTwoPlusOneContainer::kSameNS||step==AliTwoPlusOneContainer::kMixedNS){
- for(Int_t k=0; k< ind_found; k++){
- Float_t asymmetry = (part_pt-found_particle[k]->Pt())/(part_pt+found_particle[k]->Pt());
- if(step==AliTwoPlusOneContainer::kSameNS){
- fAsymmetry->Fill(asymmetry);
- }else{
- fAsymmetryMixed->Fill(asymmetry);
- }
+ //fill fTriggerPt only once, choosed kSameNS
+ if(step==AliTwoPlusOneContainer::kSameNS)
+ for(Int_t k=0; k< ind_found; k++)
+ fTriggerPt->Fill(part_pt, found_particle[k]->Pt());
+
+ //fill asymmetry only for kSameNS and kMixedNS
+ if(step==AliTwoPlusOneContainer::kSameNS||step==AliTwoPlusOneContainer::kMixedNS){
+ for(Int_t k=0; k< ind_found; k++){
+ Float_t asymmetry = (part_pt-found_particle[k]->Pt())/(part_pt+found_particle[k]->Pt());
+ if(step==AliTwoPlusOneContainer::kSameNS){
+ fAsymmetry->Fill(asymmetry);
+ }else{
+ fAsymmetryMixed->Fill(asymmetry);
}
}
-
+ }
}
-
+
//add correlated particles on the near side
for (Int_t k=0; k<assocNear->GetEntriesFast(); k++){
AliVParticle* part3 = (AliVParticle*) assocNear->UncheckedAt(k);
}//end loop to search for the first trigger particle
}
-
-
//____________________________________________________________________
AliTwoPlusOneContainer &AliTwoPlusOneContainer::operator=(const AliTwoPlusOneContainer &c)
{
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff