fTrackTypeRec(kTrackUndef),
fTrackTypeGen(kTrackUndef),
fNSkipLeadingRan(0),
+ fNSkipLeadingCone(0),
fNRandomCones(0),
fAvgTrials(1),
fExternalWeight(1),
fTrackTypeRec(kTrackUndef),
fTrackTypeGen(kTrackUndef),
fNSkipLeadingRan(0),
+ fNSkipLeadingCone(0),
fNRandomCones(0),
fAvgTrials(1),
fExternalWeight(1),
}
}
// create the branch for the random cones with the same R
- TString cName = Form("%sRandomCone",fNonStdBranch.Data());
-
+ TString cName = Form("%sRandomConeSkip%02d",fNonStdBranch.Data(),fNSkipLeadingCone);
+
if(fNRandomCones>0){
if(!AODEvent()->FindListObject(cName.Data())){
fTCARandomConesOut = new TClonesArray("AliAODJet", 0);
Double_t pC = TMath::Sqrt(pTC*pTC+pZC*pZC);
AliAODJet tmpRecC (pXC,pYC,pZC, pC);
bool skip = false;
- for(int jj = 0; jj < TMath::Min(nRec,2);jj++){// test for overlap with leading jets
+ for(int jj = 0; jj < TMath::Min(nRec,fNSkipLeadingCone);jj++){// test for overlap with leading jets
AliAODJet jet (sortedJets[jj].px(), sortedJets[jj].py(), sortedJets[jj].pz(), sortedJets[jj].E());
if(jet.DeltaR(& tmpRecC)<2.*fRparam+0.2){
skip = true;
virtual void SetBackgroundBranch(const char* c){fBackgroundBranch = c;}
virtual const char* GetBackgroundBranch(){return fBackgroundBranch.Data();}
virtual void SetNSkipLeadingRan(Int_t x){fNSkipLeadingRan = x;}
+ virtual void SetNSkipLeadingCone(Int_t x){fNSkipLeadingCone = x;}
virtual void SetNRandomCones(Int_t x){fNRandomCones = x;}
virtual void SetJetOutputBranch(const char *c){fNonStdBranch = c;}
UInt_t fFilterMask; // filter bit for slecected tracks
Int_t fTrackTypeRec; // type of tracks used for FF
Int_t fTrackTypeGen; // type of tracks used for FF
- Int_t fNSkipLeadingRan; // number of leading tracks to be skipped in the randomized event
+ Int_t fNSkipLeadingRan; // number of leading tracks to be skipped in the randomized event
+ Int_t fNSkipLeadingCone; // number of leading jets to be for the random cones
Int_t fNRandomCones; // number of generated random cones
Float_t fAvgTrials; // Average nimber of trials
Float_t fExternalWeight; // external weight
TList *fHistList; //!leading tracks to be skipped in the randomized event Output list
- ClassDef(AliAnalysisTaskJetCluster, 14)
+ ClassDef(AliAnalysisTaskJetCluster, 16)
};
#endif