AliCFMuonResTask1.C :
authorrvernet <rvernet@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 3 Apr 2009 09:50:23 +0000 (09:50 +0000)
committerrvernet <rvernet@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 3 Apr 2009 09:50:23 +0000 (09:50 +0000)
 - AliCFManager (connect input)
 - ordering  of SetParticleCutsList
AliCFMuonResTask1.cxx :
 - ordering  of CheckParticleCuts
 - Phi is the same now in TParticle and AliESDMuonTrack

CORRFW/test/muons/AliCFMuonResTask1.C
CORRFW/test/muons/AliCFMuonResTask1.cxx

index 8cd40dc..6fc6991 100644 (file)
@@ -76,6 +76,16 @@ Bool_t AliCFMuonResTask1(
     else {
       analysisChain = new TChain("esdTree");
       analysisChain->Add("/scratch/lopez/PDC08jpsi/run2-300/AliESDs.root");
+      analysisChain->Add("/scratch/lopez/PDC08jpsi/run3-1000/AliESDs.root");
+      analysisChain->Add("/scratch/lopez/PDC08jpsi/run4-1000/AliESDs.root");
+      analysisChain->Add("/scratch/lopez/PDC08jpsi/run5-1000/AliESDs.root");
+      analysisChain->Add("/scratch/lopez/PDC08jpsi/run7-1000/AliESDs.root");
+      analysisChain->Add("/scratch/lopez/PDC08jpsi/run8-1000/AliESDs.root");
+      analysisChain->Add("/scratch/lopez/PDC08jpsi/run9-1000/AliESDs.root");
+      analysisChain->Add("/scratch/lopez/PDC08jpsi/run10-1000/AliESDs.root");
+      analysisChain->Add("/scratch/lopez/PDC08jpsi/run11-1000/AliESDs.root");
+      analysisChain->Add("/scratch/lopez/PDC08jpsi/run12-1000/AliESDs.root");
+      analysisChain->Add("/scratch/lopez/PDC08jpsi/run13-1000/AliESDs.root");
    }
   }
   
