changes to run with plugin
authordsilverm <dsilverm@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 23 Sep 2010 14:58:07 +0000 (14:58 +0000)
committerdsilverm <dsilverm@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 23 Sep 2010 14:58:07 +0000 (14:58 +0000)
PWG4/totEt/AliAnalysisEtReconstructed.cxx
PWG4/totEt/macros/CreateAlienHandlerCaloEtSim.C
PWG4/totEt/macros/runCaloEt.C

index 318bf5a7cbc059e8293f009b5f0dbd78f6695646..767ffa363fe0c797a141b3b071e7eac2cbd2d84e 100644 (file)
@@ -55,12 +55,13 @@ Int_t AliAnalysisEtReconstructed::AnalyseEvent(AliVEvent* ev)
 
     Double_t protonMass = fPdgDB->GetParticle("proton")->Mass();
 
-    AliESDEvent *realEvent = dynamic_cast<AliESDEvent*>(ev);
     //for PID
     AliESDpid *pID = new AliESDpid();
-    pID->MakePID(realEvent);
-    TObjArray* list = fEsdtrackCutsTPC->GetAcceptedTracks(realEvent);
+    pID->MakePID(event);
+    TObjArray* list = fEsdtrackCutsTPC->GetAcceptedTracks(event);
     Int_t nGoodTracks = list->GetEntries();
+    // printf("nGoodTracks %d nCaloClusters %d\n", nGoodTracks, event->GetNumberOfCaloClusters());
+
     for (Int_t iTrack = 0; iTrack < nGoodTracks; iTrack++)
       {
        AliESDtrack *track = dynamic_cast<AliESDtrack*> (list->At(iTrack));
@@ -113,7 +114,7 @@ Int_t AliAnalysisEtReconstructed::AnalyseEvent(AliVEvent* ev)
         }
 
         Double_t et = track->E() * TMath::Sin(track->Theta()) + massPart;
-       // printf("Rec track: iTrack %03d eta %4.3f phi %4.3f nITSCl %d nTPCCl %d\n", iTrack, track->Eta(), track->Phi(), nItsClusters, nTPCClusters); // tmp/debug printout
+       //printf("Rec track: iTrack %03d eta %4.3f phi %4.3f nITSCl %d nTPCCl %d\n", iTrack, track->Eta(), track->Phi(), nItsClusters, nTPCClusters); // tmp/debug printout
 
         if (TMath::Abs(track->Eta()) < fCuts->GetCommonEtaCut() && CheckGoodVertex(track) && nItsClusters > fCuts->GetReconstructedNItsClustersCut() && nTPCClusters > fCuts->GetReconstructedNTpcClustersCut() )
         {
@@ -186,7 +187,7 @@ Int_t AliAnalysisEtReconstructed::AnalyseEvent(AliVEvent* ev)
 
        if (cluster->GetType() != fClusterType) continue;
        //if(cluster->GetTracksMatched() > 0) 
-       //printf("Rec Cluster: iCluster %03d E %4.3f type %d NCells %d, nmatched: %d, distance to closest: %f\n", iCluster, cluster->E(), (int)(cluster->GetType()), cluster->GetNCells(), cluster->GetNTracksMatched(), cluster->GetEmcCpvDistance()); // tmp/debug printout
+       // printf("Rec Cluster: iCluster %03d E %4.3f type %d NCells %d, nmatched: %d, distance to closest: %f\n", iCluster, cluster->E(), (int)(cluster->GetType()), cluster->GetNCells(), cluster->GetNTracksMatched(), cluster->GetEmcCpvDistance()); // tmp/debug printout
               
         
         if (cluster->E() < fClusterEnergyCut) continue;
@@ -415,7 +416,7 @@ AliAnalysisEtReconstructed::CalcTrackClusterDistance(const Float_t clsPos[3],
     }
   } // iTrack
 
