Babar/CM2 A-to-Z at Manchester

James Werner

Event reconstruction using experiment data.

recodata.cc


      1   // ------------------------------------------------------------------------
      2   //                          recodata.cc
      3   // Description:
      4   //    Class recodata - reconstruct parent particles from tracks and gammas
      5   //                       only experiment data
      6   //------------------------------------------------------------------------
      7   #include "BaBar/BaBar.hh"
      8
      9   //-----------------------
     10   // This Class's Header --
     11   //-----------------------
     12   #include "BetaMiniUser/recodata.hh"
     13
     14   //-------------
     15   // C Headers --
     16   //-------------
     17   #include <assert.h>
     18
     19   //---------------
     20   // C++ Headers --
     21   //---------------
     22   #include <iostream.h>
     23   #include <math.h>
     24
     25   //-------------------------------
     26   // Collaborating Class Headers --
     27   //-------------------------------
     28   #include "CLHEP/Alist/ConstAList.h" // produces a list
     29   #include "CLHEP/Alist/ConstAIterator.h" // tool to iterate throeugh a list
     30   #include "HepTuple/TupleManager.h"
     31   #include "HepTuple/Histogram.h"
     32   #include "HepTuple/Tuple.h"
     33   #include "AbsEvent/AbsEvent.hh"  //
     34   #include "AbsEvent/getTmpAList.hh" // for making lists
     35   #include "AbsEnv/AbsEnv.hh"
     36   #include "GenEnv/GenEnv.hh"
     37   #include "PDT/Pdt.hh"     // needed for particle data table
     38   #include "PDT/PdtPid.hh"  // the same
     39   #include "TrkBase/TrkRecoTrk.hh" // give you access to the track itself
     40                                    // There you will find DC information
     41   #include "TrkBase/TrkFit.hh"
     42   #include "PDT/PdtEntry.hh" //
     43   #include "BetaMicroAdapter/BtaCalQual.hh"  // EMC information
     44   #include "AbsParm/AbsParmIfdStrKey.hh"
     45   #include "Beta/EventInfo.hh"  // General event information
     46   #include "Beta/BtaCandidate.hh" // Candidates are what everything is about
     47   #include "BetaCoreTools/BtaMcAssoc.hh"
     48   #include "ErrLogger/ErrLog.hh"
     49
     50   //-----------------------------------------------------------------------
     51   // Local Macros, Typedefs, Structures, Unions and Forward Declarations --
     52   //-----------------------------------------------------------------------
     53
     54   static const char rcsid[] = "$Id: recodata.cc,v 1.7 2002/06/21 23:25:28 jtinslay Exp $";
     55
     56   //----------------
     57   // Constructors --
     58   //----------------
     59   recodata::recodata(
     60       const char* const theName,
     61       const char* const theDescription )
     62       : AppModule( theName, theDescription )
     63     , _eventInfoList("eventInfoList",this,"Default")       //we produce different
     64     , _btaPionList("pionCandidates",this,"piDefault")      //lists of candidates
     65     , _btaCalorList("calorCandidates",this,"CalorNeutral") // We fill them with
     66     , _btaTrackList("trackCandidates",this,"ChargedTracks")//default
     67     , _btaGamList( "btaGamList", this, "CalorNeutral")      // gamma
     68     , _btaTruthList("truthCandidates",this,"MCTruth")      //identification.
     69              // This is handeled in myAnalysis.tcl
     70              // You have to add the needed modules in AppUserBuild.cc too!
     71              // (I have deleted muons from the candidate lists. You may try
     72              // to make a muon list and loop over it)
     73   {
     74     commands()->append(& _eventInfoList);
     75     commands()->append(& _btaPionList);
     76     commands()->append(& _btaCalorList);
     77     commands()->append(& _btaTrackList);
     78     commands()->append(& _btaGamList);
     79     commands()->append(& _btaTruthList);
     80   }
     81
     82   //--------------
     83   // Destructor --
     84   //--------------
     85   recodata::~recodata( ) {}
     86
     87   //--------------
     88   // Operations --
     89   //--------------
     90
     91   // The begin(AppJob*) member function is run before any events are
     92   // processed.  In this analysis, it opens the output histogram file
     93   // and then books a number of histograms
     94
     95   AppResult
     96   recodata::beginJob( AbsEvent* anEvent )
     97   {
     98       _n_read_events = 0 ;
     99
    100       HepTupleManager* manager = gblEnv->getGen()->ntupleManager();
    101       assert(manager != 0);
    102
    103       _massaG2 = manager->histogram("Massa SD - comb 2 gamas",  12000, 0., 12. );
    104
    105       _massaC2 = manager->histogram("Massa - comb 2 Part",  12000, 0., 12. );
    106       _massaC3 = manager->histogram("Massa - comb 3 Part",  12000, 0., 12. );
    107       _massaC4 = manager->histogram("Massa - comb 4 Part",  12000, 0., 12. );
    108       _massaC5 = manager->histogram("Massa - comb 5 Part",  12000, 0., 12. );
    109       _massaC6 = manager->histogram("Massa - comb 6 Part",  12000, 0., 12. );
    110
    111       return AppResult::OK;
    112   }
    113
    114   void recodata::comb(int init, int numtermos, int itera, int tam, int *lista)
    115   {
    116   int i,j,xlista[200];
    117   float momx,momy,momz,ener,mas;
    118
    119   if((numtermos-itera)==1) {
    120     for(i=lista[numtermos-2]+1;i<tam;i++) {
    121       momx=0.0;
    122       momy=0.0;
    123       momz=0.0;
    124       ener=0.0;
    125       for(j=0;j<numtermos-1;j++)
    126         if(estado==0) {
    127           momx+=px[lista[j]];
    128           momy+=py[lista[j]];
    129           momz+=pz[lista[j]];
    130           ener+=energia[lista[j]];
    131         } else {
    132           momx+=cpx[lista[j]];
    133           momy+=cpy[lista[j]];
    134           momz+=cpz[lista[j]];
    135           ener+=cenergia[lista[j]];
    136         }
    137       if(estado==0) {
    138         momx+=px[i];
    139         momy+=py[i];
    140         momz+=pz[i];
    141         ener+=energia[i];
    142         mas=sqrt(ener*ener-(momx*momx+momy*momy+momz*momz));
    143         if((mas>0.125)&&(mas<0.155)) {
    144           cpx[contcomp]=momx;
    145           cpy[contcomp]=momy;
    146           cpz[contcomp]=momz;
    147           cenergia[contcomp]=ener;
    148         cmassa[contcomp]=mas;
    149           if(contcomp<199)
    150           contcomp++;
    151         }
    152       } else {
    153         momx+=cpx[i];
    154         momy+=cpy[i];
    155         momz+=cpz[i];
    156         ener+=cenergia[i];
    157         mas=sqrt(ener*ener-(momx*momx+momy*momy+momz*momz));
    158       }
    159       if(estado==0)
    160           _massaG2->accumulate(mas);
    161       else
    162         switch (numtermos) {
    163         case 2:
    164           _massaC2->accumulate(mas);
    165           break;
    166         case 3:
    167           _massaC3->accumulate(mas);
    168           break;
    169         case 4:
    170           _massaC4->accumulate(mas);
    171           break;
    172         case 5:
    173           _massaC5->accumulate(mas);
    174           break;
    175         case 6:
    176           _massaC6->accumulate(mas);
    177           break;
    178         }
    179     }
    180   } else {
    181     for(i=init;i<tam;i++) {
    182       for(j=0;j<itera;j++)
    183         xlista[j]=lista[j];
    184       xlista[itera]=i;
    185       comb(i+1,numtermos,itera+1,tam,xlista);
    186     }
    187   }
    188   }
    189
    190   AppResult
    191   recodata::event( AbsEvent* anEvent ) {
    192
    193     int i;
    194     _n_read_events++ ;
    195     HepPoint mcPriVtx = HepPoint(0.,0.,0.);
    196
    197     // Get the event summary info object from the EventInfo list
    198     HepAList<EventInfo>* eventList = Ifd<HepAList<EventInfo> >::get(anEvent,
    199                                                                   _eventInfoList.value());
    200     if (eventList == 0) {
    201       ErrMsg(fatal) << "Event list is null. Something is wrong in the job configuration."
    202                   << endmsg;
    203     }
    204
    205     // get list of charged tracks
    206     HepAList<BtaCandidate>* trackList = Ifd<HepAList<BtaCandidate> >::get(anEvent,
    207                                                                        _btaTrackList.value());
    208     if (trackList == 0) {
    209       ErrMsg(fatal) << "Charged track list is null. Something is wrong in the job configuration."
    210                   << endmsg;
    211     }
    212
    213     contpart=0;
    214     contcomp=0;
    215
    216     // photon list
    217
    218     HepAList<BtaCandidate>* gamList =
    219         Ifd<HepAList<BtaCandidate> >::get(anEvent, _btaGamList.value() );
    220
    221     HepAListIterator<BtaCandidate> gamIter(*gamList);
    222     BtaCandidate* theGamma;
    223     while (theGamma = gamIter() ) {
    224
    225       px[contpart]=theGamma->p4().x();
    226       py[contpart]=theGamma->p4().y();
    227       pz[contpart]=theGamma->p4().z();
    228       energia[contpart]=theGamma->energy();
    229       if(contpart<199)
    230         contpart++;
    231     } // end: while ( 0 != ( theGamma = gamIter() ) )
    232
    233     lista[0]=0;
    234     estado=0;
    235     if(contpart>1)
    236       comb(0,2,0,contpart,lista);
    237
    238   //  cout<<"Antes"<<endl;
    239   //  for(i=0;i<contpart;i++)
    240   //     cout<<i<<"  "<<cpx[i]<<"  "<<cpy[i]<<"  "<<cpz[i]<<"  "<<cenergia[i]<<"  "<<cmassa[i]<<endl;
    241   //  cout<<"Depois"<<endl;
    242
    243     // Loop over track candidates
    244
    245     HepConstAListIterator<BtaCandidate> itertrack(*trackList);
    246     const BtaCandidate* track_ptr;
    247
    248     while ( 0 != ( track_ptr = itertrack() ) ) {
    249
    250       cpx[contcomp]=track_ptr->p4().x();
    251       cpy[contcomp]=track_ptr->p4().y();
    252       cpz[contcomp]=track_ptr->p4().z();
    253       cenergia[contcomp]=track_ptr->energy();
    254       cmassa[contcomp]=sqrt(cenergia[contcomp]*cenergia[contcomp]-(cpx[contcomp]*cpx[contcomp]+
    255         cpy[contcomp]*cpy[contcomp]+cpz[contcomp]*cpz[contcomp]));
    256       if(contcomp<199)
    257         contcomp++;
    258     } // end: while ( 0 != ( track_ptr = itertrack() ) )
    259
    260   //  for(i=0;i<contcomp;i++)
    261   //      cout<<i<<"  "<<cpx[i]<<"  "<<cpy[i]<<"  "<<cpz[i]<<"  "<<cenergia[i]<<" "<<cmassa[i]<<endl;
    262     lista[0]=0;
    263     estado=1;
    264     for(i=2;(i<=4)&&(i<contcomp);i++)
    265       comb(0,i,0,contcomp,lista);
    266
    267     return AppResult::OK;
    268   }
    269
    270   AppResult
    271   recodata::endJob( AbsEvent* anEvent )
    272   {
    273        return AppResult::OK;
    274   }

