-
// ******************************************
// This task computes several jet observables like
// the fraction of energy in inner and outer coronnas,
fCentMax(100.),
fNInputTracksMin(0),
fNInputTracksMax(-1),
+fRequireITSRefit(0),
+fApplySharedClusterCut(0),
fAngStructCloseTracks(0),
fCheckMethods(0),
fDoEventMixing(0),
fFlagPhiBkg(0),
fFlagEtaBkg(0),
fFlagJetHadron(0),
-fTTLow(11),
-fTTUp(13),
+fFrac(0.8),
+fTTLowRef(11),
+fTTUpRef(13),
+fTTLowSig(15),
+fTTUpSig(19),
fHardest(0),
fFlagRandom(0),
fFlagOnlyRecoil(0),
fh2AngStructpt2C60(0x0),
fh2AngStructpt3C60(0x0),
fh2AngStructpt4C60(0x0),
+fh1TrigRef(0x0),
+fh1TrigSig(0x0),
fh2Ntriggers(0x0),
fh2Ntriggers2C10(0x0),
fh2Ntriggers2C20(0x0),
fCentMax(100.),
fNInputTracksMin(0),
fNInputTracksMax(-1),
+fRequireITSRefit(0),
+fApplySharedClusterCut(0),
fAngStructCloseTracks(0),
fCheckMethods(0),
fDoEventMixing(0),
fFlagPhiBkg(0),
fFlagEtaBkg(0),
fFlagJetHadron(0),
-fTTLow(11),
-fTTUp(13),
+fFrac(0.8),
+fTTLowRef(11),
+fTTUpRef(13),
+fTTLowSig(15),
+fTTUpSig(19),
fHardest(0),
fFlagRandom(0),
fFlagOnlyRecoil(0),
fh2AngStructpt2C60(0x0),
fh2AngStructpt3C60(0x0),
fh2AngStructpt4C60(0x0),
+fh1TrigRef(0x0),
+fh1TrigSig(0x0),
fh2Ntriggers(0x0),
fh2Ntriggers2C10(0x0),
fh2Ntriggers2C20(0x0),
fh2AngStructpt4C60 = new TH2F("Ang struct pt4 C60","",15,0.,1.5,150,0.,10.); }
-
+ fh1TrigRef=new TH1D("Trig Ref","",10,0.,10);
+ fh1TrigSig=new TH1D("Trig Sig","",10,0.,10);
fh2Ntriggers=new TH2F("# of triggers","",100,0.,100.,50,0.,50.);
fh2Ntriggers2C10=new TH2F("# of triggers2C10","",50,0.,50.,50,0.,50.);
fh2Ntriggers2C20=new TH2F("# of triggers2C20","",50,0.,50.,50,0.,50.);
-
+ fOutputList->Add(fh1TrigRef);
+ fOutputList->Add(fh1TrigSig);
fOutputList->Add(fh2Ntriggers);
fOutputList->Add(fh2Ntriggers2C10);
fOutputList->Add(fh2Ntriggers2C20);
const Int_t dimSpec = 5;
const Int_t nBinsSpec[dimSpec] = {100,6, 140, 50, fNRPBins};
const Double_t lowBinSpec[dimSpec] = {0,0,-80, 0, 0};
- const Double_t hiBinSpec[dimSpec] = {100,1, 200, 50, fNRPBins};
+ const Double_t hiBinSpec[dimSpec] = {100,1, 200, 50, static_cast<Double_t>(fNRPBins)};
fHJetSpec = new THnSparseF("fHJetSpec","Recoil jet spectrum",dimSpec,nBinsSpec,lowBinSpec,hiBinSpec);
//change binning in jet area
// -- event selection --
fHistEvtSelection->Fill(1); // number of events before event selection
- // physics selection
+
+ Bool_t selected=kTRUE;
+ selected = AliAnalysisHelperJetTasks::Selected();
+ if(!selected){
+ // no selection by the service task, we continue
+ PostData(1,fOutputList);
+ return;}
+
+
+
+ // physics selection: this is now redundant, all should appear as accepted after service task selection
AliInputEventHandler* inputHandler = (AliInputEventHandler*)
((AliAnalysisManager::GetAnalysisManager())->GetInputEventHandler());
std::cout<<inputHandler->IsEventSelected()<<" "<<fOfflineTrgMask<<std::endl;
return;
}
+
+
// vertex selection
if(!aod){
if(fDebug) Printf("%s:%d No AOD",(char*)__FILE__,__LINE__);
if(fIsPbPb){
if(fESD) {cent = fESD->GetCentrality();
if(cent) centValue = cent->GetCentralityPercentile("V0M");}
- else centValue=aod->GetHeader()->GetCentrality();
+ else centValue=((AliVAODHeader*)aod->GetHeader())->GetCentrality();
if(fDebug) printf("centrality: %f\n", centValue);
if (centValue < fCentMin || centValue > fCentMax){
if(fFlagRandom==0){
if(externalBackground)rho = externalBackground->GetBackground(0);}
- if(fFlagRandom==1){
+ if(fFlagRandom==2){
if(externalBackground)rho = externalBackground->GetBackground(2);}
-
+ if(fFlagRandom==3){
+ if(externalBackground)rho = externalBackground->GetBackground(3);}
// fetch jets
TClonesArray *aodJets[2];
aodJets[0]=0;
aodJets[1] = dynamic_cast<TClonesArray*>(fAODIn->FindListObject(fJetBranchName[1].Data())); }
- //Double_t ptsub[aodJets[0]->GetEntriesFast()];
- //Int_t inord[aodJets[0]->GetEntriesFast()];
- //for(Int_t n=0;n<aodJets[0]->GetEntriesFast();n++){
- // ptsub[n]=0;
- // inord[n]=0;}
+
Int_t nT=0;
TList ParticleList;
- if(fHardest==1) nT = GetListOfTracks(&ParticleList);
- if(fHardest==0) nT=SelectTrigger(&ParticleList);
+ Double_t minT=0;
+ Double_t maxT=0;
+ Int_t number=0;
+ Double_t dice=fRandom->Uniform(0,1);
+ if(dice>fFrac){ minT=fTTLowRef;
+ maxT=fTTUpRef;}
+ if(dice<=fFrac){minT=fTTLowSig;
+ maxT=fTTUpSig;}
+
+
+
+ if(fHardest==1 || fHardest==2) nT = GetListOfTracks(&ParticleList);
+ if(fHardest==0) nT=SelectTrigger(&ParticleList,minT,maxT,number);
+ if(nT<0){
+ PostData(1, fOutputList);
+ return;}
+
+ if(dice>fFrac) fh1TrigRef->Fill(number);
+ if(dice<=fFrac)fh1TrigSig->Fill(number);
+
for (Int_t iJetType = 0; iJetType < 2; iJetType++) {
fListJets[iJetType]->Clear();
if (!aodJets[iJetType]) continue;
Double_t phibig=0.;
Double_t etasmall=0;
Double_t ptsmall=0;
- Double_t areasmall=0;
+ // Double_t areasmall=0;
Double_t phismall=0.;
Int_t iCount=0;
Int_t trigJet=-1;
Int_t trigBBTrack=-1;
- Int_t trigInTrack=-1;
- fRPAngle = aod->GetHeader()->GetEventplane();
-
-
+ // Int_t trigInTrack=-1;
+ fRPAngle = ((AliVAODHeader*)aod->GetHeader())->GetEventplane();
+ if(fHardest==0 || fHardest==1){
AliVParticle *partback = (AliVParticle*)ParticleList.At(nT);
if(!partback){
PostData(1, fOutputList);
return;}
+ if(fSemigoodCorrect){
+ Double_t disthole=RelativePhi(partback->Phi(),fHolePos);
+ if(TMath::Abs(disthole)+fHoleWidth>TMath::Pi()-0.6){
+ PostData(1, fOutputList);
+ return;}}
+
+
+ }
+
+
+ for(Int_t tt=0;tt<ParticleList.GetEntries();tt++){
+ if(fHardest==0||fHardest==1){if(tt!=nT) continue;}
+ AliVParticle *partback = (AliVParticle*)ParticleList.At(tt);
+ if(!partback) continue;
+ if(partback->Pt()<8) continue;
- //for(Int_t tt=0;tt<ParticleList.GetEntries();tt++){
- //if(fFlagOnlyHardest!=0){if(tt!=nT) continue;}
- //AliVParticle *partback = (AliVParticle*)ParticleList.At(tt);
Double_t accep=2.*TMath::Pi()*1.8;
Int_t injet4=0;
Int_t injet=0;
+
if(fSemigoodCorrect){
Double_t disthole=RelativePhi(partback->Phi(),fHolePos);
- if(TMath::Abs(disthole)+fHoleWidth>TMath::Pi()-0.6){
- PostData(1, fOutputList);
- return;}
+ if(TMath::Abs(disthole)+fHoleWidth>TMath::Pi()-0.6) continue;}
- }
+
fh2Ntriggers->Fill(centValue,partback->Pt());
Double_t phiBinT = RelativePhi(partback->Phi(),fRPAngle);
if(phitt<0)phitt+=TMath::Pi()*2.;
Int_t phiBintt = GetPhiBin(phitt-fRPAngle);
- Double_t fillspec[] = {centValue,jetbig->EffectiveAreaCharged(),ptcorr,partback->Pt(),phiBintt};
+ Double_t fillspec[] = {centValue,jetbig->EffectiveAreaCharged(),ptcorr,partback->Pt(), static_cast<Double_t>(phiBintt)};
fHJetSpec->Fill(fillspec);
if(iCount==0){
trigJet=i;
trigBBTrack=nT;
- trigInTrack=ippt;
+ // trigInTrack=ippt;
iCount=iCount+1;}
etasmall = jetsmall->Eta();
phismall = jetsmall->Phi();
ptsmall = jetsmall->Pt();
- areasmall = jetsmall->EffectiveAreaCharged();
+ // areasmall = jetsmall->EffectiveAreaCharged();
Double_t tmpDeltaR=(phismall-phibig)*(phismall-phibig)+(etasmall-etabig)*(etasmall-etabig);
tmpDeltaR=TMath::Sqrt(tmpDeltaR);
//Fraction in the jet core
-
+ }
for(int it = 0;it < aod->GetNumberOfTracks();++it){
- AliAODTrack *tr = aod->GetTrack(it);
+ AliAODTrack *tr = dynamic_cast<AliAODTrack*>(aod->GetTrack(it));
+ if(!tr) AliFatal("Not a standard AOD");
Bool_t bGood = false;
if(fFilterType == 0)bGood = true;
else if(fFilterType == 1)bGood = tr->IsHybridTPCConstrainedGlobal();
else if(fFilterType == 2)bGood = tr->IsHybridGlobalConstrainedGlobal();
- if((fFilterMask>0)&&!(tr->TestFilterBit(fFilterMask)))continue;
+ if((fFilterMask>0)&&!(tr->TestFilterBit(fFilterMask)))continue;
+ if(fRequireITSRefit==1){if((tr->GetStatus()&AliESDtrack::kITSrefit)==0)continue;}
if(bGood==false) continue;
- if(TMath::Abs(tr->Eta())>0.9)continue;
+ if (fApplySharedClusterCut) {
+ Double_t frac = Double_t(tr->GetTPCnclsS()) /Double_t(tr->GetTPCncls());
+ if (frac > 0.4) continue;
+ }
+ if(TMath::Abs(tr->Eta())>0.9)continue;
if(tr->Pt()<0.15)continue;
list->Add(tr);
iCount++;
- Int_t AliAnalysisTaskJetCore::SelectTrigger(TList *list){
+Int_t AliAnalysisTaskJetCore::SelectTrigger(TList *list,Double_t minT,Double_t maxT,Int_t &number){
Int_t iCount = 0;
AliAODEvent *aod = 0;
if(!fESD)aod = fAODIn;
else aod = fAODOut;
if(!aod)return 0;
Int_t index=-1;
- Double_t triggers[100];
- for(Int_t cr=0;cr<100;cr++){triggers[cr]=-1;}
+ Int_t triggers[100];
+ Int_t triggers2[100];
+ for(Int_t cr=0;cr<100;cr++){triggers[cr]=-1;
+ triggers2[cr]=-1;}
Int_t im=0;
+ Int_t im2=0;
for(int it = 0;it < aod->GetNumberOfTracks();++it){
- AliAODTrack *tr = aod->GetTrack(it);
+ AliAODTrack *tr = dynamic_cast<AliAODTrack*>(aod->GetTrack(it));
+ if(!tr) AliFatal("Not a standard AOD");
Bool_t bGood = false;
if(fFilterType == 0)bGood = true;
else if(fFilterType == 1)bGood = tr->IsHybridTPCConstrainedGlobal();
- else if(fFilterType == 2)bGood =
- tr->IsHybridGlobalConstrainedGlobal();
+ else if(fFilterType == 2)bGood = tr->IsHybridGlobalConstrainedGlobal();
if((fFilterMask>0)&&!(tr->TestFilterBit(fFilterMask)))continue;
+ if(fRequireITSRefit==1){if((tr->GetStatus()&AliESDtrack::kITSrefit)==0)continue;}
if(bGood==false) continue;
+ if (fApplySharedClusterCut) {
+ Double_t frac = Double_t(tr->GetTPCnclsS()) /Double_t(tr->GetTPCncls());
+ if (frac > 0.4) continue;
+ }
if(TMath::Abs(tr->Eta())>0.9)continue;
if(tr->Pt()<0.15)continue;
list->Add(tr);
iCount++;
- if(tr->Pt()>=fTTLow && tr->Pt()<fTTUp){
+ if(tr->Pt()>=minT && tr->Pt()<maxT){
triggers[im]=iCount-1;
- im=im+1;}
+ im=im+1;
+ if(tr->Pt()>=20.){triggers2[im2]=iCount-1;
+ im2=im2+1;}
- }
+ }}
+
+ number=im;
Int_t rd=0;
- if(im==0) rd=0;
- if(im>0) rd=fRandom->Integer(im-1);
- index=triggers[rd];
+ if(im2==0) rd=0;
+ if(im2>0) rd=fRandom->Integer(im2);
+ index=triggers2[rd];
+ AliVParticle *tr1 = (AliVParticle*)list->At(index);
+
+
+
+ for(Int_t kk=0;kk<im;kk++){
+ //if(kk==rd) continue;
+ if(index==triggers[kk]) continue;
+ Int_t lab=triggers[kk];
+ AliVParticle *tr2 = (AliVParticle*)list->At(lab);
+
+ Double_t detat=tr1->Eta()-tr2->Eta();
+ Double_t dphit=RelativePhi(tr1->Phi(),tr2->Phi());
+ Double_t deltaRt=TMath::Sqrt(detat*detat+dphit*dphit);
+
+ if(deltaRt>0.4) number=number-1;}
Int_t index=-1;
Double_t ptmax=-10;
Double_t dphi=0;
- Double_t dif=0;
+ // Double_t dif=0;
Int_t iCount=0;
for(int it = 0;it < aod->GetNumberOfTracks();++it){
- AliAODTrack *tr = aod->GetTrack(it);
+ AliAODTrack *tr = dynamic_cast<AliAODTrack*>(aod->GetTrack(it));
+ if(!tr) AliFatal("Not a standard AOD");
if((fFilterMask>0)&&!(tr->TestFilterBit(fFilterMask)))continue;
if(TMath::Abs(tr->Eta())>0.9)continue;
if(tr->Pt()<0.15)continue;
if(TMath::Abs(dphi)<TMath::Pi()-0.6) continue;
if(tr->Pt()>ptmax){ ptmax=tr->Pt();
index=iCount-1;
- dif=dphi; }}
+ // dif=dphi;
+ }}
return index;
else aod = fAODOut;
for(int it = 0;it < aod->GetNumberOfTracks();++it){
- AliAODTrack *tr = aod->GetTrack(it);
+ AliAODTrack *tr = dynamic_cast<AliAODTrack*>(aod->GetTrack(it));
+ if(!tr) AliFatal("Not a standard AOD");
if((fFilterMask>0)&&!(tr->TestFilterBit(fFilterMask)))continue;
if(TMath::Abs(tr->Eta())>0.9)continue;
if(tr->Pt()<0.15)continue;