--- PamelaLevel2/src/PamLevel2.cpp 2006/12/11 17:49:03 1.11 +++ PamelaLevel2/src/PamLevel2.cpp 2007/03/21 09:07:32 1.33 @@ -7,9 +7,10 @@ * Default constructor */ PamTrack::PamTrack(){ - trk_track = this->TrkTrack::GetTrkTrack(); - calo_track = this->CaloTrkVar::GetCaloTrkVar(); - tof_track = this->ToFTrkVar::GetToFTrkVar(); + trk_track = 0; + calo_track = 0; + tof_track = 0; + candeleteobj = 0; }; //-------------------------------------- // @@ -19,14 +20,66 @@ * Constructor */ PamTrack::PamTrack(TrkTrack* t, CaloTrkVar* c, ToFTrkVar* o){ - trk_track = this->TrkTrack::GetTrkTrack(); - calo_track = this->CaloTrkVar::GetCaloTrkVar(); - tof_track = this->ToFTrkVar::GetToFTrkVar(); - if(t) *trk_track = *t; - if(c) *calo_track = *c; - if(o) *tof_track = *o; + + trk_track = 0; + calo_track = 0; + tof_track = 0; +// if(t)trk_track = new TrkTrack(*t); +// if(c)calo_track = new CaloTrkVar(*c); +// if(o)tof_track = new ToFTrkVar(*o); + if(t)trk_track = t; + if(c)calo_track = c; + if(o)tof_track = o; + candeleteobj = 0; + }; +PamTrack::PamTrack(const PamTrack& track){ + + TrkTrack *t = track.trk_track; + CaloTrkVar *c = track.calo_track; + ToFTrkVar *o = track.tof_track; + + trk_track = 0; + calo_track = 0; + tof_track = 0; + if(t)trk_track = new TrkTrack(*t); + if(c)calo_track = new CaloTrkVar(*c); + if(o)tof_track = new ToFTrkVar(*o); + candeleteobj = 1; + +} +void PamTrack::Clear(){ +// cout << "PamTrack::Clear() "<TrkTrack::Clear(); + if(calo_track) calo_track->CaloTrkVar::Clear();//??? + if(tof_track) tof_track->ToFTrkVar::Clear();//??? + }else{ + trk_track = 0; + calo_track = 0; + tof_track = 0; + } +} +void PamTrack::Delete(){ +// cout << "PamTrack::Delete() "<TrkTrack::Clear(); + delete trk_track; + } + if(calo_track){ + calo_track->CaloTrkVar::Clear();//??? + delete calo_track; + } + if(tof_track){ + tof_track->ToFTrkVar::Clear();//??? + delete tof_track; + } + }else{ + Clear(); + } +} //-------------------------------------- // // @@ -36,95 +89,184 @@ */ PamLevel2::PamLevel2(){ - trk2_obj = TrkLevel2::GetTrkLevel2(); - trk1_obj = TrkLevel1::GetTrkLevel1(); - trkh_obj = TrkHough::GetTrkHough(); - calo1_obj = CaloLevel1::GetCaloLevel1(); - calo2_obj = CaloLevel2::GetCaloLevel2(); - tof_obj = ToFLevel2::GetToFLevel2(); - trig_obj = TrigLevel2::GetTrigLevel2(); - s4_obj = S4Level2::GetS4Level2(); - nd_obj = NDLevel2::GetNDLevel2(); - ac_obj = AcLevel2::GetAcLevel2(); - orb_obj = OrbitalInfo::GetOrbitalInfo(); - - run_obj = new GL_RUN(); - -// trk1_clone = 0; -// trk2_clone = 0; -// trkh_clone = 0; -// calo1_clone = 0; -// calo2_clone = 0; -// tof_clone = 0; -// trig_clone = 0; -// s4_clone = 0; -// nd_clone = 0; -// ac_clone = 0; -// orb_clone = 0; - -// T_clone = 0; -// C_clone = 0; -// O_clone = 0; -// R_clone = 0; -// S_clone = 0; -// N_clone = 0; -// A_clone = 0; -// B_clone = 0; +// trk2_obj = TrkLevel2::GetTrkLevel2(); +// trk1_obj = TrkLevel1::GetTrkLevel1(); +// trkh_obj = TrkHough::GetTrkHough(); +// calo1_obj = CaloLevel1::GetCaloLevel1(); +// calo2_obj = CaloLevel2::GetCaloLevel2(); +// tof_obj = ToFLevel2::GetToFLevel2(); +// trig_obj = TrigLevel2::GetTrigLevel2(); +// s4_obj = S4Level2::GetS4Level2(); +// nd_obj = NDLevel2::GetNDLevel2(); +// ac_obj = AcLevel2::GetAcLevel2(); +// orb_obj = OrbitalInfo::GetOrbitalInfo(); + + h0_obj = 0; + trk0_obj = 0; + calo0_obj = 0; + + trk2_obj = 0; + trk1_obj = 0; + trkh_obj = 0; + calo1_obj = 0; + calo2_obj = 0; + tof_obj = 0; + trig_obj = 0; + s4_obj = 0; + nd_obj = 0; + ac_obj = 0; + orb_obj = 0; + + run_obj = 0;//new GL_RUN(); + soft_obj = 0;// Emiliano + irun = -1LL; + runfirstentry = 0LL; + runlastentry = 0LL; + + l0_file = NULL; + l0_tree = NULL; + iroot = -1; + dbc = 0; + + prevshift = 0; + + run_tree = NULL; + run_tree_clone = NULL; + sel_tree = NULL; + sel_tree_clone = NULL; + + irunentry = -1LL; + pam_tree = NULL; + for(Int_t i=0; iGetenv("PAM_DBHOST"); + const char *pamdbuser=gSystem->Getenv("PAM_DBUSER"); + const char *pamdbpsw=gSystem->Getenv("PAM_DBPSW"); + if ( !pamdbhost ) pamdbhost = ""; + if ( !pamdbuser ) pamdbuser = ""; + if ( !pamdbpsw ) pamdbpsw = ""; + if ( strcmp(pamdbhost,"") ) host = pamdbhost; + if ( strcmp(pamdbuser,"") ) user = pamdbuser; + if ( strcmp(pamdbpsw,"") ) psw = pamdbpsw; - for(Int_t i=0; i<8; i++ )tree_clone[i]=NULL; + - sorted_tracks = new TRefArray(); +// sorted_tracks = 0;//new TRefArray(); + CAL0 = false; CAL1 = true; CAL2 = true; TRK2 = true; TRK1 = false; TRKh = false; - TRG = true; - TOF = true; - S4 = true; - ND = true; - AC = true; - ORB = true; + TRKh = false; + TRG = true; + TOF = true; + TOF0 = false; + S4 = true; + ND = true; + AC = true; + ORB = true; + RUN = true; + + SELLI = -1; + + tsorted=0; + timage=0; }; /** - * Destructor + * Delete the event (act as Dtor) */ -PamLevel2::~PamLevel2(){ - delete run_obj; - Clear(); +void PamLevel2::Delete(){ + + if(run_obj)delete run_obj; + if(soft_obj)delete soft_obj; //Emiliano + +// cout << "void PamLevel2::Clear()"<Delete(); + delete tsorted; + } + if(timage){ + timage->Delete(); + delete timage; + } + + if(dbc){ + dbc->Close(); + delete dbc; + } + + if(l0_file)l0_file->Close(); + if(pam_tree)pam_tree->Delete();; + if(run_tree)run_tree->Delete();; + if(sel_tree)sel_tree->Delete();; + for(Int_t i=0; iDelete();; + if(run_tree_clone)run_tree_clone->Delete();; + if(sel_tree_clone)sel_tree_clone->Delete();; + + }; + /** - * Clear the event + * Clear the event (NB! does not deallocate objects) */ void PamLevel2::Clear(){ - TrkLevel1::Clear(); - - TrkLevel2::Clear(); - CaloLevel2::Clear(); - ToFLevel2::Clear(); - TrigLevel2::Clear(); - S4Level2::Clear(); - NDLevel2::Clear(); - AcLevel2::Clear(); - OrbitalInfo::Clear(); - -// if(trk1_clone)trk1_clone->Clear(); -// if(trk2_clone)trk2_clone->Clear(); -// if(trkh_clone)trkh_clone->Clear(); -// if(calo1_clone)calo1_clone->Clear(); -// if(calo2_clone)calo2_clone->Clear(); -// if(tof_clone)tof_clone->Clear(); -// if(trig_clone)trig_clone->Clear(); -// if(s4_clone)s4_clone->Clear(); -// if(nd_clone)nd_clone->Clear(); -// if(ac_clone)ac_clone->Clear(); -// if(orb_clone)orb_clone->Clear(); +// cout << "void PamLevel2::Clear()"<Clear(); // Emiliano: Do not deallocate run_obj here, it will give segmentation fault! call clear instead +// if(soft_obj) soft_obj->Clear(); + + if(h0_obj) h0_obj->Clear(); +// if(trk0_obj) trk0_obj->Clear(); + if(calo0_obj) calo0_obj->Clear(); + if(trk1_obj) trk1_obj->Clear(); + if(trk2_obj) trk2_obj->Clear(); + if(trkh_obj) trkh_obj->Clear(); + if(calo1_obj)calo1_obj->Clear(); + if(calo2_obj)calo2_obj->Clear(); + if(tof_obj) tof_obj->Clear(); + if(trig_obj) trig_obj->Clear(); + if(s4_obj) s4_obj->Clear(); + if(nd_obj) nd_obj->Clear(); + if(ac_obj) ac_obj->Clear(); + if(orb_obj) orb_obj->Clear(); - sorted_tracks->Delete(); // clean the reference array +// if(sorted_tracks)sorted_tracks->Clear(); +// sorted_tracks.Clear(); + + if(tsorted){ + tsorted->Delete(); + } + if(timage){ + timage->Delete(); + } }; @@ -136,18 +278,68 @@ void *PamLevel2::GetPointerTo(const char* c ){ TString objname = c; - if(!objname.CompareTo("TrkLevel1"))return &trk1_obj; - if(!objname.CompareTo("TrkLevel2"))return &trk2_obj; - if(!objname.CompareTo("TrkHough"))return &trkh_obj; - if(!objname.CompareTo("CaloLevel1"))return &calo1_obj; - if(!objname.CompareTo("CaloLevel2"))return &calo2_obj; - if(!objname.CompareTo("ToFLevel2"))return &tof_obj; - if(!objname.CompareTo("TrigLevel2"))return &trig_obj; - if(!objname.CompareTo("S4Level2"))return &s4_obj; - if(!objname.CompareTo("NDLevel2"))return &nd_obj; - if(!objname.CompareTo("AcLevel2"))return &ac_obj; - if(!objname.CompareTo("OrbitalInfo"))return &orb_obj; + + if(!objname.CompareTo("TrkLevel1")) { + if(!trk1_obj){ + trk1_obj = new TrkLevel1(); + trk1_obj->Set(); + } + return &trk1_obj; + }; + if(!objname.CompareTo("TrkLevel2")) { + if(!trk2_obj){ + trk2_obj = new TrkLevel2(); + trk2_obj->Set(); + } + return &trk2_obj; + }; + if(!objname.CompareTo("TrkHough")) { + if(!trkh_obj) trkh_obj = new TrkHough(); + return &trkh_obj; + }; + if(!objname.CompareTo("CaloLevel1")) { + if(!calo1_obj) calo1_obj = new CaloLevel1(); + return &calo1_obj; + }; + if(!objname.CompareTo("CaloLevel2")) { + if(!calo2_obj){ + calo2_obj = new CaloLevel2(); + calo2_obj->Set(); + }; + return &calo2_obj; + }; + if(!objname.CompareTo("ToFLevel2")) { + if(!tof_obj){ + tof_obj = new ToFLevel2(); + tof_obj->Set(); + } + return &tof_obj; + }; + if(!objname.CompareTo("TrigLevel2")) { + if(!trig_obj) trig_obj = new TrigLevel2(); + return &trig_obj; + }; + if(!objname.CompareTo("S4Level2")) { + if(!s4_obj) s4_obj = new S4Level2(); + return &s4_obj; + }; + if(!objname.CompareTo("NDLevel2")) { + if(!nd_obj) nd_obj = new NDLevel2(); + return &nd_obj; + }; + if(!objname.CompareTo("AcLevel2")) { + if(!ac_obj) ac_obj = new AcLevel2(); + return &ac_obj; + }; + if(!objname.CompareTo("OrbitalInfo")){ + if(!orb_obj) orb_obj = new OrbitalInfo(); + return &orb_obj; + }; + if(!objname.CompareTo("RunInfo"))return &run_obj; + + if(!objname.CompareTo("SoftInfo"))return &soft_obj; // Emiliano + return NULL; }; //-------------------------------------- @@ -159,8 +351,10 @@ * (If seqno = -1 retrieves the self-trigger calorimeter track) */ CaloTrkVar *PamLevel2::GetCaloStoredTrack(int seqno){ - - if( CaloLevel2::ntrk()==0 ){ + + if( !calo2_obj )return 0; + + if( calo2_obj->CaloLevel2::ntrk()==0 ){ cout << "PamLevel2::GetCaloStoredTrack(int) : requested tracker SeqNo "<< seqno <<" but no Calorimeter tracks are stored"<CaloLevel2::GetCaloTrkVar(it_calo); it_calo++; - } while( c && seqno != c->trkseqno && it_calo < CaloLevel2::ntrk()); + } while( c && seqno != c->trkseqno && it_calo < calo2_obj->CaloLevel2::ntrk()); if(!c || seqno != c->trkseqno){ c = 0; @@ -190,7 +384,9 @@ */ ToFTrkVar *PamLevel2::GetToFStoredTrack(int seqno){ - if( ToFLevel2::ntrk()==0 ){ + if( !tof_obj )return 0; + + if( tof_obj->ToFLevel2::ntrk()==0 ){ cout << "PamLevel2::GetToFStoredTrack(int) : requested tracker SeqNo "<< seqno <<" but no ToF tracks are stored"<ToFLevel2::GetToFTrkVar(it_tof); it_tof++; - } while( c && seqno != c->trkseqno && it_tof < ToFLevel2::ntrk()); + } while( c && seqno != c->trkseqno && it_tof < tof_obj->ToFLevel2::ntrk()); if(!c || seqno != c->trkseqno){ c = 0; @@ -219,7 +415,10 @@ * Give the pamela track associated to a tracker track, retrieving related calorimeter and tof track information. */ PamTrack* PamLevel2::GetPamTrackAlong(TrkTrack* t){ - + + cout <<"PamLevel2::GetPamTrackAlong(TrkTrack* t) **obsolete** "<=0 && itrk < TrkLevel2::ntrk() ){ + if( itrk >=0 && itrk < trk2_obj->TrkLevel2::ntrk() ){ - TrkTrack *t = TrkLevel2::GetStoredTrack(itrk); + TrkTrack *t = trk2_obj->TrkLevel2::GetStoredTrack(itrk); track = GetPamTrackAlong(t); + + }else{ - cout << "PamLevel2::GetStoredTrack(int) : tracker track SeqNo "<< itrk <<" does not exist (GetNTracks() = "<Delete(); //temporaneo??? - + +// if(!sorted_tracks)sorted_tracks = new TRefArray(); +// sorted_tracks->Clear(); +// sorted_tracks.Clear(); + + if(!tsorted)tsorted = new TClonesArray("PamTrack",trk2_obj->GetNTracks()); + tsorted->Delete(); + TClonesArray &ttsorted = *tsorted; + if(!timage)timage = new TClonesArray("PamTrack",trk2_obj->GetNTracks()); + timage->Delete(); + TClonesArray &ttimage = *timage; + + // loop over the tracks sorted by the tracker Bool_t use_TRK = how.Contains("TRK", TString::kIgnoreCase); Bool_t use_CAL = how.Contains("CAL", TString::kIgnoreCase); @@ -298,22 +517,27 @@ // cout << "use_CAL "<TrkLevel2::GetNTracks(); i++){ - TrkTrack *ts = 0; + TrkTrack *ts = 0; + CaloTrkVar *cs = 0; + ToFTrkVar *os = 0; // get tracker tracks - TrkTrack *tp = TrkLevel2::GetTrack(i); //tracker - CaloTrkVar *cp = 0; - ToFTrkVar *op = 0; - + TrkTrack *tp = trk2_obj->TrkLevel2::GetTrack(i); //tracker + CaloTrkVar *cp = GetCaloStoredTrack(tp->GetSeqNo()); + ToFTrkVar *op = GetToFStoredTrack(tp->GetSeqNo()); + + TrkTrack *ti = 0; //tracker (image) + CaloTrkVar *ci = 0; + ToFTrkVar *oi = 0; // cout << "trk track n. "<HasImage()){ - TrkTrack *ti = TrkLevel2::GetTrackImage(i); //tracker (image) - CaloTrkVar *ci = 0; - ToFTrkVar *oi = 0; + ti = trk2_obj->TrkLevel2::GetTrackImage(i); //tracker (image) + ci = GetCaloStoredTrack(ti->GetSeqNo()); + oi = GetToFStoredTrack(ti->GetSeqNo()); // cout << "its image "< 15 && //no. of fit planes on Y view - varcfit[1] < 1000. && //fit variance on Y view + calo2_obj->npcfit[1] > 15 && //no. of fit planes on Y view + calo2_obj->varcfit[1] < 1000. && //fit variance on Y view + cp && ci && true){ - cp = GetCaloStoredTrack(tp->GetSeqNo()); - ci = GetCaloStoredTrack(ti->GetSeqNo()); - Float_t resy_p = cp->tbar[0][1] - cbar[0][1]; if(resy_p < 0)resy_p= - resy_p; - Float_t resy_i = ci->tbar[0][1] - cbar[0][1]; if(resy_i < 0)resy_i= - resy_i; + Float_t resy_p = cp->tbar[0][1] - calo2_obj->cbar[0][1]; if(resy_p < 0)resy_p= - resy_p; + Float_t resy_i = ci->tbar[0][1] - calo2_obj->cbar[0][1]; if(resy_i < 0)resy_i= - resy_i; if(resy_p <= resy_i) tp_score++; else ti_score++; @@ -349,13 +576,15 @@ // ------------------------ // check the number of hit pmts along the track // on S12 S21 and S32, where paddles are parallel to Y axis - if( use_TOF ){ + if( use_TOF && !tof_obj ){ + cout << "void PamLevel2::SortTracks(TString how): how= "<GetSeqNo()); - oi = GetToFStoredTrack(ti->GetSeqNo()); /* cout << "track: npmtadc "<< op->npmtadc << endl; cout << "track: npmttdc "<< op->npmttdc << endl; @@ -363,12 +592,12 @@ cout << "image: npmttdc "<< oi->npmttdc << endl;*/ for (Int_t ih=0; ih < op->npmtadc; ih++){ - Int_t pl = GetPlaneIndex( (op->pmtadc).At(ih) ); + Int_t pl = tof_obj->GetPlaneIndex( (op->pmtadc).At(ih) ); if(pl == 1 || pl == 2 || pl == 5)nphit_p++; }; for (Int_t ih=0; ih < oi->npmtadc; ih++){ - Int_t pl = GetPlaneIndex( (oi->pmtadc).At(ih) ); + Int_t pl = tof_obj->GetPlaneIndex( (oi->pmtadc).At(ih) ); if(pl == 1 || pl == 2 || pl == 5)nphit_i++; }; @@ -398,27 +627,59 @@ // ------------------------ // the winner is.... // ------------------------ - if (tp_score > ti_score) ts = tp;//the track sorted by the tracker!! - else if (tp_score < ti_score) ts = ti;//its image!! - else { - ts = tp; + if (tp_score > ti_score) { +// ts = tp;//the track sorted by the tracker!! +// cs = cp; +// os = op; + }else if (tp_score < ti_score) { + ts = ti;//its image!! + cs = ci; + os = oi; + + ti = tp;//its image!! + ci = cp; + oi = op; + + tp = ts;//its image!! + cp = cs; + op = os; + + + }else { +// ts = tp; +// cs = cp; +// os = op; // cout << "Warning - track image ambiguity not solved" << endl; // cout << ts->GetNtot() << " "<< ts->chi2 << " " << npcfit[1] << " "<< nphit_p << endl; }; }else{ - ts = tp; +// ts = tp; +// cs = cp; +// os = op; }; // cout <<" SortTracks() "<Add(ts);//save the track in the sorted array +// sorted_tracks->Add(ts);//save the track in the sorted array +// sorted_tracks.Add(ts);//save the track in the sorted array +// sorted_tracks.Add(tp);//save the track in the sorted array // cout << "SortTracks:: sorted_tracks->Add(it) "<GetEntries() != trk2_obj->GetNTracks() ){ + cout << "void PamLevel2::SortTracks(TString how): tsorted->GetEntries() "<GetEntries()<<" != trk2_obj->GetNTracks() = "<GetNTracks() <Delete(); tsorted=0; + timage->Delete(); timage=0; + } + //Restore Object count //To save space in the table keeping track of all referenced objects - //we assume that our events do not address each other. We reset the - //object count to what it was at the beginning of the event. + //We reset the object count to what it was at the beginning of the event. TProcessID::SetObjectCount(ObjectNumber); }; @@ -430,13 +691,24 @@ * This method overrides TrkLevel2::GetTracks(), where sorting is done by decreasing number of fit points and increasing chi^2. * PamLevel2::GetTracks() keeps the same track order given by TrkLevel2::GetTracks(), but checks image selection by using calorimeter and ToF tracking information. */ -TRefArray *PamLevel2::GetTracks(){ +// TRefArray *PamLevel2::GetTracks(){ + +// // *-*-*-*-*-*-*-*-*-*-*-*-* +// SortTracks("+CAL+TOF"); +// // *-*-*-*-*-*-*-*-*-*-*-*-* + +// // return sorted_tracks; +// return &sorted_tracks; + +// }; +TClonesArray *PamLevel2::GetTracks(){ // *-*-*-*-*-*-*-*-*-*-*-*-* SortTracks("+CAL+TOF"); // *-*-*-*-*-*-*-*-*-*-*-*-* - return sorted_tracks; + return tsorted; + }; //-------------------------------------- // @@ -448,18 +720,43 @@ * @param it Track number, ranging from 0 to GetNTracks(). */ PamTrack *PamLevel2::GetTrack(int it){ + +// if(!trk2_obj) return 0; +// // *-*-*-*-*-*-*-*-*-*-*-*-* +// SortTracks("+CAL+TOF"); +// // *-*-*-*-*-*-*-*-*-*-*-*-* +// // if(!sorted_tracks)return 0; +// if(sorted_tracks.IsEmpty())return 0; + +// PamTrack *track = 0; + +// if( it >=0 && it < trk2_obj->TrkLevel2::GetNTracks() && itGetEntries())return 0; + PamTrack *track = 0; - if( it >=0 && it < TrkLevel2::GetNTracks() && itGetEntriesFast() ){ - TrkTrack *t = (TrkTrack*)sorted_tracks->At(it); - track = GetPamTrackAlong(t); + if( it >=0 && it < trk2_obj->TrkLevel2::GetNTracks()){ +// TrkTrack *t = (TrkTrack*)sorted_tracks.At(it); +// track = GetPamTrackAlong(t); + TClonesArray &t = *(tsorted); + track = (PamTrack*)t[it]; + }else{ - cout << "PamLevel2::GetTrack(int) : tracker track SeqNo "<< it <<" does not exist (GetNTracks() = "<=0 && it < trk2_obj->TrkLevel2::GetNTracks() ){ +// TrkTrack *temp = (TrkTrack*)sorted_tracks.At(it); +// if( temp->HasImage() ){ +// TrkTrack *t = trk2_obj->TrkLevel2::GetStoredTrack(temp->GetImageSeqNo()); +// image = GetPamTrackAlong(t); +// }else{ +// cout <<"PamLevel2::GetTrackImage(int) : Track SeqNo "<GetEntries())return 0; PamTrack *image = 0; - if( it >=0 && it < TrkLevel2::GetNTracks() ){ - TrkTrack *temp = (TrkTrack*)sorted_tracks->At(it); - if( temp->HasImage() ){ - TrkTrack *t = TrkLevel2::GetStoredTrack(temp->GetImageSeqNo()); - image = GetPamTrackAlong(t); + if( it >=0 && it < trk2_obj->TrkLevel2::GetNTracks() ){ + TClonesArray &t = *(tsorted); + PamTrack *temp = (PamTrack*)t[it]; + if( temp->GetTrkTrack()->HasImage() ){ + TClonesArray &t = *(timage); + image = (PamTrack*)t[it]; }else{ - cout <<"PamLevel2::GetTrackImage(int) : Track SeqNo "<GetName(); if(!CheckLevel2File(fname))return NULL; - UInt_t *found=0; +// UInt_t *found=0; + + cout<< "GetPamTree(TFile*,TString): detector list --> "; + if(TRK1)cout<<"TRK1 "; + if(TRK2)cout<<"TRK2 "; + if(TRKh)cout<<"TRKH "; + if(CAL1)cout<<"CAL1 "; + if(CAL2)cout<<"CAL2 "; + if(TOF)cout<<"TOF "; + if(TRG)cout<<"TRG "; + if(AC)cout<<"AC "; + if(ND)cout<<"ND "; + if(S4)cout<<"S4 "; + if(ORB)cout<<"ORB "; + cout << endl; + if(SELLI)cout<<">>> Found selection-list <<<"<cd(); // Tracker TTree *T = (TTree*)f->Get("Tracker"); if(T && (TRK2||TRK1||TRKh)) { if(TRK2)T->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2")); - else T->SetBranchStatus("TrkLevel2",0,found); +// else T->SetBranchStatus("TrkLevel2",0,found); if(TRK2)cout << "Tracker : set branch address TrkLevel2"<SetBranchAddress("TrkLevel1", GetPointerTo("TrkLevel1")); - else T->SetBranchStatus("TrkLevel1",0,found); +// else T->SetBranchStatus("TrkLevel1",0,found); if(TRK1)cout << "Tracker : set branch address TrkLevel1"<SetBranchAddress("TrkHough", GetPointerTo("TrkHough")); - else T->SetBranchStatus("TrkHough",0,found); +// else T->SetBranchStatus("TrkHough",0,found); if(TRKh)cout << "Tracker : set branch address TrkHough"<AddFriend(T); @@ -540,16 +886,17 @@ TTree *C = (TTree*)f->Get("Calorimeter"); if(C && (CAL2||CAL1)) { if(CAL2)C->SetBranchAddress("CaloLevel2", GetPointerTo("CaloLevel2")); - else C->SetBranchStatus("CaloLevel2",0,found); +// else C->SetBranchStatus("CaloLevel2",0,found); if(CAL2)cout << "Calorimeter : set branch address CaloLevel2"<SetBranchAddress("CaloLevel1", GetPointerTo("CaloLevel1")); - else C->SetBranchStatus("CaloLevel1",0,found); +// else C->SetBranchStatus("CaloLevel1",0,found); if(CAL1)cout << "Calorimeter : set branch address CaloLevel1"<AddFriend(C); }else{ cout << "Calorimeter : missing tree"<Get("ToF"); if(O && TOF) { @@ -610,8 +957,24 @@ }else{ cout << "OrbitalInfo : missing tree"<Get("SelectionList"); + if(L && SELLI==1) { +// L->SetBranchAddress("RunEntry",&irun); +// cout << "SelectionList: set branch address RunEntry"<SetBranchAddress("EventEntry",&irunentry); +// cout << "SelectionList: set branch address EventEntry"<AddFriend("SelectionList"); + cout << " TTree *PamLevel2::GetPamTree(TFile, TString) >>> SelectionList not implemente!!!!"<GetEntries()<SetOwner(); - char *fullpath; +// char *fullpath; +// const char *fullpath; // if no input file list is given: if ( flisttxt != "" ){ - if( !gSystem->IsFileInIncludePath(flisttxt,&fullpath) ){ - cout <<"File "<ConcatFileName(gSystem->DirName(flisttxt),gSystem->BaseName(flisttxt)); - flisttxt = fullpath; - +// if( !gSystem->IsFileInIncludePath(flisttxt,&fullpath) ){ +// if( !(fullpath = gSystem->FindFile("./",flisttxt)) ){ +// cout <<"File "<ConcatFileName(gSystem->DirName(flisttxt),gSystem->BaseName(flisttxt)); + if( !gSystem->ChangeDirectory(ddir) ){ cout << "Cannot change directory : "<ChangeDirectory(wdir); // back to the working directory + return 0; + } int line=0; while (1) { TString file; @@ -671,16 +1043,23 @@ } if(file.Contains("#"))file = file(0,file.First("#")); // cout <<"(2) " << file << endl; - if( gSystem->IsFileInIncludePath(file,&fullpath) ){ - +// if( gSystem->IsFileInIncludePath(file,&fullpath) ){ +// if( (fullpath = gSystem->FindFile(ddir,file)) ){ + if( file.EndsWith(".root") ){ + char *fullpath = gSystem->ConcatFileName(gSystem->DirName(ddir),gSystem->BaseName(file)); contents->Add(new TSystemFile(fullpath,gSystem->DirName(fullpath)));// add file to the list - - }else{ + delete fullpath; + } +// }else{ // if(file.Data()!="")cout << "File: "< missing "<< endl; - }; +// }; }; in.close(); - + } else { + char *fullpath = gSystem->ConcatFileName(gSystem->DirName(flisttxt),gSystem->BaseName(flisttxt)); + contents->Add(new TSystemFile(fullpath,gSystem->DirName(fullpath)));// add file to the list + delete fullpath; + }; }else{ cout << "No input file list given."< GetName(); if( name.EndsWith(".root") ){ - char *fullpath; - gSystem->IsFileInIncludePath(name,&fullpath); +// const char *fullpath = gSystem->FindFile(ddir,name); +// char *fullpath; +// gSystem->IsFileInIncludePath(name,&fullpath); + char *fullpath = gSystem->ConcatFileName(gSystem->DirName(ddir),gSystem->BaseName(name)); contents->Add(new TSystemFile(fullpath,gSystem->DirName(fullpath))); + delete fullpath; }; } delete temp; @@ -711,7 +1093,7 @@ gSystem->ChangeDirectory(wdir); // back to the working directory // cout << endl << "Selected files:" << endl; // contents->Print(); - cout << contents->GetEntries()<<" files selected\n"; + cout << contents->GetEntries()<<" files \n"; // cout << endl; // cout << "Working directory: "<< gSystem->WorkingDirectory()<< endl; return contents; @@ -730,9 +1112,23 @@ TChain *Trout =0; - if( !detlist.IsNull() )SetWhichTrees(detlist); - +// if( !detlist.IsNull() )SetWhichTrees(detlist); + SetWhichTrees(detlist); + cout<< "GetPamTree(TList*,TString): input detector list --> "; + if(TRK1)cout<<"TRK1 "; + if(TRK2)cout<<"TRK2 "; + if(TRKh)cout<<"TRKH "; + if(CAL1)cout<<"CAL1 "; + if(CAL2)cout<<"CAL2 "; + if(TOF)cout<<"TOF "; + if(TRG)cout<<"TRG "; + if(AC)cout<<"AC "; + if(ND)cout<<"ND "; + if(S4)cout<<"S4 "; + if(ORB)cout<<"ORB "; + cout << endl; + TChain *T = 0; TChain *C = 0; TChain *O = 0; @@ -741,6 +1137,8 @@ TChain *N = 0; TChain *A = 0; TChain *B = 0; + + TChain *L = 0; if(TRK2||TRK1||TRKh) T = new TChain("Tracker"); if(CAL2||CAL1) C = new TChain("Calorimeter"); @@ -750,13 +1148,15 @@ if(ND) N = new TChain("NeutronD"); if(AC) A = new TChain("Anticounter"); if(ORB) B = new TChain("OrbitalInfo"); + + L = new TChain("SelectionList"); // loop over files and create chains TIter next(fl); TSystemFile *questo = 0; while ( (questo = (TSystemFile*) next()) ) { TString name = questo->GetName(); -// cout << "File: "<< name << endl; + cout << "File: "<< name << endl; if( CheckLevel2File(name) ){ if(TRK2||TRK1||TRKh) T->Add(name); if(CAL1||CAL2) C->Add(name); @@ -766,22 +1166,23 @@ if(ND) N->Add(name); if(AC) A->Add(name); if(ORB) B->Add(name); + if(SELLI==1) L->Add(name); }; } cout << "done chain \n"; - UInt_t *found=0; +// UInt_t *found=0; // Tracker if(T && (TRK2||TRK1||TRKh)) { if(TRK2)T->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2")); - else T->SetBranchStatus("TrkLevel2",0,found); +// else T->SetBranchStatus("TrkLevel2",0,found); if(TRK2)cout << "Tracker : set branch address TrkLevel2"<SetBranchAddress("TrkLevel1", GetPointerTo("TrkLevel1")); - else T->SetBranchStatus("TrkLevel1",0,found); +// else T->SetBranchStatus("TrkLevel1",0,found); if(TRK1)cout << "Tracker : set branch address TrkLevel1"<SetBranchAddress("TrkHough", GetPointerTo("TrkHough")); - else T->SetBranchStatus("TrkHough",0,found); +// else T->SetBranchStatus("TrkHough",0,found); if(TRKh)cout << "Tracker : set branch address TrkHough"<AddFriend("Tracker"); @@ -791,10 +1192,10 @@ // Calorimeter if(C && (CAL2||CAL1)) { if(CAL2)C->SetBranchAddress("CaloLevel2", GetPointerTo("CaloLevel2")); - else C->SetBranchStatus("CaloLevel2",0,found); +// else C->SetBranchStatus("CaloLevel2",0,found); if(CAL2)cout << "Calorimeter : set branch address CaloLevel2"<SetBranchAddress("CaloLevel1", GetPointerTo("CaloLevel1")); - else C->SetBranchStatus("CaloLevel1",0,found); +// else C->SetBranchStatus("CaloLevel1",0,found); if(CAL1)cout << "Calorimeter : set branch address CaloLevel1"<AddFriend("Calorimeter"); @@ -855,55 +1256,26 @@ }else{ cout << "OrbitalInfo : missing tree"<>> Found selection-list <<<"<SetBranchAddress("RunEntry",&irun); + cout << "SelectionList: set branch address RunEntry"<SetBranchAddress("EventEntry",&irunentry); + cout << "SelectionList: set branch address EventEntry"<AddFriend("SelectionList"); + }else{ +// cout << "SelectionList : missing tree"<Delete(); + }; // cout<GetEntries()<AddFriend("Tracker"); -// } -// // Calorimeter -// if(CAL2||CAL1) { -// if(!Tout)Tout=C; -// else Tout->AddFriend("Calorimeter"); -// }; -// // ToF -// if(TOF) { -// if(!Tout)Tout=O; -// else Tout->AddFriend("ToF"); -// }; -// // Trigger -// if(TRG) { -// if(!Tout)Tout=R; -// else Tout->AddFriend("Trigger"); -// }; -// // S4 -// if(S4) { -// if(!Tout)Tout=S; -// else Tout->AddFriend("S4"); -// }; -// // Neutron Detector -// if(ND) { -// if(!Tout)Tout=N; -// else Tout->AddFriend("NeutronD"); -// }; -// // Anticounters -// if(AC) { -// if(!Tout)Tout=A; -// else Tout->AddFriend("Anticounter"); -// }; -// // OrbitalInfo -// if(ORB) { -// if(!Tout)Tout=B; -// else Tout->AddFriend("OrbitalInfo"); -// }; - -// cout<GetEntries()<GetEntries() )PamLevel2::SetBranchAddress( Tout ); - + pam_tree = Trout; + return Trout; } @@ -985,6 +1357,13 @@ t->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo")); cout << "OrbitalInfo : set branch address OrbitalInfo"<SetBranchAddress("RunEntry", &irun); + cout << "SelectionList: set branch address RunEntry"<SetBranchAddress("EventEntry", &irunentry); + cout << "SelectionList: set branch address EventEntry"<SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo")); cout << "OrbitalInfo : set branch address OrbitalInfo"<SetBranchAddress("RunEntry", &irun); + cout << "SelectionList: set branch address RunEntry"<SetBranchAddress("EventEntry", &irunentry); + cout << "SelectionList: set branch address EventEntry"<Add(name); }; } + + if(R->GetNtrees()){ - R->SetBranchAddress("RunInfo", GetPointerTo("RunInfo")); - cout << "Run : set branch address RunInfo"<SetBranchAddress("RunInfo", GetPointerTo("RunInfo")); + cout << "Run : set branch address RunInfo"<SetBranchAddress("SoftInfo", GetPointerTo("SoftInfo")); // Emiliano + cout << "Software : set branch address SoftInfo"<Get("Run"); - R->SetBranchAddress("RunInfo", GetPointerTo("RunInfo")); - cout << "Run : set branch address RunInfo"<SetBranchAddress("RunInfo", GetPointerTo("RunInfo")); + cout << "Run : set branch address RunInfo"<SetBranchAddress("SoftInfo", GetPointerTo("SoftInfo")); // Emiliano + cout << "Software : set branch address SoftInfo"<GetEntries() <= 0 ) return(false); + // + +// Int_t oldrun = irun; + Long64_t oldrun = irun; + // -------------------------------------- + // if it is a full file (not preselected) + // -------------------------------------- + if(SELLI==0){ + + if( !GetOrbitalInfo() ){ + cout << "Bool_t PamLevel2::UpdateRunInfo(TChain *run, ULong64_t iev) -- ERROR -- missing OrbitalInfo "<GetEntry(irun); + runfirstentry = 0LL; + runlastentry += (Long64_t)(this->GetRunInfo()->NEVENTS); + if ( (Long64_t)(this->GetRunInfo()->NEVENTS) > 0LL ) runlastentry -= 1LL; +// cout << " ))))) UPDATE RUN INFO ((((( @iev "<ID<<" irun "<absTime >= GetRunInfo()->RUNHEADER_TIME && GetOrbitalInfo()->absTime <= GetRunInfo()->RUNTRAILER_TIME) && irun < run->GetEntries()-1 ){ +// while ( iev > (runfirstentry+(ULong64_t)(this->GetRunInfo()->NEVENTS-1)) && irun < run->GetEntries() ){ + irun++; + run->GetEntry(irun); + runfirstentry = runlastentry; + if ( (Long64_t)(this->GetRunInfo()->NEVENTS) > 0LL ) runfirstentry += 1LL; + runlastentry += (Long64_t)(this->GetRunInfo()->NEVENTS); +// cout << " ))))) UPDATE RUN INFO ((((( @iev "<ID<<" irun "<GetRunInfo()->NEVENTS,(ULong64_t)(runfirstentry+(ULong64_t)(this->GetRunInfo()->NEVENTS))); + // printf(" abstime %u trailertime %u \n",GetOrbitalInfo()->absTime,GetRunInfo()->RUNTRAILER_TIME); + // printf(" IDRUN %u \n",GetRunInfo()->ID); + // +// prevshift = 0; + // + if ( irun == run->GetEntries()-1 && fromfirst && !(GetOrbitalInfo()->absTime >= GetRunInfo()->RUNHEADER_TIME && GetOrbitalInfo()->absTime <= GetRunInfo()->RUNTRAILER_TIME)){ + printf(" resetting irun (it should NOT happen!!!)\n"); + fromfirst = false; + irun = 0; + run->GetEntry(irun); + runfirstentry = 0ULL; + runlastentry += (Long64_t)(this->GetRunInfo()->NEVENTS); + if ( (Long64_t)(this->GetRunInfo()->NEVENTS) > 0LL ) runlastentry -= 1LL; + }; + // + }; + // + if ( !fromfirst && + !(GetOrbitalInfo()->absTime >= GetRunInfo()->RUNHEADER_TIME && + GetOrbitalInfo()->absTime <= GetRunInfo()->RUNTRAILER_TIME) + ) { + printf(" Something very wrong here: cannot find RUN containing absolute time %u \n",GetOrbitalInfo()->absTime); + return false; + } + // + if ( irun == oldrun || irun >= run->GetEntries() ) return(false); + // + // printf(" iev %llu irun %i nevents %u 1st %llu last %llu \n",iev,irun,this->GetRunInfo()->NEVENTS,(ULong64_t)runfirstentry,(ULong64_t)runlastentry); + // + prevshift = 0; + cout << " ))))) UPDATE RUN INFO ((((( @iev "<ID<<" irun "<GetEntry(iev); + if(irun != oldrun){ + run->GetEntry(irun); + cout << " ))))) UPDATE RUN INFO ((((( @iev "<ID<<" irun "<GetEntries() <= 0 ) { +// cout << " Bool_t PamLevel2::UpdateRunInfo(ULong64_t iev) -- ERROR -- run tree is empty"<GetEntry(irun); +// runfirstentry = iev; +// if( iev != 0){ +// cout << "Bool_t PamLevel2::UpdateRunInfo(TChain *run, ULong64_t iev) -- WARNING -- "<absTime >= GetRunInfo()->RUNHEADER_TIME && GetOrbitalInfo()->absTime <= GetRunInfo()->RUNTRAILER_TIME) && irun < run->GetEntries() ){ +// irun++; +// if( irun == run->GetEntries() ){ +// if(!fromfirst){ +// cout << "Bool_t PamLevel2::UpdateRunInfo(TChain *run, ULong64_t iev) -- ERROR -- event entry #"<GetEntry(irun); +// runfirstentry = iev; +// }; +// // +// if ( irun == oldrun ) return(false); +// // +// cout << " ))))) UPDATE RUN INFO ((((( @iev "<ID<<" irun "<GetEntry(iev); +// if(irun != oldrun){ +// run->GetEntry(irun); +// cout << " ))))) UPDATE RUN INFO ((((( @iev "<ID<<" irun "< "; - if(TRK1)cout<<"TRK1 "; - if(TRK2)cout<<"TRK2 "; - if(TRKh)cout<<"TRKH "; - if(CAL1)cout<<"CAL1 "; - if(CAL2)cout<<"CAL2 "; - if(TOF)cout<<"TOF "; - if(TRG)cout<<"TRG "; - if(AC)cout<<"AC "; - if(ND)cout<<"ND "; - if(S4)cout<<"S4 "; - if(ORB)cout<<"ORB "; - cout << endl; +// cout<< "Set detector list --> "; +// if(TRK1)cout<<"TRK1 "; +// if(TRK2)cout<<"TRK2 "; +// if(TRKh)cout<<"TRKH "; +// if(CAL1)cout<<"CAL1 "; +// if(CAL2)cout<<"CAL2 "; +// if(TOF)cout<<"TOF "; +// if(TRG)cout<<"TRG "; +// if(AC)cout<<"AC "; +// if(ND)cout<<"ND "; +// if(S4)cout<<"S4 "; +// if(ORB)cout<<"ORB "; +// cout << endl; }; @@ -1227,24 +1837,26 @@ * Set tree/branch detector flags from the content of a tree */ void PamLevel2::GetWhichTrees(TFile* f){ - + + + cout << "void PamLevel2::GetWhichTrees(TFile* f) --- WARNING!! --- ...potrebbe non funzionare "<GetName()<IsZombie() ){ @@ -1354,7 +1966,21 @@ }; - lk->Delete(); + delete lk; + +// cout<< "Get detector list from input file --> "; +// if(TRK1)cout<<"TRK1 "; +// if(TRK2)cout<<"TRK2 "; +// if(TRKh)cout<<"TRKH "; +// if(CAL1)cout<<"CAL1 "; +// if(CAL2)cout<<"CAL2 "; +// if(TOF)cout<<"TOF "; +// if(TRG)cout<<"TRG "; +// if(AC)cout<<"AC "; +// if(ND)cout<<"ND "; +// if(S4)cout<<"S4 "; +// if(ORB)cout<<"ORB "; +// cout << endl; return ; @@ -1386,7 +2012,7 @@ Bool_t RUN__ok = false; - + Bool_t SELLI__ok = false; cout << "Checking file: "<GetName() << endl; -// cout << " Get tree: " << f->Get(key->GetName())<GetName() << ""<GetClassName()<Get(key->GetName())<GetClassName()=="TTree" && nev_previous && nev != nev_previous ){ @@ -1413,9 +2040,23 @@ // return false; // }; + if( !strcmp(key->GetName(),"Run" ) )RUN__ok = true; //========================================================= + if( !strcmp(key->GetName(),"SelectionList" ) ){ + SELLI__ok = true; + if(SELLI==1){ + Int_t nevt = ((TTree*)f->Get("SelectionList"))->GetEntries(); + if( nev && nevt!=nev){ + cout << "File: "<< f->GetName() <<" discarded ---- SelectionList tree has "<GetName(),"Trigger" ) ){ TRG__ok = true; if(TRG){ @@ -1527,7 +2168,41 @@ }; }; + + if( SELLI==-1 )SELLI = (Int_t)SELLI__ok; + if( SELLI==0 && SELLI__ok ){ + cout << "File: "<< f->GetName() <<" discarded ---- found SelectionList (it is not a full-event file)" << endl; + return false; + } + if( SELLI==1 && !SELLI__ok ){ + cout << "File: "<< f->GetName() <<" discarded ---- SelectionList missing" << endl; + return false; + } +// cout << "SELLI "< "; +// if(TRK1__ok)cout<<"TRK1 "; +// if(TRK2__ok)cout<<"TRK2 "; +// if(TRKh__ok)cout<<"TRKH "; +// if(CAL1__ok)cout<<"CAL1 "; +// if(CAL2__ok)cout<<"CAL2 "; +// if(TOF__ok)cout<<"TOF "; +// if(TRG__ok)cout<<"TRG "; +// if(AC__ok)cout<<"AC "; +// if(ND__ok)cout<<"ND "; +// if(S4__ok)cout<<"S4 "; +// if(ORB__ok)cout<<"ORB "; +// cout << endl; + + + if(TRK2 && TRK1__ok)TRK1=1; +// ---------------------------------------------------------------------------- +// NOTA +// se c'e` il level1, lo devo necessarimente leggere. +// infatti (non ho capito perche`) i cluster vengono letti e allocati in memoria +// comunque, ma non vengono disallocati da PamLevel2::Clear() +// ---------------------------------------------------------------------------- if(!RUN__ok) { @@ -1572,21 +2247,33 @@ return false; }; - if(ND && !ND__ok){ - cout << "File: "<< f->GetName() <<" discarded ---- Missing ND tree"<< endl; - return false; - }; - if(TRG && !TRG__ok){ - cout << "File: "<< f->GetName() <<" discarded ---- Missing Trigger tree"<< endl; - return false; - }; - - - - + if(ND && !ND__ok){ + cout << "File: "<< f->GetName() <<" discarded ---- Missing ND tree"<< endl; + return false; + }; + if(TRG && !TRG__ok){ + cout << "File: "<< f->GetName() <<" discarded ---- Missing Trigger tree"<< endl; + return false; + }; + - lk->Delete(); +// lk->Delete(); +// delete lk; f->Close(); + +// cout<< "CheckLevel2File(TString): detector list --> "; +// if(TRK1)cout<<"TRK1 "; +// if(TRK2)cout<<"TRK2 "; +// if(TRKh)cout<<"TRKH "; +// if(CAL1)cout<<"CAL1 "; +// if(CAL2)cout<<"CAL2 "; +// if(TOF)cout<<"TOF "; +// if(TRG)cout<<"TRG "; +// if(AC)cout<<"AC "; +// if(ND)cout<<"ND "; +// if(S4)cout<<"S4 "; +// if(ORB)cout<<"ORB "; +// cout << endl; return true; @@ -1602,11 +2289,11 @@ cout << "Create clones of PAMELA trees "<GetTree()->CloneTree(0); - TString name = tree_clone[i]->GetName(); + pam_tree_clone[i] = fChain->GetTree()->CloneTree(0); + TString name = pam_tree_clone[i]->GetName(); name.Append("_clone"); -// tree_clone[i]->SetName(name.Data()); - cout << tree_clone[i]->GetName() <SetName(name.Data()); + cout << pam_tree_clone[i]->GetName() <GetListOfFriends(); @@ -1616,14 +2303,16 @@ while( (T_friend = (TFriendElement*)next()) ){ // cout<IsA()->GetName()<<" "<GetName()<GetTree() << dec<GetTree()->GetName()<< endl; - tree_clone[i] = T_friend->GetTree()->CloneTree(0); - tree_clone[i]->SetAutoSave(1000000); - name = tree_clone[i]->GetName(); + pam_tree_clone[i] = T_friend->GetTree()->CloneTree(0); + pam_tree_clone[i]->SetAutoSave(1000000); + name = pam_tree_clone[i]->GetName(); name.Append("_clone"); -// tree_clone[i]->SetName(name.Data()); - cout << tree_clone[i]->GetName() << endl; +// pam_tree_clone[i]->SetName(name.Data()); + cout << pam_tree_clone[i]->GetName() << endl; i++; } + + delete li; cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <Branch("RunInfo","GL_RUN",GetPointerTo("RunInfo")); + cout << "Run : branch RunInfo"<Branch("SoftInfo","SoftInfo",GetPointerTo("SoftInfo")); + cout << "Run : branch SoftInfo"<Branch("RunEntry",&irun,"runentry/I"); + sel_tree_clone->Branch("EventEntry",&irunentry,"eventry/I"); + + + Int_t i=0; if(TRK1||TRK2||TRKh){ - tree_clone[i] = new TTree("Tracker","PAMELA tracker level2 data "); + pam_tree_clone[i] = new TTree("Tracker","PAMELA tracker level2 data "); if(TRK1) { - tree_clone[i]->Branch("TrkLevel1","TrkLevel1", GetPointerTo("TrkLevel1")); - tree_clone[i]->BranchRef(); + pam_tree_clone[i]->Branch("TrkLevel1","TrkLevel1", GetPointerTo("TrkLevel1")); + pam_tree_clone[i]->BranchRef(); cout << "Tracker : branch TrkLevel1"<Branch("TrkLevel2", "TrkLevel2",GetPointerTo("TrkLevel2")); + pam_tree_clone[i]->Branch("TrkLevel2", "TrkLevel2",GetPointerTo("TrkLevel2")); cout << "Tracker : branch TrkLevel2"<Branch("TrkHough","TrkHough", GetPointerTo("TrkHough")); + pam_tree_clone[i]->Branch("TrkHough","TrkHough", GetPointerTo("TrkHough")); cout << "Tracker : branch TrkHough"<Branch("CaloLevel1", "CaloLevel1", GetPointerTo("CaloLevel1")); + pam_tree_clone[i]->Branch("CaloLevel1", "CaloLevel1", GetPointerTo("CaloLevel1")); cout << "Calorimeter : branch CaloLevel1"<Branch("CaloLevel2","CaloLevel2", GetPointerTo("CaloLevel2")); + pam_tree_clone[i]->Branch("CaloLevel2","CaloLevel2", GetPointerTo("CaloLevel2")); cout << "Calorimeter : branch CaloLevel2"<Branch("ToFLevel2","ToFLevel2", GetPointerTo("ToFLevel2")); + pam_tree_clone[i] = new TTree("ToF","PAMELA ToF level2 data "); + pam_tree_clone[i]->Branch("ToFLevel2","ToFLevel2", GetPointerTo("ToFLevel2")); cout << "ToF : branch ToFLevel2"<Branch("TrigLevel2","TrigLevel2", GetPointerTo("TrigLevel2")); + pam_tree_clone[i] = new TTree("Trigger","PAMELA trigger level2 data "); + pam_tree_clone[i]->Branch("TrigLevel2","TrigLevel2", GetPointerTo("TrigLevel2")); cout << "Trigger : branch TrigLevel2"<Branch("S4Level2","S4Level2", GetPointerTo("S4Level2")); + pam_tree_clone[i] = new TTree("S4","PAMELA S4 level2 data "); + pam_tree_clone[i]->Branch("S4Level2","S4Level2", GetPointerTo("S4Level2")); cout << "S4 : branch S4Level2"<Branch("NDLevel2","NDLevel2", GetPointerTo("NDLevel2")); + pam_tree_clone[i] = new TTree("NeutronD","PAMELA neutron detector level2 data "); + pam_tree_clone[i]->Branch("NDLevel2","NDLevel2", GetPointerTo("NDLevel2")); cout << "NeutronD : branch NDLevel2"<Branch("AcLevel2","AcLevel2", GetPointerTo("AcLevel2")); + pam_tree_clone[i] = new TTree("Anticounter","PAMELA anticounter detector level2 data "); + pam_tree_clone[i]->Branch("AcLevel2","AcLevel2", GetPointerTo("AcLevel2")); cout << "Anticounter : branch AcLevel2"<Branch("OrbitalInfo","OrbitalInfo", GetPointerTo("OrbitalInfo")); + pam_tree_clone[i] = new TTree("OrbitalInfo","PAMELA oribital info "); + pam_tree_clone[i]->Branch("OrbitalInfo","OrbitalInfo", GetPointerTo("OrbitalInfo")); cout << "OrbitalInfo : branch OrbitalInfo"<Branch("TrkLevel1","TrkLevel1", &trk1_clone); -// T_clone->BranchRef(); -// cout << "Tracker : branch TrkLevel1"<Branch("TrkLevel2", "TrkLevel2",&trk2_clone); -// cout << "Tracker : branch TrkLevel2"<Branch("TrkHough","TrkHough", &trkh_clone); -// cout << "Tracker : branch TrkHough"<Branch("CaloLevel1", "CaloLevel1", &calo1_clone); -// cout << "Calorimeter : branch CaloLevel1"<Branch("CaloLevel2","CaloLevel2", &calo2_clone); -// cout << "Calorimeter : branch CaloLevel2"<Branch("ToFLevel2","ToFLevel2", &tof_clone); -// cout << "ToF : branch ToFLevel2"<Branch("TrigLevel2","TrigLevel2", &trig_clone); -// cout << "Trigger : branch TrigLevel2"<Branch("S4Level2","S4Level2", &s4_clone); -// cout << "S4 : branch S4Level2"<Branch("NDLevel2","NDLevel2", &nd_clone); -// cout << "NeutronD : branch NDLevel2"<Branch("AcLevel2","AcLevel2", &ac_clone); -// cout << "Anticounter : branch AcLevel2"<Branch("OrbitalInfo","OrbitalInfo", &orb_clone); -// cout << "OrbitalInfo : branch OrbitalInfo"<Clear(); -// trk2_obj->Copy(*trk2_clone); -// // *trk2_clone = *trk2_obj; -// } -// if(trkh_clone) *trkh_clone = *trkh_obj; -// if(calo1_clone){ -// // *calo1_clone = *calo1_obj; -// calo1_clone->Clear(); -// calo1_obj->Copy(*calo1_clone); -// } -// if(calo2_clone){ -// // *calo2_clone = *calo2_obj; -// calo2_clone->Clear(); -// calo2_obj->Copy(*calo2_clone); -// } -// if(tof_clone) *tof_clone = *tof_obj; -// if(trig_clone) *trig_clone = *trig_obj; -// if(s4_clone) *s4_clone = *s4_obj; -// if(nd_clone) *nd_clone = *nd_obj; -// if(ac_clone) *ac_clone = *ac_obj; -// if(orb_clone) *orb_clone = *orb_obj; - -// if(T_clone)T_clone->Fill(); -// if(C_clone)C_clone->Fill(); -// if(O_clone)O_clone->Fill(); -// if(R_clone)R_clone->Fill(); -// if(S_clone)S_clone->Fill(); -// if(N_clone)N_clone->Fill(); -// if(A_clone)A_clone->Fill(); -// if(O_clone)O_clone->Fill(); + +// cout << "PamLevel2::FillCloneTrees()" << irunentry << endl; - for(Int_t i=0; i<8; i++){ - if(tree_clone[i])tree_clone[i]->Fill(); + for(Int_t i=0; iFill(); } + if(sel_tree_clone)sel_tree_clone->Fill(); + } TTree* PamLevel2::GetCloneTree(TString name){ - for(Int_t i=0; i<8; i++){ - if(tree_clone[i]){ - TString na = tree_clone[i]->GetName(); - if(!name.CompareTo(na))return tree_clone[i]; + for(Int_t i=0; iGetName(); + if(!name.CompareTo(na))return pam_tree_clone[i]; }; } + if(run_tree_clone){ + TString na = run_tree_clone->GetName(); + if(!name.CompareTo(na))return run_tree_clone; + } + if(sel_tree_clone){ + TString na = sel_tree_clone->GetName(); + if(!name.CompareTo(na))return sel_tree_clone; + } return NULL; } void PamLevel2::WriteCloneTrees(){ cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <GetName()<Write(); + cout << run_tree_clone->GetName()<Write(); + cout << sel_tree_clone->GetName()<Write(); + for(Int_t i=0; iGetName()<Write(); }; } cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <GetEntry(ii) ){ + cout << " Int_t PamLevel2::GetEntry(Int_t) -- ERROR -- error reading pam tree"<Fill(); + +// cout << "PamLevel2::GetEntry("<ID_ROOT_L0; +// cout << "iroot "<IsConnected())){ + cout<<"~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~"<IsConnected() )return NULL; +// cout<<"...done"<Close(); + l0_file->Delete(); + } + cout << "Opening LEVEL0 file: "<< filename << endl; + FileStat_t t; + if( gSystem->GetPathInfo(filename.Data(),t) ){ + cout << " PamLevel2::GetYodaTree() -- ERROR opening file "<Get("Physics"); + if(!h0_obj)h0_obj = new EventHeader(); + l0_tree->SetBranchAddress("Header" ,&h0_obj); + prevshift = 0; + //--------------------------------------------------- + // TRACKER: + if(TRK0){ + if(!trk0_obj){ + trk0_obj = new TrkLevel0(); + trk0_obj->Set(); + }; + l0_tree->SetBranchAddress("Tracker" ,trk0_obj->GetPointerToTrackerEvent()); + } + //--------------------------------------------------- + // CALORIMETER: + if(CAL0){ + if(!calo0_obj){ + calo0_obj = new CaloLevel0(); + calo0_obj->Set(); + }; + l0_tree->SetBranchAddress("Calorimeter" ,calo0_obj->GetPointerToCalorimeterEvent()); + // cout << "PamLevel2::GetYodaTree() --- level0 calorimeter not implemented "<IsConnected())){ + cout << " TTree* PamLevel2::GetYodaTree( ) -- no DB connected... hai fatto qualche cazzata "<EV_FROM); -// if(TRK1||TRK2||TRKh){ -// TTree *T = new TTree("Tracker_clone","PAMELA tracker level2 data "); -// if(TRK1) { -// trk1_clone = new TrkLevel1(); -// T->Branch("TrkLevel1","TrkLevel1", &trk1_clone); -// T->BranchRef(); -// cout << "Tracker : branch TrkLevel1"<Branch("TrkLevel2", "TrkLevel2",&trk2_clone); -// cout << "Tracker : branch TrkLevel2"<Branch("TrkHough","TrkHough", &trkh_clone); -// cout << "Tracker : branch TrkHough"<AddFriend("Tracker_clone"); -// } -// // Calorimeter -// if(CAL1||CAL2){ -// TTree *C = new TTree("Calorimeter_clone","PAMELA calorimeter level2 data "); -// if(CAL1) { -// calo1_clone = new CaloLevel1(); -// C->Branch("CaloLevel1", "CaloLevel1", &calo1_clone); -// cout << "Calorimeter : branch CaloLevel1"<Branch("CaloLevel2","CaloLevel2", &calo2_clone); -// cout << "Calorimeter : branch CaloLevel2"<AddFriend("Calorimeter_clone"); -// } - -// // ToF -// if(TOF) { -// TTree *O = new TTree("ToF_clone","PAMELA ToF level2 data "); -// tof_clone = new ToFLevel2(); -// O->Branch("ToFLevel2","ToFLevel2", &tof_clone); -// cout << "ToF : branch ToFLevel2"<AddFriend("ToF_clone"); -// }; -// // Trigger -// if(TRG) { -// TTree *R = new TTree("Trigger_clone","PAMELA trigger level2 data "); -// trig_clone = new TrigLevel2(); -// R->Branch("TrigLevel2","TrigLevel2", &trig_clone); -// cout << "Trigger : branch TrigLevel2"<AddFriend("Trigger_clone"); -// }; -// // S4 -// if(S4) { -// TTree *S = new TTree("S4_clone","PAMELA S4 level2 data "); -// s4_clone = new S4Level2(); -// S->Branch("S4Level2","S4Level2", &s4_clone); -// cout << "S4 : branch S4Level2"<AddFriend("S4_clone"); -// }; -// // Neutron Detector -// if(ND) { -// TTree *N = new TTree("NeutronD_clone","PAMELA neutron detector level2 data "); -// nd_clone = new NDLevel2(); -// N->Branch("NDLevel2","NDLevel2", &nd_clone); -// cout << "NeutronD : branch NDLevel2"<AddFriend("NeutronD_clone"); -// }; -// // Anticounters -// if(AC) { -// TTree *A = new TTree("Anticounter_clone","PAMELA anticounter detector level2 data "); -// ac_clone = new AcLevel2(); -// A->Branch("AcLevel2","AcLevel2", &ac_clone); -// cout << "Anticounter : branch AcLevel2"<AddFriend("Anticounter_clone"); -// }; -// // OrbitalInfo -// if(ORB) { -// TTree *B = new TTree("OrbitalInfo_clone","PAMELA oribital info "); -// orb_clone = new OrbitalInfo(); -// B->Branch("OrbitalInfo","OrbitalInfo", &orb_clone); -// cout << "OrbitalInfo : branch OrbitalInfo"<AddFriend("OrbitalInfo_clone"); -// }; -// cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <EV_FROM "<EV_FROM <<" quella giusta "<ID_ROOT_L0 "<ID_ROOT_L0<absTime << endl; +// cout << " trk_calib_used "<TRK_CALIB_USED<< endl; + + if( !GetOrbitalInfo() ){ + cout << "Int_t PamLevel2::GetYodaEntry() -- ERROR -- missing OrbitalInfo "<0){ + cout << " level0 <--> level2 mismatch ( irun "<GetEntry(quellagiusta+(Long64_t)shift+(Long64_t)prevshift); + shift++; + if( !GetEventHeader() ){ + cout << "Int_t PamLevel2::GetYodaEntry() -- ERROR -- missing EventHeader "<pkt_num << " --- L0 --- "<< GetEventHeader()->GetPscuHeader()->GetCounter()<OBT << " --- L0 --- "<< GetEventHeader()->GetPscuHeader()->GetOrbitalTime()<GetEntries()+1 )cout << ">>> end of level0 tree <<<"<OBT "<< GetOrbitalInfo()->OBT << endl; +// cout << " GetEventHeader()->GetPscuHeader()->GetOrbitalTime() "<< GetEventHeader()->GetPscuHeader()->GetOrbitalTime() << endl; +// cout << " GetOrbitalInfo()->pkt_num "<< GetOrbitalInfo()->pkt_num << endl; +// cout << " GetEventHeader()->GetPscuHeader()->GetCounter() "<< GetEventHeader()->GetPscuHeader()->GetCounter() << endl; +// printf(" IDRUN %u \n",GetRunInfo()->ID); +// + if ( prevshift != 0 && (quellagiusta+(Long64_t)shift) == GetYodaTree()->GetEntries() ){ + prevshift = 0; + shift = -1; + }; -// tree_clone = Tout; -// tree_clone->SetDirectory(0); + }while( ( GetOrbitalInfo()->OBT != (UInt_t)(GetEventHeader()->GetPscuHeader()->GetOrbitalTime()) || GetOrbitalInfo()->pkt_num != (UInt_t)(GetEventHeader()->GetPscuHeader()->GetCounter())) && (quellagiusta+(Long64_t)shift) < GetYodaTree()->GetEntries() ); -// return Tout; -// } -// /** -// * Fill a tree (created with GetNewPamTree) -// * -// */ -// //void PamLevel2::FillNewPamTree(TTree *T){ -// void PamLevel2::FillNewPamTree(){ - + if ( (quellagiusta+(Long64_t)shift) == GetYodaTree()->GetEntries() ) cout << " Big trouble here, no such event in Level0 data! " <GetEntry(quellagiusta); + if ( shift > 1 ) prevshift += (shift-1); + + return answer; + +} -// if(trk1_clone) *trk1_clone = *trk1_obj; -// if(trk2_clone){ -// trk2_clone->Clear(); -// trk2_obj->Copy(*trk2_clone); -// // *trk2_clone = *trk2_obj; -// } -// if(trkh_clone) *trkh_clone = *trkh_obj; -// if(calo1_clone){ -// // *calo1_clone = *calo1_obj; -// calo1_clone->Clear(); -// calo1_obj->Copy(*calo1_clone); -// } -// if(calo2_clone){ -// // *calo2_clone = *calo2_obj; -// calo2_clone->Clear(); -// calo2_obj->Copy(*calo2_clone); -// } -// if(tof_clone) *tof_clone = *tof_obj; -// if(trig_clone) *trig_clone = *trig_obj; -// if(s4_clone) *s4_clone = *s4_obj; -// if(nd_clone) *nd_clone = *nd_obj; -// if(ac_clone) *ac_clone = *ac_obj; -// if(orb_clone) *orb_clone = *orb_obj; - -// TTree *T = tree_clone; - -// T->Fill(); //fill main tree -// // cout<IsA()->GetName()<<" "<GetName()<GetListOfFriends(); -// TIter next(li); -// TFriendElement* T_friend=0; -// while( (T_friend = (TFriendElement*)next()) ){ -// // cout<IsA()->GetName()<<" "<GetName()<GetTree() << dec<GetTree()->Fill();//fill friends -// } -// }