/** * Packet Scan * author Marcelli * version 1.0 - 2 February 2006 * * Description: This script creates a Multigraph to show packetID variable (for all packets) vs. OBT with an axplicative legend * ** Version 1.25 - 13 January 2005 * Changed Int_t to Float because of variable range size * (UInt_t has been excluded beacuse of uncompatibility with TGraph) * * version 1.3 - 22 February 2005 - Nagni * For compatibility with batch mode excution: * 1) Added "include " and "using namespace std" * 2) Removed gROOT->Reset() * * Version 1.4 * Date 08 March 2005 - Nagni * Added "format" parameter for saving the produced image in various formats * (for a complete list of types refer to TPad::SaveAs method) * * Version 1.5 * Date 09 February 2006 - Marcelli * Update to work with new Yoda output * * * Version 1.6 * Date 27 February 2006 - Marcelli * For compilation: * Added function "int main(int argc, char* argv[])" * * * Description: To show packet type recorded by PAMELA (Packet type vs. OBT) * * Parameters: * TSTring base - the path to the root directory for the specific Pamela unpack session * There is no default value, without this input the program will not run * TString outDir - the path where to save the output image (Default = ./) * TString format - the format which will be used fo rsave the produced images (Default = "jpg") * */ #include #include #include #include "TString.h" #include "TFile.h" #include "TGraph.h" #include "TTree.h" #include "TMultiGraph.h" #include "TCanvas.h" #include "TLegend.h" #include "TObjString.h" #include "EventHeader.h" #include "PscuHeader.h" using namespace std; void PacketScan(TString base, TString outDir, TString format) { TFile *file = new TFile(base.Data()); if (!file){ printf("No such file in the directory has been found"); return; } if (outDir == "" ) outDir = "."; TList *list=file->GetListOfKeys(); Int_t numkey = file->GetNkeys(); pamela::EventHeader *eh = 0; pamela::PscuHeader *ph = 0; Float_t obt[1]={0}; Float_t id[1]={0}; Int_t cont[48]={0}; std::stringstream oss, oss1, oss2, oss3; TString filename = ((TObjString*)base.Tokenize('/')->Last())->GetString(); filename = ((TObjString*)filename.Tokenize('.')->First())->GetString(); oss.str(""); oss << "PACKET ID - OBT: " << filename.Data(); TCanvas *finalCanv1 = new TCanvas("Packet_1", "PacketScan1", 1280, 1024); finalCanv1->SetFillColor(10); TMultiGraph *mg1 = new TMultiGraph(); mg1->SetTitle(oss.str().c_str()); TLegend *leg1b = new TLegend(0.80,0.67,0.99,0.99, ""); TCanvas *finalCanv2 = new TCanvas("Packet_2", "PacketScan2", 1280, 1024); finalCanv2->SetFillColor(10); TMultiGraph *mg2 = new TMultiGraph(); mg2->SetTitle(oss.str().c_str()); TLegend *leg2 = new TLegend(0.87,0.70,0.99,0.99, ""); TCanvas *finalCanv3 = new TCanvas("Packet_3", "PacketScan3", 1280, 1024); finalCanv3->SetFillColor(10); TMultiGraph *mg3 = new TMultiGraph(); mg3->SetTitle(oss.str().c_str()); TLegend *leg3 = new TLegend(0.87,0.65,0.99,0.99, ""); for (Int_t i=0; iAt(i); char *name=(char *)(key->GetName()); TTree* tr = (TTree*)file->Get(name); while((tr->IsZombie())){ i++; TObject *key = list->At(i); char *name=(char *)(key->GetName()); TTree* tr = (TTree*)file->Get(name); } Long64_t nevents = tr->GetEntries(); const Int_t size = nevents; tr->SetBranchAddress("Header", &eh); for (Int_t j = 0; j < size; j++){ tr->GetEntry(j); ph = eh->GetPscuHeader(); obt[0]= ph->GetOrbitalTime(); id[0]=(ph->GetPacketId1()); Int_t type = (int)id[0]; switch(type){ //physic block case 7: { TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(2); graph->SetMarkerStyle(21); graph->SetMarkerSize(0.7); mg1->Add(graph); if (cont[0]==0) leg1b->AddEntry(graph,"PhysEndRun (0x07)","p"); cont[0]++; break; } case 16:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(2); graph->SetMarkerStyle(21); graph->SetMarkerSize(0.6); mg1->Add(graph); if (cont[1]==0) leg1b->AddEntry(graph,"Physics (0x10)","p"); cont[1]++; break; } case 32:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(4); graph->SetMarkerStyle(21); graph->SetMarkerSize(1); mg1->Add(graph); if (cont[2]==0) leg1b->AddEntry(graph,"RunHeader/Trailer (0x20/21)","p"); cont[2]++; break; } case 33:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(4); graph->SetMarkerStyle(21); graph->SetMarkerSize(1); mg1->Add(graph); if (cont[2]==0) leg1b->AddEntry(graph,"RunHeader/Trailer (0x20/21)","p"); cont[2]++; break; } //calibration block case 8:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(2); graph->SetMarkerStyle(23); graph->SetMarkerSize(1.3); mg1->Add(graph); if (cont[23]==0) leg1b->AddEntry(graph,"CalibCalPulse1 (0x08)","p"); cont[23]++; break; } case 9:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(3); graph->SetMarkerStyle(23); graph->SetMarkerSize(1.3); mg1->Add(graph); if (cont[24]==0) leg1b->AddEntry(graph,"CalibCalPulse2 (0x09)","p"); cont[24]++; break; } case 17:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(4); graph->SetMarkerStyle(23); graph->SetMarkerSize(1.3); mg1->Add(graph); if (cont[25]==0) leg1b->AddEntry(graph,"CalibTrkBoth (0x11)","p"); cont[25]++; break; } case 18:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(5); graph->SetMarkerStyle(23); graph->SetMarkerSize(1.3); mg1->Add(graph); if (cont[26]==0) leg1b->AddEntry(graph,"CalibTrk1 (0x12)","p"); cont[26]++; break; } case 19:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(6); graph->SetMarkerStyle(23); graph->SetMarkerSize(1.3); mg1->Add(graph); if (cont[27]==0) leg1b->AddEntry(graph,"CalibTrk2 (0x13)","p"); cont[27]++; break; } case 29:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(30); graph->SetMarkerStyle(23); graph->SetMarkerSize(1.3); mg1->Add(graph); if (cont[28]==0) leg1b->AddEntry(graph,"Calib_Cal (0x1d)","p"); cont[28]++; break; } case 22:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(9); graph->SetMarkerStyle(23); graph->SetMarkerSize(1.3); mg1->Add(graph); if (cont[29]==0) leg1b->AddEntry(graph,"CalibTof (0x16)","p"); cont[29]++; break; } case 23:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(7); graph->SetMarkerStyle(23); graph->SetMarkerSize(1.3); mg1->Add(graph); if (cont[30]==0) leg1b->AddEntry(graph,"CalibS4 (0x17)","p"); cont[30]++; break; } case 24:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(49); graph->SetMarkerStyle(23); graph->SetMarkerSize(1.3); mg1->Add(graph); if (cont[31]==0) leg1b->AddEntry(graph,"CalibCalPed (0x18)","p"); cont[31]++; break; } case 25:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(33); graph->SetMarkerStyle(23); graph->SetMarkerSize(1.3); mg1->Add(graph); if (cont[32]==0) leg1b->AddEntry(graph,"Calib1_Ac1 (0x19)","p"); cont[32]++; break; } case 26:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(21); graph->SetMarkerStyle(23); graph->SetMarkerSize(1.3); mg1->Add(graph); if (cont[33]==0) leg1b->AddEntry(graph,"Calib2_Ac1 (0x1a)","p"); cont[33]++; break; } case 27:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(16); graph->SetMarkerStyle(23); graph->SetMarkerSize(1.3); mg1->Add(graph); if (cont[34]==0) leg1b->AddEntry(graph,"Calib1_Ac2 (0x1b)","p"); cont[34]++; break; } case 28:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(37); graph->SetMarkerStyle(23); graph->SetMarkerSize(1.3); mg1->Add(graph); if (cont[35]==0) leg1b->AddEntry(graph,"Calib2_Ac2 (0x1c)","p"); cont[35]++; break; } case 34:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(41); graph->SetMarkerStyle(23); graph->SetMarkerSize(1.3); mg1->Add(graph); if (cont[36]==0) leg1b->AddEntry(graph,"CalibHeader/Trailer (0x22/23)","p"); cont[36]++; break; } case 35:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(41); graph->SetMarkerStyle(23); graph->SetMarkerSize(1.3); mg1->Add(graph); if (cont[36]==0) leg1b->AddEntry(graph,"CalibHeader/Trailer (0x22/23)","p"); cont[36]++; break; } //initialization block case 36:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(20); graph->SetMarkerStyle(22); graph->SetMarkerSize(1.3); mg1->Add(graph); if (cont[4]==0) leg1b->AddEntry(graph,"InitHeader/Trailer (0x24/25)","p"); cont[4]++; break; } case 37:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(20); graph->SetMarkerStyle(22); graph->SetMarkerSize(1.3); mg1->Add(graph); if (cont[4]==0) leg1b->AddEntry(graph,"InitHeader/Trailer (0x24/25)","p"); cont[4]++; break; } case 48:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(3); graph->SetMarkerStyle(21); graph->SetMarkerSize(1); mg2->Add(graph); if (cont[38]==0) leg2->AddEntry(graph,"EventTrk (0x30)","p"); cont[38]++; break; } //test block case 64:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(2); graph->SetMarkerStyle(22); graph->SetMarkerSize(1); mg2->Add(graph); if (cont[39]==0) leg2->AddEntry(graph,"Test_trk (0x40)","p"); cont[39]++; break; } case 65:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(3); graph->SetMarkerStyle(22); graph->SetMarkerSize(1); mg2->Add(graph); if (cont[40]==0) leg2->AddEntry(graph,"Test_Tof (0x41)","p"); cont[40]++; break; } //Log block case 80:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(3); graph->SetMarkerStyle(21); graph->SetMarkerSize(1); mg2->Add(graph); if (cont[41]==0) leg2->AddEntry(graph,"Log (0x50)","p"); cont[41]++; break; } //Arrdump-Vardump-Tabdump block case 81:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(2); graph->SetMarkerStyle(23); graph->SetMarkerSize(1); mg2->Add(graph); if (cont[42]==0) leg2->AddEntry(graph,"VarDump (0X51)","p"); cont[42]++; break; } case 82:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(3); graph->SetMarkerStyle(23); graph->SetMarkerSize(1); mg2->Add(graph); if (cont[43]==0) leg2->AddEntry(graph,"ArrDump (x52)","p"); cont[43]++; break; } case 83:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(4); graph->SetMarkerStyle(23); graph->SetMarkerSize(1); mg2->Add(graph); if (cont[44]==0) leg2->AddEntry(graph,"TabDump (0x53)","p"); cont[44]++; break; } //Tmtc Block case 84:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(2); graph->SetMarkerStyle(21); graph->SetMarkerSize(1); mg2->Add(graph); if (cont[45]==0) leg2->AddEntry(graph,"Tmtc (0x54)","p"); cont[45]++; break; } //Mcmd Block case 85:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(4); graph->SetMarkerStyle(21); graph->SetMarkerSize(1); mg2->Add(graph); if (cont[46]==0) leg2->AddEntry(graph,"Mcmd (0x55)","p"); cont[46]++; break; } case 96:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(20); graph->SetMarkerStyle(21); graph->SetMarkerSize(1); mg2->Add(graph); if (cont[22]==0) leg2->AddEntry(graph,"ForceFECmd (0x60)","p"); cont[22]++; break; } case 112:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(1); graph->SetMarkerStyle(22); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[6]==0) leg3->AddEntry(graph,"Ac1Init (0x70)","p"); cont[6]++; break; } case 113:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(2); graph->SetMarkerStyle(22); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[7]==0) leg3->AddEntry(graph,"CalInit (0x71)","p"); cont[7]++; break; } case 114:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(3); graph->SetMarkerStyle(22); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[8]==0) leg3->AddEntry(graph,"TrkInit (0x72)","p"); cont[8]++; break; } case 115:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(4); graph->SetMarkerStyle(22); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[9]==0) leg3->AddEntry(graph,"TofInit (0x73)","p"); cont[9]++; break; } case 116:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(5); graph->SetMarkerStyle(22); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[10]==0) leg3->AddEntry(graph,"TrgInit (0x74)","p"); cont[10]++; break; } case 117:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(6); graph->SetMarkerStyle(22); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[11]==0) leg3->AddEntry(graph,"NdInit (0x75)","p"); cont[11]++; break; } case 118:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(7); graph->SetMarkerStyle(22); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[12]==0) leg3->AddEntry(graph,"S4Init (0x76)","p"); cont[12]++; break; } case 119:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(32); graph->SetMarkerStyle(22); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[13]==0) leg3->AddEntry(graph,"Ac2Init (0x77)","p"); cont[13]++; break; } //Alarm block case 129:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(1); graph->SetMarkerStyle(21); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[14]==0) leg3->AddEntry(graph,"CalAlarm (0x81)","p"); cont[14]++; break; } case 130:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(42); graph->SetMarkerStyle(21); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[15]==0) leg3->AddEntry(graph,"Ac1Alarm (0x82)","p"); cont[15]++; break; } case 131:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(3); graph->SetMarkerStyle(21); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[16]==0) leg3->AddEntry(graph,"TrkAlarm (0x83)","p"); cont[16]++; break; } case 132:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(4); graph->SetMarkerStyle(21); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[17]==0) leg3->AddEntry(graph,"TrgAlarm (0x84)","p"); cont[17]++; break; } case 133:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(5); graph->SetMarkerStyle(21); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[18]==0) leg3->AddEntry(graph,"TofAlarm (0x85)","p"); cont[18]++; break; } case 134:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(6); graph->SetMarkerStyle(21); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[19]==0) leg3->AddEntry(graph,"S4Alarm (0x86)","p"); cont[19]++; break; } case 136:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(7); graph->SetMarkerStyle(21); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[20]==0) leg3->AddEntry(graph,"Ac2Alarm (0x88)","p"); cont[20]++; break; } /*case 137:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(32); graph->SetMarkerStyle(21); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[21]==0) leg3->AddEntry(graph,"ac2_alarm (137=0x89)","p"); cont[21]++; break; } case 138:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(20); graph->SetMarkerStyle(21); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[22]==0) leg3->AddEntry(graph,"s4_adc_ac_2 (138=0x8a)","p"); cont[22]++; break; }*/ case 161:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(2); graph->SetMarkerStyle(3); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[3]==0) leg3->AddEntry(graph,"TsbT (0xA1)","p"); cont[3]++; break; } case 171:{ TGraph *graph = new TGraph(1, obt, id); graph->SetMarkerColor(4); graph->SetMarkerStyle(3); graph->SetMarkerSize(1.3); mg3->Add(graph); if (cont[5]==0) leg3->AddEntry(graph,"TsbB (0xAB)","p"); cont[5]++; break; } } } } //*******************************TO DRAW AND SAVE*****************************************************/// finalCanv1->cd(); mg1->SetMinimum(0); mg1->SetMaximum(55); mg1->Draw("AP"); mg1->GetXaxis()->SetTitle("OBT (ms)"); mg1->GetXaxis()->CenterTitle(); mg1->GetYaxis()->SetTitle("Packet ID"); mg1->GetYaxis()->CenterTitle(); leg1b->Draw(); oss1.str(""); oss1 << outDir.Data() << filename.Data() << "_PacketScan_1." << format.Data(); finalCanv1->Update(); finalCanv1->SaveAs(oss1.str().c_str()); finalCanv2->cd(); mg2->SetMinimum(60); mg2->SetMaximum(100); mg2->Draw("AP"); mg2->GetXaxis()->SetTitle("OBT (ms)"); mg2->GetXaxis()->CenterTitle(); mg2->GetYaxis()->SetTitle("Packet ID"); mg2->GetYaxis()->CenterTitle(); leg2->Draw(); oss2.str(""); oss2 << outDir.Data() << filename.Data() << "_PacketScan_2." << format.Data(); finalCanv2->Update(); finalCanv2->SaveAs(oss2.str().c_str()); finalCanv3->cd(); mg3->SetMinimum(110); mg3->SetMaximum(180); mg3->Draw("AP"); mg3->GetXaxis()->SetTitle("OBT (ms)"); mg3->GetXaxis()->CenterTitle(); mg3->GetYaxis()->SetTitle("Packet ID"); mg3->GetYaxis()->CenterTitle(); leg3->Draw(); oss3.str(""); oss3 << outDir.Data() << filename.Data() << "_PacketScan_3." << format.Data(); finalCanv3->Update(); finalCanv3->SaveAs(oss3.str().c_str()); file->Close(); } int main(int argc, char* argv[]){ TString path; TString outDir = "./"; TString format = "jpg"; if (argc < 2){ printf("You have to insert at least the file to analyze \n"); printf("Try '--help' for more information. \n"); exit(1); } if (!strcmp(argv[1], "--help")){ printf( "Usage: PacketScan FILE [OPTION] \n"); printf( "\t --help Print this help and exit \n"); printf( "\t -outDir[path] Path where to put the output [default ./] \n"); printf( "\t -format[ps|gif|jpg] Format for output files [default 'jpg'] \n"); exit(1); } path=argv[1]; for (int i = 2; i < argc; i++){ if (!strcmp(argv[i], "-outDir")){ if (++i >= argc){ printf( "-outDir needs arguments. \n"); printf( "Try '--help' for more information. \n"); exit(1); } else{ outDir = argv[i]; continue; } } if (!strcmp(argv[i], "-format")){ if (++i >= argc){ printf( "-format needs arguments. \n"); printf( "Try '--help' for more information. \n"); exit(1); } else{ format = argv[i]; continue; } } } PacketScan(argv[1], outDir, format); }