[Download Source Code.]

recodata.hh


      1   //--------------------------------------------------------------------------
      2   //                          recodata.hh
      3   //
      4   // Description:
      5   //    Class recodata - reconstruct the parent particle from traks and gammas
      6   //                       Only data from experiment
      7   //------------------------------------------------------------------------
      8
      9   #ifndef RECODATA_HH
     10   #define RECODATA_HH
     11
     12   //----------------------
     13   // Base Class Headers --
     14   //----------------------
     15   #include "Framework/APPModule.hh"
     16   #include "Framework/AbsParmBool.hh"
     17   #include "AbsParm/AbsParmIfdStrKey.hh"
     18
     19   //------------------------------------
     20   // Collaborating Class Declarations --
     21   //------------------------------------
     22   class HepTupleManager;
     23   class HepHistogram;
     24   class HepTuple;
     25
     26   //            ---------------------
     27   //            -- Class Interface --
     28   //            ---------------------
     29
     30   class recodata : public AppModule {
     31
     32   //--------------------
     33   // Instance Members --
     34   //--------------------
     35
     36   public:
     37
     38       // Constructors
     39       recodata( const char* const theName, const char* const theDescription );
     40
     41       // Destructor
     42       virtual ~recodata( );
     43
     44       // Operations
     45
     46       virtual AppResult beginJob( AbsEvent* anEvent);
     47       virtual AppResult event( AbsEvent* anEvent );
     48       virtual AppResult endJob  ( AbsEvent* anEvent );
     49       virtual void comb(int init, int numtermos, int itera, int tam, int *lista);
     50
     51   protected:
     52
     53       AbsParmIfdStrKey _eventInfoList;
     54       AbsParmIfdStrKey _btaCalorList;
     55       AbsParmIfdStrKey _btaTrackList;
     56       AbsParmIfdStrKey _btaGamList;
     57       AbsParmIfdStrKey _btaTruthList;
     58       AbsParmIfdStrKey _btaPionList;
     59
     60   private:
     61     int  _n_read_events; // this is our event counter
     62     int contpart,contcomp,estado,lista[200];
     63     float px[200],py[200],pz[200],energia[200],massa[200];
     64     float cpx[200],cpy[200],cpz[200],cenergia[200],cmassa[200];
     65
     66     HepHistogram* _massaG2;
     67
     68     HepHistogram* _massaC2;
     69     HepHistogram* _massaC3;
     70     HepHistogram* _massaC4;
     71     HepHistogram* _massaC5;
     72     HepHistogram* _massaC6;
     73
     74   };
     75
     76   #endif