@@ -187,12 +197,13 @@ Bool_t AliCFMuonResTask1(
   recList->AddLast(recKineCuts);
 
   // CREATE THE INTERFACE TO CORRECTION FRAMEWORK USED IN THE TASK
-  printf("CREATE INTERFACE AND CUTS\n");
+
+   printf("CREATE INTERFACE AND CUTS\n");
   AliCFManager* man = new AliCFManager() ;
   man->SetParticleContainer     (container);
-  man->SetParticleCutsList(AliCFManager::kPartGenCuts,mcList);
-  man->SetParticleCutsList(AliCFManager::kPartRecCuts,recList);
 
+  man->SetParticleCutsList(AliCFManager::kPartGenCuts,mcList);
+  man->SetParticleCutsList(AliCFManager::kPartAccCuts,recList);
 
   //CREATE THE TASK
   printf("CREATE TASK\n");
@@ -237,7 +248,9 @@ Bool_t AliCFMuonResTask1(
 
   cinput0->SetData(analysisChain);
   mgr->AddTask(task);
-  mgr->ConnectInput(task,0,cinput0);
+
+  mgr->ConnectInput(task,0,mgr->GetCommonInputContainer());
+
   mgr->ConnectOutput(task,1,coutput1);
   mgr->ConnectOutput(task,2,coutput2);
 
@@ -258,10 +271,9 @@ void Load() {
 
   //load the required aliroot libraries
   gSystem->Load("libANALYSIS") ;
-
   gSystem->Load("libANALYSISalice") ;
 
-  //gSystem->Load("libCORRFW.so") ;
+//  gSystem->Load("libCORRFW.so") ;
   gSystem->Load("$ALICE_ROOT/lib/tgt_linux/libCORRFW.so") ;
 
   //compile online the task class
index 83ce494..19f9147 100644 (file)
@@ -175,14 +175,14 @@ void AliCFMuonResTask1::UserExec(Option_t *)
     part0 = stack->Particle(p0); \r
    // selection of the rapidity for first muon\r
     AliMCParticle *mcpart0 = new AliMCParticle(part0);\r
-    if (!fCFManager->CheckParticleCuts(AliCFManager::kPartRecCuts,mcpart0)) continue;\r
+    if (!fCFManager->CheckParticleCuts(AliCFManager::kPartAccCuts,mcpart0)) continue;\r
     Int_t pdg0 = part0->GetPdgCode();\r
 \r
     Int_t p1 = part->GetDaughter(1);\r
     part1 = stack->Particle(p1);\r
    // selection of the rapidity for second muon\r
     AliMCParticle *mcpart1 = new AliMCParticle(part1);\r
-    if (!fCFManager->CheckParticleCuts(AliCFManager::kPartRecCuts,mcpart1)) continue;\r
+    if (!fCFManager->CheckParticleCuts(AliCFManager::kPartAccCuts,mcpart1)) continue;\r
     Int_t pdg1 = part1->GetPdgCode();\r
  \r
    // 0 mu- 1 mu+\r
@@ -195,7 +195,7 @@ void AliCFMuonResTask1::UserExec(Option_t *)
        pz0 = part0->Pz();\r
        py0 = part0->Py();\r
        px0 = part0->Px();\r
-       phi0 = part0->Phi(); // Warning in TParticle Phi = pi + ATan2(Py,Px) = [0,2pi] \r
+       phi0 = part0->Phi(); \r
        phi0 = Phideg(phi0);    \r
        rapmc0=Rap(e0,pz0);\r
 \r
@@ -258,9 +258,9 @@ void AliCFMuonResTask1::UserExec(Option_t *)
   AliESDInputHandler *esdH = dynamic_cast<AliESDInputHandler*>\r
       (AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());\r
   fESD = esdH->GetEvent();\r
-\r
   Int_t mult1 = fESD->GetNumberOfMuonTracks() ;\r
 \r
\r
     for (Int_t j = 0; j<mult1; j++) { \r
        AliESDMuonTrack* mu1 = new AliESDMuonTrack(*(fESD->GetMuonTrack(j)));\r
        Float_t zr1 = mu1->Charge();\r
@@ -269,12 +269,12 @@ void AliCFMuonResTask1::UserExec(Option_t *)
            Float_t pxr1 = mu1->Px();\r
            Float_t pyr1 = mu1->Py();\r
            Float_t pzr1 = mu1->Pz();\r
-           Float_t phir1 = mu1->Phi(); // Warning in AliESDMuonTrack Phi = ATan2(Py,Px) = [-pi,pi]\r
-           phir1 = phir1 * 57.296;\r
+           Float_t phir1 = mu1->Phi(); \r
+           phir1=Phideg(phir1);\r
            Float_t er1 = mu1->E();\r
            Float_t rap1=Rap(er1,pzr1);\r
            // selection of the rapidity for first muon\r
-           if (!fCFManager->CheckParticleCuts(AliCFManager::kPartRecCuts,mu1)) continue;\r
+           if (!fCFManager->CheckParticleCuts(AliCFManager::kPartAccCuts,mu1)) continue;\r
 \r
            for (Int_t jj = 0; jj<mult1; jj++) {\r
                AliESDMuonTrack* mu2 = new AliESDMuonTrack(*(fESD->GetMuonTrack(jj)));\r
@@ -285,11 +285,11 @@ void AliCFMuonResTask1::UserExec(Option_t *)
                    Float_t pyr2 = mu2->Py();\r
                    Float_t pzr2 = mu2->Pz();\r
                    Float_t phir2 = mu2->Phi();\r
-                   phir2 = phir2 * 57.296;\r
+                   phir2 = Phideg(phir2);\r
                    Float_t er2 = mu2->E();\r
                    Float_t rap2=Rap(er2,pzr2);\r
                    // selection of the rapidity for second muon\r
-                   if (!fCFManager->CheckParticleCuts(AliCFManager::kPartRecCuts,mu2)) continue;\r
+                   if (!fCFManager->CheckParticleCuts(AliCFManager::kPartAccCuts,mu2)) continue;\r
 \r
                    Float_t prec = TMath::Sqrt((pxr1+pxr2)*(pxr1+pxr2)+(pyr1+pyr2)*(pyr1+pyr2)+\r
                                                (pzr1+pzr2)*(pzr1+pzr2));\r
@@ -314,7 +314,7 @@ void AliCFMuonResTask1::UserExec(Option_t *)
            }      // second mu Loop\r
        }          // mu- Selection\r
     }        \r
-\r
\r
 //  ----------\r
   fHistEventsProcessed->Fill(0);\r
   PostData(1,fHistEventsProcessed) ;\r
@@ -346,7 +346,7 @@ const Float_t AliCFMuonResTask1::Rap(Float_t e, Float_t pz)
 //________________________________________________________________________\r
 const Float_t AliCFMuonResTask1::Phideg(Float_t phi) \r
 {\r
-// calculate Phi from TParticle in range [-180,180] \r
+// calculate Phi in range [-180,180] \r
     Float_t phideg;\r
     \r
        phideg = phi-TMath::Pi();\r