-  //  printf("CalcTrackClusterDistance: bestTrkMatch %d origTrkMatch %d distance %f\n", bestTrkMatchId, *trkMatchId, distance);
+  // printf("CalcTrackClusterDistance: bestTrkMatch %d origTrkMatch %d distance %f\n", bestTrkMatchId, *trkMatchId, distance);
   *trkMatchId = bestTrkMatchId;
   return distance;
 }
index 36c5717bf661e468f862f1820b54e821303f3a5a..50d859339b088b7c5859a5aa464a9ad2029cb891 100644 (file)
@@ -15,7 +15,7 @@ AliAnalysisGrid* CreateAlienHandlerCaloEtSim(TString outputDir, TString outputNa
 
   // Set versions of used packages
   plugin->SetAPIVersion("V1.1x");
-  plugin->SetROOTVersion("v5-27-05");
+  plugin->SetROOTVersion("v5-27-05-build1");
   plugin->SetAliROOTVersion("v4-20-09-AN-build1");
   // Declare input data to be processed.
 
@@ -36,7 +36,7 @@ AliAnalysisGrid* CreateAlienHandlerCaloEtSim(TString outputDir, TString outputNa
   // XML collections added via this method can be combined with the first method if
   // the content is compatible (using or not tags)
   plugin->AddDataFile("tag.xml");
-  //  plugin->AddDataFile("wn.xml"); // test
+  // plugin->AddDataFile("wn.xml"); // test
   // file generated with:  find -x tag /alice/sim/LHC10d1/117222/* AliESDs.root > tag.xml
 
   // Define alien work directory where all files will be copied. Relative to alien $HOME.
@@ -49,6 +49,9 @@ AliAnalysisGrid* CreateAlienHandlerCaloEtSim(TString outputDir, TString outputNa
   // Declare all libraries (other than the default ones for the framework. These will be
   // loaded by the generated analysis macro. Add all extra files (task .cxx/.h) here.
   plugin->SetAdditionalLibs("AliAnalysisEtCuts.h AliAnalysisEt.h AliAnalysisEtMonteCarlo.h AliAnalysisEtMonteCarloPhos.h AliAnalysisEtMonteCarloEmcal.h AliAnalysisEtReconstructed.h AliAnalysisEtReconstructedPhos.h AliAnalysisEtReconstructedEmcal.h AliAnalysisTaskTotEt.h AliAnalysisEtCuts.cxx AliAnalysisEt.cxx AliAnalysisEtMonteCarlo.cxx AliAnalysisEtMonteCarloPhos.cxx AliAnalysisEtMonteCarloEmcal.cxx AliAnalysisEtReconstructed.cxx AliAnalysisEtReconstructedPhos.cxx AliAnalysisEtReconstructedEmcal.cxx AliAnalysisTaskTotEt.cxx");
+  plugin->SetExecutableCommand("aliroot -b -q");
+  // add extra include files/path
+  plugin->AddIncludePath("-I. -I$ALICE_ROOT/EMCAL -I$ALICE_ROOT/ANALYSIS");
 
   // No need for output file names. Procedure is automatic. <-- not true
   plugin->SetDefaultOutputs(kFALSE);
@@ -74,5 +77,6 @@ AliAnalysisGrid* CreateAlienHandlerCaloEtSim(TString outputDir, TString outputNa
   plugin->SetPrice(1); 
   // Optionally modify split mode (default 'se')    
   plugin->SetSplitMode("se");
+
   return plugin;
 } 
index 9c78553b8d7d0359550aa9c8156daa14b3860f88..b88a21a6adbc3e6aa1976a93fbe30ec2972dd0c6 100644 (file)
@@ -19,7 +19,8 @@ void runCaloEt(bool submit = false, // true or false
   gSystem->Load("libMinuit");
 
   gSystem->AddIncludePath("-I$ALICE_ROOT/include");
-  
+  gSystem->AddIncludePath("-I. -I$ALICE_ROOT/EMCAL -I$ALICE_ROOT/ANALYSIS");
+
   gSystem->Load("libSTEERBase");
   gSystem->Load("libESD");
   gSystem->Load("libAOD");