/[PAMELA software]/eventviewer/flight/src/FEVdetector.cpp
ViewVC logotype

Diff of /eventviewer/flight/src/FEVdetector.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.17 by mocchiut, Tue Mar 20 15:22:00 2007 UTC revision 1.18 by mocchiut, Fri Mar 30 11:16:43 2007 UTC
# Line 134  void FEVdetector::SetDDEC(TString de){ Line 134  void FEVdetector::SetDDEC(TString de){
134  };  };
135    
136  void FEVdetector::GetWindow(){  void FEVdetector::GetWindow(){
   printf(" ah \n");  
137    //    //
138    if ( level.file == 2 ){    if ( level.file == 2 ){
139      L2->GetEntry(maxevent);      L2->GetEntry(minevent);
140    } else {    } else {
141      otr->GetEntry(maxevent);      otr->GetEntry(minevent);
142    }    }
   //  
143    if ( level.file < 2 ){    if ( level.file < 2 ){
144      ph = eh->GetPscuHeader();      ph = eh->GetPscuHeader();
145      lastevno = (int)ph->Counter;      firstevno = (int)ph->Counter;
146    } else {    } else {
147      if ( var.ORB ){      if ( var.ORB ){
148        lastevno = L2->GetOrbitalInfo()->pkt_num;        firstevno = L2->GetOrbitalInfo()->pkt_num;
149      } else {      } else {
150        lastevno = 0;        firstevno = 0;
151      };      };
152    };    };
153    printf(" mmh minev %u \n",minevent);    //
154    if ( level.file == 2 ){    if ( level.file == 2 ){
155      L2->GetEntry(minevent);      L2->GetEntry(maxevent);
156    } else {    } else {
157      otr->GetEntry(minevent);      otr->GetEntry(maxevent);
158    }    }
159      //
160    if ( level.file < 2 ){    if ( level.file < 2 ){
161      ph = eh->GetPscuHeader();      ph = eh->GetPscuHeader();
162      firstevno = (int)ph->Counter;      lastevno = (int)ph->Counter;
163    } else {    } else {
164      if ( var.ORB ){      if ( var.ORB ){
165        firstevno = L2->GetOrbitalInfo()->pkt_num;        lastevno = L2->GetOrbitalInfo()->pkt_num;
166      } else {      } else {
167        firstevno = 0;        lastevno = 0;
168      };      };
169    };    };
170    printf(" aaah \n");    //
171  }  }
172    
173  void FEVdetector::ClearVariables(){  void FEVdetector::ClearVariables(){
# Line 193  void FEVdetector::ClearVariables(){ Line 192  void FEVdetector::ClearVariables(){
192  void FEVdetector::GetEntry(Int_t i){  void FEVdetector::GetEntry(Int_t i){
193    thisentry = i;      thisentry = i;  
194    if ( level.file == 2 ){    if ( level.file == 2 ){
195      printf("qui \n");      L2->Clear();
196      L2->GetEntry(i);      L2->GetEntry(i);
     printf("qua \n");  
197      ShowInfo("RunInfo");      ShowInfo("RunInfo");
198      ShowInfo("SoftInfo");      ShowInfo("SoftInfo");
199      ShowInfo("OrbitalInfo");      ShowInfo("OrbitalInfo");
# Line 346  void FEVdetector::ShowInfo(TString detec Line 344  void FEVdetector::ShowInfo(TString detec
344    //  TObjArray *leaf_array = 0;                    //  TObjArray *leaf_array = 0;                
345    //    //
346    if ( !strcmp(detector.Data(),"RunInfo") || !strcmp(detector.Data(),"SoftInfo") ){    if ( !strcmp(detector.Data(),"RunInfo") || !strcmp(detector.Data(),"SoftInfo") ){
347      if ( runtree ) b1 = runtree->FindBranch(detector.Data());      if ( L2->GetRunTree() ) b1 = L2->GetRunTree()->FindBranch(detector.Data());
348    } else {    } else {
349      b1 = otr->FindBranch(detector.Data());      b1 = otr->FindBranch(detector.Data());
350    };    };
# Line 855  void FEVdetector::GetGeneralInfo(){ Line 853  void FEVdetector::GetGeneralInfo(){
853      trcs = (TString)trc.str().c_str();      trcs = (TString)trc.str().c_str();
854    } else {    } else {
855      headcold = headc;      headcold = headc;
856        if ( var.RUN ){
857          var.idrun = L2->GetRunInfo()->ID;
858          var.idroot = L2->GetRunInfo()->ID_ROOT_L0;
859          var.EVF = L2->GetRunInfo()->EV_FROM;
860        };
861      if ( var.ORB ){      if ( var.ORB ){
862        headc = L2->GetOrbitalInfo()->pkt_num;        headc = L2->GetOrbitalInfo()->pkt_num;
863        OBT = L2->GetOrbitalInfo()->OBT;        OBT = L2->GetOrbitalInfo()->OBT;
# Line 1012  TChain* FEVdetector::Load(TFile &mainfil Line 1015  TChain* FEVdetector::Load(TFile &mainfil
1015    //    //
1016    file = &mainfile;    file = &mainfile;
1017    //    //
   if ( !L2 ) L2 = new PamLevel2();  
   //  
1018    if ( level.file == 2 ){    if ( level.file == 2 ){
1019      //      //
1020      printf(" %s \n",ddec.Data());      printf(" %s \n",ddec.Data());
1021      otr = L2->GetPamTree(gSystem->DirName(file->GetName()),file->GetName(),ddec.Data());      if ( !L2 ){
1022      runtree = L2->GetRunTree(gSystem->DirName(file->GetName()),file->GetName());        L2 = new PamLevel2("",file->GetName(),ddec.Data());
1023          otr = L2->GetPamTree();
1024        } else {
1025          //
1026          otr = L2->GetPamTree(gSystem->DirName(file->GetName()),file->GetName(),ddec.Data());
1027          L2->GetRunTree(gSystem->DirName(file->GetName()),file->GetName());
1028          //
1029        };
1030        //
1031      //      //
1032      if ( !otr ) return(0);      if ( !otr ) return(0);
1033      //      //
1034      otr->CanDeleteRefs(true);      //    otr->CanDeleteRefs(true);
1035      if ( runtree ) runtree->CanDeleteRefs(true);      //    if ( runtree ) runtree->CanDeleteRefs(true);
1036      //      //
1037    } else {    } else {
1038      //      //
# Line 5507  void FEVdetector::ShowTRK(Bool_t upd){ Line 5516  void FEVdetector::ShowTRK(Bool_t upd){
5516    
5517  void FEVdetector::ShowTrack(){  void FEVdetector::ShowTrack(){
5518    //    //
5519      if ( level.file != 2 ) return;
5520      //
5521    //    //
5522    //    //
5523    if ( var.CALO && var.CALOT ){    if ( var.CALO && var.CALOT ){
# Line 5604  void FEVdetector::ShowTrack(){ Line 5615  void FEVdetector::ShowTrack(){
5615      };      };
5616    };    };
5617    //    //
5618    if ( !var.TRK || level.file != 2 ) return;    if ( !var.TRK ) return;
5619    //    //
5620    if ( L2->GetTrkLevel2()->GetNTracks() > 0 ){    if ( L2->GetTrkLevel2()->GetNTracks() > 0 ){
5621      for (Int_t nt = 0; nt < L2->GetTrkLevel2()->GetNTracks(); nt++){      for (Int_t nt = 0; nt < L2->GetTrkLevel2()->GetNTracks(); nt++){
5622    
5623        //      TrkTrack *track = L2->GetTrkLevel2()->GetTrack(nt);        //      TrkTrack *track = L2->GetTrkLevel2()->GetTrack(nt);
5624          //      printf("qui \n");
5625        PamTrack *ptrack = L2->GetTrack(nt);        PamTrack *ptrack = L2->GetTrack(nt);
5626        TrkTrack *track = ptrack->GetTrkTrack();        TrkTrack *track = ptrack->GetTrkTrack();
5627          //      printf("qua \n");
5628        Int_t npoint = 100;        Int_t npoint = 100;
5629        Float_t zin[100];        Float_t zin[100];
5630        Double_t xout[100];        Double_t xout[100];
# Line 5819  void FEVdetector::ShowCalo(Bool_t upd){ Line 5832  void FEVdetector::ShowCalo(Bool_t upd){
5832          Int_t se = 5;          Int_t se = 5;
5833          if ( !l && m%2 == 0 ) se = 3;          if ( !l && m%2 == 0 ) se = 3;
5834          if ( !l && m%2 != 0 ) se = 2;          if ( !l && m%2 != 0 ) se = 2;
5835          if (  l && m%2 == 0 ) se = 1;          if (  l && m%2 != 0 ) se = 1;
5836          if (  l && m%2 != 0 ) se = 0;          if (  l && m%2 == 0 ) se = 0;
5837          //          //
5838          // determine what kind of event we are going to analyze          // determine what kind of event we are going to analyze
5839          //          //
# Line 6027  void FEVdetector::ShowCalo(Bool_t upd){ Line 6040  void FEVdetector::ShowCalo(Bool_t upd){
6040          Int_t se = 5;          Int_t se = 5;
6041          if (l == 0 && m%2 == 0) se = 3;          if (l == 0 && m%2 == 0) se = 3;
6042          if (l == 0 && m%2 != 0) se = 2;          if (l == 0 && m%2 != 0) se = 2;
6043          if (l == 1 && m%2 == 0) se = 1;          if (l == 1 && m%2 != 0) se = 1;
6044          if (l == 1 && m%2 != 0) se = 0;                  if (l == 1 && m%2 == 0) se = 0;        
6045          //          //
6046          // determine what kind of event we are going to analyze          // determine what kind of event we are going to analyze
6047          //          //
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.23