[Download Source Code.]

recodata.tcl


      1   #------------------------------------------------------------------------------
      2   #  MomEnH.tcl
      3   #------------------------------------------------------------------------------
      4   # always source the error logger early in your main tcl script
      5   sourceFoundFile ErrLogger/ErrLog.tcl
      6   sourceFoundFile FrameScripts/FwkCfgVar.tcl
      7   sourceFoundFile FrameScripts/talkto.tcl
      8   # Disable the use of envvars
      9   set ProdTclOnly true
     10
     11   # set the error logging level to 'warning'.  If you encounter a configuration
     12   # error you can get more information using 'trace'
     13   ErrLoggingLevel warning
     14
     15   #
     16   #  Turn off some specialty items
     17   #
     18   module disable NTrk
     19   module disable MomEnH
     20   module disable PEntp
     21   module disable colPEntp
     22   module disable EnPBip
     23   module disable LisEvt
     24   module disable verMCTruth
     25   module disable recoevt
     26
     27   ## allowed values of BetaMiniReadPersistence are (currently) "Kan", "Bdb"
     28   ##
     29   set BetaMiniReadPersistence Kan
     30
     31   ## allowed (non-expert) values of levelOfDetail are "micro", "cache", "extend"
     32   ## or "refit"
     33   ##
     34   FwkCfgVar levelOfDetail "cache"
     35
     36   ## allowed values of ConfigPatch are "Run1", "Run2" or "MC".  This MUST be set
     37   ## consistent ## with your input data type or you will get INCONSISTENT OR
     38   ## INCORRECT RESULTS
     39   ##
     40   FwkCfgVar ConfigPatch "Run2"
     41
     42   ##
     43   ##  You can enter input collections two ways: either append them to a list, or
     44   ##  explicitly enter them in the input module. Do one or the other, BUT NOT
     45   ##  BOTH.
     46   ##  If inputList is set before executing btaMini.tcl, that will automatically
     47   ##  add the collections to the appropriate input module, otherwise make sure you
     48   ##  talk to the right one.
     49   ##
     50   ## lappend inputList collection1 collection2 ...
     51   ##
     52   ##  OR THE FOLLOWING (choose the correct one based on persistence)
     53   ##
     54   ## talkto BdbEventInput {
     55   ## talkto KanEventInput {
     56   ##    input add collection1
     57   ##    input add collection2
     58   ##    ...
     59   ## }
     60
     61
     62   ##
     63   ## Set the number of events to run. If this isn't set, all events in the
     64   ## input collections will be processed.
     65   ##
     66   FwkCfgVar NEvent
     67
     68   ## choose the flavor of ntuple to write (hbook or root) and the file name
     69   ##
     70   FwkCfgVar BetaMiniTuple "root"
     71   FwkCfgVar histFileName "enphist.roo"
     72
     73   ## create Everything path and add core sequences to it. btaMiniPhysics is the
     74   ## same as btaMini, just appending a few standard list generating modules. For
     75   ## reading data with stored composites, you may have a conflict running
     76   ## btaMiniPhyscs.tcl
     77   ##
     78   ## You can also run (most of) the PhysProdSequence, complete with its 3 gamma
     79   ## conversion finders, etc. Consider disabling the portion of this sequence
     80   ## that you do not need to save yourself some time.  The BetaLumiSequence
     81   ## and TagProd sequences are left off, as they otherwise cause problems.
     82   ##
     83   sourceFoundFile BetaMiniUser/btaMini.tcl
     84   #sourceFoundFile BetaMiniUser/btaMiniPhysics.tcl
     85   #sourceFoundFile BetaMiniUser/btaMiniPhysProdSequence.tcl
     86   ## Add Analysis module
     87   ##
     88   path append Everything recodata
     89
     90   ##
     91   ##  If your job has a tag-level filter, here is how you should run it
     92   ##  so as to avoid wasting time reading the mini when the tag filter fails
     93   ##  Here's a simple example that restricts to just multi-hadron events
     94   ##  on Kan input
     95
     96   module clone TagFilterByName TagBGFMultiHadron
     97   module talk TagBGFMultiHadron
     98     andList set BGFMultiHadron
     99     assertIfMissing set true
    100   exit
    101
    102   module talk EvtCounter
    103     printFreq set 100000
    104   exit
    105
    106   #sequence append BetaMiniReadSequence -a KanEventUpdateTag TagBGFMultiHadron
    107
    108   path list
    109   if [info exists NEvent] {
    110     ev begin -nev $NEvent
    111   } else {
    112     ev begin
    113   }
    114
    115   ErrMsg trace "completed OK"
    116   exit
    117   

