Multiple high pT triggers study
authorlcunquei <lcunquei@cern.ch>
Sat, 29 Nov 2014 13:59:43 +0000 (14:59 +0100)
committerlcunquei <lcunquei@cern.ch>
Sat, 29 Nov 2014 13:59:43 +0000 (14:59 +0100)
PWGJE/AliAnalysisTaskJetCore.cxx
PWGJE/AliAnalysisTaskJetCore.h

index a403ea6..044112c 100644 (file)
@@ -89,6 +89,7 @@ fDoEventMixing(0),
 fFlagPhiBkg(0),
 fFlagEtaBkg(0),
 fFlagJetHadron(0),
+fDodiHadron(0),
 fFrac(0.8),
 fTTLowRef(11),
 fTTUpRef(13),
@@ -152,6 +153,8 @@ fh2AngStructpt3C60(0x0),
 fh2AngStructpt4C60(0x0),
 fh1TrigRef(0x0),
 fh1TrigSig(0x0),
+fh1TrackPhiDistance(0x0),
+fh1TrackRDistance(0x0), 
 fh2Ntriggers(0x0),
 fh2Ntriggers2C10(0x0),
 fh2Ntriggers2C20(0x0), 
@@ -220,6 +223,7 @@ fDoEventMixing(0),
 fFlagPhiBkg(0),
 fFlagEtaBkg(0),
 fFlagJetHadron(0),
+fDodiHadron(0),
 fFrac(0.8),
 fTTLowRef(11),
 fTTUpRef(13),
@@ -283,6 +287,8 @@ fh2AngStructpt3C60(0x0),
 fh2AngStructpt4C60(0x0),    
 fh1TrigRef(0x0),
 fh1TrigSig(0x0),
+fh1TrackPhiDistance(0x0),
+fh1TrackRDistance(0x0), 
 fh2Ntriggers(0x0),
 fh2Ntriggers2C10(0x0),
 fh2Ntriggers2C20(0x0),
@@ -436,6 +442,8 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects()
    
     fh1TrigRef=new TH1D("Trig Ref","",10,0.,10);
     fh1TrigSig=new TH1D("Trig Sig","",10,0.,10);  
+    fh1TrackPhiDistance=new TH1D("PhiDistance","",35,0.,3.5);
+    fh1TrackRDistance=new TH1D("RDistance","",10,0.,1); 
     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.);
@@ -498,6 +506,9 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects()
 
         fOutputList->Add(fh1TrigRef);
         fOutputList->Add(fh1TrigSig); 
+        fOutputList->Add(fh1TrackPhiDistance);
+        fOutputList->Add(fh1TrackRDistance);
+        fOutputList->Add(fh1TrigSig); 
        fOutputList->Add(fh2Ntriggers);
         fOutputList->Add(fh2Ntriggers2C10);
         fOutputList->Add(fh2Ntriggers2C20); 
@@ -1278,8 +1289,12 @@ Int_t  AliAnalysisTaskJetCore::SelectTrigger(TList *list,Double_t minT,Double_t
        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;}      
+       fh1TrackPhiDistance->Fill(TMath::Abs(dphit));
+       fh1TrackRDistance->Fill(deltaRt);
+       
+       if(fDodiHadron==1)  if(deltaRt>0.4) number=number-1;
+       if(fDodiHadron==2) if((deltaRt>0.4) && (TMath::Abs(dphit)>TMath::Pi()-0.6)) number=number-1;}
+
 
      
   
index 2d68707..667020d 100644 (file)
@@ -86,6 +86,7 @@ public:
    virtual void     SetFlagPhiBkg(Int_t yesno){fFlagPhiBkg=yesno;}
    virtual void     SetFlagEtaBkg(Int_t yesno){fFlagEtaBkg=yesno;}
    virtual void     SetFlagJetHadron(Int_t yesno){fFlagJetHadron=yesno;}
+   virtual void     SetFlagDodiHadron(Int_t yesno){fDodiHadron=yesno;}
    virtual void     SetFlagFrac(Float_t fraction){fFrac=fraction;} 
    virtual void     SetTTLowRef(Float_t ttlow){fTTLowRef=ttlow;}
    virtual void     SetTTUpRef(Float_t ttup){fTTUpRef=ttup;}
@@ -152,6 +153,7 @@ private:
    Int_t   fFlagPhiBkg;
    Int_t   fFlagEtaBkg;
    Int_t   fFlagJetHadron;
+   Int_t   fDodiHadron;
    Float_t fFrac;               //fraction of events to evaluate signal
    Float_t fTTLowRef;
    Float_t fTTUpRef;
@@ -221,6 +223,8 @@ private:
   
      TH1D*      fh1TrigRef;                 //ref multiple triggers
      TH1D*      fh1TrigSig;                 //sig multiple triggers    
+     TH1D*      fh1TrackPhiDistance;   //distance between high pT triggers
+     TH1D*      fh1TrackRDistance;     //distance between high pT triggers
      TH2F*      fh2Ntriggers;              //triggers
      TH2F*      fh2Ntriggers2C10;             //centrality bias of triggers 
      TH2F*      fh2Ntriggers2C20;             //centrality bias of triggers