]> git.uio.no Git - u/mrichter/AliRoot.git/commitdiff
remove debugging printouts
authormverweij <marta.verweij@cern.ch>
Mon, 24 Feb 2014 09:57:39 +0000 (10:57 +0100)
committermverweij <marta.verweij@cern.ch>
Mon, 24 Feb 2014 10:07:30 +0000 (11:07 +0100)
remove debugging printouts

remove debugging printouts

PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskEmcalJetMass.cxx
PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskEmcalJetMassResponse.cxx

index 983a2695989c4abeda6951ab5d47c094c8daa4b4..fc998d62b101521e546c9881105b39fb5f37b87f 100644 (file)
@@ -254,10 +254,7 @@ Bool_t AliAnalysisTaskEmcalJetMass::FillHistograms()
 {
   // Fill histograms.
 
-  AliInfo(Form("%s",GetName()));
-
   AliEmcalJet* jet1 = NULL;
-
   AliJetContainer *jetCont = GetJetContainer(fContainerBase);
   if(jetCont) {
     jetCont->ResetCurrentID();
@@ -265,7 +262,7 @@ Bool_t AliAnalysisTaskEmcalJetMass::FillHistograms()
       Double_t fraction = jetCont->GetFractionSharedPt(jet1);
       if(fMinFractionShared>0. && fraction<fMinFractionShared) continue;
 
-      Double_t ptJet1 = jet1->Pt() - GetRhoVal(fContainerBase)*jet1->Area();//jetCont->GetJetPtCorr(jetCont->GetCurrentID());//
+      Double_t ptJet1 = jet1->Pt() - GetRhoVal(fContainerBase)*jet1->Area();
       fh2PtJet1VsLeadPtAllSel[fCentBin]->Fill(ptJet1,jet1->MaxTrackPt());
       fh2PtVsMassJet1All[fCentBin]->Fill(ptJet1,jet1->M());
       fpPtVsMassJet1All[fCentBin]->Fill(ptJet1,jet1->M());
@@ -287,43 +284,16 @@ Bool_t AliAnalysisTaskEmcalJetMass::FillHistograms()
       if((Et*jetCont->GetJetRadius())>0.) 
        massOverEtR = jet1->M()/(Et*jetCont->GetJetRadius());
       fh2EtMassOverEtRSq[fCentBin]->Fill(Et,massOverEtR*massOverEtR);
-
-      //check if matched gen level jet exists and do analysis
-      AliEmcalJet *jetGen = jet1->ClosestJet();
-
-      Printf("AA jet");
-      Printf("jet 4-vector: %f,%f,%f,%f",jet1->Px(),jet1->Py(),jet1->Pz(),jet1->E());
-      Printf("pT: %f  pTcorr: %f  M: %f",jet1->Pt(),ptJet1,jet1->M());
-      Printf("eta: %f  phi: %f",jet1->Eta(),jet1->Phi());
-      if(jetGen) {
-       Printf("gen jet");
-       Printf("jet 4-vector: %f,%f,%f,%f",jetGen->Px(),jetGen->Py(),jetGen->Pz(),jetGen->E());
-       Printf("pT: %f  M: %f",jetGen->Pt(),jetGen->M());
-       Printf("eta: %f  phi: %f",jetGen->Eta(),jetGen->Phi());
-      }
-
     }
   }
-
+  
   return kTRUE;
 }
 
 //________________________________________________________________________
 Double_t AliAnalysisTaskEmcalJetMass::GetJetMass(AliEmcalJet *jet) {
-
-  Double_t deltaM = jet->M() - fJetMassAvg;
-  Double_t scale = deltaM / jet->M();
-
-  Double_t mt2 = jet->E()*jet->E() - jet->Pt()*jet->Pt();
-  Double_t et2 = jet->M()*jet->M() + jet->Pt()*jet->Pt();
-
-  Double_t pxScale = jet->Px()*scale;
-  Double_t pyScale = jet->Py()*scale;
-  Double_t pzScale = jet->Pz()*scale;
-  
-  // Printf("scaled jet 4-vector: %f-%f-%f-%f",pxScale,pyScale,pzScale);
-
-  return deltaM;
+  //calc subtracted jet mass, not implemented
+  return jet->M();
 
 }
 
index 724185a180faa7744ee02b369f9bf66b49b773d9..e6856a988549590874e9cd8ae4109156549e8b0a 100644 (file)
@@ -223,10 +223,7 @@ Bool_t AliAnalysisTaskEmcalJetMassResponse::FillHistograms()
 {
   // Fill histograms.
 
-  AliInfo(Form("%s",GetName()));
-
   AliEmcalJet* jet1 = NULL;
-
   AliJetContainer *jetCont = GetJetContainer(fContainerBase);
   if(jetCont) {
     jetCont->ResetCurrentID();