[Download Source Code.]


Gamma - Gamma reconstruction.

This example shows the use of grid and recursive reconstruction algorithm to obtain Pi0 from the combination of 2 gammas from EMC. Pi0 decays 98.78% of time in 2 photons. The dataset to obtain this histograms are Tau11-Run3-OnPeak-R14 with 80,000,000 events.

This is the overview of the energy histogram between 0 and 1 GeV showing Pi0 peak (134 Mev).

A closer look with energies between 0 and 200 MeV. There is a second peak around 50 MeV wich is background residuals due high quantity of low energy gammas.

This is the Pi0 peak, reconstructed from 2 gammas.

This is the region where eta (mass=547 MeV) should be. However, due the litle branch rate (0.00014% * 39%) it is impossible see the peak over the huge amount of background. The number of etas in this peak should be around 43 events.

This is the region where eta-prime should appear (mass=957 MeV). The number of events should be 0.000074% * 2% = 1 event !!!

Particle - Particle reconstruction.

This example shows the use of grid and recursive reconstruction algorithm to obtain the mass through 2 particles combination with Drift Chamber data. Is there a peak around 750 MeV ?

3 Particles reconstruction.

This example shows the use of grid and recursive reconstruction algorithm to obtain the mass through 3 particles combination with Drift Chamber data. At this point, the effect of background is so intense, that nothing can be seen.

