if(fTCARandomConesOut)fTCARandomConesOut->Delete();
delete fTCARandomConesOut;
- //if(fTCARandomConesOutRan)fTCARandomConesOutRan->Delete();
- //delete fTCARandomConesOutRan;
+ if(fTCARandomConesOutRan)fTCARandomConesOutRan->Delete();
+ delete fTCARandomConesOutRan;
}
fGhostEtamax(1.5),
background(0),
fTCARandomConesOut(0x0),
- //fTCARandomConesOutRan(0x0),
+ fTCARandomConesOutRan(0x0),
fRandom(0),
fh1Xsec(0x0),
fh1Trials(0x0),
fh1DeltapT(0x0),
fh1Rho(0x0),
fh1PtRandCone(0x0),
- //fh1Area(0x0),
+ fh1Area(0x0),
fh1efficiencyPt(0x0),
fh2efficiencyPhi(0x0),
fh1ZPhySel(0x0),
fGhostEtamax(1.5),
background(0),
fTCARandomConesOut(0x0),
- //fTCARandomConesOutRan(0x0),
+ fTCARandomConesOutRan(0x0),
fRandom(0),
fh1Xsec(0x0),
fh1Trials(0x0),
fh1DeltapT(0x0),
fh1Rho(0x0),
fh1PtRandCone(0x0),
- //fh1Area(0x0),
+ fh1Area(0x0),
fh1efficiencyPt(0x0),
fh2efficiencyPhi(0x0),
fh1ZPhySel(0x0),
fTCARandomConesOut->SetName(cName.Data());
AddAODBranch("TClonesArray",&fTCARandomConesOut,fNonStdFile.Data());
}
- /*
// create the branch with the random for the random cones on the random event
cName = Form("%sRandomCone_random",fNonStdBranch.Data());
if(!AODEvent()->FindListObject(cName.Data())){
fTCARandomConesOutRan->SetName(cName.Data());
AddAODBranch("TClonesArray",&fTCARandomConesOutRan,fNonStdFile.Data());
}
- */
+
if(fNonStdFile.Length()!=0){
//
// case that we have an AOD extension we need to fetch the jets from the extended output
fh1DeltapT = new TH1F("fh1DeltapT","DeltapT",100,-50,50);
fh1Rho = new TH1F("fh1Rho","Rho",100,0,200);
fh1PtRandCone = new TH1F("fh1PtRandCone","pt",100,0,200);
- //fh1Area = new TH1F("fh1Area","area",100,0,1);
+ fh1Area = new TH1F("fh1Area","area",100,0,1);
const Int_t saveLevel = 3; // large save level more histos
if(saveLevel>0){
fHistList->Add(fh1DeltapT);
fHistList->Add(fh1Rho);
fHistList->Add(fh1PtRandCone);
- //fHistList->Add(fh1Area);
+ fHistList->Add(fh1Area);
}
// =========== Switch on Sumw2 for all histos ===========
// handle and reset the output jet branch
if(fTCARandomConesOut)fTCARandomConesOut->Delete();
- //if(fTCARandomConesOutRan)fTCARandomConesOutRan->Delete();
+ if(fTCARandomConesOutRan)fTCARandomConesOutRan->Delete();
//
// Execute analysis for current event
// create a random jet within the acceptance
Double_t etaMax = fTrackEtaWindow - fRparam;//0.9 - 0.4
Int_t nCone = 0;
- // Int_t nConeRan = 0;
+ Int_t nConeRan = 0;
Double_t pTC = 1; // small number
Double_t etaC = etaMax*2.*(fRandom->Rndm()-0.5); // +- etamax
Double_t phiC = fRandom->Rndm()*2.*TMath::Pi(); // 0 - 2pi
tmpRecC.SetBgEnergy(0,0); // this is use as temporary storage of the summed p_T below
if(fTCARandomConesOut)new ((*fTCARandomConesOut)[nCone++]) AliAODJet(tmpRecC);
- //if(fTCARandomConesOutRan)new ((*fTCARandomConesOutRan)[nConeRan++]) AliAODJet(tmpRecC);
+ if(fTCARandomConesOutRan)new ((*fTCARandomConesOutRan)[nConeRan++]) AliAODJet(tmpRecC);
// loop over the reconstructed particles and add up the pT in the random cones
// maybe better to loop over randomized particles not in the real jets...
jC->SetBgEnergy(jC->ChargedBgEnergy()+vp->Pt(),0);
}
}// add up energy in cone
- /*
+
// the randomized input changes eta and phi, but keeps the p_T
if(fTCARandomConesOutRan){
Double_t pTR = vp->Pt();
jC->SetBgEnergy(jC->ChargedBgEnergy()+vTmpRanR.Pt(),0);
}
}
- */
}// loop over recparticles
} //fTCARandomConesOut
Float_t jetArea = fRparam*fRparam*TMath::Pi();
rC->SetEffArea(jetArea,0);
}
}
- /*
if(fTCARandomConesOutRan){
AliAODJet* rC = (AliAODJet*)fTCARandomConesOutRan->At(0);
// same with random
rC->SetEffArea(jetArea,0);
}
}//find the random jets
- */
}//inclusive Jets > 0
//Calculate delta pT
fh1DeltapT->Fill(ptSub);// delta pT
fh1Rho->Fill(background);// background rho
fh1PtRandCone->Fill(randCone->Pt());// pT of random cone
- //fh1Area->Fill(randCone->EffectiveAreaCharged()); // area of random cone; should always be pi*0.4^2 = 0.5
+ fh1Area->Fill(randCone->EffectiveAreaCharged()); // area of random cone; should always be pi*0.4^2 = 0.5
}else{
if(fDebug)Printf("%s:%d No random Cone found",(char*)__FILE__,__LINE__);
}
if (fDebug > 2){
if(fTCARandomConesOut)Printf("%s:%d RC %d",(char*)__FILE__,__LINE__,fTCARandomConesOut->GetEntriesFast());
- //if(fTCARandomConesOutRan)Printf("%s:%d RC Ran %d",(char*)__FILE__,__LINE__,fTCARandomConesOutRan->GetEntriesFast());
+ if(fTCARandomConesOutRan)Printf("%s:%d RC Ran %d",(char*)__FILE__,__LINE__,fTCARandomConesOutRan->GetEntriesFast());
}
PostData(1, fHistList);
}