/[PAMELA software]/DarthVader/TrackerLevel2/src/TrkStruct.cpp
ViewVC logotype

Diff of /DarthVader/TrackerLevel2/src/TrkStruct.cpp

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

revision 1.5 by pam-fi, Thu Nov 30 17:05:09 2006 UTC revision 1.7 by pam-fi, Thu May 24 16:45:48 2007 UTC
# Line 4  Line 4 
4   * \date   * \date
5   */   */
6  #include <TrkStruct.h>  #include <TrkStruct.h>
7    //#include <TrkCalib.h>
8    
9  using namespace pamela;  using namespace pamela;
10  using namespace std;  using namespace std;
# Line 22  float anint(float f){ Line 22  float anint(float f){
22      ret1=(float)ret;      ret1=(float)ret;
23      return ret1;              return ret1;        
24  }  }
 /**  
  * Method to fill the tracker calibration-parameter struct from on-line calibrations  
  */  
 void cTrkCalib::FillACalibFrom(TFile* f0, Int_t ev_reg1, Int_t ev_reg2){  
   
     TTree *tr1 = (TTree*)f0->Get("CalibTrk1");    
     TTree *tr2 = (TTree*)f0->Get("CalibTrk2");  
   
     CalibTrk1Event *calibdata1 = 0;  
     CalibTrk2Event *calibdata2 = 0;  
     tr1->SetBranchAddress("CalibTrk1",&calibdata1);  
     tr2->SetBranchAddress("CalibTrk2",&calibdata2);  
     tr1->GetEntry(ev_reg1);  
     tr2->GetEntry(ev_reg2);  
     Int_t dspn1(0),dspn2(0);  
           
     for(Int_t i=0;i<6;i++){  
         dspn1=calibdata1->DSPnumber[i]-1;  
         dspn2=calibdata2->DSPnumber[i]-1;  
         /*     cout<<"dspn1= "<<dspn1<<"   dspn2= "<<dspn2<<endl; */  
         for(Int_t ii=0;ii<3072;ii++){  
             Int_t j=(ii/128);  
             Int_t jj=ii-j*128;  
             pedestal[jj][j][dspn1]=calibdata1->DSPped_par[i][ii];  
             pedestal[jj][j][dspn2]=calibdata2->DSPped_par[i][ii];  
             pedestal_t[jj][j][dspn1]=aint(calibdata1->DSPped_par[i][ii]);  
             pedestal_t[jj][j][dspn2]=aint(calibdata2->DSPped_par[i][ii]);  
             sigma[jj][j][dspn1]=calibdata1->DSPsig_par[i][ii];  
             sigma[jj][j][dspn2]=calibdata2->DSPsig_par[i][ii];  
             sigma_t[jj][j][dspn1]=anint(calibdata1->DSPsig_par[i][ii]);  
             sigma_t[jj][j][dspn2]=anint(calibdata2->DSPsig_par[i][ii]);  
             //        cout<<"sigma trk1= "<<sigma[jj][j][dspn1]<<endl;  
             //        cout<<"sigma_t trk1= "<<sigma_t[jj][j][dspn1]<<endl;  
             //        cout<<"sigma trk2= "<<sigma[jj][j][dspn2]<<endl;  
             //        cout<<"sigma_t trk2= "<<sigma_t[jj][j][dspn2]<<endl;  
           
             //  
             // invert the bad  
             //  
             if(calibdata1->DSPbad_par[i][ii]==0)  
                 bad[jj][j][dspn1]=1;  
             else if(calibdata1->DSPbad_par[i][ii]==1)  
                 bad[jj][j][dspn1]=0;  
             if(calibdata2->DSPbad_par[i][ii]==0)  
                 bad[jj][j][dspn2]=1;  
             else if(calibdata2->DSPbad_par[i][ii]==1)  
                 bad[jj][j][dspn2]=0;  
         };  
     };  
           
     tr1->Delete();  
     tr2->Delete();  
 };  
   
 /**  
  * Method to fill the tracker calibration-parameter struct from on-line calibrations  
  */  
 void cTrkCalib::FillFCalibFrom(TFile* f0, Int_t ev_reg1, Int_t ev_reg2){  
   
     TTree *tr1 = (TTree*)f0->Get("CalibTrk1");    
     TTree *tr2 = (TTree*)f0->Get("CalibTrk2");  
   
     CalibTrk1Event *calibdata1 = 0;  
     CalibTrk2Event *calibdata2 = 0;  
     tr1->SetBranchAddress("CalibTrk1",&calibdata1);  
     tr2->SetBranchAddress("CalibTrk2",&calibdata2);  
     tr1->GetEntry(ev_reg1);  
     tr2->GetEntry(ev_reg2);  
     Int_t dspn1(0),dspn2(0);  
           
     for(Int_t i=0;i<6;i++){  
         dspn1=calibdata1->DSPnumber[i]-1;  
         dspn2=calibdata2->DSPnumber[i]-1;  
         /*     cout<<"dspn1= "<<dspn1<<"   dspn2= "<<dspn2<<endl; */  
         for(Int_t ii=0;ii<3072;ii++){  
             Int_t j=(ii/128);  
             Int_t jj=ii-j*128;  
             pedestal[jj][j][dspn1]=calibdata1->DSPped_par[i][ii];  
             pedestal[jj][j][dspn2]=calibdata2->DSPped_par[i][ii];  
             sigma[jj][j][dspn1]=calibdata1->DSPsig_par[i][ii];  
             sigma[jj][j][dspn2]=calibdata2->DSPsig_par[i][ii];  
             //        cout<<"sigma trk1= "<<sigma[jj][j][dspn1]<<endl;  
             //        cout<<"sigma_t trk1= "<<sigma_t[jj][j][dspn1]<<endl;  
             //        cout<<"sigma trk2= "<<sigma[jj][j][dspn2]<<endl;  
             //        cout<<"sigma_t trk2= "<<sigma_t[jj][j][dspn2]<<endl;  
           
             //  
             // invert the bad  
             //  
             if(calibdata1->DSPbad_par[i][ii]==0)  
                 bad[jj][j][dspn1]=1;  
             else if(calibdata1->DSPbad_par[i][ii]==1)  
                 bad[jj][j][dspn1]=0;  
             if(calibdata2->DSPbad_par[i][ii]==0)  
                 bad[jj][j][dspn2]=1;  
             else if(calibdata2->DSPbad_par[i][ii]==1)  
                 bad[jj][j][dspn2]=0;  
         };  
     };  
           
     tr1->Delete();  
     tr2->Delete();  
 };  
 /**  
  * Method to fill the tracker calibration-parameter struct from on-line calibrations  
  */  
 void cTrkCalib::FillTCalibFrom(TFile* f0, Int_t ev_reg1, Int_t ev_reg2){  
   
     TTree *tr1 = (TTree*)f0->Get("CalibTrk1");    
     TTree *tr2 = (TTree*)f0->Get("CalibTrk2");  
   
     CalibTrk1Event *calibdata1 = 0;  
     CalibTrk2Event *calibdata2 = 0;  
     tr1->SetBranchAddress("CalibTrk1",&calibdata1);  
     tr2->SetBranchAddress("CalibTrk2",&calibdata2);  
     tr1->GetEntry(ev_reg1);  
     tr2->GetEntry(ev_reg2);  
     Int_t dspn1(0),dspn2(0);  
           
     for(Int_t i=0;i<6;i++){  
         dspn1=calibdata1->DSPnumber[i]-1;  
         dspn2=calibdata2->DSPnumber[i]-1;  
         /*     cout<<"dspn1= "<<dspn1<<"   dspn2= "<<dspn2<<endl; */  
         for(Int_t ii=0;ii<3072;ii++){  
             Int_t j=(ii/128);  
             Int_t jj=ii-j*128;  
             pedestal_t[jj][j][dspn1]=aint(calibdata1->DSPped_par[i][ii]);  
             pedestal_t[jj][j][dspn2]=aint(calibdata2->DSPped_par[i][ii]);  
             sigma_t[jj][j][dspn1]=anint(calibdata1->DSPsig_par[i][ii]);  
             sigma_t[jj][j][dspn2]=anint(calibdata2->DSPsig_par[i][ii]);  
         };  
     };  
           
     tr1->Delete();  
     tr2->Delete();  
 };  
 /**  
  * Method to fill the tracker calibration-parameter struct from default calibration  
  */  
 void cTrkCalib::FillTCalibFrom(TString path){  
       
     path_.FillWith(path);  
     fillpedsigfromdefault_();  
     if(path_.error) throw -216;  
   
 };  
   
 /**  
  * Method to fill the tracker VA1-mask struct from calibration  
  */  
 void cTrkMask::Set(TFile* f0, Int_t ev_reg1, Int_t ev_reg2){  
       
     TrkCalib *calib = new TrkCalib();  
   
     for(Int_t i=0; i<2; i++){  
         TTree *tr=0;  
         if(i==0){  
             //-------------------------  
             //read calibration packet 1  
             //-------------------------  
             tr = (TTree*)f0->Get("CalibTrk1");    
             tr->SetBranchAddress("CalibTrk1",calib->GetPointerTo());  
             tr->GetEntry(ev_reg1);  
         }  
         if(i==1){  
             //-------------------------  
             //read calibration packet 2  
             //-------------------------  
             tr = (TTree*)f0->Get("CalibTrk2");  
             tr->SetBranchAddress("CalibTrk2",calib->GetPointerTo());  
             tr->GetEntry(ev_reg2);  
         }  
         for(Int_t ipkt=0; ipkt<NPLANE; ipkt++){  
             for(Int_t ivk=0; ivk<NVK; ivk++){  
                 Int_t view   = calib->GetView(ipkt);  
                 Float_t sigm = calib->GetMean("SIG-BAD",ipkt,ivk+1);  
                 Float_t sigmin = 0.;  
                 if(view%2)sigmin = 6.5;  
                 else      sigmin = 2.5;  
                 Int_t mask = 1;  
                 if( sigm<sigmin )mask = 0;  
               
                 mask_vk_run[ivk][view-1] = mask;  
               
 //              if(mask==0)cout<< " VIEW "<<view<<" VK "<<ivk+1<<" -- <SIG> "<<sigm<<" ---> MASKED!"<<endl;  
             }    
         }  
         if(tr)tr->Delete();  
     }  
   
     if(calib)calib->Delete();  
   
 };  
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.23