Running data events reconstruction in grid

Program: recodata.cc/hh/tcl
         recodataJob.tcl

[jamwer@bfb gridrecodata]$ ./easygrid Tau11-Run3-OnPeak-R14
Searching pre selected skimdata.
skimData files not found ! Running BbkDatasetTcl.
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-1.tcl (3587123 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-2.tcl (3586623 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-3.tcl (3599537 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-4.tcl (3714806 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-5.tcl (3437617 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-6.tcl (3400520 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-7.tcl (3709093 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-8.tcl (3736506 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-9.tcl (3621644 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-10.tcl (3720442 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-11.tcl (3597749 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-12.tcl (3743655 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-13.tcl (3742591 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-14.tcl (3735212 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-15.tcl (3750987 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-16.tcl (3612365 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-17.tcl (3690077 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-18.tcl (3771847 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-19.tcl (3744435 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-20.tcl (1501758 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-21.tcl (1923806 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-22.tcl (564550 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-23.tcl (367724 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-24.tcl (3689414 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-25.tcl (2375655 events)
BbkDatasetTcl: wrote Tau11-Run3-OnPeak-R14-26.tcl (1162866 events)
Selected 26 collections, 81088602/444539210 events, ~30773.0/pb
Searching previous handlers.
Handlers not found. Submiting to GRID . Wait end of process...

bash-2.05b$ ssh -l root bfa.tier2.hep.man.ac.uk
root@bfa.tier2.hep.man.ac.uk's password:
Last login: Tue Nov 30 11:11:38 2004 from pc73.hep.man.ac.uk
[root@bfa root]# qstat
Job id           Name             User             Time Use S Queue
---------------- ---------------- ---------------- -------- - -----
12736.bfa        STDIN            babar003         00:18:23 R long
12737.bfa        STDIN            babar003         00:18:22 R long
12738.bfa        STDIN            babar003         00:18:21 R long
12739.bfa        STDIN            babar003         00:18:22 R long
12740.bfa        STDIN            babar003         00:18:22 R long
12741.bfa        STDIN            babar003         00:18:05 R long
12743.bfa        STDIN            babar003         00:16:50 R long
12744.bfa        STDIN            babar003         00:16:50 R long
12745.bfa        STDIN            babar003         00:16:50 R long
12746.bfa        STDIN            babar003         00:16:47 R long
12747.bfa        STDIN            babar003                0 Q long
12748.bfa        STDIN            babar003                0 Q long
12749.bfa        STDIN            babar003                0 Q long
12750.bfa        STDIN            babar003                0 Q long
12751.bfa        STDIN            babar003                0 Q long
12752.bfa        STDIN            babar003                0 Q long
12753.bfa        STDIN            babar003                0 Q long
12754.bfa        STDIN            babar003                0 Q long
12755.bfa        STDIN            babar003                0 Q long
12756.bfa        STDIN            babar003                0 Q long
12757.bfa        STDIN            babar003                0 Q long
12758.bfa        STDIN            babar003                0 Q long
12759.bfa        STDIN            babar003                0 Q long
12760.bfa        STDIN            babar003                0 Q long
12761.bfa        STDIN            babar003                0 Q long
12762.bfa        STDIN            babar003                0 Q long

Grid Results.

Date:1-7 December/2004 Number job submitted: 26 Jobs Number Successful results: 15 Jobs Number of Aborts: 11 42% 8 errors (computer locked for 24 hours each one of 3 attempts!): Cannot read JobWrapper output, both from Condor and from Maradona. Could it be the processing time and queue time limit ? 2,375,655 evts took 15 hours => 2,640 evts/minute One of failed programs had 3,587,123 evts => 22.6 hours < 24 hours limit ! 3 errors: Cannot download fullboot.sh from gsiftp://lcgrb01.gridpp.rl.ac.uk/var/edgwl/SandboxDir/Gj/https_3a_2f_2flxn1188.cern.ch_3a9000_2fGjjrlSsHnclaQKJ5CHkHdg/input/ Resubmission using bfa.tier2.hep.man.ac.uk:2119/jobmanager-lcgpbs-infinite queue: Date: 14-15 December/2004 Number of jobs: 11 Number of Successful results: 10 Jobs Number of aborts: 1 Cannot download fullboot.sh from ... REresubmission [jamwer@bfb gridrecodata]$ ./resub Tau11-Run3-OnPeak-R14 Job Resubmission. [jamwer@bfb gridrecodata]$ ./easygrid Tau11-Run3-OnPeak-R14 Searching pre selected skimdata. Searching previous handlers. Checking if jobs finished. ### Handle -> https://lxn1188.cern.ch:9000/z0rfVYJWxfddnyMrZf9Dgw Current Status: Ready https://lxn1188.cern.ch:9000/z0rfVYJWxfddnyMrZf9Dgw still pendent. 1 jobs did not finished ! Try again later. Submission OK! Selected Virtual Organisation name (from --vo option): babar Connecting to host lcgrb01.gridpp.rl.ac.uk, port 7772 Logging to host lcgrb01.gridpp.rl.ac.uk, port 9002 ********************************************************************************************* JOB SUBMIT OUTCOME The job has been successfully submitted to the Network Server. Use edg-job-status command to check job current status. Your job identifier (edg_jobId) is: - https://lxn1188.cern.ch:9000/z0rfVYJWxfddnyMrZf9Dgw ********************************************************************************************* [jamwer@bfb gridrecodata]$ edg-job-get-logging-info -v 1 https://lxn1188.cern.ch:9000/z0rfVYJWxfddnyMrZf9Dgw **** Error: API_NATIVE_ERROR **** Error while calling the "Eve::getEvents" native api Unable to retrieve the Job Events for: https://lxn1188.cern.ch:9000/z0rfVYJWxfddnyMrZf9Dgw edg_wll_JobLog: No such file or directory: no matching events found The job disapeared after 2 days!!! ReReReSubmission...
[Complete errors report]
Top

Last modified:
Copyright 2004 Manchester University
Feedback to: jamwer@hep.man.ac.uk